Donsker 建立了随机游动的不变原理,即其标准化后的极限为布朗运动. 具体地, 设{ξn:n≥0} 是直线上一列独立同分布的随机变量, 满足对于任意的n,ξn的期望 Eξn=0及方差 ξn=σ2,其中 σ>0. 定义随机游动{Sn:n≥0}为 Sn:=n∑i=1ξi,n≥1,S0=0, 及其线性内插S(t): t≥ 0为 S(t):=S[t]+(t−[t])(S[t]+1−S[t]),∀t≥0. 记相应的尺度过程为 Sn(t)=S(nt)√n, ∀t≥0, Donsker 证明了过程 {Sn(t):t≥0} 弱收敛到布朗运动 {W(t):t≥0}, 且方差 W(t)=σ2t.
受Donsker 不变原理启发,关于随机游动局部时的不变原理引起了人们的广泛关注. 记随机游动的局部时过程为 phi(x,n),即 phi(x,n) 为 {Sn} 前 n 步中到过 x 处的次数. 对于紧邻随机游动,即 ξn 分别以 1/2的概率取 ±1, Borodin[2] 借助于不变嵌入方法建立了局部时的不变原理,即: 两指标的过程 n−1/2phi([x√n],[nt]) 收敛到布朗运动 {W(t):t≥0} 在时刻 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} 为带停留的紧邻随机游动是指
定义反射过程 |S| 的局部时为 L(j,n):=♯{r:1≤r≤n,|Sr|=j},j≥0,n≥1. 令 τ0≡0,对于任意 m≥1, 设 τm:=inf{n>τm−1:Sn=0,Sn−1≠0}. 用 D(E) 表示空间E上左极右连函数全体并赋予Skorohod 拓扑. 对任意自然数 N≥1,通过局部时做尺度变换,我们定义空间 D[0,∞) 上的过程
本文目标是证明当 N→∞ 时,lN(⋅) 在空间 D[0,∞) 上收敛. 即
定理1.1 作为空间 D[0,∞) 上的一个元素,当 N→∞ 时, lN(⋅)⇒l∗W(⋅,τ), 其中 l∗W(x,t)=lW(x,t)+lW(−x,t),τ=inf{t>0:lW(0,t)>σ−2}. ``⇒" 表``弱收敛".
注 1.1 当 r=0 时,σ2=1,对应的正是Rogers[1]的结果.
Kesten 等人[8]和Wang[19]揭示了带停留的紧邻随机游动的内蕴分支机制; 该机制对应了一个两物种分支过程,其均值矩阵是非正则的. 注意到, 与停时相关的局部时可表示成这一两物种分支过程的泛函. 本文将借助于这一表达来证明停时相关的局部时的收敛性. 关于多物种分支过程的收敛性,可参见文献[7, 10], 但文献中都是假定均值矩阵是正则的. 本文要处理的是非正则的情形.
定理 1.1 的证明将在下一节展开.
首先,我们来回顾带停留随机游动的内蕴分支机制,具体可参见文献[]. 设 (UN(0),ZN(0))=(N,0). 对于任意 i≥1, 定义
Zeitouni[20] 指出过程 |S| 的每个从 i 到 i+1 或者到 i 的步子必然出现在某个连续的从 i−1 到 i 的步子之间. 并且在每个连续的从 i−1 到 i 的步子之间的这两种步子服从多项分布. 因此,在(UN(i−1),ZN(i−1)),(UN(i−2),ZN(i−2)),⋯,(UN(0),ZN(0)) 生成的σ-代数条件下,随机向量 (UN(i),ZN(i)) 同分布于 UN(i−1) 个独立同分布的随机变量之和, 这些随机变量服从如下分布 P(V1=(a,b))=(a+ba)parbp, a,b≥0. 也就是说,(UN(n),ZN(n))n≥0 是一个两物种分支过程,初始分布为 (UN(0),ZN(0))=(N,0). 下面我们精确刻画这一结论.
引理2.1 (i) (2.1) 式定义的过程 (UN(n),ZN(n))n≥0 是一个两物种分支过程. 其后代分布为
(ii) 过程 (UN(n))n≥0 是一个单物种分支过程,初始分布为 UN(0)=N,分支机制为 P(UN(i)=a|UN(i−1)=1)=(12)a+1, a≥0.
证 关于第一个结论,读者可以通过该引理前的一段话加以理解,严格证明请参见文献[]. 由(2.2)式,我们可以得到边际分布 P(UN(i)=a|UN(i−1)=1)=∑b≥0(a+ba)parbp=(12)a+1, a≥0. 因此 UN 是一个单物种分支过程.
从0 到 τN 这段时间内,反射随机游动 |S| 的第一步是从 0 到 1, 最后一步是从 1 到 0. 所以在这段时间内,如果游动 |S| 在某个时间从 j 处到了 j+1 处,那么它必然要在将来某个时间回到 j. 因此, 可用过程 (UN,ZN) 来表示局部时过程 L(j;τN),j≥1. 我们有
定理2.1 (Lamperti[11],Lindvall[12]) 作为空间 D([0,∞)) 上的一个过程, 当n→∞,UN⇒U, 其中 U 是随机微分方程
由(2.3)式有,对于任意 [Nx]≥1, lN(x)=UN([Nx])+UN([Nx]−1)+ZN([Nx])N=UN(x)+UN(x−1N)+ZN(x). 因此,欲证定理 1.1,只需证 {2UN(x)+ZN(x):x≥0}⇒{l∗W(x,τ):x≥0}.
定义 XN(k):=2UN(k)+ZN(k),以及 XN(x):=2UN(x)+ZN(x).
本文将根据文献[1,定理15.6]证明结论. 为此, 我们给出如下三个命题. 前两个是关于有限维分布收敛的,最后一个是关于胎紧性的.
命题2.1 对于任意 α,β>0 和 x≥0,当 N→∞ 时, αUN(x)+βZN(x)⇒(α+βrp)U(x).
证 由引理 2.1,(UN(n),ZN(n))n≥0 是一个两物种分支过程. 根据其分支性和马氏性, E(e−λ(αUN(n)+βZN(n))|UN(n−1)=1)=∞∑m=0∞∑j=0e−λ(αm+βj)P(UN(n)=m,ZN(n)=j|UN(n−1)=1)=∞∑m=0∞∑j=0e−λ(αm+βj)Cmm+jpmrjp=∞∑m=0∞∑j=0Cmm+j(e−λαp)m(e−λβr)jp=p1−e−λαp−e−λβr. 因此,对于固定的 x≥0, FN(x,λ;α,β):=Ee−λ(αUN([Nx])/N+βZN([Nx])/N)=E(E(e−λ(αUN([Nx])/N+βZN([Nx])/N)|UN([Nx]−1)))=E(p1−e−λα/Np−e−λβ/Nr)UN([Nx]−1). 由定理2.1, UN(x)→U(x); 所以 EsUN(x)→EsU(x). 再结合过程 (UN([Nx])) 的马氏性,有 EsUN([Nx])/N=E[(12−s1/N)UN([Nx]−1)]→EsU(x). 注意到 sN:=(2−1−e−λα/Nq−e−λβ/Nrp)N→e−λ(α+rpβ), 于是,当 N→∞ 时, FN(x,λ;α,β)=E(12−s1/Nn)UN([Nx]−1)→Ee−λ(α+rpβ)U(x). 因此,当 N→∞ 时, αUN(x)+βZN(x)→(α+rpβ)U(x). 证毕.
命题2.2 对于任意 k≥1,设 λi,αi,βi,i=1,⋯,k 均为正实数. 那么对任意 0≤x1≤x2≤⋯≤xk,当 N→∞ 时,
引理2.2 设 {Y(t)}t≥0 是一个连续状态分支过程,Φ(t,λ) 是 Y(t) 的拉普拉斯变换. 那么 E(exp{−λY(t)}|Y(t1))=ΦY(t1)(t−t1,λ),t1≤t. 一般地,对任意 λi≥0,i=1,⋯,k,0≤t1≤t2≤⋯≤tk, E(exp{−k∑i=1λiY(ti)})=E(exp{−k−1∑i=1λiY(ti)}ΦY(tk−1)(tk−tk−1,λk)).
命题 2.2的证明 对于任意 λ,x≥0,α,β≥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))k≥1 的分支性,有 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). 因此,
在(2.5)式中,对于所有 i,令 αi=2,βi=1,可得 (2.6)式. 证毕.
命题2.3 对于固定的 T>0,存在 α>1,β≥0, 以及 C>0 使得对于任意 0≤x1≤x≤x2≤T, E(|XN(x)−XN(x1)|β|XN(x2)−XN(x)|β)≤C|x2−x1|α.
为证该命题,我们首先证明如下引理.
引理2.3 假设 i,j,k 均为非负整数满足 i≤j≤k. 设 β=43,则存在一个常数 L<∞ 使得 E(|XN(k)−XN(j)|β|XN(j)−XN(i)|β)≤LN3/4β(k+2−i)βjβ/4.
证 证明中,我们用 K 来表示某一常数,它独立于其他参数,且其具体取值在不同的地方可能会不同. 设 FNk 为由 {(UN(i),UN(i)):i≤k} 所生成的 σ -代数. 因为当 k≥1 时,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))k≥1 的马氏性和分支性,有 E[(UN(k)−UN(j))2|FNj]=UN(j)(U1(k−j))=(k−j)σ2UN(j). 注意到分支过程 (UN(k))k≥1 只有其第一个物种产生后代,所以 E[(ZN(k)−ZN(j))2|FNj]=E[(UN(j)∑i=1(Z1i(k−j)−rp)+(rpUN(j)−ZN(j)))2|FNj]=[σ21+(rp)2(k−j−1)σ2]UN(j)+[rpUN(j)−ZN(j)]2, 其中 Z1i(k−j),i≥1 是 Z1(k−j) 的独立同分布的复制品,并且它们均与 UN(j) 相互独立. 所以存在常数 K>0 使得
命题 2.3 的证明 令 α=β=43, 设 0≤x1≤x≤x2≤T,N≥1.
当 x2−x1≤1/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β)=0≤C|x2−x1|β, 当 x2−x1≥1/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|x2−x1|β, 这里 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{t≥0: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=l∗W(x,τ). 定理得证.