6.15

Задачи

[57/100%]
Показать
LaTeX
Задача 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+\vectϵ)=16\mathbb {P}\left(\mathbf{X}_{n+1} = \mathbf{X}_{n}+\vect {\epsilon }\right) = \frac{1}{6} для \vectϵ=(±1,0,0),(0,±1,0),(0,0,±1)\vect {\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) Покажите, что мера \vectv\vect {v} удовлетворяет \vectvG=0\vect {v} \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) и единственным стационарным распределением \vectπ\vect {\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 имеет конечное пространство состояний и стационарное распределение \vectπ\vect {\pi }, то относительная энтропия D(a(t);\vectπ)D(a(t) ; \vect {\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}, где \vectπ\vect {\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}?