王焕许， 张博洋，柴　畅，王　烁

(绥化学院 信息工程学院，黑龙江 绥化 152061)



线性分段连续型随机延迟微分方程Euler方法的收敛性

王焕许， 张博洋，柴畅，王烁

(绥化学院 信息工程学院，黑龙江 绥化 152061)

摘要：研究线性分段连续型随机延迟微分方程的数值解的收敛性，采用的是Euler方法，在处理线性项的矩阵时，证明的方法主要应用了矩阵欧几里得范数，从而达到要研究线性分段连续型随机延迟微分方程数值解的收敛性的目的，这也是本文解决问题的关键。

关键词：随机延迟微分方程；Euler方法；收敛性；数值解；矩阵欧几里得范数；

1引言与预备

考虑n维的线性分段连续型随机延迟微分方程

(1)

其中B(t)是一个布朗运动过程。

由随机微分的定义， 此方程等价于下面的积分方程

(2)

yn+1=yn+(Ayn+Byh([nh]))h+(Cyn+Dyh([nh]))ΔBn

(3)

其中 n=0,1,2,…,并且ΔBn=B(tn)-B(tn-1)，yn是x(tn)的近似解，而yh([nh])是x([nh])的近似解，为了计算简便易懂，令n=km+l(k=0,1,2,…,l=0,1,2,…,m-1)，因此(3)可以简化为下面的格式，

ykm+l+1=ykm+l+(Aykm+l+Bykm)h+(Cykm+l+Dykm)ΔBkm+l

(4)

其中ΔBkm+l=B(tkm+l)-B(tkm+l-1)，ykm+l是x(tkm+l)的近似解，而ykm是x([tkm+l])的近似解。把(4)进行连续化得到下面的式子，

(5)

其中对于t∈[tkm+l,tkm+l+1),z(t)=ykm+l,z([t])=ykm.

2Euler方法的收敛性

下面给出线性分段连续型随机延迟微分方程数值解的收敛性。首先给出两个引理。

引理1存在一个正数C1使得方程(1.2) 的解和连续的欧拉方法的数值解(5)满足

(6)

其中C1=C1(T,L2)是与h无关的常数。

证明 由(5)可得

由H¨older 不等式可得

对于所有0≤t≤T，得到

通过Doob's 鞅不等式, 得到

因此，有

由Gronwall不等式得到

其中C1=3E|x0|2exp{6(T|A|2+T|B|2+4|C|2+4|D|2)}.

同理可得

其中C1=3E|x0|2exp{6(T|A|2+T|B|2+4|C|2+4|D|2)}. 证毕。

引理2 存在一个正数C2使得下面的式子成立

(7)

其中C2=C2(T,L2)是与h 无关的常数。

证明 对于t∈[0,T]，则存在两个常数k,l使得t∈[tkm+l,tkm+l+1)，那么

y(t)=ykm+l+(Aykm+l+Bykm)(t-tkm+l)+(Cykm+l+Dykm)(B(t)-B(tkm+l))

因此，

|y(t)-z(t)|2=|ykm+l+(Aykm+l+Bykm)(t-tkm+l)

+(Cykm+l+Dykm)(B(t)-B(tkm+l))-ykm+l|2

≤|(Aykm+l+Bykm)(t-tkm+l)+(Cykm+l+Dykm)(B(t)-B(tkm+l))|2

对上面的式子取期望得到……