Дискретный Марковский процесс с дискретным временем.

advertisement
Дискретный Марковский процесс с дискретным временем.
Марковский случайный дискретный процесс, протекающий в системе S, характеризуется
состояниями и моментами времени, в которые происходит переход системы из одного состояния в другое. Такие моменты времени могут быть заранее известными или случайными.
Определение 1. Случайный процесс, протекающий в некоторой системе S, называется
процессом с дискретным временем, если переходы системы системы из одного состояния в
другое происходит в заранее известные моменты времени t0, t1, …, tk, которые называют шагами или этапами процесса. В промежутки времени между смежными шагами состояние системы не изменяется.
Определение 2. Случайный процесс, протекающий в некоторой системе S, называется
процессом с непрерывным временем, если переходы системы из одного состояния в другое
возможны в любые, заранее неизвестное, случайные моменты времени.
Рассмотрим Марковский случайный процесс с дискретным временем.
Пусть S1, S2, …, Sk – возможные состояния системы S. Переходы системы из одного состояния в другое может происходить только в моменты времени t0, t1, …, tk. В момент времени t [tk, tk+1), k=1, 2, 3…, система находится в состоянии S(t)=S(tk) и процесс можно рассматривать как случайную функцию шагов tk или номеров шагов.
Обозначим Si(k) (i=1, 2, …n, k= 1, 2, …) событие, состоящее в том, что система S с k-го
шага до (k+1)-го находится в состоянии Si. Тогда случайный процесс с дискретным временем
можно представлять случайной последовательностью (по индексу k) случайных событий
Si(k), i=1, 2, …,n; k=1, 2 …, называемой цепью.
Определение 3. Случайная последовательность называется Марковской цепью, если для
каждого шага вероятность перехода из любого состояния Si в любое состояние Sj не зависит
от того когда и как система S оказалась в состоянии Si.
В любой момент времени t, система S может пребывать только в одном из состояний S1,
S2, …, Sn, то при каждом k= 1, 2,… события S1(k), S2(k), …, Sn(k) несовместны и образуют
полную группу событий.
Реализация дискретного случайного процесса с дискретным временем: за любой конечный промежуток времени можно представить неслучайной конечной последовательностью
(цепью) по индексу k событий Si(k), (i=1, 2, …,n; k=1, 2 …)
Основными
характеристиками
Марковских
цепей
является
вероятности
Pi(k) = P(Si(k)), (i=1, 2, …,n; k=1, 2 …) событий Si(k).
Определение 4. Вероятности Рi(k), (i=1, 2, …,n; k=1, 2 …) называются вероятностями состояний.
Для вычисления вероятностей состояния Рi(k) используются переходные вероятности, которые определяются следующим образом.
Определение 5. Переходной вероятностью Pij(k) из i-го состояния в j-e состояние для kго шага (i=1, 2, …,n; j= 1, 2, …,n ; k=1, 2 …) называют вероятность непосредственного перехода системы S в момент времени tk из состояния Si в состояние Sj.
Если i=j, то переходная вероятность Pij(k)=Pii(k) называют вероятностью задержки системы S в состоянии Si.
Определение 6. Если переходные вероятности не зависят от шагов k, то Марковская цепь
называется однородной, Pij(k)=Pij.
Переходные вероятности для однородной Марковской цепи можно представить в виде
квадратной матрицы n-го порядка.
P  Pij i , j1
n
 P11

P
  21
...

P
 n1
P21
P22
...
Pn 2
... P1n 

... P2 n 
... ... 

... Pnn 
Элементы матрицы Р не зависят от номера шага k и ее порядок определяется числом возможных состояний системы S, на главной диагонали находятся вероятности задержки состояния системы S.
Переходная вероятность как условная вероятность
Pij  PS j (k ) / Si (k  1) 
События S1(k), …, Sn(k) несовместны, образуют полную группу событий, и,
n
P
j1
ij
 1 i  1, 2, ..., n ,
что означает равенство 1 суммы элементов каждой из строк матрицы Р.
Определение 7. Матрица, каждый элемент которой неотрицателен и сумма элементов
каждой строки равен 1, называется стохастической.
Если стохастическая матрица обладает тем свойством, что сумма элементов каждого из
столбцов равна 1, то такая матрица называется двоякостохастической.
Определение 8. Граф состояний системы S с указанием переходных вероятностей называется размеченным.
Пример размеченного графа состояний системы.
P12
S2
S1
P23
P31
P24
P42
S4
S3
P43
Вероятности переходов равные 0 и вероятности задержек на размеченном графе не обозначаются.
В силу свойств стохастической матрицы Р вероятности задержки можно вычислить по
формуле:
n
Pii  1   Pij ,
i  1,..., n
j1
j i
Определение 9. Вектор- строка вероятностей состояния (P1(0), P1(0), …, Pn(0)) в начальный момент времени t=0, непосредственно предшествующий первому шагу, называется вектором начального распределения вероятностей.
Свойства:
1.
P1(0)+P2(0)+…+Pn(0)=1
2.
Если в начальный момент времени t=0 система S находилась в состоянии Sm, то Pm(0)=1 и начальное распределение вероятностей P1(0)=0, P2(0)=0, …,
Pm(0)=1, Pm+1(0)=0,…, Pn(0)=0.
3.
Если известно начальное распределение вероятностей и матрица переходных вероятностей, то можно вычислить вероятности состояния системы для
любого k-го шага.
Теорема 1.
Для однородной Марковской цепи вектор- строка вероятностей состояния от k-го до
(k+1)-го шага равна произведению вектор- строки вероятностей состояния от (k-1)-го до k-го
шага на матрицу переходных вероятностей.
(1)
P1 (k), P2 (k),..., Pn (k)  P1 (k  1), P2 (k  1),..., Pn (k  1)  P
Доказательство
Произведение вектор- строки вероятностей состояния размерности (1n) на матрицу размерности (nn) определяет вектор строку размерности (1n) вероятностей состояний системы на k-м шаге.
Для каждого шага k=1, 2, … рассмотрим n гипотез Hi(k-1), i= 1, 2, …,n состоящих в том,
что от (k-1)-го шага до k-го система S находилась в состоянии Si. Эти гипотезы для каждого
шага несовместны и образуют полную группу событий, с вероятностями состояния P(Hi(k1))=Pi(k-1).
Условная вероятность P(Sj(k)/Si(k-1)) есть переходная вероятность Pij и по формуле полной вероятности получим:
Pj (k )   PH i (k  1)   PS j (k ) / Si (k  1)    Pi (k  1)  Pij
n
n
i 1
i 1
j  1, 2, ..., n
что и доказывает справедливость формулы (2.1), которая является рекуррентной и позволяет вычислять вероятности состояния системы для любого k-го шага, и для однородной
Марковской цепи справедлива формула
(2)
P1 (k), P2 (k),..., Pn (k)  P1 (0), P2 (0),..., Pn (0)  P k .....k  1,2,...
Пример дискретного Марковского процесса с дискретным временем для экономической
системы.
Состояние коммерческого банка характеризуется одной из процентных ставок 2%, 3%,
4% которые устанавливаются в начале каждого квартала и могут изменяться только в начале
следующего квартала.
Состояние системы S1- процентная ставка 2%, S2- ставка 3%, S3-ставка 4%.
Анализ работы в предыдущие годы показал, что изменение переходных вероятностей с
течением времени пренебрежимо мало. Определить вероятности состояний банка на конец
текущего года, если в конце прошедшего процентная ставка составляла 3%, размеченный
граф состояний имел вид:
0,4
S1
0,3
0,2
S2
0,3
S3
0,2
0,1
В данном случае смоделировать поведение банка на конец года можно в виде однородного Марковского дискретного случайного процесса с дискретным временем, или в виде однородной Марковской цепи.
По размеченному графу состояний составим матрицу переходных вероятностей.
 0.4 0.4 0.2 
P   0.2 0.5 0.3 


 0.1 0.3 0.6 
 0.207 0.404 0.389 
P   0.202 0.402 0.396 


 0.195 0.396 0.409 
4
(P1(4), P2(4), P3(4))= (0, 1, 0)*P4= (0.2020, 0.4015, 0.3965)
P1(4)=0.2020,
P2(4)=0.4015,
P3(4)=0.3965
Наиболее вероятное состояние процентной ставки на конец года P2(4)=0.4015 т.е. ставка
3%.
Download