- •Глава 1. Принципы моделирования случайных элементов 6
- •Глава 3. Метод Монте-Карло и его применения 51
- •3.4. Задания 63
- •Глава 4. Система моделирования gpss 64
- •Глава 5. Моделирование сложных систем на основе gpss 127
- •5.2.1. Иммитационная модель системы 136
- •Глава 6. Задания лабораторного практикума по gpss 152
- •Глава 7. Статистическое моделирование в финансовом анализе. 181
- •Введение
- •Глава 1.Принципы моделирования случайных элементов
- •1.1. Моделирование случайных величин
- •1.1.1.Моделирование базовых случайных величин (бсв) Понятие бсв
- •Датчики бсв
- •Метод Макларена-Марсальи
- •Специальные тесты проверки точности моделирования бсв
- •Тест «совпадения моментов»
- •Тест «ковариация»
- •Тест «равномерность двухмерного распределения»
- •1.1.2.Статический анализ точности моделирования
- •1.1.3.Задания
- •1.2. Моделирование дискретных случайных величин (дсв)
- •1.2.1.Понятие дсв. Общая схема алгоритма моделирования
- •1.2.3.Графический анализ точности моделирования
- •1.2.4.Задания
- •Дискретное равномерное распределение
- •1.3. Моделирование непрерывных случайных величин (нсв)
- •1.3.1.Понятие нсв. Универсальные методы моделирования нсв. Понятие нсв.
- •Метод обратной функции.
- •Метод исключения
- •1.3.2.Тесты проверки точности моделирования нсв Критерии согласия
- •- Критерий согласия Пирсона
- •Критерий согласия Колмогорова.
- •Критерий серий
- •Критерий “нисходящих” и ”восходящих” серий
- •1.3.3.Алгоритмы моделирования для основных непрерывных распределений Равномерное распределение
- •Алгоритм моделирования
- •Одномерное нормальное распределение
- •Алгоритм моделирования
- •Экспоненциальное распределение
- •Алгоритм моделирования
- •Распределение Лапласа
- •Алгоритм моделирования
- •Алгоритмы моделирования:
- •Бета-распределение.
- •Связь с другими распределениями
- •Алгоритмы моделирования:
- •Алгоритмы моделирования
- •Алгоритм моделирования
- •Моделирование случайных величин с заданной гистограммой
- •Алгоритм моделирования
- •Моделирование случайных величин с заданным полигоном частот
- •1.3.4.Задания Равномерное распределение
- •Одномерное нормальное распределение
- •Распределение Лапласа
- •Распределение Вейбулла-Гнеденко
- •Гамма распределение
- •Бета-распределение
- •Распределение Коши
- •Логистическое распределение
- •Распределение Стьюдента
- •Распределение Фишера
- •3.2. Вычисление определенного интеграла методом Монте-Карло
- •Пример.
- •Текст программы:
- •3.3. Решение системы алгебраических уравнений методом Монте-Карло
- •3.4. Задания
- •Глава 4.Система моделирования gpss
- •4.1. Краткие сведения из теории
- •4.2. Основные правила и операторы языка gpss
- •4.2.1.Структура операторов gpss
- •4.2.2.Основные операторы языка gpss
- •4.3. Основные блоки gpss/pc и связанные с ними объекты
- •4.3.1.Блоки, связанные с транзактами
- •4.3.2.Блоки, связанные с аппаратными объектами
- •4.3.3.Блоки для сбора статистических данных
- •4.3.4.Блоки, изменяющие маршруты транзактов
- •4.3.5.Блоки, работающие с памятью
- •4.3.6.Блоки для работы со списками пользователя
- •4.4. Управляющие операторы gpss/pc
- •4.5. Некоторые приемы конструирования gpss-моделей
- •4.5.1.Косвенная адресация
- •4.5.2.Обработка одновременных событий
- •4.6. Команды gpss/pc и технология работы с пакетом
- •4.6.1.Загрузка интегрированной среды
- •4.6.2.Ввод новой модели
- •4.6.3.Редактирование текста модели
- •4.6.4.Запись и считывание модели с диска
- •4.6.5.Прогон модели и наблюдение за моделированием
- •4.6.6.Получение и интерпретация стандартного отчета
- •5.1.3.Анализ результатов иммитационного моделирования
- •5.2. Моделирование работы участка цеха
- •5.2.1.Иммитационная модель системы
- •110 Generate 30,5 Первый сегмент модели
- •280 Generate 20,5 Второй сегмент модели
- •450 Generate 480 Третий сегмент модели (таймер)
- •5.2.2.Анализ результатов иммитационного моделирования
- •Обсуждение результатов моделирования
- •105 A2 storage 2 Второй станок моделируется как накопитель
- •5.3. Порядок оформления отчета
- •5.3.1.Содержание отчета
- •Задания для самопроверки
- •5.3.2.Создание выходного файла
- •Задание для домашней подготовки
- •5.3.3.Описание элементов выходного файла
- •Глава 6.Задания лабораторного практикума по gpss
- •Глава 7.Статистическое моделирование в финансовом анализе.
- •7.1. Моделирование наращенных сумм и современных величин платежа.
- •Задание 4.
- •7.2. Моделирование величины реальной наращенной суммы платежей с учетом инфляции.
- •7.3. Моделирование обменного курса валют при двойной конверсии валют.
- •Приложение
- •Примечание к таблице п.3
- •Литература
Глава 7.Статистическое моделирование в финансовом анализе.
В детерминированном финансовом анализе параметры финансовой операции такие как, первоначальная сумма, продолжительность ссуды, ставка процентов, считаются заданными. Теория таких расчетов хорошо развита [ ].
Однако в реальных ситуациях часть из этих параметров, а то и все параметры в целом, могут быть не полностью определены в силу различных причин. Чтобы сделать обоснованные выводы о ходе таких финансовых операций, необходимы дополнительные априорные предположения о поведении таких не вполне определенных параметров. В данном разделе мы будем считать, что не вполне определенные параметры финансовой операции являются случайными величинами с заданными, либо спрогнозированными законами распределения вероятностей. Тогда и результат финансовой операции (наращенная сумма, современная величина) будут также случайными величинами и для их оценки можно успешно применять аппарат статистического моделирования. Исследователя будут интересовать такие числовые характеристики, например, наращенной суммы, как минимальное и максимальное значение, среднее значение, дисперсия и вероятность попадания наращенной суммы в заданный интервал. По результатам статистического моделирования эти числовые характеристики могут быть оценены.
В финансовой математике, в силу ряда причин [ ], наибольшее применение находят два типа распределений случайных величин: дискретное распределение и “трапецеидальное” распределение Tr(a,b;c,d), a<=c<=d<=b. Распределение Tr (a, b; c, d)-это абсолютно-непрерывное распределение с плотностью распределения заданной на интервале [a, b] и принимающей постоянное значение 2/(b+d-a-c) на интервале [c, d], линейно возврастающее на интервале [a, c] и линейно убывающее на интервале [d, b]. Если a = c, d = b, то ‘’трапецеидальное’’ распределение обращается в равномерное распределение на отрезке [a, b]. Если c = d, то "трапецеидальное’’ распределение трансформируется в “треугольное” распределение. Наконец, если c = d = (a+b)/2, то имеем стандартное “треугольное” распределение .
Эти типы распределений, имеющие хорошую экономическую интерпретацию, легко моделируются.
7.1. Моделирование наращенных сумм и современных величин платежа.
Приведем формулы для вычисления наращенных сумм и современных величин платежей для ставок протых, сложных процентов и банковской учетной ставки. Все величины, стоящие в правых частях ниже приведенных формул, являются параметрами финансовой операции. Часть из этих параметров может быть не полностью заданной. В этих случаях не полностью определенные параметры будем рассматривать, как случайные величины, и будем моделировать их числовые значения. После того как все значения неизвестных параметров будут смоделированы, можно по приведенным ниже формулам, вычислить результирующее значение финансовой операции. А по совокупности результирующих значений можно оценить все числовые характеристики финансовой операции, используя методы математической статистики.
Формула наращения по ставке простых процентов, изменяющейся кусочно - линейно от интервала к интервалу, имеет вид:
, (1)
где S - наращенная
сумма, P - первоначальная
сумма,
,
- ставка простых процентов, определенная
на интервале ns,
n1+n2+…+nk
= n. Каждая
из величин, стоящих в правой части (1),
может быть случайной величиной.
Если уравнение (1) разрешить относительно P, то получим современное значение Р наращенной суммы S:
(2)
Если наращение производится по ставкам сложных процентов, то формула (1) трансформируется в формулу
,
(3)
а формула (2) – в формулу :
.
(4)
Учет векселя с суммой S, поглощаемого за t дней до его погашения по годовой банковской ставке d, осуществляется по формуле
.
(5)
Задание 1.
В начале года на депозит положено 560
тысяч рублей сроком на один год. По
условиям финансового соглашения банк
имеет право изменить процентную ставку
в начале каждого квартала. Первоначальная
ставка процентов
известна
и равна 70 % годовых, ставки процентов
,соответственно,
во втором, третьем и четвертом кварталах
неизвестны заранее. Есть основания
считать, что
будет
равномерно распределена на интервале
[73 %, 78 %] ;
будет иметь “трапецеидальное”
распределение на интервале [73 %, 78 %] с
максимальным значением плотности
распределения на интервале [74 %, 76 %], а
ставка процентов
в четвертом квартале примет значение
79 % с вероятностью 0,4 и 80 % с вероятностью
0,6. причем
- независимые случайные величины.
По 10 реализациям ставок процентов
оценить числовые характеристики
наращенной суммы S:
Smin,
Smax,
E{S},
D{S},
P{ 97600
97900 }.
Решение. Обозначим через
= (
)
трехмерный вектор. Пусть в результате
моделирования получены 10 реализаций
этого вектора:
(1)
= (0,723; 0,748; 0,8),
(2)
= (0,746; 0,733; 0,79).
(3)
= (0,758; 0,756; 0,8),
(4)
= (0,734; 0.779; 0,8),
(5)
= (0,734; 0,741; 0,79),
(6)
= (0,729; 0,757; 0,8),
(7)
= (0,755; 0,736; 0,8),
(8)
= (0,739; 0,747; 0,8),
(9)
= (0,753; 0,758; 0,79),
(10)
= (0,742; 0,762; 0,8).
Используя формулу (1) с 1 = 0,7; n1 = n2 = n3 = n4 =0,25, получаем: S(1) = 975,94; S(2) = 975,66; S(3) = 981,96; S(4) = 981,82; S(5) = 975,1; S(6) = 978,04; S(7) = 978,74; S(8) = 978,04; S(9) = 980,14; S(10) = 980,56.
По выборочным значениям S(I),
=
,
оценим следующие числовые характеристики
S:
Smin
= 975,1 тыс.; Smax
= 980,56тыс.; E{S}
=
=978,6
тыс.; D{S}
= 6,276 тыс.; P{97600
97900} = 0,5.
Задание2.
По 1000 реализаций S
вычисленных по формуле (1), оценить
числовые характеристики S:
Smin,
Smax,
E{S},
D{S},
P{
},для
следующих значений параметров, входящих
в правую часть формулы (1):
=
4, n1 = n2
= n3 =
n4 = 0,25
года, P = 1230, S1
= 1335, S2
= 1340;
,
,
,
– независимые случайные величины,
1
R (8,2 %; 9 %),
2
R (8,4 %; 9,1 %),
-
3
8,6 %
8,7 %
8,8 %
9,2 %
P
0,1
0,3
0,5
0,1
-
4
8,7 %
8,9 %
9,3 %
P
0,1
0,6
0,3
2) = 3,P = 15800; n = 1 год, S1 = 16800,S2 = 17000, n1, n2 -независимые дискретные случайные величины:
-
n1, лет
1/12
1/4
P
0,7
0,3
-
n2, лет
1/4
1/3
1/2
P
0,5
0,4
0,1
Если n1 = 1/12, то 1 R (3 %; 4,2 %), если n1 = 1/4, то 1 R (5 %; 6,2 %), если n2 = 1/4, то 2 R (5 %; 6,2 %); n2=1/3, то 2 R (6 %; 7 %);
если n2 =
1/2, то
2
R (6,5 %; 7,8 %). Если n3
0,5, то
3
R
(5 %; 7,8 %); если n3
> 0,5, то
3
R (6 %; 8 %).
Задание3.
По 1000 реализаций S вычисленных по формуле (3), оценить числовые характеристики S: Smin, Smax, E{S}, D{S}, P{ }, для следующих значений параметров, входящих в правую часть формулы (3):
= 2, n1 = n2 = 2 года, P = 10000, S1 = 13900,S2 = 14100; ,– независимые случайные величины, 1 R(9 %; 10,2 %), 2 Tr (8 %,12 %; 9 %,10 %);
2) = 3, P = 20800; n = 6 лет, n1 = 2 года, 1 = 9 %, n2 - дискретная случайная величина:
-
n2, годы
1
2
P
0,8
0,2
если n2 = 1, то
-
2
9 %
9,1 %
9,3 %
P
0,7
0,2
0,1
если n2 = 2, то
-
2
9,2 %
9,3 %
9,4 %
P
0,6
0,3
0,1
если n3 = 3, то 3=10 %, если n3 = 2, то 3 = 9,2 %.
