Простая цепь Маркова полностью характеризуется матрицей переходных вероятностей { pij }.
Будем рассматривать вероятности Рij матрицы переходных вероятностей, как условные вероятности наступления в данном испытании события Aj , при условии, что результатом предыдущего испытания было событие Аi .
Моделирование цепи Маркова в этом случае состоит в последовательном выборе событий Аj в соответствии с заданными переходными вероятностями. Сначала выбирается некоторое начальное состояние цепи Маркова, которое может быть задано априорно, либо выбирается по жребию.
Пусть начальным состоянием будет состояние Аl , которое задается вероятностями pl1, pl2, . . ., plK , образующими l - строку матрицы переходных вероятностей.
Так как сумма вероятностей в l-ой строке тождественно равна единице, то формируется число g , равномерно распределенное на интервале [0, 1] и проверяется в какой из интервалов plj = Di это число попадает:
(1)
Пусть число g попало в интервал, соответствующий m-му состоянию. Это означает, что следующим событием данной реализации цепи Маркова будет событие Аm . Выбирается m-ая строка матрицы {Рij}, вновь формируется число g и проверяется условие (1) , но уже для m-ой строки. Аналогичная процедура повторяется и далее.
Метод генерирования процессов Маркова с дискретным множеством состояний.
Процессы Маркова с дискретным множеством состояний задаются: дискретным множеством состояний и интенсивностями переходов из одного состояния в другое состояние , . Интенсивность выхода из состояния определяется через интенсивности переходов следующим образом
= .
Тогда моделирование процесса Маркова с дискретным множеством состояний осуществляется следующим образом: генерируется по вышеизложенному случайная величина, имеющая экспоненциальное распределение с параметром , а затем осуществляется переход в , , с вероятностями