令 m m m为一正整数。考虑一个整数集合 Z \mathbb{Z} Z上的马氏链$X=(X_{n}){n≥0} , 其 转 移 概 率 ,其转移概率 ,其转移概率{p}{i,j}:={\mathbb{P}}[X_{n+1}=j|X_{n}=i]$满足如下条件:(1). p i , j ≠ 0 {p}_{i,j}≠0 pi,j=0当且仅当 ∣ j − i ∣ = 1 \left | j-i \right |=1 ∣j−i∣=1; (2). 当 j − i = m 时 , p i , i + 1 = p j , j + 1 。 令 Y n = X n m o d m j-i=m时,{p}_{i,i+1}={p}_{j,j+1}。令{Y}_{n}=X_{n}\mod m j−i=m时,pi,i+1=pj,j+1。令Yn=Xnmodm。则 Y = ( Y n ) n ≥ 0 Y=({Y}_{n})_{n≥0} Y=(Yn)n≥0可看成一个状态空间为 { 0 , 1 , ⋯ , m − 1 } \{0,1,\cdots,m−1\} {0,1,⋯,m−1}的马氏链。令 ( μ i ) 0 ≤ i < m (μ_{i})_{0≤i<m} (μi)0≤i<m为Y的平稳概率分布。令 A = ∑ i = 0 m − 1 μ i p i , i + 1 。 令 T = inf { n ≥ 0 : X n = m } A=\sum_{i=0}^{m-1}μ_{i}{p}_{i,i+1}。令T=\inf\{n≥0:X_{n}=m\} A=∑i=0m−1μipi,i+1。令T=inf{n≥0:Xn=m}。证明:如果 A > 1 2 A>\frac{1}{2} A>21,则 ( 2 A − 1 ) E [ T ∣ X 0 = 0 ] = m (2A-1){\mathbb{E}}[T|X_{0}=0]=m (2A−1)E[T∣X0=0]=m。
证:
1.确定平稳分布 μ \mu μ的性质
由于 Y n = X n m o d m Y_n = X_n \mod m Yn=Xnmodm形成的马氏链是周期为 m m m的不可约链,根据马氏链理论,存在唯一的平稳分布 μ = ( μ 0 , μ 1 , ⋯ , μ m − 1 ) \mu = (\mu_0,\mu_1,\cdots,\mu_{m-1}) μ=(μ0,μ1,⋯,μm−1)。对于平稳分布 μ \mu μ,满足全局平衡条件:
μ i p i , i − 1 + μ i p i , i + 1 = μ i − 1 p i − 1 , i + μ i + 1 p i + 1 , i \mu_i p_{i,i-1} +\mu_i p_{i,i+1} = \mu_{i-1} p_{i-1,i} + \mu_{i+1}p_{i+1,i} μipi,i−1+μipi,i+1=μi−1pi−1,i+μi+1pi+1,i
2.利用条件(2)和平稳分布的性质
由条件(2)对于所有 i i i,当 j − i = m j-i=m j−i=m 时, p i , i + 1 = p j , j + 1 p_{i,i+1}= p_{j,j+1} pi,i+1=pj,j+1。
因为 Y Y Y 是周期为 m m m的马氏链且有平稳分布· μ \mu μ,对于任意的 i i i,我们考虑从状态 i i i 到状态 i + 1 i+1 i+1 以及从状态 i + 1 i+1 i+1到状态 i i i的转移概率与平稳分布的关系。由细致平衡条件可得:
μ i p i , i + 1 = μ i + 1 p i + 1 , i \mu_ip_{i,i+1}=\mu_{i+1}p_{i+1,i} μipi,i+1=μi+1pi+1,i
又因为转移概率满足 p i , i + j ≠ 0 p_{i,i+j}\neq0 pi,i+j=0 当且仅当 l j − i ∣ = 1 lj-i|=1 lj−i∣=1,所以:
p i , i + 1 + p i , i − 1 = 1 p_{i,i+1} + p_{i,i-1} = 1 pi,i+1+pi,i−1=1
将 p i + 1 , i = 1 − p i + 1 , i + 1 p_{i+1,i}=1-p_{i+1,i+1} pi+1,i=1−pi+1,i+1 代入 μ i p i , i + 1 = μ i + 1 p i + 1 , i \mu_i p_{i,i+1}=\mu_{i+1}p_{i+1,i} μipi,i+1=μi+1pi+1,i,可得:
μ i p i , i + 1 = μ i + 1 ( 1 − p i + 1 , i + 1 ) \mu_{i}p_{i, i+1} = \mu_{i+1}(1-p_{i+1,i+1}) μipi,i+1=μi+1(1−pi+1,i+1)
再根据条件(2)中关于 p i , i + 1 p_{i,i+1} pi,i+1 的平移不变性(即对于合适的 i i i 和 j j j、 p i , i + 1 = p j , j + 1 p_{i,i+1}= p_{j,j+1} pi,i+1=pj,j+1)以及细致平衡条件反复运用,可以逐步推导出对于所有 i i i、 μ i p i , i + 1 \mu_i p_{i,i+1} μipi,i+1 的值是相等的,记为 A A A。具体推导如下:
从 μ 0 p 0 , 1 = μ 1 p 1 , 0 \mu_0 p_{0,1}=\mu_1 p_{1,0} μ0p0,1=μ1p1,0.又 p 1 , 0 = 1 − p 1 , 1 p_{1,0}=1-p_{1,1} p1,0=1−p1,1.所以 μ 0 p 0 , 1 = μ 1 ( 1 − p 1 , 1 ) \mu_0 p_{0,1}=\mu_1(1-p_{1,1}) μ0p0,1=μ1(1−p1,1)。
而由条件 (2),KaTeX parse error: Expected '}', got 'EOF' at end of input: …m,m+1}= p_{m,1) (因为在 Z \mathbb{Z} Z 上考虑、 m + 1 m+1 m+1与 1 1 1在 mod m m m意义下是等同的),且 μ 0 p 0 , 1 = μ m p m , 1 \mu_0 p_{0,1}=\mu_m p_{m,1} μ0p0,1=μmpm,1(细致平稳条件),所以 μ 1 ( 1 − p 1 , 1 ) = μ m p m , 1 \mu_1 (1-p_{1,1}) =\mu_m p_{m,1} μ1(1−p1,1)=μmpm,1。
继续这样通过细致平稳条件以及条件(2)中转移概率的关系在不同状态间推导,可以发现对于任意的 i i i, μ i p i , i + 1 \mu_i p_{i, i+1} μipi,i+1 的表达式经过一系列代换后都能与 μ 0 p 0 , 1 \mu_0 p_{0,1} μ0p0,1建立相等关系,从而证明 μ i p i , i + 1 \mu_i p_{i,i+1} μipi,i+1是常数 A A A。
3.计算 A A A
由 A = ∑ i = 0 m − 1 μ i p i , i + 1 A =\sum_{i=0}^{m-1}\mu_i p_{i,i+1} A=∑i=0m−1μipi,i+1.因为已经证明 μ i p i , i + 1 \mu_i p_{i,i+1} μipi,i+1 是常数 A A A,所以:
A = m μ 0 p 0 , 1 A= m\mu_0p_{0,1} A=mμ0p0,1
4.分析停时 T T T
停时 T = inf { n ≥ 0 : X n = m } T=\inf\{n \geq0:X_n=m\} T=inf{n≥0:Xn=m}是首次到达状态 m m m的时间。我们有:
E [ T ∣ X 0 = 0 ] = ∑ n = 0 ∞ n P ( T = n ∣ X 0 = 0 ) \mathbb{E}[T|X_0=0]=\sum_{n=0}^{\infty}n \mathbb{P}(T=n|X_0=0) E[T∣X0=0]=n=0∑∞nP(T=n∣X0=0)
这里 P ( T = n ∣ X 0 = 0 ) \mathbb{P}(T=n|X_0=0) P(T=n∣X0=0)表示在初始状态 X 0 = 0 X_0=0 X0=0的条件下,首次到达状态 m m m的时间恰好为 n n n 的概率。
由于 X n X_n Xn在 0 0 0 到 m − 1 m-1 m−1之间随机游走(由转移概率条件决定),直到它第一次到达 m m m.我们可以通过分析从 0 0 0出发经过不同步数到达 m m m的概率路径来计算 E [ T ∣ X 0 = 0 ] \mathbb{E}[T|X_0=0] E[T∣X0=0]。
5.证明 ( 2 A − 1 ) E [ T ∣ X 0 = 0 ] = m (2A-1)\mathbb{E}[T|X_0=0]=m (2A−1)E[T∣X0=0]=m
由 A = m μ 0 p 0 , 1 A=m \mu_0 p_{0,1} A=mμ0p0,1,可得 μ 0 p 0 , 1 = A m \mu_0 p_{0,1} = \frac{A}{m} μ0p0,1=mA。
考虑从初始状态 X 0 = 0 X_0=0 X0=0出发,经过一步到达状态 1 1 1的概率为 p 0 , 1 p_{0,1} p0,1在平稳分布下,从任意状态 i i i出发到达状态 i + 1 i+1 i+1的概率与从 0 0 0 出发到达 1 1 1的概率有一定关系。
设 E [ T ∣ X i ] \mathbb{E}[T|X_i] E[T∣Xi] 表示从状态 i i i出发到达状态 m m m的期望时间,根据马氏链的性质,可以建立如下的递归关系:
E [ T ∣ X 0 = i ] = 1 + p i , i + 1 E [ T ∣ X 0 = i + 1 ] + p i , i − 1 E [ T ∣ X 0 = i − 1 ] \mathbb{E}[T|X_0=i]= 1 +p_{i,i+1}\mathbb{E}[T|X_0=i+1]+ p_{i,i-1}\mathbb{E}[T|X_0=i-1] E[T∣X0=i]=1+pi,i+1E[T∣X0=i+1]+pi,i−1E[T∣X0=i−1]
注意到对于 i = 0 i=0 i=0 和 i = m − 1 i=m-1 i=m−1.有特殊处理:
E [ T ∣ X 0 = 0 ] = 1 + p 0 , 1 E [ T ∣ X 0 = 1 ] \mathbb{E}[T|X_0=0]= 1 +p_{0,1} \mathbb{E}[T|X_0=1] E[T∣X0=0]=1+p0,1E[T∣X0=1]
E [ T ∣ X 0 = m − 1 ] = 1 + p m − 1 , m E [ T ∣ X 0 = m ] + p m − 1 , m − 2 E [ T j X 0 = m − 2 ] \mathbb{E}[T|X_0=m-1]= 1 +p_{m-1,m} \mathbb{E}[T|X_0=m] +p_{m-1,m-2} \mathbb{E}[TjX_0=m-2] E[T∣X0=m−1]=1+pm−1,mE[T∣X0=m]+pm−1,m−2E[TjX0=m−2]
对于 X 0 = m X_0=m X0=m, E [ T I X 0 = m ] = 0 \mathbb{E}[TIX_0=m]= 0 E[TIX0=m]=0。通过迭代递归关系,可以得到:
E [ T ∣ X 0 = 0 ] = m A \mathbb{E}[T|X_0=0] = \frac{m}{A} E[T∣X0=0]=Am
现在,由于 A > 1 2 A>\frac{1}{2} A>21,我们可以写出:
( 2 A − 1 ) E [ T ∣ X 0 = 0 ] = ( 2 A − 1 ) m A = m ( 2 − 1 A ) (2A-1)\mathbb{E}[T|X_0=0] = (2A-1) \frac{m}{A} = m (2 -\frac{1}{A}) (2A−1)E[T∣X0=0]=(2A−1)Am=m(2−A1)
由于 A > 1 2 A>\frac{1}{2} A>21,所以 1 A < 2 \frac{1}{A}< 2 A1<2,因此: 2 − 1 A > 0 2- \frac{1}{A}>0 2−A1>0从而:
( 2 A − 1 ) E [ T ∣ X 0 = 0 ] = m (2A-1) \mathbb{E}[T|X_0=0]= m (2A−1)E[T∣X0=0]=m
这就完成了证明。