- •Глава 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
- •Литература
Тест «ковариация»
Ковариационной
функцией случайной последовательности
называется функция целочисленной
переменной
:
Если
--
независимые, одинаково распределённые
по закону R(0,1) случайные
величины, то
и
независимы для любого
и следовательно :
(7)
Пусть
--
оценка
по выборке
,
полученной в результате
-
кратного обращения к исследуемому
датчику :
где 1<t<<
--
выборочное среднее. Заметим , что
-
выборочная дисперсия).
Тест «ковариация» позволяет проверить свойство (7) (гипотезу ) для последовательности и описывается следующим решающим правилом :
принимается
(8)
где :
для
Δ - порог, определённый для заданного
уровня значимости
по формуле :
Δ = Ф-1(1 - ).
В пакете СТАТМОД предполагается использование эквивалентной формы правила (8) :
принимается
где
Представим (8) в виде:
принимается
(9)
Здесь
- верхняя и нижняя доверительные
границы для
,
соответствующие доверительной вероятности
Решающее правило в форме (9) удобно использовать при визуальном анализе графиков ковариационной функции и её оценки .
Тест «равномерность двухмерного распределения»
Из выборочных
значений
полученных в результате
-
кратного обращения к датчику БСВ построим
([
]-целая
часть числа
)
векторов :
(10)
Единичный
квадрат
с
центром в точке
разобьём на
ячеек
:
где
-- произвольные вещественные числа.
Вычислим
частоты
попадания
точек с координатами (10) в
ячеек
гиперкуба.
Если
-- независимые, одинаково распределённые
по закону R(0,1) случайные
величины, то теоретические вероятности
попадания
точек с координатами
в ячейки
равны
площадям этих ячеек :
Описываемый тест используется для
проверки гипотезы Н0 о равномерности
двухмерного распределения векторов
и представляет собой следующее решающее
правило:
принимается
(11)
где в случае истинной гипотезы H0
и
статистика
(12)
имеет x2 – распределение с k
– 1 степенями свободы, а порог Δ
определяется как квантиль этого
распределения:
,
где ε – заданный уровень значимости.
В пакете СТАТМОД предполагается использование эквивалентной формы правила (11):
принимается
где Р - значение вычисляется по формуле: P=1-F(x2), здесь F(.)- функция распределения статистики (12).
1.1.2.Статический анализ точности моделирования
Тесты проверки точности моделирования
БСВ допускают графическую интерпретацию,
удобную для быстрого визуального анализа
качества смоделированной выборки
реализаций БСВ.
В пакете СТАТМОД реализованы следующие графики:
диаграмма рассеяния (иллюстрирует зависимость между ai и at-1
;
для проверки факта наличия зависимости
используется тест «ковариация»);
гистограмма и плотность распределения (позволяет осуществить сравнительный анализ теоретического и эмпирического распределения выборки А; для проверки согласия распределения А с равномерным законом используются критерии согласия);
график корреляционной функции с указанием доверительных границ (служит для графической поддержки теста «ковариация»);
визуализация выборки (иллюстрирует зависимость ai от t
и может использоваться для графической
поддержки критериев серий).
