36

Скачкообразные марковские процессы

[17/71%]
Показать
LaTeX
Задача 36.1

Доказать, что процесс Пуассона с интенсивностью λ\lambda является однородным скачкообразным марковским процессом. Найти переходную функцию этого процесса.

?
Примечание.
?

Пусть случайный процесс ξ(t)\xi (t) принимает не более чем счётное число значений {x1,x2,…}\{ x_{1}, x_{2}, \ldots \} и имеет кусочно-постоянные, непрерывные справа траектории. Пусть σ(t)\sigma (t) — σ\sigma-алгебра, порождённая процессом ξ\xi до момента времени tt включительно. Процесс ξ(t)\xi (t) называется скачкообразным марковским процессом, если для любых состояний xi,xjx_{i}, x_{j}, для любых t,h>0t, h > 0 и для любого события B∈σ(t)B \in \sigma (t) выполняется P(ξ(t+h)=xj∣ξ(t)=xi,B)=P(ξ(t+h)=xj∣ξ(t)=xi)\mathbb {P}\left(\xi (t+h)=x_{j} \mid \xi (t)=x_{i}, B\right) = \mathbb {P}\left(\xi (t+h)=x_{j} \mid \xi (t)=x_{i}\right).

Скачкообразный марковский процесс ξ(t)\xi (t) однороден по времени, если для любых состояний xix_{i} и xjx_{j} найдётся такая переходная функция pij(h)p_{ij}(h), h≥0h \geq 0, что для всех tt выполняется P(ξ(t+h)=xj∣ξ(t)=xi)=pij(h)\mathbb {P}\left(\xi (t+h)=x_{j} \mid \xi (t)=x_{i}\right) = p_{ij}(h), при этом для любых ii и h≥0h \geq 0 имеет место ∑jpij(h)=1\sum_{j} p_{ij}(h) = 1. Известно, что распределение однородного скачкообразного марковского процесса ξ(t)\xi (t), t≥0t \geq 0, однозначно определяется начальным распределением P(ξ(0)=xi)\mathbb {P}\left(\xi (0)=x_{i}\right) и матрицей переходных функций ∥pij(⋅)∥\left\| p_{ij}(\cdot )\right\|.

Отметим, что pii(0)=1p_{ii}(0) = 1 и pij(0)=0p_{ij}(0) = 0 при i≠ji \neq j. Если функции pij(h)p_{ij}(h) дифференцируемы, то P(ξ(t+h)=xj∣ξ(t)=xi)=pij(0)+qijh+o(h)\mathbb {P}\left(\xi (t+h)=x_{j} \mid \xi (t)=x_{i}\right) = p_{ij}(0) + q_{ij}h + o(h) при h→0h \to 0, где qij=pij′(0)q_{ij} = p'_{ij}(0). При i≠ji \neq j величина qij≥0q_{ij} \geq 0 есть интенсивность перехода процесса из состояния xix_{i} в состояние xjx_{j}. Матрица Q=∥qij∥Q = \left\| q_{ij}\right\| называется матрицей переходных интенсивностей (или генератором, или инфинитезимальным оператором); в случае конечного числа состояний для любого ii имеет место равенство ∑jqij=0\sum_{j} q_{ij} = 0.

Задача 36.2

Доказать, что переходные вероятности однородного скачкообразного марковского процесса удовлетворяют уравнению Колмогорова — Чепмена:

pij(t+h)=∑kpik(t)pkj(h). p_{ij}(t+h) = \sum _{k} p_{ik}(t)p_{kj}(h).
?
Задача 36.3

Доказать, что в случае дифференцируемости переходные вероятности однородного скачкообразного марковского процесса с конечным числом состояний удовлетворяют прямым уравнениям Колмогорова

pij′(t)=∑kpik(t)qkj, p'_{ij}(t) = \sum _{k} p_{ik}(t)q_{kj},

а также обратным уравнениям Колмогорова

pij′(t)=∑kqikpkj(t). p'_{ij}(t) = \sum _{k} q_{ik}p_{kj}(t).
?
Задача 36.4

Пусть однородный скачкообразный марковский процесс имеет два состояния 0 и 1, где, к примеру, 0 означает, что некоторый субъект здоров, а 1 — что он болен. Найти переходные вероятности, если матрица переходных интенсивностей равна

Q=(−λλμ−μ) Q = \begin{pmatrix} -\lambda & \lambda \\ \mu & -\mu \end{pmatrix}

при некоторых λ,μ>0\lambda , \mu > 0. Найти предел переходных вероятностей pij(t)p_{ij}(t) при t→∞t \to \infty.

?
Задача 36.5

Пусть ξ(t)\xi (t) — однородный скачкообразный марковский процесс с матрицей интенсивностей перехода Q=∥qij∥Q = \left\| q_{ij}\right\|. Доказать, что если в момент времени tt процесс ξ(⋅)\xi (\cdot ) находится в состоянии xix_{i}, то распределение времени до следующего скачка процесса является экспоненциальным с параметром qiiq_{ii}, независимо от предыстории процесса до момента tt, а в момент следующего скачка процесс переходит в состояние xj≠xix_{j} \neq x_{i} с вероятностью −qij/qii-q_{ij}/q_{ii} независимо от времени скачка и от предыстории процесса до него.

?
Задача 36.6

В условиях предыдущей задачи положим T0=0T_{0} = 0, TnT_{n} — время nn-го скачка процесса ξ(⋅)\xi (\cdot ). Доказать, что при условии ξ(Tn)=xi\xi (T_{n}) = x_{i} величина Tn+1−TnT_{n+1} - T_{n} имеет экспоненциальное распределение с параметром qiiq_{ii} и не зависит от предыстории процесса до момента времени TnT_{n}.

?
Задача 36.7

Доказать, что в условиях предыдущей задачи последовательность ξn∗=ξ(Tn)\xi_{n}^{*} = \xi (T_{n}) является цепью Маркова. Найти вероятности перехода этой вложенной цепи Маркова.

?
Задача 36.8

Доказать, что для вложенной цепи Маркова ξn∗\xi_{n}^{*} состояния xix_{i} и xjx_{j} являются сообщающемися тогда и только тогда, когда pij(t)>0p_{ij}(t) > 0 при всех t>0t > 0.

?
Задача 36.9

Пусть скачкообразный марковский процесс со значениями 1, 2, 3, 4 и 5 определяется матрицей переходных интенсивностей

Q=(−410032−631000−550004−4000000). Q = \begin{pmatrix} -4 & 1 & 0 & 0 & 3 \\ 2 & -6 & 3 & 1 & 0 \\ 0 & 0 & -5 & 5 & 0 \\ 0 & 0 & 4 & -4 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{pmatrix}.

Найти матрицу переходных вероятностей и классифицировать состояния.

?
Задача 36.10

Некоторая частица перемещается между тремя состояниями как однородный скачкообразный марковский процесс с матрицей переходных интенсивностей

Q=(−2λλλλ−2λλλλ−2λ). Q = \begin{pmatrix} -2\lambda & \lambda & \lambda \\ \lambda & -2\lambda & \lambda \\ \lambda & \lambda & -2\lambda \end{pmatrix}.

Найти вероятности перехода за время tt и их пределы при t→∞t \to \infty.

?
Задача 36.11

Однородный скачкообразный марковский процесс с тремя значениями задаётся матрицей переходных интенсивностей

Q=(−(λ+ν)λνμ−(μ+ν)ν000). Q = \begin{pmatrix} -(\lambda +\nu ) & \lambda & \nu \\ \mu & -(\mu +\nu ) & \nu \\ 0 & 0 & 0 \end{pmatrix}.
?
(а)

Найти вероятности перехода за время tt и их пределы при t→∞t \to \infty.

(б)

Найти матрицу переходных вероятностей вложенной цепи Маркова и найти её стационарное распределение.

Задача 36.12

Однородный скачкообразный марковский процесс со значениями 1, 2, 3, 4, 5 и 6 задан матрицей переходных интенсивностей

Q=(−5111110−4130020−42000000000000−220003−30). Q = \begin{pmatrix} -5 & 1 & 1 & 1 & 1 & 1 \\ 0 & -4 & 1 & 3 & 0 & 0 \\ 2 & 0 & -4 & 2 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & -2 & 2 \\ 0 & 0 & 0 & 3 & -3 & 0 \end{pmatrix}.
?
(а)

Найти матрицу переходных вероятностей вложенной цепи Маркова и с её помощью провести классификацию состояний и найти вероятности поглощения каждым из существенных классов состояний.

(б)

Найти все стационарные распределения для марковского процесса и для вложенной цепи Маркова. Совпадают ли они? Почему?

Задача 36.13

Динамика системы задаётся скачкообразным марковским процессом с переходными интенсивностями

Q=(−3123−4142−6). Q = \begin{pmatrix} -3 & 1 & 2 \\ 3 & -4 & 1 \\ 4 & 2 & -6 \end{pmatrix}.

Показать, что все состояния процесса являются сообщающими, и найти его стационарное распределение.

?
Задача 36.14

Пусть ξ(t)\xi (t) — целочисленное случайное блуждание в непрерывном времени с задерживающим барьером в состоянии 0, т.е. скачкообразный марковский процесс с переходными интенсивностями qi,i+1=λ>0q_{i,i+1} = \lambda > 0 при i≥0i \geq 0, qi,i−1=μ>0q_{i,i-1} = \mu > 0 при i≥1i \geq 1, qii=−(λ+μ)q_{ii} = -(\lambda +\mu ) при i≥1i \geq 1, q00=−λq_{00} = -\lambda, qij=0q_{ij} = 0 для всех других j≠ij \neq i. Показать, что стационарное распределение процесса ξ(t)\xi (t) существует тогда и только тогда, когда λ<μ\lambda < \mu, и найти его в этом случае.

?
Задача 36.15

Пусть ξ(t)\xi (t) — процесс рождения и гибели в непрерывном времени в пространстве состояний Z+\mathbb {Z}^{+}, т.е. однородный скачкообразный марковский процесс с переходными интенсивностями qi,i+1=λiq_{i,i+1} = \lambda_{i} при i≥0i \geq 0 (вероятность рождения в состоянии ii), qi,i−1=μiq_{i,i-1} = \mu_{i} при i≥1i \geq 1 (вероятность гибели), qii=−(λ+μ)iq_{ii} = -(\lambda +\mu )_{i} при i≥1i \geq 1, q00=0q_{00} = 0, qij=0q_{ij} = 0 для всех других j≠ij \neq i. Под λ\lambda можно понимать интенсивность рождения потомка одним индивидуумом, а под μ\mu — интенсивность гибели одного индивидуума. С помощью вложенной цепи Маркова показать, что 0 является поглощающим состоянием процесса, и найти вероятность поглощения в нуле процесса, стартующего из произвольного состояния i>0i > 0.

?
Задача 36.16

На фабрике MM машин, каждая из работающих в момент времени tt машин с вероятностью μh+o(h)\mu h + o(h) может сломаться в интервале времени (t,t+h)(t, t+h) — независимо от того, что происходило до момента времени tt. Инженер ремонтирует каждую из сломанных машин. Если в момент tt есть хотя бы одна сломанная машина, то вероятность того, что в интервале времени (t,t+h)(t, t+h) ровно одна из них будет отремонтирована, есть λh+o(h)\lambda h + o(h) — опять независимо от истории процесса до момента времени tt и от работающих машин. Пусть ξ(t)\xi (t) — количество машин в рабочем состоянии в момент времени tt.

?
(а)

Найти матрицу переходных интенсивностей QQ процесса ξ(t)\xi (t).

(б)

Найти стационарное распределение процесса.

(в)

Найти слабый предел распределения числа работающих машин при M→∞M \to \infty.

Задача 36.17

Однородный скачкообразный марковский процесс ξ(t)\xi (t) со значениями 1, 2, 3, 4, 5 и 6 задаётся матрицей переходных интенсивностей

Q=(−2101001−4110100−4040000−2020040−4000200−2). Q = \begin{pmatrix} -2 & 1 & 0 & 1 & 0 & 0 \\ 1 & -4 & 1 & 1 & 0 & 1 \\ 0 & 0 & -4 & 0 & 4 & 0 \\ 0 & 0 & 0 & -2 & 0 & 2 \\ 0 & 0 & 4 & 0 & -4 & 0 \\ 0 & 0 & 2 & 0 & 0 & -2 \end{pmatrix}.
?
(а)

Найти матрицу переходных вероятностей вложенной цепи Маркова ξn∗\xi_{n}^{*}. С её помощью провести классификацию состояний и найти вероятности поглощения каждым из существенных классов.

(б)

Найти все стационарные распределения процесса ξ(t)\xi (t) и вложенной цепи Маркова ξn∗\xi_{n}^{*}. Совпадают ли они? Почему?

(в)

Найти предельные распределения для марковского процесса и для цепи Маркова, если ξ(0)=2\xi (0) = 2.

(г)

Найти предельные распределения для марковского процесса и для цепи Маркова, если ξ(0)=3\xi (0) = 3.