Марковские цепи / TVIMS
.DOCЭргодическая теорема для цепей Маркова
] известно (???)
a) выполняется это уравнение
б) Всегда выполняется 1 или 2:
-
все j=0
-
j=1
-
у цепи имеется единственное стационарное распределение тогда существует положительный y не единственный возвратный класс состояний(pj)
Доказательство
-вектор
означает что p=p т.е. (I-p)=0 – условие для стационарного распределения
j-стационарное распределение т.е. цепь стремится к своему старшему распределению.
Если какое-то является опис. Распределением м/у (то легко показать qj(m)j при n) т.к. qj(n)= qpj(n)j q=j)
то различные отрезка времени от 0 до N – (достаточно большой) т.е. до выхода на стационарный режим.
Лемма Фату:
] существуют положительные числа
- ряд сходится
рассмотрим нижний предел
а) (1) по лемме Фату
(1) ур. Ч.К.
Следовательно
Хотим доказать равенство что меньше не может быть
] При j0 < возьмем и сложим эти неравенства по всем j это <
поменяем порядок суммирования i и j местами
одно и тоже
Получили противоречие
Следовательно имеет место равенство.
б) j=по уравнению из а) p=
pn=
следовательно
для любого j j(1-I)=0
] существует м/у Xn и есть подмножество Xnk, n1<n2 … <nk
Xnk – ж/у (но она не однородная)
Напишем Xt, где tR+ (0,)
И этот процесс можно было бы считать м/у.
Он является: м/у с непрерывным временем, если при любом выборе t1…tk…последовательность X+k – это м/у. И будем считать что она однородная м/у, если Xt*k является однородным м/у для всякого t (сколь угодно маленького) Установим соответствие между м/у (с непрерывным вектором и нет).
Вектор ??? Q=(q1,q2,…) можно pin(t)=для всех интервалов времени длинны t задать.
Переходим En-1, En+1, En.
Матрица перех. Вероятностей за 1 шаг у нас нет т.к. t – мало. Можно определить p(k)
Pij(t)=p(Xt+=j/X=i) (pk у нас нет)
<t+ Условная вероятность составляет E в момент t+ при условии что в момент это сост. EI
т.к. мы не можем возводить матрицу в нецелую степень)
(зависит от продолжения t)
Уравнение Ч.К. выполняется
(ничего не можем сказать) и
тоже будет выполнятся как
Конечномерное распред. Мн. Можно определить м/у (не однородная)
P(xt1=i1,X2=i2,…,Xtk=ik)=pik-1,ik(tk-1,tk)*pik-2,ik-1*(tk-2,tk-1)…pi0i1(t0,t1)q+0 pij(t) – могут менятся со временем указ. конца пром. времени.
м/у с непр. вр. однозначно определяется вектором начального распределения и набором матриц
p(t)=(pij(t))
но p(t+)=p(t)*p() – обладают этим свойством t,>0
эти матрицы над полугруппами (если обладают такими св-ми)
p(0)=E
i=j pii(0)=1
i<>j pij(0)=0
? – как определить эту полугруппу матриц? Как найти pij(t)? Какую хар-ку выбрать через которую эта подгруппа определится?
Достаточно знать полугруппу на маленьком промежутке времени и знаем всю гр. Pn(t)=p{Z(t)=n}
Действительно t=nt0+t1, t1<t0
P(t)=p(nt0)*p(t1)=p(t0)*p(t1) Нужно знать значение p(t) близких к нулю будем знать хар-ку Маркова инфинитизимальная х-ка
] p(t) – диф-ал в нуле, т.е. pij(t) в точке 0
Покажем, что она будет диф-ма и в любой другой точке
(p(t+t)-p(t))/t=(p(t)p(t)-p(t))/t=(p(t)(p(t)-E))/t=p(t)(p(t)-p(0))/tp(t)p’(0) при t0
p’(0) – диф-л в нуле значит, что: если ф-я дифф. pii(е)=(значит в нуле)1+iit+0(t)
1 – диагональный элемент матрицы
i<>j pij(t)=ij=0
Выведем диф-ое ур-ие с помощью которого pij можно находить
Введём матр. Производную
=(ij) – сумма эл-ов по строкам = 0
Теорема: В этих условиях имеет место системы уравнений – системы Колмогорова
-
Обратная систему уравнений Колмогорова
pij’(t)=
-
Прямая СУ Колмогорова
pij’=
Начальные условия
Эти системы при i<>j решаются однозначно имеется момент однозначно востановить p(t)
Доказательство:
перейдем к пределу
- обратное уравнение системы Колмогорова
2-й способ представления
- прямая СУ Колмогорова / Что и требовалось доказать
Q и - определяют м/у однозначно (остальные формулы для дискретных цепей используются для м/у непрер. врем.)
Пример:
Расс. произв, р. со временем возрастает т.е. Xt имеет возрастающую траекторию т.е. pij(t)=0, если j<i
] ij= задаём матрицу
Возьмём СУ Колмогорова
Обратная СУК:
(2) (*)
(<>0, j=I, j=I+1 2 слагаемых)
(*) = -pij(t)+pi+1,j i<>j
(Есть один особый случай когда i=j)
-
pij’(t)=-pii(t)
pi+1,j не будет т.к. j>=ш, а так от i+1 до I (а н ас ??? возрастает)
Решим эти уравнения
(1) ln pii(t)=-t+C pii(t)=e-t+C1 t=0C1=1
pii(t)=e-t
-
fi,i+1(t)=etpi,i+j(t) fi,i+1’(t)=et(pi,i+j(t)+pi,i+j’(t))=fi+1,i+j(t)
Сист ема диф ур для f
fi,i+j’(t)=fi+1,i+j(t)
fii(t)=f00(t)=1
Получили fi,i+1’(t)
(чтобы доказать по ММИ некую формулу, для того чтбы найти формулу рассмотрим частный случай)
fi,i+1’=fi+1,i+1(t)=
fi,i+1(t)=t+C
pi,i+1(0)=0 | t=0
fi,i+1(0)=0 |
fi,i+1(t)=t