扩散过程平稳分布的存在唯一性与稳定性

应用数学

MATHEMATICAAPPLICATA

2010,23(4):713—718

扩散过程平稳分布的存在唯一性与稳定性

徐立峰

(湖北师范学院数学与统计学院,湖北黄石435002)

摘要:用耦合方法研究扩散过程稳定性.利用KRW概率距离的对偶表示,获得平稳

分布的存在唯一性、遍历性和依分布稳定性并且给出了收敛速度估计。

关键词:扩散过程;随机微分方程;平稳分布;依分布稳定;耦合

中图分类号:0211.63

文献标识码:AAMS(2000)主题分类:60j60;60H10文章编号:1001—9847(2010)04—0713—06

1.预备知识与引理

设X:是R。上的扩散过程,它满足如下It6型随机微分方程

dX,一6(X。)dt+仃(X,)dB,,(1.1)

式中B,是某概率空间(n,易P)上d维Weiner过程,6(・):R。一R。是d维Borel可测函数,盯(・):R。一R以。是d×dBorel可测函数矩阵.对z∈R。和盯∈R积。定义范数

.Zd2(蚤h|z)l/2H小(i善1%12)1/2・

用P(t,z,dy)表示扩散过程X,的转移概率,简记为P。.(1.1)对应的扩散算子

(卯(z)=pd(z)矗弛)+芝1一i暑d“z)为i弛w∈a(㈣.(1.2)

简记为以“(z),6(z)),其中n(z)=d(z)T(z).

设P・,Pz∈缎R。)(R。上全体概率测度),p∈缄R。×R。),称∥是弘,与岸。的耦合,如果满足边缘性:

P(Al×R。)一卢l(A1)(A1∈∥),卢(R。×A2)一P2(A2)(A2∈∥).

用K(户-,∥z)表示P-与弘z的全体耦合概率,不致混淆时简记为K.对R。上的距离函数p(x,y)一lz—Yl(z,y∈R。)定义P,与户2的KRW概率距离:

w(∥l,户2)=i,n!Ip(x,y)产(dz,dy).

若存在皿∈K(p。,卢2),使

w(pl,胆)一Ip(x,y)fi(dx,dy),

・收稿日期:2009~09—08

基金项目:教育部重点科研项目(209078)作者简介:徐立峰.男.汉族,上海人,讲师。研究方向:概率论与数理统计.

714应用数学2010称乒是.£‘。与户。的最优耦合,最优耦合总是存在的,且此结论可推广到户-,卢z为转移概率的情形‘¨.

对丌∈缄掣),r,P。(-):=In(dx)P(t,z,・).若V£>o,7cP。=硝,则称丌是X,的平稳分布.若对Vz∈R。,当t—+。。时P(t,z,dy)弱收敛于某平稳分布丌(ay),称X依分布稳定.‘

KRW概率距离W(卢。,肛2)一般地强于弱收敛意义下的概率距离ct,(p,,pz)(即Levy—Pro—horov距离),事实上我们有W(卢,,p。)≥叫2(∥z,脾)‘21,因此耦合方法成为研究过程依分布稳定的一个可供选择的工具.

任给厂∈C(R。),定义w)=磐岭掣,驴=(r

Ec(妁,w)<oo)-

L(,)叫做,的Lipschitz常数.

对参考点z。∈R4,记

绲=(卢∈苏掣),lp(x,zo)卢(出)<o。).

注意,若Ⅳ∈玩,则VYo∈R4,(1.3)

Ip(z,Y。)卢(dLc)≤I[P(z,zo)+p(xo,yo)]产(dz)

=lp(z,zo)p(dz)+p(xo,Yo)<∞,

因此z。是任意的.

若7r是X,的平稳分布,且对V户∈玩,W(∥,(・),丌(・))一O(t—co),称X,在KRW概率距离意义下是遍历的.本文中简称为遍历.若存在A>0,使V卢∈玩,Ⅳ(胪。(・),丌(・))≤e--知W(p(・),丌(・)),则称X,指数遍历.

引理1.1C23Vpl,p2∈琨,V矿(胁,脾)一sup{lld(胁一re)f(Y)f,厂∈乡,L(D≤1).

引理1.1被称为KRW距离的对偶表示.

对于一般Markov过程平稳分布的存在唯一性及遍历性有如下一般性准则.

引理1.2[妇设x,(f≥o)是距离可测空间(E,lD,o上的Markov过程,P(t,z,dy)是其转移概率,满足

(i)VT>0,存在0≤Cl(T),c2(T)<∞,使

IP(t,z,dy)p(y,zo)≤c1(丁)+cz(T)lD(z,zo),£∈[o,T3,

其中z。是E中任意固定的一个参考点;(1-4)

(ii)对t∈[o,T3,存在0≤c(£)<1,使

w(P(t,zl,・),P(t,X2,・))≤c(t)p(xl,z2),zl,X2∈E.

则X。存在唯一平稳分布7r∈玩,且X。是遍历的.

本文中我们将用耦合方法着重讨论(1.1)依分布稳定性,不同的是我们采用了对偶原理,它的好处是将研究对象从测度与算子转化为普通实函数,相对简单.本文利用扩散半群P:对Lipschitz常数的压缩性刻划了x,的极限性质,同时得到X。平稳分布的存在唯一性、遍历性及依分布稳定性,并且给出了收敛速度估计,而同类结果中有些需要预先假定平稳分布的存在性,定理条件也弱于某些已有的研究.(1.5)

第4期徐立峰:扩散过程平稳分布的存在唯一性与稳定性7152.主要结果与应用

定理2.1若方程(1.1)的系数满足线性增长条件

6(z)I2+I盯(z)I2≤K(1+l

且对Vt>0,存在0≤r(£)<1,使

L(P。厂)≤,.(£)L(,),厂∈g。,(2.2)z2),z∈RJ,(2.1)

(i)扩散过程X。有唯一的平稳分布丌∈绲,且X,遍历,因而也依分布稳定;

。(ii)V∥∈甥,

W(zP(£,z,・),丌(・))≤r(t)W(F(・),7r(.)).

证首先证明(i),这只需验证引理1.2的两个条件成立.

广(2.3)定义由方程(1.1)确定的期望算子

EⅡ,(x,)]一Jf(y)p(t,z,dy)=P,厂(z).

由T(1.1)的系数满足线性增长条件(2.1),依文[4]引理21.2有矩估计:对VT>0,存在c(丁)>o,使

EzX≤,L1+Z\,酽nzJ∈R£∈[O明,L24

由TJenSen不等式

E21Ⅸrl≤LPI叉r

≤[(1+fz2J“2≤L(1+Iz2)e“7’]“2(2.5)f)2・e。‘D]172一(1+jz})・ec(n/2.

在(1.4)中取参考点z。一0,则

j.P(£,z,dy)P(y,o)一P妒(z,o)一Ez

即此(1.4),再证(1.5).

W(P(t,zl,・),P(t,z2,・))x。I≤(1+lz1)ef‘n/2:一c1(丁)+c2(T)J0(2,。).由引理1.1,并注意到五是Feller过程(以(・)表示z处的点测度)。

一sup{lJ[PQ,zz,dy)一P(t,锄,d了)]厂(y)I,,∈箩,L(厂)≤1)

一sup{Ij以。Pt(d3,)厂(y)一以:Pr(dy)厂(y’I,厂∈驴,L(厂)≤1}

一sup{lJ.乱。cdz,[JIP(t,x,dy),cy,]~.f%cdz,P(t,x,dy)厂cy,]I,厂∈箩,Lc厂,≤-)2sup{lj良,(dz)Pr厂(z)一j以。(dx)Pr,(z)I,,∈g。,L(厂)≤1)

一sup{fJ文。(dx)Pr,(z)一j以。(dy)P。厂(y)f,fE箩,L(厂)≤l}.

记卢是艿1与屯的联合分布,则五(刚,A)=以,(A),互(A,R。)=如(A),于是

上式=sup

≤sup”厂(嘶(dx,W)--P,,(y)压(Ru,dy)f,,∈箩,L(p≤1},JP,f(x)--P,厂(了)f卢(如,dy),厂∈妙,L(厂)≤1}

supIfL(Ⅳ)p(x,y)fi(dx,dy),厂∈矽,L(厂)≤1}

sup{fr(£)L(厂)ID(z,y)fi(dx,dy),fE驴,L(厂)≤ll

716应用数学20lO

≤r(£)Ip(x,y)犀(dx,dy)=r(t)p(xl,z2).(2.6)

由引理1.2,X。有唯一平稳分布丌,且对V卢∈绲,w(∥。(・),zr(・))一o(当t一∞),特别地,取肛=乱,则

W(&P,(・),7f(・))一W(P(£,X,・),丌(・))一0

由于按KRW距离收敛强于弱收敛,故X。依分布稳定.

再证(2.3)式,由X,对丌的不变性,

W(肛P(£,z,・),zr(・))一w(£一∞),(f,工,・),护(z,o,・))

=Sup…矿。江,dy),(Y)一J护(∽,dy)厂(Y)l,,∈驴,L(,)≤l}

一p㈣卜(dz)P(t,z,dy)p一『[』rr(dz)P(t,z,dy)p|',∈驴郴≤・)一suP{lI户(dz)P。,(z)~l7f(dz)P。厂(z)l,,∈g’,L(厂)≤11.1J√

设P(・,・)为岸(・)与丌(・)的最优耦合,与(2.6)类似地可证

w(pP(£,z,・),丌(・))≤,.(f)Ip(z,了)P(dz,dy)一r(t)W(ff,丌).

注2.1设X,为Wiener过程,其转移概率为P(t,士,dy),则线性增长条件(2.1)显然成立(6—0,d==l,rx。),且由文[2]知L(只,)≤L(厂),但熟知Wiener过程没有平稳分布,这说明定理条件(2.2)中的r(£)<1不可放宽为,.(£)≤1,即条件(2.2)是精确的.

注2.2(2.3)反映了扩散半群P,关于Lipschitz常数压缩性与关于KRW距离压缩性的密切联系,由于测度的弱收敛距离锄(卢-,P2)≤W172硼(户1,P2),且w(以(・),丌(・))=Ip(x,y)丌(dy)<f)o(丌(・)∈绲).因此(2.3)也给出了弱收敛速度的估计r172(f)(当r(£)一o).特别地当r(£)一e咄(A>o)时,X,指数遍历,因而也依分布指数稳定.

记妖a(x,y),b(x,y))是颤n(z),6(z))与自身的耦合算子,即

以为妒2,、k,≯“0))’/“(z),c(x,y)、

其中““y)一(c一(而y”以。使矩阵d。,y)非负定,灰z,y)一(::i;).

推论2.1若定理2.1中(2.1)式成立,且存在g与自身的耦合算子勿和常数口>0,使

埤(z,y)≤~qD(z,y),z,Y∈R。,

则X,有唯一平稳分布,r(・)且指数遍历,因而也依分布指数稳定.

可推出(2.7)我们只需验证耦合算子压缩条件(2.7)强于(2.2).事实上根据文[73定理2.3,由(2.7)

E“’ID(X,,Y,)≤e--'p(z,y),

其中yf是由Y出发的解过程.记T—inf{t,X,-=Y,).

P。厂(z)~P,I厂(y)

z~Y(2.8):!.芝[』!墨!!]二旦[丛茎业一!二卫!羔上二f!∑!!!Iz—Y}z—yIT>t]

姑’yJ一』鬻趔-I:Y,匕IX“Yu・1鲁掣YI‰z~l[TM。

≤w炉j-I等并bd≤e气(n

第4期徐立峰:扩散过程平稳分布的存在唯一性与稳定性717于是L(P,厂)≤e1L(厂),因口>0,由定理2.1,X。有唯一平稳分布丌,且由(2.3)X,指数遍历,特别地,(2.3)中取/z一以,则Ⅳ(P(t,X,・),丌(・))≤e1W(艿:(・),丌(・)),即X,依分布指数稳定,指数阶为口/2.

推论2.1明显地优于文[5]中的结果,后者需事先假定平稳分布存在,也优于文f7]中定理2,这里证明了指数遍历性,因而也给出了收敛速度估计.而文[7]在相同条件下,仅得到遍历性.

注2.3在注2.1中我们指出条件L(P,厂)≤L(厂)一般地不能得到平稳分布的存在性,但如果我们假设平稳分布7l"存在唯一,类似于(2.3)的证明仍然可以得到W(P(t,z,・),丌(・))≤w(以(・),丌(・)),此即文[5]中的平均稳定性.

由于Lipschitz常数的压缩性在许多领域有重要应用,因此也为定理2.1的应用提供了一些现成的结果,Herbst和Pitt在研究扩散半群的单调性时给出了如下结果.

引理2.1【83如果b,(z)∈C1(R。)且Ⅱ。(z)仅依赖于zi和zJ(1≤i,J≤d),则

L(P,厂)≤eKtL(,),f∈驴,t≥0,

其中

K翟p罡慧{掣+萎I掣h

考虑如下Hopfield随机神经网络

dX,一[一AX,+Bg(X。)]出+盯(X。)dB,,(2.9)

其中A—diag(aj,d2,…,(1d),c£i>0,1≤i≤d,B一(%)以d,g(z)=(gl(A),92(zz),…,gJ(zd))T,仃(z)一(盯i(z))‘,×d.假定(2.9)有唯一解.

定理2.2在方程(2.9)中,若gi(z)可导,盯(z)=diag(a1(z1),…,O'd(zd)),gf(z),吼(z)满足线性增长条件.

K=sup…ma…x(一ai+b.g,7(zi)+∑J1≤i≤d、:Z』bjigi7(zi)I),J

则当K<0时,(2.9)有唯一平稳分布且依分布指数稳定.

证这足引理2.1与定理2.1的直接结果,且此时c(£)一e触,故指数阶为一K/2.

也可通过构造耦合计算(2.2),利用文[9]中构造的反射耦合我们有

H--=-s。≠up,南{l

论成立.定理2.3在随机神经网络(2.9)中,令d(z)一盯(y)I2一L鱼鱼堕-_鼍望笔掣+<一A(x--y),X--y>+(Big(x)~g(y)],z—y>},则当H<0时,(2.9)有唯一平稳分布且依分布指数稳定(其中(-,・>表示R。上的内积).由文[9]定理1并经简单计算,此时有L(P,厂)≤eHtL(厂),f∈2。,故由定理2.1结

例2.1考虑三维随机神经网络

rdX-(£)=[一X1(£)一0.49l(X1(£))+0.592(Xz(£))+0.69。(K(£))]出+2.4Xl(t)dBl(£),Jd%(£)一[一1.1x2(£)+1.19l(x1(£))一0.692(x2(£))+0.29。(x3(£))]出+2.5托(t)dB2(£),Idx。(f)=[一0.8x。(£)+0.39l(x1(£))+0.792(X2(f))一0.59s(x。(£))]出+2.4X。(t)clB。(£),其中Ⅱ向应函数g。(“)=薹{芝三,i=l,2,3.由于此时o<g7(“)=1一eu_e-U)2<1.易知

718应用数学20lO方程有唯一解,且不难得到K=max(一1+0.7,一1.1+0.7,一0.8+0.5)一一0.3.由定理2.2,X(£)依分布指数稳定,指数阶为0.15.

例2.2考虑二维随机神经网络

,dXl(£)=[一0.9Xl(£)+1.29(X2(£))]d£+0.2XI(t)dWl(£),

其中响应函数g(“)满足Lipschitz条件,如g(摊)=L型_』掣或g(“)一雨eu_e-u.1dx2(£)一[一1.1x2(f)+0.49(Xl(t))]dt+0.25X2(£)dW2(£).

通过计算不难得到上界估计:H<0.252—0.22/2—0.9+1.6/2一一0.0575<0,由定理2.3,X(£)依分布指数稳定,指数阶为o.028.

参考文献:

[13张绍义,徐侃.转移概率最优可测耦合的存在性[J].数学学报,1997,40(1):5-13.

[z]CHENMufa.FromMarkov

ic,1996.ChainstONon-equilibriumParticleSystems[M].Singapore:WorldSeientif-

[3]张绍义.最优可测耦合的存在性与Markov过程的遍历性[J].中国科学(A辑),1998,28(11):999—1008.[4]黄志远.随机分析学基础[M].2版.jE京:科学出版社,2001.

[5]张绍义.耦合方法与随机微分方程平均稳定性[J].北京师范大学学报:自然科学版。1999,35(3):318—

322.

[6]CHENMufa,LIShaofu.Couplingmethodsformultidimensionaldiffusionprocesses[J].Ann.Prob.,

1989。17(1):151-177.

[7]徐侃,张绍义.Markov耦合与Markov过程的遍历性[J].数学杂志,2001,21(3):315—318.

[8]HERBSTI。DITTL.Diffusionequationtechniquesinstochasticmonotonicityandpositivecorrebtions

TheoryRelatFielbs,1991(87):275—312.[J].Probability

[9]王风雨.扩散半群Lipschitz常数的压缩性[J].北京师范大学学报:自然科学版,1996,32(4):447—451.[10]徐侃,徐立峰.一般状态空间跳过程不变测度的存在性和唯一性[J].数学杂志。2004,24(5):561—564.

ExistenceandUniquenessoftheStationaryDistributionand

StabilityinDistributionforSomeDiffusionProcesses

XL厂Lifeng

(College

na)ofMathematicsandStatistics,HubeiNormalUniversity,Huangshi435002.Chi—

Abstract:Couplingmethodisused

thedualityexpressionoftostudythestabilityofsomediffusionprocesses.ByofstationaryKRW—metricweprovetheexistenceanduniqueness

ratedistribution,andtheestimationofconvergenceisobtainedtoo.

Keywords:Diffusionprocesses;Stochasticdifferentialequation;Stationarydistribution;Stabilityindistribution;Coupling

扩散过程平稳分布的存在唯一性与稳定性

作者:

作者单位:

刊名:

英文刊名:

年,卷(期):徐立峰, XU Lifeng湖北师范学院数学与统计学院,湖北,黄石,435002应用数学MATHEMATICA APPLICATA2010,23(4)

参考文献(10条)

1. 王风雨 扩散半群Lipschitz常数的压缩性 1996(04)

2. HERBST I;DITT L Diffusion equation techniques in stochastic monotonicity and positive correbtions1991(87)

3. 徐侃;张绍义 Markov耦合与Markov过程的遍历性[期刊论文]-数学杂志 2001(03)

4. 徐侃;徐立峰 一般状态空间跳过程不变测度的存在性和唯一性[期刊论文]-数学杂志 2004(05)

5. CHEN Mufa;LI Shaofu Coupling methods for multidimensional diffusion processes[外文期刊] 1989(01)

6. 张绍义 耦合方法与随机微分方程平均稳定性 1999(03)

7. 黄志远 随机分析学基础 2001

8. 张绍义 最优可测耦合的存在性与Markov过程的遍历性 1998(11)

9. CHEN Mufa From Markov Chains to Non-equilibrium Particle Systems 1996

10. 张绍义;徐侃 转移概率最优可测耦合的存在性 1997(01)

引证文献(1条)

1. 徐立峰 Markov调制的随机系统平稳分布的存在与唯一性[期刊论文]-数学杂志 2012(1)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_yingysx201004004.aspx

应用数学

MATHEMATICAAPPLICATA

2010,23(4):713—718

扩散过程平稳分布的存在唯一性与稳定性

徐立峰

(湖北师范学院数学与统计学院,湖北黄石435002)

摘要:用耦合方法研究扩散过程稳定性.利用KRW概率距离的对偶表示,获得平稳

分布的存在唯一性、遍历性和依分布稳定性并且给出了收敛速度估计。

关键词:扩散过程;随机微分方程;平稳分布;依分布稳定;耦合

中图分类号:0211.63

文献标识码:AAMS(2000)主题分类:60j60;60H10文章编号:1001—9847(2010)04—0713—06

1.预备知识与引理

设X:是R。上的扩散过程,它满足如下It6型随机微分方程

dX,一6(X。)dt+仃(X,)dB,,(1.1)

式中B,是某概率空间(n,易P)上d维Weiner过程,6(・):R。一R。是d维Borel可测函数,盯(・):R。一R以。是d×dBorel可测函数矩阵.对z∈R。和盯∈R积。定义范数

.Zd2(蚤h|z)l/2H小(i善1%12)1/2・

用P(t,z,dy)表示扩散过程X,的转移概率,简记为P。.(1.1)对应的扩散算子

(卯(z)=pd(z)矗弛)+芝1一i暑d“z)为i弛w∈a(㈣.(1.2)

简记为以“(z),6(z)),其中n(z)=d(z)T(z).

设P・,Pz∈缎R。)(R。上全体概率测度),p∈缄R。×R。),称∥是弘,与岸。的耦合,如果满足边缘性:

P(Al×R。)一卢l(A1)(A1∈∥),卢(R。×A2)一P2(A2)(A2∈∥).

用K(户-,∥z)表示P-与弘z的全体耦合概率,不致混淆时简记为K.对R。上的距离函数p(x,y)一lz—Yl(z,y∈R。)定义P,与户2的KRW概率距离:

w(∥l,户2)=i,n!Ip(x,y)产(dz,dy).

若存在皿∈K(p。,卢2),使

w(pl,胆)一Ip(x,y)fi(dx,dy),

・收稿日期:2009~09—08

基金项目:教育部重点科研项目(209078)作者简介:徐立峰.男.汉族,上海人,讲师。研究方向:概率论与数理统计.

714应用数学2010称乒是.£‘。与户。的最优耦合,最优耦合总是存在的,且此结论可推广到户-,卢z为转移概率的情形‘¨.

对丌∈缄掣),r,P。(-):=In(dx)P(t,z,・).若V£>o,7cP。=硝,则称丌是X,的平稳分布.若对Vz∈R。,当t—+。。时P(t,z,dy)弱收敛于某平稳分布丌(ay),称X依分布稳定.‘

KRW概率距离W(卢。,肛2)一般地强于弱收敛意义下的概率距离ct,(p,,pz)(即Levy—Pro—horov距离),事实上我们有W(卢,,p。)≥叫2(∥z,脾)‘21,因此耦合方法成为研究过程依分布稳定的一个可供选择的工具.

任给厂∈C(R。),定义w)=磐岭掣,驴=(r

Ec(妁,w)<oo)-

L(,)叫做,的Lipschitz常数.

对参考点z。∈R4,记

绲=(卢∈苏掣),lp(x,zo)卢(出)<o。).

注意,若Ⅳ∈玩,则VYo∈R4,(1.3)

Ip(z,Y。)卢(dLc)≤I[P(z,zo)+p(xo,yo)]产(dz)

=lp(z,zo)p(dz)+p(xo,Yo)<∞,

因此z。是任意的.

若7r是X,的平稳分布,且对V户∈玩,W(∥,(・),丌(・))一O(t—co),称X,在KRW概率距离意义下是遍历的.本文中简称为遍历.若存在A>0,使V卢∈玩,Ⅳ(胪。(・),丌(・))≤e--知W(p(・),丌(・)),则称X,指数遍历.

引理1.1C23Vpl,p2∈琨,V矿(胁,脾)一sup{lld(胁一re)f(Y)f,厂∈乡,L(D≤1).

引理1.1被称为KRW距离的对偶表示.

对于一般Markov过程平稳分布的存在唯一性及遍历性有如下一般性准则.

引理1.2[妇设x,(f≥o)是距离可测空间(E,lD,o上的Markov过程,P(t,z,dy)是其转移概率,满足

(i)VT>0,存在0≤Cl(T),c2(T)<∞,使

IP(t,z,dy)p(y,zo)≤c1(丁)+cz(T)lD(z,zo),£∈[o,T3,

其中z。是E中任意固定的一个参考点;(1-4)

(ii)对t∈[o,T3,存在0≤c(£)<1,使

w(P(t,zl,・),P(t,X2,・))≤c(t)p(xl,z2),zl,X2∈E.

则X。存在唯一平稳分布7r∈玩,且X。是遍历的.

本文中我们将用耦合方法着重讨论(1.1)依分布稳定性,不同的是我们采用了对偶原理,它的好处是将研究对象从测度与算子转化为普通实函数,相对简单.本文利用扩散半群P:对Lipschitz常数的压缩性刻划了x,的极限性质,同时得到X。平稳分布的存在唯一性、遍历性及依分布稳定性,并且给出了收敛速度估计,而同类结果中有些需要预先假定平稳分布的存在性,定理条件也弱于某些已有的研究.(1.5)

第4期徐立峰:扩散过程平稳分布的存在唯一性与稳定性7152.主要结果与应用

定理2.1若方程(1.1)的系数满足线性增长条件

6(z)I2+I盯(z)I2≤K(1+l

且对Vt>0,存在0≤r(£)<1,使

L(P。厂)≤,.(£)L(,),厂∈g。,(2.2)z2),z∈RJ,(2.1)

(i)扩散过程X。有唯一的平稳分布丌∈绲,且X,遍历,因而也依分布稳定;

。(ii)V∥∈甥,

W(zP(£,z,・),丌(・))≤r(t)W(F(・),7r(.)).

证首先证明(i),这只需验证引理1.2的两个条件成立.

广(2.3)定义由方程(1.1)确定的期望算子

EⅡ,(x,)]一Jf(y)p(t,z,dy)=P,厂(z).

由T(1.1)的系数满足线性增长条件(2.1),依文[4]引理21.2有矩估计:对VT>0,存在c(丁)>o,使

EzX≤,L1+Z\,酽nzJ∈R£∈[O明,L24

由TJenSen不等式

E21Ⅸrl≤LPI叉r

≤[(1+fz2J“2≤L(1+Iz2)e“7’]“2(2.5)f)2・e。‘D]172一(1+jz})・ec(n/2.

在(1.4)中取参考点z。一0,则

j.P(£,z,dy)P(y,o)一P妒(z,o)一Ez

即此(1.4),再证(1.5).

W(P(t,zl,・),P(t,z2,・))x。I≤(1+lz1)ef‘n/2:一c1(丁)+c2(T)J0(2,。).由引理1.1,并注意到五是Feller过程(以(・)表示z处的点测度)。

一sup{lJ[PQ,zz,dy)一P(t,锄,d了)]厂(y)I,,∈箩,L(厂)≤1)

一sup{Ij以。Pt(d3,)厂(y)一以:Pr(dy)厂(y’I,厂∈驴,L(厂)≤1}

一sup{lJ.乱。cdz,[JIP(t,x,dy),cy,]~.f%cdz,P(t,x,dy)厂cy,]I,厂∈箩,Lc厂,≤-)2sup{lj良,(dz)Pr厂(z)一j以。(dx)Pr,(z)I,,∈g。,L(厂)≤1)

一sup{fJ文。(dx)Pr,(z)一j以。(dy)P。厂(y)f,fE箩,L(厂)≤l}.

记卢是艿1与屯的联合分布,则五(刚,A)=以,(A),互(A,R。)=如(A),于是

上式=sup

≤sup”厂(嘶(dx,W)--P,,(y)压(Ru,dy)f,,∈箩,L(p≤1},JP,f(x)--P,厂(了)f卢(如,dy),厂∈妙,L(厂)≤1}

supIfL(Ⅳ)p(x,y)fi(dx,dy),厂∈矽,L(厂)≤1}

sup{fr(£)L(厂)ID(z,y)fi(dx,dy),fE驴,L(厂)≤ll

716应用数学20lO

≤r(£)Ip(x,y)犀(dx,dy)=r(t)p(xl,z2).(2.6)

由引理1.2,X。有唯一平稳分布丌,且对V卢∈绲,w(∥。(・),zr(・))一o(当t一∞),特别地,取肛=乱,则

W(&P,(・),7f(・))一W(P(£,X,・),丌(・))一0

由于按KRW距离收敛强于弱收敛,故X。依分布稳定.

再证(2.3)式,由X,对丌的不变性,

W(肛P(£,z,・),zr(・))一w(£一∞),(f,工,・),护(z,o,・))

=Sup…矿。江,dy),(Y)一J护(∽,dy)厂(Y)l,,∈驴,L(,)≤l}

一p㈣卜(dz)P(t,z,dy)p一『[』rr(dz)P(t,z,dy)p|',∈驴郴≤・)一suP{lI户(dz)P。,(z)~l7f(dz)P。厂(z)l,,∈g’,L(厂)≤11.1J√

设P(・,・)为岸(・)与丌(・)的最优耦合,与(2.6)类似地可证

w(pP(£,z,・),丌(・))≤,.(f)Ip(z,了)P(dz,dy)一r(t)W(ff,丌).

注2.1设X,为Wiener过程,其转移概率为P(t,士,dy),则线性增长条件(2.1)显然成立(6—0,d==l,rx。),且由文[2]知L(只,)≤L(厂),但熟知Wiener过程没有平稳分布,这说明定理条件(2.2)中的r(£)<1不可放宽为,.(£)≤1,即条件(2.2)是精确的.

注2.2(2.3)反映了扩散半群P,关于Lipschitz常数压缩性与关于KRW距离压缩性的密切联系,由于测度的弱收敛距离锄(卢-,P2)≤W172硼(户1,P2),且w(以(・),丌(・))=Ip(x,y)丌(dy)<f)o(丌(・)∈绲).因此(2.3)也给出了弱收敛速度的估计r172(f)(当r(£)一o).特别地当r(£)一e咄(A>o)时,X,指数遍历,因而也依分布指数稳定.

记妖a(x,y),b(x,y))是颤n(z),6(z))与自身的耦合算子,即

以为妒2,、k,≯“0))’/“(z),c(x,y)、

其中““y)一(c一(而y”以。使矩阵d。,y)非负定,灰z,y)一(::i;).

推论2.1若定理2.1中(2.1)式成立,且存在g与自身的耦合算子勿和常数口>0,使

埤(z,y)≤~qD(z,y),z,Y∈R。,

则X,有唯一平稳分布,r(・)且指数遍历,因而也依分布指数稳定.

可推出(2.7)我们只需验证耦合算子压缩条件(2.7)强于(2.2).事实上根据文[73定理2.3,由(2.7)

E“’ID(X,,Y,)≤e--'p(z,y),

其中yf是由Y出发的解过程.记T—inf{t,X,-=Y,).

P。厂(z)~P,I厂(y)

z~Y(2.8):!.芝[』!墨!!]二旦[丛茎业一!二卫!羔上二f!∑!!!Iz—Y}z—yIT>t]

姑’yJ一』鬻趔-I:Y,匕IX“Yu・1鲁掣YI‰z~l[TM。

≤w炉j-I等并bd≤e气(n

第4期徐立峰:扩散过程平稳分布的存在唯一性与稳定性717于是L(P,厂)≤e1L(厂),因口>0,由定理2.1,X。有唯一平稳分布丌,且由(2.3)X,指数遍历,特别地,(2.3)中取/z一以,则Ⅳ(P(t,X,・),丌(・))≤e1W(艿:(・),丌(・)),即X,依分布指数稳定,指数阶为口/2.

推论2.1明显地优于文[5]中的结果,后者需事先假定平稳分布存在,也优于文f7]中定理2,这里证明了指数遍历性,因而也给出了收敛速度估计.而文[7]在相同条件下,仅得到遍历性.

注2.3在注2.1中我们指出条件L(P,厂)≤L(厂)一般地不能得到平稳分布的存在性,但如果我们假设平稳分布7l"存在唯一,类似于(2.3)的证明仍然可以得到W(P(t,z,・),丌(・))≤w(以(・),丌(・)),此即文[5]中的平均稳定性.

由于Lipschitz常数的压缩性在许多领域有重要应用,因此也为定理2.1的应用提供了一些现成的结果,Herbst和Pitt在研究扩散半群的单调性时给出了如下结果.

引理2.1【83如果b,(z)∈C1(R。)且Ⅱ。(z)仅依赖于zi和zJ(1≤i,J≤d),则

L(P,厂)≤eKtL(,),f∈驴,t≥0,

其中

K翟p罡慧{掣+萎I掣h

考虑如下Hopfield随机神经网络

dX,一[一AX,+Bg(X。)]出+盯(X。)dB,,(2.9)

其中A—diag(aj,d2,…,(1d),c£i>0,1≤i≤d,B一(%)以d,g(z)=(gl(A),92(zz),…,gJ(zd))T,仃(z)一(盯i(z))‘,×d.假定(2.9)有唯一解.

定理2.2在方程(2.9)中,若gi(z)可导,盯(z)=diag(a1(z1),…,O'd(zd)),gf(z),吼(z)满足线性增长条件.

K=sup…ma…x(一ai+b.g,7(zi)+∑J1≤i≤d、:Z』bjigi7(zi)I),J

则当K<0时,(2.9)有唯一平稳分布且依分布指数稳定.

证这足引理2.1与定理2.1的直接结果,且此时c(£)一e触,故指数阶为一K/2.

也可通过构造耦合计算(2.2),利用文[9]中构造的反射耦合我们有

H--=-s。≠up,南{l

论成立.定理2.3在随机神经网络(2.9)中,令d(z)一盯(y)I2一L鱼鱼堕-_鼍望笔掣+<一A(x--y),X--y>+(Big(x)~g(y)],z—y>},则当H<0时,(2.9)有唯一平稳分布且依分布指数稳定(其中(-,・>表示R。上的内积).由文[9]定理1并经简单计算,此时有L(P,厂)≤eHtL(厂),f∈2。,故由定理2.1结

例2.1考虑三维随机神经网络

rdX-(£)=[一X1(£)一0.49l(X1(£))+0.592(Xz(£))+0.69。(K(£))]出+2.4Xl(t)dBl(£),Jd%(£)一[一1.1x2(£)+1.19l(x1(£))一0.692(x2(£))+0.29。(x3(£))]出+2.5托(t)dB2(£),Idx。(f)=[一0.8x。(£)+0.39l(x1(£))+0.792(X2(f))一0.59s(x。(£))]出+2.4X。(t)clB。(£),其中Ⅱ向应函数g。(“)=薹{芝三,i=l,2,3.由于此时o<g7(“)=1一eu_e-U)2<1.易知

718应用数学20lO方程有唯一解,且不难得到K=max(一1+0.7,一1.1+0.7,一0.8+0.5)一一0.3.由定理2.2,X(£)依分布指数稳定,指数阶为0.15.

例2.2考虑二维随机神经网络

,dXl(£)=[一0.9Xl(£)+1.29(X2(£))]d£+0.2XI(t)dWl(£),

其中响应函数g(“)满足Lipschitz条件,如g(摊)=L型_』掣或g(“)一雨eu_e-u.1dx2(£)一[一1.1x2(f)+0.49(Xl(t))]dt+0.25X2(£)dW2(£).

通过计算不难得到上界估计:H<0.252—0.22/2—0.9+1.6/2一一0.0575<0,由定理2.3,X(£)依分布指数稳定,指数阶为o.028.

参考文献:

[13张绍义,徐侃.转移概率最优可测耦合的存在性[J].数学学报,1997,40(1):5-13.

[z]CHENMufa.FromMarkov

ic,1996.ChainstONon-equilibriumParticleSystems[M].Singapore:WorldSeientif-

[3]张绍义.最优可测耦合的存在性与Markov过程的遍历性[J].中国科学(A辑),1998,28(11):999—1008.[4]黄志远.随机分析学基础[M].2版.jE京:科学出版社,2001.

[5]张绍义.耦合方法与随机微分方程平均稳定性[J].北京师范大学学报:自然科学版。1999,35(3):318—

322.

[6]CHENMufa,LIShaofu.Couplingmethodsformultidimensionaldiffusionprocesses[J].Ann.Prob.,

1989。17(1):151-177.

[7]徐侃,张绍义.Markov耦合与Markov过程的遍历性[J].数学杂志,2001,21(3):315—318.

[8]HERBSTI。DITTL.Diffusionequationtechniquesinstochasticmonotonicityandpositivecorrebtions

TheoryRelatFielbs,1991(87):275—312.[J].Probability

[9]王风雨.扩散半群Lipschitz常数的压缩性[J].北京师范大学学报:自然科学版,1996,32(4):447—451.[10]徐侃,徐立峰.一般状态空间跳过程不变测度的存在性和唯一性[J].数学杂志。2004,24(5):561—564.

ExistenceandUniquenessoftheStationaryDistributionand

StabilityinDistributionforSomeDiffusionProcesses

XL厂Lifeng

(College

na)ofMathematicsandStatistics,HubeiNormalUniversity,Huangshi435002.Chi—

Abstract:Couplingmethodisused

thedualityexpressionoftostudythestabilityofsomediffusionprocesses.ByofstationaryKRW—metricweprovetheexistenceanduniqueness

ratedistribution,andtheestimationofconvergenceisobtainedtoo.

Keywords:Diffusionprocesses;Stochasticdifferentialequation;Stationarydistribution;Stabilityindistribution;Coupling

扩散过程平稳分布的存在唯一性与稳定性

作者:

作者单位:

刊名:

英文刊名:

年,卷(期):徐立峰, XU Lifeng湖北师范学院数学与统计学院,湖北,黄石,435002应用数学MATHEMATICA APPLICATA2010,23(4)

参考文献(10条)

1. 王风雨 扩散半群Lipschitz常数的压缩性 1996(04)

2. HERBST I;DITT L Diffusion equation techniques in stochastic monotonicity and positive correbtions1991(87)

3. 徐侃;张绍义 Markov耦合与Markov过程的遍历性[期刊论文]-数学杂志 2001(03)

4. 徐侃;徐立峰 一般状态空间跳过程不变测度的存在性和唯一性[期刊论文]-数学杂志 2004(05)

5. CHEN Mufa;LI Shaofu Coupling methods for multidimensional diffusion processes[外文期刊] 1989(01)

6. 张绍义 耦合方法与随机微分方程平均稳定性 1999(03)

7. 黄志远 随机分析学基础 2001

8. 张绍义 最优可测耦合的存在性与Markov过程的遍历性 1998(11)

9. CHEN Mufa From Markov Chains to Non-equilibrium Particle Systems 1996

10. 张绍义;徐侃 转移概率最优可测耦合的存在性 1997(01)

引证文献(1条)

1. 徐立峰 Markov调制的随机系统平稳分布的存在与唯一性[期刊论文]-数学杂志 2012(1)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_yingysx201004004.aspx


相关内容

  • 大气污染物扩散模式的应用研究综述
  • 环境污染与防治 第29卷 第5期 2007年5 月 大气污染物扩散模式的应用研究综述3 迟妍妍1 张惠远2# (1.中国环境科学研究院,北京100012;2.国家环保总局环境规划院,北京100012) 摘要应用大气污染物扩散模式可以模拟不同尺度.气象.地形条件下工业污染物在大气中的输送与扩散特征,为 ...

  • 基于MCMC方法的贝叶斯统计推断
  • 基础及前沿研究 Fundamentalandfrontierresearch 中国科技信息2012年第10期 CHINASCIENCE AND TEC帆OGY INFORMATIONM_y.2012 DOI:10.3969/j.issn.1001-8972.2012.10.024 基于MCMC方法的 ...

  • 大气边界层湍流复杂系统的观测试验分析与层次相似理论的探讨
  • 湍流是流体在特定条件下所表现出的一种特殊的运动现象,是一种具有强烈涡旋性的不规则运动.湍流问题的研究具有十分重要的科学意义,湍流尤其是大气边界层物理研究的核心问题,因为大气边界层内的流体运动几乎总是处于湍流状态.与一般的流体力学的湍流相比,大气边界层湍流是一个受到更多内外因素影响的非常复杂的多尺度时 ...

  • 晶体生长的有关问题
  • 晶体生长的有关问题 关于退火与热应力及位错 向大家请教两个问题: 1.晶体硅在生长完成后,退火时间越长,热应力释 晶体位错:位错是晶体中局部滑移区域的边界线,即是晶体中的一种线缺陷:它是决定属等晶体力学性质的基本因素,也对晶体的其他许多性质(包括晶体生长)有着严重的影响.通过化学腐蚀可在晶体表面上观 ...

  • 锂离子电池正极材料LiFePO4的结构和电化学反应机理
  • 锂离子电池正极材料LiFePO 4的结构和电化学反应机理 王连亮1, 2 马培华1 李法强1 诸葛芹1 ( 中国科学院青海盐湖研究所 西宁 810008: 中国科学院研究生院 北京 100039) 12 摘 要 十年来的研究并没有对LiFePO 4的电化学反应机理形成准确一致的认识.复合阴离子(PO ...

  • 人文地理学复习资料
  • 第一节 人文地理学的研究对象和特性 一 .地理学与人文地理学 1.当代地理学:研究地球表层地理环境的结构.分布及其发展变化的规律性以及人地关系的科学 2.人文地理学定义:是研究地球表层人类活动或人与地理环境相互关系所形成的现象分布和变化规律的科学 3.人文地理学的研究对象:是关于人类活动的空间差异和 ...

  • 第二节 细胞膜
  • 教学内容 第三章 生命的基本单位--细胞 第二节 细胞膜 一.概念 细胞膜(cell membrane) 又称细胞质膜(plasmamembrane),是指围绕在细胞最外层.将原生质与周围环境分隔开,由膜脂和膜蛋白组成的一层极薄的膜,它是生物膜的重要组成部分. 在光学显微镜下,细胞膜只不过是细胞与周 ...

  • 中南大学生物化学知识点
  • 绪论 1.掌握生物化学的概念和基本研究内容; 概念:应用化学的理论和方法来研究生命现象,在分子水平上解释和阐明生命现象的化学本质 生物化学的基本研究内容: 生命体的分子基础 生物分子的化学组成.结构.性质: 生物分子的结构.功能与生命现象的关系, 如:遗传信息的储存.传递和表达. 生物机体中生物分子 ...

  • 高中生物必修1笔记汇集
  • 第一章走近细胞 第一节从生物圈到细胞 知识梳理: 1病毒没有细胞结构,但必须依赖(活细胞)才能生存. 2生命活动离不开细胞,细胞是生物体结构和功能的(基本单位). 3生命系统的结构层次:(细胞).(组织).(器官).(系统).(个体).(种群)(群落).(生态系统).(生物圈). 4 血液属于(组织 ...