Скачиваний:
53
Добавлен:
01.05.2014
Размер:
93.18 Кб
Скачать

Эргодическая теорема для цепей Маркова

] известно (???) 

a) выполняется это уравнение

б) Всегда выполняется 1 или 2:

  1. все j=0

  2. j=1

  3. у цепи имеется единственное стационарное распределение  тогда существует положительный 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, где tR+  (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))/tp(t)p’(0) при t0

p’(0) – диф-л в нуле значит, что: если ф-я дифф. pii(е)=(значит в нуле)1+iit+0(t)

1 – диагональный элемент матрицы

i<>j pij(t)=ij=0

Выведем диф-ое ур-ие с помощью которого pij можно находить

Введём матр. Производную

=(ij) – сумма эл-ов по строкам = 0

Теорема: В этих условиях имеет место системы уравнений – системы Колмогорова

  1. Обратная систему уравнений Колмогорова

pij’(t)=

  1. Прямая СУ Колмогорова

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)

  1. pij’(t)=-pii(t)

pi+1,j не будет т.к. j>=ш, а так от i+1 до I (а н ас ??? возрастает)

Решим эти уравнения

(1) ln pii(t)=-t+C pii(t)=e-t+C1 t=0C1=1

pii(t)=e-t

  1. 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

Соседние файлы в папке Марковские цепи