- •Упражнения
- •Классификация состояний
- •Упражнения
- •Классификация состояний в неразложимой цепи Маркова
- •Классификация по асимптотическим свойствам вероятностей Рjj(n)
- •Классификация по арифметическим свойствам вероятностей Рjj(n)
- •Упражнения
- •Симметричное случайное блуждание
- •Упражнения
- •28. Шары.
- •29. Прямая.
- •30. Плоскость.
- •32. Частицы.
- •Теорема солидарности
- •Упражнения
- •Эргодичность и стационарные распределения
- •Алгебраический метод нахождения стационарного распределения
- •Упражнения
- •Дополнение (без доказательства)
- •Литература
- •Оглавление
- •Часть I 1
Упражнения
33. Доказать, что в неразложимой конечной цепи Маркова не может быть нулевых состояний.
34. Доказать, что в неразложимой конечной цепи Маркова не может быть невозвратных состояний.
35. Могут ли все состояния цепи Маркова со счетным числом состояний быть невозвратными?
36. Доказать, что для конечной цепи Маркова состояние возвратно тогда и только тогда, когда оно существенно.
37. Доказать, что, если в цепи Маркова состояние εj возвратно, то система с вероятностью 1 за бесконечное число шагов бесконечно много раз возвратится в εj. Если состояние εj невозвратно, то система с вероятностью 1 за бесконечное число шагов лишь конечное число раз возвратится в εj. Другими словами, после некоторого конечного числа шагов система никогда больше не возвратится в εj.
38. Пусть ξt – номер состояния в цепи Маркова в момент времени t; P(ξ0 = 1) = 1,
матрица
переходных вероятностей имеет вид
.
Положим
Показать, что последовательность ηt является цепью Маркова, найти ее матрицу вероятностей перехода.
39 – 41. Пусть ξ1, ξ2,… – независимые случайные величины, P(ξt = 1) = 1 – P(ξt = –1) = p.
Является ли цепью Маркова последовательность ηt, если :
39. ηt = ξt∙ ξt+1.
40. ηt = ξ1∙ ξ2∙… ∙ξt.
41. ηt = φ(ξt,ξt+1), где φ(–1, –1) = 1, φ(–1, 1) = 2, φ(1, –1) = 3, φ(1, 1) = 4.
Эргодичность и стационарные распределения
Определение
11.
Состояние
εj,
для которого существует
,
не зависящий от i, называется эргодическим.
Замечание. Позднее будет показано, что не зависящий от i существует для любого возвратного состояния εj, не являющегося ни нулевым, ни периодическим. Таким образом, термин «эргодическое состояние» – это синоним для «возвратное, не нулевое, не периодическое состояние».
Смысл эргодичности. Существуют вероятности попадания системы в состояние εj через большой промежуток времени, причем они не зависят от начального состояния системы.
Определение
12. Цепь
Маркова называется эргодической,
если для любых i,
j
существует
(все
состояния эргодические).
Определение
13.
Распределение вероятностей {ak}
называется
стационарным распределением цепи
Маркова, если для любого n
справедливо:
(n).
Пример
4.
Пусть дана однородная цепь Маркова с
матрицей переходных вероятностей P=
Стационарным
распределением для него является (a1,
a2)
= (2/3, 1/3). Проверим, что
(1).
(2/3,
1/3) ∙
Это верно и для n = 2, 3, 4:
P2
=
(2/3,
1/3) ∙
=
P3
=
(2/3,
1/3) ∙
=
P4
=
(2/3,
1/3) ∙
=
Упражнение.
Найдите общий вид матрицы P(n)
= Pn
=
и
убедитесь в стационарности распределения.
Найдите пределы элементов матрицы Pn
и убедитесь, что
Замечание. Если для однородной цепи Маркова справедливо: , то {ak} – стационарное распределение.
Пример 5. Пусть дана однородная эргодическая цепь Маркова с матрицей переходных вероятностей
P=
Найти стационарное распределение.
Решение. Ищем решение уравнения (a1, a2)∙P = (a1, a2).
(a1,
a2)∙
(a1,
a2),
1/3a1 + 1/4 a2 = a1, 2/3a1 + 3/4 a2 = a2, а также a1 + a2 = 1.
Решая систему, находим: a1 = 3/11, a2 = 8/11.
Выясним общие условия существования стационарного распределения.
Теорема 5. Неразложимая непериодическая цепь Маркова принадлежит одному из следующих классов:
а)
или все состояния невозвратные, или
все нулевые. В этом случае
и не существует стационарного
распределения;
б)
или все состояния эргодические, т.е.,
.
В этом случае {Uk}–
стационарное распределение, и не
существует никаких других стационарных
распределений.
Доказательство:
Рассмотрим случай б). Заметим, что
поэтому «обрезанная» сумма не превышает
1, и это верно и в пределе. Так как
,
N
получаем, что U1+…+UN
≤ 1.
Далее,
для любых j,
k,
m,
n
справедливо:
.
Положим n
= 1:
.
Если мы сужаем множество индексов в
правой части, то справедливо неравенство:
,
где А = {1, 2,…, N}.
Пусть теперь m
→ ∞. По теореме о мажорируемой сходимости
можно переходить к пределу. Поскольку
,
,
то
,
что верно для любого А, то есть, для
любого N,
и в пределе при N→
∞ получим
.
Докажем, что имеет место равенство.
Если
для каких–то k
,
то, просуммировав, получим:
,
однако они равны, так как состоят из
одних и тех же слагаемых. Т.о.,
Следовательно, это стационарное распределение.
Докажем
единственность. Пусть {Vk}
удовлетворяет (*),
.
Умножим на pkj
= pkj(1).:
,
просуммируем по k.
.
Но
=Vj,
а
.
Получили:
.
Опять умножим обе части на pjk=
pjk(1),
,
просуммируем по j.
,
(3)
и т.д. по индукции. Пусть
покажем, что
,
просуммируем по k.
Заменим индекс j
на k,
получим:
Пусть
n
→ ∞, с учетом того, что
в пределе получаем Vk
= (V1+…)Uk
= Uk
(так как сумма вероятностей стационарного
распределения равна единице).
Таким образом, стационарное распределение {Vk} совпадает с {Uk}, то есть, стационарное распределение единственно.
В
случае а), если все состояния или
невозвратные, или нулевые, предположим,
что распределение {Vk}
стационарно. Тогда, как показано выше,
имеет место уравнение
.
Но
противоречие. Поэтому стационарное
распределение для неразложимой цепи
может существовать только в эргодическом
случае.
