Далее необходимо определить максимальное и минимальное значение ключевых факторов и задать характер распределения вероятностей. Значение данных интервалов по всем годам проекта было выдано экспертами. За характер распределения было взято нормальное, поскольку, практика риск-анализа показывает, что в подавляющем большинстве случаев при оценки риска пользуется нормальное распределение /9/.
На основе выбранного распределения проводится имитация ключевых факторов, причем количество имитаций должно быть таким, чтобы совокупность случайных пробных значений могла считаться репрезентативной. В данном случае это количество составило 500 имитаций.
Для генерации двух случайных чисел, распределенных по нормальному закону, использовался метод Морсальи-Брея. При этом математическое ожидание высчитывалось как середина интервала разброса числа, которое было получено, как уже говорилось, от экспертов. Дисперсия рассчитывалась как середина длины отрезка данного интервала.
Из результата имитационных попыток определялись необходимые значения величин математических ожиданий и дисперсий денежных потоков для каждого периода проекта, а также их ковариации.
Математическое ожидание чистого приведенного эффекта рассчитывалось по формуле (26)
(26)где I - первоначальные капиталовложения, денежные единицы,
N - число планово-учетных периодов проекта,
I - номер планово-учетного периода
E(Si) -математическое ожидание чистого денежного потоки в i-й учетный период проекта, денежные единицы,
к - коэффициент дисконтирования, доли единицы.
В предложенной модели в качестве критерия для учета неопределенности (риска) стратегического инвестиционного проекта было выбрано среднеквадратическое отклонение чистого приведенного эффекта (NPV) от его математического ожидания. Если по результатам расчетов будет выявлено, что дисперсия достаточно большая, то риск того, что ожидаемое значение NPV примет запланированное значение, будет также достаточно большим.
Понятие дисперсии, так же как и математического ожидания, функции распределения и ковариации не определено для нечетких чисел, поэтому использование формулы (22) вызывает определенные трудности. В результате для учета риска была рассмотрена разработанная Недосекиным А.О. оценка возможности того, что по результатам инвестиционного процесса значение NPV окажется ниже предустановленного граничного уровня /6/. Таким образом, первый критерий оценки риска проекта будет вычисляться по следующей формуле
(27)где
, (28) . (29)где Risk(G) – вероятность того, что значение чистого приведенного эффекта окажется ниже предустановленного граничного уровня, доли единицы,
G – уровень эффективности проекта, денежные единицы,
NPVmin – минимальное значение чистого приведенного эффекта из заданного диапазона, денежные единицы,
NPVmах – максимальное значение чистого приведенного эффекта из заданного диапазона, денежные единицы,
NPVav – среднее значение чистого приведенного эффекта из заданного диапазона, денежные единицы,
α1 – функция принадлежности нечеткого числа NPV, доли единицы.
Степень риска Risk(G) принимает значение от 0 до 1. Каждый инвестор, исходя из своих инвестиционных предпочтений, может классифицировать значение Risk(G), выделив для себя отрезок неприемлемых значений риска. Возможна также более подробная градация степени риска, которая будет рассмотрена позже.
Формула для оценки ликвидности стратегического инвестиционного проекта, в силу введения нечетких чисел в модель и описанных в литературном обзоре операций, примет следующий вид
(30)
где Rt1– минимальное значение коэффициента ликвидности в момент времени t, доли единицы,
Rt2– среднее (наиболее ожидаемое) значение коэффициента ликвидности в момент времени t, доли единицы,
Rt3– максимальное значение коэффициента ликвидности в момент времени t, доли единицы,
Si1– минимальное значение чистого денежного потока в i-й планово-учетный период, денежные единицы,
Si2– среднее значение чистого денежного потока в i-й планово-учетный период, денежные единицы,
Si3– максимальное значение чистого денежного потока в i-й планово-учетный период, денежные единицы,
a1 – минимальное значение безрисковой ставки дисконтирования, доли единицы,
a2 – среднее значение безрисковой ставки дисконтирования, доли единицы,
a3 – максимальное значение безрисковой ставки дисконтирования, доли единицы,
n – число планово-учетных периодов проекта,
I – номер планово-учетных периодов,
j – номер планово-учетного периода на момент реализации стратегического инвестиционного проекта,
NPVt - фактически полученная стоимость стратегического инвестиционного проекта (денежные потоки, полученные на момент времени t), денежные единицы,
I1 –минимальное значение первоначальных капиталовложений, денежные единицы,
I2 – максимальное значение первоначальных капиталовложений, денежные единицы,
I3 – среднее значение первоначальных капиталовложений, денежные единицы.
Критерий покрытия в рамках данной модели оценивается по следующей формуле:
(31)
где Сt1 – минимальное значение критерия покрытия в момент времени t, доли единицы,
Сt2 – максимальное значение критерия покрытия в момент времени t, доли единицы,
Аt1 - минимальное значение собственного капитала субъекта в момент времени t, денежные единицы,
Аt2 - максимальное значение собственного капитала субъекта в момент времени t, денежные единицы,
Zt1 - минимальное значение заемного капитала в момент времени t, денежные единицы,
Zt2 - максимальное значение заемного капитала в момент времени t, денежные единицы.
Для получения однозначного ответа об уровне риска инвестиционного проекта на основе полученных данных необходимо ввести процедуру принятия решения. Такая процедура описана Смирновым А.П. в /17/.
Оператор оценивает входные данные с помощью субъективных качественных понятий типа “много”, ”мало” и т.п. Эти качественные оценки значений переменных uформализуются с помощью так называемых лингвистических переменных /11/.
Модель управления в рассматриваемом случае есть модель связи между входными переменными u и выходной переменной v. Механизм этой связи включает суждения оператора о значениях переменных. В результате на основе численного значения каждой из входных переменных оператор присваивает им качественные (нечеткие) значения. Свое решение он также принимает на основе нечеткого значения выходной переменной. Это означает, что оператор интуитивно пользуется нечеткой логикой, а конкретно – правилами нечеткого вывода.
Правила вначале формулируются с помощью термов (словесных описаний значений входных переменных). Каждое правило представляет собой текст, определяющий некоторое нечеткое отношение R между входными переменными u и выходной переменной v. Обозначим порядковый номер правила через L.
Для превращения текста правила в формальную процедуру нужно установить вид правила композиционного вывода и форму нечеткой импликации.
В качестве правила композиционного вывода для рассматриваемого класса систем может быть принята максиминная композиция, а в качестве нечеткой импликации – правило минимума (пересечение нечетких множеств предпосылки и заключения).
Нечеткое отношение R для L-го правила между j-й входной переменной ujи выходной переменной v в соответствии с принятым правилом минимума выражено следующей функцией принадлежности
(32)Здесь индекс i(L) означает индекс i-го терма в L-м правиле вывода. Функция принадлежности (31) отображает отношение связи между числовыми значениями в паре (uj,v). Чем больше ее значение, тем теснее эта связь.
Результаты измерения (наблюдения) входных переменных могут быть выражены как обычными числовыми (четкими) значениями, так и качественными или размытыми значениями (нечеткими числами).
Пусть входные переменные uj представлены нечеткими числами fj с функциями принадлежности mfj(uj). Заметим, что эти функции есть результат работы системы наблюдения (измерения) в отличие от ранее введенных функций mji(uj), которые выражают мнение эксперта-оператора по поводу конкретных значений uj. Тогда в соответствии с принятым правилом композиционного вывода можно записать связь между выходной переменной v и входной переменной uj следующим образом
(33)Здесь Mj(v) есть функция принадлежности, устанавливающая локальную связь между нечеткой входной переменной uj и нечеткой выходной переменной v.
Если система наблюдения дает конкретные числовые значения uj=Ej, то формула (33) преобразуется к следующему виду