Практическая работа № 4. Моделирование случайных процессов в экономических системах с дискретными состояниями и непрерывным временем
Задания практической работы рассчитаны на 2 часа аудиторных занятий. В конце работы приведены индивидуальные задания для самостоятельной работы, выполняемые по указанию преподавателя, рассчитанные на 2 часа самостоятельной работы. Практическая работа опирается на материал лекции 5.
Оглавление.
Цель работы
-
Имитация простейшего потока событий с заданной интенсивностью;
-
моделирование случайных величин с заданным законом распределения;
-
имитация функционирования экономической системы с простейшими потоками событий.
Перед началом работы
Создайте на своем сетевом диске в папке ИМЭП/Практика папку Занятие4.
По результатам занятия в папке будут сохранены два проекта.
Модель системы «Гараж» с непрерывным потоками отказов и восстановлений.
Постановка задачи
Постановка задачи выполнена на лекции 5.
Построить модель, реализующую описанный процесс.
Определить финальные вероятности установившегося режима функционирования системы.
Указанную задачу можно реализовать двумя способами:
-
с помощью системы дифференциальных уравнений Колмогорова;
-
с помощью имитации простейших потоков.
Цель моделирования
Построить модель процесса двумя способами, сравнить результаты имитационной и теоретической моделей и сделать вывод о правильности подхода к имитации простейшего потока.
Задание 4.1. ДУ Колмогорова
Разработать модель системы «Гараж», основанную на ДУ Колмогорова.
Технология работы
-
Создайте новый проект и сохраните его с именем ГаражНепрДУ. Выберите вид модели – Непрерывный элементарный объект.
-
Опишите параметры класса – интенсивности потоков отказов и восстановлений (рис. 1).
-
Опишите внутренние переменные класса – вероятности быть исправной или неисправной (рис. 1).
-
Присвойте этим переменным некоторые начальные значения. В дальнейшем вы можете провести эксперименты с моделью и сделать вывод, как влияют начальные вероятности на финальные вероятности состояний.
Примечание. Этим переменным при описании присвоены некоторые значения (на рисунке 0.7 и 0.3) – это начальные значения процесса. Можно было бы задать и другие начальные значения, например 1 и 0. Это означает, что в начальный момент все машины были исправны.
-
Опишите дифференциальные уравнения Колмогорова (рис. 1).
-
Запустите модель.
-
Проведите эксперименты и определите:
-
финальные вероятности процесса;
-
влияние начальных вероятностей состояний на финальные вероятности.
Рис. 1
Задание 4.2. Имитация простейших потоков
Создать имитационную модель системы «Гараж», в которой функционируют два простейших потока – поток отказов и поток восстановлений. Найти финальные вероятности. Сравнить результаты с моделью, построенной в задании 4.1.
Технология работы
-
Создайте новый проект и сохраните его с именем ГаражНепрИмит. Выберите вид модели – Элементарный гибридный объект.
-
Опишите параметры класса:
-
интенсивности потоков отказов и восстановлений;
-
начальную вероятность быть исправным (рис. 2).
-
Опишите внутренние переменные класса:
-
Индикатор нахождения в состоянии k (может принимать значения 1 или 2);
-
Время нахождения в состоянии (t 1, t 2);
-
Переменную для равномерно распределенного случайного числа (ksi).
Рис. 2
-
Создайте карту поведения и в ней опишите состояния Начальник, Исправен, Неисправен (рис. 3).
-
Опишите входные действия для состояния Начало (рис. 3).– генерируется случайное равномерно распределенное число ksi для определения направления перехода по начальной вероятности либо в состояние Исправен, либо в состояние Неисправен.
Рис. 3
6. Опишите входные действия для состояния Исправен (рис. 4):
-
индикатору нахождения в состояние k присвоить значение 1;
-
сгенерировать случайное равномерно распределенное число ksi для вычисления случайной величины, распределенной по показательному закону;
-
вычислить случайный интервал времени t 1 перехода в состояние Неисправен по формуле, полученной на лекции 5:
.
Рис. 4
-
Аналогично опишите входные действия для состояния Неисправен.
-
Опишите условия перехода из состояния Неисправен в состояние Исправен.
-
Скомпилируйте модель.
-
Запустите модель. Проследите по временной диаграмме переход из состояния в состояние (рис. 5).
Рис. 5
-
Включите в отчет описание модели и диаграмму процесса.
Задание 4.3. Накопление статистики и вычисление финальных вероятностей
При моделировании цепи Маркова (пр. раб. 3) мы вычисляли вероятность нахождения в том или ином состоянии как отношение количества переходов в данное состояние к общему числу переходов, т. к. эти переходы происходили через одинаковые промежутки времени. В непрерывной цепи Маркова так поступать нельзя, т. к. длительность нахождения в каждом состоянии является случайным числом. Поэтому используется следующий алгоритм вычисления вероятности нахождения в состояниях:
-
Вычисляется суммарное время пребывания в состоянии Исправен.
-
Вычисляется суммарное время пребывания в состоянии Неисправен.
-
Вычисляется суммарное время процесса.
-
Вероятность нахождения в состоянии вычисляется как отношение времени нахождения в состоянии к общей длительности процесса.
Технология работы
-
Добавьте внутренние переменные (рис. 6).
Рис. 6
2. Опишите входные действия в состоянии (рис. 7).
Рис. 7
-
Скомпилируйте модель и запустите. Постройте временную диаграмму для наблюдения за изменением вероятностей.
-
Проведите эксперименты и определите:
-
финальные вероятности процесса;
-
влияние начальных вероятностей состояний на финальные вероятности.
Рис. 8. Переходный режим
Рис. 9. Установившийся режим
-
Проведите эксперименты и сравните результаты расчетов по теоретической и имитационной моделям. Включите в отчет скриншоты построенных моделей, диаграмму процесса в переходный период и диаграмму установившегося процесса (рис. 8, 9), таблицу сравнительных результатов.
Вид модели
|
Через сколько отсчетов времени система переходит в стационарный режим?
|
Результаты расчета вероятностей состояний
|
P1
Исправен
|
P2
Ремонт
|
ГаражНепрДУ
|
|
|
|
ГаражНепрИмитация
|
|
|
|
Задание 4.4. Для самостоятельной работы
Граф экономической системы имеет вид, как показано на рисунке. Заданы интенсивности переходов из состояния в состояние. Записать дифференциальные уравнения Колмогорова. Найти финальные вероятности теоретически.
Построить две модели процесса функционирования системы:
Модель 1 – путем интегрирования системы ДУ;
Модель 2 – путем составления имитационной модели.
Определить финальные вероятности на моделях.
Составить отчет о проделанной работе в электронном виде. Включить в отчет:
-
исходный граф системы;
-
исходные данные;
-
систему ДифУр-й;
-
систему для определения финальных вероятностей теоретически;
-
скриншоты описания класса, карты поведения, временную диаграмму установившегося процесса;
-
результаты моделирования в виде таблицы:
Состояние
|
Вероятность нахождения в состоянии Si
|
Модель 1
|
Модель 2
|
Время, начиная с которого можно считать вероятность установившейся
|
Значение
|
Количество тактов имитации для достижения установившегося процесса
|
Значение
|
S1
|
|
|
|
|
S2
|
|
|
…
|
|
|
Вариант 1
Вариант 2
Вариант 3
Вариант 4
Вариант 5

Достарыңызбен бөлісу: |