本系列题集将在本学期长期更新,收集老师课上布置的部分课堂习题、“黑脸”(难度较大的)作业题及部分小测题。

预备知识(概率论回顾与拓展)

X1,X2,X3X_1,X_2,X_3是独立同分布的随机变量,都服从正态分布N(1,1)N(1,1),令Z=X1+X2X3,W=2X1X2+X3Z=X_1+X_2-X_3,W=2X_1-X_2+X_3.

(1)求P(Z>W)P(Z>W)(用标准正态分布的分布函数Φ\Phi表示)

点击查看答案

Z,WZ,WX1,X2,X3X_1,X_2,X_3表示,即求P(X12X2+2X3<0)P(X_1-2X_2+2X_3\lt 0)。又因为X1,X2,X3X_1,X_2,X_3均为正态分布且独立同分布,故X12X2+2X3N(1,9)X_1-2X_2+2X_3\sim N(1,9)
所以设X12X2+2X3=YX_1-2X_2+2X_3=Y,则Y13N(0,1)\frac{Y-1}{3}\sim N(0,1),故

P(Y<0)=P(Y13<13)=1Φ(13)P(Y< 0)=P\left(\frac{Y-1}{3}< -\frac{1}{3}\right)=1-\Phi\left(\frac{1}{3}\right)

(2)验证随机变量ZZWW是独立的。

点击查看答案

证明:因为

Cov(Z,W)=Cov(X1+X2X3,2X1X2+X3)=2VarX1VarX2VarX3=0\mathrm{Cov}(Z,W)=\mathrm{Cov}(X_1+X_2-X_3,2X_1-X_2+X_3)=2\mathrm{Var}X_1-\mathrm{Var}X_2-\mathrm{Var}X_3=0

所以Corr(Z,W)=0\mathrm{Corr}(Z,W)=0ZZWW不相关。又因为ZZWW均服从一维正态分布,所以由不相关可以推出独立。

甲、乙、丙三人在玩下面的游戏。每个人开始时都有且仅有11元钱。
66秒钟就会响铃,铃响时每个有钱的玩家随机选择另外两名玩家中的一名并给予该玩家11
(例如:甲和乙可能会同时决定给丙,丙可能会决定把她的钱给甲,这样第一轮游戏结束时丙22元,甲11元,乙00元。
在第二轮比赛中,乙没有钱给,但甲和丙可能会选择对方,那么在第二轮游戏结束时仍然是丙22元,甲11元,乙00元)。
那么,600600次铃响后,

(1)每个玩家都有11元钱的概率是多少?

点击查看答案

设第nn轮的状态:
AnA_n:每个玩家都有11元钱(11种可能);
BnB_n:一个玩家有11元,一个玩家有22元,剩下一个玩家没有钱(66种可能)。
满足P(An)+P(Bn)=1P(A_n)+P(B_n)=1(没有其他可能)
若上一轮状态为An1A_{n-1},则下一轮有23=82^3=8种可能,其中22种可能结果为AnA_n(甲给乙、乙给丙、丙给甲;甲给丙、丙给乙、乙给甲)
若上一轮状态为Bn1B_{n-1},则下一轮有22=42^2=4种可能(因为没钱的人无法给钱),其中只有11种可能结果为AnA_n22元的玩家给11元,11元的玩家给00元)
综上,P(An)=28P(An1)+14P(Bn1)=14P(A_n)=\frac{2}{8}P(A_{n-1})+\frac{1}{4}P(B_{n-1})=\frac{1}{4},即每轮每个玩家有11元钱的概率均为14\frac{1}{4}

(2)玩家甲恰有22元的概率?

点击查看答案

设第nn轮每个玩家都有11元钱的状态为AnA_n11种可能),玩家甲有22元(22种可能)的状态为BnB_n,有11元(非AnA_n情况)的状态为CnC_n22种可能),有00元的状态为DnD_n22种可能)。则

P(Bn)=14P(An1)+14P(Bn1)+14P(Cn1)+14P(Dn1)=14P(B_n)=\frac{1}{4}P(A_{n-1})+\frac{1}{4}P(B_{n-1})+\frac{1}{4}P(C_{n-1})+\frac{1}{4}P(D_{n-1})=\frac{1}{4}

当然也可以用对称性理解:Bn,Cn,DnB_n,C_n,D_n出现的概率是均等的(都有两种可能),且概率和为34\frac{3}{4},故每个概率分别为14\frac{1}{4}

XX为取值在[a,b][a,b]上的一个随机变量,证明:Var(X)(ba)24\mathrm{Var}(X)\leq\dfrac{(b-a)^2}{4}
(此证明思想可用于证明马尔科夫/切比雪夫不等式)

点击查看答案

证明:

\begin{aligned}
   \mathrm{Var}X&=\mathbb{E}[(X-\mathbb{E}X)^2]\\
   &=\mathbb{E}[(X-\mathbb{E}X)^2\cdot 1_{X>\mathbb{E}X}]+\mathbb{E}[(X-\mathbb{E}X)^2\cdot 1_{X\leq \mathbb{E}X}]\\
   &\leq \mathbb{E}[(b-\mathbb{E}X)^2\cdot 1_{X>\mathbb{E}X}]+\mathbb{E}[(a-\mathbb{E}X)^2\cdot 1_{X\leq \mathbb{E}X}]\\
   &=(b-\mathbb{E}X)^2\cdot\underbrace{P(X>\mathbb{E}X)}_ y+(a-\mathbb{E}X)^2\cdot\underbrace{P(X\leq \mathbb{E}X)}_{1-y}\\
   &\leq\frac{(b-a)^2}{4}
\end{aligned}

当且仅当EX=a+b2,y=12\mathbb{E}X=\frac{a+b}{2},y=\frac{1}{2}时取等号。(最后一个不等号可以用二元函数求极值的方法求解)


这里我们直接证明一个一般的不等式:若随机变量XXpp阶矩存在,则对任意x>0x>0,有

   P(|X|>x)\leq\frac{\mathbb{E}(|X|^p)}{x^p}

马尔科夫不等式与切比雪夫不等式分别对应p=1p=1p=2p=2的情形。
证明:令Y=Xp,z=xp>0Y=|X|^p,z=x^p>0,则

\mathbb{E}Y=\mathbb{E}[Y\cdot 1_{Y>z}]+\mathbb{E}[Y\cdot 1_{Y\leq z}]

又因为E[Y1Y>z]E[x1Y>z]=zP(Y>z)\mathbb{E}[Y\cdot 1_{Y>z}]\geq \mathbb{E}[x\cdot 1_{Y>z}]=zP(Y>z),所以

   P(Y>z)=P(X^p>x^p)\leq\frac{\mathbb{E}Y}{z}=\frac{|X|^p}{x^p}

又因为Xp>xp|X|^p>x^p等价于X>x|X|>x,所以不等式成立。

{Xn,n1}\{X_n,n\geq 1\}是独立同分布的随机变量序列,证明当nn\to\inftyXnn\dfrac{X_n}{n}几乎处处收敛到00当且仅当E(X1)\mathbb{E}(X_1)存在。(由此可构造依概率收敛无法推出几乎处处收敛的反例)

点击查看答案
ϵ>0\forall \epsilon>0,设Ak{Xkk>ϵ}A_k\doteq\left\{\left|\dfrac{X_k}{k}\right|>\epsilon\right\},则{Ak,k1}\{A_k,k\geq 1\}相互独立,且$
\begin{aligned}
   \sum_{k=1}^\infty P(A_k)&=\sum_{k=1}^\infty P(|X_k|>k\epsilon)\\
   &=\sum_{k=1}^\infty kP(k\epsilon< |X_k|< (k+1)\epsilon)\\
   &=\frac{1}{\epsilon}\sum_{k=1}^nk\epsilon P(k\epsilon< |X_k|< (k+1)\epsilon)\\
   &=\frac{1}{\epsilon}\sum_{k=1}^\infty\int_{k\epsilon}^{(k+1)\epsilon}k\epsilon p_{X_k}(x)dx+\frac{1}{\epsilon}\sum_{k=1}^\infty\int_{-(k+1)\epsilon}^{-k\epsilon}k\epsilon p_{X_k}(x)dx\\
   &\leq\frac{2}{\epsilon}\sum_{k=1}^\infty\int_{k\epsilon}^{(k+1)\epsilon}|x| p_{|X_k|}(x)dx\\
   &=\frac{2}{\epsilon}\int_{\epsilon}^\infty |x|p_{|X_k|}(x)dx\\
   &\leq\frac{2}{\epsilon}\mathbb{E}|X_k|< +\infty
\end{aligned}

最后一个小于是因为EXk=EX1\mathbb{E}|X_k|=\mathbb{E}|X_1|EX1\mathbb{E}X_1存在(积分绝对收敛)。于是由Borel–Cantelli引理,AnA_n可以无限次发生,limnXnna.s.0\lim_{n\to\infty}\frac{X_n}{n}\stackrel{a.s.}{\rightarrow}0
实际上,这里用了一个思想:

   \begin{aligned}
   \mathbb{E}|X|&=\int_{0}^\infty xdF(x)\\
   &=-\int_{0}^\infty xd(1-F(x))\\
   &=-x(1-F(x))|_{0}^{+\infty}+\int_{0}^\infty (1-F(x))dx\\
   &=\int_{0}^\infty P(|X|>x)dx\approx \sum_{k=1}^{\infty}\epsilon P(|X|>k\epsilon)
   \end{aligned}

无论XnX_n服从何种分布,limnP(Xnn>ϵ)=limnP(Xn>nϵ)=0\lim_{n\to\infty}P\left(\left|\dfrac{X_n}{n}\right|>\epsilon\right)=\lim_{n\to\infty}P(|X_n|>n\epsilon)=0,所以

limnXnnp0\lim_{n\to\infty}\frac{X_n}{n}\stackrel{p}{\rightarrow}0

只要EXn\mathbb{E}X_n不存在(如柯西分布),即为依概率收敛无法推出几乎处处收敛的反例。

(秘书问题)某岗位有nn个应聘者,假设应聘者的面试顺序是完全随机的。每次只面试一人,只有通过和不通过两种结果,且面试结束后立即知晓,没有通过的应聘者,面试官就失去再次招聘他的机会,一旦某位应聘者面试通过,后面的面试就自然取消。假定面试官是公正的,面试通过的人在面试官面试过的应聘者中最优。为了使得面试官能以最大的概率录用到应聘者中的最佳人选,问面试最少要面试几个人?

点击查看答案

最优策略是:
先拒绝前rr个人,然后选择第一个比前rr个人都优秀的人。对rr的计算如下:设最优秀的人顺序为第kk个(概率为1n\dfrac{1}{n}),则若krk\leq r则选到最优的概率为00(被拒绝),若k>rk>r,则要保证最优秀的人被选中需要前k1k-1个候选人中最优秀的人处于前rr个候选人中,概率为rk1\dfrac{r}{k-1}。所以,总概率为

   P(r) = \sum_{k=1}^r\frac{1}{n}\cdot 0+\sum_{k=r+1}^n \frac{1}{n} \cdot \frac{r}{k-1}

nn\to\infty时,

   \sum_{k=r+1}^n \frac{1}{n} \cdot \frac{r}{k-1}\approx\frac{r}{n}(\ln n-\ln r)=\frac{r}{n}\ln\frac{n}{r}

对关于rr的函数rnlnnr\dfrac{r}{n}\ln\dfrac{n}{r}求导,得极大值为1e\dfrac{1}{e},当且仅当r=ner=\dfrac{n}{e},即先拒绝前1e37%\dfrac{1}{e}\approx 37\%的候选人。

关于为什么这个策略是最优的,这里给出两个链接:算法趣谈(1)——秘书问题与37%法则秘书问题——知乎,涉及到动态规划的知识,在此不赘述。实则我也没大理解ww

一个算法的伪代码如下:

Initialize T in {1,2,...,10} uniformly, S=0
While T < 6
   For k in {1,2,...,T}
      S=S+k
   EndFor
   Take T in {1,2,...,10} uniformly
EndWhile
Print S

问:

(1)若第6行被注释掉,最后输出的S值平均是多少?

点击查看答案

在这一条件下,算法只对TT赋了一次值,那么While内语句就最多被执行一次。(而且只有T<6T< 6时执行)
利用重期望公式,

ES=E(E(ST))=i=110E(ST=i)P(T=i)=110i=15E(ST=i)=110i=15i(i+1)2=3.5\begin{aligned} \mathbb{E}S&=\mathbb{E}(\mathbb{E}(S|T))=\sum_{i=1}^{10}\mathbb{E}(S|T=i)P(T=i)\\ &=\frac{1}{10}\sum_{i=1}^5\mathbb{E}(S|T=i)\\ &=\frac{1}{10}\sum_{i=1}^5\frac{i(i+1)}{2}=3.5 \end{aligned}

(2)若第6行没有被注释掉,最后输出的S值平均是多少?

点击查看答案

SS分为两个部分(设第一次T<6T< 6):第一次For循环S的增量N1N_1,与第一次For循环结束到整个While结束S的增量S1S_1
那么有

ES=E(E(N1+S1T))=110i=15E(N1+S1T=i)=110i=15(E(N1T=i)+E(S1T=i))=110i=15(i(i+1)2+ES)=3.5+12ES\begin{aligned} \mathbb{E}S&=\mathbb{E}(\mathbb{E}(N_1+S_1|T))=\frac{1}{10}\sum_{i=1}^5\mathbb{E}(N_1+S_1|T=i)\\ &=\frac{1}{10}\sum_{i=1}^5(\mathbb{E}(N_1|T=i)+\mathbb{E}(S_1|T=i))\\ &=\frac{1}{10}\sum_{i=1}^5\left(\frac{i(i+1)}{2}+\mathbb{E}S\right)\\ &=3.5+\frac{1}{2}\mathbb{E}S \end{aligned}

所以ES=7\mathbb{E}S=7
此题已初见随机过程雏形

设仪器寿命服从参数为λ\lambda的指数分布,假定仪器损坏后能及时更换新仪器,而且仪器之间的寿命相互没有影响。
那么从某个仪器装上后开始算,长为tt的时间内更换的仪器数量的平均数是多少?

点击查看答案

N(t)N(t)0t0\sim t时间内仪器的更换数量(给定tt后属于随机变量),TT为仪器的寿命(服从指数分布,设其概率密度函数为f(t)f(t)),那么

E(N(t))=E(E(N(t)T))=0E(N(t)T=u)f(u)du\mathbb{E}(N(t))=\mathbb{E}(\mathbb{E}(N(t)|T))=\int_0^\infty \mathbb{E}(N(t)|T=u)f(u)du

又因为

E(N(t)T=u)={E(1+N(u,t)T=u),u<t,0,otherwise\mathbb{E}(N(t)|T=u)=\begin{cases} \mathbb{E}(1+N(u,t)|T=u),&u< t,\\ 0,&otherwise \end{cases}

其中N(u,t)N(u,t)表示utu\sim t时间内仪器的更换数量,且E(N(u,t)T=u)=E(N(tu))\mathbb{E}(N(u,t)|T=u)=\mathbb{E}(N(t-u))
所以

E(N(t))=0t(1+E(N(tu)))f(u)du\mathbb{E}(N(t))=\int_0^t (1+\mathbb{E}(N(t-u)))f(u)du

g(t)=E(N(t))g(t)=\mathbb{E}(N(t))

g(t)=0t(1+g(tu))λeλudu=0t(1+g(u))λeλ(tu)du\begin{aligned} g(t)&=\int_0^t (1+g(t-u))\lambda e^{-\lambda u}du\\ &=\int_0^t (1+g(u))\lambda e^{-\lambda (t-u)}du \end{aligned}

再令φ(t)=eλtg(t)\varphi(t)=e^{\lambda t}g(t),得

φ(t)=0tλeλu+λφ(u)du\varphi(t)=\int_0^t\lambda e^{\lambda u}+\lambda\varphi(u)du

两边求导:

φ(t)=λeλt+λφ(t)\varphi'(t)=\lambda e^{\lambda t}+\lambda\varphi(t)

利用一阶线性非齐次微分方程的知识,解得φ(t)=λteλt\varphi(t)=\lambda te^{\lambda t},即E(N(t))=g(t)=λt\mathbb{E}(N(t))=g(t)=\lambda t.

事实上,这里N(t)N(t)表示的就是随机过程里面的泊松过程。

设随机变量X,YX,Y满足EX,EY,VarX,VarY\mathbb{E}X,\mathbb{E}Y,\mathrm{Var}X,\mathrm{Var}Y均存在。
(1)试确定参数a,ba,b,使得E[(YaXb)2]\mathbb{E}[(Y-aX-b)^2]取到最小值。(由此可得到一般的线性逼近)

点击查看答案

f(a,b)=E[(YaXb)2]=EY2+a2EX2+b22aE[XY]2bEY+2abEXf(a,b)=\mathbb{E}[(Y-aX-b)^2]=\mathbb{E}Y^2+a^2\mathbb{E}X^2+b^2-2a\mathbb{E}[XY]-2b\mathbb{E}Y+2ab\mathbb{E}X

分别对aabb求偏导,得

{fa(a,b)=2aEX2+2bEX2E[XY]=0fb(a,b)=2aEX+2b2EY=0\begin{cases} f_a(a,b)=2a\mathbb{E}X^2+2b\mathbb{E}X-2\mathbb{E}[XY]=0\\ f_b(a,b)=2a\mathbb{E}X+2b-2\mathbb{E}Y=0 \end{cases}

解得

a=E[XY]EXEYEX2(EX)2=Cov(X,Y)VarX,b=EYaEXa=\frac{\mathbb{E}[XY]-\mathbb{E}X\mathbb{E}Y}{\mathbb{E}X^2-(\mathbb{E}X)^2}=\frac{\mathrm{Cov}(X,Y)}{\mathrm{Var}X},b=\mathbb{E}Y-a\mathbb{E}X

由此可得

YaX+b=Cov(X,Y)VarX(XEX)+EYY\approx aX+b=\frac{\mathrm{Cov}(X,Y)}{\mathrm{Var}X}(X-\mathbb{E}X)+\mathbb{E}Y

此为最佳线性逼近(会在线性模型中出现)

(2)证明:设g(X)g(X)为关于XX的函数,则E[(Yg(X))2]\mathbb{E}[(Y-g(X))^2]渠道最小值当且仅当g(X)=E(YX)g(X)=\mathbb{E}(Y|X)。(即最佳逼近)

点击查看答案

证明:

E[(Yg(X))2]=E[(YE(YX)+E(YX)g(X))2]=E[(YE(YX))2]+E[(E(YX)g(X))2]=E[(YE(YX))2]+E[(E(YX)g(X))2]\begin{aligned} \mathbb{E}[(Y-g(X))^2]&=\mathbb{E}[(Y-\mathbb{E}(Y|X)+\mathbb{E}(Y|X)-g(X))^2]\\ &=\mathbb{E}[(Y-\mathbb{E}(Y|X))^2]+\mathbb{E}[(\mathbb{E}(Y|X)-g(X))^2]\\ &=\mathbb{E}[(Y-\mathbb{E}(Y|X))^2]+\mathbb{E}[(\mathbb{E}(Y|X)-g(X))^2] \end{aligned}

其中

E[(YE(YX))(E(YX)g(X))]=E[E[(YE(YX))(E(YX)g(X))X]]=E[(E(YX)g(X))E[(YE(YX))X]]=E[(E(YX)g(X))(E[YX]E(YX))]=0\begin{aligned} \mathbb{E}[(Y-\mathbb{E}(Y|X))(\mathbb{E}(Y|X)-g(X))]&=\mathbb{E}[\mathbb{E}[(Y-\mathbb{E}(Y|X))(\mathbb{E}(Y|X)-g(X))|X]]\\ &=\mathbb{E}[(\mathbb{E}(Y|X)-g(X))\mathbb{E}[(Y-\mathbb{E}(Y|X))|X]]\\ &=\mathbb{E}[(\mathbb{E}(Y|X)-g(X))(\mathbb{E}[Y|X]-\mathbb{E}(Y|X))]=0 \end{aligned}

所以

E[(Yg(X))2]E[(YE(YX))2]\mathbb{E}[(Y-g(X))^2]\geq \mathbb{E}[(Y-\mathbb{E}(Y|X))^2]

去等号当且仅当g(X)=E(YX)g(X)=\mathbb{E}(Y|X)


补充:在随机向量典型题目第10题中,我们得到二元正态分布对应的条件概率
E(XY)=EX+ρVarXVarY(YEY)\mathbb{E}(X\mid Y)=\mathbb{E}X+\rho\dfrac{\sqrt{\mathrm{Var}X}}{\sqrt{\mathrm{Var}Y}}(Y-\mathbb{E}Y),这和上面第一问的线性逼近得到的结果几乎相同(除了X,YX,Y顺序不一样)
实际上,这说明二元正态分布下最优逼近就是线性逼近。

(Sanov不等式):假设X1,X2,,XmX_1,X_2,\cdots,X_m是取值在[0,1][0,1]区间上的独立同分布随机变量,记E(Xi)=p\mathbb{E}\left(X_i\right)=p.
证明:对任意q[0,1]q\in[0,1]

P(1mi=1mXiq)emD(qp)P\left(\frac{1}{m}\sum_{i=1}^{m}X_i\geq q\right)\leq e^{-mD(q||p)}

其中D(qp)=qlogqp+(1q)log1q1pD(q||p)=q\log\frac{q}{p}+(1-q)\log\frac{1-q}{1-p}.

点击查看答案

我们利用概率母函数(设为MX(s)M_X(s))证明:

P(1mi=1mXiq)=P(i=1mXimq)=P(esi=1mXiesmq)esmqE[esi=1mXi]=esmqi=1mE[esXi]=esmq[MX(s)]m=em(sqlogMX(s))\begin{aligned} P\left(\frac{1}{m}\sum_{i=1}^{m}X_i\geq q\right)&=P\left(\sum_{i=1}^{m}X_i\geq mq\right)\\ &=P\left(e^{s\sum_{i=1}^{m}X_i}\geq e^{smq}\right)\\ &\leq e^{-smq}\mathbb{E}[e^{s\sum_{i=1}^{m}X_i}]=e^{-smq}\sum_{i=1}^{m}\mathbb{E}[e^{sX_i}]\\ &=e^{-smq}[M_X(s)]^m=e^{-m(sq-\log M_X(s))} \end{aligned}

下面对m(sqlogMX(s))-m(sq-\log M_X(s))进行估计:
因为Xi[0,1]X_i\in[0,1]EXi=p\mathbb{E}X_i=p,由函数凸性,

esXi(1Xi)e0+Xies=1Xi+Xiese^{sX_i}\leq (1-X_i)e^0+X_ie^s=1-X_i+X_ie^s

两边取期望得

MX(s)1p+pesM_X(s)\leq 1-p+pe^s

代入上式:m(sqlog(1p+pes))-m(sq-\log(1-p+pe^s))ss求导,

m(qpes1p+pes)=0-m\left(q-\frac{pe^s}{1-p+pe^s}\right)=0

解得

es=q(1p)p(1q),s=logq(1p)p(1q)e^s=\frac{q(1-p)}{p(1-q)},s=\log\frac{q(1-p)}{p(1-q)}

代回原不等式,即得

m(sqlogMX(s))m(qlogq(1p)p(1q)log(1p+pq(1p)p(1q)))=mD(qp)\begin{aligned} -m(sq-\log M_X(s))&\leq -m\left(q\log\frac{q(1-p)}{p(1-q)}-\log\left(1-p+p\frac{q(1-p)}{p(1-q)}\right)\right)\\ &=-mD(q||p) \end{aligned}

证明完成。


注:这里的D(qp)D(q||p)也称为KL散度,而这个不等式说明样本均值的偏离程度会随样本增加呈指数级下降。

(Chernoff bound) 若随机变量X满足:对任意bs0b\geq s\geq 0ϕ(s)=E(esX)<+\phi\left(s\right)=\mathbb{E}\left(e^{sX}\right)< +\infty.
证明:对任意x>0x>0

P(X>x)inf0sbesxϕ(s)P\left(X>x\right)\le\inf_{0\le s\le b}{e^{-sx}\phi\left(s\right)}
点击查看答案

证明:因为

P(X>x)=P(esX>esx)esxE(esX)=esxϕ(s)\begin{aligned} P(X>x)&=P(e^{sX}>e^{sx})\\ &\leq e^{-sx}\mathbb{E}(e^{sX})\\ &=e^{-sx}\phi(s) \end{aligned}

又因为上式对任意0sb0\leq s\leq b均成立,所以P(X>x)inf0sbesxϕ(s)P(X>x)\leq\inf_{0\le s\le b}{e^{-sx}\phi\left(s\right)}。(如果ϕ(s)\phi(s)连续则可改为最小值)
(可用于估计一些非典型分布的概率)

简单随机模型

p12p\neq\dfrac{1}{2}时, 证明Bernoulli过程不是独立增量过程.

点击查看答案

独立增量过程的定义:{X(t+1)X(t),t1},X(1)\{X(t+1)-X(t),t\geq 1\},X(1)相互独立

证明:
因为P(X2X1=0X1=1)=P(X2=1)=pP(X_2-X_1=0|X_1=1)=P(X_2=1)=p,而P(X2X1=0)=P(X2=1,X1=1)+P(X2=0,X1=0)=p2+(1p)2P(X_2-X_1=0)=P(X_2=1,X_1=1)+P(X_2=0,X_1=0)=p^2+(1-p)^2
所以当p12p\neq\frac{1}{2}时,pp2+(1p)2p\neq p^2+(1-p)^2P(X2X1=0X1=1)P(X2X1=0)P(X_2-X_1=0|X_1=1)\neq P(X_2-X_1=0)P(X2X1=0)P(X_2-X_1=0)P(X1)P(X_1)不独立,得证。

已知ΘU(0,2π)\Theta\sim U\left(0,2\pi\right),对任意ttYt=Acos(wt+Θ)Y_t=A\cos\left(wt+\Theta\right)
Y={Yt;t0}Y=\{Y_t;t\geq0\}为正弦波随机过程,其中AA为振幅,ww为角频率都是常数,Θ\Theta是随机相位. 证明YY是严平稳过程。

点击查看答案

证明:
任意时刻t1,t2,,tnt_1,t_2,\cdots,t_n,与任意非负整数hh,考查联合分布(Yt1,Yt2,,Ytn)(Y_{t_1},Y_{t_2},\cdots,Y_{t_n})(Yt1+h,Yt2+h,,Ytn+h)(Y_{t_1+h},Y_{t_2+h},\cdots,Y_{t_n+h})
因为

Yt1+h=Acos(w(t1+h)+Θ)=Acos(wt1+Θ+wh)Y_{t_1+h}=A\cos(w(t_1+h)+\Theta)=A\cos(wt_1+\Theta+wh)

(Θ+wh)mod2πU(0,2π)(\Theta+wh)\mod 2\pi\sim U(0,2\pi),所以(Θ+wh)mod2π(\Theta+wh)\mod 2\piΘ\Theta同分布,
那么
Acos(wt1+Θ+wh)A\cos(wt_1+\Theta+wh)Acos(wt1+Θ)A\cos(wt_1+\Theta)同分布,进而得到
(Yt1,Yt2,,Ytn)(Y_{t_1},Y_{t_2},\cdots,Y_{t_n})(Yt1+h,Yt2+h,,Ytn+h)(Y_{t_1+h},Y_{t_2+h},\cdots,Y_{t_n+h})同分布,YY的严平稳性得证。

定义随机过程Z={Zt;t0}Z=\{Z_t;t\geq0\},其中Zt=Xcos(2πt)+Ysin(2πt)Z_t=X\cos\left(2\pi t\right)+Y\sin\left(2\pi t\right), 其中XXYY是两个不相关的随机变量,分布不同但E[X]=E[Y]=0,D[X]=D[Y]=1\mathbb{E}[X]=\mathbb{E}[Y]=0,D[X]=D[Y]=1.
(1)证明ZZ是平稳过程.
(2)若进一步假设XXYY独立且都服从正态分布, 证明ZZ是严平稳过程.
(3)求ZZ的二维分布概率密度函数族.

点击查看答案

(1)证明:
因为

E[Zt]=E[X]cos(2πt)+E[Y]sin(2πt)=0\mathbb{E}[Z_t]=\mathbb{E}[X]\cos(2\pi t)+\mathbb{E}[Y]\sin(2\pi t)=0

tt无关;又因为E[X2]=D[X]+(E[X])2=1=E[Y2],E[XY]=Cov(X,Y)=0\mathbb{E}[X^2]=D[X]+(\mathbb{E}[X])^2=1=\mathbb{E}[Y^2],\mathbb{E}[XY]=\mathrm{Cov}(X,Y)=0X,YX,Y不相关),所以

Cov(Zs,Zt)=E[ZsZt]=E[(Xcos(2πs)+Ysin(2πs))(Xcos(2πt)+Ysin(2πt))]=E[X2]=1cos(2πs)cos(2πt)+E[Y2]=1sin(2πs)sin(2πt)+E[XY]=0()=cos(2πs)cos(2πt)+sin(2πs)sin(2πt)=cos(2π(st))\begin{aligned} \mathrm{Cov}(Z_s,Z_t)&=\mathbb{E}[Z_sZ_t]\\ &=\mathbb{E}[(X\cos(2\pi s)+Y\sin(2\pi s))(X\cos(2\pi t)+Y\sin(2\pi t))]\\ &=\underbrace{\mathbb{E}[X^2]}_{=1}\cos(2\pi s)\cos(2\pi t)+\underbrace{\mathbb{E}[Y^2]}_{=1}\sin(2\pi s)\sin(2\pi t)+\underbrace{\mathbb{E}[XY]}_{=0}(\dots)\\ &=\cos(2\pi s)\cos(2\pi t)+\sin(2\pi s)\sin(2\pi t)=\cos(2\pi(s-t)) \end{aligned}

s,ts,t无关,所以ZZ为宽平稳过程。
(2)证明:
因为X,YX,Y相互独立且均服从正态分布N(0,1)N(0,1),所以(X,Y)(X,Y)服从二维标准正态分布

N((00),(1001))\mathcal{N} \left( \begin{pmatrix} 0 \\ 0 \end{pmatrix}, \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} \right)

Zt+h=Xcos(2π(t+h))+Ycos(2π(t+h))=X(cos2πtcos2πhsin2πtsin2πh)+Y(cos2πtcos2πhsin2πtsin2πh)=(Xcos2πh+Ycos2πh)cos2πt(Xsin2πh+Ysin2πh)sin2πt\begin{aligned} Z_{t+h}&=X\cos(2\pi(t+h))+Y\cos(2\pi(t+h))\\ &=X(\cos2\pi t\cdot\cos2\pi h-\sin2\pi t\cdot\sin2\pi h)+Y(\cos2\pi t\cdot\cos2\pi h-\sin2\pi t\cdot\sin2\pi h)\\ &=(X\cos2\pi h+Y\cos2\pi h)\cos2\pi t-(-X\sin2\pi h+Y\sin2\pi h)\sin2\pi t \end{aligned}

X=Xcos2πh+Ycos2πh,Y=Xsin2πh+Ysin2πhX'=X\cos2\pi h+Y\cos2\pi h,Y'=-X\sin2\pi h+Y\sin2\pi h,则(X,Y)(X,Y)(X,Y)(X',Y')是单位正交变换,所以(X,Y)(X',Y')同样满足二维标准正态分布。
所以Zt+hZ_{t+h}ZtZ_t同分布,进而任意t1<t2<tnt_1< t_2\cdots< t_nhh
(Zt1,Zt2,,Ztn)(Z_{t_1},Z_{t_2},\cdots,Z_{t_n})(Zt1+h,Zt2+h,,Ztn+h)(Z_{t_{1+h}},Z_{t_{2+h}},\cdots,Z_{t_{n+h}})同分布,ZZ为严平稳过程。
(3)即求(Zs,Zt)(0s<t)(Z_s,Z_t)(\forall 0\leq s< t)的概率密度函数。由前两问可知,(Zs,Zt)(Z_s,Z_t)服从二维正态分布N(0,0,1,1,cos(2π(st)))N(0,0,1,1,\cos(2\pi(s-t)))(具体概率密度函数略)

S={Sn;n1}S=\{S_n;n\geq 1\}(q,p)(q, p)-简单随机游动, 其中q=1p,p>12q=1-p,p > \frac{1}{2}τ1=inf{n>0:Sn=1}\tau_1=\inf\{n > 0: S_n = 1\}
(1)求τ1\tau_1的分布(默认初始位置为00)。

点击查看答案

因为SS是简单随机游动且右连续,所以

P(τ1=nS0=0)=1nP(Sn=1S0=0)={1n(nn+12)pn+12qn12,n+12Z0,otherwise.\begin{aligned} P(\tau_1=n|S_0=0)&=\frac{1}{n}P(S_n=1|S_0=0)\\ &=\begin{cases} \frac{1}{n}\binom{n}{\frac{n+1}{2}}p^{\frac{n+1}{2}}q^{\frac{n-1}{2}},&\frac{n+1}{2}\in\Z\\ 0,&otherwise. \end{cases} \end{aligned}

(2)求E(τ1S0=0)\mathbb{E}(\tau_1|S_0 = 0).

点击查看答案
E(τ1S0=0)=E[E(τ1S1,S0=0)S0=0]=E(τ1S1=1,S0=0)P(S1=1S0=0)+E(τ1S1=1,S0=0)P(S1=1S0=0)=1P(S1=1S0=0)+E(1+τ1+τ1S1=1,S0=0)P(S1=1S0=0)\begin{aligned} \mathbb{E}(\tau_1|S_0=0)&=\mathbb{E}[\mathbb{E}(\tau_1|S_1,S_0=0)|S_0=0]\\ &=\mathbb{E}(\tau_1|S_1=1,S_0=0)P(S_1=1|S_0=0)+\mathbb{E}(\tau_1|S_1=-1,S_0=0)P(S_1=-1|S_0=0)\\ &=1\cdot P(S_1=1|S_0=0)+\mathbb{E}(1+\tau'_1+\tau''_1|S_1=-1,S_0=0)P(S_1=-1|S_0=0) \end{aligned}

其中τ1\tau'_1τ1\tau''_1分别为从1-1出发首次到00的时间与首次到00后,从00出发到11的时间.
然后

E(1+τ1+τ1S1=1,S0=0)P(S1=1S0=0)=P(S1=1S0=0)+E(τ1)P(S1=1S0=0)+E(τ1)P(S1=1S0=0)\begin{aligned} &\mathbb{E}(1+\tau'_1+\tau''_1|S_1=-1,S_0=0)P(S_1=-1|S_0=0)\\ =&P(S_1=-1|S_0=0)+\mathbb{E}(\tau'_1)P(S_1=-1|S_0=0)+\mathbb{E}(\tau''_1)P(S_1=-1|S_0=0) \end{aligned}

E(τ1)=E(τ1)=E(τ1S0=0)\mathbb{E}(\tau'_1)=\mathbb{E}(\tau''_1)=\mathbb{E}(\tau_1|S_0=0),所以

E(τ1S0=0)=p+q+2qE(τ1S0=0)=1+2qE(τ1S0=0)\mathbb{E}(\tau_1|S_0=0)=p+q+2q\mathbb{E}(\tau_1|S_0=0)=1+2q\mathbb{E}(\tau_1|S_0=0)

E(τ1S0=0)=112q\mathbb{E}(\tau_1|S_0=0)=\dfrac{1}{1-2q}.(这也是为什么要求q<12q< \dfrac{1}{2}

(3)证明P(τ1<)=1P(\tau_1< \infty)=1.(即从任意位置出发,一定能经过有限次游动回到位置11)。

点击查看答案

证明思路1:设f(i)=Pi(τ1<)f(i)=P_i(\tau_1< \infty)表示从位置ii出发,最终经过有限次到达11的概率。则

f(i)=pf(i+1)+qf(i1)f(i)=pf(i+1)+qf(i-1)

解得

f(i)=A+B(qp)if(i)=A+B\left(\frac{q}{p}\right)^i

其中A,BA,B为常数。又因为

{f(0)=1f(i)1,iZqp<1\begin{cases} f(0)=1\\ f(i)\leq 1,\forall i\in\Z\\ \dfrac{q}{p}< 1 \end{cases}

所以

B=0,A=1f(i)1B=0,A=1\Longrightarrow f(i)\equiv 1

所以P(τ1<)=1P(\tau_1< \infty)=1.

注:上面用到的差分方程在处理简单随机游动问题时非常有用,务必牢记。

证明思路2(使用条件期望):

P(τ1<)=E[E(1{τ1<}W1)]=E(1{τ1<}W1=1)P(W1=1)+E(1{τ1<}W1=1)P(W1=1)=p+qE(1{τ1+τ1<}W1=1)=p+qE(1{τ1<})E(1{τ1<})=p+q[P(τ1<)]2\begin{aligned} P(\tau_1< \infty) &= \mathbb{E}[\mathbb{E}(\mathbf{1}_{\{\tau_1 <\infty\}}|W_1)]\\ &=\mathbb{E}(\mathbf{1}_{\{\tau_1 <\infty\}}|W_1=1)P(W_1=1)+\mathbb{E}(\mathbf{1}_{\{\tau_1 <\infty\}}|W_1=-1)P(W_1=-1)\\ &=p+q\mathbb{E}(\mathrm{1}_{\{\tau'_1+\tau''_1<\infty\}}|W_1=-1)\\ &=p+q\mathbb{E}(\mathrm{1}_{\{\tau'_1<\infty\}})\cdot \mathbb{E}(\mathrm{1}_{\{\tau''_1<\infty\}})\\ &=p+q[P(\tau_1<\infty)]^2 \end{aligned}

因此P(τ1<)=1±14pq2q=1±pq2qP(\tau_1< \infty)=\dfrac{1\pm\sqrt{1-4pq}}{2q}=\dfrac{1\pm|p-q|}{2q}
又因为p>qp>q,且P(τ1<)1P(\tau_1<\infty)\leq 1,所以只能

P(τ1<)=1(pq)2q=1P(\tau_1<\infty)=\frac{1-(p-q)}{2q}=1

pqp\leq q时,P(τ1<)=pqP(\tau_1<\infty)=\dfrac{p}{q}.

S={Sn;n0}S=\{S_n;n\geq 0\}是初值为00(q,p)(q, p)-简单随机游动,用TT表示SS在首次到达11之前回到00的次数。
已知qpq\leq p, 对任意k0k\geq 0, 证明P(T=k)=qkpP(T = k) = q^kp.

点击查看答案

证明:
设事件AiA_i表示第ii次从00出发,并在到达11前回到00
那么P(Ai)=p1(0,1)=qP(A_i)=p_1(0,-1)=q(因为第一步到达1-1后,P1(τ0<)=1P_{-1}(\tau_0< \infty)=1,此为上题第3问结论),所以

P(T=k)=P(A1A2AkAk+1)=qk(1q)=qkp.P(T=k)=P(A_1A_2\cdots A_k\overline{A_{k+1}})=q^k(1-q)=q^kp.

对初值为00的简单对称随机游动WW,求P(W10,W20,,W2n10,W2n=0)=P(W_1\neq 0,W_2\neq 0,\cdots,W_{2n-1}\neq 0,W_{2n}=0)=(2nn)(2n1)22n\frac{\binom{2n}{n}}{(2n-1)2^{2n}}

这一过程等价于从00出发,经过2n2n步后第一次回到00
将情况分为22种(另一种对称):第一步移动到11,之后位置一直大于00,到2n12n-1步回到11。下面讨论这种情况的路径数。
首先,所有从11开始经过2n22n-2步后回到11的路径数为(2n2n1)\displaystyle\binom{2n-2}{n-1}(即向上n1n-1步+向下n1n-1步)。
然后中间回到过00的路径:根据反射原理,

P(W2n1=1,min2k2n1Wk1S1=1)=P(W2n2=1,min1k2n2Wk0S0=1)=P(W2n2=1S0=1)\begin{aligned} &P(W_{2n-1}=1,\min_{2\leq k\leq 2n-1}W_k\leq 1|S_1=1)\\ =&P(W_{2n-2}=1,\min_{1\leq k\leq 2n-2}W_k\leq 0|S_0=1)\\ =&P(W_{2n-2}=1|S_0=-1) \end{aligned}

因此路径数为(2n2n)\displaystyle\binom{2n-2}{n}nn条上+n2n-2条下)。
所以最终的路径数量为

2((2n2n1)(2n2n))=2n(2n2n1)=12n1(2nn)2\left(\binom{2n-2}{n-1}-\binom{2n-2}{n}\right)=\frac{2}{n}\binom{2n-2}{n-1}=\frac{1}{2n-1}\binom{2n}{n}

而总路径数为22n2^{2n},故概率为(2nn)(2n1)22n\dfrac{\binom{2n}{n}}{(2n-1)2^{2n}}

S={Sn;n0}S=\{S_n;n\geq 0\}是右连续随机游动,步长XX的分布为P(X=k)=pk,k=1,0,1,2,P(X=-k)=p_k,k = -1,0,1,2,\cdots。记m=E(X)m=\mathbb{E}(|X|),令τ1=inf{n>0;Sn=1}\tau_1 = \inf\{n > 0; S_n= 1\},求E(τ1S0=0)\mathbb{E}(\tau_1|S_0= 0).

点击查看答案

E(τ1S0=0)=a\mathbb{E}(\tau_1 | S_0 = 0)=a,则由条件期望的平滑性及全期望公式:

E(τ1S0=0)=E[E(τ1S1,S0)S0=0]=1P(S1=1S0=0)+E(τ1S1=0,S0=0)P(S1=0S0=0)+k=1E(τ1S1=k,S0=0)P(S1=kS0=0)\begin{aligned} \mathbb{E}(\tau_1 | S_0 = 0) &= \mathbb{E}[\mathbb{E}(\tau_1 | S_1, S_0) | S_0 = 0]\\ &= 1\cdot P(S_1=-1|S_0=0)+\mathbb{E}(\tau_1 | S_1=0, S_0=0)P(S_1=0|S_0=0)+\sum_{k=1}^\infty \mathbb{E}(\tau_1 | S_1=k, S_0=0)P(S_1=k|S_0=0) \end{aligned}

其中E(τ1S1=0,S0=0)=1+E(τ1S0=0)=1+a\mathbb{E}(\tau_1 | S_1=0, S_0=0)=1+\mathbb{E}(\tau_1 | S_0 = 0) = 1+aE(τ1S1=k,S0=0)=1+i=0kE(τ1S0=0)=(k+1)a\mathbb{E}(\tau_1 | S_1=k, S_0=0)=1+\sum_{i=0}^k\mathbb{E}(\tau_1 | S_0=0)=(k+1)a(即从k-k11经过k+1k+1个首回时),因此

a=p1+p0(1+a)+k=1pk[1+(k+1)a]=1+ap0+ak=1pk(k+1)=1+a(1p1)+ak=1kpk\begin{aligned} a&=p_{-1}+p_0(1+a)+\sum_{k=1}^\infty p_k\cdot[1+(k+1)a]\\ &=1+ap_0+a\sum_{k=1}^\infty p_k\cdot(k+1)\\ &=1+a(1-p_{-1})+a\sum_{k=1}^\infty kp_k \end{aligned}

因此a=1p1k=1kpka=\dfrac{1}{p_{-1}-\sum_{k=1}^\infty kp_k}。又因为

m=E(X)=p1+k=1kpkm=\mathbb{E}(|X|)=p_{-1}+\sum_{k=1}^\infty kp_k

所以

a=12p1m.a=\frac{1}{2p_{-1}-m}.

N1,N2N_1,N_2分别是强度为λ1,λ2\lambda_1, \lambda_2的泊松过程且相互独立. 证明N=N1+N2N = N_1 + N_2是强度为λ1+λ2\lambda_1 + \lambda_2的泊松过程.

点击查看答案

N1(t)N_1(t)N2(t)N_2(t)是两个独立的泊松过程,强度分别为 λ1\lambda_1λ2\lambda_2
根据泊松过程的定义,任意t>0t > 0

N1(t)Poisson(λ1t),N2(t)Poisson(λ2t)N_1(t)\sim \text{Poisson}(\lambda_1 t),N_2(t) \sim \text{Poisson}(\lambda_2 t)

且相互独立。因此N(t)=N1(t)+N2(t)Poisson((λ1+λ2)t)N(t)=N_1(t)+N_2(t)\sim\text{Poisson}((\lambda_1+\lambda_2) t)
又因为N1N_1N2N_2均为独立平稳增量过程,所以NN也是独立平稳增量过程。因此NN是强度为λ1+λ2\lambda_1 + \lambda_2的泊松过程。

TkT_k表示NN的第kk次随机事件发生的时间.对任意t>x>0t>x>0,求:
(1)P(tTk>xN(t)=k)P(t-T_k>x|N(t)=k)

点击查看答案

k=0k=0时,P(tT0>xN(t)=0)=P(t>xN(t)=0)=1P(t-T_0>x|N(t)=0)=P(t>x|N(t)=0)=1
下面讨论k1k\geq 1情况:
因为给定N(t)=kN(t)=k时,TkT_k与均匀分布U(0,t)U(0,t)的最大次序统计量同分布,故TkT_k的概率密度函数为k(xt)k11t\displaystyle k\left(\frac{x}{t}\right)^{k-1}\frac{1}{t},因此

P(tTk>xN(t)=k)=P(Tk<txN(t)=k)=0txk(ut)k11tdu=ktk0txuk1du=(txt)k\begin{aligned} P(t-T_k>x|N(t)=k)&=P(T_k< t-x|N(t)=k)\\ &=\int_{0}^{t-x}k\left(\frac{u}{t}\right)^{k-1}\frac{1}{t}du\\ &=\frac{k}{t^k}\int_{0}^{t-x}u^{k-1}du=\left(\frac{t-x}{t}\right)^k \end{aligned}

综上,P(tT0>xN(t)=k)=(1xt)k\displaystyle P(t-T_0>x|N(t)=k)=\left(1-\frac{x}{t}\right)^k.

(2)P(tTN(t)>x)P(t-T_{N(t)}>x)P(TN(t)+1t>x)P(T_{N(t)+1}-t>x)

点击查看答案
P(tTN(t)>x)=E[P(tTN(t)>x)N(t)]=k=0P(TN(t)<txN(t)=k)P(N(t)=k)=k=0(1xt)k(λt)kk!eλt=eλti=0(λ(tx))kk!=eλteλ(tx)=eλx\begin{aligned} P(t-T_{N(t)}>x)&=\mathbb{E}[P(t-T_{N(t)}>x)|N(t)]\\ &=\sum_{k=0}^\infty P(T_{N(t)}< t-x|N(t)=k)P(N(t)=k)\\ &=\sum_{k=0}^\infty \left(1-\frac{x}{t}\right)^k\cdot\frac{(\lambda t)^k}{k!}e^{-\lambda t}\\ &=e^{-\lambda t}\sum_{i=0}^\infty\frac{(\lambda (t-x))^k}{k!}\\ &=e^{-\lambda t}\cdot e^{\lambda(t-x)}=e^{-\lambda x} \end{aligned}P(TN(t)+1t>x)=P(TN(t)+1>t+x)=k=0P(N(t+x)N(t)=0N(t)=k)P(N(t)=k)=P(N(t+x)N(t)=0)=P(N(x)=0)=eλx\begin{aligned} P(T_{N(t)+1}-t>x)&=P(T_{N(t)+1}>t+x)\\ &=\sum_{k=0}^\infty P(N(t+x)-N(t)=0|N(t)=k)P(N(t)=k)\\ &=P(N(t+x)-N(t)=0)\\ &=P(N(x)=0)=e^{-\lambda x} \end{aligned}

可以发现这两个概率相等,这是由指数分布的无记忆性以及泊松过程的平稳增量性决定的。

(3)以WkW_k表示NN的第kk次事件的间隔时间,对任意x>0x > 0,求P(WN(t)+1>x)P(W_{N(t)+1}>x)

点击查看答案
P(WN(t)+1>x)=P(TN(t)+1TN(t)>x)=k=0P(TN(t)+1TN(t)>x,N(t)=k)\begin{aligned} P(W_{N(t)+1}>x)&=P(T_{N(t)+1}-T_{N(t)}>x)\\ &=\sum_{k=0}^\infty P(T_{N(t)+1}-T_{N(t)}>x,N(t)=k) \end{aligned}

k=0k=0时,

P(T1>x,N(t)=0)=P(T1>x,T1>t)=eλmax(x,t)\begin{aligned} P(T_1>x,N(t)=0)&=P(T_1>x,T_1>t)\\ &=e^{-\lambda\max(x,t)} \end{aligned}

k>0k>0时,

P(TN(t)+1TN(t)>x,N(t)=k)=P(Tkt,Tk+1>t,Tk+1Tk>x)=P(Tkt,Wk+1>max(x,tTk))\begin{aligned} &P(T_{N(t)+1}-T_{N(t)}>x,N(t)=k)\\ &=P(T_k\leq t,T_{k+1}>t,T_{k+1}-T_k>x)\\ &=P(T_k\leq t,W_{k+1}>\max(x,t-T_k))\\ \end{aligned}

又因为TkT_kWk+1W_{k+1}相互独立,TkGa(k,λ),Wk+1Exp(λ)T_k\sim Ga(k,\lambda),W_{k+1}\sim\mathrm{Exp}(\lambda),所以

P(Tkt,Wk+1>max(x,tTk))=0tλkuk1eλu(k1)!dumax(x,tu)λeλvdv=0tλkuk1eλu(k1)!eλmax(x,tu)du\begin{aligned} &P(T_k\leq t,W_{k+1}>\max(x,t-T_k))\\ &=\int_{0}^t\frac{\lambda^ku^{k-1}e^{-\lambda u}}{(k-1)!}du\int_{\max(x,t-u)}^\infty\lambda e^{-\lambda v}dv\\ &=\int_{0}^t\frac{\lambda^ku^{k-1}e^{-\lambda u}}{(k-1)!}\cdot e^{-\lambda\max(x,t-u)}du \end{aligned}

因此

k=1P(Tkt,Wk+1>max(x,tTk))=k=10tλkuk1eλu(k1)!eλmax(x,tu)du=0tλk=1(λu)k1eλu(k1)!eλmax(x,tu)du=0tλeλmax(x,tu)du=λ0teλmax(x,u)du={λteλx,xt,(1+λx)eλxeλt,x<t\begin{aligned} &\sum_{k=1}^\infty P(T_k\leq t,W_{k+1}>\max(x,t-T_k))\\ &=\sum_{k=1}^\infty\int_{0}^t\frac{\lambda^ku^{k-1}e^{-\lambda u}}{(k-1)!}\cdot e^{-\lambda\max(x,t-u)}du\\ &=\int_{0}^t\lambda\sum_{k=1}^\infty\frac{(\lambda u)^{k-1}e^{-\lambda u}}{(k-1)!}\cdot e^{-\lambda\max(x,t-u)}du\\ &=\int_0^t\lambda e^{-\lambda\max(x,t-u)}du\\ &=\lambda\int_0^te^{-\lambda\max(x,u)}du\\ &=\begin{cases} \lambda t e^{-\lambda x}, & x \geq t, \\ (1 + \lambda x)e^{-\lambda x} - e^{-\lambda t}, & x < t \end{cases} \end{aligned}

综上,

P(WN(t)+1>x)=k=0P(WN(t)+1>x,N(t)=k)={eλx+λteλx,xt,eλt+(1+λx)eλxeλt,x<t=(1+λmax(x,t))eλx.\begin{aligned} P(W_{N(t)+1}>x)&=\sum_{k=0}^\infty P(W_{N(t)+1}>x,N(t)=k)\\ &=\begin{cases} e^{-\lambda x}+\lambda t e^{-\lambda x}, & x \geq t, \\ e^{-\lambda t}+(1 + \lambda x)e^{-\lambda x} - e^{-\lambda t}, & x < t \end{cases}\\ &=(1+\lambda\max(x,t))e^{-\lambda x}. \end{aligned}

离散时间随机游动假定了每一次游动都在整数时刻点发生,如果游动可以在任意时刻点发生,我们可以怎样建模?试着用泊松过程和简单(q,p)(q,p)-随机游动构造一个随机过程模型,并在p>qp>q时计算从00位置出发首次到达位置11的平均时间。

提示:连续时间随机游动

点击查看答案

模型建立:设(q,p)(q,p)-随机游动{Xt,t0}\{X_t,t\geq 0\}的游动时刻服从强度为λ\lambda的泊松过程,即在[0,t][0,t]时间内发生N(t)N(t)次游动,相邻游动时间间隔服从指数分布Exp(λ)\mathrm{Exp}(\lambda),则可以定义连续时间随机游动:

S(t)=XN(t)S(t)=X_{N(t)}

下面求τ=inf{t0;S(t)=1}\tau=\inf\{t\geq 0;S(t)=1\}
因为离散时间下E(τ1S0=0)=112q\mathbb{E}(\tau_1|S_0=0)=\dfrac{1}{1-2q}(因为p>qp>q),即首次到达位置11游动步数的期望为112q\dfrac{1}{1-2q};又因为游动时间间隔的期望为1λ\dfrac{1}{\lambda},所以

E(τ)=112q1λ=1λ(12q)\mathbb{E}(\tau)=\frac{1}{1-2q}\cdot\frac{1}{\lambda}=\frac{1}{\lambda(1-2q)}

某个零件在运行过程中收到撞击,用N(t)N(t)表示(0,t](0,t]时间段内该零件受到的撞击次数。经验表明N=(N(t):t0)N=(N(t):t\geq 0)为参数λ\lambda的泊松过程。
每次撞击会给零件带来一定的磨损,假定每次磨损量为ξ1,ξ2,\xi_1,\xi_2,\cdots均服从参数为β\beta的指数分布且相互独立。若磨损总量达到或超过α>0\alpha>0就要更换零件,问该零件的平均寿命。

点击查看答案

设零件的寿命为XX,则

X=inf{t;i=1N(t)ξiα}X=\inf\left\{t;\sum_{i=1}^{N(t)}\xi_i\geq \alpha\right\}

EX=0tdF(t)=0P(X>t)dt=0P(i=1N(t)ξi<α)dt=0k=0P(i=1N(t)ξi<αN(t)=k)P(N(t)=k)dt=0k=0P(i=1kξi<α)P(N(t)=k)dt\begin{aligned} \mathbb{E}X&=\int_{0}^\infty tdF(t)=\int_{0}^\infty P(X>t)\mathrm{d}t\\ &=\int_{0}^\infty P\left(\sum_{i=1}^{N(t)}\xi_i< \alpha\right)\mathrm{d}t\\ &=\int_{0}^\infty\sum_{k=0}^\infty P\left(\sum_{i=1}^{N(t)}\xi_i< \alpha|N(t)=k\right)P(N(t)=k)\mathrm{d}t\\ &=\int_{0}^\infty\sum_{k=0}^\infty P\left(\sum_{i=1}^k\xi_i< \alpha\right)P(N(t)=k)\mathrm{d}t \end{aligned}

{ξi,i1}\{\xi_i,i\geq 1\}均服从Exp(β)\mathrm{Exp}(\beta),所以可以将Tk={i=1kξi,k1}\displaystyle T_k'=\left\{\sum_{i=1}^k\xi_i,k\geq 1\right\}看作强度为β\beta的泊松过程{N(t),t0}\{N'(t),t\geq 0\}的到达时间。于是

P(i=1kξi<α)=P(N(α)k)=1n=0keβα(βα)nn!=eβαn=k+1(βα)nn!P\left(\sum_{i=1}^k\xi_i< \alpha\right)=P(N'(\alpha)\geq k)=1 - \sum_{n=0}^k e^{-\beta\alpha} \frac{(\beta\alpha)^n}{n!} = e^{-\beta \alpha} \sum_{n=k+1}^\infty \frac{(\beta \alpha)^n}{n!}

因此

EX=0k=0[eβαn=k+1(βα)nn!](λt)kk!eλtdt=0eβαn=1(βα)nn![k=0n1(λt)kk!]eλtdt=0eβαn=1(βα)nn!P(N(t)n1)dt=eβαn=1(βα)nn!0P(N(t)n1)dt=eβαn=1(βα)nn!nλ=1λeβαn=1n(βα)nn!=1λeβα(βα)eβα=βαλ.\begin{aligned} \mathbb{E}X&=\int_{0}^\infty\sum_{k=0}^\infty\left[e^{-\beta \alpha} \sum_{n=k+1}^\infty\frac{(\beta \alpha)^n}{n!}\right]\cdot\frac{(\lambda t)^k}{k!}e^{-\lambda t}\mathrm{d}t\\ &=\int_{0}^{\infty} e^{-\beta\alpha} \sum_{n=1}^{\infty} \frac{(\beta\alpha)^n}{n!} \left[ \sum_{k=0}^{n-1} \frac{(\lambda t)^k}{k!} \right] e^{-\lambda t}\mathrm{d}t\\ &=\int_{0}^{\infty} e^{-\beta\alpha} \sum_{n=1}^{\infty} \frac{(\beta\alpha)^n}{n!} P(N(t) \leq n-1)\mathrm{d}t\\ &=e^{-\beta\alpha} \sum_{n=1}^{\infty} \frac{(\beta\alpha)^n}{n!} \int_{0}^{\infty} P(N(t) \leq n-1)\mathrm{d}t\\ &= e^{-\beta \alpha} \sum_{n=1}^\infty \frac{(\beta \alpha)^n}{n!} \cdot \frac{n}{\lambda}\\ &=\frac{1}{\lambda} e^{-\beta \alpha} \sum_{n=1}^\infty \frac{n(\beta \alpha)^n}{n!}\\ &=\frac{1}{\lambda} e^{-\beta \alpha} \cdot (\beta \alpha) e^{\beta \alpha}\\ &=\frac{\beta\alpha}{\lambda}. \end{aligned}

离散时间马氏链

(存储模型)设ξi,i=1,2,\xi_i, i = 1, 2, \dots为独立同分布的非负整数值随机变量,令

X0=K,Xn={0(Xn1ξn),k<Xn1K,0(Kξn),Xn1k,n1,X_0 = K, \quad X_n = \begin{cases} 0 \lor (X_{n-1} - \xi_n), & k < X_{n-1} \leq K, \\ 0 \lor (K - \xi_n), & X_{n-1} \leq k, \end{cases} \quad n \geq 1,

其中k<Kk < K为给定的两个正整数。
(1)证明X={Xn,n0}X = \{X_n, n \geq 0\}为离散时间马氏链。

点击查看答案

(1)证明:显然XX的状态空间S={0,1,,K}S = \{0, 1, \cdots, K\}。对任意n1n \geq 1以及i1,,in+1Si_1, \cdots, i_{n+1} \in S

P(Xn+1=in+1Xn=in,,X1=i1,X0=K)={P((Kξn+1)0=in+1Xn=in,,X0=K),ink,P((inξn+1)0=in+1Xn=in,,X0=K),k<inK,\begin{aligned} &P(X_{n+1} = i_{n+1} | X_n = i_n, \cdots, X_1 = i_1, X_0 = K) \\ &=\begin{cases} P((K - \xi_{n+1}) \lor 0 = i_{n+1} | X_n = i_n, \cdots, X_0 = K), & i_n \leq k, \\ P((i_n - \xi_{n+1}) \lor 0 = i_{n+1} | X_n = i_n, \cdots, X_0 = K), & k < i_n \leq K, \end{cases} \end{aligned}

对任意i1i \geq 1, 由XiX_i的定义可知XiX_iXi1X_{i-1}ξi\xi_i的函数,如此递归代入后可知XiX_i是随机变量ξ1,,ξi\xi_1, \cdots, \xi_i的函数。这表明事件

{Xn=in,,X1=i1,X0=K}\{X_n = i_n, \cdots, X_1 = i_1, X_0 = K\}

由随机变量ξ1,,ξn\xi_1, \cdots, \xi_n确定。再由{ξi,i1}\{ \xi_i, i \geq 1 \}的独立性,

P(Xn+1=in+1Xn=in,,X1=i1,X0=K)={P((Kξn+1)0=in+1),ink,P((inξn+1)0=in+1),k<inK,={P((Kξn+1)0=in+1Xn=in),ink,P((inξn+1)0=in+1Xn=in),k<inK,=P(Xn+1=in+1Xn=in).\begin{aligned} &P(X_{n+1} = i_{n+1} | X_n = i_n, \cdots, X_1 = i_1, X_0 = K) \\ &=\begin{cases} P((K - \xi_{n+1}) \lor 0 = i_{n+1}), & i_n \leq k, \\ P((i_n - \xi_{n+1}) \lor 0 = i_{n+1}), & k < i_n \leq K, \end{cases}\\ &=\begin{cases} P((K - \xi_{n+1}) \lor 0 = i_{n+1} | X_n = i_n), & i_n \leq k, \\ P((i_n - \xi_{n+1}) \lor 0 = i_{n+1} | X_n = i_n), & k < i_n \leq K, \end{cases}\\ &= P(X_{n+1} = i_{n+1} | X_n = i_n). \end{aligned}

所以XX为离散时间马氏链。

(2) 设P(ξn=0)=P(ξn=2)=1/4, P(ξn=1)=1/2P(\xi_n = 0) = P(\xi_n = 2) = 1/4,\ P(\xi_n = 1) = 1/2。当k=1,K=3k = 1, K = 3时,写出XX的转移概率矩阵。

点击查看答案

(2)当Xn1=0X_{n-1}=0时,Xn=3ξnX_{n}=3-\xi_n,故矩阵第一行为

(0141214)\left(0\quad\frac{1}{4}\quad\frac{1}{2}\quad\frac{1}{4}\right)

类似可得到剩下三行转移概率(请读者自行推导),最终得到转移概率矩阵为

P=(0141214014121414121400141214).P=\begin{pmatrix} 0&\dfrac{1}{4}&\dfrac{1}{2}&\dfrac{1}{4}\\[10pt] 0&\dfrac{1}{4}&\dfrac{1}{2}&\dfrac{1}{4}\\[10pt] \dfrac{1}{4}&\dfrac{1}{2}&\dfrac{1}{4}&0\\[10pt] 0&\dfrac{1}{4}&\dfrac{1}{2}&\dfrac{1}{4} \end{pmatrix}.

(3) 在(2)的条件下求X3X_3的分布。

点击查看答案

(3)初始分布为

π(0)=(0001)\pi^{(0)}=(0\quad 0\quad 0\quad 1)

X3X_3分布为

π(3)=π(0)P3=(332113238316).\pi^{(3)}=\pi^{(0)}P^3=\left(\frac{3}{32}\quad\frac{11}{32}\quad\frac{3}{8}\quad\frac{3}{16}\right).

XX是不可约dd周期马氏链,状态空间为SS。对任意iSi \in S,令

Sk={jS;pi,j(nd+k)>0,n0},k=0,,d1.S_k = \{ j \in S; \, p_{i,j}^{(nd+k)} > 0, \, n \geq 0\}, \quad k = 0, \dots, d-1.

证明:

(1)Sk,0kd1S_k, 0 \leq k \leq d-1,互不相交而且S=k=0d1Sk\displaystyle S = \bigcup_{k=0}^{d-1} S_k

点击查看答案

(1)先证明S=k=0d1Sk\displaystyle S = \bigcup_{k=0}^{d-1} S_k:显然k=0d1SkS\displaystyle \bigcup_{k=0}^{d-1} S_k\subseteq S,而对任意jSj\in S,因为XX不可约,所以存在正整数mm使pi,j(m)>0p_{i,j}^{(m)}>0.
那么取km(modd)k\equiv m(\mathrm{mod}\hspace{0.2em}d),则jSkj\in S_k,故Sk=0d1Sk\displaystyle S\subseteq\bigcup_{k=0}^{d-1} S_k,所以S=k=0d1Sk\displaystyle S = \bigcup_{k=0}^{d-1} S_k
再证明SkS_k互不相交:若存在k1k2k_1\neq k_2使Sk1Sk1S_{k_1}\cap S_{k_1}\neq\varnothing,则存在jSkSlj\in S_k\cap S_l,于是存在n1,n2Nn_1,n_2\in\N,使

pi,j(n1d+k1)>0,pi,j(n2d+k2)>0.p_{i,j}^{(n_1d+k_1)}>0,\,\quad p_{i,j}^{(n_2d+k_2)}>0.

又因为XX不可约,所以jij\to i,存在lNl\in N使pj,i(l)>0p_{j,i}^{(l)}>0,那么

pii(n1d+k1+l)pi,j(n1d+k1)pj,i(l)>0,pii(n2d+k2+l)pi,j(n2d+k2)pj,i(l)>0\begin{aligned} p_{ii}^{(n_1d+k_1+l)}\geq p_{i,j}^{(n_1d+k_1)}p_{j,i}^{(l)}>0,\\ p_{ii}^{(n_2d+k_2+l)}\geq p_{i,j}^{(n_2d+k_2)}p_{j,i}^{(l)}>0 \end{aligned}

所以dn1d+k1+l,dn2d+k2+ld(n1n2)d+(k1k2)d\mid n_1d+k_1+l,d\mid n_2d+k_2+l\Longrightarrow d\mid (n_1-n_2)d+(k_1-k_2),于是d(k1k2)d\mid (k_1-k_2)
k1,k2{0,1,,d1}k_1,k_2\in\{0,1,\cdots,d-1\},故k1k2<d|k_1-k_2|< d,矛盾。故SkS_k互不相交。

(2)令Yn=XndY_n = X_{nd},那么Y={Yn;n0}Y = \{ Y_n; \, n \geq 0 \}是以P(d)P^{(d)}为一步转移概率矩阵的非周期马氏链。特别地,若X0SkX_0 \in S_k,那么YY是状态空间为SkS_k的非周期马氏链。

点击查看答案

(2)证明分三步:

  1. 证明YY是马氏链,且转移概率矩阵为P(d)P^{(d)}
    XX的马氏性,
P{Yn+1=lY0,Y1,,Yn=j}=P{X(n+1)d=lX0,Xd,,Xnd=j}=P{Xnd+d=lXnd=j}=pjl(d).\begin{aligned} P\{Y_{n+1} = l \mid Y_0, Y_1, \dots, Y_n = j\} &= P\{X_{(n+1)d} = l \mid X_0, X_d, \dots, X_{nd} = j\}\\ &= P\{X_{nd+d} = l \mid X_{nd} = j\}\\ &= p_{jl}^{(d)}. \end{aligned}

所以YY是马氏链,且一步转移概率矩阵为P(d)P^{(d)}
2. 证明YY是非周期马氏链
对任意状态jSj\in S,记YYjj的周期为dY(j)=gcd{m1:pjj(md)>0}d_Y(j)=\gcd\{m\geq1:p_{jj}^{(md)}>0\},而在XX中如果jj存在周期,则一定存在sNs\in\N,使pjj(sd)>0p_{jj}^{(sd)}>0,所以dY(j)=1d_Y(j)=1YY是非周期马氏链。
3. 证明若X0SkX_0 \in S_k,那么YY是状态空间为SkS_k的非周期马氏链
uSau \in S_a,且 puv(t)>0p_{uv}^{(t)} > 0t=sd+bt = sd + b,则 vS(a+b)moddv \in S_{(a+b) \bmod d}
YnY_n的一步转移对应原链走 dd步:t=d=1d+0t = d = 1 \cdot d + 0,余数 b=0b = 0
若初态 X0=iSkX_0 = i \in S_k,则

uSk,puv(d)>0    vS(k+0)modd=Sk.u \in S_k,\quad p_{uv}^{(d)} > 0 \implies v \in S_{(k+0) \bmod d} = S_k.

即转移矩阵 P(d)P^{(d)}仅在 SkS_k内部转移,不会跑到其余 Sl(lk)S_l (l \ne k)
又原链不可约,故对任意 u,vSku, v \in S_k,存在 mm使得 puv(md)>0p_{uv}^{(md)} > 0,因此 P(d)P^{(d)}SkS_k上不可约。
所以YY是状态空间为SkS_k的非周期马氏链。

XX是马氏链,X0=i,τi=inf{n>0;Xn=i},τi(2)=inf{n>τi;Xn=i}X_0 = i, \tau_i =\inf\{ n > 0; X_n = i\},\tau_i^{(2)} = \inf\{n > \tau_i; X_n=i\}。证明τi\tau_iτi(2)τi\tau_i^{(2)}-\tau_i独立且同分布。

点击查看答案

对任意正整数m,nm,n

P(τi=m,τi(2)τi=n)=P(X1i,X2i,,Xm1i,Xm=i,Xm+1i,,Xm+n1i,Xm+n=i)=P(Xm+1i,,Xm+n1i,Xm+n=iX1i,X2i,,Xm1i,Xm=i)P(X1i,X2i,,Xm1i,Xm=i)=P(Xm+1i,,Xm+n1i,Xm+n=iXm=i)P(τi=m)=P(X1i,,Xn1i,Xn=iX0=i)P(τi=m)=Pi(τi=n)P(τi=m)\begin{aligned} P(\tau_i=m,\tau_i^{(2)}-\tau_i=n)&=P(X_1\neq i,X_2\neq i,\cdots,X_{m-1}\neq i,X_{m}=i,X_{m+1}\neq i,\cdots,X_{m+n-1}\neq i,X_{m+n}=i)\\ &=P(X_{m+1}\neq i,\cdots,X_{m+n-1}\neq i,X_{m+n}=i|X_1\neq i,X_2\neq i,\cdots,X_{m-1}\neq i,X_{m}=i)P(X_1\neq i,X_2\neq i,\cdots,X_{m-1}\neq i,X_{m}=i)\\ &=P(X_{m+1}\neq i,\cdots,X_{m+n-1}\neq i,X_{m+n}=i|X_{m}=i)P(\tau_i=m)\\ &=P(X_{1}\neq i,\cdots,X_{n-1}\neq i,X_{n}=i|X_{0}=i)P(\tau_i=m)\\ &=P_i(\tau_i=n)P(\tau_i=m) \end{aligned}

其中Pi(τi=n)P_i(\tau_i=n)表示从ii出发首次返回ii的时间为nn的概率,恰等于P(τi(2)τi=n)P(\tau_i^{(2)}-\tau_i=n),所以τi\tau_iτi(2)τi\tau_i^{(2)}-\tau_i独立同分布。

在马氏链中,如果ii常返,且iji\rightarrow j,则jj常返。

证明:若iji\rightarrow j,则存在正整数mm,使pi,j(m)>0p_{i,j}^{(m)}>0。又因为ii常返,所以k=1pi,i(k)=\displaystyle\sum_{k=1}^\infty p_{i,i}^{(k)}=\infty,且jij\rightarrow i(否则从jj出发无法到达ii,说明ii非常返)。
所以存在正整数nn,使pj,i(n)>0p_{j,i}^{(n)}>0。于是

k=0pj,j(k)k=m+npj,i(n)pi,i(knm)pi,j(m)=pj,i(n)pi,j(m)k=0pi,i(k)=+,\sum_{k=0}^\infty p_{j,j}^{(k)} \geq \sum_{k=m+n}^\infty p_{j,i}^{(n)} p_{i,i}^{(k-n-m)} p_{i,j}^{(m)} = p_{j,i}^{(n)} p_{i,j}^{(m)} \sum_{k=0}^\infty p_{i,i}^{(k)} = +\infty,

所以jj常返。

在马氏链中,如果ii非常返,且iji\rightarrow j,则jj非常返。

一个简单的反例:状态空间{i,j}\{i,j\},转移概率矩阵(0.50.501)\begin{pmatrix}0.5 &0.5\\0 &1\end{pmatrix},则ii非常返且iji\rightarrow j,但jj是常返状态(吸收态)。

假设马氏链X={Xn,n0}X = \{X_n, n \geq 0\}的状态空间为{1,2,3}\{1, 2, 3\},转移概率矩阵为

(0.70.20.10.30.50.2001)\begin{pmatrix} 0.7 &0.2 &0.1\\ 0.3 &0.5 &0.2\\ 0 &0 &1 \end{pmatrix}

τ3=inf{n>0,Xn=3}\tau_3 = \inf\{n > 0, \, X_n = 3\},求P(X3=1X0=1,τ3>3)P(X_3 = 1 | X_0 = 1, \tau_3 > 3)

点击查看答案

注意到状态33是吸收态,所以

P(X3=1X0=1,τ3>3)=P(X3=1,X0=1,τ3>3)P(X0=1,τ3>3)=P(X3=1,X0=1)P(X0=1,X33)=p1,1(3)p1,1(3)+p1,2(3)=0.4570.457+0.230.665.\begin{aligned} P(X_3 = 1 | X_0 = 1, \tau_3 > 3)&=\frac{P(X_3 = 1,X_0 = 1, \tau_3 > 3)}{P(X_0 = 1, \tau_3 > 3)}\\ &=\frac{P(X_3 = 1,X_0 = 1)}{P(X_0 = 1, X_3\neq 3)}\\ &=\frac{p_{1,1}^{(3)}}{p_{1,1}^{(3)}+p_{1,2}^{(3)}}\\ &=\frac{0.457}{0.457+0.23}\approx 0.665. \end{aligned}

此题关键在于注意到吸收态,由此简化首访条件。

记马氏链XX的状态空间为S={1,2,3,4}S=\{1,2,3,4\},转移概率矩阵为

(14141414121316001316121601213)\begin{pmatrix} \dfrac{1}{4} & \dfrac{1}{4} & \dfrac{1}{4} & \dfrac{1}{4} \\[10pt] \dfrac{1}{2} & \dfrac{1}{3} & \dfrac{1}{6} & 0\\[10pt] 0 & \dfrac{1}{3} & \dfrac{1}{6} & \dfrac{1}{2} \\[10pt] \dfrac{1}{6} & 0 & \dfrac{1}{2} & \dfrac{1}{3} \end{pmatrix}

τi=inf{n1;Xn=i}\tau_i = \inf\{n \geq 1; \, X_n = i\},求P1(τ3<τ4)P_1(\tau_3 < \tau_4)

点击查看答案

进行前进一步分解:设X0,X1X_0,X_1分别表示出发状态与第一步到达状态,则

P1(τ3<τ4)=P(τ3<τ4X0=1)=k=14P1(τ3<τ4X1=k)P1(X1=k)=14P1(τ3<τ4)+14P2(τ3<τ4)+141+140\begin{aligned} P_1(\tau_3 < \tau_4)&=P(\tau_3 < \tau_4|X_0=1)\\ &=\sum_{k=1}^4P_1(\tau_3 < \tau_4|X_1=k)P_1(X_1=k)\\ &=\frac{1}{4}P_1(\tau_3 < \tau_4)+\frac{1}{4}P_2(\tau_3 < \tau_4)+\frac{1}{4}\cdot 1+\frac{1}{4}\cdot 0 \end{aligned}

于是我们需要求P2(τ3<τ4)P_2(\tau_3< \tau_4)。同理,

P2(τ3<τ4)=k=14P2(τ3<τ4X1=k)P2(X1=k)=12P1(τ3<τ4)+13P2(τ3<τ4)+161\begin{aligned} P_2(\tau_3 < \tau_4)&=\sum_{k=1}^4P_2(\tau_3 < \tau_4|X_1=k)P_2(X_1=k)\\ &=\frac{1}{2}P_1(\tau_3 < \tau_4)+\frac{1}{3}P_2(\tau_3 < \tau_4)+\frac{1}{6}\cdot 1 \end{aligned}

联立这两个方程,解得

{P1(τ3<τ4)=59P2(τ3<τ4)=23.\begin{cases} \boxed{\displaystyle P_1(\tau_3 < \tau_4)=\frac{5}{9}}\\[10pt] \displaystyle P_2(\tau_3 < \tau_4)=\frac{2}{3}. \end{cases}

设不可约正常返马氏链XX的转移概率矩阵为(pi,j)i,jS,π=(πi,iS)(p_{i,j})_{i,j \in S}, \pi = (\pi_i, i \in S)为平稳分布。
Yn=(Xn,Xn+1)Y_n = (X_n, X_{n+1}),证明Y={Yn;n0}Y = \{Y_n; n \geq 0\}是状态空间为M={(i,j);i,jS,pi,j>0}\mathcal{M} = \{(i,j); i,j \in S, p_{i,j} > 0\}的不可约正常返马氏链,平稳分布为(πipi,j,(i,j)M)(\pi_ip_{i,j}, (i,j) \in \mathcal{M})

点击查看答案
  1. YY是马氏链
    对任意i,j,k,lSi,j,k,l\in S满足(i,j),(k,l)M(i,j),(k,l)\in\mathcal{M}
P(Yn+1=(k,l)Yn=(i,j))={pj,l,k=j0,kjP(Y_{n+1}=(k,l)|Y_n=(i,j))=\begin{cases} p_{j,l},&k=j\\ 0,&k\neq j \end{cases}

因此转移概率q(i,j),(k,l)=1{k=j}pj,lq_{(i,j),(k,l)}=\mathbf{1}_{\{k=j\}}p_{j,l}YnY_n之前的状态无关,故YY是马氏链。
2. YY是不可约的
因为XX是不可约的,所以任意两个状态j,kSj,k\in S,有jkj\rightarrow k,即存在一条路径j=x0x1xm1xm=kj=x_0\rightarrow x_1\rightarrow\cdots\rightarrow x_{m-1}\rightarrow x_m=k使pj,x1,,pxm1,k>0p_{j,x_1},\cdots,p_{x_{m-1},k}>0
因此对任意i,j,k,lS,(i,j),(k,l)Mi,j,k,l\in S,(i,j),(k,l)\in\mathcal{M}q(i,j)(k,l)(m+1)q(i,j)(j,x1)q(xm1,k)(k,l)>0q_{(i,j)(k,l)}^{(m+1)}\geq q_{(i,j)(j,x_1)}\cdots q_{(x_{m-1},k)(k,l)}>0,所以(i,j)(k,l)(i,j)\rightarrow (k,l)YY是不可约的。
3. YY存在平稳分布
定义π~(i,j)πipi,j,(i,j)M\widetilde{\pi}_{(i,j)}\coloneqq \pi_i p_{i,j}, (i,j) \in \mathcal{M},则

(i,j)Mπ~(i,j)=iSjSπipi,j=iSπijSpi,j=iSπi=1\sum_{(i,j) \in \mathcal{M}} \widetilde{\pi}_{(i,j)} = \sum_{i \in S} \sum_{j \in S} \pi_i p_{i,j} = \sum_{i \in S} \pi_i \sum_{j \in S} p_{i,j} = \sum_{i \in S} \pi_i = 1

所以π~\widetilde{\pi}是概率分布。又因为对任意(k,l)M(k,l) \in \mathcal{M},

(i,j)Mπ~(i,j)q(i,j),(k,l)=i,jπipi,j1{j=k}pk,l=(iπipi,k)pk,l.\sum_{(i,j) \in \mathcal{M}} \widetilde{\pi}_{(i,j)} q_{(i,j),(k,l)} = \sum_{i,j} \pi_i p_{i,j} \mathbf{1}_{\{j=k\}} p_{k,l} = \left( \sum_i \pi_i p_{i,k} \right) p_{k,l}.

由于π\piXX的平稳分布,有

iπipi,k=πk\sum_i \pi_i p_{i,k}=\pi_k

于是

(i,j)Mπ~(i,j)q(i,j),(k,l)=πkpk,l=π~(k,l).\sum_{(i,j) \in \mathcal{M}} \widetilde{\pi}_{(i,j)} q_{(i,j),(k,l)} = \pi_k p_{k,l} = \widetilde{\pi}_{(k,l)}.

π~\widetilde{\pi}YY的平稳分布。
4. YY是正常返的
因为YY是不可约马氏链,且有平稳分布,所以YY一定是不可约正常返马氏链。

综上,命题得证。

证明:不可约马尔可夫链XX正常返的充要条件是存在非负数列{zi;iS}\{z_i; i \in S\}及一个状态jj使得对任意iji \neq jkSpi,kzkzi1\displaystyle\sum_{k \in S} p_{i,k} z_k \leq z_i - 1而且kSpj,kzk<\displaystyle\sum_{k \in S} p_{j,k} z_k < \infty

点击查看答案

证明:

  1. 必要性:若 XX 不可约正常返,对任意 iji \neq j,令 zi=mi,j,zj=0z_i = m_{i,j}, z_j = 0,那么任意 iji \neq j
zi=mi,j=kjpi,kmk,j+1=kSpi,kzk+1z_i = m_{i,j} = \sum_{k \neq j} p_{i,k} m_{k,j} + 1 = \sum_{k \in S} p_{i,k} z_k + 1

kSpj,kzk=kjpj,kmk,jmj,j<.\sum_{k \in S} p_{j,k} z_k = \sum_{k \neq j} p_{j,k} m_{k,j} \leq m_{j,j} < \infty.
  1. 充分性:记δ=kSpj,kzk\displaystyle\delta = \sum_{k \in S} p_{j,k} z_k,则0<δ<0 < \delta < \infty。任取iji \neq j,对任意n2n \geq 2,迭代得,
0kpi,k(n)zk=kSlSpi,l(n1)pl,kzk=lSpi,l(n1)kSpl,kzkljpi,l(n1)(zl1)+pi,j(n1)δlSpi,l(n1)zl+pi,j(n1)(1+δ)1zi+k=1n1pi,j(k)(1+δ)n.\begin{aligned} 0 \leq \sum_k p_{i,k}^{(n)} z_k &= \sum_{k \in S} \sum_{l \in S} p_{i,l}^{(n-1)} p_{l,k} z_k = \sum_{l \in S} p_{i,l}^{(n-1)} \sum_{k \in S} p_{l,k} z_k\\ &\leq \sum_{l \neq j} p_{i,l}^{(n-1)} (z_l - 1) + p_{i,j}^{(n-1)} \delta\\ &\leq \sum_{l \in S} p_{i,l}^{(n-1)} z_l + p_{i,j}^{(n-1)} (1 + \delta) - 1\\ &\leq\cdots\leq z_i+\sum_{k=1}^{n-1}p_{i,j}^{(k)}(1+\delta)-n. \end{aligned}

这表明

k=1n1pi,j(k)(1+δ)nzi.\sum_{k=1}^{n-1} p_{i,j}^{(k)} (1 + \delta) \geq n - z_i.

两边除以nn并令nn \to \infty

limn1nk=1npi,j(k)11+δ>0.\lim_{n \to \infty} \frac{1}{n} \sum_{k=1}^n p_{i,j}^{(k)} \geq \frac{1}{1 + \delta} > 0.

因此jj为正常返状态,从而XX是正常返的。

实际上,这道题就是Foster-Lyapunov定理的一个特殊形式,其在稳定性理论中起着重要作用。(Lyapunov这个名字似乎在概率论中也出现过……)

(1)设不可约非常返或零常返马氏链XX的状态空间S={0,1,2,}S = \{0,1,2,\cdots\},证明:若数列{yk}\{y_k\}满足

limk+yk=+\lim_{k \to +\infty} y_k = +\infty

那么对任意iSi \in S

limnkSpi,k(n)yk=+.\lim_{n \to \infty} \sum_{k \in S} p_{i,k}^{(n)} y_k = +\infty.
点击查看答案

证明:任意kSk \in S,由 X 非常返或零常返可知,pi,k(n)0p_{i,k}^{(n)} \to 0。注意到

kSpi,k(n)yk=k<Npi,k(n)yk+kNpi,k(n)yk(1k<Npi,k(n))minsNys.\sum_{k \in S} p_{i,k}^{(n)} y_k = \sum_{k < N} p_{i,k}^{(n)} y_k + \sum_{k \geq N} p_{i,k}^{(n)} y_k \geq \left(1 - \sum_{k < N} p_{i,k}^{(n)}\right) \min_{s \geq N} y_s.

先令nn \to \infty,可得limnk<Npi,k(n)=0\displaystyle\lim_{n \to \infty} \sum_{k < N} p_{i,k}^{(n)} = 0,再令NN \to \infty,可知此时结论成立。

(2)设不可约链XX的状态空间S={0,1,2,}S = \{0,1,2,\cdots\},若存在正数{yi;iS}\{y_i; i \in S\}及某个状态jj使得对任意iji \neq j,

kSpi,kykyi\sum_{k \in S} p_{i,k} y_k \leq y_i

而且limi+yi=+\displaystyle\lim_{i \to +\infty} y_i = +\infty,证明XX常返。

点击查看答案

证明:构造一个新马氏链YY使其转移概率为p~i,k=pi,k,ij,kS\tilde{p}_{i,k} = p_{i,k}, i \neq j, k \in S,

p~j,j=δ+(1δ)pj,j,p~j,k=(1δ)pj,k, kj.\tilde{p}_{j,j} = \delta + (1 - \delta)p_{j,j},\quad \tilde{p}_{j,k} = (1 - \delta)p_{j,k},\ k \neq j.

由于XX是不可约的,因此YY也是不可约的。对任意iji \neq jkSp~i,kykyi\displaystyle\sum_{k \in S} \tilde{p}_{i,k} y_k \leq y_i,而且Δ=kSp~j,kyk>0\displaystyle\Delta = \sum_{k \in S} \tilde{p}_{j,k} y_k > 0。注意到

kSp~i,k(n)yk=kSlSp~i,l(n1)p~l,kyk=lSp~i,l(n1)kSp~l,kyk=ljp~i,l(n1)kSpl,kyk+Δp~i,j(n1)ljp~i,l(n1)yl+Δp~i,j(n1)yi+Δk=1n1p~i,j(k).\begin{aligned} \sum_{k \in S} \tilde{p}_{i,k}^{(n)} y_k &= \sum_{k \in S} \sum_{l \in S} \tilde{p}_{i,l}^{(n-1)} \tilde{p}_{l,k} y_k \\ &= \sum_{l \in S} \tilde{p}_{i,l}^{(n-1)} \sum_{k \in S} \tilde{p}_{l,k} y_k\\ &=\sum_{l \neq j} \tilde{p}_{i,l}^{(n-1)} \sum_{k \in S} p_{l,k} y_k + \Delta \tilde{p}_{i,j}^{(n-1)}\\ &\leq \sum_{l \neq j} \tilde{p}_{i,l}^{(n-1)} y_l + \Delta \tilde{p}_{i,j}^{(n-1)} \\ &\leq \cdots \leq y_i + \Delta \sum_{k=1}^{n-1} \tilde{p}_{i,j}^{(k)}. \end{aligned}

YY是非常返的,那么令nn \to \infty,由上一问结论知k=1n1p~i,j(k)+\displaystyle\sum_{k=1}^{n-1} \tilde{p}_{i,j}^{(k)} \to +\infty,因此k=0p~i,j(k)=+\displaystyle\sum_{k=0}^{\infty} \tilde{p}_{i,j}^{(k)} = +\infty,这表明YY的状态jj是常返的,从而YY是常返的,矛盾。
因此YY一定是常返的,故方程组

zi=ljp~i,lzl+p~i,j=ljpi,lzl+pi,j,ijz_i = \sum_{l \neq j} \tilde{p}_{i,l} z_l + \tilde{p}_{i,j} = \sum_{l \neq j} p_{i,l} z_l + p_{i,j},\quad i \neq j

的最小非负解恒为11,由此可知XX是常返的。

当非负整数值马氏链X={Xn;n0}X = \{X_n; n \geq 0\}的转移概率满足:
(1)

{p0,j=pj,j0pi,j=pji+1,i1,ji1j=0pj=1,j=0jpj=m<1,0<p0,p0+p1<1\begin{cases} p_{0,j} = p_j, \quad j \geq 0\\ p_{i,j} = p_{j-i+1}, \quad i \geq 1,\quad j \geq i-1\\ \displaystyle\sum_{j=0}^{\infty} p_j = 1, \sum_{j=0}^{\infty} j p_j = m < 1, \\ 0 < p_0, p_0 + p_1 < 1 \end{cases}

(2)

{pi,0=1k=0ipk,i0pi,j=pi+1j,j1,ij1j=0pj=1,j=0jpj=m>1,0<p0<1\begin{cases} \displaystyle p_{i,0} = 1-\sum_{k=0}^i p_k, \quad i \geq 0\\ p_{i,j} = p_{i+1-j}, \quad j \geq 1,\quad i \geq j-1\\ \displaystyle\sum_{j=0}^{\infty} p_j = 1, \sum_{j=0}^{\infty} j p_j = m > 1, \\ 0 < p_0 < 1 \end{cases}

分别讨论XX的常返性(正常返/零常返/非常返)。

点击查看第1问答案

(1)由条件0<p0,p0+p1<10 < p_0, p_0 + p_1 < 1可知XX是不可约的。转移概率矩阵形式为

(p0p1p2p3p4p0p1p2p3p40p0p1p2p300p0p1p2000p0p1)\begin{pmatrix} p_0 & p_1 & p_2 & p_3 & p_4 & \cdots \\ p_0 & p_1 & p_2 & p_3 & p_4 & \cdots \\ 0 & p_0 & p_1 & p_2 & p_3 & \cdots \\ 0 & 0 & p_0 & p_1 & p_2 & \cdots \\ 0 & 0 & 0 & p_0 & p_1 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}

解法1:注意到这个马氏链(看作随机游动)是左连续的,设di=Ei[τi1],i1d_i=\mathbb{E}_i[\tau_{i-1}],i\geq 1,即从状态ii出发最终到达状态i1i-1的平均时间。那么可得到关系式

di=p01+k=1pk(1+di1+k+di2+k++di)d_i=p_0\cdot 1+\sum_{k=1}^\infty p_k(1+d_{i-1+k}+d_{i-2+k}+\cdots+d_i)

又因为对每个状态i1i\geq 1,其向下一步的结构相同,故所有did_i相同(记为dd),代入得

d=k=0pk+dk=1kpk=1+mdd=11md=\sum_{k=0}^\infty p_k+d\sum_{k=1}^\infty kp_k=1+md\Longrightarrow d=\frac{1}{1-m}

然后再计算E0[τ0]\mathbb{E}_0[\tau_0],即从00出发最终回到00的时间期望:

E0[τ0]=p01+k=1kdpk=1+dm=11m<.\mathbb{E}_0[\tau_0]=p_0\cdot 1+\sum_{k=1}^\infty kdp_k=1+dm=\frac{1}{1-m}< \infty.

所以00是正常返状态。又XX不可约,故XX是正常返的。


解法2:对任意k0k \geq 0,构造zk=k1mz_k = \dfrac{k}{1-m}。那么对任意i0i \neq 0

kSpi,kzk=11mkSpi,kk=11m(i1+ki1pi,k(k(i1)))=11m(i1+m)=i1m1=zi1,kSp0,kzk=11mkSp0,kk=m1m<+.\begin{aligned} \sum_{k \in S} p_{i,k} z_k &= \frac{1}{1-m} \sum_{k \in S} p_{i,k} k \\ &= \frac{1}{1-m} \left( i - 1 + \sum_{k \geq i-1} p_{i,k} \left( k - (i-1) \right) \right)\\ &= \frac{1}{1-m} \left( i - 1 + m \right) = \frac{i}{1-m} - 1 = z_i - 1,\\ \sum_{k \in S} p_{0,k} z_k &= \frac{1}{1-m} \sum_{k \in S} p_{0,k} k = \frac{m}{1-m} < +\infty. \end{aligned}

由上上题结论(Foster-Lyapunov定理)可知XX是正常返的。

点击查看第2问答案

(2)转移概率矩阵形式为

(1p0p00001p0p1p1p0001k=02pkp2p1p001k=03pkp3p2p1p0)\begin{pmatrix} 1-p_0 & p_0 & 0 & 0 & 0 & \cdots \\ 1-p_0-p_1 & p_1 & p_0 & 0 & 0 & \cdots \\ \displaystyle 1-\sum_{k=0}^{2}p_k & p_2 & p_1 & p_0 & 0 & \cdots \\ \displaystyle 1-\sum_{k=0}^{3}p_k & p_3 & p_2 & p_1 & p_0 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}

η(s)=j=0pjsj\displaystyle \eta(s) = \sum_{j=0}^{\infty} p_j s^j,则由j=0jpj=m>1\displaystyle \sum_{j=0}^{\infty} j p_j = m > 1可知η(s)=s\eta(s) = s(0,1)(0,1)内存在唯一解qq。对任意i0i \geq 0,令πi=(1q)qi>0\pi_i = (1-q)q^i > 0,那么

i=0πi=1,i=0πipi,0=i=0(1k=0ipk)πi=1k=0pki=k(1q)qi=1k=0pkqk=1η(q)=1q=π0,\begin{aligned} \sum_{i=0}^{\infty} \pi_i &= 1,\\ \sum_{i=0}^{\infty} \pi_i p_{i,0}&= \sum_{i=0}^{\infty} \left(1 - \sum_{k=0}^{i} p_k\right) \pi_i\\ &= 1 - \sum_{k=0}^{\infty} p_k \sum_{i=k}^{\infty} (1-q)q^i= 1 - \sum_{k=0}^{\infty} p_k q^k\\ &= 1 - \eta(q)= 1 - q= \pi_0, \end{aligned}

且对任意j0j \neq 0

i=0πipi,j=i=j1pij+1(1q)qi=(1q)qj1i=j1pij+1qij+1=(1q)qj1k=0pkqk=(1q)qj1η(q)=(1q)qj=πj.\begin{aligned} \sum_{i=0}^{\infty} \pi_i p_{i,j}&= \sum_{i=j-1}^{\infty} p_{i-j+1} (1-q)q^i\\ &= (1-q)q^{j-1} \sum_{i=j-1}^{\infty} p_{i-j+1} q^{i-j+1}= (1-q)q^{j-1} \sum_{k=0}^{\infty} p_k q^k\\ &= (1-q)q^{j-1} \eta(q)= (1-q)q^j= \pi_j. \end{aligned}

综上可知XX存在平稳分布,又由条件0<p0<10 < p_0 < 1可知XX是不可约的,因此XX是正常返的。

那么这一问里面平稳分布πi=(1q)qi\pi_i = (1-q)q^i是如何想到的呢?根据条件,平稳分布需满足

πj=i=0πipi,j=i=j1πipi,j=i=j1πipi+1j=k=0πk+j1pk\pi_j=\sum_{i=0}^\infty\pi_ip_{i,j}=\sum_{i=j-1}^\infty\pi_ip_{i,j}=\sum_{i=j-1}^\infty\pi_ip_{i+1-j}=\sum_{k=0}^\infty\pi_{k+j-1}p_{k}

为了将变量kkjj分离,考虑πk+j1=sk+j1\pi_{k+j-1}=s^{k+j-1},则πj=sj1k=0skpk=sj1η(s)\displaystyle\pi_j=s^{j-1}\sum_{k=0}^\infty s_{k}p_k=s^{j-1}\eta(s),因此只需要η(s)=s\eta(s)=s即可。【又因为m>1m>1,所以ss一定有解】

假设XX是一个不可约马氏链,状态空间为(0,1,2,,N)(0,1,2,\cdots,N),转移概率pi,jp_{i,j}满足条件jSjpi,j=i\displaystyle\sum_{j \in S} j p_{i,j} = i,令qiq_iXXii出发到达00之前先到达NN的概率,求qiq_i

点击查看答案

容易验证{qi}\{q_i\}满足方程组

qi=jSpi,jqji>0q_i = \sum_{j \in S} p_{i,j}q_j \quad i > 0

而且是该方程组的最小非负解。此外由条件可知{iN}\displaystyle\left\{\frac{i}{N}\right\}也是该方程组的非负解。
qiiNq_i \leq \dfrac{i}{N}且存在某个i0i_0使得qi0<i0Nq_{i_0} < \dfrac{i_0}{N},那么由不可约性,存在k1k \geq 1使得pN,i0(k)>0p_{N,i_0}^{(k)} > 0。由此可得

1=qN=jSqjpN,j(k)<jSjNpN,j(k)=1NjSjpN,j(k)=1NiSpN,i(k1)jSjpi,j=11 = q_N = \sum_{j \in S} q_j p_{N,j}^{(k)}< \sum_{j \in S} \frac{j}{N} p_{N,j}^{(k)} = \frac{1}{N} \sum_{j \in S} j p_{N,j}^{(k)} = \frac{1}{N} \sum_{i \in S} p_{N,i}^{(k-1)} \sum_{j \in S} j p_{i,j} = 1

矛盾表明qi=iNq_i = \dfrac{i}{N}对所有iSi \in S都成立。


关于上述推导中第二个等号和最后一个等号的解释:

qN=j1qj1pN,j1=j1j2qj2pj1,j2pN,j1=j2pN,j2(2)==jpN,j(k),iSpN,i(k1)jSjpi,j=iSipN,i(k1)=iSi(lSpN,l(k2)pl,i)=lSpN,l(k2)(iSipl,i)=lSlpN,l(k2)==lSlpN,l(0)=N\begin{aligned} \begin{aligned} q_N&=\sum_{j_1}q_{j_1}p_{N,j_1}=\sum_{j_1}\sum_{j_2}q_{j_2}p_{j_1,j_2}p_{N,j_1}\\ &=\sum_{j_2}p_{N,j_2}^{(2)}=\cdots=\sum_{j}p_{N,j}^{(k)}, \end{aligned}\\ \begin{aligned} \sum_{i \in S} p_{N,i}^{(k-1)} \sum_{j \in S} j p_{i,j}&=\sum_{i \in S} ip_{N,i}^{(k-1)}=\sum_{i \in S}i\left(\sum_{l\in S}p_{N,l}^{(k-2)}\cdot p_{l,i}\right)\\ &=\sum_{l \in S}p_{N,l}^{(k-2)}\left(\sum_{i\in S}i\cdot p_{l,i}\right)\\ &=\sum_{l \in S}lp_{N,l}^{(k-2)}=\cdots=\sum_{l \in S}lp_{N,l}^{(0)}=N \end{aligned} \end{aligned}

马氏链简单应用

设 G-W 分支过程X={Xn;n0}X=\{X_n; n \geq 0\}的后代均值为m<1m < 1,后代方差为σ2\sigma^2,假定X0=1X_0 = 1,求存活过粒子数的方差。

点击查看答案

设存活过粒子数为S=n=0Xn\displaystyle S=\sum_{n=0}^\infty X_n,则因为m<1m< 1,所以ES=11m\mathbb{E}S=\dfrac{1}{1-m}。而根据条件方差公式:

Var(S)=Var[E(SY1,1)]+E[Var(SY1,1)]\mathrm{Var}(S)=\mathrm{Var}[\mathbb{E}(S|Y_{1,1})]+\mathbb{E}[\mathrm{Var}(S|Y_{1,1})]

又因为S=X0+i=1Y1,1Si=1+i=1Y1,1Si\displaystyle S=X_0+\sum_{i=1}^{Y_{1,1}} S_i=1+\sum_{i=1}^{Y_{1,1}} S_i,故

Var[E(SY1,1)]=Var[E(1+i=1Y1,1SiY1,1)]=Var(1+Y1,1ES)=Var(Y1,1)(1m)2=σ2(1m)2,E[Var(SY1,1)]=E[Var(1+i=1Y1,1SiY1,1)]=E[i=1Y1,1Var(Si)]=E(Y1,1)Var(S)=mVar(S).\begin{aligned} \mathrm{Var}[\mathbb{E}(S|Y_{1,1})]&=\mathrm{Var}\left[\mathbb{E}\left.\left(1+\sum_{i=1}^{Y_{1,1}} S_i\right|Y_{1,1}\right)\right]\\ &=\mathrm{Var}(1+Y_{1,1}\mathbb{E}S)\\ &=\frac{\mathrm{Var}(Y_{1,1})}{(1-m)^2}=\frac{\sigma^2}{(1-m)^2},\\[10pt] \mathbb{E}[\mathrm{Var}(S|Y_{1,1})]&=\mathbb{E}\left[\mathrm{Var}\left.\left(1+\sum_{i=1}^{Y_{1,1}} S_i\right|Y_{1,1}\right)\right]\\ &=\mathbb{E}\left[\sum_{i=1}^{Y_{1,1}}\mathrm{Var}(S_i)\right]\\ &=\mathbb{E}(Y_{1,1})\mathrm{Var}(S)=m\mathrm{Var}(S). \end{aligned}

综上可得Var(S)=σ2(1m)2+mVar(S)\displaystyle\mathrm{Var}(S)=\frac{\sigma^2}{(1-m)^2}+m\mathrm{Var}(S),解得Var(S)=σ2(1m)3\mathrm{Var}(S)=\dfrac{\sigma^2}{(1-m)^3}

进一步可得到二阶矩E(S2)=σ2+1m(1m)3\mathbb{E}(S^2)=\dfrac{\sigma^2+1-m}{(1-m)^3}

布朗运动

对任意a>0a > 0,参考以下结论,求(B(t),τa)(B(t), \tau_a)的联合密度函数f(x,s)f(x, s),其中xR,0<stx\in\R, 0 < s\leq t

对任意t,s>0t, s > 0(B(t),B(t+s))(B^*(t), B(t+s))的联合概率密度函数g(z,x)=z2zuπttsexp{(2zu)22t(xu)22s}du,z>0,xR.\displaystyle g(z, x) = \int_{-\infty}^z \frac{2z - u}{\pi t \sqrt{ts}}\exp\left\{-\frac{(2z - u)^2}{2t} - \frac{(x - u)^2}{2s}\right\}\mathrm{d}u, \quad z > 0, x \in \mathbb{R}.

点击查看答案

因为{τas}\{\tau_a\leq s\}{B(s)a}\{B^*(s) \geq a\}等价,所以(B(t),τa)(B(t), \tau_a)的分布函数

F(x,s)=axg(z,y)dydz.F(x,s)=\int_a^\infty\int_{-\infty}^x g(z,y)\mathrm{d}y\mathrm{d}z.

又已知(B(t),B(t+s))(B^*(t), B(t+s))的联合概率密度函数g(z,x)g(z,x),令z=az=a

g(a,x)=a2auπttsexp{(2au)22t(xu)22s}du.g(a,x)=\int_{-\infty}^a \frac{2a - u}{\pi t \sqrt{ts}}\exp\left\{-\frac{(2a - u)^2}{2t} - \frac{(x - u)^2}{2s}\right\}\mathrm{d}u.

u=2avu=2a-v

g(a,x)=avπttsexp{v22t(x2a+v)22s}dv=avπttsexp{v22t(x2a)22sv22s+v(x2a)s}dv=1πttsexp{(x2a)22s}avexp{v22tv22s+v(x2a)s}dv=a2πs3/2tsexp(a22s(xa)22(ts)).\begin{aligned} g(a,x)&=\int_a^\infty \frac{v}{\pi t \sqrt{ts}}\exp\left\{-\frac{v^2}{2t} - \frac{(x - 2a + v)^2}{2s}\right\}\mathrm{d}v\\ &=\int_a^\infty \frac{v}{\pi t \sqrt{ts}}\exp\left\{-\frac{v^2}{2t} - \frac{(x - 2a)^2}{2s} - \frac{v^2}{2s} + \frac{v(x - 2a)}{s}\right\}\mathrm{d}v\\ &=\frac{1}{\pi t \sqrt{ts}}\exp\left\{-\frac{(x - 2a)^2}{2s}\right\}\int_a^\infty v \exp\left\{-\frac{v^2}{2t} - \frac{v^2}{2s} + \frac{v(x - 2a)}{s}\right\}\mathrm{d}v\\ &=\frac{a}{2\pi s^{3/2} \sqrt{t - s}} \exp\left(-\frac{a^2}{2s} - \frac{(x - a)^2}{2(t - s)}\right). \end{aligned}

因此(B(t),τa)(B(t), \tau_a)的联合密度函数

f(x,s)=a2πs3/2tsexp(a22s(xa)22(ts)).f(x,s)=\frac{a}{2\pi s^{3/2} \sqrt{t - s}} \exp\left(-\frac{a^2}{2s} - \frac{(x - a)^2}{2(t - s)}\right).

不过本题还有一个相对简单的证明方法(不使用上面的结论):
τa\tau_a的密度函数为g(s)=a2πs3ea22s\displaystyle g(s) = \frac{a}{\sqrt{2\pi s^3}} e^{-\frac{a^2}{2s}}。由布朗运动的马氏性和独立增量性可知,在τa\tau_a给定的条件下,B(t)B(t)的分布

P(B(t)xτa=s)=P(B(t)xB(s)=a,B(u)<a,0<u<s)=P(B(t)xB(s)=a)=P(B(ts)xa).\begin{aligned} P(B(t) \leq x | \tau_a = s) &= P(B(t) \leq x | B(s) = a, B(u) < a, 0 < u < s)\\ &= P(B(t) \leq x | B(s) = a) = P(B(t-s) \leq x - a). \end{aligned}

因此B(t)B(t)的条件密度为fB(t)τa(xs)=12π(ts)e(xa)22(ts)\displaystyle f_{B(t)|\tau_a}(x|s) = \frac{1}{\sqrt{2\pi(t-s)}} e^{-\frac{(x-a)^2}{2(t-s)}},其中t=st=s时退化为1{x=a}1_{\{x=a\}}
(B(t),τa)(B(t), \tau_a)的联合密度函数f(x,s)=g(s)fB(t)τa(xs)=a2πs3(ts)ea22s+(xa)22(ts)\displaystyle f(x, s) = g(s) f_{B(t)|\tau_a}(x|s) = \frac{a}{2\pi \sqrt{s^3(t-s)}} e^{-\frac{a^2}{2s}+\frac{(x-a)^2}{2(t-s)}}

对任意aRa \in \mathbb{R},记τa\tau_a为标准布朗运动B(t)B(t)首次访问aa的时间,以及La(λ)=E(eλτa),λ>0L_a(\lambda) = \mathbb{E}(e^{-\lambda \tau_a}), \lambda > 0
(1)证明:La(λ)=L1(λ)a=L1(a2λ)L_a(\lambda)=L_1(\lambda)^{|a|}=L_1(a^2\lambda)
(2)计算La(λ)L_a(\lambda)

点击查看答案

(1)当a=0a=0时,τ0=0\tau_0=0,因此L0(λ)=1=L1(λ)0=L1(0)L_0(\lambda)=1=L_1(\lambda)^{0}=L_1(0)。下面先证明a>0a>0的情形:
利用布朗运动的独立增量性容易证明对任意正整数nnτn\tau_n可以看成是nn个独立同分布的τ1\tau_1之和,从而Ln(λ)=L1(λ)nL_n(\lambda) = L_1(\lambda)^n;类似,对任意正整数mmτ1\tau_1可以看作是mm个独立同分布的τ1m\tau_{\frac{1}{m}}之和,因此L1(λ)=L1m(λ)mL_1(\lambda) = L_{\frac{1}{m}}(\lambda)^m,从而L1m(λ)=L1(λ)1mL_{\frac{1}{m}}(\lambda) = L_1(\lambda)^{\frac{1}{m}}。因此对任意有理数nm,Lnm(λ)=L1(λ)nm\dfrac{n}{m},L_{\frac{n}{m}}(\lambda) = L_1(\lambda)^{\frac{n}{m}}.
由布朗运动轨道连续性可知τa\tau_aaa增加单调不降,因此La(λ)L_a(\lambda)aa单调下降。对任意无理数aa,存在有理数列bn,cnb_n, c_n使得bn<a<cnb_n < a < c_nbna,cnab_n \to a,\, c_n \to a。从而

L1(λ)a=limnLbn(λ)La(λ)limnLcn(λ)=L1(λ)a.L_1(\lambda)^a = \lim_{n \to \infty} L_{b_n}(\lambda) \ge L_a(\lambda) \ge \lim_{n \to \infty} L_{c_n}(\lambda) = L_1(\lambda)^a.

由此可得La(λ)=L1(λ)aL_a(\lambda) = L_1(\lambda)^a对任意a>0a > 0成立。
又因为1aB(a2t)\dfrac{1}{a}B(a^2t)也是标准布朗运动(定义为WtW_t),则其首次访问11的时间τ1W=τaa2\tau_1^W=\dfrac{\tau_a}{a^2}(与τ1\tau_1同分布)。所以

La(λ)=E(eλτa)=E(eλa2τ1)=L1(a2λ).L_a(\lambda)=\mathbb{E}(e^{-\lambda \tau_a})=\mathbb{E}(e^{-\lambda a^2\tau_1})=L_1(a^2\lambda).

对于a<0a< 0的情形,因为τa\tau_{a}τa\tau_{-a}同分布,所以La(λ)=La(λ)L_a(\lambda)=L_{-a}(\lambda)
综上,La(λ)=L1(λ)a=L1(a2λ)L_a(\lambda)=L_1(\lambda)^{|a|}=L_1(a^2\lambda)
(2)因为τa\tau_a的密度函数为f(t)=a2πt3ea22tf(t) = \frac{|a|}{\sqrt{2\pi t^3}} e^{-\frac{a^2}{2t}},因此

La(λ)=E(eλτa)=0eλtf(t)dt=0a2πt3ea2+2λt22tdt=e2λa0a2πt3e12(at2λt)2dt=e2λa022πe12(u2λau)2du(u=at)=e2λa+12πe12w2(1+ww2+42λa)dw(w=u2λau)=e2λa.\begin{aligned} L_a(\lambda)=\mathbb{E}(e^{-\lambda \tau_a}) &= \int_{0}^{\infty} e^{-\lambda t} f(t) \mathrm{d}t = \int_{0}^{\infty} \frac{|a|}{\sqrt{2\pi t^3}} e^{-\frac{a^2 + 2\lambda t^2}{2t}} \mathrm{d}t \\ &= e^{-\sqrt{2\lambda}|a|} \int_{0}^{\infty} \frac{|a|}{\sqrt{2\pi t^3}} e^{-\frac{1}{2} \left( \frac{|a|}{\sqrt{t}} - \sqrt{2\lambda t} \right)^2} \mathrm{d}t \\ &= e^{-\sqrt{2\lambda}|a|} \int_{0}^{\infty} \frac{2}{\sqrt{2\pi}} e^{-\frac{1}{2} \left( u - \frac{\sqrt{2\lambda}|a|}{u} \right)^2} du \quad \left( u = \frac{|a|}{\sqrt{t}} \right) \\ &= e^{-\sqrt{2\lambda}|a|} \int_{-\infty}^{+\infty} \frac{1}{\sqrt{2\pi}} e^{-\frac{1}{2} w^2} \left( 1 + \frac{w}{\sqrt{w^2 + 4\sqrt{2\lambda}|a|}} \right) \mathrm{d}w \quad \left( w = u - \frac{\sqrt{2\lambda}|a|}{u} \right) \\ &= e^{-\sqrt{2\lambda}|a|}. \end{aligned}

当然也可以用鞅的停时定理计算:可以验证{Mn}={eλBnnλ22,n0}\{M_n\}=\{e^{\lambda B_n-\frac{n\lambda^2}{2}},n\geq 0\}为鞅,且代入τa\tau_a后满足停时定理条件,故

1=E[M0]=E[Mτa]=E[eλaτaλ22]=eλaE[eτaλ22]1=\mathbb{E}[M_0]=\mathbb{E}[M_{\tau_a}]=\mathbb{E}[e^{\lambda a-\frac{\tau_a\lambda^2}{2}}]=e^{\lambda a}\mathbb{E}[e^{-\frac{\tau_a\lambda^2}{2}}]

因此

La(λ22)=eλaLa(λ)=ea2λ.L_a\left(\frac{\lambda^2}{2}\right)=e^{-\lambda a}\Longrightarrow L_a(\lambda)=e^{|a|\sqrt{2\lambda}}.

M={M(t);t0}M = \{M(t); t \geq 0\}是指数布朗运动,即M(t)=eμt+σB(t)M(t) = e^{\mu t + \sigma B(t)},求E(M(s)M(t)=x)\mathbb{E}(M(s)|M(t) = x)

点击查看答案

因为E(M(s)M(t)=x)=eμsE(eσB(s)B(t)=lnxμσ)\displaystyle \mathbb{E}(M(s)|M(t) = x) = e^{\mu s} \mathbb{E}\left(e^{\sigma B(s)} | B(t) = \frac{\ln x - \mu}{\sigma}\right),下面进行分类讨论:
sts \leq t时,B(s)N(s(lnxμt)σt,(ts)st)\displaystyle B(s) \sim N \left( \frac{s(\ln x - \mu t)}{\sigma t}, \frac{(t-s)s}{t} \right),因此

E(eσB(s)B(t)=(lnxμ)/σ)=es(lnxμt)t+σ2(ts)s2t.\mathbb{E}(e^{\sigma B(s)} | B(t) = (\ln x - \mu)/\sigma) = e^{\frac{s(\ln x - \mu t)}{t} + \frac{\sigma^2 (t-s)s}{2t}}.

进而E(M(s)M(t)=x)=eslnxt+σ2(ts)s2t\displaystyle \mathbb{E}(M(s)|M(t) = x) = e^{\frac{s \ln x}{t} + \frac{\sigma^2 (t-s)s}{2t}}
s>ts > t时,B(s)N(lnxμtσ,st)\displaystyle B(s) \sim N \left( \frac{\ln x - \mu t}{\sigma}, s - t \right),因此

E(eσB(s)B(t)=(lnxμt)/σ)=elnxμt+σ2(st)2.\mathbb{E}(e^{\sigma B(s)} | B(t) = (\ln x - \mu t)/\sigma) = e^{\ln x - \mu t + \frac{\sigma^2 (s-t)}{2}}.

进而E(M(s)M(t)=x)=elnx+μ(st)+σ2(st)2\displaystyle \mathbb{E}(M(s)|M(t) = x) = e^{\ln x + \mu (s-t) + \frac{\sigma^2 (s-t)}{2}}

τ2\tau_2是标准布朗运动{B(t):t0}\{B(t): t \geq 0\}首次到达22的时间,求aa使得τ2\tau_2a[B(2)]2\dfrac{a}{[B(2)]^2}同分布。$
点击查看答案

方法一:由于T2T_2非负,aa非负。对任意λ0\lambda \geq 0,

E(eλT2)=e22λ,E(eλaB2(2))=2014πeλax2ex24dx=2014πe14(x22ax)2dxe2a=e2a=P(B(t)2)=P(t2B(2)2).\begin{aligned} \mathbb{E}(e^{-\lambda T_2}) &= e^{-2\sqrt{2\lambda}},\\ \mathbb{E}(e^{-\lambda \frac{a}{B^2(2)}}) &= 2 \int_0^\infty \frac{1}{\sqrt{4\pi}} e^{-\lambda \frac{a}{x^2}} e^{-\frac{x^2}{4}} \mathrm{d}x\\ &= 2 \int_0^\infty \frac{1}{\sqrt{4\pi}} e^{-\frac{1}{4}\left(x - \frac{2\sqrt{2a}}{x}\right)^2} \mathrm{d}x e^{-\sqrt{2a}} = e^{-\sqrt{2a}}\\ &= P(|B(t)| \geq 2) = P(\sqrt{\frac{t}{2}}|B(2)| \geq 2). \end{aligned}

令两式相等,可得a=8a = 8
方法二:对任意t0t \geq 0

P(T2t)=P(B(t)2)=P(B(t)2)=P(t2B(t)2)=P(t2B(2)24)=P(8B(2)2t).\begin{aligned} P(T_2 \leq t) &= P(B^*(t) \geq 2)\\ &=P(|B(t)|\geq 2)=P\left(\sqrt{\frac{t}{2}}|B(t)|\geq 2\right)\\ &=P\left(\frac{t}{2}|B(2)|^2 \geq 4\right) = P\left(\frac{8}{|B(2)|^2} \leq t\right). \end{aligned}

由此可知a=8a = 8

离散时间鞅

XX是鞅,则对任意nmn\geq mE(XnXm)=E(Xm2)\mathbb{E}(X_nX_m)=\mathbb{E}(X_m^2)

推导如下:

E(XnXm)=E(E(XnXmX0,,Xn1))=E(XmE(XnX0,,Xn1))=E(XmXn1)==E(Xm2).\begin{aligned} \mathbb{E}(X_nX_m)&=\mathbb{E}(\mathbb{E}(X_nX_m|X_0,\cdots,X_{n-1}))\\ &=\mathbb{E}(X_m\mathbb{E}(X_n|X_0,\cdots,X_{n-1}))\\ &=\mathbb{E}(X_mX_{n-1})\\ &=\cdots=\mathbb{E}(X_m^2). \end{aligned}

Y={Yn;n0}Y=\{Y_n;n\geq 0\}关于X={Xn;n0}X=\{X_n;n\geq 0\}是鞅,TT为一停时,证明:
(1){YnT;n0}\{Y_{n\land T};n\geq 0\}关于XX也是鞅;
(2)E(YnT)=E(Y0)\mathbb{E}(Y_{n\land T})=\mathbb{E}(Y_0),且E(YnT)\mathbb{E}(|Y_{n\land T}|)单调不降。
(3)当n=1E(YnYn11{Tn})<+\displaystyle\sum_{n=1}^\infty \mathbb{E}\left( |Y_n - Y_{n-1}| \mathbf{1}_{\{T \ge n\}} \right) < +\infty时,YnT=E(YTX0,,Xn)Y_{n\land T}=\mathbb{E}(Y_T|X_0,\cdots,X_n)

点击查看答案

(1)因为E(YnT)i=1nE(Yi)<\displaystyle \mathbb{E}(|Y_{n\land T}|)\leq\sum_{i=1}^n\mathbb{E}(|Y_i|)< \infty,故YnTY_{n\land T}可积;
又因为

E(Y(n+1)TX1,,Xn)=E(Y(n+1)T1{Tn+1}X1,,Xn)+E(Y(n+1)T1{Tn}X1,,Xn)=1{Tn+1}E(Yn+1X1,,Xn)+k=1nE(Yk1{T=k}X1,,Xn)=1{Tn+1}Yn+k=1nYk1{T=k}=YnT\begin{aligned} &\mathbb{E}(|Y_{(n+1)\land T}|\mid X_1,\cdots,X_n)\\ =&\mathbb{E}(|Y_{(n+1)\land T}|\mathbf{1}_{\{T\geq n+1\}}\mid X_1,\cdots,X_n)+\mathbb{E}(|Y_{(n+1)\land T}|\mathbf{1}_{\{T\leq n\}}\mid X_1,\cdots,X_n)\\ =&\mathbf{1}_{\{T\geq n+1\}}\mathbb{E}(|Y_{n+1}\mid X_1,\cdots,X_n)+\sum_{k=1}^n\mathbb{E}(|Y_{k}|\mathbf{1}_{\{T=k\}}\mid X_1,\cdots,X_n)\\ =&\mathbf{1}_{\{T\geq n+1\}}Y_n+\sum_{k=1}^n|Y_{k}|\mathbf{1}_{\{T=k\}}\\ =&Y_{n\land T} \end{aligned}

因此{YnT;n0}\{Y_{n\land T};n\geq 0\}关于XX也是鞅。
(2)因为nTn\land T是有界停时,故满足有界停时定理,因此E(YnT)=E(Y0)\mathbb{E}(Y_{n\land T})=\mathbb{E}(Y_0)

注:E(YT)=E(Y0)\mathbb{E}(Y_{T})=\mathbb{E}(Y_0)不一定成立(需满足一般停时定理条件)

又因为{YnT;n0}\{Y_{n\land T};n\geq 0\}关于XX是鞅,所以{YnT;n0}\{|Y_{n\land T}|;n\geq 0\}关于XX是下鞅,故E(YnT)\mathbb{E}(|Y_{n\land T}|)单调不降。
(3)证明:记D0=M0D_0 = M_0,Dm=MmMm1D_m = M_m - M_{m-1},m1m \geq 1。那么

E(YTX0,,Xn)=E(k=0nYk1{T=k}X0,,Xn)=E(k=0nm=0kDm1{T=k}X0,,Xn)=E(m=0nDm1{Tm}X0,,Xn)=m=0nDm1{Tm}+E(m=n+1Dm1{Tm}X0,,Xn).\begin{aligned} \mathbb{E}(Y_T | X_0, \cdots, X_n) &= \mathbb{E} \left( \sum_{k=0}^n Y_k 1_{\{T=k\}} \mid X_0, \cdots, X_n \right)\\ &= \mathbb{E} \left( \sum_{k=0}^n \sum_{m=0}^k D_m 1_{\{T=k\}} \mid X_0, \cdots, X_n \right)\\ &= \mathbb{E} \left( \sum_{m=0}^n D_m 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right)\\ &= \sum_{m=0}^n D_m 1_{\{T \geq m\}} + \mathbb{E} \left( \sum_{m=n+1}^\infty D_m 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right). \end{aligned}

E(Dm+1X0,,Xm)=0\mathbb{E}(D_{m+1} | X_0, \cdots, X_m) = 0及重期望公式可知对任意N>nN > n,

E(m=n+1Dm1{Tm}X0,,Xn)=E(m=N+1Dm1{Tm}X0,,Xn).\mathbb{E} \left( \sum_{m=n+1}^\infty D_m 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right) = \mathbb{E} \left( \sum_{m=N+1}^\infty D_m 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right).

注意到

E(m=N+1Dm1{Tm}X0,,Xn)E(m=N+1Dm1{Tm}X0,,Xn)\left| \mathbb{E} \left( \sum_{m=N+1}^\infty D_m 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right) \right| \leq \mathbb{E} \left( \sum_{m=N+1}^\infty |D_m| 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right)

上式右端非负,且随NN的增加单调不增,而且由条件可知

E(E(m=N+1Dm1{Tm}X0,,Xn))=E(m=N+1Dm1{Tm})0\mathbb{E} \left( \mathbb{E} \left( \sum_{m=N+1}^\infty |D_m| 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right) \right)= \mathbb{E} \left( \sum_{m=N+1}^\infty |D_m| 1_{\{T \geq m\}} \right) \to 0

这表明

E(m=N+1Dm1{Tm}X0,,Xn)0a.s.\mathbb{E} \left( \sum_{m=N+1}^\infty |D_m| 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right) \to 0 \quad \text{a.s.}

NN的任意性可知

E(m=n+1Dm1{Tm}X0,,Xn)=0\mathbb{E} \left( \sum_{m=n+1}^\infty D_m 1_{\{T \geq m\}} \mid X_0, \cdots, X_n \right) = 0

因此

E(YTX0,,Xn)=m=0nDm1{Tm}=m=0nTDm=XnT.\mathbb{E}(Y_T | X_0, \cdots, X_n) = \sum_{m=0}^n D_m 1_{\{T \geq m\}} = \sum_{m=0}^{n\land T} D_m = X_{n\land T}.

已知Yk,k1Y_k, k \geq 1独立且都服从参数为θ\theta的泊松分布。令Sn=k=1n(1Yk),n1\displaystyle S_n = \sum_{k=1}^n (1 - Y_k), n \geq 1
(1)θ>1\theta > 1时求P(maxn1Sn1)\displaystyle P\left(\max_{n \geq 1} S_n \geq 1\right)
(2)θ1\theta \leq 1时,令τ=inf{n1;Sn1}\tau = \inf\{n \geq 1; S_n \geq 1\},对任意u>0u > 0,求E(euτ)\mathbb{E}(e^{-u\tau})

点击查看答案

τ=inf{n1;Sn1}\tau = \inf\{n \geq 1; S_n \geq 1\}。由于SnS_n是右连续的随机游动,因此若τ<\tau < \infty,则Sτ=1S_\tau = 1
注意到,对任意λ>0\lambda > 0E(eλY1)=eλθ+eλθ\mathbb{E}(e^{\lambda Y_1}) = e^{\lambda - \theta + e^{-\lambda \theta}},记Φ(θ,λ)=λθ+eλθ\Phi(\theta, \lambda) = \lambda - \theta + e^{-\lambda \theta}
(1)即求P(τ<+)P(\tau < +\infty)。当θ>1\theta > 1时,存在λ0>0\lambda_0 > 0使得Φ(θ,λ0)=0\Phi(\theta, \lambda_0) = 0,且对任意λ>λ0\lambda > \lambda_0Φ(θ,λ)>0\Phi(\theta, \lambda) > 0。容易验证eλSnΦ(θ,λ)ne^{\lambda S_n - \Phi(\theta, \lambda) n}关于SS是鞅。由有界停时定理

1=E(eλS0)=E(eλSnτΦ(θ,λ)(nτ))=E(eλΦ(θ,λ)τ1{τ<n})+E(eλSnΦ(θ,λ)n1{τ>n})\begin{aligned} 1 = \mathbb{E}(e^{\lambda S_0}) &= \mathbb{E}(e^{\lambda S_n \wedge \tau - \Phi(\theta, \lambda)(n \wedge \tau)})\\ &= \mathbb{E}(e^{\lambda - \Phi(\theta, \lambda) \tau 1_{\{\tau < n\}}}) + \mathbb{E}(e^{\lambda S_n - \Phi(\theta, \lambda) n 1_{\{\tau > n\}}}) \end{aligned}

注意到τ>n\tau > n时,Sn<1S_n < 1。令nn \to \infty1=E(eλΦ(θ,λ)τ1{τ<+})1 = \mathbb{E}(e^{\lambda - \Phi(\theta, \lambda) \tau 1_{\{\tau < +\infty\}}})
再令λλ0\lambda \to \lambda_0,得1=E(eλ01{τ<+})1 = \mathbb{E}(e^{\lambda_0 1_{\{\tau < +\infty\}}})。由此可得P(τ<+)=eλ0P(\tau < +\infty) = e^{-\lambda_0}
(2)当θ1\theta \leq 1时,重复上面讨论。注意到Φ(θ,0)=0\Phi(\theta, 0) = 0,而且对任意λ>0\lambda > 0Φ(θ,λ)>0\Phi(\theta, \lambda) > 0。即此时上面的λ0=0\lambda_0 = 0,因此P(τ<+)=1P(\tau < +\infty) = 1。由此可得eλ=E(eΦ(θ,λ)τ)e^{-\lambda} = \mathbb{E}(e^{-\Phi(\theta, \lambda) \tau})。取λu>0\lambda_u > 0使得

Φ(θ,λu)=λuθ+eλuθ=u,\Phi(\theta, \lambda_u) = \lambda_u - \theta + e^{-\lambda_u \theta} = u,

那么E(euτ)=eλu\mathbb{E}(e^{-u \tau}) = e^{-\lambda_u}