
- •Аннотация
- •1.Основные понятия теории моделирования
- •1.2 Определение абстрактной системы по р. Калману
- •Где –упорядоченное по возрастанию множество моментов времени,
- •1.3 Классификация систем по типу поведения
- •1.4 Принципы моделирования
- •2. Основы технологии имитационного моделирования
- •1.2 Аналитическое и имитационное моделирование
- •2.2 Понятие статистического эксперимента
- •2 .3 Область применения имитационных моделей
- •2.4 Управление модельным временем
- •2.5 Оценка качества имитационной модели
- •3. Типовые математические схемы моделей систем
- •3.1.2. Системы разностных уравнений.
- •3.1.3. Примеры систем, описываемых в классе d – схем
- •3.2 Дискретно-детерминированная модель. Конечные автоматы
- •3.3 Дискретно-стохастическая модель. Вероятностный автомат
- •3.4 Марковский процесс с дискретным временем перехода. Уравнения для вероятностей состояний.
- •3.5 Непрерывно-стохастическая модель. Потоки случайных событий.
- •3.6. Марковские процессы с непрерывным временем. Дифференциальные уравнения относительно вероятностей состояний.
- •3.5. Процессы в системах массового обслуживания как непрерывные
- •3.6, Непрерывно-стохастическая модель. Системы массового обслуживания.
- •3.7 Обобщенные модели
- •4. Описание программных средств визуального моделирования системы matlab
- •4.1 Некоторые сведения о пакете matlab
- •4.2 Краткая характеристика расширения matlab – simulink.
- •4.3 Описание программного средства stateflow.
- •5. Разработка моделей средствами simulink-stateflow
- •5.1 Вероятностный автомат
- •5.2 Моделирование процессов функционирования систем массового обслуживания как марковских процессов.
- •5.2.1 Смо замкнутого типа
- •5.2.2 Смо циклического типа с ветвлением переходов
- •5.2.3 Смо смешанного типа
- •5.2.4 Смо смешанного типа с различными характеристиками каналов
- •5.3 Общая схема моделирования смо.
- •Заключение
3.4 Марковский процесс с дискретным временем перехода. Уравнения для вероятностей состояний.
Пусть
система имеет конечное или счетное
пространство состояний
и
вероятность
того, что если система в момент времени
была в состоянии
,
то в момент времени
она окажется в состоянии
(переходная вероятность);
принимает значения на множестве
.
X(t)
называют
марковским процессом, если вероятность
перехода зависит только от
и. не зависит от того, как этот процесс
проходил в прошлом. Будущее зависит от
прошлого только через настоящее. Если
не зависит от
,
то процесс является стационарным.
Марковский процесс называется эргодическим, если по истечении достаточно большого промежутка времени вероятности состояний системы становятся независимыми от начального состояния системы.
Чтобы найти условие эргодичности, вводится понятие транзитивности
Процесс называется транзитивным, если для любой пары состояний
Это
условие не означает, что
и
имеют непосредственную связь, но на
графе можно найти путь из
в
.
Это эквивалентно тому, что нет состояний
без выхода и без входа. Например,
электрическая лампочка, которая может
находиться в состояниях включена,
выключена, перегорела (см. рис.), не
является транзитивной системой.
Требование
транзитивности нужно, чтобы процесс
при t®
µ
не «застрял» в одном состоянии или
группе состояний. Теорема Маркова
утверждает, что транзитивный стационарный
марковский процесс с конечным числомсостояний являетсяэргодическим.Режим или решение,
прикоторомназывается
стационарным или установившимся.
Для марковского процесса вероятность находиться в данном состоянии в данный момент времени (в данном такте) можно получить из сведений о предшествующем состоянии. Действительно, пусть время изменяется по тактам или по шагам с постоянным интервалом и пусть:
–шаг
приращения времени;
–вероятность
того, что в момент
система находилась в состоянии
,
=1,2,…,
;
–вероятность
перехода из состояния
в состояние
в момент времени
.
Тогда, в соответствии с формулой полной вероятности, будем иметь
(1)
(здесь
выступают как вероятностиn
взаимоисключающих гипотез или
предположений о состоянии на текущий
момент времени
,
а
-
условные вероятности получить состояние
при этих гипотезах в следующий момент).
Заметим,
что сумма вероятностей
j=1,2..n
как сумма вероятностей несовместных
событий, составляющих полную группу
(теорема сложения вероятностей), равна
1, т.е.
.
(2)
Перепишем (1) с учетом (2) в виде:
или:
(3)
.
(4)
Это по существу система разностных уравнений, которая может быть записана в виде:
(5)
или, если система стационарна, в виде:
,
(6)
где
–
номер такта в заданном интервале (в
частности
,
).
Изменение
вероятностей
ищется решением системы разностных
уравнений (4), (5), (6) при заданных начальных
условиях
при выполнении нормирующего условия:
для
любых
в заданном
интервале.
(7)
Обозначив
,
(8)
уравнение (6) можно записать в виде:
.
(9)
Таким образом, если задана матрица переходных вероятностей и начальное состояние процесса, решение разностной системы (6) позволяет найти численное решение поставленной задачи – задачи построения процесса изменения вероятностей нахождения системы в тех или иных состояниях в зависимости от номера такта.
И в этом смысле система (6) представляет собой аналитическую модель для определения основных характеристик процесса перехода системы из одного дискретного состояния в другое по времени (или по тактам) при заданной матрице переходной вероятности и начальному состоянию.
Если
вероятности
не зависят от номера шага, марковский
процесс является стационарным (или
однородным). В этом случае приk
система (6) может иметь установившееся
решение, которое определяется из условия
,
или из условия
.
Установившиеся значения искомых
вероятностей могут быть найдены из
решения системы алгебраических уравнений,
полученных из (9) приравниванием левой
части к нулю.
Рассмотрим пример простой марковской цепи, заданной таблицей 9 вероятностей переходов и таблицей 10 выходов. Верхняя строка и левый столбец таблицы 9 соответствуют начальному состоянию. Граф переходов представлен на Рис.3.
0 |
0.5 |
0 |
0 |
0.5 |
0 |
0 |
0 |
1 |
0 |
0 |
0 |
0.75 |
0 |
0.25 |
0 |
0 |
0.4 |
0 |
0.6 |
0 |
1.0 |
0 |
0 |
0 |

Рис.3. Граф переходов состояний .
таблица 10
X |
|
|
|
|
|
Y |
0 |
0 |
1 |
1 |
0 |
В
соответствии с (9) можно подсчитать
финальные вероятности нахождения
системы в одном из 4-х состояний (т.е.
вероятности нахождения в определенном
состоянии по прошествии длительного
времени или так называемые финальные
вероятности). Пусть эти вероятности
.
Соответственно графу переходов получим:
аналогично:
и условие нормировки:
;
В
результате получим
.
Естественно, аппарат марковских цепей находит широкое применение при исследовании характеристик реальных систем, моделируемых в классе P-схем.
Схемы вероятностных автоматов (P-схемы) широко используются для разработки методов проектирования систем с дискретными состояниями, проявляющих статистически закономерное случайное поведение. Их также можно использовать в качестве генераторов марковских последовательностей случайных событий, которые необходимы для построения моделей процессов функционирования реальных систем и воздействий внешней среды.