Таблица 1.3. 1 кривая Верхняя граница 0.Нижняя граница -0.2 кривая Верхняя граница 0.Нижняя граница -0.3 кривая Верхняя граница 0.Нижняя граница -0.Опираясь на метод имитационного моделирования, построена выборка, каждый элемент которой представляет собой вектор.
Координата вектора - это доля вложений в одну из сфер деятельности. Таблица 1.3.4 представляет собой основную рабочую таблицу задачи. Ее последний столбец содержит ожидаемое значение дохода для соответствующего варианта распределения, вычисляемое по формуле (1.3.10).
Таблица 1.3.Фрагмент выборки для метода Монте-Карло Распределение средств между 3 видами деятельности x1 x2 x3 всего доход 1- доход 2- доход 3- Валовый вложен го го го доход ий направле направл направле ния ения ния 337 777 525 1639 625.0 664.6 1138.7 2428.1842 113 44 1999 871.4 515.1 26.9 1413.456 1320 13 1789 651.5 664.8 14.6 1331.135 1667 70 1872 577.9 664.8 44.5 1287.1339 61 421 1821 811.7 447.2 1117.2 2376.580 777 551 1908 678.1 664.6 1140.0 2482.1829 74 3 1906 870.1 465.7 12.0 1347.1394 421 80 1895 819.3 655.8 53.9 1529.1345 215 217 1777 812.5 601.2 495.7 1909.448 1133 269 1850 649.8 664.8 781.7 2096.1732 232 30 1994 860.2 610.2 20.4 1490.81 984 787 1852 565.0 664.8 1141.8 2371.1558 310 126 1994 840.3 638.4 126.2 1605.860 236 234 1330 733.6 612.2 592.3 1938.1374 294 103 1771 816.6 634.1 83.1 1533.659 648 343 1650 694.5 663.9 1033.5 2391.1180 591 72 1843 788.2 663.2 46.2 1497.1863 46 69 1978 873.4 425.0 43.7 1342.1247 669 47 1963 798.4 664.1 28.5 1490.Максимальное значение дохода составило 2559.1 и соответствует распределению суммы 1974 ед. В следующей пропорции (1 112; 418; 444). Следующий наиболее выгодный вариант вложений (доход - 2550.2 ед.) предполагает распределение 1990 ед. между тремя направлениями в соотношении (1219, 369, 402), т.е. каждому направлению и в первом и во втором случае соответствует примерно одинаковая доля общей суммы.
Второй подход. Поскольку второй вариант постановки стохастической задачи основан на двухкритериальной задаче, имеющей следующие целевые функции доход max, (1.3.11) риск min, то обратим внимание на возможности количественной оценки указанных критериев и подходы к решению, связанные с векторным характером задачи.
Особенностью этой задачи, как и любой задачи векторной оптимизации, является принятие решения, последствия которого станут до конца ясными лишь в будущем. Как и в любой жизненной ситуации, эти последствия не могут быть объективно оценены при помощи математических расчетов. Каждый получаемый, эффективный по Парето, вариант решения имеет оценку по двум противоречивым критериям, расчеты не могут определить компромисс между этими критериями. Но поскольку выбор одного из вариантов должен быть осуществлен в любом случае, то неопределенность, связанная с наличием двух критериев, может быть устранена только на основе информации, полученной от экспертов или же с помощью внешнего критерия, не входящего в состав целевых функций задачи.
Количественная оценка критерия дохода определяется, согласно описанному выше подходу, как сумма прибылей, полученных во всех сферах деятельности. Этот подход, отличный от классической оценки доходности портфеля вложений как средневзвешенной доходностей его составляющих, кажется разумным для моделирования зависимости затраты -выпуск, поскольку развитие производства в любой отрасли в некоторой степени соответствует жизненному циклу, при этом существуют пороговые значения сумм инвестиций, при переходе через которые изменяется величина прибыли на единицу вложений.
При определении оценки риска портфеля вложений кроме классического метода измерения риска с помощью ковариаций n n n n n 2 Cij X X = X + Cij X X i i i j i j i=1 j =1 i=1 i=1 j=1, ji могут быть использованы и альтернативные подходы. Например, определение риска с использованием понятий теории информации, позволяющие оценить содержание, полезность полученных данных для их потребителя. Применительно к проблеме распределения средств: информация, полученная от экспертов, может быть названа полезной, если после того, как управляющий учел ее при принятии решения, были достигнуты более высокие результаты. При анализе распределения инвестиций будем использовать понятия теории информации: неопределенность и энтропия. Величину n n * N( X / X ) = - xi ln x*, где xi > 0 и xi =1, i i=1 i=(если брать log2, то N(X*/X) измеряется в битах) называют неопределенностью для ситуации, когда истинное распределение есть X*, а человек, решающий задачу, думает, что оно есть X.
Изменение неопределенности можно интерпретировать как процесс запасания полезной информации. За нулевой уровень часто удобно принимать запас полезной информации при xi =1 n, i =1,n.
Если мы ничего не знаем о распределении X*, то в некотором смысле гипотеза xi =1 n является оптимальной: если мы не имеем сведений о событиях, объектах, то будем считать их равновероятными, имеющими одинаковые характеристики. При таком подходе запас полезной информации, содержащийся в гипотезе X относительно задачи с распределением ответа X*, дается соотношением * In = log( n ) - N( X / X ).
Т.о. полезная информация зависит не от точности, с которой известны X, а от того, насколько X соответствует X*.
Нижней границей неопределенности является энтропия распределения X*, равная n * H( X ) = - x* ln x* i i i=Модифицируем определение энтропии, нормируем ее:
" H ( X ) =1 - H( X ) / ln( n ) Величина H(X) принимает значения из интервала [0, 1), H(X) = 0 при xi =1 n для i = 1,n. Ее значение можно рассматривать как меру риска портфеля вложений. H(X) показывает степень отклонения распределения X от равномерного распределения (1/n, 1/n,..., 1/n). При равномерном распределении средств велика вероятность того, что выбор сделан неправильно, и чем меньше H(X), то есть чем ближе распределение к равномерному, значит, тем больше возможность ошибиться при наборе портфеля, и тем больше риск; т.е. неопределенность рассматривается как мера риска, возможности сделать неправильный выбор.
Если мы хотим при измерении риска учитывать только вероятность нежелательных результатов инвестирования, то простейшим способом его измерения является вероятность недобора, измеряющая шансы на то, что доходность направления окажется ниже ожидаемой доходности. По существу, это доля вероятностного распределения, лежащая справа от ожидаемой доходности.
Более сложные способы измерения риска получения доходности ниже ожидаемой производятся с помощью семейства статистических данных, известных как частичные моменты низких порядков. Например, средний недобор измеряет среднее отклонение доходности ценной бумаги вниз от ожидаемой доходности. Он является более полезным, чем вероятность недобора, т.к. принимает во внимание величину каждого отрицательного отклонения и показывает какова может быть величина уменьшения относительно ожидаемой доходности.
Мерой риска, учитывающей лишь случаи снижения доходности по отношению к среднему значению, является полувариация. Она рассчитывается как обычная вариация кроме тех случаев, когда доходность выше ожидаемой.
Еще одним альтернативным методом измерения риска является определение величины потерь, такой, что потери в доходности портфеля за определенный период времени с заданной вероятностью не превысят этой величины (стоимостная оценка риска VAR).
Перейдем к формулировке задачи, в которой второй критерий - оценка рискованности предлагаемого варианта распределения средств - может определяться одним из описанных выше способом.
Постановка интервальной задачи. Приводимая ниже задача отличается от классических стохастических моделей планирования производства тем, что кроме областей возможного изменения параметров (оценок доходности каждой сферы деятельности;
величин, характеризующих рискованность вложений средств в каждое из рассматриваемых направлений) содержит еще и два критерия. Т.е. отражает как бы два вида неопределенности:
неточность, неполноту, частичную недоступность информации о внешней среде и противоречивость внутренних целей любой производственной деятельности. Указанное обстоятельство приводит к некоторым изменениям в применении классических критериев принятия решений.
Рассмотрим математическую постановку задачи. В реальной ситуации точные значения внешних переменных обычно неизвестны, поэтому возникает необходимость привлечения других методов, которые способны дать представление о значениях и закономерностях изменения переменных в течение того периода, на который строится прогноз. Далее предлагается подход к решению задачи распределения средств на основе модели, экзогенные параметры которой предварительно оценены по статистической информации.
Возьмем в качестве интервальной оценки будущей доходности сферы деятельности отрезок, задаваемый минимальным и максимальным ее значением в предыдущих периодах. При этом необходимо быть уверенным, что в следующем периоде существенных изменений на рынке не произойдет. Т. о.
сформируем вектор верхних и нижних границ min Eit Ei max Eit (1.3.12) t t или.
E E E Аналогичным образом определим границы для матрицы ковариаций:
C C C или cij cij cij, (1.3.13) t t где cij = min cij, cij = max cij.
t t Приходим к задаче, о коэффициентах которой известно лишь то, что они находятся в некотором множестве, отражающем все реальные ситуации:
T min X CT X, C C C,(1.3.14) X max ET X, E E E, (1.3.15) X n X = 1, X 0, i =1, 2,..., n. (1.3.16) i i i=Возможные алгоритмы решения. Так как ограничения (1.3.12)-(1.3.13) задают бесконечное число возможных реализаций будущей ситуации, исследовать которые все не представляется возможным, то упростим задачу (1.3.14) - (1.3.16), введя параметр [ 0,1], определяющий характеристики ситуаций следующим образом:
k k C = C + (1 - ) C, E = E + (1 - ) E, k k k k где = + h, h =, = 0, k =1, N. (1.3.17) k k -1 N Задача (1.3.14) - (1.3.16) решается N раз при соответствующих k k коэффициентах целевых функций C и E. При отсутствии экспертных процедур выбора единственного решения из получаемых на каждом шаге множеств эффективных точек может быть использована процедура выбора недоминируемых точек, заключающаяся в определении альтернативы, наиболее близкой к УоптимальнойФ по всем критериям. Т.е. этот подход позволяет сузить множество Парето-оптимальных точек до единственной.
Пусть множество критериев задачи векторной оптимизации имеет вид f = ( f1, f2,..., fm ), а - множество допустимых значений аргументов, то лидеальная точка определяется следующим образом. Каждая отдельная компонента F(X) имеет максимум при некотором X, предположим, fi ( X ) достигает своего экстремума при X. Можно записать, что * * * extr fi ( X ) = fi ( X ) = fi*. Тогда вектор f =(f1*, f2,..., fm ) есть лидеальная точка, т.е. вектор всех экстремальных допустимых значений, достигаемых отдельными целевыми функциями на множестве. Но такое идеальное решение чаще всего невозможно. Будем предполагать, что ЛПР стремится найти такое решение, которое было бы как можно ближе расположено к идеальной точке.
Выбор такого Укомпромиссного вариантаФ возможно и не является оптимальным, но оказывается наиболее предпочтительным по совокупности критериев.
Алгоритм Для двухкритериальной задачи (1.3.14)-(1.3.16) определяем лидеальную точку, решая две оптимизационные задачи:
T k min X C X n 1* из нее находим точку X и соответствующее xi =1, xi i=значение целевой функции f1* ;
kT max E X n 2* из нее находим точку X и соответствующее xi =1, xi i=* значение целевой функции f2.
Идеальная точка в пространстве критериев имеет * следующие координаты F* = ( f1*, f2 ).
Ищем решение задачи d( F( X ), F* ) min n, xi =1, xi i=где d( F( X ), F* ) расстояние между двумя точками в пространстве критериев, определяемое одним из следующих способов * d( F,F* ) =[(f1 - f1* ) + (f2 - f2 )] или * d( F,F* ) = f1 - f1* + f2 - f2.
После N повторений процедур (1)-(2) получаем серию точек k* {(, X )}, соответствующих различным реализациям внешних и k k k внутренних условий ( C, E ). Следующим сложным вопросом является процедура выбора определенной программы из этого множества точек.
Пример. В приведенном далее примере используются следующие исходные данные: таблица 1.3.5 содержит верхние и нижние границы координат вектора доходностей направлений деятельности; таблица 1.3.6 - содержит матрицы, элементы которых представляют собой максимальное и минимальное значение матрицы ковариаций доходностей.
Таблица 1.3.Е Е 10,7 0,20,4 0,31,2 0,Таблица 1.3.С С 1 2 31 2 1 0,25 0,15 0,5 1 0,5 0,6 0,2 0,15 0,3 0,2 2 0,6 0,65 0,30,5 0,2 0,7 3 0,89 0,45 0,Интервал для каждого нечетко заданного параметра разбит шестью внутренними точками (в формуле (1.3.17) выбран шаг = 0,125 ). Первоначальная выборка состояла из 111 наборов. Для каждого варианта распределения вычисляют значения дохода и рискованности вложений (4 и 5 столбцы рабочей таблицы).
Фрагмент рабочей таблицы Идеальное значение ф-ции цели max E min R min D 0,974217 0,44617 0,х1 х2 х3 Доходность Риск Расстояние 1 0,38 0,236 0,3810,741496 0,61478 0,2 0,10068 0,400 0,4990,732163 0,59126 0,3 0,59648 0,354 0,0480,578156 0,49446 0,4 0,89910 0,042 0,0580,683144 0,48018 0,5 0,8846 0,107 0,0080,643281 0,45354 0,6 0,95846 0,033 0,0070,666943 0,44617 0,7 0,01449 0,976 0,0080,360866 0,5574 0,8 0,40742 0,390 0,2010,623576 0,54158 0,9 0,86324 0,043 0,0930,696084 0,50243 0,10 0,13858 0,0929 0,7680,932991 0,77854 0,Определив максимальное значение дохода и минимальное значение рискованности по всем вариантам, получим лидеальную точку (0,97; 0,44). Далее при каждом для каждого варианта к распределения (X1, X2, X3) определяется точка, самая близкая к лидеальной. Расстояние между точками определяется стандартной евклидовой мерой. В результате получена таблица 1.3.7.
Таблица 1.3.№ точки № х1 х2 х3 в ER к п/п выборке 10 0,015 0,323 0,66178 0,933 0,2 0,125 0,015 0,323 0,66178 0,875 0,3 0,25 0,015 0,323 0,66178 0,817 0,4 0,5 0,863 0,043 0,0939 0,658 0,5 0,625 0,863 0,043 0,0939 0,639 0,6 0,75 0,958 0,033 0,0076 0,613 0,7 0,875 0,958 0,033 0,00776 0,599 0,81 0,958 0,033 0,00776 0,586 0,Описанные возможности моделирования факторов неопределенности реальных ситуаций перспективного планирования развития предприятия, основанные на сочетании традиционных и стохастических моделей с методами принятия решений, предоставляют широкие возможности поиска компромиссного варианта распределения имеющихся ресурсов.
Pages: | 1 | 2 | 3 | 4 | 5 | ... | 6 | Книги по разным темам