Курс лекций для специальности Прикладная математика и информатика
Вид материала | Курс лекций |
СодержаниеМарковское свойство |
- Программа по курсу "Математика. Алгебра и геометрия" для специальности 080801 (351400), 143.45kb.
- Программа вступительного экзамена по математике подготовки магистров по направлению, 86.94kb.
- Программа дисциплины дс. 08 «Информационная безопасность» для студентов специальности, 149.66kb.
- Программа дисциплины ф дифференциальные уравнения для студентов специальности 010501, 101.63kb.
- Программа дисциплины математический анализ и обыкновенные дифференциальные уравнения., 139.76kb.
- Рабочая программа по дисциплине «принятие управленческих решений» (по выбору) для специальности, 89.25kb.
- Программа дисциплины Современная прикладная алгебра для направления 010500 Прикладная, 214.78kb.
- Программа дисциплины ф. 8 Общая физика Разделы «Механика», «Колебания и волны», «Молекулярная, 113.79kb.
- «Прикладная математика и информатика», 3781.56kb.
- Основная образовательная программа специальности высшего профессионального образования, 68.9kb.
Пусть
![](images/5855-nomer-6132cc61.gif)
![](images/5855-nomer-46cd7dff.gif)
![](images/5855-nomer-46cd7dff.gif)
![](images/5855-nomer-m5c0a2202.gif)
![](images/5855-nomer-1b2cc633.gif)
Определение. Последовательность случайных величин
![](images/5855-nomer-b4c5d90.gif)
![](images/5855-nomer-m68ef6d0a.gif)
![](images/5855-nomer-13063254.gif)
![](images/5855-nomer-m46c83c88.gif)
![](images/5855-nomer-6452189c.gif)
![](images/5855-nomer-1b2cc633.gif)
![](images/5855-nomer-18149a9d.gif)
Множество
![](images/5855-nomer-46cd7dff.gif)
![](images/5855-nomer-2417dbe2.gif)
![](images/5855-nomer-6fa6833e.gif)
Марковское свойство:
![](images/5855-nomer-372c7819.gif)
Определим
![](images/5855-nomer-m5bcb6d7a.gif)
![](images/5855-nomer-23d507d4.gif)
![](images/5855-nomer-314873ed.gif)
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-2417dbe2.gif)
Если
![](images/5855-nomer-m186b1470.gif)
![](images/5855-nomer-7ee8069a.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-6132cc61.gif)
Чтобы знать все вероятности
![](images/5855-nomer-56e447dd.gif)
.
Например,
![](images/5855-nomer-292d7522.gif)
![](images/5855-nomer-79cc16e9.gif)
Рассмотрим матрицу, которую обозначим
![](images/5855-nomer-m4a31b088.gif)
![](images/5855-nomer-71670e8e.gif)
Для матрицы
![](images/5855-nomer-m7e943538.gif)
![](images/5855-nomer-m265abdda.gif)
Любая матрица с такими совйствами называется стохастической.
Итак, чтобы уметь описывать цепь Маркова, необходимо знать
- Вектор начальных состояний
и матрицу
.
Задачи.
- Пусть
- цепь Маркова. Показать, что
тоже цепь Маркова.
-
независимые одинаково распределенные целочисленные случайные величины с известными вероятностями
Какой вид имеет матрица
Лекция 14 (7.12.10)
Зададимся следующим вопросом: верно ли, что для однородных марковских цепей справедливо равенство
![](images/5855-nomer-m5a85a08f.gif)
при всех натуральных k?
Утверждение:
![](images/5855-nomer-6132cc61.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-m5a85a08f.gif)
![](images/5855-nomer-m46c83c88.gif)
![](images/5855-nomer-mc5f3267.gif)
![](images/5855-nomer-m46c83c88.gif)
![](images/5855-nomer-69e9f430.gif)
![](images/5855-nomer-52432088.gif)
Доказательство: индукцией по
![](images/5855-nomer-2417dbe2.gif)
База:
![](images/5855-nomer-m35ddeea4.gif)
Пусть верно для
![](images/5855-nomer-m4e0b9eea.gif)
![](images/5855-nomer-2417dbe2.gif)
Рассмотрим
![](images/5855-nomer-4bb9b0e8.gif)
=
![](images/5855-nomer-m2af482f4.gif)
![](images/5855-nomer-3f403ba3.gif)
![](images/5855-nomer-m53687454.gif)
![](images/5855-nomer-mf508e86.gif)
![](images/5855-nomer-3e20f3ee.gif)
![](images/5855-nomer-m17e8deac.gif)
![](images/5855-nomer-m2e486d7d.gif)
![](images/5855-nomer-4cf99077.gif)
![](images/5855-nomer-6a7f82af.gif)
Введем обозначение:
![](images/5855-nomer-2696e589.gif)
![](images/5855-nomer-b9337b.gif)
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-314873ed.gif)
![](images/5855-nomer-2417dbe2.gif)
![](images/5855-nomer-md311bfc.gif)
Пусть известна
![](images/5855-nomer-m71880c1b.gif)
Как найти
![](images/5855-nomer-7d28b85e.gif)
![](images/5855-nomer-2417dbe2.gif)
Теорема (Маркова – Чепмена - Колмогорова)
![](images/5855-nomer-m46c83c88.gif)
![](images/5855-nomer-52432088.gif)
![](images/5855-nomer-m513f49de.gif)
![](images/5855-nomer-m12dcab1b.gif)
![](images/5855-nomer-m4b7f71d0.gif)
![](images/5855-nomer-m71880c1b.gif)
![](images/5855-nomer-m71880c1b.gif)
![](images/5855-nomer-m4b7f71d0.gif)
![](images/5855-nomer-m12dcab1b.gif)
Следствия:
=
;
.
Доказательство (теоремы). Достаточно доказать равенство:
![](images/5855-nomer-68ed75eb.gif)
Имеем
![](images/5855-nomer-414c7293.gif)
![](images/5855-nomer-1cb49a85.gif)
![](images/5855-nomer-6d36da8b.gif)
![](images/5855-nomer-3e8c788d.gif)
Доказали следующее равенство:
![](images/5855-nomer-m513f49de.gif)
![](images/5855-nomer-m71880c1b.gif)
![](images/5855-nomer-m4b7f71d0.gif)
![](images/5855-nomer-m1596ee09.gif)
Замечание (к теореме и утверждению): Умножать и делить можно лишь тогда, когда в знаменателе не 0.
Пусть
![](images/5855-nomer-m48413e13.gif)
![](images/5855-nomer-m590ef7c4.gif)
![](images/5855-nomer-6bada1a1.gif)
![](images/5855-nomer-mc6dd90.png)
![](images/5855-nomer-2dd6cdde.gif)
![](images/5855-nomer-m30698612.gif)
Можно и наоборот, по
![](images/5855-nomer-m418a747e.gif)
Задача: Знаем
![](images/5855-nomer-m7054b1b.gif)
![](images/5855-nomer-m418a747e.gif)
§25. Классификация состояний.
Пусть
![](images/5855-nomer-mc5f3267.gif)
Определения.
Состояние
![](images/5855-nomer-314873ed.gif)
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-m61eaac98.gif)
![](images/5855-nomer-7a0044fb.gif)
![](images/5855-nomer-2417dbe2.gif)
![](images/5855-nomer-m138dd33d.gif)
Состояния
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-314873ed.gif)
![](images/5855-nomer-m1459eb7d.gif)
![](images/5855-nomer-7a0044fb.gif)
![](images/5855-nomer-m61eaac98.gif)
![](images/5855-nomer-m60031c6e.gif)
![](images/5855-nomer-7a0044fb.gif)
![](images/5855-nomer-23416f5.gif)
![](images/5855-nomer-69d19978.gif)
![](images/5855-nomer-77ff7e5b.gif)
![](images/5855-nomer-2417dbe2.gif)
![](images/5855-nomer-m532294e.gif)
Состояние
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-52770f7a.gif)
Утверждение:
и
Доказательство: Существует
![](images/5855-nomer-2417dbe2.gif)
![](images/5855-nomer-m138dd33d.gif)
![](images/5855-nomer-m532294e.gif)
![](images/5855-nomer-m204701d5.gif)
Рассмотрим
![](images/5855-nomer-2d6ee77d.gif)
- Пусть
- существенное состояние и
. Тогда
тоже существенное состояние.
Доказательство: Предположим противное: пусть
![](images/5855-nomer-314873ed.gif)
![](images/5855-nomer-7a0044fb.gif)
![](images/5855-nomer-m7af21693.gif)
![](images/5855-nomer-32477add.gif)
![](images/5855-nomer-m119da8c9.gif)
![](images/5855-nomer-m61eaac98.gif)
![](images/5855-nomer-32477add.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-5fbb79d8.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-m1ad1b04b.gif)
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-32477add.gif)
![](images/5855-nomer-m1ad1b04b.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-m60031c6e.gif)
![](images/5855-nomer-m61eaac98.gif)
- Пусть
- существенные состояния и
. Пусть
- множество состояний, сообщающихся с
,
- множество состояний, сообщающихся с
. Тогда
.
Доказательство: следует непосредственно из предыдущих утверждений.
![](images/5855-nomer-m64046af3.png)
![](images/5855-nomer-1be76372.png)
Введем обозначение:
![](images/5855-nomer-m36657394.gif)
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-2417dbe2.gif)
![](images/5855-nomer-46510956.gif)
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-2417dbe2.gif)
![](images/5855-nomer-m1f95291d.gif)
![](images/5855-nomer-bd4efe5.gif)
Определение.
Состояние
![](images/5855-nomer-bd4efe5.gif)
![](images/5855-nomer-7a0044fb.gif)
![](images/5855-nomer-mf765bee.gif)
Если
![](images/5855-nomer-13f2808b.gif)
![](images/5855-nomer-bd4efe5.gif)
Примеры:
k – состояние (все целые числа на вещественной прямой, их счетное количество)
![](images/5855-nomer-m400ddc15.gif)
Здесь возвратных состояний нет.
-
,
,
,
.
Т.о., состояние 3 возвратное.
Теорема (Критерий возвратности)
- возвратное состояние
(ряд расходится);
- невозвратное состояние
.
Доказательство теоремы.
Пусть имеется последовательность чисел
![](images/5855-nomer-m108f427c.gif)
![](images/5855-nomer-m6e145e2.gif)
![](images/5855-nomer-m66897ff9.gif)
![](images/5855-nomer-m244910c3.gif)
![](images/5855-nomer-40dca844.gif)
![](images/5855-nomer-m108f427c.gif)
Рассмотрим
![](images/5855-nomer-m3fd6f5f7.gif)
![](images/5855-nomer-5a870a29.gif)
![](images/5855-nomer-m149640c1.gif)
![](images/5855-nomer-7e59f660.gif)
Имеем по формуле полной вероятности:
![](images/5855-nomer-m52dbf767.gif)
![](images/5855-nomer-59088fd.gif)
![](images/5855-nomer-m573db612.gif)
![](images/5855-nomer-7adb2d00.gif)
![](images/5855-nomer-5e84a5bf.gif)
![](images/5855-nomer-me701580.gif)
![](images/5855-nomer-7adb2d00.gif)
![](images/5855-nomer-76fc346a.gif)
![](images/5855-nomer-m73f7d0dd.gif)
Заметим, что
![](images/5855-nomer-m66b345ca.gif)
Из (*) следует
![](images/5855-nomer-m7ecd82ed.gif)
![](images/5855-nomer-11b9d9f2.gif)
![](images/5855-nomer-bd4efe5.gif)
Если ряд расходится, то при
![](images/5855-nomer-m6c32fba6.gif)
![](images/5855-nomer-m320ae59.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-3fe5d627.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-m66b345ca.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-bd4efe5.gif)
Лекция 15 (14.12.10)
Пусть
![](images/5855-nomer-75422c90.gif)
![](images/5855-nomer-mecffa3d.gif)
![](images/5855-nomer-53175ce5.gif)
Тогда
![](images/5855-nomer-5bfd2521.gif)
![](images/5855-nomer-53175ce5.gif)
![](images/5855-nomer-m17ed43dd.gif)
![](images/5855-nomer-1355fa29.gif)
![](images/5855-nomer-m17ed43dd.gif)
…
Понятно, что эти классы не пересекаются, т.е.
![](images/5855-nomer-m40301f1.gif)
Перенумеруем состояния. Сначала занумеруем состояния из класса S0, затем из класса S1 и т.д.
Тогда матрица переходных вероятностей за n шагов будет иметь следующий вид:
| ![]() | ![]() | ![]() | … |
![]() | | | | |
![]() | 0 | | 0 | 0 |
![]() | 0 | 0 | | 0 |
![]() | | | | |
![](images/5855-nomer-m4fe1be16.gif)
Каждый блок в отдельности – стохастическая матрица.
Определение. Цепь Маркова называется неразложимой, если она состоит из одного класса сообщающихся между собой состояний.
Теорема. (О солидарности для возвратности)
Пусть ЦМ неразложимая
![](images/5855-nomer-300113d7.gif)
Доказательство.
Пусть
![](images/5855-nomer-3493f7b5.gif)
![](images/5855-nomer-566f319f.gif)
Пусть
![](images/5855-nomer-b6f05cc.gif)
Напомним, что по критерию возвратности
![](images/5855-nomer-3493f7b5.gif)
![](images/5855-nomer-7a0044fb.gif)
![](images/5855-nomer-m7c91d287.gif)
Рассмотрим
![](images/5855-nomer-6dc510f8.gif)
![](images/5855-nomer-335a91d.gif)
![](images/5855-nomer-300113d7.gif)
![](images/5855-nomer-m484ddd93.gif)
Случайные блуждания по целым точкам в
![](images/5855-nomer-b78446c.gif)
![](images/5855-nomer-111bf1c7.png)
![](images/5855-nomer-m3906870e.gif)
![](images/5855-nomer-75422c90.gif)
![](images/5855-nomer-d9f1516.gif)
![](images/5855-nomer-m5bb1d4af.gif)
![](images/5855-nomer-m2121923f.gif)
![](images/5855-nomer-22263dcc.gif)
Проверим, будут ли состояния этой ЦМ возвратными. Понятно, что достаточно провести проверку только для одного состояния, т.к. все состояния сообщаются между собой.
Очевидно, что
![](images/5855-nomer-m4be65b61.gif)
Рассмотрим
![](images/5855-nomer-m4aea7fcd.gif)
![](images/5855-nomer-2b871937.gif)
Используем формулу Стирлинга (
![](images/5855-nomer-m3cd04de1.gif)
![](images/5855-nomer-m60d04a4b.gif)
![](images/5855-nomer-5d8fd01a.gif)
![](images/5855-nomer-7539e7d1.gif)
![](images/5855-nomer-2b317973.gif)
![](images/5855-nomer-m3906870e.gif)
![](images/5855-nomer-2a10e20c.gif)
![](images/5855-nomer-60519b69.gif)
![](images/5855-nomer-31798b96.gif)
Итак, если случайное блуждание по целым точкам на вещественной прямой симметрично (т.е. p=q), то оно возвратное и наоборот.
Случайное блуждание по решётке «целочисленных» точек в
![](images/5855-nomer-4ab3bfd3.gif)
![](images/5855-nomer-3e49746c.png)
Симметричные случайные блуждания на плоскости тоже будут возвратными.
Рассмотрим событие, которое означает возврат в начальное состояние через
![](images/5855-nomer-3ff74ac.gif)
Каждая из таких ситуаций может быть записана цепочкой длины
![](images/5855-nomer-3ff74ac.gif)
![](images/5855-nomer-3ff74ac.gif)
Здесь П означает шаг вправо, Л – шаг влево, В – шаг вверх и Н – шаг вниз, причем для возврата в начальное состояние число шагов в каждом направлении должно быть следующим:
![](images/5855-nomer-m4d17405a.gif)
![](images/5855-nomer-3ff74ac.gif)
Тогда
![](images/5855-nomer-m425f082b.gif)
![](images/5855-nomer-1d81d6f8.gif)
![](images/5855-nomer-477475a2.gif)
Ряд с такими членами расходится, и следовательно, такое случайное блуждание тоже возвратное.
Случайное блуждание в
![](images/5855-nomer-m34f60884.gif)
Рассмотрим симметричные случайные блуждания по точкам с целочисленными координатами в трехмерном пространстве. Здесь шаг вверх – ВВ, шаг вниз – ВН, шаг вправо – П , шаг влево – Л, шаг вперед – ВП и шаг назад – Н и пусть
![](images/5855-nomer-3456a75.gif)
Для возврата в начальное состояние число шагов в каждом направлении должно быть
следующим:
![](images/5855-nomer-265e3a08.gif)
![](images/5855-nomer-3ff74ac.gif)
Тогда
![](images/5855-nomer-75000c49.gif)
![](images/5855-nomer-m8fddbab.gif)
![](images/5855-nomer-m630a216a.gif)
Заметим, что
![](images/5855-nomer-m1fbb572d.gif)
![](images/5855-nomer-78645bb2.gif)
Поскольку
![](images/5855-nomer-m5628c7fb.gif)
то начальное состояние невозвратное, а следовательно и все состояния невозвратные.