Loading [MathJax]/jax/output/HTML-CSS/jax.js
  数学物理学报  2015, Vol. 35 Issue (6): 1180-1189   PDF (326 KB)    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
孙鸿雁
么亚楠
带停留对称随机游动局部时的渐进行为
孙鸿雁1, 么亚楠2    
1 中国地质大学(北京) 北京 100081;
2 北京师范大学 北京 100875
摘要: 该文研究带停留对称随机游动,模型的局部时可写成一个非退化两物种分支过程的泛函.由这个表示建立模型局部时的收敛性质.
关键词: 布朗局部时     随机游动     分支过程    
Asymptotic Behavior of the Local Time for the Symmetric Random Walk with Stay
Sun Hongyan1, Yao Yanan2     
1 China University of Geosciences, Beijing 100081;
1 School of Mathematical Sciences, Beijing Normal University, Beijing 100875
Abstract: We consider the symmetric random walk with stay in this paper, whose local time can be represented by a function of a non-primitive two-type branching process. Using this representation, we establish the convergence of the local time of the random walk.
Key words: Brownian local time     Random walk     Branching processes    
1引言

Donsker 建立了随机游动的不变原理,即其标准化后的极限为布朗运动. 具体地, 设{ξn:n0} 是直线上一列独立同分布的随机变量, 满足对于任意的n,ξn的期望 Eξn=0及方差 ξn=σ2,其中 σ>0. 定义随机游动{Sn:n0}Sn:=ni=1ξi,n1,S0=0, 及其线性内插S(t): t 0为 S(t):=S[t]+(t[t])(S[t]+1S[t]),t0. 记相应的尺度过程为 Sn(t)=S(nt)n, t0, Donsker 证明了过程 {Sn(t):t0} 弱收敛到布朗运动 {W(t):t0}, 且方差 W(t)=σ2t.

受Donsker 不变原理启发,关于随机游动局部时的不变原理引起了人们的广泛关注. 记随机游动的局部时过程为 phi(x,n),即 phi(x,n){Sn}n 步中到过 x 处的次数. 对于紧邻随机游动,即 ξn 分别以 1/2的概率取 ±1, Borodin[2] 借助于不变嵌入方法建立了局部时的不变原理,即: 两指标的过程 n1/2phi([xn],[nt]) 收敛到布朗运动 {W(t):t0} 在时刻 t 之前于 x 处的局部时. 本文中,[x] 表示 x 向下取整. Perkins[15] 利用组合方法也证明了这一不变原理. 对于一般随机游动,其局部时的不变原理由 Borodin 和 Ibragimov[3] 所建立. Csörgö 和Revész[5] 则基于文献[4] 中随机游动局部时的表示证明了局部时强收敛意义下的不变原理. 进一步, Perkins[16]借助于Tanaka 关于布朗运动局部时的刻画(参见文献 [14]) 获得了一般随机游动与其局部时构成的联合过程的不变原理.

另一方面,Rogers[18]考察了与停时相关的局部时过程的极限行为,即过程 phi([xN],τN)/N的极限行为,其中 τN 为游动第 N 次到达 0 的时间. 事实上,对于紧邻随机游动,Rogers[18] 借助于随机游动的内蕴分支机制证明 了与停时相关的局部时过程弱收敛于布朗运动的与停时相关的局部时. 关于紧邻随机游动的内蕴分支机制参见文献Dwass[6]及Kesten 等人[8]. 但是对于一般随机游动,目前还没有关于与停时相关的局部时过程的极限行为的结论.

本文的主要目的是: 对于带停留的紧邻随机游动,考察与停时相关的局部时过程的极限行为. 称 {Sn} 为带停留的紧邻随机游动是指

P(ξn=1)=p=P(ξn=1),P(ξn=0)=r=12p,0<p<1/2. (1.1)
注意到,当 r=0 时,{Sn} 为紧邻随机游动,不允许粒子在原地停留; 此时与停时相关的局部时过程的极限行为已被 Rogers[18] 所建立. 因此我们感兴趣的是 0<r<1.

定义反射过程 |S| 的局部时为 L(j,n):={r:1rn,|Sr|=j},j0,n1.τ00,对于任意 m1,τm:=inf{n>τm1:Sn=0,Sn10}.D(E) 表示空间E上左极右连函数全体并赋予Skorohod 拓扑. 对任意自然数 N1,通过局部时做尺度变换,我们定义空间 D[0,) 上的过程

lN(x)={L([Nx],τN)N,Nx1,2(1+r1r),Nx<1. (1.2)
另外,记号 lW(x;t) 表示布朗运动 {W(t):t0} 在时刻 t 之前 x 点处的局部时,即: 对于任意Borel可测函数 f, t0f(W(s))ds=f(x)lW(x,t)dx.

本文目标是证明当 N 时,lN() 在空间 D[0,) 上收敛. 即

定理1.1 作为空间 D[0,) 上的一个元素,当 N 时, lN()lW(,τ), 其中 lW(x,t)=lW(x,t)+lW(x,t),τ=inf{t>0:lW(0,t)>σ2}. ``" 表``弱收敛".

注 1.1r=0 时,σ2=1,对应的正是Rogers[1]的结果.

Kesten 等人[8]和Wang[19]揭示了带停留的紧邻随机游动的内蕴分支机制; 该机制对应了一个两物种分支过程,其均值矩阵是非正则的. 注意到, 与停时相关的局部时可表示成这一两物种分支过程的泛函. 本文将借助于这一表达来证明停时相关的局部时的收敛性. 关于多物种分支过程的收敛性,可参见文献[7, 10], 但文献中都是假定均值矩阵是正则的. 本文要处理的是非正则的情形.

定理 1.1 的证明将在下一节展开.

2 定理1.1 的证明

首先,我们来回顾带停留随机游动的内蕴分支机制,具体可参见文献[]. 设 (UN(0),ZN(0))=(N,0). 对于任意 i1, 定义

UN(i):={0nτN:|Sn|=i,|Sn+1|=i+1},ZN(i):={0nτN:|Sn|=i,|Sn+1|=i} (2.1)
分别来表示反射过程 |S| 在时刻 τN 之前,从 i 处到 i+1 以及从 ii 的步数.

Zeitouni[20] 指出过程 |S| 的每个从 ii+1 或者到 i 的步子必然出现在某个连续的从 i1i 的步子之间. 并且在每个连续的从 i1i 的步子之间的这两种步子服从多项分布. 因此,在(UN(i1),ZN(i1)),(UN(i2),ZN(i2)),,(UN(0),ZN(0)) 生成的σ-代数条件下,随机向量 (UN(i),ZN(i)) 同分布于 UN(i1) 个独立同分布的随机变量之和, 这些随机变量服从如下分布 P(V1=(a,b))=(a+ba)parbp, a,b0. 也就是说,(UN(n),ZN(n))n0 是一个两物种分支过程,初始分布为 (UN(0),ZN(0))=(N,0). 下面我们精确刻画这一结论.

引理2.1 (i) (2.1) 式定义的过程 (UN(n),ZN(n))n0 是一个两物种分支过程. 其后代分布为

P((UN(i),ZN(i))=(a,b)|(UN(i1),ZN(i1))=(1,0))=(a+ba)parbp, a,b0,P((UN(i),ZN(i))=(0,0)|(UN(i1),ZN(i1))=(0,1))=1. (2.2)

(ii) 过程 (UN(n))n0 是一个单物种分支过程,初始分布为 UN(0)=N,分支机制为 P(UN(i)=a|UN(i1)=1)=(12)a+1, a0.

关于第一个结论,读者可以通过该引理前的一段话加以理解,严格证明请参见文献[]. 由(2.2)式,我们可以得到边际分布 P(UN(i)=a|UN(i1)=1)=b0(a+ba)parbp=(12)a+1, a0. 因此 UN 是一个单物种分支过程.

从0 到 τN 这段时间内,反射随机游动 |S| 的第一步是从 0 到 1, 最后一步是从 1 到 0. 所以在这段时间内,如果游动 |S| 在某个时间从 j 处到了 j+1 处,那么它必然要在将来某个时间回到 j. 因此, 可用过程 (UN,ZN) 来表示局部时过程 L(j;τN),j1. 我们有

{L(j,τN):j1}={UN(j1)+UN(j)+ZN(j):j1}. (2.3)
UN(x):=UN([Nx])/N, t0,ZN(x):=ZN([Nx])/N, x0. Lamperti[11] 证明 (UN)N1 的有限维分布收敛到连续状态分支过程,而Lindvall[12] 证明了过程(UN)N1 的胎紧性. 所以综合两者得到如下结论.

定理2.1 (Lamperti[11],Lindvall[12]) 作为空间 D([0,)) 上的一个过程, n,UNU, 其中 U 是随机微分方程

dzt=2z+tdBt,z0=1 (2.4)
的唯一解.

由(2.3)式有,对于任意 [Nx]1, lN(x)=UN([Nx])+UN([Nx]1)+ZN([Nx])N=UN(x)+UN(x1N)+ZN(x). 因此,欲证定理 1.1,只需证 {2UN(x)+ZN(x):x0}{lW(x,τ):x0}.

定义 XN(k):=2UN(k)+ZN(k),以及 XN(x):=2UN(x)+ZN(x).

本文将根据文献[1,定理15.6]证明结论. 为此, 我们给出如下三个命题. 前两个是关于有限维分布收敛的,最后一个是关于胎紧性的.

命题2.1 对于任意 α,β>0x0,当 N 时, αUN(x)+βZN(x)(α+βrp)U(x).

由引理 2.1,(UN(n),ZN(n))n0 是一个两物种分支过程. 根据其分支性和马氏性, E(eλ(αUN(n)+βZN(n))|UN(n1)=1)=m=0j=0eλ(αm+βj)P(UN(n)=m,ZN(n)=j|UN(n1)=1)=m=0j=0eλ(αm+βj)Cmm+jpmrjp=m=0j=0Cmm+j(eλαp)m(eλβr)jp=p1eλαpeλβr. 因此,对于固定的 x0, FN(x,λ;α,β):=Eeλ(αUN([Nx])/N+βZN([Nx])/N)=E(E(eλ(αUN([Nx])/N+βZN([Nx])/N)|UN([Nx]1)))=E(p1eλα/Npeλβ/Nr)UN([Nx]1). 由定理2.1, UN(x)U(x); 所以 EsUN(x)EsU(x). 再结合过程 (UN([Nx])) 的马氏性,有 EsUN([Nx])/N=E[(12s1/N)UN([Nx]1)]EsU(x). 注意到 sN:=(21eλα/Nqeλβ/Nrp)Neλ(α+rpβ), 于是,当 N 时, FN(x,λ;α,β)=E(12s1/Nn)UN([Nx]1)Eeλ(α+rpβ)U(x). 因此,当 N 时, αUN(x)+βZN(x)(α+rpβ)U(x). 证毕.

命题2.2 对于任意 k1,设 λi,αi,βi,i=1,,k 均为正实数. 那么对任意 0x1x2xk,当 N 时,

E(exp{ki=1λi(αiUN(x)+βiZN(x))})E(exp{ki=1λi(αi+βirp)U(xi)}). (2.5)
特别地,当 N 时,
E(exp{ki=1λiXN(x)})E(exp{ki=1λi(2+rp)U(xi)}). (2.6)
为了证明这一命题,我们需要用到连续状态分支过程的分支性,为了方便读者参考,我们将之陈述如下

引理2.2 {Y(t)}t0 是一个连续状态分支过程,Φ(t,λ)Y(t) 的拉普拉斯变换. 那么 E(exp{λY(t)}|Y(t1))=ΦY(t1)(tt1,λ),t1t. 一般地,对任意 λi0,i=1,,k,0t1t2tk, E(exp{ki=1λiY(ti)})=E(exp{k1i=1λiY(ti)}ΦY(tk1)(tktk1,λk)).

命题 2.2的证明 对于任意 λ,x0,α,β0, 令 FN(x,λ;α,β):=E(exp{λ(αUN(x)+βZN(x))}) 以及 F(x,λ;α,β):=E(exp{λ(α+βrp)U(x)}).

我们采用数学归纳法证明. 事实上,当 k=1 时对应的结论已在命题2.1中被证明.

假设当 k=m 时,(2.5) 式成立. 下证 k=m+1 时结论成立. 注意到 (UN(0),ZN(0))=(N,0). 由过程 (UN(k),ZN(k))k1 的分支性,有 E(exp{λm+1(αm+1UN(xm+1)+βm+1ZN(xm+1))}|F[Nxm])=FUN(xm)N([Nxm+1][Nxm]N,λm+1;αm+1,βm+1). 因此,

E[exp{m+1i=1λi(αiUN(xi)+βiZN(xi))}]=E[E(exp{m+1i=1λi(αiUN(xi)+βiZN(xi))}|F[Nxm])]=E[exp{mi=1λi(αiUN(xi)+βiZN(xi))}FUN(xm)N([Nxm+1][Nxm]N,λm+1;αm+1,βm+1)]=E[exp{mi=1λi(αiUN(xi)+βiZN(xi))}{FUN(xm)N(xm+1xm,λm+1;αm+1,βm+1)FUN(xm)(xm+1xm,λm+1;αm+1,βm+1)}]+E[exp{mi=1λi(αiUN(xi)+βiZN(xi))}FUN(xm)(xm+1xm,λm+1;αm+1,βm+1)]. (2.7)
由命题 2.1,可见当 N 时, FN([Nxm+1][Nxm]N,λm+1;αm+1,βm+1)F(xm+1xm,λm+1;αm+1,βm+1). 易见对于 y0 一致地有
FyN([Nxm+1][Nxm]N,λm+1;αm+1,βm+1)Fy(xm+1xm,λm+1;αm+1,βm+1). (2.8)
所以当 N 时,(2.7) 式的右端趋于 0,由归纳假设,第二项收敛到 E(exp{mi=1λi(αi+βirp)U(xi)+lnF(xm+1xm,λm+1;αm+1,βm+1)U(xm)}), 由引理 2.2,它等于 E(exp{m+1i=1λi(αi+βirp)U(xi)}). 所以 k=m+1 时结论成立. 因此 (2.5) 式成立.

在(2.5)式中,对于所有 i,令 αi=2,βi=1,可得 (2.6)式. 证毕.

命题2.3 对于固定的 T>0,存在 α>1,β0, 以及 C>0 使得对于任意 0x1xx2T, E(|XN(x)XN(x1)|β|XN(x2)XN(x)|β)C|x2x1|α.

为证该命题,我们首先证明如下引理.

引理2.3 假设 i,j,k 均为非负整数满足 ijk. 设 β=43,则存在一个常数 L< 使得 E(|XN(k)XN(j)|β|XN(j)XN(i)|β)LN3/4β(k+2i)βjβ/4.

证明中,我们用 K 来表示某一常数,它独立于其他参数,且其具体取值在不同的地方可能会不同. 设 FNk 为由 {(UN(i),UN(i)):ik} 所生成的 σ -代数. 因为当 k1 时,XN(k)=2UN(k)+ZN(k),所以 E(|XN(k)XN(j)|2|FNj)=E(|2UN(k)2UN(j)+ZN(k)ZN(j)|2|FNj)8E[(UN(k)UN(j))2|Fj]+2E[(ZN(k)ZN(j))2|FNj].

σ2:=U1(1),以及 σ21:=Z1(1). 由过程 (UN(k))k1 的马氏性和分支性,有 E[(UN(k)UN(j))2|FNj]=UN(j)(U1(kj))=(kj)σ2UN(j). 注意到分支过程 (UN(k))k1 只有其第一个物种产生后代,所以 E[(ZN(k)ZN(j))2|FNj]=E[(UN(j)i=1(Z1i(kj)rp)+(rpUN(j)ZN(j)))2|FNj]=[σ21+(rp)2(kj1)σ2]UN(j)+[rpUN(j)ZN(j)]2, 其中 Z1i(kj),i1Z1(kj) 的独立同分布的复制品,并且它们均与 UN(j) 相互独立. 所以存在常数 K>0 使得

E(|XN(k)XN(j)|2|FNj)K(kj)UN(j)+(rpUN(j)ZN(j))2. (2.9)
由 Hölder 不等式, E(|XN(k)XN(j)|β|FNj)Eβ/2(|XN(k)XN(j)|2|FNj). 所以结合(2.9)式,有 E(|XN(k)XN(j)|β|XN(j)XN(i)|β)=E[|XN(j)XN(i)|βE(|XN(k)XN(j)|β|FNj)]E[|XN(j)XN(i)|βEβ/2(|XN(k)XN(j)|2|FNj)]E{|XN(j)XN(i)|β[K(kj)UN(j)+(rpUN(j)ZN(j))2]β/2}. 现在,令 q1=32,q2=3, 那么1q1+1q2=1,1q1=β2,1q2=β4. 由 Hölder 不等式,有
E{|XN(j)XN(i)|β[K(kj)UN(j)+(rpUN(j)ZN(j))2]β/2}Eβ2(|XN(j)XN(i)|2)Eβ4[K(kj)UN(j)+(rpUN(j)ZN(j))2]2. (2.10)
由(2.9)式有, E(|XN(k)XN(j)|2)E[K(kj)UN(j)+(rpUN(j)ZN(j))2]. 因为
E[(rpUN(i)ZN(i))2]=E[(rpUN(i1)m=1(U1m(1)1))+E(UN(i1)l=1(Z1l(1)rp))]22E(rpUN(i1)m=1(U1m(1)1))2+2E(UN(i1)l=1(Z1l(1)rp))2, (2.11)
其中 U1m(1),m1U1(1) 的独立同分布的复制品,且均与 UN(i1) 相互独立; Z1l(1),l1Z1(1) 的独立同分布的复制品, 且均与 UN(j1) 相互独立. 于是 E[(rpUN(i)ZN(i))2]2σ2(rp)2E(UN(i1))+2σ21E(UN(i1))KN. 所以,
E[K(ji)UN(i)+(rpUN(i)ZN(i))2]K(ji+1)N. (2.12)
类似于 (2.11)式的推导,可得 E[K(kj)UN(j)+(rpUN(j)ZN(j))2]22K2(kj)2E(UN(j))2+2E(rpUN(j)ZN(j))42NK2(kj)2jσ2+KNjKN(kj+1)2j, 最后两个式子中 K 的值需要做适当调整. 结合 (2.10) 和 (2.12)式,可得存在常数 L>0 使得 E(|XN(k)XN(j)|β|XN(j)XN(i)|β)LN3/4β(k+2i)βjβ/4. 证毕.

命题 2.3 的证明α=β=43, 设 0x1xx2T,N1.

x2x11/N 时,[Nx2][Nx1]1,所以 [Nx]=[Nx1] 或者 [Nx]=[Nx2],所以此时 E(|XN(x)XN(x1)|β|XN(x2)XN(x)|β)=E((|XN([Nx])XN([Nx1])|β|XN([Nx2])XN([Nx])|β)N2β)=0C|x2x1|β,x2x11/N 时,由引理 2.3,有 E(|XN(x)XN(x1)|β|XN(x2)XN(x)|β)=E((|XN([Nx])XN([Nx1])|β|XN([Nx2])XN([Nx])|β)N2β)LN3/4r([Nx2][Nx1])β([Nx])β4N2βC|x2x1|β, 这里 C 依赖于常数 T. 所以命题 2.3 得证.

定理 1.1 的证明 注意到过程 U 轨道连续, 又有命题2.2和命题2.3,于是根据文献[1,定理15.6] 可知在空间 D([0,T])2UN(x)+ZN(x)(2+rp)U(x). 再由 文献[13,定理3'],这一收敛亦在空间 D([0,)) 上成立. 由 Ray-Knight 定理,文献[18] 得到 2U(x)=l(x,ρ)=l(x,ρ)+l(x,ρ), 其中 l(x,t) 是标准布朗运动的局部时,ρ=inf{t0:l(0,t)>1}. 注意到 (2+rp)=2σ2. 于是 (2+rp)U(x)=σ2l(x,ρ). 由局部时的定义,易知 lW(x,t)d=σ2l(x,σ2t). 再注意到上面 ρ 的定义, 以及定理 1.1 中 τ 的定义,有 (lW(x,t),τ)d=(σ2l(x,σ2t),σ2ρ). 所以,(2+rp)U(x)d=lW(x,τ). 定理得证.

参考文献
[1] Billingsley P. Convergence of Probability Measures. New York:John Wiley Sons, 1968
[2] Borodin A N. The asymptotic behavior of local time of recurrent random walk with finite variance. Theory Prob Appl, 1981, 26:758-772
[3] Borodin A N, Ibragimow I A. Limit theorems for functionals of random walks. Proceeding of the Steklov Institute of Mathematics, 1995, 195(2):259 pages
[4] Csörgö M, Revész P. Three strong approximations of the local time of a Wiener process and their applications to invariance. Limit Theorems in Probability and Statistics, 1982, 1:223-254
[5] Csörgö M, Revséz P. On strong invariance for local time of partial sums. Stoch Proc Their Appl, 1985, 20:59-84
[6] Dwass, M. Branching processes in simple random walk. Proc Amer Math Soc, 1975, 51:270-274
[7] Ispány M, Pap G. Asymptotic behavior of critical primitive multi-type branching processes with immigration. 2012, arXiv:1401.3440v1[math.PR]
[8] Kesten H, Kozlov M V, Spitzer F. A limit law for random walk in a random environment. Comp Math, 1975, 30:145-168
[9] Knight F B. Random walks and a sojourn density of Brownian motion. Trans Amer Math Soc, 1963, 109:56-86
[10] Kurtz T G. Diffusion approximations for branching processes. Branching Processes Adv Prob, 1979, 5:262-292
[11] Lamperti J. The limit of a sequence of branching processes. Z Wahr Schein-lichkeisth, 1967, 7:271-288
[12] Lindvall T. Convergence of critical Galton-Watson processes. J Appl Prob, 1972, 9:445-450
[13] Lindvall T. Weak convergence of probability measure and random functions in the function space D([0;1)). J Appl Prob, 1973, 10:109-121
[14] McKean H P. Hölder condition for Brownian local time. J Math Kyoto Univ, 1962, 1:196-201
[15] Perkins E. A global intrinsic characterization of local time. Ann Prob, 1981, 9:800-817
[16] Perkins E. Weak invariance principle for local time. Prob Theory Rel Fields, 1982, 60:437-451
[17] Ray D B. Sojourn times of diffusion processes. Ill J Math, 1963, 7:615-630
[18] Rogers L C G. Brownian local times and branching processes//Picard J. Lecture Notes in Math 1059. Berlin:Springer, 1982:42-55
[19] Wang H M. Mean and variance of first passage time of non-homogeneous random walk. Front Math China, 2012, 7(3):551-559
[20] Zeitouni O. Random walks in random environment//Picard J. Lecture Notes in Math 1837. Berlin:Springer, 2004:189-312
带停留对称随机游动局部时的渐进行为
孙鸿雁, 么亚楠