Глава 6

Цепи Маркова

[190/100%]
Показать
LaTeX
§
Задача 6.1.1

Покажите, что любая последовательность независимых случайных величин, принимающих значения в счётном множестве SS, является марковской цепью. При каком условии эта цепь однородна?

?
Задача 6.1.2

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

?
(a)

Наибольшее число XnX_{n}, выпавшее к nn-му броску.

(b)

Число NnN_{n} шестёрок среди nn бросков.

(c)

В момент rr — время CrC_{r}, прошедшее с последней выпавшей шестёрки.

(d)

В момент rr — время BrB_{r} до следующей шестёрки.

Задача 6.1.3

Пусть {Sn:n≥0}\left\{ S_{n}: n \geq 0\right\} — простое случайное блуждание с S0=0S_{0} = 0; покажите, что Xn=∣Sn∣X_{n} = \left|S_{n}\right| определяет марковскую цепь, и найдите переходные вероятности этой цепи. Пусть Mn=max⁡{Sk:0≤k≤n}M_{n} = \max \left\{ S_{k}: 0 \leq k \leq n\right\}; покажите, что Yn=Mn−SnY_{n} = M_{n}-S_{n} определяет марковскую цепь. Что произойдёт, если S0≠0S_{0} \neq 0?

?
Задача 6.1.4

Пусть XX — марковская цепь, и пусть {nr:r≥0}\left\{ n_{r}: r \geq 0\right\} — неограниченная возрастающая последовательность натуральных чисел. Покажите, что Yr=XnrY_{r} = X_{n_{r}} образует (возможно, неоднородную) марковскую цепь. Найдите матрицу переходных вероятностей YY, когда nr=2rn_{r} = 2 r, а XX является:

?
(a)

простым случайным блужданием, и

(b)

ветвящимся процессом.

Задача 6.1.5

Пусть XX — марковская цепь на SS, и пусть I:Sn→{0,1}I: S^{n} \rightarrow \left\{ 0,1\right\}. Покажите, что распределение Xn,Xn+1,…X_{n}, X_{n+1}, \ldots при условии {I(X1,…,Xn)=1}∩{Xn=i}\left\{ I\left(X_{1}, \ldots , X_{n}\right) = 1\right\} \cap \left\{ X_{n} = i\right\} совпадает с распределением Xn,Xn+1,…X_{n}, X_{n+1}, \ldots при условии {Xn=i}\left\{ X_{n} = i\right\}.

?
Задача 6.1.6

Пусть XX — марковская цепь на SS, и пусть TT — случайная величина, принимающая значения в {0,1,2,…}\left\{ 0,1,2, \ldots \right\}, обладающая тем свойством, что индикаторная функция I{T=n}I_{\left\{ T = n\right\} } события T=nT = n является функцией величин X1,X2,…,XnX_{1}, X_{2}, \ldots , X_{n}. Такая случайная величина TT называется моментом остановки, и приведённое выше определение требует, чтобы вопрос о том, выполняется ли T=nT = n, был разрешим при знании лишь прошлого и настоящего, X0,X1,…,XnX_{0}, X_{1}, \ldots , X_{n}, без какой-либо дополнительной информации о будущем.

Покажите, что

P(XT+m=j∣Xk=xk для 0≤k<T,XT=i)=P(XT+m=j∣XT=i) \mathbb {P}\left(X_{T+m} = j \mid X_{k} = x_{k} \text{ для } 0 \leq k < T, X_{T} = i\right) = \mathbb {P}\left(X_{T+m} = j \mid X_{T} = i\right)

при m≥0,i,j∈Sm \geq 0, i, j \in S и любых последовательностях состояний (xk)\left(x_{k}\right).

?
Задача 6.1.7

Пусть XX — марковская цепь с пространством состояний SS, и предположим, что h:S→Th: S \rightarrow T — взаимно однозначное отображение. Покажите, что Yn=h(Xn)Y_{n} = h\left(X_{n}\right) определяет марковскую цепь на TT. Обязано ли это быть так, если hh не является взаимно однозначным?

?
Задача 6.1.8

Пусть XX и YY — марковские цепи на множестве Z\mathbb {Z} целых чисел.

?
(a)

Обязательно ли последовательность Zn=Xn+YnZ_{n} = X_{n}+Y_{n} является марковской цепью?

(b)

Является ли ZZ марковской цепью, если XX и YY — независимые цепи? Приведите доказательство или контрпример.

(c)

Покажите, что ZZ — марковская цепь, если XX и YY независимы друг от друга и имеют независимые приращения.

Задача 6.1.9

Пусть XX — марковская цепь. Какие из следующих последовательностей являются марковскими цепями?

?
(a)

Xm+rX_{m+r} при r≥0r \geq 0.

(b)

X2mX_{2 m} при m≥0m \geq 0.

(c)

Последовательность пар (Xn,Xn+1)\left(X_{n}, X_{n+1}\right) при n≥0n \geq 0.

Задача 6.1.10

Пусть XX — марковская цепь. Покажите, что при 1<r<n1 < r < n

P(Xr=k∣Xi=xi при i=1,2,…,r−1,r+1…,n)=P(Xr=k∣Xr−1=xr−1,Xr+1=xr+1) \begin{aligned} & \mathbb {P}\left(X_{r} = k \mid X_{i} = x_{i} \text{ при } i = 1,2, \ldots , r-1, r+1 \ldots , n\right) \\ & = \mathbb {P}\left(X_{r} = k \mid X_{r-1} = x_{r-1}, X_{r+1} = x_{r+1}\right) \end{aligned}
?
Задача 6.1.11

Пусть {Xn:n≥1}\left\{ X_{n}: n \geq 1\right\} — независимые одинаково распределённые целочисленные случайные величины. Пусть Sn=∑r=1nXrS_{n} = \sum_{r = 1}^{n} X_{r}, причём S0=0S_{0} = 0, Yn=Xn+Xn−1Y_{n} = X_{n}+X_{n-1} с X0=0X_{0} = 0, и Zn=∑r=0nSrZ_{n} = \sum_{r = 0}^{n} S_{r}. Какие из следующих последовательностей образуют марковские цепи:

?
(a)

SnS_{n},

(b)

YnY_{n},

(c)

ZnZ_{n},

(d)

последовательность пар (Sn,Zn)\left(S_{n}, Z_{n}\right)?

Задача 6.1.12

Стохастическая матрица P\mathbf{P} называется дважды стохастической, если ∑ipij=1\sum_{i} p_{i j} = 1 для всех jj. Она называется субстохастической, если ∑ipij≤1\sum_{i} p_{i j} \leq 1 для всех jj. Покажите, что если P\mathbf{P} стохастична (соответственно, дважды стохастична, субстохастична), то Pn\mathbf{P}^{n} стохастична (соответственно, дважды стохастична, субстохастична) при всех nn.

?
Задача 6.1.13

Пусть XX — марковская цепь на конечном пространстве состояний SS с матрицей переходных вероятностей P\mathbf{P}, и пусть C={Cj:j∈J}\mathcal{C} = \left\{ C_{j}: j \in J\right\} — разбиение SS. Положим Yn=jY_{n} = j, если Xn∈CjX_{n} \in C_{j}. Цепь XX называется C\mathcal{C}-укрупняемой, если YY является марковской цепью.

Покажите, что XX является C\mathcal{C}-укрупняемой тогда и только тогда, когда при a,b∈Ja, b \in J величина P(Xn+1∈Cb∣Xn=i)\mathbb {P}\left(X_{n+1} \in C_{b} \mid X_{n} = i\right) постоянна для i∈Cai \in C_{a}.

?
§
Задача 6.2.1

Пусть lij(n)=P(Xn=j,Xk≠i при 1≤k<n∣X0=i)l_{i j}(n) = \mathbb {P}\left(X_{n} = j, X_{k} \neq i \text{ при } 1 \leq k < n \mid X_{0} = i\right) — вероятность того, что цепь переходит из ii в jj за nn шагов, ни разу не вернувшись в ii. Положив

Lij(s)=∑n=1∞snlij(n) L_{i j}(s) = \sum _{n = 1}^{\infty } s^{n} l_{i j}(n)

покажите, что Pij(s)=Pii(s)Lij(s)P_{i j}(s) = P_{i i}(s) L_{i j}(s) при i≠ji \neq j. Выведите отсюда, что времена первого достижения и времена последнего выхода имеют одинаковое распределение для любой марковской цепи, для которой Pii(s)=Pjj(s)P_{i i}(s) = P_{j j}(s) при всех ii и jj. Приведите пример такой цепи.

?
Задача 6.2.2

Пусть XX — марковская цепь, содержащая поглощающее состояние ss, с которым сообщаются все остальные состояния ii в том смысле, что pis(n)>0p_{i s}(n) > 0 при некотором n=n(i)n = n(i). Покажите, что все состояния, кроме ss, невозвратны.

?
Задача 6.2.3

Покажите, что состояние ii возвратно тогда и только тогда, когда среднее число посещений цепью состояния ii при старте из ii бесконечно. Иными словами, ii возвратно тогда и только тогда, когда ∑npii(n)=∞\sum_{n} p_{i i}(n) = \infty.

?
Задача 6.2.4

Пусть Vj=∣{n≥1:Xn=j}∣V_{j} = \left|\left\{ n \geq 1: X_{n} = j\right\} \right| — число посещений марковской цепью XX состояния jj, и определим ηij=Pi(Vj=∞)\eta_{i j} = \mathbb {P}_{i}\left(V_{j} = \infty \right). Покажите, что:

?
(a)

ηii={1 если i возвратно, 0 если i невозвратно, \eta_{i i} = \begin{cases} 1 & \text{ если } i \text{ возвратно, } \\ 0 & \text{ если } i \text{ невозвратно, }\end{cases}

(b)

ηij={Pi(Tj<∞) если j возвратно, 0 если j невозвратно, \eta_{i j} = \begin{cases} \mathbb {P}_{i}\left(T_{j} < \infty \right) & \text{ если } j \text{ возвратно, } \\ 0 & \text{ если } j \text{ невозвратно, }\end{cases} где Tj=min⁡{n≥1:Xn=j}T_{j} = \min \left\{ n \geq 1: X_{n} = j\right\}.

Задача 6.2.5

Различные состояния i,ji, j марковской цепи называются симметричными, если

Pi(Tj<Ti)=Pj(Ti<Tj) \mathbb {P}_{i}\left(T_{j} < T_{i}\right) = \mathbb {P}_{j}\left(T_{i} < T_{j}\right)

где Ti=min⁡{n≥1:Xn=i}T_{i} = \min \left\{ n \geq 1: X_{n} = i\right\}. Покажите, что если X0=iX_{0} = i и пара i,ji, j симметрична, то ожидаемое число посещений jj до того, как цепь снова попадёт в ii, равно 1. [Ср. с цитатой, следующей за теоремой (3.10.18).]

?
Задача 6.2.6

Пусть XX — марковская цепь, и пусть TT — геометрическая случайная величина с P(T>n)=sn\mathbb {P}\left(T > n\right) = s^{n} при n≥0n \geq 0, независимая от XX. Рассматривая ожидаемое число посещений цепью XX заданного состояния до момента TT, докажите теорему (6.2.3).

?
Задача 6.2.7

Пусть XX — эргодическая марковская цепь, начинающаяся из aa, и предположим, что XX неприводима (в том смысле, что для состояний i,ji, j найдётся m≥0m \geq 0 такое, что pij(m)>0p_{i j}(m) > 0). Пусть a,b,ca, b, c — различные состояния, и пусть T(a,b,¬c)T(a, b, \neg c) — время до первого посещения XX состояния bb без промежуточного посещения cc. То есть, если XX посещает bb раньше cc, то T(a,b,¬c)T(a, b, \neg c) равно этому времени, а если XX посещает cc раньше bb, то полагаем T(a,b,¬c)=∞T(a, b, \neg c) = \infty. Пусть

G(a,b,¬c;s)=∑n=1∞snPa(T(a,b,¬c)=n) G(a, b, \neg c ; s) = \sum _{n = 1}^{\infty } s^{n} \mathbb {P}_{a}\left(T(a, b, \neg c) = n\right)

Покажите, что

G(a,b,¬c;s)=Fab−FacFcb1−FbcFcb G(a, b, \neg c ; s) = \frac{F_{a b}-F_{a c} F_{c b}}{1-F_{b c} F_{c b}}

где, например, Fab(s)F_{a b}(s) — производящая функция вероятностей времени первого достижения TabT_{a b} из aa в bb независимо от промежуточных посещений cc. Покажите далее, что

Pa(T(a,b,¬c)<∞)=μac+μcb−μabμbc+μcb \mathbb {P}_{a}\left(T(a, b, \neg c) < \infty \right) = \frac{\mu _{a c}+\mu _{c b}-\mu _{a b}}{\mu _{b c}+\mu _{c b}}

где, например, μab=Ea[Tab]\mu_{a b} = \mathbb {E}_{a}\left[T_{a b}\right].

?
§
Задача 6.3.1

Пусть XX — марковская цепь на {0,1,2,…}\left\{ 0,1,2, \ldots \right\} с матрицей переходных вероятностей, заданной как p0j=ajp_{0 j} = a_{j} при j≥0j \geq 0, pii=rp_{i i} = r и pi,i−1=1−rp_{i, i-1} = 1-r при i≥1i \geq 1. Классифицируйте состояния цепи и найдите их средние времена возврата.

?
Задача 6.3.2

Определите, является ли возвратным случайное блуждание по целым числам с переходными вероятностями pi,i+2=p,pi,i−1=1−pp_{i, i+2} = p, p_{i, i-1} = 1-p при всех ii.

?
Задача 6.3.3

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

?
(a)
[1−2p2p0p1−2pp02p1−2p] \left[\begin{smallmatrix} 1-2 p & 2 p & 0 \\ p & 1-2 p & p \\ 0 & 2 p & 1-2 p \end{smallmatrix}\right]
(b)
[0p01−p1−p0p001−p0pp01−p0] \left[\begin{smallmatrix} 0 & p & 0 & 1-p \\ 1-p & 0 & p & 0 \\ 0 & 1-p & 0 & p \\ p & 0 & 1-p & 0 \end{smallmatrix}\right]

В каждом случае вычислите pij(n)p_{i j}(n) и средние времена возврата состояний.

Задача 6.3.4

Частица совершает случайное блуждание по вершинам куба. На каждом шаге она остаётся на месте с вероятностью 14\frac{1}{4} либо переходит в одну из соседних вершин, каждая с вероятностью 14\frac{1}{4}. Пусть vv и ww — две диаметрально противоположные вершины. Если блуждание начинается в vv, найдите:

?
(a)

среднее число шагов до первого возвращения в vv,

(b)

среднее число шагов до первого посещения ww,

(c)

среднее число посещений ww до первого возвращения в vv.

Задача 6.3.5

В обозначениях упражнения (6.2.4) покажите, что

?
(a)

если i→ji \rightarrow j и ii возвратно, то ηij=ηji=1\eta_{i j} = \eta_{j i} = 1,

(b)

ηij=1\eta_{i j} = 1 тогда и только тогда, когда Pi(Tj<∞)=Pj(Tj<∞)=1\mathbb {P}_{i}\left(T_{j} < \infty \right) = \mathbb {P}_{j}\left(T_{j} < \infty \right) = 1.

Задача 6.3.6

Пусть TA=min⁡{n≥0:Xn∈A}T_{A} = \min \left\{ n \geq 0: X_{n} \in A\right\}, где XX — марковская цепь, а AA — подмножество пространства состояний SS, и пусть ηj=Pj(TA<∞)\eta_{j} = \mathbb {P}_{j}\left(T_{A} < \infty \right). Покажите, что

ηj={1 если j∈A∑k∈Spjkηk если j∉A \eta _{j} = \begin{cases} 1 & \text{ если } j \in A \\ \sum _{k \in S} p_{j k} \eta _{k} & \text{ если } j \notin A\end{cases}

Покажите далее, что если x=(xj:j∈S)\mathbf{x} = \left(x_{j}: j \in S\right) — произвольное неотрицательное решение этих уравнений, то xj≥ηjx_{j} \geq \eta_{j} при всех jj.

?
Задача 6.3.7

В обозначениях упражнения (6.3.6) положим ρj=Ej[TA]\rho_{j} = \mathbb {E}_{j}\left[T_{A}\right]. Покажите, что

ρj={0 если j∈A1+∑k∈Spjkρk если j∉A \rho _{j} = \begin{cases} 0 & \text{ если } j \in A \\ 1+\sum _{k \in S} p_{j k} \rho _{k} & \text{ если } j \notin A\end{cases}

и что если x=(xj:j∈S)\mathbf{x} = \left(x_{j}: j \in S\right) — произвольное неотрицательное решение этих уравнений, то xj≥ρjx_{j} \geq \rho_{j} при всех jj.

?
Задача 6.3.8

Пусть XX — неприводимая марковская цепь, и пусть AA — подмножество пространства состояний. Пусть SrS_{r} и TrT_{r} — последовательные моменты, в которые цепь входит в AA и посещает AA соответственно. Являются ли последовательности {XSr:r≥1},{XTr:r≥1}\left\{ X_{S_{r}}: r \geq 1\right\} ,\left\{ X_{T_{r}}: r \geq 1\right\} марковскими цепями? Что можно сказать о моментах, в которые цепь покидает AA?

?
Задача 6.3.9
?
(a)

Покажите, что для каждой пары состояний i,ji, j неприводимой апериодической цепи существует N=N(i,j)N = N(i, j) такое, что pij(r)>0p_{i j}(r) > 0 при всех r≥Nr \geq N.

(b)

Покажите, что существует функция ff такая, что если P\mathbf{P} — матрица переходных вероятностей неприводимой апериодической марковской цепи с nn состояниями, то pij(r)>0p_{i j}(r) > 0 для всех состояний i,ji, j и всех r≥f(n)r \geq f(n).

(c)

Покажите далее, что f(4)≥6f(4) \geq 6 и f(n)≥(n−1)(n−2)f(n) \geq (n-1)(n-2). [Указание: лемма о почтовой марке утверждает, что для взаимно простых a,ba, b наименьшее nn, такое что все целые числа, строго превосходящие nn, представимы в виде αa+βb\alpha a+\beta b при некоторых целых α,β≥0\alpha , \beta \geq 0, равно (a−1)(b−1)(a-1)(b-1).]

Задача 6.3.10

Урна первоначально содержит nn зелёных шаров и n+2n+2 красных шаров. Наугад выбирается шар: если он зелёный, то дополнительно удаляется красный шар, и оба они выбрасываются; если он красный, то он возвращается в урну вместе с ещё одним красным и ещё одним зелёным шаром. Это повторяется, пока в урне не останется зелёных шаров. Покажите, что вероятность того, что процесс завершится, равна 1/(n+1)1 /(n+1).

Теперь поменяем правила местами: если шар зелёный, он возвращается вместе с ещё одним зелёным и ещё одним красным шаром; если он красный, он выбрасывается вместе с зелёным шаром. Покажите, что ожидаемое число итераций до тех пор, пока не останется зелёных шаров, равно ∑j=1n(2j+1)=n(n+2)\sum_{j = 1}^{n}(2 j+1) = n(n+2). [Таким образом, небольшое возмущение простого симметричного случайного блуждания может быть положительно возвратным, тогда как исходное блуждание нуль-возвратно.]

?
§
Задача 6.4.1

Корректурный экземпляр книги читает бесконечная последовательность редакторов, проверяющих его на наличие ошибок. При каждом прочтении каждая ошибка обнаруживается с вероятностью pp; между прочтениями типография исправляет обнаруженные ошибки, но вносит случайное число новых ошибок (ошибки могут вноситься, даже если ни одна ошибка не была обнаружена). Предполагая обычную независимость и то, что числа новых ошибок после разных прочтений одинаково распределены, найдите выражение для производящей функции вероятностей стационарного распределения числа XnX_{n} ошибок после nn-го цикла «редактор--типография», если оно существует. Найдите его в явном виде, когда типография вносит на каждом этапе пуассоновски распределённое число ошибок.

?
Задача 6.4.2

Проделайте заново соответствующие части упражнений (6.3.1)--(6.3.4), используя новые доступные вам методы. В частности, для упражнения (6.3.3):

?
(a)

проделайте заново часть (a);

(b)

проделайте заново часть (b).

Задача 6.4.3

Пусть XnX_{n} — количество воды в водохранилище в полдень дня nn. В течение суток, начинающихся в этот момент, в водохранилище поступает количество воды YnY_{n}, а непосредственно перед полуднем каждого дня из него забирается ровно одна единица воды (если такое количество может быть найдено). Максимальная ёмкость водохранилища равна KK, а избыточный приток воды переливается и теряется. Предположим, что YnY_{n} — независимые одинаково распределённые случайные величины, и что при округлении до какой-то смехотворно малой единицы объёма все числа в этом упражнении являются неотрицательными целыми. Покажите, что (Xn)\left(X_{n}\right) — марковская цепь, и найдите её матрицу переходных вероятностей и выражение для её стационарного распределения через производящую функцию вероятностей GG величин YnY_{n}.

Найдите стационарное распределение, когда YY имеет производящую функцию вероятностей G(s)=p(1−qs)−1G(s) = p(1-q s)^{-1}.

?
Задача 6.4.4

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

?
Задача 6.4.5

Пусть (xi(n):i,n≥1)\left(x_{i}(n): i, n \geq 1\right) — ограниченное семейство вещественных чисел.

?
(a)

Покажите, что существует возрастающая последовательность натуральных чисел n1,n2,…n_{1}, n_{2}, \ldots, такая что lim⁡r→∞xi(nr)\lim_{r \rightarrow \infty } x_{i}\left(n_{r}\right) существует при всех ii.

(b)

Используя этот результат, докажите, что для неприводимой марковской цепи, если неверно, что pij(n)→0p_{i j}(n) \rightarrow 0 при n→∞n \rightarrow \infty для всех ii и jj, то существует последовательность (nr:r≥1)(n_{r}: r \geq 1) и вектор α(≠0)\mathbf{\alpha }( \neq \mathbf{0}), такие что pij(nr)→αjp_{i j}\left(n_{r}\right) \rightarrow \alpha_{j} при r→∞r \rightarrow \infty для всех ii и jj.

Задача 6.4.6

Случайное блуждание на графе. Частица совершает случайное блуждание по множеству вершин связного графа GG, который для простоты мы считаем не имеющим ни петель, ни кратных рёбер. На каждом шаге она переходит к соседу своей текущей позиции, причём каждый такой сосед выбирается с равной вероятностью. Если GG имеет η(<∞)\eta ( < \infty ) рёбер, покажите, что стационарное распределение задаётся формулой πv=dv/(2η)\pi_{v} = d_{v} /(2 \eta ), где dvd_{v} — степень вершины vv.

?
Задача 6.4.7

Покажите, что случайное блуждание на бесконечном бинарном дереве невозвратно.

?
Задача 6.4.8

В каждый момент времени n=0,1,2,…n = 0,1,2, \ldots в камеру попадает YnY_{n} частиц, где {Yn:n≥0}\left\{ Y_{n}: n \geq 0\right\} независимы и имеют распределение Пуассона с параметром λ\lambda. Времена жизни частиц независимы и имеют геометрическое распределение с параметром pp. Пусть XnX_{n} — число частиц в камере в момент времени nn. Покажите, что XX — марковская цепь, и найдите её стационарное распределение.

?
Задача 6.4.9

Случайная последовательность выпуклых многоугольников строится следующим образом: наугад выбираются два ребра текущего многоугольника, их середины соединяются, и один из двух получившихся меньших многоугольников наугад выбирается в качестве следующего члена последовательности. Пусть Xn+3X_{n}+3 — число рёбер nn-го построенного таким образом многоугольника. Найдите E[Xn]\mathbb {E}\left[X_{n}\right] через X0X_{0} и найдите стационарное распределение марковской цепи XX.

?
Задача 6.4.10

Пусть ss — состояние неприводимой марковской цепи на неотрицательных целых числах. Покажите, что цепь возвратна, если существует решение y\mathbf{y} уравнений yi≥∑j:j≠spijyj,i≠sy_{i} \geq \sum_{j: j \neq s} p_{i j} y_{j}, i \neq s, удовлетворяющее условию yi→∞y_{i} \rightarrow \infty.

?
Задача 6.4.11

Частица совершает случайное блуждание по «галстуку-бабочке» ABCDE, изображённому ниже слева, где C — узел. Из любой вершины её следующий шаг с равной вероятностью ведёт в любую соседнюю вершину. Первоначально она находится в A. Найдите ожидаемое значение:

?
(a)

момента первого возвращения в A,

(b)

числа посещений D до возвращения в A,

(c)

числа посещений C до возвращения в A,

(d)

момента первого возвращения в A, при условии, что частица не посещала E,

(e)

числа посещений D до возвращения в A, при условии, что частица не посещала E.

Задача 6.4.12

Частица стартует из AA и совершает симметричное случайное блуждание по графу, изображённому выше справа. Найдите ожидаемое число посещений BB до возвращения в AA.

?
Задача 6.4.13

Колода содержит 52 карты с метками 1,2,…,521,2, \ldots , 52, и первоначально они расположены в возрастающем порядке сверху вниз. На каждом шаге тасования верхняя карта перекладывается на одно из 52 возможных мест, определяемых остальными 51 картами, причём это место выбирается равномерно случайно, независимо от всех предыдущих шагов. Найдите среднее число шагов до того момента, когда карта 52 впервые окажется наверху.

Покажите, что после момента, в который карта с меткой 52 случайно вкладывается сверху, порядок карт в колоде равномерно распределён по 52! возможностям.

?
Задача 6.4.14

Колода содержит 52 карты с метками 1,2,…,521,2, \ldots , 52, и первоначально они расположены в возрастающем порядке сверху вниз. На каждом шаге из колоды равномерно случайно выбирается карта и кладётся наверх, независимо от всех предыдущих шагов. Найдите среднее число шагов до того момента, когда каждая карта была выбрана хотя бы один раз.

Покажите, что после момента, в который последняя выбранная случайным образом карта кладётся наверх, порядок карт в колоде равномерно распределён по 52! возможностям.

?
Задача 6.4.15

Дик и Джим по очереди пишут упражнения для включения в учебник. Дик их пишет, а Джим проверяет. Каждое упражнение содержит ошибку с вероятностью pp, независимо от остальных упражнений. У Джима есть два режима работы. В режиме A он проверяет каждое упражнение по мере его написания. В режиме B он проверяет каждое упражнение с вероятностью 1/r1 / r, где r>1r > 1, независимо от всех прочих событий.

Пусть N≥1N \geq 1. Джим работает в режиме A до тех пор, пока не обнаружит NN подряд идущих упражнений без ошибок, после чего переходит в режим B. В режиме B он работает до тех пор, пока не будет найдено первое упражнение с ошибкой, после чего возвращается в режим A.

Пусть XX — марковская цепь, находящаяся в состоянии ii, если Джим работает в режиме A и последние ii подряд идущих упражнений с момента перехода в режим A оказались без ошибок, и находящаяся в состоянии NN, если Джим находится в режиме B.

?
(a)

Запишите переходные вероятности XX и найдите её стационарное распределение.

(b)

Покажите, что доля проверяемых упражнений в долгосрочной перспективе равна 1/[1+(r−1)(1−pN)]1 /\left[1+(r-1)\left(1-p^{N}\right)\right].

(c)

Найдите выражение для доли содержащих ошибку упражнений, которые Джим в долгосрочной перспективе не обнаруживает.

Задача 6.4.16

11.39), продолжение. Частица совершает случайное блуждание по неотрицательным целым числам следующим образом. Находясь в позиции k≥0k \geq 0, она переходит в следующую позицию, равномерно распределённую на множестве {0,1,…,k,k+1}\left\{ 0,1, \ldots , k, k+1\right\}. Покажите, что последовательность позиций образует апериодическую положительно возвратную марковскую цепь, и найдите её стационарное распределение.

Найдите среднее число шагов μ\mu, необходимых для первого достижения позиции 0, начиная с позиции 1.

?
Задача 6.4.17

Пусть {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} — неприводимая марковская цепь с пространством состояний SS и матрицей переходных вероятностей P\mathbf{P} (цепь может быть как невозвратной, так и возвратной). Пусть k∈Sk \in S, и пусть x\mathbf{x} — стационарная мера, такая что xk=1x_{k} = 1. Докажите, что x≥ρ(k)\mathbf{x} \geq \rho (k), где ρ(k)\rho (k) задаётся уравнением (6.4.5). Если цепь возвратна, покажите, что x=ρ(k)\mathbf{x} = \rho (k).

?
§
Задача 6.5.1

Случайное блуждание на множестве {0,1,2,…,b}\left\{ 0,1,2, \ldots , b\right\} имеет матрицу переходных вероятностей, заданную как p00=1−λ0,pbb=1−μb,pi,i+1=λip_{00} = 1-\lambda_{0}, p_{b b} = 1-\mu_{b}, p_{i, i+1} = \lambda_{i} и pi+1,i=μi+1p_{i+1, i} = \mu_{i+1} для 0≤i<b0 \leq i < b, где 0<λi,μi<10 < \lambda_{i}, \mu_{i} < 1 для всех ii, и λi+μi=1\lambda_{i}+\mu_{i} = 1 для 1≤i<b1 \leq i < b. Покажите, что этот процесс обратим в равновесии.

?
Задача 6.5.2

Пусть XX — неприводимая, положительно возвратная, апериодическая марковская цепь на пространстве состояний SS.

?
(a)

Критерий обратимости Колмогорова. Покажите, что XX обратима в равновесии тогда и только тогда, когда

pj1,j2pj2,j3⋯pjn−1,jnpjn,j1=pj1,jnpjn,jn−1⋯pj2,j1 p_{j_{1}, j_{2}} p_{j_{2}, j_{3}} \cdots p_{j_{n-1}, j_{n}} p_{j_{n}, j_{1}} = p_{j_{1}, j_{n}} p_{j_{n}, j_{n-1}} \cdots p_{j_{2}, j_{1}}

для всех nn и всех конечных последовательностей состояний j1,j2,…,jnj_{1}, j_{2}, \ldots , j_{n}.

(b)

Условие обратимости Келли. Покажите, что XX обратима в равновесии, если для всех различных троек i,j,k∈Si, j, k \in S

pijpjkpki=pikpkjpji p_{i j} p_{j k} p_{k i} = p_{i k} p_{k j} p_{j i}

и, кроме того, существует c∈Sc \in S такое, что pic>0p_{i c} > 0 для всех i≠ci \neq c.

(c)

Рассмотрим цепь с n≥3n \geq 3 состояниями. Покажите, что критерий Колмогорова в приведённой выше форме может потребовать проверки вплоть до 12∑r=3n(nr)(r−1)\frac{1}{2} \sum_{r = 3}^{n}\binom {n}{r}(r-1) ! уравнений, тогда как условие Келли, если оно применимо, требует не более (n−13)\binom {n-1}{3}.

(d)

Покажите, что случайное блуждание на конечном дереве обратимо в равновесии.

Задача 6.5.3

Пусть XX — обратимая марковская цепь, и пусть CC — непустое подмножество пространства состояний SS. Определим марковскую цепь YY на SS матрицей переходных вероятностей Q=(qij)\mathbf{Q} = \left(q_{i j}\right), где

qij={βpij если i∈C и j∉C,pij в противном случае,  q_{i j} = \begin{cases} \beta p_{i j} & \text{ если } i \in C \text{ и } j \notin C, \\ p_{i j} & \text{ в противном случае, }\end{cases}

для i≠ji \neq j, и где β\beta — постоянная, удовлетворяющая 0<β<10 < \beta < 1. Диагональные элементы qiiq_{i i} подобраны так, что Q\mathbf{Q} является стохастической матрицей. Покажите, что YY обратима в равновесии, и найдите её стационарное распределение. Опишите ситуацию в пределе при β↓0\beta \downarrow 0.

?
Задача 6.5.4

Может ли обратимая цепь быть периодической?

?
Задача 6.5.5

Модель «собака и блохи» из примера (6.5.5) — это марковская цепь XX на пространстве состояний {0,1,…,m}\left\{ 0,1, \ldots , m\right\} с переходными вероятностями

pi,i+1=1−im,pi,i−1=im, для 0≤i≤m p_{i, i+1} = 1-\frac{i}{m}, \quad p_{i, i-1} = \frac{i}{m}, \quad \text{ для } \quad 0 \leq i \leq m

Покажите, что если X0=iX_{0} = i, то

E[Xn−m2]=(i−m2)(1−2m)n→0 при n→∞ \mathbb {E}\left[X_{n}-\frac{m}{2}\right] = \left(i-\frac{m}{2}\right)\left(1-\frac{2}{m}\right)^{n} \rightarrow 0 \quad \text{ при } n \rightarrow \infty
?
Задача 6.5.6

Какие из следующих (стационарных) цепей являются обратимыми марковскими цепями?

?
(a)

Цепь X={Xn}X = \left\{ X_{n}\right\} с матрицей переходных вероятностей P=(1−ααβ1−β)\mathbf{P} = \left(\begin{array}{cc}1-\alpha & \alpha \\ \beta & 1-\beta \end{array}\right), где α+β>0\alpha +\beta > 0.

(b)

Цепь Y={Yn}Y = \left\{ Y_{n}\right\} с матрицей переходных вероятностей P=(0p1−p1−p0pp1−p0)\mathbf{P} = \left(\begin{array}{ccc}0 & p & 1-p \\ 1-p & 0 & p \\ p & 1-p & 0\end{array}\right), где 0<p<10 < p < 1.

(c)

Zn=(Xn,Yn)Z_{n} = \left(X_{n}, Y_{n}\right), где XnX_{n} и YnY_{n} независимы и удовлетворяют (a) и (b).

Задача 6.5.7

Пусть Xn,YnX_{n}, Y_{n} — независимые простые случайные блуждания. Пусть ZnZ_{n} — пара (Xn,Yn)(X_{n}, Y_{n}), урезанная так, чтобы лежать в области Xn≥0,Yn≥0,Xn+Yn≤aX_{n} \geq 0, Y_{n} \geq 0, X_{n}+Y_{n} \leq a, где aa — целое число. Найдите стационарное распределение ZnZ_{n}.

?
Задача 6.5.8

Покажите, что неприводимая марковская цепь с конечным пространством состояний и матрицей переходных вероятностей P\mathbf{P} обратима в равновесии тогда и только тогда, когда P=DS\mathbf{P} = \mathbf{D S} для некоторой симметричной матрицы S\mathbf{S} и диагональной матрицы D\mathbf{D} со строго положительными диагональными элементами. Покажите далее, что для обратимости в равновесии необходимо, но не достаточно, чтобы P\mathbf{P} имела вещественные собственные значения.

?
Задача 6.5.9

Случайное блуждание на графе. Пусть GG — конечный связный граф без петель и кратных рёбер, и пусть XX — случайное блуждание на GG, как в упражнении (6.4.6). Покажите, что XX обратимо в равновесии.

?
Задача 6.5.10

Рассмотрим случайное блуждание по строго положительным целым числам с переходными вероятностями

pi,i−1=12⋅i+2i+1,pi,i+1=12⋅ii+1,i≥2 p_{i, i-1} = \frac{1}{2} \cdot \frac{i+2}{i+1}, \quad p_{i, i+1} = \frac{1}{2} \cdot \frac{i}{i+1}, \quad i \geq 2

и p11=34,p12=14p_{11} = \frac{3}{4}, p_{12} = \frac{1}{4}. Покажите, что это блуждание положительно возвратно, и найдите среднее время возврата состояния ii.

?
Задача 6.5.11

Случайный жук совершает случайное блуждание по пяти вершинам, состоящим из главных точек компаса (обозначенных n,e,s,wn, e, s, w) и центра (обозначенного cc), с переходными вероятностями

pen=pws=14,pcn=pne=pec=pce=psw=pwn=18,pnc=pes=psc=pcs=pse=pwc=116,pcw=pnw=132. \begin{aligned} p_{e n} & = p_{w s} = \frac{1}{4}, \\ p_{c n} = p_{n e} = p_{e c} & = p_{c e} = p_{s w} = p_{w n} = \frac{1}{8}, \\ p_{n c} = p_{e s} = p_{s c} & = p_{c s} = p_{s e} = p_{w c} = \frac{1}{16}, \\ p_{c w} & = p_{n w} = \frac{1}{32}. \end{aligned}

Прочие переходы имеют вероятность 0. Покажите, что среднее время возврата центра равно μc=112\mu_{c} = \frac{11}{2}.

?
Задача 6.5.12

Пусть XX — неприводимая (но не обязательно апериодическая) марковская цепь на счётном пространстве состояний SS с матрицей переходных вероятностей P\mathbf{P} и инвариантным распределением π\pi. Пусть a∈(0,1)a \in (0,1), и пусть L=aP+(1−a)I\mathbf{L} = a \mathbf{P}+(1-a) \mathbf{I}, где I\mathbf{I} — единичная матрица.

?
(a)

Покажите, что L\mathbf{L} является матрицей переходных вероятностей неприводимой апериодической марковской цепи YY с инвариантным распределением π\pi.

(b)

Покажите, что если XX обратима в равновесии, то и YY обратима.

§
Задача 6.6.1

Теорема Маркова--Какутани утверждает, что для любого выпуклого компактного подмножества CC пространства Rn\mathbb {R}^{n} и любого линейного непрерывного отображения TT множества CC в себя TT имеет неподвижную точку (в том смысле, что T(x)=xT(x) = x для некоторого x∈Cx \in C). Используя это, докажите, что конечная стохастическая матрица обладает неотрицательным ненулевым левым собственным вектором, соответствующим собственному значению 1.

?
Задача 6.6.2

Пусть T\mathbf{T} — матрица размера m×nm \times n, и пусть v∈Rn\mathbf{v} \in \mathbb {R}^{n}. Теорема Фаркаша утверждает, что выполняется ровно одно из следующих условий:

(i) существует x∈Rm\mathbf{x} \in \mathbb {R}^{m}, такое что x≥0\mathbf{x} \geq \mathbf{0} и xT=v\mathbf{x} \mathbf{T} = \mathbf{v},

(ii) существует y∈Rn\mathbf{y} \in \mathbb {R}^{n}, такое что yv′<0\mathbf{y v}^{\prime } < 0 и Ty′≥0\mathbf{T y}^{\prime } \geq \mathbf{0}.

Используя это, докажите, что конечная стохастическая матрица обладает неотрицательным ненулевым левым собственным вектором, соответствующим собственному значению 1.

?
Задача 6.6.3

Предположим, что вы делаете ставки на скачки с mm возможными исходами. Имеется nn букмекеров, и единичная ставка у ii го букмекера приносит tijt_{i j}, если наступает jj й исход скачек. Вектор x=(x1,x2,…,xn)\mathbf{x} = \left(x_{1}, x_{2}, \ldots , x_{n}\right), где xr∈(−∞,∞)x_{r} \in (-\infty , \infty ) — ваша ставка у rr го букмекера, называется схемой ставок. Покажите, что выполняется ровно одно из (a) и (b):

(a) существует функция вероятностей p=(p1,p2,…,pm)\mathbf{p} = \left(p_{1}, p_{2}, \ldots , p_{m}\right), такая что ∑j=1mtijpj=0\sum_{j = 1}^{m} t_{i j} p_{j} = 0 при всех значениях ii,

(b) существует схема ставок x\mathbf{x}, при которой вы наверняка выигрываете, то есть ∑i=1nxitij>0\sum_{i = 1}^{n} x_{i} t_{i j} > 0 при всех jj.

?
Задача 6.6.4

Пусть XX — марковская цепь с пространством состояний S={1,2,3}S = \left\{ 1,2,3\right\} и матрицей переходных вероятностей

P=[1−pp001−ppp01−p] \mathbf{P} =\left[\begin{smallmatrix} 1-p & p & 0 \\ 0 & 1-p & p \\ p & 0 & 1-p \end{smallmatrix}\right]

где 0<p<10 < p < 1. Докажите, что

Pn=[a1na2na3na3na1na2na2na3na1n] \mathbf{P}^{n} =\left[\begin{smallmatrix} a_{1 n} & a_{2 n} & a_{3 n} \\ a_{3 n} & a_{1 n} & a_{2 n} \\ a_{2 n} & a_{3 n} & a_{1 n} \end{smallmatrix}\right]

где a1n+ωa2n+ω2a3n=(1−p+pω)na_{1 n}+\omega a_{2 n}+\omega^{2} a_{3 n} = (1-p+p \omega )^{n}, а ω\omega — комплексный кубический корень из 1.

?
Задача 6.6.5

Пусть P\mathbf{P} — матрица переходных вероятностей марковской цепи с конечным пространством состояний. Пусть I\mathbf{I} — единичная матрица, U\mathbf{U} — матрица размера ∣S∣×∣S∣\left|S\right| \times \left|S\right|, все элементы которой равны единице, а 1\mathbf{1} — вектор-строка длины ∣S∣\left|S\right|, все элементы которой равны единице. Пусть π\mathbf{\pi } — неотрицательный вектор с ∑iπi=1\sum_{i} \pi_{i} = 1. Покажите, что πP=π\mathbf{\pi } \mathbf{P} = \mathbf{\pi } тогда и только тогда, когда π(I−P+U)=1\mathbf{\pi }(\mathbf{I}-\mathbf{P}+\mathbf{U}) = \mathbf{1}. Выведите отсюда, что если P\mathbf{P} неприводима, то π=1(I−P+U)−1\mathbf{\pi } = \mathbf{1}(\mathbf{I}-\mathbf{P}+\mathbf{U})^{-1}.

?
Задача 6.6.6

Шахматная фигура совершает случайное блуждание по шахматной доске; на каждом шаге она с равной вероятностью делает любой из доступных ходов. Чему равно среднее время возврата угловой клетки, если фигура — это:

?
(a)

король?

(b)

ферзь?

(c)

слон?

(d)

конь?

(e)

ладья?

Задача 6.6.7

Ладья и слон совершают независимые симметричные случайные блуждания с синхронными шагами по доске 4×44 \times 4 (16 клеток). Если они начинают вместе в угловой клетке, покажите, что ожидаемое число шагов до их повторной встречи в той же угловой клетке равно 448/3448 / 3.

?
Задача 6.6.8

Найдите nn-шаговые переходные вероятности pij(n)p_{i j}(n) для цепи XX с матрицей переходных вероятностей

P=[0121213145122314112] \mathbf{P} =\left[\begin{smallmatrix} 0 & \frac{1}{2} & \frac{1}{2} \\ \frac{1}{3} & \frac{1}{4} & \frac{5}{12} \\ \frac{2}{3} & \frac{1}{4} & \frac{1}{12} \end{smallmatrix}\right]
?
Задача 6.6.9

Страницы всемирной сети образуют ориентированный граф WW с nn вершинами (представляющими страницы), соединёнными ориентированными рёбрами (представляющими ссылки). Наличие ссылки из ii в jj обозначается i→ji \rightarrow j, и граф задаётся своей матрицей смежности L=(lij)L = \left(l_{i j}\right), где lij=1l_{i j} = 1, если i→ji \rightarrow j, и lij=0l_{i j} = 0 в противном случае. Исходящая степень did_{i} (соответственно, входящая степень cic_{i}) вершины ii — это число ссылок, исходящих из ii (соответственно, ведущих в ii). Говорят, что вершина ii является тупиковой, если di=0d_{i} = 0.

Поведение быстро скучающего веб-сёрфера моделируется случайным блужданием по WW. Пусть b∈(0,1)b \in (0,1). Из любой тупиковой вершины блуждание переходит в случайно выбранную вершину WW, причём каждая вершина имеет вероятность 1/n1 / n. Находясь в нетупиковой вершине ii, блуждание с вероятностью b<1b < 1 переходит в случайную связанную вершину (каждая с вероятностью 1/di1 / d_{i}), а с вероятностью 1−b1-b переходит в случайную вершину WW (каждая с вероятностью 1/n1 / n).

?
(a)

Покажите, что матрицу переходных вероятностей P\mathbf{P} можно записать в виде P=bQ+v′e\mathbf{P} = b \mathbf{Q}+\mathbf{v}^{\prime } \mathbf{e}, где Q=(qij)\mathbf{Q} = \left(q_{i j}\right) с

qij={1/di если i→j1/n если i тупиковая0 в противном случае q_{i j} = \begin{cases} 1 / d_{i} & \text{ если } i \rightarrow j \\ 1 / n & \text{ если } i \text{ тупиковая} \\ 0 & \text{ в противном случае}\end{cases}

а v=(vi)\mathbf{v} = \left(v_{i}\right) — вектор-строка с vi=(1−b)/nv_{i} = (1-b) / n, и e\mathbf{e} — вектор-строка, все элементы которой равны 1.

(b)

Выведите отсюда, что стационарное распределение π\mathbf{\pi } задаётся формулой π={(1−b)/n}e(I−bQ)−1\mathbf{\pi } = \left\{ (1-b) / n\right\} \mathbf{e}(\mathbf{I}-b \mathbf{Q})^{-1}, где I\mathbf{I} — единичная матрица.

(c)

Объясните, почему элементы π\mathbf{\pi }, расположенные в порядке убывания, дают описание относительной популярности веб-страниц (называемой Google «PageRank» — это их товарный знак для запатентованного алгоритма).

Задача 6.6.10

Пусть P\mathbf{P} — матрица переходных вероятностей неприводимой марковской цепи на конечном пространстве состояний, и пусть π\mathbf{\pi } — левый собственный вектор P\mathbf{P}, соответствующий собственному значению 1. Покажите непосредственно из уравнения π=πP\mathbf{\pi } = \mathbf{\pi } \mathbf{P}, что элементы π\pi либо все положительны, либо все отрицательны, и тем самым докажите теорему (6.6.1d): существует единственное распределение π\mathbf{\pi }, удовлетворяющее πP=π\mathbf{\pi } \mathbf{P} = \mathbf{\pi }, причём все компоненты π\mathbf{\pi } строго положительны.

?
§
Задача 6.7.1

Пусть ZnZ_{n} — размер nn-го поколения ветвящегося процесса с Z0=1Z_{0} = 1 и P(Z1=k)=2−k\mathbb {P}\left(Z_{1} = k\right) = 2^{-k} при k≥0k \geq 0. Покажите непосредственно, что при n→∞,P(Zn≤2yn∣Zn>0)→1−e−2y,y>0n \rightarrow \infty , \mathbb {P}\left(Z_{n} \leq 2 y n \mid Z_{n} > 0\right) \rightarrow 1-e^{-2 y}, y > 0, в согласии с теоремой (6.7.8).

?
Задача 6.7.2

Пусть ZZ — надкритический ветвящийся процесс с Z0=1Z_{0} = 1 и производящей функцией размера семьи GG. Предположим, что вероятность вырождения η\eta удовлетворяет 0<η<10 < \eta < 1. Найдите способ описания процесса ZZ при условии его конечного вырождения.

?
Задача 6.7.3

Пусть ZnZ_{n} — размер nn-го поколения ветвящегося процесса с Z0=1Z_{0} = 1 и P(Z1=k)=qpk\mathbb {P}\left(Z_{1} = k\right) = q p^{k} при k≥0k \geq 0, где p+q=1p+q = 1 и p>12p > \frac{1}{2}. Используя ответ к упражнению (6.7.2), покажите, что при условии конечного вырождения ZZ процесс растёт подобно ветвящемуся процессу с размерами поколений Z~n\widetilde{Z}_{n}, удовлетворяющими Z~0=1\widetilde{Z}_{0} = 1 и P(Z~1=k)=pqk\mathbb {P}\left(\widetilde{Z}_{1} = k\right) = p q^{k} при k≥0k \geq 0.

?
Задача 6.7.4
?
(a)

Покажите, что E[X∣X>0]≤E[X2]/E[X]\mathbb {E}\left[X \mid X > 0\right] \leq \mathbb {E}\left[X^{2}\right] / \mathbb {E}\left[X\right] для любой неотрицательной случайной величины XX.

(b)

Пусть ZnZ_{n} — размер nn-го поколения ветвящегося процесса с Z0=1Z_{0} = 1 и P(Z1=k)=qpk\mathbb {P}\left(Z_{1} = k\right) = q p^{k} при k≥0k \geq 0, где p>12p > \frac{1}{2}. Используя пункт (a), покажите, что E[Zn/μn∣Zn>0]≤2p/(p−q)\mathbb {E}\left[Z_{n} / \mu^{n} \mid Z_{n} > 0\right] \leq 2 p /(p-q), где μ=p/q\mu = p / q.

(c)

Покажите, что в обозначениях пункта (b) E[Zn/μn∣Zn>0]→p/(p−q)\mathbb {E}\left[Z_{n} / \mu^{n} \mid Z_{n} > 0\right] \rightarrow p /(p-q) при n→∞n \rightarrow \infty.

§
Задача 6.8.1

Мухи и осы садятся на вашу тарелку с едой в соответствии с независимыми пуассоновскими процессами с интенсивностями λ\lambda и μ\mu соответственно. Покажите, что появления летающих объектов образуют пуассоновский процесс с интенсивностью λ+μ\lambda +\mu.

?
Задача 6.8.2

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

?
Задача 6.8.3

Пусть TnT_{n} — момент nn-го поступления в пуассоновском процессе NN с интенсивностью λ\lambda, и определим процесс избыточного времени жизни E[t]=TN(t)+1−t\mathbb {E}\left[t\right] = T_{N(t)+1}-t — время, которое нужно ждать после момента tt до следующего поступления. Покажите, обусловливая по T1T_{1}, что

P(E[t]>x)=e−λ(t+x)+∫0tP(E[t−u]>x)λe−λudu \mathbb {P}\left(\mathbb {E}\left[t\right] > x\right) = e^{-\lambda (t+x)}+\int _{0}^{t} \mathbb {P}\left(\mathbb {E}\left[t-u\right] > x\right) \lambda e^{-\lambda u} d u

Решите это интегральное уравнение, чтобы найти функцию распределения E[t]\mathbb {E}\left[t\right]. Объясните полученный результат.

?
Задача 6.8.4

Пусть BB — простой процесс рождения из пункта (6.8.15 b)(6.8.15 \mathrm{~ b}) с B(0)=IB(0) = I; интенсивности рождения равны λn=nλ\lambda_{n} = n \lambda. Запишите прямую систему уравнений для процесса и выведите отсюда, что

P(B(t)=k)=(k−1I−1)e−Iλt(1−e−λt)k−I,k≥I \mathbb {P}\left(B(t\right) = k) = \binom {k-1}{I-1} e^{-I \lambda t}\left(1-e^{-\lambda t}\right)^{k-I}, \quad k \geq I

Покажите также, что E[B(t])=Ieλt\mathbb {E}\left[B(t\right]) = I e^{\lambda t} и Var⁡[B(t])=Ie2λt(1−e−λt)\operatorname {Var}\left[B(t\right]) = I e^{2 \lambda t}\left(1-e^{-\lambda t}\right).

?
Задача 6.8.5

Пусть BB — процесс простого рождения с иммиграцией (6.8.11в) с параметрами λ\lambda и ν\nu и с B(0)=0B(0) = 0; интенсивности рождения равны λn=nλ+ν\lambda_{n} = n \lambda +\nu. Запишите последовательность дифференциально-разностных уравнений для pn(t)=P(B(t)=n)p_{n}(t) = \mathbb {P}\left(B(t) = n\right). Не решая эти уравнения, используйте их, чтобы показать, что m(t)=E[B(t)]m(t) = \mathbb {E}\left[B(t)\right] удовлетворяет уравнению m′(t)=λm(t)+νm^{\prime }(t) = \lambda m(t)+\nu, и решите его относительно m(t)m(t).

?
Задача 6.8.6

Пусть NN — процесс рождения с интенсивностями λ0,λ1,…\lambda_{0}, \lambda_{1}, \ldots, и пусть N(0)=0N(0) = 0. Покажите, что pn(t)=P(N(t)=n)p_{n}(t) = \mathbb {P}\left(N(t\right) = n) задаётся формулой

pn(t)=1λn∑i=0nλie−λit∏j=0j≠inλjλj−λi p_{n}(t) = \frac{1}{\lambda _{n}} \sum _{i = 0}^{n} \lambda _{i} e^{-\lambda _{i} t} \prod _{\substack {j = 0 \\ j \neq i}}^{n} \frac{\lambda _{j}}{\lambda _{j}-\lambda _{i}}

при условии, что λi≠λj\lambda_{i} \neq \lambda_{j} при i≠ji \neq j.

?
Задача 6.8.7

Предположим, что общий процесс рождения из предыдущего упражнения таков, что ∑nλn−1<∞\sum_{n} \lambda_{n}^{-1} < \infty. Покажите, что λnpn(t)→f(t)\lambda_{n} p_{n}(t) \rightarrow f(t) при n→∞n \rightarrow \infty, где ff — плотность случайной величины T=sup⁡{t:N(t)<∞}T = \sup \left\{ t: N(t) < \infty \right\}. Выведите отсюда, что E[N(t]∣N(t)<∞)\mathbb {E}\left[N(t\right] \mid N(t) < \infty ) конечно или бесконечно в зависимости от сходимости или расходимости ∑nnλn−1\sum_{n} n \lambda_{n}^{-1}.

Найдите преобразование Лапласа ff в замкнутой форме для случая, когда λn=(n+12)2\lambda_{n} = \left(n+\frac{1}{2}\right)^{2}, и выведите отсюда выражение для ff.

?
Задача 6.8.8

Светофор горит зелёным в момент времени 00 и впоследствии переключается между зелёным и красным в моменты пуассоновского процесса с интенсивностью λ\lambda. Начиная с момента времени x>0x > 0, пусть W(x)W(x) — время ожидания до первого включения зелёного света. Найдите распределение W(x)W(x).

?
Задача 6.8.9

Условное свойство простого процесса рождения. Пусть X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} — простой процесс рождения с интенсивностью λ\lambda, и пусть X(0)=1X(0) = 1. Пусть b≥1b \geq 1. Покажите, что при условии {X(t)=b+1}\left\{ X(t) = b+1\right\} моменты bb рождений имеют то же распределение, что и вариационный ряд случайной выборки объёма bb из плотности

f(x)=λe−λ(t−x)1−e−λt,0≤x≤t f(x) = \frac{\lambda e^{-\lambda (t-x)}}{1-e^{-\lambda t}}, \quad 0 \leq x \leq t
?
Задача 6.8.10

Валуны падают по жёлобу (кулуару) в моменты пуассоновского процесса с интенсивностью λ\lambda, а альпинисты поднимаются по жёлобу в моменты пуассоновского процесса с интенсивностью μ\mu (оба процесса независимы друг от друга). Если падение и подъём происходят в течение интервала длины cc или меньше, говорят, что произошло совпадение. Покажите, что время TT до первого совпадения имеет среднее

E[T]=1λ+μ{1+λ2+μ2λμ+2e−(λ+μ)c−e−2(λ+μ)c}/{1−e−2(λ+μ)c} \mathbb {E}\left[T\right] = \frac{1}{\lambda +\mu }\left\{ 1+\frac{\lambda ^{2}+\mu ^{2}}{\lambda \mu }+2 e^{-(\lambda +\mu ) c}-e^{-2(\lambda +\mu ) c}\right\} /\left\{ 1-e^{-2(\lambda +\mu ) c}\right\}

Покажите, что cE[T]→(2λμ)−1c \mathbb {E}\left[T\right] \rightarrow (2 \lambda \mu )^{-1} при c↓0c \downarrow 0. Можете ли вы доказать этот последний результат напрямую?

?
Задача 6.8.11

Пусть S1,S2S_{1}, S_{2} — моменты первых двух поступлений в пуассоновском процессе с интенсивностью λ\lambda, в порядке поступления. Покажите, что

P(s<S1≤t<S2)=λ(t−s)e−λt,0<s<t<∞ \mathbb {P}\left(s < S_{1} \leq t < S_{2}\right) = \lambda (t-s) e^{-\lambda t}, \quad 0 < s < t < \infty

и выведите отсюда совместную плотность S1S_{1} и S2S_{2}.

?
Задача 6.8.12

Внештатному продавцу платят RR единиц за каждую продажу, и комиссионные поступают в моменты пуассоновского процесса интенсивности λ\lambda. Расходы на жизнь расходуют его ресурсы с единичной скоростью. Если его начальное состояние равно S>0S > 0, покажите, что вероятность того, что он когда-либо обанкротится, равна aS/Ra^{S / R}, где aa — наименьшее x>0x > 0, такое что x=eλR(x−1)x = e^{\lambda R(x-1)}.

Докажите, что a∈(0,1)a \in (0,1), если λR>1\lambda R > 1, тогда как a=1a = 1, если λR<1\lambda R < 1.

?
§
Задача 6.9.1

Пусть λμ>0\lambda \mu > 0, и пусть XX — марковская цепь на {1,2}\left\{ 1,2\right\} с генератором

G=[−μμλ−λ] \mathbf{G} =\left[\begin{smallmatrix} -\mu & \mu \\ \lambda & -\lambda \end{smallmatrix}\right]
?
(a)

Запишите прямые уравнения и решите их относительно переходных вероятностей pij(t),i,j=p_{i j}(t), i, j = 1,2.

(b)

Вычислите Gn\mathbf{G}^{n} и с его помощью найдите ∑n=0∞(tn/n!)Gn\sum_{n = 0}^{\infty }\left(t^{n} / n!\right) \mathbf{G}^{n}. Сравните ваш ответ с ответом к пункту (a).

(c)

Решите уравнение πG=0\pi \mathbf{G} = \mathbf{0}, чтобы найти стационарное распределение. Проверьте, что pij(t)→πjp_{i j}(t) \rightarrow \pi_{j} при t→∞t \rightarrow \infty.

Задача 6.9.2

В продолжение предыдущего упражнения найдите:

?
(a)

P(X(t)=2∣X(0)=1,X(3t)=1)\mathbb {P}\left(X(t\right) = 2 \mid X(0) = 1, X(3 t) = 1),

(b)

P(X(t)=2∣X(0)=1,X(3t)=1,X(4t)=1)\mathbb {P}\left(X(t\right) = 2 \mid X(0) = 1, X(3 t) = 1, X(4 t) = 1).

Задача 6.9.3

Задания поступают в компьютерную очередь в соответствии с пуассоновским процессом интенсивности λ\lambda. Центральный процессор обрабатывает их одно за другим в порядке поступления, и время выполнения каждого имеет показательное распределение с параметром μ\mu, причём времена выполнения разных заданий независимы друг от друга и от процесса поступления. Пусть X(t)X(t) — число заданий в системе (выполняющихся или ожидающих) в момент времени tt, где X(0)=0X(0) = 0. Объясните, почему XX — марковская цепь, и запишите её генератор. Покажите, что стационарное распределение существует тогда и только тогда, когда λ<μ\lambda < \mu, и найдите его в этом случае.

?
Задача 6.9.4

Пусть X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} — марковская цепь со стационарным распределением π\mathbf{\pi }. Мы можем выбирать значения XX в моменты пуассоновского процесса: пусть NN — пуассоновский процесс с интенсивностью λ\lambda, независимый от XX, и определим Yn=X(Tn)Y_{n} = X\left(T_{n}\right). Покажите, что Y={Yn:n≥0}Y = \left\{ Y_{n}: n \geq 0\right\} — дискретная марковская цепь с тем же стационарным распределением, что и XX. (Это иллюстрирует свойство PASTA: пуассоновские поступления видят усреднённые по времени характеристики.) [Полное предположение о независимости NN и XX не является необходимым для этого вывода. Достаточно, чтобы {N(s):s≥t}\left\{ N(s): s \geq t\right\} было независимо от {X(s):s≤t}\left\{ X(s): s \leq t\right\} — свойство, известное как «отсутствие предвидения». Не требуется даже, чтобы XX была марковской; свойство PASTA выполняется для многих подходящих эргодических процессов.]

?
Задача 6.9.5

Пусть XX — марковская цепь с непрерывным временем с генератором G\mathbf{G}, удовлетворяющим gi=−gii>0g_{i} = -g_{i i} > 0 при всех ii. Пусть HA=inf⁡{t≥0:X(t)∈A}H_{A} = \inf \left\{ t \geq 0: X(t) \in A\right\} — время достижения множества состояний AA, и пусть ηj=Pj(HA<∞)\eta_{j} = \mathbb {P}_{j}\left(H_{A} < \infty \right) — вероятность когда-либо достичь AA, начав из jj. Используя свойства цепи скачков, которую можно считать «благополучной», покажите, что ∑kgjkηk=0\sum_{k} g_{j k} \eta_{k} = 0 при j∉Aj \notin A.

?
Задача 6.9.6

В продолжение предыдущего упражнения пусть μj=Ej[HA]\mu_{j} = \mathbb {E}_{j}\left[H_{A}\right]. Покажите, что вектор μ\mathbf{\mu } является минимальным неотрицательным решением уравнений

μj=0 если j∈A,1+∑k∈Sgjkμk=0 если j∉A \mu _{j} = 0 \quad \text{ если } j \in A, \quad 1+\sum _{k \in S} g_{j k} \mu _{k} = 0 \quad \text{ если } j \notin A
?
Задача 6.9.7

Пусть XX — марковская цепь с непрерывным временем и переходными вероятностями pij(t)p_{i j}(t), и определим Fi=inf⁡{t>T1:X(t)=i}F_{i} = \inf \left\{ t > T_{1}: X(t) = i\right\}, где T1T_{1} — момент первого скачка XX. Покажите, что если gii≠0g_{i i} \neq 0, то Pi(Fi<∞)=1\mathbb {P}_{i}\left(F_{i} < \infty \right) = 1 тогда и только тогда, когда ii возвратно.

?
Задача 6.9.8

Пусть XX — простое симметричное случайное блуждание по целым числам в непрерывном времени, так что

pi,i+1(h)=pi,i−1(h)=12λh+o(h) p_{i, i+1}(h) = p_{i, i-1}(h) = \frac{1}{2} \lambda h+\mathrm{o}(h)

Покажите, что это блуждание возвратно. Пусть TT — время, проведённое в mm за время экскурсии из 0. Найдите распределение TT.

?
Задача 6.9.9

Пусть ii — невозвратное состояние марковской цепи с непрерывным временем XX с X(0)=iX(0) = i. Покажите, что суммарное время, проведённое в состоянии ii, имеет показательное распределение.

?
Задача 6.9.10

Пусть XX — асимметричное простое случайное блуждание в непрерывном времени по неотрицательным целым числам с удержанием в 00, так что

pij(h)={λh+o(h) если j=i+1,i≥0μh+o(h) если j=i−1,i≥1 p_{i j}(h) = \begin{cases} \lambda h+\mathrm{o}(h) & \text{ если } j = i+1, i \geq 0 \\ \mu h+\mathrm{o}(h) & \text{ если } j = i-1, i \geq 1\end{cases}

Предположим, что X(0)=0X(0) = 0 и λ>μ\lambda > \mu. Покажите, что суммарное время VrV_{r}, проведённое в состоянии rr, имеет показательное распределение с параметром λ−μ\lambda -\mu.

Предположим теперь, что X(0)X(0) имеет некоторое общее распределение с производящей функцией вероятностей GG. Найдите ожидаемое количество времени, проведённого в 0, через GG.

?
Задача 6.9.11

Пусть XX — марковская цепь с непрерывным временем на конечном пространстве состояний SS с генератором G=(gij)\mathbf{G} = \left(g_{i j}\right). Покажите из первых принципов, что переходные вероятности удовлетворяют

pij(h)={1+giih+o(h) если i=jgijh+o(h) если i≠j p_{i j}(h) = \begin{cases} 1+g_{i i} h+\mathrm{o}(h) & \text{ если } i = j \\ g_{i j} h+\mathrm{o}(h) & \text{ если } i \neq j\end{cases}
?
Задача 6.9.12

Пусть XX — марковская цепь на целых числах Z\mathbb {Z} с генератором, удовлетворяющим gi,i−1=gi,i+1=2ig_{i, i-1} = g_{i, i+1} = 2^{i} при i∈Zi \in \mathbb {Z}, и gi,j=0g_{i, j} = 0 для остальных пар (i,j)(i, j) с i≠ji \neq j. Взрывается ли XX?

?
Задача 6.9.13

Популяция растёт под угрозой полного уничтожения. Она моделируется марковской цепью с генератором G=(gij)\mathbf{G} = \left(g_{i j}\right), удовлетворяющим

gi,i+1=1i+2,i≥0gi,0=1(i+1)(i+2),i≥1 \begin{aligned} g_{i, i+1} & = \frac{1}{i+2}, & & i \geq 0 \\ g_{i, 0} & = \frac{1}{(i+1)(i+2)}, & & i \geq 1 \end{aligned}

причём остальные внедиагональные элементы G\mathbf{G} равны 0. Покажите, что эта цепь нуль-возвратна.

?
Задача 6.9.14

Пусть Zn=X(nh)Z_{n} = X(n h), где XX — марковская цепь и h>0h > 0. Покажите, что ii возвратно для ZZ тогда и только тогда, когда оно возвратно для XX. Покажите, что ZZ неприводима тогда и только тогда, когда XX неприводима.

?
§
Задача 6.10.1

Пусть N<∞N < \infty, и пусть Q\mathcal{Q} — пространство матриц размера N×NN \times N с вещественными элементами, с нормой

∣Q∣=sup⁡x≠0∣Qx∣∣x∣,Q∈Q |Q| = \sup _{\mathbf{x} \neq \mathbf{0}} \frac{|Q \mathbf{x}|}{|\mathbf{x}|}, \quad Q \in \mathcal{Q}

где ∣y∣\left|\mathbf{y}\right| — евклидова норма вектора y\mathbf{y}, а супремум берётся по всем ненулевым векторам-столбцам.

?
(a)

Покажите для Q1,Q2∈QQ_{1}, Q_{2} \in \mathcal{Q}, что

∣Q1+Q2∣≤∣Q1∣+∣Q2∣,∣Q1Q2∣≤∣Q1∣⋅∣Q2∣ \left|Q_{1}+Q_{2}\right| \leq \left|Q_{1}\right|+\left|Q_{2}\right|, \quad \left|Q_{1} Q_{2}\right| \leq \left|Q_{1}\right| \cdot \left|Q_{2}\right|
(b)

Покажите для Q∈QQ \in \mathcal{Q}, что En:=∑k=0nQk/k!E_{n} := \sum_{k = 0}^{n} Q^{k} / k! сходится по норме ∣⋅∣\left|\cdot \right| к пределу, который мы обозначаем E=E(Q)E = E(Q).

(c)

Покажите, что E(Q1+Q2)=E(Q1)E(Q2)E(Q_{1}+Q_{2}) = E(Q_{1}) E(Q_{2}), если Q1Q_{1} и Q2Q_{2} — коммутирующие элементы Q\mathcal{Q}.

Задача 6.10.2

Пусть YY — неприводимая дискретная марковская цепь на счётно-бесконечном пространстве состояний SS с матрицей переходных вероятностей Y=(yij)\mathbf{Y} = \left(y_{i j}\right), удовлетворяющей yii=0y_{i i} = 0 для всех состояний ii, и со стационарным распределением ν\mathbf{\nu }. Постройте процесс XX с непрерывным временем на SS, для которого YY является цепью скачков, такой что XX не имеет стационарного распределения.

?
Задача 6.10.3

Пусть X=(X(t):t≥0)X = (X(t): t \geq 0) — марковская цепь на Z\mathbb {Z} с генератором G=(gij)\mathbf{G} = \left(g_{i j}\right), заданным как

gi,i−1=i2+1,gi,i=−2(i2+1),gi,i+1=i2+1,i∈Z g_{i, i-1} = i^{2}+1, \quad g_{i, i} = -2\left(i^{2}+1\right), \quad g_{i, i+1} = i^{2}+1, \quad i \in \mathbb {Z}

Покажите, что XX возвратна. Является ли XX положительно возвратной?

?
Задача 6.10.4

Пусть XX — марковская цепь на Z\mathbb {Z} с генератором G=(gij)G = \left(g_{i j}\right), удовлетворяющим

gi,i−1=3∣i∣,gi,i=−3∣i∣+1,gi,i+1=2⋅3∣i∣,i∈Z g_{i, i-1} = 3^{\left|i\right|}, \quad g_{i, i} = -3^{\left|i\right|+1}, \quad g_{i, i+1} = 2 \cdot 3^{\left|i\right|}, \quad i \in \mathbb {Z}

Покажите, что XX невозвратна, но обладает инвариантным распределением. Объясните это.

?
Задача 6.10.5

Пусть X=(Xt:t≥0)X = \left(X_{t}: t \geq 0\right) — марковская цепь с генератором G=(gij)\mathbf{G} = \left(g_{i j}\right) на конечном пространстве состояний SS, и пусть f:S→Rf: S \rightarrow \mathbb {R} — функция, которую мы отождествляем с вектором f=(f(i):i∈S)f = (f(i): i \in S). Покажите, что

Gf(i)=∑j∈Sgij(f(j)−f(i)) \mathbf{G} f(i) = \sum _{j \in S} g_{i j}(f(j)-f(i))

где Gf\mathbf{G} f обозначает обычное матричное умножение. Покажите, что

Gf(i)=lim⁡t→01t[Ei[f(Xt)]−f(i)],i∈S \mathbf{G} f(i) = \lim _{t \rightarrow 0} \frac{1}{t}\left[\mathbb {E}_{i}\left[f\left(X_{t}\right)\right]-f(i)\right], \quad i \in S

и выведите отсюда, что

Ei[f(Xt)]=f(i)+∫0tEi[Gf(Xs)]ds \mathbb {E}_{i}\left[f\left(X_{t}\right)\right] = f(i)+\int _{0}^{t} \mathbb {E}_{i}\left[\mathbf{G} f\left(X_{s}\right)\right] d s
?
§
Задача 6.11.1

Опишите цепь скачков для процесса рождения и гибели с интенсивностями λn\lambda_{n} и μn\mu_{n}.

?
Задача 6.11.2

Рассмотрим процесс иммиграции-гибели XX — процесс рождения и гибели с интенсивностями рождения λn=λ\lambda_{n} = \lambda и интенсивностями гибели μn=nμ\mu_{n} = n \mu. Найдите матрицу переходных вероятностей цепи скачков YY и покажите, что её стационарное распределение равно

πn=12(n!)(1+nρ)ρne−ρ \pi _{n} = \frac{1}{2(n!)}\left(1+\frac{n}{\rho }\right) \rho ^{n} e^{-\rho }

где ρ=λ/μ\rho = \lambda / \mu. Объясните, почему оно отличается от стационарного распределения XX.

?
Задача 6.11.3

Рассмотрим процесс рождения и гибели XX с λn=nλ\lambda_{n} = n \lambda и μn=nμ\mu_{n} = n \mu при всех n≥0n \geq 0. Предположим, что X(0)=1X(0) = 1, и пусть η(t)=P1(X(t)=0)\eta (t) = \mathbb {P}_{1}\left(X(t) = 0\right). Покажите, что η\eta удовлетворяет дифференциальному уравнению

η′(t)+(λ+μ)η(t)=μ+λη(t)2 \eta ^{\prime }(t)+(\lambda +\mu ) \eta (t) = \mu +\lambda \eta (t)^{2}

Отсюда найдите η(t)\eta (t) и вычислите P1(X(t)=0∣X(u)=0)\mathbb {P}_{1}\left( X(t) = 0 \mid X(u) = 0 \right) при 0<t<u0 < t < u.

?
Задача 6.11.4

Для процесса рождения и гибели из предыдущего упражнения с λ<μ\lambda < \mu покажите, что распределение X(t)X(t) при условии {X(t)>0}\left\{ X(t) > 0\right\} сходится при t→∞t \rightarrow \infty к геометрическому распределению.

?
Задача 6.11.5

Пусть XX — процесс рождения и гибели с λn=nλ\lambda_{n} = n \lambda и μn=nμ\mu_{n} = n \mu, и предположим, что X(0)=1X(0) = 1. Покажите, что момент времени TT, в который X(t)X(t) впервые принимает значение 0, удовлетворяет

E[T∣T<∞]={1λlog⁡(μμ−λ) если λ<μ1μlog⁡(λλ−μ) если λ>μ \mathbb {E}\left[T \mid T < \infty \right] = \begin{cases} \frac{1}{\lambda } \log \left(\frac{\mu }{\mu -\lambda }\right) & \text{ если } \lambda < \mu \\ \frac{1}{\mu } \log \left(\frac{\lambda }{\lambda -\mu }\right) & \text{ если } \lambda > \mu \end{cases}

Что происходит при λ=μ\lambda = \mu?

?
Задача 6.11.6

Пусть XX — процесс рождения и гибели из упражнения (6.11.5) с λ≠μ\lambda \neq \mu, и пусть Vr(t)V_{r}(t) — суммарное время, проведённое процессом в состоянии r≥0r \geq 0 до момента tt. Найдите распределение V1(∞)V_{1}(\infty ) и производящую функцию ∑rsrE[Vr(t)]\sum_{r} s^{r} \mathbb {E}\left[V_{r}(t)\right]. С их помощью покажите двумя способами, что E[V1(∞)]=[max⁡{λ,μ}]−1\mathbb {E}\left[V_{1}(\infty )\right] = [\max \left\{ \lambda , \mu \right\} ]^{-1}. Покажите далее, что E[Vr(∞)]=λr−1r−1[max⁡{λ,μ}]−r\mathbb {E}\left[V_{r}(\infty )\right] = \lambda^{r-1} r^{-1}[\max \left\{ \lambda , \mu \right\} ]^{-r}.

?
Задача 6.11.7

Повторите вычисления упражнения (6.11.6) в случае λ=μ\lambda = \mu.

?
Задача 6.11.8

Рассмотрим процесс рождения и гибели XX с интенсивностями рождения λn>0\lambda_{n} > 0 при n≥0n \geq 0, интенсивностями гибели μn>0\mu_{n} > 0 при n>0n > 0 и μ0=0\mu_{0} = 0. Пусть X(0)=n>0X(0) = n > 0, и пусть DnD_{n} — время до первого момента, когда процесс примет значение n−1n-1.

?
(a)

Покажите, что dn=E[Dn]d_{n} = \mathbb {E}\left[D_{n}\right] удовлетворяет

λndn+1=μndn−1,n≥1 \lambda _{n} d_{n+1} = \mu _{n} d_{n}-1, \quad n \geq 1
(b)

Покажите, что производящая функция моментов Mn(θ)=E[eθDn]M_{n}(\theta ) = \mathbb {E}\left[e^{\theta D_{n}}\right] удовлетворяет

(λn+μn−θ)Mn(θ)=μn+λnMn(θ)Mn+1(θ),n≥1 \left(\lambda _{n}+\mu _{n}-\theta \right) M_{n}(\theta ) = \mu _{n}+\lambda _{n} M_{n}(\theta ) M_{n+1}(\theta ), \quad n \geq 1
Задача 6.11.9

В модели популяции биоплёнки предположим, что имеется nn доступных для колонизации «ниш» (или «источников пищи»). Пусть X(t)X(t) — число занятых ниш в момент времени tt, и предположим, что XX — марковская цепь, эволюционирующая следующим образом. Время жизни любой колонии имеет показательное распределение с параметром μ\mu; если X(t)=iX(t) = i, то интенсивность образования новой колонии в пустой нише равна λi(n−i)\lambda i(n-i). Можно предполагать обычную независимость.

При X(0)=1,2,…X(0) = 1,2, \ldots найдите среднее время до вымирания популяции, то есть до момента, когда ни одна ниша не занята. Обсудите последствия для случая большого nn.

?
§
Задача 6.12.1

Клиенты, заходящие в магазин, обслуживаются единственным продавцом в порядке их прибытия. Они прибывают в соответствии с пуассоновским процессом с интенсивностью λ\lambda, а времена их обслуживания — независимые показательно распределённые случайные величины с параметром μ\mu. Рассматривая цепь скачков, покажите, что ожидаемая продолжительность периода занятости BB продавца равна (μ−λ)−1(\mu -\lambda )^{-1} при λ<μ\lambda < \mu. (Период занятости длится с момента, когда клиент прибывает и застаёт продавца свободным, до ближайшего последующего момента, когда продавец снова свободен.)

?
Задача 6.12.2

Иммигранты прибывают в моменты пуассоновского процесса интенсивности ν\nu, и каждый независимо основывает простой процесс рождения интенсивности λ\lambda. В моменты независимого пуассоновского процесса интенсивности δ\delta популяция полностью уничтожается. Найдите производящую функцию вероятностей популяции X(t)X(t) при условии X(0)=0X(0) = 0.

?
Задача 6.12.3

В рамках упражнения (6.12.2) предположим, что каждый иммигрант порождает простой процесс рождения и гибели с интенсивностями λ\lambda и μ\mu. Покажите, что среднее значение размера популяции остаётся ограниченным тогда и только тогда, когда δ>λ−μ\delta > \lambda -\mu.

?
Задача 6.12.4

Очередь M/G/∞M / G / \infty. FTP-сервер принимает клиентов в моменты пуассоновского процесса с параметром λ\lambda, начиная с момента 0. ii-й клиент остаётся подключённым в течение времени SiS_{i}, где SiS_{i} — независимые одинаково распределённые случайные величины, независимые от процесса поступлений. Предполагая, что сервер обладает бесконечной ёмкостью, покажите, что число клиентов, обслуживаемых в момент времени tt, имеет распределение Пуассона с параметром λ∫0t[1−G(x)]dx\lambda \int_{0}^{t}[1-G(x)] d x, где GG — общая функция распределения величин SiS_{i}. Покажите, что среднее этого распределения сходится к λE[S]\lambda \mathbb {E}\left[S\right] при t→∞t \rightarrow \infty.

?
§
Задача 6.13.1

В некотором городе в момент времени t=0t = 0 медведей нет. Бурые медведи и медведи гризли прибывают в соответствии с независимыми пуассоновскими процессами BB и GG с интенсивностями β\beta и γ\gamma соответственно.

?
(a)

Покажите, что первый медведь окажется бурым с вероятностью β/(β+γ)\beta /(\beta +\gamma ).

(b)

Найдите вероятность того, что между двумя последовательными бурыми медведями прибудет ровно rr медведей гризли.

(c)

При условии B(1)=1B(1) = 1 найдите ожидаемое значение момента прибытия первого медведя.

Задача 6.13.2

Пусть Π\Pi — точки неоднородного пуассоновского процесса на Rd\mathbb {R}^{d} с функцией интенсивности λ\lambda. Пусть S=∑x∈Πg(x)S = \sum_{\mathbf{x} \in \Pi } g(\mathbf{x}), где gg — (измеримая) функция, которую мы для удобства считаем неотрицательной.

?
(a)

Покажите непосредственно, что E[S]=∫Rdg(x)λ(x)dx\mathbb {E}\left[S\right] = \int_{\mathbb {R}^{d}} g(\mathbf{x}) \lambda (\mathbf{x}) d \mathbf{x} и Var⁡[S]=∫Rdg(x)2λ(x)dx\operatorname {Var}\left[S\right] = \int_{\mathbb {R}^{d}} g(\mathbf{x})^{2} \lambda (\mathbf{x}) d \mathbf{x}, при условии сходимости этих интегралов.

(b)

Покажите, что

E[e−tS]=exp⁡{−∫Rd(1−e−tg(x))λ(x)dx},t>0 \mathbb {E}\left[e^{-t S}\right] = \exp \left\{ -\int _{\mathbb {R}^{d}}\left(1-e^{-t g(\mathbf{x})}\right) \lambda (\mathbf{x}) d \mathbf{x}\right\} , \quad t > 0

и выведите отсюда, что P(S<∞)=1\mathbb {P}\left(S < \infty \right) = 1, если ∫Rdmin⁡{1,g(x)}λ(x)dx<∞\int_{\mathbb {R}^{d}} \min \left\{ 1, g(\mathbf{x})\right\} \lambda (\mathbf{x}) d \mathbf{x} < \infty.

(c)

Если выполнено интегральное условие пункта (b), покажите, что характеристическая функция ϕ\phi величины SS удовлетворяет

ϕ(t)=exp⁡{−∫Rd(1−eitg(x))λ(x)dx},t∈R \phi (t) = \exp \left\{ -\int _{\mathbb {R}^{d}}\left(1-e^{i t g(\mathbf{x})}\right) \lambda (\mathbf{x}) d \mathbf{x}\right\} , \quad t \in \mathbb {R}
Задача 6.13.3

Пусть Π\Pi — пуассоновский процесс с постоянной интенсивностью λ\lambda на поверхности сферы радиуса 1 в R3\mathbb {R}^{3}. Пусть PP — процесс, задаваемый координатами (X,Y)(X, Y) точек, спроецированных на плоскость, проходящую через центр сферы. Покажите, что PP — пуассоновский процесс, и найдите его функцию интенсивности.

?
Задача 6.13.4

Повторите упражнение (6.13.3) для случая, когда Π\Pi — однородный пуассоновский процесс на шаре {(x1,x2,x3):x12+x22+x32≤1}\left\{ \left(x_{1}, x_{2}, x_{3}\right) : x_{1}^{2}+x_{2}^{2}+x_{3}^{2} \leq 1\right\}.

?
Задача 6.13.5

Вы втыкаете булавки в карту Земли в проекции Меркатора в соответствии с пуассоновским процессом постоянной интенсивности λ\lambda. Какова функция интенсивности соответствующего процесса на глобусе? Какой была бы функция интенсивности на карте, если бы вы образовали пуассоновский процесс постоянной интенсивности λ\lambda падений метеоритов на поверхность Земли?

?
Задача 6.13.6

rr-я точка TrT_{r} пуассоновского процесса NN постоянной интенсивности λ\lambda на R+\mathbb {R}_{+}порождает эффект Xre−α(t−Tr)X_{r} e^{-\alpha \left(t-T_{r}\right)} в момент времени t≥Trt \geq T_{r}, где XrX_{r} независимы, одинаково распределены и имеют конечную дисперсию. Найдите среднее и дисперсию суммарного эффекта S(t)=∑r=1N(t)Xre−α(t−Tr)S(t) = \sum_{r = 1}^{N(t)} X_{r} e^{-\alpha \left(t-T_{r}\right)} через первые два момента XrX_{r}, и вычислите Cov⁡[S(s),S(t)]\operatorname {Cov}\left[S(s), S(t)\right].

Каково поведение корреляции ρ(S(s),S(t))\rho (S(s), S(t)) при s→∞s \rightarrow \infty с фиксированным t−st-s?

?
Задача 6.13.7

Пусть NN — неоднородный пуассоновский процесс на R+\mathbb {R}_{+} с функцией интенсивности λ\lambda. Найдите совместную плотность первых двух интервалов между событиями и выведите отсюда, что в общем случае они не являются независимыми.

?
Задача 6.13.8

Пусть {Nr(t):r≥1}\left\{ N_{r}(t): r \geq 1\right\} — семейство независимых пуассоновских процессов на R+\mathbb {R}_{+} с постоянными интенсивностями {λr:r≥1}\left\{ \lambda_{r}: r \geq 1\right\} соответственно, такое что ∑rλr=λ<∞\sum_{r} \lambda_{r} = \lambda < \infty. Положим N(t)=∑rNr(t)N(t) = \sum_{r} N_{r}(t), и пусть II обозначает индекс процесса, дающего первую точку в NN, наступающую в момент времени TT. Покажите, что

P(I=i,T≥t)=P(I=i)P(T≥t)=λiλe−λt,i≥1 \mathbb {P}\left(I = i, T \geq t\right) = \mathbb {P}\left(I = i\right) \mathbb {P}\left(T \geq t\right) = \frac{\lambda _{i}}{\lambda } e^{-\lambda t}, \quad i \geq 1
?
Задача 6.13.9

Пусть Π\Pi — пуассоновский процесс на R2\0\mathbb {R}^{2} \backslash \mathbf{0} с функцией интенсивности λ(u,v)=(u2+v2)−3/2\lambda (u, v) = \left(u^{2}+v^{2}\right)^{-3 / 2}. Каждая точка (U,V)∈Π(U, V) \in \Pi порождает прямую Ux+Vy=1U x+V y = 1 на плоскости x/y,Lx / y, \mathbb {L}.

?
(a)

Пусть (p,θ)(p, \theta ) — полярные координаты основания перпендикуляра, опущенного из начала координат на прямую ux+vy=1u x+v y = 1 на плоскости x/y,Lx / y, \mathbb {L}. Выразите (p,θ)(p, \theta ) через (u,v)(u, v).

(b)

Покажите, что отображение из пункта (a) переводит пуассоновский процесс прямых в равномерный пуассоновский процесс на полосе S⊥=R×[0,π)S^{\perp } = \mathbb {R} \times [0, \pi ).

(c)

Покажите, что процесс прямых на L\mathbb {L} инвариантен относительно сдвигов и вращений L\mathbb {L}.

Задача 6.13.10

Большой континент пересекает дважды бесконечное прямое шоссе, вдоль которого грузовики припаркованы в точках пуассоновского процесса с постоянной интенсивностью 1. Массы грузовиков — независимые одинаково распределённые случайные величины, независимые от мест парковки. Пусть GG — гравитационное притяжение, оказываемое грузовиками на пешехода единичной массы, стоящего рядом с шоссе. Гравитационную постоянную можно считать равной 1.

Покажите, что GG имеет характеристическую функцию вида ϕ(t)=exp⁡(−c∣t∣1/2)\phi (t) = \exp \left(-c|t|^{1 / 2}\right), где c>0c > 0. Выразите cc через среднее типичной массы MM.

?
§
Задача 6.14.1

Пусть P\mathbf{P} — стохастическая матрица на конечном множестве Θ\Theta со стационарным распределением π\mathbf{\pi }. Определим скалярное произведение ⟨x,y⟩=∑k∈Θxkykπk\langle \mathbf{x}, \mathbf{y}\rangle = \sum_{k \in \Theta } x_{k} y_{k} \pi_{k} и пусть l2(π)={x∈RΘ:⟨x,x⟩<∞}l^{2}(\pi ) = \left\{ \mathbf{x} \in \mathbb {R}^{\Theta }:\langle \mathbf{x}, \mathbf{x}\rangle < \infty \right\}. Покажите, в очевидных обозначениях, что P\mathbf{P} обратима относительно π\pi тогда и только тогда, когда ⟨x,Py⟩=⟨Px,y⟩\langle \mathbf{x}, \mathbf{P y}\rangle = \langle \mathbf{P x}, \mathbf{y}\rangle для всех x,y∈l2(π)\mathbf{x}, \mathbf{y} \in l^{2}(\pi ).

?
Задача 6.14.2

Покажите, что возможным выбором вероятностей принятия в общем алгоритме Гастингса являются

bij=πjgjiπigij+πjgji b_{i j} = \frac{\pi _{j} g_{j i}}{\pi _{i} g_{i j}+\pi _{j} g_{j i}}

где G=(gij)\mathbf{G} = \left(g_{i j}\right) — матрица предложений.

?
Задача 6.14.3

Пусть SS — счётное множество. Для каждого j∈Sj \in S множества Ajk,k∈SA_{j k}, k \in S, образуют разбиение интервала [0,1][0,1]. Пусть g:S×[0,1]→Sg: S \times [0,1] \rightarrow S задана как g(j,u)=kg(j, u) = k, если u∈Ajku \in A_{j k}. Последовательность случайных величин {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} строится рекурсивно как Xn+1=g(Xn,Un+1),n≥0X_{n+1} = g\left(X_{n}, U_{n+1}\right), n \geq 0, где {Un:n≥1}\left\{ U_{n}: n \geq 1\right\} — независимые случайные величины, равномерно распределённые на [0,1][0,1]. Покажите, что XX — марковская цепь, и найдите её матрицу переходных вероятностей.

?
Задача 6.14.4

Пусть U=(ust)\mathbf{U} = \left(u_{s t}\right) — конечная стохастическая матрица размера ∣S∣×∣T∣\left|S\right| \times \left|T\right|. Эргодическим коэффициентом Добрушина называется

d(U)=12sup⁡i,j∈S∑t∈T∣uit−ujt∣ d(\mathbf{U}) = \frac{1}{2} \sup _{i, j \in S} \sum _{t \in T}\left|u_{i t}-u_{j t}\right|
?
(a)

Покажите, что если V\mathbf{V} — конечная стохастическая матрица размера ∣T∣×∣U∣\left|T\right| \times \left|U\right|, то d(UV)≤d(U)d(V)d(\mathbf{U V}) \leq d(\mathbf{U}) d(\mathbf{V}).

(b)

Пусть XX и YY — дискретные марковские цепи с одной и той же матрицей переходных вероятностей P\mathbf{P}, и покажите, что

∑k∣P(Xn=k)−P(Yn=k)∣≤d(P)n∑k∣P(X0=k)−P(Y0=k)∣ \sum _{k}\left|\mathbb {P}\left(X_{n} = k\right)-\mathbb {P}\left(Y_{n} = k\right)\right| \leq d(\mathbf{P})^{n} \sum _{k}\left|\mathbb {P}\left(X_{0} = k\right)-\mathbb {P}\left(Y_{0} = k\right)\right|
Задача 6.14.5

Пусть π\mathbf{\pi } — положительная функция вероятностей на конечном множестве Θ\Theta, и пусть P\mathbf{P} — матрица переходных вероятностей неприводимой апериодической марковской цепи со стационарным распределением π\mathbf{\pi }. Пусть W=(W(i):i∈Θ)W = (W(i): i \in \Theta ) — вектор случайных величин, такой что P(W(i)=j)=pij\mathbb {P}\left(W(i\right) = j) = p_{i j} при i,j∈Θi, j \in \Theta, и используем WW в качестве правила обновления в алгоритме «склейки из прошлого» для выборки из π\pi.

?
(a)

Если W(i),i∈ΘW(i), i \in \Theta, независимы, покажите, что время слияния почти наверное конечно.

(b)

Приведите два примера ситуаций, в которых время слияния почти наверное бесконечно.

Задача 6.14.6

Покажите, что распределение Изинга из примера (6.14.2) удовлетворяет решёточному условию FKG (6.14.20).

?
§
Задача 6.15.1

Классифицируйте состояния дискретных марковских цепей с пространством состояний S={1,2,3,4}S = \left\{ 1,2,3,4\right\} и следующими матрицами переходных вероятностей:

?
(a)
[13230012120014014120001] \left[\begin{smallmatrix} \frac{1}{3} & \frac{2}{3} & 0 & 0 \\ \frac{1}{2} & \frac{1}{2} & 0 & 0 \\ \frac{1}{4} & 0 & \frac{1}{4} & \frac{1}{2} \\ 0 & 0 & 0 & 1 \end{smallmatrix}\right]

Вычислите f34(n)f_{34}(n) и выведите отсюда, что вероятность окончательного поглощения в состоянии 44, начиная с 33, равна 23\frac{2}{3}.

(b)
[01212013002310000010] \left[\begin{smallmatrix} 0 & \frac{1}{2} & \frac{1}{2} & 0 \\ \frac{1}{3} & 0 & 0 & \frac{2}{3} \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 \end{smallmatrix}\right]

Найдите средние времена возврата состояний.

Задача 6.15.2

Матрица переходных вероятностей называется дважды стохастической, если суммы всех её столбцов равны 11, то есть если ∑ipij=1\sum_{i} p_{i j} = 1 для всех j∈Sj \in S.

?
(a)

Покажите, что если конечная цепь имеет дважды стохастическую матрицу переходных вероятностей, то все её состояния положительно возвратны, и что если, кроме того, цепь неприводима и апериодична, то pij(n)→N−1p_{i j}(n) \rightarrow N^{-1} при n→∞n \rightarrow \infty, где NN — число состояний.

(b)

Покажите, что если бесконечная неприводимая цепь имеет дважды стохастическую матрицу переходных вероятностей, то её состояния либо все нуль-возвратны, либо все невозвратны.

Задача 6.15.3

Докажите, что сообщающиеся между собой состояния марковской цепи имеют одинаковый период.

?
Задача 6.15.4
?
(a)

Покажите, что для каждой пары состояний i,ji, j неприводимой апериодической цепи существует N=N(i,j)N = N(i, j), такое что pij(n)>0p_{i j}(n) > 0 при всех n≥Nn \geq N.

(b)

Пусть XX и YY — независимые неприводимые апериодические цепи с одним и тем же пространством состояний SS и матрицей переходных вероятностей P\mathbf{P}. Покажите, что двумерная цепь Zn=(Xn,Yn),n≥0Z_{n} = \left(X_{n}, Y_{n}\right), n \geq 0, неприводима и апериодична.

(c)

Покажите, что двумерная цепь ZZ может быть приводимой, если XX и YY периодичны.

Задача 6.15.5

Предположим, что {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} — дискретная марковская цепь с X0=iX_{0} = i. Пусть NN — общее число последующих посещений цепью состояния jj. Покажите, что

P(N=n)={1−fij если n=0fij(fjj)n−1(1−fjj) если n≥1 \mathbb {P}\left(N = n\right) = \begin{cases} 1-f_{i j} & \text{ если } n = 0 \\ f_{i j}\left(f_{j j}\right)^{n-1}\left(1-f_{j j}\right) & \text{ если } n \geq 1\end{cases}

и выведите отсюда, что P(N=∞)=1\mathbb {P}\left(N = \infty \right) = 1 тогда и только тогда, когда fij=fjj=1f_{i j} = f_{j j} = 1.

?
Задача 6.15.6

Пусть ii и jj — два состояния дискретной марковской цепи. Покажите, что если ii сообщается с jj, то существует положительная вероятность достичь jj из ii, ни разу не вернувшись в ii по пути. Выведите отсюда, что если цепь неприводима и возвратна, то вероятность fijf_{i j} когда-либо достичь jj из ii равна 1 для всех ii и jj.

?
Задача 6.15.7

Пусть {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} — возвратная неприводимая марковская цепь на пространстве состояний SS с матрицей переходных вероятностей P\mathbf{P}, и пусть x\mathbf{x} — положительное решение уравнения x=xP\mathbf{x} = \mathbf{x P}.

?
(a)

Покажите, что

qij(n)=xjxipji(n),i,j∈S,n≥1 q_{i j}(n) = \frac{x_{j}}{x_{i}} p_{j i}(n), \quad i, j \in S, n \geq 1

задаёт nn-шаговые переходные вероятности возвратной неприводимой марковской цепи на SS, вероятности первого перехода которой задаются как

gij(n)=xjxilji(n),i≠j,n≥1 g_{i j}(n) = \frac{x_{j}}{x_{i}} l_{j i}(n), \quad i \neq j, n \geq 1

где lji(n)=Pj(Xn=i,T>n)l_{j i}(n) = \mathbb {P}_{j}\left(X_{n} = i, T > n\right) и T=min⁡{m>0:Xm=j}T = \min \left\{ m > 0: X_{m} = j\right\}.

(b)

Покажите, что x\mathbf{x} единственно с точностью до мультипликативной постоянной.

(c)

Пусть Tj=min⁡{n≥1:Xn=j}T_{j} = \min \left\{ n \geq 1: X_{n} = j\right\}, и определим hij=Pi(Tj≤Ti)h_{i j} = \mathbb {P}_{i}\left(T_{j} \leq T_{i}\right). Покажите, что xihij=xjhjix_{i} h_{i j} = x_{j} h_{j i} для всех i,j∈Si, j \in S.

Задача 6.15.8

Последовательность u={un:n≥0}u = \left\{ u_{n}: n \geq 0\right\} называется «последовательностью восстановления», если

u0=1,un=∑i=1nfiun−i при n≥1 u_{0} = 1, \quad u_{n} = \sum _{i = 1}^{n} f_{i} u_{n-i} \quad \text{ при } n \geq 1

для некоторого семейства f={fn:n≥1}f = \left\{ f_{n}: n \geq 1\right\} неотрицательных чисел, суммирующихся в 1.

?
(a)

Покажите, что uu является последовательностью восстановления тогда и только тогда, когда существует марковская цепь XX на счётном пространстве состояний SS, такая что un=P(Xn=s∣X0=s)u_{n} = \mathbb {P}\left(X_{n} = s \mid X_{0} = s\right) для некоторого возвратного s∈Ss \in S и всех n≥1n \geq 1.

(b)

Покажите, что если uu и vv — последовательности восстановления, то таковой является и {unvn:n≥0}\left\{ u_{n} v_{n}: n \geq 0\right\}.

Задача 6.15.9

Рассмотрим симметричное случайное блуждание в трёх измерениях по множеству точек {(x,y,z):x,y,z=0,±1,±2,…}\left\{ (x, y, z): x, y, z = 0, \pm 1, \pm 2, \ldots \right\}; этот процесс представляет собой последовательность точек {Xn:n≥0}\left\{ \mathbf{X}_{n}: n \geq 0\right\}, такую что P(Xn+1=Xn+ϵ)=16\mathbb {P}\left(\mathbf{X}_{n+1} = \mathbf{X}_{n}+\mathbf{\epsilon }\right) = \frac{1}{6} для ϵ=(±1,0,0),(0,±1,0),(0,0,±1)\mathbf{\epsilon } = (\pm 1,0,0),(0, \pm 1,0),(0,0, \pm 1). Предположим, что X0=(0,0,0)\mathbf{X}_{0} = (0,0,0). Покажите, что

P(X2n=(0,0,0))=(16)2n∑i+j+k=n(2n)!(i!j!k!)2=(12)2n(2nn)∑i+j+k=n(n!3ni!j!k!)2 \mathbb {P}\left(\mathbf{X}_{2 n} = (0,0,0)\right) = \left(\frac{1}{6}\right)^{2 n} \sum _{i+j+k = n} \frac{(2 n)!}{(i!j!k!)^{2}} = \left(\frac{1}{2}\right)^{2 n}\binom {2 n}{n} \sum _{i+j+k = n}\left(\frac{n!}{3^{n} i!j!k!}\right)^{2}

и с помощью формулы Стирлинга выведите отсюда, что начало координат — невозвратное состояние.

?
Задача 6.15.10

Рассмотрим трёхмерную версию модели рака (6.12.12). Если κ=1\kappa = 1, неизбежны ли в этом случае империи из теоремы (6.12.14)?

?
Задача 6.15.11

Пусть XX — дискретная марковская цепь с пространством состояний S={1,2}S = \left\{ 1,2\right\} и матрицей переходных вероятностей

P=[1−ααβ1−β] \mathbf{P} =\left[\begin{smallmatrix} 1-\alpha & \alpha \\ \beta & 1-\beta \end{smallmatrix}\right]

Классифицируйте состояния цепи. Предположим, что αβ>0\alpha \beta > 0 и αβ≠1\alpha \beta \neq 1. Найдите nn-шаговые переходные вероятности и покажите непосредственно, что они сходятся к единственному стационарному распределению при n→∞n \rightarrow \infty. При каких значениях α\alpha и β\beta цепь обратима в равновесии?

?
Задача 6.15.12

NN чёрных шаров и NN белых шаров размещаются в двух урнах так, что каждая содержит NN шаров. После каждой единицы времени из каждой урны наугад выбирается по одному шару, и эти два выбранных шара меняются местами. Пусть состоянием системы обозначается число чёрных шаров в первой урне. Запишите матрицу переходных вероятностей этой марковской цепи и найдите единственное стационарное распределение. Обратима ли цепь в равновесии?

?
Задача 6.15.13

Рассмотрим марковскую цепь на множестве S={0,1,2,…}S = \left\{ 0,1,2, \ldots \right\} с переходными вероятностями pi,i+1=aip_{i, i+1} = a_{i}, pi,0=1−ai,i≥0p_{i, 0} = 1-a_{i}, i \geq 0, где (ai:i≥0)\left(a_{i}: i \geq 0\right) — последовательность постоянных, удовлетворяющих 0<ai<10 < a_{i} < 1 при всех ii. Пусть b0=1,bi=a0a1⋯ai−1b_{0} = 1, b_{i} = a_{0} a_{1} \cdots a_{i-1} при i≥1i \geq 1. Покажите, что цепь

?
(a)

возвратна тогда и только тогда, когда bi→0b_{i} \rightarrow 0 при i→∞i \rightarrow \infty,

(b)

положительно возвратна тогда и только тогда, когда ∑ibi<∞\sum_{i} b_{i} < \infty, и запишите стационарное распределение, если последнее условие выполнено. Пусть AA и β\beta — положительные постоянные, и предположим, что ai=1−Ai−βa_{i} = 1-A i^{-\beta } при всех достаточно больших ii. Покажите, что цепь

(c)

невозвратна, если β>1\beta > 1,

(d)

положительно возвратна, если β<1\beta < 1.

Наконец, если β=1\beta = 1, покажите, что цепь

(e)

положительно возвратна, если A>1A > 1,

(f)

нуль-возвратна, если A≤1A \leq 1.

Задача 6.15.14

Пусть XX — марковская цепь с непрерывным временем, счётным пространством состояний SS и полугруппой {Pt}\left\{ \mathbf{P}_{t}\right\}. Покажите, что pij(t)p_{i j}(t) — непрерывная функция tt. Пусть g(t)=−log⁡pii(t)g(t) = -\log p_{i i}(t); покажите, что gg — непрерывная функция, g(0)=0g(0) = 0, и g(s+t)≤g(s)+g(t)g(s+t) \leq g(s)+g(t). Говорят, что gg «субаддитивна», и хорошо известная теорема даёт результат, что

lim⁡t↓0g(t)t=λ существует и λ=sup⁡t>0g(t)t≤∞ \lim _{t \downarrow 0} \frac{g(t)}{t} = \lambda \quad \text{ существует и } \quad \lambda = \sup _{t > 0} \frac{g(t)}{t} \leq \infty

Выведите отсюда, что предел gii=lim⁡t↓0t−1{pii(t)−1}g_{i i} = \lim_{t \downarrow 0} t^{-1}\left\{ p_{i i}(t)-1\right\} существует.

?
Задача 6.15.15

Пусть XX — марковская цепь с непрерывным временем и генератором G=(gij)\mathbf{G} = \left(g_{i j}\right). Покажите, что XX неприводима тогда и только тогда, когда для любой пары различных состояний i,ji, j существует последовательность различных состояний i,k1,k2,…,kn,ji, k_{1}, k_{2}, \ldots , k_{n}, j, такая что gi,k1gk1,k2⋯gkn,j>0g_{i, k_{1}} g_{k_{1}, k_{2}} \cdots g_{k_{n}, j} > 0.

?
Задача 6.15.16
?
(a)

Пусть T>0T > 0, и пусть X={X(t):0≤t≤T}X = \left\{ X(t): 0 \leq t \leq T\right\} — неприводимая, невзрывающаяся марковская цепь со стационарным распределением π\pi, и предположим, что X(0)X(0) имеет распределение π\pi. Пусть Y(t)=X(T−t)Y(t) = X(T-t) при 0≤t≤T0 \leq t \leq T. Мы называем XX обратимой (в равновесии), если XX и YY имеют одинаковые совместные распределения.

(i) Покажите, что YY — (непрерывная слева) марковская цепь с переходными вероятностями p^ij(t)=(πj/πi)pji(t)\widehat{p}_{i j}(t) = \left(\pi_{j} / \pi_{i}\right) p_{j i}(t) и генератором G^\widehat{\mathbf{G}}, удовлетворяющим πjg^ji=πigij\pi_{j} \hat{g}_{j i} = \pi_{i} g_{i j}, где pji(t)p_{j i}(t) и G=(gij)\mathbf{G} = \left(g_{i j}\right) относятся к XX. Покажите, что YY неприводима и невзрывающаяся со стационарным распределением π\pi.

(ii) Покажите, что XX обратима в равновесии тогда и только тогда, когда выполняются уравнения детального баланса πigij=πjgji\pi_{i} g_{i j} = \pi_{j} g_{j i} (для всех ii и jj).

(iii) Покажите, что мера ν\mathbf{\nu } удовлетворяет νG=0\mathbf{\nu } \mathbf{G} = \mathbf{0}, если она удовлетворяет уравнениям детального баланса.

(b)

Пусть XX неприводима и невзрывающаяся со стационарным распределением π\pi, и предположим, что X(0)X(0) имеет распределение π\pi.

(i) Критерий Колмогорова. Покажите, что XX обратима тогда и только тогда, когда для всех nn и всех конечных последовательностей состояний k1,k2,…,knk_{1}, k_{2}, \ldots , k_{n}

gk1,k2gk2,k3⋯gkn−1,kngkn,k1=gk1,kngkn,kn−1⋯gk2,k1 g_{k_{1}, k_{2}} g_{k_{2}, k_{3}} \cdots g_{k_{n-1}, k_{n}} g_{k_{n}, k_{1}} = g_{k_{1}, k_{n}} g_{k_{n}, k_{n-1}} \cdots g_{k_{2}, k_{1}}

(ii) Критерий Келли. Покажите, что XX обратима, если для всех различных троек i,j,k∈Si, j, k \in S выполнено gijgjkgki=gikgkjgjig_{i j} g_{j k} g_{k i} = g_{i k} g_{k j} g_{j i}, и, кроме того, существует c∈Sc \in S, такое что gic>0g_{i c} > 0 для всех i≠ci \neq c.

(c)

Покажите, что любая неприводимая цепь XX ровно с двумя состояниями обратима в равновесии.

(d)

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

Задача 6.15.17

Покажите, что не всякую дискретную марковскую цепь можно вложить в цепь с непрерывным временем. Точнее, пусть

P=[α1−α1−αα] при некотором 0<α<1 \mathbf{P} =\left[\begin{smallmatrix} \alpha & 1-\alpha \\ 1-\alpha & \alpha \end{smallmatrix}\right] \quad \text{ при некотором } 0 < \alpha < 1

— матрица переходных вероятностей. Покажите, что полугруппа {Pt}\left\{ \mathbf{P}_{t}\right\} переходных вероятностей в непрерывном времени, такая что P1=P\mathbf{P}_{1} = \mathbf{P}, существует тогда и только тогда, когда 12<α<1\frac{1}{2} < \alpha < 1. В этом случае покажите, что {Pt}\left\{ \mathbf{P}_{t}\right\} единственна, и вычислите её через α\alpha.

?
Задача 6.15.18

Рассмотрим процесс иммиграции-гибели X(t)X(t) — процесс рождения и гибели с интенсивностями λn=λ\lambda_{n} = \lambda, μn=nμ\mu_{n} = n \mu. Покажите, что его производящая функция G(s,t)=E[sX(t)]G(s, t) = \mathbb {E}\left[s^{X(t)}\right] задаётся формулой

G(s,t)={1+(s−1)e−μt}Iexp⁡{ρ(s−1)(1−e−μt)} G(s, t) = \left\{ 1+(s-1) e^{-\mu t}\right\} ^{I} \exp \left\{ \rho (s-1)\left(1-e^{-\mu t}\right)\right\}

где ρ=λ/μ\rho = \lambda / \mu и X(0)=IX(0) = I. Выведите отсюда предельное распределение X(t)X(t) при t→∞t \rightarrow \infty.

?
Задача 6.15.19

Пусть NN — неоднородный пуассоновский процесс на R+=[0,∞)\mathbb {R}_{+} = [0, \infty ) с функцией интенсивности λ\lambda.

?
(a)

Запишите прямые и обратные уравнения для NN и решите их.

(b)

Пусть N(0)=0N(0) = 0; найдите плотность времени TT до первого поступления в процессе. Если λ(t)=c/(1+t)\lambda (t) = c /(1+t), покажите, что E[T]<∞\mathbb {E}\left[T\right] < \infty тогда и только тогда, когда c>1c > 1.

Задача 6.15.20

Последовательные предложения за мой дом — независимые одинаково распределённые случайные величины X1X_{1}, X2,…X_{2}, \ldots с плотностью ff и функцией распределения FF. Пусть Y1=X1Y_{1} = X_{1}, пусть Y2Y_{2} — первое предложение, превышающее Y1Y_{1}, и вообще пусть Yn+1Y_{n+1} — первое предложение, превышающее YnY_{n}. Покажите, что Y1,Y2,…Y_{1}, Y_{2}, \ldots являются моментами поступлений в неоднородном пуассоновском процессе с функцией интенсивности λ(t)=f(t)/(1−F(t))\lambda (t) = f(t) /(1-F(t)). Величины YiY_{i} называются «рекордными значениями».

Теперь пусть Z1Z_{1} — первое полученное предложение, являющееся на данный момент вторым по величине, и пусть Z2Z_{2} — второе такое предложение, и так далее. Покажите, что ZiZ_{i} являются моментами поступлений неоднородного пуассоновского процесса с функцией интенсивности λ\lambda.

?
Задача 6.15.21

Пусть NN — пуассоновский процесс с постоянной интенсивностью λ\lambda, и пусть Y1,Y2,…Y_{1}, Y_{2}, \ldots — независимые случайные величины с общей характеристической функцией ϕ\phi и плотностью ff. Процесс N∗(t)=Y1+Y2+⋯+YN(t)N^{*}(t) = Y_{1}+Y_{2}+\cdots +Y_{N(t)} называется сложным пуассоновским процессом. YnY_{n} — изменение значения N∗N^{*} при nn-м поступлении пуассоновского процесса NN. Представьте это так. «Случайный будильник» звонит в моменты поступлений пуассоновского процесса. При nn-м звонке процесс N∗N^{*} накапливает дополнительную величину YnY_{n}. Запишите прямое уравнение для N∗N^{*} и с его помощью найдите характеристическую функцию N∗(t)N^{*}(t). Видите ли вы непосредственно, почему она имеет найденную вами форму?

?
Задача 6.15.22

Если функция интенсивности λ\lambda неоднородного пуассоновского процесса NN сама является случайным процессом, то NN называется дважды стохастическим пуассоновским процессом (или процессом Кокса).

?
(a)

Рассмотрим случай, когда λ(t)=Λ\lambda (t) = \Lambda при всех tt, а Λ\Lambda — случайная величина, принимающая одно из двух значений λ1\lambda_{1} или λ2\lambda_{2}, каждое с равной вероятностью 12\frac{1}{2}. Найдите производящую функцию вероятностей N(t)N(t) и выведите отсюда её среднее и дисперсию.

(b)

Для дважды стохастического пуассоновского процесса NN покажите, что Var⁡[N(t])≥E[N(t])\operatorname {Var}\left[N(t\right]) \geq \mathbb {E}\left[N(t\right]).

(c)

Пусть MM — обычный пуассоновский процесс на временном интервале [0,∞)[0, \infty ) с постоянной интенсивностью 1. Пусть M∗M^{*} получен из MM удалением kk-го поступления для каждого нечётного значения kk. Является ли M∗M^{*}: (i) пуассоновским процессом, или (ii) дважды стохастическим пуассоновским процессом?

Задача 6.15.23

Покажите, что простой процесс рождения XX с параметром λ\lambda является дважды стохастическим пуассоновским процессом с функцией интенсивности λ(t)=λX(t)\lambda (t) = \lambda X(t).

?
Задача 6.15.24

Марковская цепь X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} — это процесс рождения, интенсивности λk(t)\lambda_{k}(t) которого зависят также от времени tt и задаются как

P(X(t+h)=k+1∣X(t)=k)=1+μk1+μth+o(h) \mathbb {P}\left(X(t+h\right) = k+1 \mid X(t) = k) = \frac{1+\mu k}{1+\mu t} h+\mathrm{o}(h)

при h↓0h \downarrow 0. Покажите, что производящая функция вероятностей G(s,t)=E[sX(t)]G(s, t) = \mathbb {E}\left[s^{X(t)}\right] удовлетворяет

∂G∂t=s−11+μt{G+μs∂G∂s},0<s<1 \frac{\partial G}{\partial t} = \frac{s-1}{1+\mu t}\left\{ G+\mu s \frac{\partial G}{\partial s}\right\} , \quad 0 < s < 1

Отсюда найдите среднее и дисперсию X(t)X(t), когда X(0)=IX(0) = I.

?
Задача 6.15.25
?
(a)

Пусть XX — процесс рождения и гибели со строго положительными интенсивностями рождения λ0,λ1,…\lambda_{0}, \lambda_{1}, \ldots и интенсивностями гибели μ1,μ2,…\mu_{1}, \mu_{2}, \ldots Пусть ηi\eta_{i} — вероятность того, что X(t)X(t) когда-либо примет значение 0, начиная с X(0)=iX(0) = i. Покажите, что

λjηj+1−(λj+μj)ηj+μjηj−1=0,j≥1 \lambda _{j} \eta _{j+1}-\left(\lambda _{j}+\mu _{j}\right) \eta _{j}+\mu _{j} \eta _{j-1} = 0, \quad j \geq 1

и выведите отсюда, что ηi=1\eta_{i} = 1 для всех ii, если ∑1∞ej=∞\sum_{1}^{\infty } e_{j} = \infty, где ej=μ1μ2⋯μj/(λ1λ2⋯λj)e_{j} = \mu_{1} \mu_{2} \cdots \mu_{j} /\left(\lambda_{1} \lambda_{2} \cdots \lambda_{j}\right).

(b)

Для дискретной цепи на неотрицательных целых числах с

pj,j+1=(j+1)2j2+(j+1)2 и pj,j−1=j2j2+(j+1)2 p_{j, j+1} = \frac{(j+1)^{2}}{j^{2}+(j+1)^{2}} \quad \text{ и } \quad p_{j, j-1} = \frac{j^{2}}{j^{2}+(j+1)^{2}}

найдите вероятность того, что цепь когда-либо посетит 00, начиная с 1.

Задача 6.15.26

Найдите хорошее необходимое условие и хорошее достаточное условие для того, чтобы процесс рождения и гибели XX из задачи (6.15.25а) был честным.

?
Задача 6.15.27

Пусть XX — простой симметричный процесс рождения и гибели с λn=μn=nλ\lambda_{n} = \mu_{n} = n \lambda, и пусть TT — время до вырождения. Покажите, что

P(T≤x∣X(0)=I)=(λx1+λx)I \mathbb {P}\left(T \leq x \mid X(0\right) = I) = \left(\frac{\lambda x}{1+\lambda x}\right)^{I}

и выведите отсюда, что вырождение достоверно, если P(X(0)<∞)=1\mathbb {P}\left(X(0\right) < \infty ) = 1. Покажите, что P(λT/I≤x∣X(0)=I)→e−1/x\mathbb {P}\left(\lambda T / I \leq x \mid X(0\right) = I) \rightarrow e^{-1 / x} при I→∞I \rightarrow \infty.

?
Задача 6.15.28

Пусть XX — процесс иммиграции-гибели-катастроф, то есть процесс рождения и гибели с параметрами λi=λ,μi=iμ\lambda_{i} = \lambda , \mu_{i} = i \mu, с дополнительной возможностью «катастроф», сводящих популяцию к 0. Катастрофы происходят в моменты пуассоновского процесса интенсивности δ\delta, независимо от всех предшествующих рождений и гибелей.

?
(a)

Покажите, что XX обладает стационарным распределением, и найдите выражение для производящей функции этого распределения.

(b)

Покажите, что в равновесии среднее X(t)X(t) равно λ/(δ+μ)\lambda /(\delta +\mu ).

Задача 6.15.29

С каждым достаточно «хорошим» (скажем, измеримым по Лебегу) подмножеством BB вещественной прямой R\mathbb {R} связана случайная величина X(B)X(B), такая что

(a) X(B)X(B) принимает значения в {0,1,2,…}\left\{ 0,1,2, \ldots \right\},

(b) если B1,B2,…,BnB_{1}, B_{2}, \ldots , B_{n} не пересекаются, то X(B1),X(B2),…,X(Bn)X\left(B_{1}\right), X\left(B_{2}\right), \ldots , X\left(B_{n}\right) независимы, и, кроме того, X(B1∪B2)=X(B1)+X(B2)X\left(B_{1} \cup B_{2}\right) = X\left(B_{1}\right)+X\left(B_{2}\right)

(c) распределение X(B)X(B) зависит от BB только через её меру Лебега («длину») ∣B∣\left|B\right|, и

P(X(B)≥1)P(X(B)=1)→1 при ∣B∣→0 \frac{\mathbb {P}\left(X(B) \geq 1\right)}{\mathbb {P}\left(X(B) = 1\right)} \rightarrow 1 \quad \text{ при }|B| \rightarrow 0

Покажите, что XX — пуассоновский процесс.

?
Задача 6.15.30

Пусть NN — пуассоновский процесс на R2\mathbb {R}^{2} с постоянной интенсивностью λ\lambda, и пусть R(1)<R(2)<R_{(1)} < R_{(2)} <... — упорядоченные расстояния от начала координат до точек процесса.

?
(a)

Покажите, что R(1)2,R(2)2,…R_{(1)}^{2}, R_{(2)}^{2}, \ldots — точки пуассоновского процесса на R+=[0,∞)\mathbb {R}_{+} = [0, \infty ) с интенсивностью λπ\lambda \pi.

(b)

Покажите, что R(k)R_{(k)} имеет плотность

f(r)=2πλr(λπr2)k−1e−λπr2(k−1)!,r>0 f(r) = \frac{2 \pi \lambda r\left(\lambda \pi r^{2}\right)^{k-1} e^{-\lambda \pi r^{2}}}{(k-1)!}, \quad r > 0
Задача 6.15.31

Пусть XX — nn-мерный пуассоновский процесс с постоянной интенсивностью λ\lambda. Покажите, что объём наибольшего (nn-мерного) шара с центром в начале координат, не содержащего ни одной точки XX, имеет показательное распределение. Выведите отсюда плотность расстояния RR от начала координат до ближайшей точки XX. Покажите, что E[R]=Γ(1/n)/{n(λc)1/n}\mathbb {E}\left[R\right] = \Gamma (1 / n) /\left\{ n(\lambda c)^{1 / n}\right\}, где cc — объём единичного шара в Rn\mathbb {R}^{n}, а Γ\Gamma — гамма-функция.

?
Задача 6.15.32

Деревня из N+1N+1 жителей охвачена эпидемией. Пусть X(t)X(t) — число заболевших в момент времени tt, и предположим, что X(0)=1X(0) = 1 и XX — процесс рождения с интенсивностями λi=λi(N+1−i)\lambda_{i} = \lambda i(N+1-i). Пусть TT — время, необходимое для того, чтобы заболели все члены популяции. Покажите, что

E[T]=1λ∑k=1N1k(N+1−k) \mathbb {E}\left[T\right] = \frac{1}{\lambda } \sum _{k = 1}^{N} \frac{1}{k(N+1-k)}

и выведите отсюда, что

E[T]=2(log⁡N+γ)λ(N+1)+O(N−2) \mathbb {E}\left[T\right] = \frac{2(\log N+\gamma )}{\lambda (N+1)}+\mathrm{O}\left(N^{-2}\right)

где γ\gamma — постоянная Эйлера. Примечательно, что E[T]\mathbb {E}\left[T\right] убывает с ростом NN при больших NN.

?
Задача 6.15.33

Частица имеет скорость V(t)V(t) в момент времени tt, причём предполагается, что V(t)V(t) принимает значения в {n+12:n≥0}\left\{ n+\frac{1}{2}: n \geq 0\right\}. Переходы в течение (t,t+h)(t, t+h) возможны следующим образом:

P(V(t+h)=w∣V(t)=v)={(v+12)h+o(h) если w=v+11−2vh+o(h) если w=v(v−12)h+o(h) если w=v−1 \mathbb {P}\left(V(t+h\right) = w \mid V(t) = v) = \begin{cases} \left(v+\frac{1}{2}\right) h+\mathrm{o}(h) & \text{ если } w = v+1 \\ 1-2 v h+\mathrm{o}(h) & \text{ если } w = v \\ \left(v-\frac{1}{2}\right) h+\mathrm{o}(h) & \text{ если } w = v-1\end{cases}

Первоначально V(0)=12V(0) = \frac{1}{2}. Пусть

G(s,t)=∑n=0∞snP(V(t)=n+12) G(s, t) = \sum _{n = 0}^{\infty } s^{n} \mathbb {P}\left(V(t) = n+\frac{1}{2}\right)
?
(a)

Покажите, что

∂G∂t=(1−s)2∂G∂s−(1−s)G \frac{\partial G}{\partial t} = (1-s)^{2} \frac{\partial G}{\partial s}-(1-s) G

и выведите отсюда, что G(s,t)={1+(1−s)t}−1G(s, t) = \left\{ 1+(1-s) t\right\}^{-1}.

(b)

Покажите, что ожидаемая длина mn(T)m_{n}(T) времени, в течение которого V=n+12V = n+\frac{1}{2} на временном интервале [0,T][0, T], задаётся формулой

mn(T)=∫0TP(V(t)=n+12)dt m_{n}(T) = \int _{0}^{T} \mathbb {P}\left(V(t) = n+\frac{1}{2}\right) d t

и что при фиксированном k,mk(T)−log⁡T→−∑i=1ki−1k, m_{k}(T)-\log T \rightarrow -\sum_{i = 1}^{k} i^{-1} при T→∞T \rightarrow \infty.

(c)

Чему равна ожидаемая скорость частицы в момент времени tt?

Задача 6.15.34

Последовательность случайных целых чисел X0,X1,…X_{0}, X_{1}, \ldots строится следующим образом. Сначала X0=0,X1=1X_{0} = 0, X_{1} = 1. При n≥1n \geq 1, при условии X0,X1,…,XnX_{0}, X_{1}, \ldots , X_{n}, следующее значение Xn+1X_{n+1} с равной вероятностью равно либо Xn+Xn−1X_{n}+X_{n-1}, либо ∣Xn−Xn−1∣\left|X_{n}-X_{n-1}\right|.

?
(a)

Является ли XX марковской цепью?

(b)

Используя марковскую цепь Yn=(Xn−1,Xn)Y_{n} = \left(X_{n-1}, X_{n}\right), найдите вероятность того, что XX достигнет значения 3 раньше, чем вновь посетит 0.

(c)

Покажите, что вероятность того, что YY когда-либо достигнет состояния (1,1)(1,1), начав из (1,2)(1,2), равна 12(3−5)\frac{1}{2}(3-\sqrt{5}).

Задача 6.15.35

Возьмём правильный шестиугольник и соединим противоположные углы прямыми линиями, пересекающимися в точке C. Частица совершает симметричное случайное блуждание по этим 7 вершинам, начиная из A(≠C)\mathrm{A}( \neq \mathrm{C}). Найдите:

?
(a)

вероятность возвращения в A без посещения C,

(b)

ожидаемое время возвращения в A,

(c)

ожидаемое число посещений C до возвращения в A,

(d)

ожидаемое время возвращения в A при условии отсутствия предшествующего посещения C.

Задача 6.15.36

Марковские цепи определяются следующими процедурами в произвольный момент времени nn :

?
(a)

Модель Бернулли. Два соседних сосуда A и B содержат каждый по mm частиц; mm частиц типа I и mm частиц типа II. В каждом сосуде наугад выбирается по частице. Если они разных типов, они меняются местами с вероятностью α\alpha, если частица типа I находится в A, либо с вероятностью β\beta, если частица типа I находится в B. Пусть XnX_{n} — число частиц типа I в A в момент времени nn.

(b)

Модель Эренфеста «собака и блохи». Два соседних сосуда содержат в сумме mm частиц. Наугад выбирается частица. Если она в A, она перемещается в B с вероятностью α\alpha, если она в B, она перемещается в A с вероятностью β\beta. Пусть YnY_{n} — число частиц в A в момент времени nn. В каждом случае найдите матрицу переходных вероятностей и стационарное распределение цепи.

Задача 6.15.37

Пусть XX — неприводимая марковская цепь с непрерывным временем на пространстве состояний SS с переходными вероятностями pjk(t)p_{j k}(t) и единственным стационарным распределением π\mathbf{\pi }, и запишем P(X(t)=j)=aj(t)\mathbb {P}\left(X(t\right) = j) = a_{j}(t). Если c(x)c(x) — вогнутая функция, покажите, что функция d(t)=∑j∈Sπjc(aj(t)/πj)d(t) = \sum_{j \in S} \pi_{j} c\left(a_{j}(t) / \pi_{j}\right) возрастает до c(1)c(1) при t→∞t \rightarrow \infty.

Относительная энтропия (или дивергенция Кульбака--Лейблера) двух строго положительных функций вероятностей f,gf, g на подмножестве SS целых чисел определяется как

D(f;g)=∑i∈Sf(i)log⁡(f(i)/g(i)) D(f ; g) = \sum _{i \in S} f(i) \log (f(i) / g(i))

Докажите, что если XX имеет конечное пространство состояний и стационарное распределение π\mathbf{\pi }, то относительная энтропия D(a(t);π)D(a(t) ; \mathbf{\pi }) монотонно убывает до 0 при t→∞t \rightarrow \infty.

?
Задача 6.15.38

В обозначениях предыдущей задачи пусть uk(t)=P(X(t)=k∣X(0)=0)u_{k}(t) = \mathbb {P}\left(X(t\right) = k \mid X(0) = 0), и предположим, что цепь обратима в равновесии (см. задачу (6.15.16)). Покажите, что u0(2t)=∑j(π0/πj)uj(t)2u_{0}(2 t) = \sum_{j}\left(\pi_{0} / \pi_{j}\right) u_{j}(t)^{2}, и выведите отсюда, что u0(t)u_{0}(t) убывает до π0\pi_{0} при t→∞t \rightarrow \infty.

?
Задача 6.15.39

Пусть Π\Pi — множество точек пуассоновского процесса на Rd\mathbb {R}^{d} с постоянной интенсивностью λ\lambda. Каждая точка смещается, причём смещения независимы и одинаково распределены. Покажите, что получившийся точечный процесс является пуассоновским процессом с интенсивностью λ\lambda.

?
Задача 6.15.40

Для удобства предположим в задаче (6.15.39), что смещения имеют непрерывную функцию распределения и конечное среднее, и что d=1d = 1. Предположим также, что первоначально вы находитесь в начале координат, а в возмущённом процессе перемещаетесь в точку aa. Пусть LR\mathrm{L}_{\mathrm{R}} — число точек, ранее находившихся слева от вас, которые теперь находятся справа, а RL\mathrm{R}_{\mathrm{L}} — число точек, ранее находившихся справа от вас, которые теперь находятся слева. Покажите, что E[LR]=E[RL]\mathbb {E}\left[\mathrm{L}_{\mathrm{R}}\right] = \mathbb {E}\left[\mathrm{R}_{\mathrm{L}}\right] тогда и только тогда, когда a=μa = \mu, где μ\mu — среднее смещение частицы.

Выведите отсюда, что если автомобили въезжают в начало длинной дороги в моменты пуассоновского процесса, имея независимые одинаково распределённые скорости, то, если вы двигаетесь со средней скоростью, в долгосрочной перспективе частота, с которой вас обгоняют другие автомобили, равна частоте, с которой вы обгоняете другие автомобили.

?
Задача 6.15.41

Муравьи заходят на кухню в моменты пуассоновского процесса NN интенсивности λ\lambda; каждый из них посещает кладовую, а затем раковину, и уходит. rr-й муравей проводит время XrX_{r} в кладовой и YrY_{r} у раковины (и Xr+YrX_{r}+Y_{r} на кухне в целом), причём векторы Vr=(Xr,Yr)V_{r} = \left(X_{r}, Y_{r}\right) и VsV_{s} независимы при r≠sr \neq s. В момент времени t=0t = 0 на кухне нет муравьёв. Найдите совместное распределение чисел A(t)A(t) муравьёв в кладовой и B(t)B(t) муравьёв у раковины в момент времени tt.

Покажите, что при t→∞t \rightarrow \infty число муравьёв на кухне сходится по распределению, при условии E[Xr+Yr]<∞\mathbb {E}\left[X_{r}+Y_{r}\right] < \infty.

Теперь предположим, что муравьи прибывают парами в моменты пуассоновского процесса, но затем разделяются и ведут себя независимо, как описано выше. Найдите совместное распределение чисел муравьёв в двух местах.

?
Задача 6.15.42

Пусть {Xr:r≥1}\left\{ X_{r}: r \geq 1\right\} — независимые показательные случайные величины с параметром λ\lambda, и положим Sn=∑r=1nXrS_{n} = \sum_{r = 1}^{n} X_{r}. Покажите, что:

?
(a)

Yk=Sk/Sn,1≤k≤n−1Y_{k} = S_{k} / S_{n}, 1 \leq k \leq n-1, имеют то же распределение, что и вариационный ряд независимых величин {Uk:1≤k≤n−1}\left\{ U_{k}: 1 \leq k \leq n-1\right\}, равномерно распределённых на (0,1)(0,1),

(b)

Zk=Xk/Sn,1≤k≤nZ_{k} = X_{k} / S_{n}, 1 \leq k \leq n, имеют то же совместное распределение, что и координаты точки (U1,…,Un)\left(U_{1}, \ldots , U_{n}\right), выбранной равномерно случайно на симплексе ∑r=1nur=1,ur≥0\sum_{r = 1}^{n} u_{r} = 1, u_{r} \geq 0 для всех rr.

Задача 6.15.43

Пусть XX — дискретная марковская цепь с конечным числом состояний и матрицей переходных вероятностей P=(pij)\mathbf{P} = \left(p_{i j}\right), где pij>0p_{i j} > 0 для всех i,ji, j. Покажите, что существует λ∈(0,1)\lambda \in (0,1), такое что ∣pij(n)−πj∣<λn\left|p_{i j}(n)-\pi_{j}\right| < \lambda^{n}, где π\mathbf{\pi } — стационарное распределение.

?
Задача 6.15.44

В условиях задачи (6.15.43) пусть Vi(n)=∑r=0n−1I{Xr=i}V_{i}(n) = \sum_{r = 0}^{n-1} I_{\left\{ X_{r} = i\right\} } — число посещений цепью состояния ii до момента nn. Покажите, что

E[∣1nVi(n)−πi∣2]→0 при n→∞ \mathbb {E}\left[\left|\frac{1}{n} V_{i}(n)-\pi _{i}\right|^{2}\right] \rightarrow 0 \quad \text{ при } n \rightarrow \infty

Покажите далее, что если ff — произвольная ограниченная функция на пространстве состояний, то

E[∣1n∑r=0n−1f(Xr)−∑i∈Sf(i)πi∣2]→0 \mathbb {E}\left[\left|\frac{1}{n} \sum _{r = 0}^{n-1} f\left(X_{r}\right)-\sum _{i \in S} f(i) \pi _{i}\right|^{2}\right] \rightarrow 0
?
Задача 6.15.45

Пусть AA и B=(B0,B1,…,Bn)\mathbf{B} = \left(B_{0}, B_{1}, \ldots , B_{n}\right) — дискретная случайная величина и вектор соответственно. Условная энтропия AA относительно B\mathbf{B} определяется как H(A∣B)=E[E[−log⁡f(A∣B]∣B])H(A \mid \mathbf{B}) = \mathbb {E}\left[\mathbb {E}\left[-\log f(A \mid \mathbf{B}\right] \mid \mathbf{B}\right]), где f(a∣b)=P(A=a∣B=b)f(a \mid \mathbf{b}) = \mathbb {P}\left(A = a \mid \mathbf{B} = \mathbf{b}\right). Пусть XX — апериодическая марковская цепь на конечном пространстве состояний. Покажите, что

H(Xn+1∣X0,X1…,Xn)=H(Xn+1∣Xn) H\left(X_{n+1} \mid X_{0}, X_{1} \ldots , X_{n}\right) = H\left(X_{n+1} \mid X_{n}\right)

и что

H(Xn+1∣Xn)→−∑iπi∑jpijlog⁡pij при n→∞ H\left(X_{n+1} \mid X_{n}\right) \rightarrow -\sum _{i} \pi _{i} \sum _{j} p_{i j} \log p_{i j} \quad \text{ при } n \rightarrow \infty

если XX апериодична с единственным стационарным распределением π\pi.

?
Задача 6.15.46

Пусть XX и YY — независимые возвратные процессы рождения и гибели с одинаковыми параметрами (и без взрывов). Не предполагается, что X0=Y0X_{0} = Y_{0}. Покажите, что:

?
(a)

для любого A⊆R,∣P(Xt∈A)−P(Yt∈A)∣→0A \subseteq \mathbb {R},\left|\mathbb {P}\left(X_{t} \in A\right)-\mathbb {P}\left(Y_{t} \in A\right)\right| \rightarrow 0 при t→∞t \rightarrow \infty,

(b)

если P(X0≤Y0)=1\mathbb {P}\left(X_{0} \leq Y_{0}\right) = 1, то E[g(Xt)]≤E[g(Yt)]\mathbb {E}\left[g\left(X_{t}\right)\right] \leq \mathbb {E}\left[g\left(Y_{t}\right)\right] для любой возрастающей функции gg.

Задача 6.15.47

Число птиц в лесу в момент времени tt — марковский процесс с непрерывным временем XX. Пищевые ресурсы накладывают ограничение 0≤X(t)≤n0 \leq X(t) \leq n. Конкуренция приводит к тому, что переходные вероятности подчиняются

pk,k+1(h)=λ(n−k)h+o(h),pk,k−1(h)=μkh+o(h) p_{k, k+1}(h) = \lambda (n-k) h+\mathrm{o}(h), \quad p_{k, k-1}(h) = \mu k h+\mathrm{o}(h)

Найдите E[sX(t)]\mathbb {E}\left[s^{X(t)}\right], а также среднее и дисперсию X(t)X(t), когда X(0)=rX(0) = r. Что происходит при t→∞t \rightarrow \infty?

?
Задача 6.15.48

Счётчик совершает неприводимое случайное блуждание по вершинам 0,1,20,1,2 треугольника на рисунке ниже, с матрицей переходных вероятностей

P=[0p0q0q10p1p2q20] \mathbf{P} =\left[\begin{smallmatrix} 0 & p_{0} & q_{0} \\ q_{1} & 0 & p_{1} \\ p_{2} & q_{2} & 0 \end{smallmatrix}\right]

где pi+qi=1p_{i}+q_{i} = 1 при всех ii. Покажите, что стационарное распределение π\pi имеет

π0=1−q2p13−q1p0−q2p1−q0p2 \pi _{0} = \frac{1-q_{2} p_{1}}{3-q_{1} p_{0}-q_{2} p_{1}-q_{0} p_{2}}

с соответствующими формулами для π1,π2\pi_{1}, \pi_{2}.

Предположим, что вы выигрываете одну песету за каждый шаг блуждания по часовой стрелке и теряете одну песету за каждый шаг против часовой стрелки. Покажите, что в равновесии средний выигрыш за шаг равен

γ=∑i(2pi−1)πi=3(2p0p1p2−p0p1−p1p2−p2p0+p0+p1+p2−1)3−q1p0−q2p1−q0p2 \gamma = \sum _{i}\left(2 p_{i}-1\right) \pi _{i} = \frac{3\left(2 p_{0} p_{1} p_{2}-p_{0} p_{1}-p_{1} p_{2}-p_{2} p_{0}+p_{0}+p_{1}+p_{2}-1\right)}{3-q_{1} p_{0}-q_{2} p_{1}-q_{0} p_{2}}

Рассмотрим теперь три случая этого процесса: A. Пусть pi=12−ap_{i} = \frac{1}{2}-a для каждого ii, где a>0a > 0. Покажите, что средний выигрыш за шаг удовлетворяет γA<0\gamma_{\mathrm{A}} < 0. B. Пусть p0=110−a,p1=p2=34−ap_{0} = \frac{1}{10}-a, p_{1} = p_{2} = \frac{3}{4}-a, где a>0a > 0. Покажите, что γB<0\gamma_{\mathrm{B}} < 0 при достаточно малых aa. C. На каждом шаге счётчик с равной вероятностью движется в соответствии с переходными вероятностями случая A или случая B, причём выбор делается независимо на каждом шаге. Покажите, что в этом случае p0=310−a,p1=p2=58−ap_{0} = \frac{3}{10}-a, p_{1} = p_{2} = \frac{5}{8}-a. Покажите, что γC>0\gamma_{\mathrm{C}} > 0 при достаточно малых aa. Тот факт, что две систематически невыгодные игры можно объединить в выгодную игру, называется парадоксом Парронда. Такие ставки в казино недоступны.

?
Задача 6.15.49

Автомобили въезжают в начало длинной дороги пуассоновским потоком интенсивности λ\lambda, начиная с момента времени t=0t = 0. Автомобиль имеет постоянную скорость V>0V > 0, являющуюся случайной величиной. Скорости автомобилей независимы, одинаково распределены и независимы от процесса въезда. Автомобили могут свободно обгонять друг друга. Покажите, что число автомобилей на первых xx милях дороги в момент времени tt имеет распределение Пуассона с параметром λE[V−1min⁡{x,Vt}]\lambda \mathbb {E}\left[V^{-1} \min \left\{ x, V t\right\} \right].

?
Задача 6.15.50

События происходят в моменты пуассоновского процесса интенсивности λ\lambda, и вам предлагается пари, основанное на этом процессе. Пусть t>0t > 0. Вам нужно произнести слово «сейчас» сразу после события, которое, как вы думаете, окажется последним, наступившим до момента tt. Вы выигрываете, если угадали, иначе проигрываете. Если до tt не произошло ни одного события, вы проигрываете. Если вы не выбрали событие до момента времени tt, вы проигрываете.

Рассмотрим стратегию, при которой вы выбираете первое событие, произошедшее после заданного момента времени ss, где 0<s<t0 < s < t.

?
(a)

Вычислите выражение для вероятности выигрыша при использовании этой стратегии.

(b)

При каком значении ss эта вероятность максимальна?

(c)

Если λt≥1\lambda t \geq 1, покажите, что вероятность выигрыша при использовании этого значения ss равна e−1e^{-1}.

Задача 6.15.51

Новый профессор Оксбриджа хочет купить дом и может позволить себе потратить до одного миллиона фунтов. Отказавшись от услуг обычных агентов по недвижимости, она обращается к своей любимой интернет-странице объявлений о недвижимости, на которой дома появляются в моменты пуассоновского процесса интенсивности λ\lambda в день. Можно считать, что цены на дома — независимые случайные величины, равномерно распределённые на интервале (800,000,2,000,000)(800,000,2,000,000). Она решает осмотреть каждый доступный по цене дом, объявленный в течение следующих 30 дней. Время, затрачиваемое на осмотр любого данного дома, равномерно распределено на промежутке (1,2)(1,2) часа. Чему равна производящая функция моментов суммарного времени, затраченного на осмотр домов?

?
Задача 6.15.52

Пусть X={Xn:n≥0}X = \left\{ X_{n}: n \geq 0\right\} — неприводимая апериодическая марковская цепь на конечном пространстве состояний SS, и пусть hij=Ei[min⁡{n≥0:Xn=j}]h_{i j} = \mathbb {E}_{i}\left[\min \left\{ n \geq 0: X_{n} = j\right\} \right] обозначает среднее время достижения jj (заметим, что hii=0h_{i i} = 0). Пусть Ki=∑jhijπjK_{i} = \sum_{j} h_{i j} \pi_{j} — среднее время достижения состояния ZZ, выбранного случайно в соответствии со стационарным распределением π\pi. Покажите, что KiK_{i} не зависит от выбора ii.

?
Задача 6.15.53

Профессор ходит пешком между домом и работой. У неё есть в общей сложности rr зонтиков, распределённых между домом и работой. Если идёт дождь, когда она выходит из дома или с работы, она берёт с собой зонт (если он доступен). Предположим, что в начале любой прогулки идёт дождь с вероятностью pp (при обычной независимости). Пусть XnX_{n} — число зонтиков, доступных ей в начале её nn-й прогулки.

?
(a)

Объясните, почему XX — марковская цепь, и запишите её матрицу переходных вероятностей.

(b)

Покажите, что цепь имеет стационарное распределение π\pi, заданное как

πi={1−pr+1−p если i=01r+1−p если i=1,2,…,r \pi _{i} = \begin{cases} \frac{1-p}{r+1-p} & \text{ если } i = 0 \\ \frac{1}{r+1-p} & \text{ если } i = 1,2, \ldots , r\end{cases}

Какая доля прогулок в долгосрочной перспективе приводит к тому, что она промокает?

(c)

Пусть r=1r = 1 и X1=1X_{1} = 1. Вычислите среднее число прогулок, совершённых до того, как она промокнет.

Задача 6.15.54

Паук взбирается по вертикальному водостоку высотой hh со скоростью 1. В моменты пуассоновского процесса постоянной интенсивности λ\lambda паук смывается обратно вниз водостока. После этого он возобновляет подъём. Пусть TT — время достижения верха, а NN — число промежуточных смываний. Покажите, что

E[e−θTsN]=(λ+θ)e−(λ+θ)hλ+θ−λs(1−e−(λ+θ)h),θ,s∈R \mathbb {E}\left[e^{-\theta T} s^{N}\right] = \frac{(\lambda +\theta ) e^{-(\lambda +\theta ) h}}{\lambda +\theta -\lambda s\left(1-e^{-(\lambda +\theta ) h}\right)}, \quad \theta , s \in \mathbb {R}

Вычисляя E[e−θT∣N=n]\mathbb {E}\left[e^{-\theta T} \mid N = n\right] или иным способом, найдите E[T∣N=n]\mathbb {E}\left[T \mid N = n\right] при n≥0n \geq 0.

?
Задача 6.15.55

Пусть XX — эргодическая марковская цепь с матрицей переходных вероятностей P\mathbf{P} и стационарным распределением π\pi. Покажите, что для любого множества состояний AA

∑i∈Aj∉Aπipij=∑i∉Aj∈Aπipij \sum _{\substack {i \in A \\ j \notin A}} \pi _{i} p_{i j} = \sum _{\substack {i \notin A \\ j \in A}} \pi _{i} p_{i j}
?
Задача 6.15.56

Турист единичной массы стоит в начале координат плоскости R2\mathbb {R}^{2}. Валуны с независимыми одинаково распределёнными массами M1,M2,…M_{1}, M_{2}, \ldots разбросаны по плоскости в точках пуассоновского процесса интенсивности 1. Пусть GRG_{R} — xx-компонента гравитационного притяжения, действующего на туриста со стороны валунов, находящихся на расстоянии не более RR от него. Гравитационную постоянную можно считать равной 1.

Покажите, что при R→∞,GRR \rightarrow \infty , G_{R} сходится по распределению к распределению Коши с характеристической функцией вида ϕ(t)=e−c∣t∣\phi (t) = e^{-c\left|t\right|}, и выразите cc через типичную массу MM.

?
Задача 6.15.57

Распределение Хольцмарка для звёздной гравитации. Пусть звёзды одинаковой массы mm расположены в точках пуассоновского процесса интенсивности 1 в R3\mathbb {R}^{3}. Пусть GRG_{R} — xx-компонента гравитационного притяжения со стороны звёзд, находящихся на расстоянии не более RR от начала координат, действующего на путешественника единичной массы в начале координат. Гравитационную постоянную можно считать равной 1.

?
(a)

Покажите, что при R→∞,GRR \rightarrow \infty , G_{R} сходится по распределению к симметричному распределению с характеристической функцией ϕ(t)=exp⁡{−c∣t∣3/2}\phi (t) = \exp \left\{ -c|t|^{3 / 2}\right\}, где c>0c > 0.

(b)

Каков будет ответ, если звёзды имеют независимые одинаково распределённые случайные массы MiM_{i}?