Advances in Applied Mathematics
Vol. 13  No. 02 ( 2024 ), Article ID: 82270 , 18 pages
10.12677/AAM.2024.132077

双尺度随机时滞微分方程的平均原理

贺鑫

长安大学理学院,陕西 西安

收稿日期:2024年1月28日;录用日期:2024年2月22日;发布日期:2024年2月29日

摘要

本文研究了分数布朗运动驱动的非自治双尺度随机时滞微分方程的平均原理。首先,通过广义Stieltjes积分和随机平均原理,推导了非自治双尺度系统的均方收敛定理。然后,结合均方收敛定理和停时理论,分别得到了原系统和平均系统的矩估计。最后,证明了当时间尺度参数趋于零时,慢变量方程的解过程在均方意义下收敛于平均方程的解过程。

关键词

双尺度,随机时滞微分方程,平均原理,分数布朗运动

The Averaging Principle of Two-Scale Stochastic Delay Differential Equation

Xin He

School of Sciences, Chang’an University, Xi’an Shaanxi

Received: Jan. 28th, 2024; accepted: Feb. 22nd, 2024; published: Feb. 29th, 2024

ABSTRACT

The main goal of this article is to study an average principle of a class of non-autonomous two time-scale stochastic differential delay equations driven by fractional Brownian motion. Firstly, the mean square convergence theorem for non-autonomous scale systems was derived by means of Generalized Stieltjes integral and Stochastic Average Principle. Then, combining the mean square convergence theorem and the Stopping-time theory, the moment estimates of the original system and the average system were obtained, respectively. Finally, it showed that when the time scale parameters approach zero, the solution process of the slow variable equation converges to the solution process of the mean equation in the mean square sense.

Keywords:Two-Time-Scale, Stochastic Differential Delay Equations, Averaging Principle, Fractional Brownian Motion

Copyright © 2024 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

双尺度随机时滞微分方程最早可以追溯到Itô和Gikhman创立随机微积分理论的年代。由于在化学 [1] 、力学 [2] 、电子电路 [3] 和大气海洋学 [4] 等实际应用领域中表现出来的广泛应用,双尺度随机微分方程俨然已经成为一个重要的研究领域,越来越多的学者开始重视双尺度随机微分方程的研究。对于多尺度随机微分方程的定性分析,平均原理提供了极大帮助。多尺度随机微分方程的平均原理就是在一定的条件下,将快变量过程视为随机噪声并将其消除从而得到一个极限过程,使得极限过程的解逼近原始方程中慢变量的解。平均原理广泛应用于力学、物理、控制等领域,是简化动力系统从而得到微分方程近似解的有力工具。

关于随机动力系统的平均原理有较长的历史,其奠基性工作由前苏联数学家Bogoliubov在文献 [5] 中完成。紧接着,Gikhman [6] ,Volosov [7] 和Besjes [8] 研究了非线性常微分方程的平均问题。随机平均原理首先由Stratonovich提出,此后,Khasminskii [9] 将平均原理发展到具有快–慢尺度的随机常微分方程的研究中,证明了在较弱的收敛意义下成立的平均原理。值得一提的是,Freidlin和Wentzell [10] [11] 给出了依概率收敛的平均原理。另外,Golec和Ladde [12] 以及Givon [13] 等人又进一步将结果推广到均方意义下收敛的情形。然后,Cerrai和Freidlin [14] 研究了一类无限维随机反应扩散模型的平均原理。近年来,Fu和Duan [15] 深入研究了带有高斯噪声的双尺度随机偏微分方程,得到在Lipschitz条件下的均方收敛的平均结果。随后,Fu和Liu [16] 又进一步将结果推广为:慢方程的解强收敛于平均方程的解。Xu [17] 将以布朗运动为驱动的随机系统推广到非高斯噪声(带跳)的情形并得到了强收敛结果。之后,Xu [18] 等人又进一步证得由分数布朗运动驱动的双尺度随机微分方程的平均原理。

然而,许多研究都是基于系统未来的状态与过去的状态是相互独立的假设下展开的,这与实际情况并不相符。事实上,滞后现象不可避免地出现在随机动力系统中,也就是说事物的发展趋势既依赖于当前的状态,还依赖于过去的历史状态。时滞是实际应用中许多系统的重要特征,在生物学、信号传输、随机控制等应用领域的很多数学模型中,都发挥着不可或缺的作用,许多学者都进行了研究 [19] [20] [21] 。文献 [19] 给出了随机时滞微分方程和带跳的随机时滞微分方程的平均原理,证明了平均方程的解在p阶矩意义下和依概率意义下收敛于随机微分方程的解。文献 [20] 讨论了基于Hurst参数 H ( 1 / 2 , 1 ) 的分数布朗运动驱动的随机时滞微分方程的平均原理,原始方程的解依概率收敛和均方意义收敛于相应平均方程的解。

双尺度随机微分方程的平均原理已经取得了很大的进展,但对于双尺度随机时滞微分方程的平均法的研究却很少。文献 [22] 利用随机平均原理研究了由布朗运动驱动的自治双尺度随机时滞微分方程,为一类自治双尺度随机时滞微分方程建立了一个平均原理的强极限定理。文献 [23] 则在文献 [22] 的基础上研究了Hurst参数 H ( 1 / 2 , 1 ) 的乘法分数布朗噪声驱动的自治双时间尺度随机时滞微分方程的平均原理。

受此启发,本文在文献 [20] [23] 的研究基础上,研究由分数布朗运动驱动的非自治双尺度随机时滞微分方程的平均原理,其中该非自治系统的慢变量过程由Hurst参数 H ( 1 / 2 , 1 ) 的分数布朗运动驱动,快变量过程由布朗运动驱动。借助广义Stieltjes积分和停时理论,通过随机平均原理得到非自治双尺度随机时滞微分方程的平均方程,并证得平均方程和原方程的慢分量在均方意义下解的近似等价性。

本文的其余部分组织如下:第2节预备知识。第3节是本文的主要结果及其证明过程,即证明慢方程的解是强收敛于平均方程的解。第4节总结了本文。

2. 预备知识

假设是一个带流的完备概率空间,其中 Ω = [ 0 , T ] 为实轴上的有限区间,T是正常数, 是概率测度,满足通常条件,即是右连续的且包含所有的 零测集。符号 | · | 表示 n 中的Euclidean范数, · 表示矩阵范数。 : = C ( [ τ , 0 ] ; n ) 表示所有从 [ τ , 0 ] 映到 n 的右连续且具有左极限的函数f构成的Banach空间,其范数为最大范数 · 。对于 h ( · ) C ( [ τ , ) ) t 0 ,定义一个随机过程 h : = { h ( t ) , t [ τ , ) } h t 定义为 h t ( θ ) = h ( t + θ ) θ [ τ , 0 ] W : = { W t , t [ 0 , T ] } 是定义在完备概率空间上的m维布朗运动, B H : = { B t H , t [ 0 , T ] } ( 1 / 2 < H < 1 ) 是一个定义在完备概率空间上的n维分数布朗运动。在本文中,符号C代表常数,其值在变化。

定义2.1 [23] 令 a , b 。阶数为 α ( α > 0 ) 的左侧和右侧Riemann-Liouville分数阶积分被定义为

I α + α f ( t ) = 1 Γ ( α ) a t ( t ω ) α 1 f ( ω ) d ω , I b α f ( t ) = ( 1 ) α Γ ( α ) t b ( ω x ) α 1 f ( ω ) d ω ,

其中 ( 1 ) α = e i π α Γ ( α ) = 0 r α 1 e r d r 是Euler Gamma函数。

假设 1 / 2 < H < 1 1 H < α < 1 ,记 W α , 1 ( [ 0, T ] ; n ) 表示由可测函数 f : [ 0 , T ] n ,且

f α , 1 : = 0 T s α | f ( s ) | d s + 0 T 0 s ( s ω ) α 1 | f ( s ) f ( ω ) | d ω d s < (1)

构成的空间。假设 1 / 2 < H < 1 1 H < α < 1 ,记 W 0 α , ( [ 0, T ] ; n ) 表示由可测函数 f : [ 0 , T ] n ,且

f α , : = sup t [ 0 , ] ( | f ( t ) | + s t ( r ω ) α 1 | f ( t ) f ( ω ) | d ω ) < (2)

构成的空间。对任意 λ 0 定义等价范数

f α , : = sup t [ 0 , T ] e λ t ( | f ( t ) | + s t ( r ω ) α 1 | f ( t ) f ( ω ) | d ω ) < . (3)

定义2.3 [24] 令 a , b 。阶数为 α ( α > 0 ) 的广义Riemann-Stieltjes分数阶积分被定义为

0 T f ( r ) d g ( r ) = ( 1 ) α 0 T D 0 + α f ( r ) D T 1 α g T ( r ) d r , s t f ( r ) d g ( r ) = 0 T f ( r ) 1 ( s , t ) d g ( r ) .

其中对于任意的 0 a < b T g b ( t ) : = g ( t ) g ( b ) ,且对于任意的 a < t < b ,阶数为 α ( α > 0 ) 的左侧和右侧Weyl型分数阶积分被定义为

D α + α f ( t ) : = 1 Γ ( 1 α ) [ f ( t ) ( t a ) α + α a t f ( t ) f ( s ) ( t s ) α + 1 d s ] ,

D b 1 α g b ( t ) = ( 1 ) 1 α Γ ( 1 α ) [ g ( b ) g ( t ) ( b t ) 1 α + ( 1 α ) t b g ( t ) g ( s ) ( s t ) 2 α d s ] .

其中 Γ ( · ) 是Gamma函数。

假设 1 / 2 < H < 1 1 H < α < 1 ,设 W T α , 0 ( [ 0, T ] ; n ) 是可测函数 g : [ 0 , T ] m 构成的空间,且具有范数

g α , 0 , T : = sup 0 s < t T ( | g ( t ) g ( s ) | ( t s ) 1 α + s t | g ( u ) g ( s ) | ( u s ) 2 α d u ) < .

引理2.4 [25] 假设 H ( 1 / 2 , 1 ) 且对任意的 α ( 1 H , 1 / 2 ) ,函数 f W α , 1 ( [ 0, T ] ; n ) g W 1 α , ( [ 0 , T ] ; n ) ,且对于任意 t [ 0 , T ] ,积分 0 t f d B s H 存在,则有以下不等式成立:

| 0 t f d B s H | Λ α ( B s H ) f α , 1 .

其中 Λ α ( B H ) : = 1 Γ ( 1 α ) Γ ( α ) B H α , 0 , T

进一步,根据Fernique定理,对于任何 0 < ϑ < 2 ,有 E [ e ( Λ α ( B H ) ) ϑ ] < (4)

假设 τ > 0 s , t [ τ , T ] ,记 W 0 α , ( [ s , t ] ; n ) 表示可测函数 f : [ s , t ] n 构成的空间且满足以下可积条件

f α , ( s , t ) : = sup r [ s , t ] ( | f ( r ) | + s r ( r u ) α 1 | f ( r ) f ( u ) | d u ) < . (5)

这里 f : = sup r [ s , t ] e λ r | f ( r ) | f ( r ) α ( s ) : = | f ( r ) | + s r ( r u ) α 1 | f ( r ) f ( u ) | d u

对于任意 λ 0 定义等价范数 f α , λ ( s , t ) : = sup r [ s , t ] e λ r ( | f ( r ) | + s r ( r u ) α 1 | f ( r ) f ( u ) | d u ) < (6)

对于任意 0 < λ 1 ,记 C λ ( [ s , t ] ; n ) 表示由λ-阶Hölder连续函数 f : [ s , t ] n

f α , λ ( s , t ) p : = sup r [ s , t ] e λ r ( | f ( r ) | + s r ( r u ) α 1 | f ( r ) f ( u ) | d u ) p , p 1. (7)

构成的空间。这里 f α , ( τ ) p : = f α , ( τ , T ) p f α , λ ( τ ) p : = f α , λ ( τ , T ) p f α ( τ ) p : = f α , ( τ ) p 。当 τ = 0 时,我们将省略对应范数的τ。

引理2.5 [26] 对于任意非负的a和b,且 a + b < 1 ,并且对于任意的 λ 1 ,存在正常数C,使得

0 t e λ ( t r ) ( t r ) a r b d r C λ a + b 1 .

0 a < 1 b 0 以及任意的 λ 1 时,有 0 t r b ( t r ) a r b d r Γ ( 1 a ) t b λ a 1

引理2.6 [25] 对于可测函数 f : [ 0 , T ] n ,存在正常数C,使得

0 t f ( r ) d r α C 0 t | f ( r ) | ( t r ) α d r ,

0 t f ( r ) d B r H α C Λ α ( B H ) 0 t ( ( t r ) 2 α + r α ) ( | f ( r ) | + 0 t ( r u ) α 1 | f ( r ) f ( u ) | d u ) d r

本文考虑如下形式的双尺度随机时滞微分方程:

{ d X ε ( t ) = f ( t , X t ε , Y t ε ) d t + σ 1 ( t , X ε ( t τ ) ) d B t H , t [ 0 , T ] , X 0 ε = ξ d Y ε ( t ) = 1 ε g ( X t ε , Y ε ( t ) , Y ε ( t τ ) ) d t + 1 ε σ 2 ( X t ε , Y ε ( t ) , Y ε ( t τ ) ) d W t , t [ 0 , T ] , Y 0 ε = η (8)

其中 t [ 0 , T ] τ > 0 为正常数时滞, 0 < ε 1 是一个小的正参数,代表系统中的时间尺度比率。因此,根据这个时间尺度,变量 X t ε 称为慢分量,变量 Y t ε 称为快分量。我们对系数做如下定义: f : [ 0 , T ] × × n σ 1 : [ 0 , T ] × n n × m g : × n × n n σ 2 : × n × n n × m

现在,假设方程(8)的各个系数满足以下条件:

(H1) f是可测函数,对于任意的 x i , y i ( i = 1 , 2 ) s , t [ 0 , T ] ,存在常数 θ i 0 , ( i = 1 , 2 , 3 ) , 0 < κ 1 L i 0 ( i = 1 , 2 ) ,使得

| f ( t , x 1 , y 1 ) | L 1 ( 1 + | x 1 | + | y 1 | ) ; | f ( t , x 1 , y 1 ) f ( t , x 2 , y 1 ) | L 2 ( 1 + | y 1 | θ 1 ) | x 1 x 2 | ; | f ( t , x 1 , y 1 ) f ( t , x 1 , y 2 ) | L 2 | y 1 y 2 | ( 1 + | x 1 | θ 2 + | y 1 | θ 3 + | y 2 | θ 3 ) ; | f ( t , x 1 , y 1 ) f ( s , x 1 , y 1 ) | L 2 | t s | κ ( 1 + | x 1 | θ 2 + | y 1 | θ 3 ) .

(H2) 对任意的 x 1 , x 2 n s , t [ 0 , T ] ,存在常数 β ( 0 , 1 ] L i > 0 , i = 3 , 4 , 5 , 6 使得

| x σ 1 ( t , x 1 ) | L 3 ; | σ 1 ( t , x 1 ) | L 6 ( 1 + | x | γ ) ; | x σ 1 ( t , x 1 ) x σ 1 ( t , x 2 ) | + | σ 1 ( t , x 1 ) σ 1 ( t , x 2 ) | L 4 | x 1 x 2 | ; | x σ 1 ( t , x 1 ) x σ 1 ( s , x 1 ) | + | σ 1 ( t , x 1 ) σ 1 ( s , x 1 ) | L 5 | t s | β .

(H3) 假设 g = ( ( 1 ) g , ( 2 ) g , ( 3 ) g ) σ 2 = ( ( 1 ) σ 2 , ( 2 ) σ 2 , ( 3 ) σ 2 ) 是有界的。

(H4) 对于任意的 x y 1 , y 2 , z 1 , z 2 n ,存在 λ i > 0 , i = 1 , 2 λ 1 > λ 2 > 0 且与x无关,使得

2 y 1 y 2 , g ( x , y 1 , z 1 ) g ( x , y 2 , z 2 ) + σ 2 ( x , y 1 , z 1 ) σ 2 ( x , y 2 , z 2 ) 2 λ 1 | y 1 y 2 | 2 + λ 2 | z 1 z 2 | 2 , 2 y , g ( x , y , z ) + σ 2 ( x , y , z ) λ 3 | y | 2 + λ 4 | z | 2 + C ( 1 + x 2 ) .

(H5) 对于初值 X 0 ε = ξ ,存在 λ 5 > 0 ,使得 | ξ ( t ) ξ ( s ) | λ 5 | t s | s , t [ τ , 0 ]

不难证明,在条件(H1)~(H3)和(H5)下,方程(8)存在唯一解:

{ X ε ( t ) = ξ ( 0 ) + 0 t f ( s , X s ε , Y s ε ) d s + 0 t σ 1 ( s , X ε ( s τ ) ) d B s H , Y ε ( t ) = η ( 0 ) + 1 ε 0 t g ( X s ε , Y ε ( s ) , Y ε ( s τ ) ) d s + 1 ε 0 t σ 2 ( X s ε , Y ε ( s ) , Y ε ( s τ ) ) d W s , X 0 ε = ξ , Y 0 ε = η , t [ 0 , T ] . (9)

接下来,我们推导当 ε 0 + 时,慢分量 X t ε 收敛于以下平均方程(利用随机平均原理,我们可以得到方程(8)的平均方程)的解:

d X ¯ ( t ) = f ¯ ( t , X ¯ t ) d t + σ 1 ( t , X ¯ ( t τ ) ) d B t H , t [ 0 , T ] , X ¯ 0 = ξ . (10)

其中 μ x 是相应冻结方程的转移半群 ( P t x ) t 0 的唯一不变测度。显然,在条件(H1)~(H5)下,方程(10)也有唯一解 ( X ¯ ( t ) ) t τ [27] 。

下面,我们来介绍冻结方程。对于任意 t [ 0 , T ] ,固定慢变量 x

d Y ( t ) = g ( x , Y ( t ) , Y ( t τ ) ) d t + σ 2 ( x , Y ( t ) , Y ( t τ ) ) d W t , Y 0 = η . (11)

在条件(H3)下,方程(11)有唯一的解 ( Y ( t ) ) t τ [22] 。

在本部分最后,给出以下引理,这是证明本文结论的重要工具:

定义2.7 设 F t { Y r χ , η , r t } 生成的 σ -域。对 0 ζ s T ,有

J ( s , ζ , χ , η ) = E [ f ( k δ , χ , Y s χ , η ) f ¯ ( k δ , χ ) , f ( k δ , χ , Y ζ χ , η ) f ¯ ( k δ , χ ) ] , (12)

此时,存在独立于 s , ζ 的常数 C , ρ > 0 ,使得 J ( s , ζ , χ , η ) C ( 1 + χ 2 + η 2 ) e ρ 2 ( s ζ )

证明 由式(12),调用 Y t χ , η 的马尔科夫性质,Hölder不等式、条件(H1)和文献 [22] ,可以证明。

3. 主要结果

在本节,我们给出本文的主要结论:

定理3.1 假设原始方程(8)和平均方程(10)都满足假设条件(H1)~(H5),则慢分量 X t ε 均方收敛于平均方程的解 X ¯ t 。即对任意的 t [ 0 , T ] ,对于所有的 ε ( 0 , 1 ) ,有

lim ε 0 E [ X ε X ¯ α , ( τ ) 2 ] = 0.

证明 定理3.1的证明包括以下步骤:

Step1:我们给出方程(8)的解 ( X ε , Y ε ) 的估计。

首先,我们对快慢过程有如下估计:

定理3.2 假设条件(H1)~(H5)成立,那么对任意 t [ 0 , T ] ,存在正常数C,使得当 ε ( 0 , 1 ) 时,有

E [ X ε α , ( τ ) 2 ] C .

证明 设 Λ : = Λ α ( B H ) 1 ,对任意的 λ 1 ,首先,我们估计 X ε , λ ( τ ) , t ,由条件(H1)~(H2)、引理2.4和引理2.5,有以下估计

X ε , λ ( τ ) , t : = sup s [ τ , t ] e λ s | X ε ( s ) | sup s [ τ , 0 ] e λ s | ξ ( s ) | + sup s [ 0 , t ] e λ s | X ε ( s ) | sup s [ τ , 0 ] e λ s | ξ ( s ) | + sup s [ 0 , t ] e λ s ( | 0 s f ( r , X r ε , Y r ε ) d r | + | 0 s σ 1 ( r , X ε ( r τ ) ) d B r H | ) ξ α , ( τ , 0 ) + C sup 0 s t 0 s e λ ( s r ) ( 1 + X ε , λ ( τ ) , t ) d r + Λ α ( B H ) sup 0 s t ( 0 s e λ ( s u + τ ) [ ( 1 + X ε , λ ( τ ) , t ) u α + X ε 1 , λ ( τ ) , t ] d u ) K Λ ( 1 + λ α 1 X ε , λ ( τ ) , t + λ 1 X ε 1 , λ ( τ ) , t ) . (13)

其中常数 K > 1 且与 ξ α , ( τ , 0 ) 有关。

对于 X ε 1 , λ ( τ ) , t ,有

X ε 1 , λ ( τ ) , t : = sup s [ τ , t ] e λ s τ s ( s r ) α 1 | X ε ( s ) X ε ( r ) | d r sup s [ τ , 0 ] e λ s τ s ( s r ) α 1 | ξ ( s ) ξ ( r ) | d r + sup s [ 0 , t ] e λ s τ 0 ( r ) α 1 | ξ ( 0 ) ξ ( r ) | d r + sup s [ 0 , t ] e λ s τ 0 ( s r ) α 1 | X ε ( r ) ξ ( 0 ) | d r + sup s [ 0 , t ] e λ s 0 s ( s r ) α 1 | X ε ( s ) X ε ( r ) | d r = : i = 1 4 B i . (14)

接下来,我们分别估计 B 1 B 2 B 3 B 4 。对于 B 1 B 2 ,由条件(H5),

B 1 + B 2 ξ α , ( τ , 0 ) . (15)

对于 B 3 ,由引理2.4,可以得到

B 3 sup s [ 0 , t ] e λ s s α | X ε ( s ) ξ ( 0 ) | sup s [ 0 , t ] 1 s α 0 s e λ ( s r ) ( 1 + X ε , λ ( τ ) , t ) d r + Λ α ( B H ) C α , β sup s [ 0 , t ] 1 s α 0 s e λ ( s u + τ ) [ ( 1 + X ε , λ ( τ ) , t ) u α + X ε 1 , λ ( τ ) , t ] d u Λ α ( B H ) C α , β sup s [ 0 , t ] 0 s e λ ( s u + τ ) [ ( 1 + X ε , λ ( τ ) , t ) u 2 α + ( s u ) α X ε 1 , λ ( τ ) , t ] d u K Λ ( 1 + λ 2 α 1 X ε , λ ( τ ) , t + λ α X ε 1 , λ ( τ ) , t ) . (16)

对于 B 4 ,由引理2.6,可得

B 4 sup s [ 0 , t ] e λ s 0 s ( s r ) α 1 r s f ( u , X u ε , Y u ε ) d u d r + sup s [ 0 , t ] e λ s 0 s ( s r ) α 1 | r s σ 1 ( u , X ε ( u τ ) ) d B u H | d r sup s [ 0 , t ] 0 s e λ ( s r ) ( s r ) α ( 1 + X ε , λ ( τ ) , t ) d r + C Λ α ( B H ) sup s [ 0 , t ] ( 0 s e λ ( s u + τ ) ( s u ) 2 α ( 1 + X ε , λ ( τ ) , t ) d u + 0 s e λ ( s u + τ ) ( s u ) α X ε 1 , λ ( τ ) , t d u ) . (17)

最后,将式(15)~(17)代入不等式(14),得到

X ε 1 , λ ( τ ) , t K Λ ( 1 + λ 2 α 1 X ε , λ ( τ ) , t + λ α X ε 1 , λ ( τ ) , t ) . (18)

选取 λ = ( 4 K Λ ) 1 / ( 1 α ) ,由式(13),得到

X ε , λ ( τ ) , t 4 3 K Λ ( 1 + λ 1 X ε 1 , λ ( τ ) , t ) . (19)

将上式代入不等式(18),进行简单变换,得到

X ε 1 , λ ( τ ) , t 3 2 K Λ + 2 ( K Λ ) 1 / ( 1 α ) C T , | ξ | Λ 1 / ( 1 α ) .

再将其代入不等式(19),则 X ε , λ ( τ ) , t C Λ 1 / ( 1 α )

因此,可以得到

X ε α , ( τ ) e λ T ( X ε , λ ( τ ) , T + X ε 1 , λ ( τ ) , T ) C e Λ 1 / ( 1 α ) ( 1 + Λ α 1 / ( 1 α ) ( B H ) ) .

由于 0 < 1 / ( 1 α ) < 2 和引理2.4中的 E [ e ( Λ α ( B H ) ) ϑ ] < ,可以得到

E [ X ε α , ( τ ) 2 ] C .

定理3.2得证。

此外,同理可证

X ^ ε α , ( τ ) + X ¯ α , ( τ ) C e C Λ α 1 / ( 1 α ) ( B H ) ( 1 + Λ α 1 / ( 1 α ) ( B H ) ) . (20)

E [ X ^ ε α , ( τ ) 2 + X ¯ α , ( τ ) 2 ] C . (21)

其次,慢过程有如下连续性:

定理3.3 假设条件(H1)~(H5)成立,那么任意 t [ 0 , T ] ,存在正常数C,使得对足够小的 δ ε ( 0 , 1 ) ,有

sup t [ 0 , T ] E [ X t ε X t ( δ ) ε 2 ] C δ 1 2 α .

证明 首先,我们估计 X ε ( t ) X ε ( s ) 。对于任意的 t , s [ 0 , T ] ,存在正常数C,使得

| X ε ( t ) X ε ( s ) | | s t f ( r , X r ε , Y r ε ) d r | + | s t σ 1 ( r , X ε ( r τ ) ) d B r H | C ( 1 + X ε α , ( τ ) ) | t s | + C Λ α ( B H ) ( s τ t τ ( ( 1 + X ε α , ( τ ) ) ( r s ) α d r + s τ v | X ε ( v ) X ε ( u ) | ( v u ) α 1 d u d v ) ) C ( 1 + Λ α ( B H ) ) ( 1 + X ε α , ( τ ) ) | t s | 1 α .

其次,小区间长度为 δ : = τ N < 1 (N为某一固定的正整数)。对任意的 t [ 0 , T ] θ [ τ , 0 ] ,存在常数 k , m 0 ,若 t [ k δ , ( k + 1 ) δ ] θ [ ( m + 1 ) δ , m δ ] ,则 t + θ [ ( k m 1 ) δ , ( k + 1 m ) δ ] t ( δ ) + θ [ ( k m 1 ) δ , ( k m ) δ ] 。对于 E [ X t ε X t ( δ ) ε 2 ] ,有

E [ X t ε X t ( δ ) ε 2 ] E [ m = 0 N 1 sup θ [ ( m + 1 ) δ , m δ ] | X ε ( t + θ ) X ε ( t ( δ ) + θ ) | 2 ] N max m = 0 , , N 1 E [ sup θ [ ( m + 1 ) δ , m δ ] | X ε ( t + θ ) X ε ( t ( δ ) + θ ) | 2 ] = : N max m = 0 , , N 1 J ( t , m , δ ) .

为了完成定理3.3的证明,我们只需证明 J ( t , m , δ ) C δ 2 ( 1 α )

因此,我们讨论以下三种情况:

情况1:当 m k 1 时,根据Hölder不等式,引理3.1,条件(H1)和(H2),存在正常数C,使得

J ( t , m , δ ) C E [ sup θ [ ( m + 1 ) δ , m δ ] | k δ + θ t + θ f ( r , X r ε , Y r ε ) d r + k δ + θ t + θ σ 1 ( r , X ε ( r τ ) ) d B r H | 2 ] C δ E [ k δ ( m + 1 ) δ t m δ | f ( r , X r ε , Y r ε ) | 2 d r ] + C E [ | k δ ( m + 1 ) δ t ( m + 1 ) δ σ 1 ( r , X ε ( r τ ) ) d B r H | 2 ] + C E [ sup θ [ ( m + 1 ) δ , m δ ] | t ( m + 1 ) δ t + θ σ 1 ( r , X ε ( r τ ) ) d B r H | 2 ] + C E [ sup θ [ ( m + 1 ) δ , m δ ] | k δ ( m + 1 ) δ k δ + θ σ 1 ( r , X ε ( r τ ) ) d B r H | 2 ] C δ 2 + C δ 2 ( 1 α ) E [ ( 1 + Λ α ( B H ) ) 2 ( 1 + X ε α , ( τ ) ) 2 ] C δ 2 ( 1 α ) .

情况2:当 m k 1 时,由条件(H5),存在正常数C,使得

| X ε ( t + θ ) X ε ( t ( δ ) + θ ) | 2 = | ξ ( t + θ ) ξ ( t ( δ ) + θ ) | 2 C δ 2 .

情况3:当 m = k 时,由Hölder不等式,条件(H1)和(H2),存在正常数C,使得

J ( t , m , δ ) = E [ sup θ [ ( m + 1 ) δ , m δ ] | X ε ( t + θ ) X ε ( t δ + θ ) | 2 ] C δ 2 + C E [ sup θ [ ( m + 1 ) δ , m δ ] | X ε ( t + θ ) X ε ( t ( δ ) + θ ) | 2 1 { t + θ > 0 } ] C δ 2 + C E [ sup θ [ t , k δ ] | 0 t + θ f ( r , X r ε , Y r ε ) d r | 2 ] + C E [ sup θ [ t , k δ ] | 0 t + θ σ 1 ( r , X ε ( r τ ) ) d B r H | 2 ] C δ 2 ( 1 α ) .

综上所述, sup t [ 0 , T ] E [ X t ε X t ( δ ) ε 2 ] C δ 1 2 α

定理3.3得证。

参考文献 [22] 中定理3.1和文献 [23] 中引理3.10的证明,易得如下结论:

定理3.4 假设条件(H1)~(H5)成立,那么对任意 t [ 0 , T ] ,存在正常数C,使得

sup t [ 0 , T ] E [ Y t ε 2 ] C .

Step2:根据文献 [9] ,我们先构造辅助过程 ( X ^ t ε , Y ^ t ε ) × ,再对区间 [ 0 , T ] 进行分割,令小区间长度为 δ : = τ N < 1 。对于过程 Y ^ t ε ,令它的初值 Y ^ 0 ε = Y 0 ε = η ,对 t [ k δ , min { ( k + 1 ) δ , T } ] k 0 并满足方程

Y ^ ε ( s ) = Y ^ ε ( k δ ) + 1 ε k δ t g ( X k δ ε , Y ^ ε ( s ) , Y ^ ε ( s τ ) ) d s + 1 ε k δ t σ 2 ( X k δ ε , Y ^ ε ( s ) , Y ^ ε ( s τ ) ) d W s ,

其中 X k δ ε Y k δ ε 分别为在 k δ 时的慢变量和快变量的解过程。令 · 表示取整函数,并定义 t [ 0 , T ] 时,对于过程 X ^ t ε 满足方程

X ^ ε ( t ) = ξ + 0 t f ( s ( δ ) , X s ( δ ) ε , Y ^ s ε ) d s + 0 t σ 1 ( s , X ε ( s τ ) ) d B s H . (22)

其中 s ( δ ) = s δ δ 表示处于点s右端并且距离其最近的分点。

现在,我们来估计 X t ε X ^ t ε Y t ε Y ^ t ε 这两项。

参考文献 [22] 中定理4.2和文献 [23] 中引理3.9的证明,易得如下结论:

定理3.5 假设(H1)~(H5)成立,那么对任意 t [ 0 , τ ) ,存在正常数C和独立于 ε β > 0 ,使得对足够小的 δ ε ( 0 , 1 ) ,有

sup t [ 0 , T ] E [ Y t ε Y ^ t ε 2 ] C ε 1 δ 1 2 α e β δ ε .

接下来,我们估计 X ^ ε X ε ,得到以下引理:

定理3.6 假设条件(H1)~(H5)成立,那么对任意 t [ 0 , T ] ,存在正常数C和独立于 ε β > 0 ,使得对足够小的 δ ε ( 0 , 1 ) ,有

E [ X ^ ε X ε α , λ ( τ ) 2 ] C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) .

证明 根据方程(8)和方程(22),

E [ X ^ ε X ε α , λ ( τ ) 2 ] C E [ sup t [ 0 , T ] e λ t | X ^ ε ( t ) X ε ( t ) | 2 ] + C E [ sup t [ 0 , T ] e λ t ( 0 t ( t s ) α 1 | X ^ ε ( t ) X ε ( t ) X ^ ε ( s ) + X ε ( s ) | d s ) 2 ] + C E [ sup t [ 0 , T ] e λ t ( τ 0 ( t s ) α 1 | X ^ ε ( t ) X ε ( t ) X ^ ε ( s ) + X ε ( s ) | d s ) 2 ] = : D 1 + D 2 + D 3 .

首先,对于第一项 D 1 和第二项 D 2 ,由条件(H1)~(H2),引理2.6,定理3.3和定理3.5可得

D 1 + D 2 C E [ sup t [ 0 , T ] e λ t 0 t ( f ( s , X s ε , Y s ε ) f ( s , X s ε , Y ^ s ε ) ) d s α 2 ] + C E [ sup t [ 0 , T ] e λ t 0 t ( f ( s , X s ε , Y ^ s ε ) f ( s , X s ( δ ) ε , Y ^ s ε ) ) d s α 2 ] + C E [ sup t [ 0 , T ] e λ t 0 t ( f ( s , X s ( δ ) ε , Y ^ s ε ) f ( s ( δ ) , X s ( δ ) ε , Y ^ s ε ) ) d s α 2 ] C 0 T ( E [ Y ε Y ^ s ( δ ) ε 2 + X ε X s ( δ ) ε 2 ] ) d s + C δ 2 C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) .

其次,对于第三项 D 3 ,有

D 3 C E [ sup t [ 0 , T ] e λ t ( τ 0 ( t s ) α 1 | X ^ ε ( t ) X ε ( t ) X ^ ε ( 0 ) + X ε ( 0 ) | d s ) 2 ] + C E [ sup t [ 0 , T ] e λ t ( τ 0 ( s ) α 1 | X ^ ε ( 0 ) X ε ( 0 ) ξ ( 0 ) + ξ ( 0 ) | d s ) 2 ] C E [ sup t [ 0 , T ] e λ t ( τ 0 ( t s ) α 1 | X ^ ε ( t ) X ε ( t ) | d s ) 2 ]

C E [ sup t [ 0 , T ] t 2 α e λ t | 0 t ( f ( s , X s ε , Y s ε ) f ( s , X s ε , Y ^ s ε ) ) d s | 2 ] + C E [ sup t [ 0 , T ] t 2 α e λ t | 0 t ( f ( s , X s ε , Y ^ s ε ) f ( s , X s ( δ ) ε , Y ^ s ε ) ) d s | 2 ] + C E [ sup t [ 0 , T ] t 2 α e λ t | 0 t ( f ( s , X s ( δ ) ε , Y ^ s ε ) f ( s ( δ ) , X s ( δ ) ε , Y ^ s ε ) ) d s | 2 ] C E [ sup t [ 0 , T ] t 1 2 α e λ t 0 t ( X s ε X s ( δ ) ε 2 + Y s ε Y ^ s ε 2 ) d s ] + C α , T δ 2 C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) .

综上所述,

E [ X ^ s ε X ε α , λ ( τ ) 2 ] C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) .

定理3.6得证。

Step3:对 X ^ ε X ¯ 进行估计。

R > 1 ,我们定义以下停时 τ R

τ R : = i n f { t 0 : B H | α , 0 , t R } T . (17)

定理3.7 假设(H1)~(H5)成立,那么对任意 t [ 0 , T ] ,存在正常数C和独立于 ε β > 0 ,使得对足够小的 δ ε ( 0 , 1 ) ,有

E [ X ^ ε X ¯ α , λ ( τ ) 2 ] C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) + C ( δ + ε δ 1 ) + C R 1 E [ B H α , 0 , T 2 ] .

证明由方程(8)和方程(22),有

E [ X ^ ε X ¯ α , λ ( τ ) 2 ] E [ X ^ ε X ¯ α , λ ( τ ) 2 1 { τ R < T } ] + E [ X ^ ε X ¯ α , λ ( τ ) 2 1 { τ R T } ] . (23)

对于不等式(23)右侧的第一项,根据切比雪夫不等式,可得

E [ X ^ ε X ¯ α , λ ( τ ) 2 1 { τ R < T } ] ( E [ X ^ ε X ¯ α , λ ( τ ) 4 ] ) 1 / 2 ( ( τ R < T ) ) 1 / 2 .

很容易得到 ( τ R < T ) ( B H α , 0 , T R ) R 1 E [ B H α , 0 , T 2 ]

根据文献 [25] 中的引理7.5中 lim R R 1 E [ B H α , 0 , T 2 ] = 0 ,可以得到

E [ X ^ ε X ¯ α , λ ( τ ) 2 1 { τ R < T } ] C R 1 E [ B H α , 0 , T 2 ] .

对于不等式(23)右侧的第二项,

E [ X ^ ε X ¯ α , λ ( τ ) 2 1 D ] C E [ sup t [ 0 , T ] e λ t | X ^ ε ( t ) X ¯ ( t ) | 2 1 D ] + C E [ sup t [ 0 , T ] e λ t ( 0 t ( t s ) α 1 | X ^ ε ( t ) X ¯ ( t ) X ^ ε ( s ) + X ¯ ( s ) | d s ) 2 1 D ] + C E [ sup t [ 0 , T ] e λ t ( τ 0 ( t s ) α + 1 | X ^ ε ( t ) X ¯ ( t ) X ^ ε ( s ) + X ¯ ( s ) | d s ) 2 1 D ] = : B 1 + B 2 + B 3 ,

其中 D : = { B H α , 0 , T R }

对于前两项 B 1 B 2 ,有

B 1 + B 2 C E [ sup t [ 0 , T ] e λ t 0 t ( f ( s ( δ ) , X s ( δ ) ε , Y ^ s ε ) f ¯ ( s ( δ ) , X s ( δ ) ε ) ) d s α 2 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( f ¯ ( s ( δ ) , X s ( δ ) ε ) f ¯ ( s , X s ( δ ) ε ) ) d s α 2 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( f ¯ ( s , X s ( δ ) ε ) f ¯ ( s , X s ε ) ) d s α 2 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( f ¯ ( s , X s ε ) f ¯ ( s , X ^ s ε ) ) d s α 2 1 D ]

+ C E [ sup t [ 0 , T ] e λ t 0 t ( f ¯ ( s , X ^ s ε ) f ¯ ( s , X ¯ s ε ) ) d s α 2 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( σ 1 ( s , X ^ ε ( s τ ) ) σ 1 ( s , X ¯ ( s τ ) ) ) d B s H α 2 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( σ 1 ( s , X ε ( s τ ) ) σ 1 ( s , X ^ ε ( s τ ) ) ) d B s H α 2 1 D ] = : i = 1 7 C i .

第一项 C 1 可写成

C 1 C E [ sup t [ 0 , T ] | 0 t ( f ( s ( δ ) , X s ( δ ) ε , Y ^ s ε ) f ¯ ( s ( δ ) , X s ( δ ) ε ) ) d s | 2 ] + C E [ sup t [ 0 , T ] ( 0 t ( t s ) α 1 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ s ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | d s ) 2 ] = : C 11 + C 12 .

对于 C 11 ,由条件(H1)和基本不等式,有

C 11 C E [ sup t [ 0 , T ] | k = 0 t δ 1 k δ ( k + 1 ) δ ( f ( k δ , X k δ ε , Y ^ s ε ) f ¯ ( k δ , X k δ ε ) ) d s | 2 ] + C E [ sup t [ 0 , T ] | t δ δ t ( f ( s ( δ ) , X s ( δ ) ε , Y ^ s ε ) f ¯ ( s ( δ ) , X s ( δ ) ε ) ) d s | 2 ]

C δ 2 + C E [ sup t [ 0 , T ] t δ k = 0 t δ 1 | k δ ( k + 1 ) δ ( f ( k δ , X k δ ε , Y ^ s ε ) f ¯ ( k δ , X k δ ε ) ) d s | 2 ] C δ 2 + C δ 2 max 0 k T δ 1 E [ | k δ ( k + 1 ) δ ( f ( k δ , X k δ ε , Y ^ s ε ) f ¯ ( k δ , X k δ ε ) ) d s | 2 ] C δ 2 + C ε 2 δ 2 max 0 k T δ 1 0 δ ε ζ δ ε J k ( s , ζ ) d s d ζ .

对于 C 12 ,根据Hölder不等式和 α ( 1 H , 1 / 2 ) ,可以得到

C 12 C E [ sup t [ 0 , T ] ( 0 t ( t s ) α 1 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ s ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | d s ) 2 ] C E [ sup t [ 0 , T ] 0 t ( t s ) α 1 2 d s × 0 t ( t s ) α 3 2 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 d s ] C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l c d s ] + C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l d s ] = : C 121 + C 122 .

其中 1 是一个指示函数, l : = { t < ( s δ + 2 ) δ } l c : = { t ( s δ + 2 ) δ }

由(H1)和条件 λ 1 λ 2 λ 1 λ 2 + 1 ,当 λ 1 λ 2 0 时,则

C 121 C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 | s ( s δ + 1 ) δ ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l c d s ] + C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 | t ( δ ) t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l c d s ] + C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 ( t δ s δ 1 ) k = s δ + 1 t δ 1 | k δ ( k + 1 ) δ ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l c d s ]

C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 ( ( s δ + 1 ) δ s ) | s ( s δ + 1 ) δ ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l c d s ] + C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 ( t t δ δ ) | t ( δ ) t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l c d s ] + C δ 1 E [ sup t [ 0 , T ] 0 t ( t t ( δ ) ) 1 2 α k = s δ + 1 t δ 1 | k δ ( k + 1 ) δ ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 l c d s ] C δ + C ε 2 δ 2 max 0 k T δ 1 0 δ ε ζ δ ε J k ( s , ζ ) d s d ζ .

对于 C 122 ,设 j : = { t δ > 1 } j c : = { t δ 1 } ,由条件(H1)和 t s < s δ δ s + 2 δ 2 δ ,有

C 122 C E [ sup t [ 0 , T ] 0 ( t δ 1 ) δ ( t s ) α 3 2 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 j l d s ] + C E [ sup t [ 0 , T ] ( t δ 1 ) δ t ( t s ) α 3 2 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 j l d s ] + C E [ sup t [ 0 , T ] 0 t ( t s ) α 3 2 | s t ( f ( r ( δ ) , X r ( δ ) ε , Y ^ r ε ) f ¯ ( r ( δ ) , X r ( δ ) ε ) ) d r | 2 1 j c l d s ] C δ 2 sup t [ 0 , T ] ( 0 ( t δ 1 ) δ ( t s ) α + 1 2 1 j l d s ) + C sup t [ 0 , T ] ( ( t δ 1 ) δ t ( t s ) α + 1 2 1 j l d s ) + C [ sup t [ 0 , T ] 0 t ( t s ) α + 1 2 1 j c l d s ] C δ 3 2 α .

因此, C 1 C α , T ε 2 δ 2 max 0 k T δ 1 0 δ ε ζ δ ε J k ( s , ζ ) d s d ζ

其中 0 ζ s δ ε J k ( s , ζ ) = E [ f ( k δ , X k δ ε , Y ^ s ε + k δ ε ) f ¯ ( k δ , X k δ ε ) , f ( k δ , X k δ ε , Y ^ ζ ε + k δ ε ) f ¯ ( k δ , X k δ ε ) ] (24)

下面,我们来估计 J k ( s , ζ ) :对 Y ^ ε ( t ) 作时间尺度变换,即对固定的k和 s [ 0 , δ ]

Y ^ ε ( s + k δ ) = Y ^ ε ( k δ ) + k δ s + k δ g ( X k δ ε , Y ^ ε ( r ) , Y ^ ε ( r τ ) ) d r + 0 s / ε σ 2 ( X k δ ε , Y ^ ε ( r ) , Y ^ ε ( r τ ) ) d W r = Y ^ ε ( k δ ) + 1 ε 0 s g ( X k δ ε , Y ^ ε ( r + k δ ) , Y ^ ε ( r + k δ τ ) ) d r + 1 ε 0 s σ 2 ( X k δ ε , Y ^ ε ( r + k δ ) , Y ^ ε ( r + k δ τ ) ) d W r , (25)

其中 W t * = W t + k δ W k δ 。对于固定的 ε > 0 r 0 ,设 Y X k δ ε , Y ^ k δ ε ( r / ε + θ ) = Y ^ ε ( r + k δ + θ ) θ [ τ , 0 ] 。令 W ¯ t 是布朗运动且与 W t 独立。构建过程 Y X k δ ε , Y ^ k δ ε

Y X k δ ε , Y ^ k δ ε ( s / ε ) = Y ^ ε ( k δ ) + 0 s / ε g ( X k δ ε , Y X k δ ε , Y ^ k δ ε ( r ) , Y X k δ ε , Y ^ k δ ε ( r τ ) ) d r + 0 s / ε σ 2 ( X k δ ε , Y X k δ ε , Y ^ k δ ε ( r ) , Y X k δ ε , Y ^ k δ ε ( r τ ) ) d W ¯ r = Y ^ ε ( k δ ) + 1 ε 0 s g ( X k δ ε , Y X k δ ε , Y ^ k δ ε ( r / ε ) , Y X k δ ε , Y ^ k δ ε ( r / ε τ ) ) d r + 1 ε 0 s σ 2 ( X k δ ε , Y X k δ ε , Y ^ k δ ε ( r / ε ) , Y X k δ ε , Y ^ k δ ε ( r / ε τ ) ) d W ˜ r ε , (26)

其中 W ˜ t ε = ε W ¯ t / ε

结合(25)和(26),可得 ( X k δ ε , Y ^ s + k δ ε ) ~ ( X k δ ε , Y s / ε k δ ε , Y ^ k δ ε ) s [ 0 , δ ) ,这里∼表示同分布。

对于当 s [ 0 , δ ) 时,由不等式(24),有

J k ( s , ζ ) = E [ f ( k δ , X k δ ε , Y s X k δ ε , Y ^ k δ ε ) f ¯ ( k δ , X k δ ε ) , f ( k δ , X k δ ε , Y ζ X k δ ε , Y ^ k δ ε ) f ¯ ( k δ , X k δ ε ) ] .

M k δ ε 是由 X k δ ε , Y ^ k δ ε 生成的σ-代数,它与 { Y r χ , η : r 0 } 独立。由文献 [28] 中的定理7.1.2,再结合引理2.7,可以得到 J k ( s , ζ ) C ( 1 + E [ X k δ ε 2 ] + E [ Y ^ k δ ε 2 ] ) e ρ 2 ( s ζ )

因此,对于 C 1 C 1 C ( ε δ 1 + δ )

通过定理3.6,估计 i = 2 5 C i ,可以得到

i = 2 5 C i C E [ sup t [ 0 , T ] e λ t 0 t ( t s ) α | f ¯ ( s ( δ ) , X s ( δ ) ε ) f ¯ ( s , X s ( δ ) ε ) | 2 d s 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( t s ) α | f ¯ ( s , X s ( δ ) ε ) f ¯ ( s , X s ε ) | 2 d s 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( t s ) α | f ¯ ( s , X s ε ) f ¯ ( s , X ^ s ε ) | 2 d s 1 D ] + C E [ sup t [ 0 , T ] e λ t 0 t ( t s ) 2 α | f ¯ ( s , X ^ s ) f ¯ ( s , X ¯ s ε ) | 2 d s 1 D ] C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) + C λ α 1 E [ X ^ ε X ¯ α , λ ( τ ) 2 1 D ] .

对于 C 6 ,有

C 6 C E [ sup t [ 0 , T ] 0 t e λ t [ ( t r ) 2 α + r α ] σ ^ 1 ( r , X ^ ε ( r τ ) ) σ 1 ( r , X ¯ ( r τ ) ) α 2 1 D d r ] C E [ sup t [ 0 , T ] 0 t e λ t [ ( t r ) 2 α + r α ] | σ ^ 1 ( r , X ^ ε ( r τ ) ) σ 1 ( r , X ¯ ( r τ ) ) | 2 1 D d r ] + C E [ sup t [ 0 , T ] 0 t e λ t [ ( t r ) 2 α + r α ] ( 0 r ( t u ) α 1 × | σ 1 ( r , X ^ ε ( r τ ) ) σ 1 ( r , X ¯ ε ( r τ ) ) σ 1 ( u , X ^ ε ( u τ ) ) + σ 1 ( u , X ¯ ( u τ ) ) | d u ) 2 1 D d r ] = : G 61 + G 62 .

对于 C 61 ,通过Hölder不等式,

G 61 C E [ sup t [ 0 , T ] 0 t e λ ( t r ) [ ( t r ) 2 α + r α ] sup q [ τ , r ] e λ r | X ^ ε ( u ) X ¯ ( u ) | 2 1 D d r ] C λ 2 α 1 E [ X ^ ε X ¯ α , λ ( τ ) 2 1 D ] .

根据(H2)和文献 [25] 的引理7.1,存在正常数C,使得

| σ ( t 1 , x 1 ) σ ( t 2 , x 2 ) σ ( t 1 , x 3 ) + σ ( t 2 , x 4 ) | C | x 1 x 2 x 3 + x 4 | + C | x 1 x 3 | | t 2 t 1 | β + C | x 1 x 3 | ( | x 1 x 2 | γ + | x 3 x 4 | γ ) .

对于 C 62 ,有

C 62 C E [ sup t [ 0 , T ] τ t τ e λ t [ ( t s τ ) 2 α + ( s + τ ) α ] × ( τ s | X ^ ε ( s ) X ¯ ( s ) X ^ ε ( u ) + X ¯ ( u ) | ( s u ) α + 1 d u ) 2 1 D d s ] + C E [ sup t [ 0 , T ] τ t τ e λ t [ ( t s τ ) 2 α + ( s + τ ) α ] × ( | X ^ ε ( s ) X ¯ ( s ) | τ s | X ^ ε ( s ) X ^ ε ( u ) | γ ( s u ) 1 + α d u ) 2 1 D d s ]

+ C E [ sup t [ 0 , T ] τ t τ e λ t [ ( t s τ ) 2 α + ( s + τ ) α ] × ( | X ^ ε ( s ) X ¯ ( s ) | τ s | X ¯ ( s ) X ¯ ( u ) | γ ( s u ) 1 + α d u ) 2 1 D d s ] C E τ t τ e λ t [ ( t s τ ) 2 α + ( s + τ ) α ] × ( 1 + Δ ( X ^ ε ) + Δ ( X ¯ ) ) 2 X ^ ε X ¯ α , λ ( τ ) 2 1 D d s ] ,

其中 Δ ( X ^ ε ) = sup s [ τ , T ] τ s | X ^ ε ( s ) X ^ ( u ) | γ ( s u ) α + 1 d u Δ ( X ¯ ) = sup s [ τ , T ] τ s | X ¯ ( s ) X ¯ ( u ) | γ ( s u ) α + 1 d u

由引理3.4,在 B H α , 0 , T R 的条件下,存在正常数C,使得

Δ ( X ^ ε ) + Δ ( X ¯ ) C Λ α γ ( B H ) [ ( 1 + X ^ ε α , ( τ ) ) γ + ( 1 + X ¯ α , ( τ ) ) γ ] sup s [ τ , T ] τ s ( s r ) ( 1 α ) γ 1 α d r C Λ α γ ( B H ) ( 1 + X ^ ε α , ( τ ) γ + X ¯ α , ( τ ) γ ) C .

其中 ( 1 α ) γ α > 0 α ( 0 , γ / 2 )

因此,可得 C 6 C 61 + C 62 C λ 2 α 1 E [ X ^ ε X ¯ α , λ ( τ ) 2 1 D ]

同理, C 7 的估计如下 C 7 C λ 2 α 1 E [ X ^ ε X α , λ ( τ ) 2 1 D ]

同理, B 3 的估计为

B 3 C ( ε δ 1 + δ ) + C λ 2 α 1 E [ X ^ ε X ¯ α , λ ( τ ) 2 1 D ] + C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) .

综上所述,

E [ X ^ ε X ¯ α , λ ( τ ) 2 ] C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) + C ( δ + ε δ 1 ) + C R 1 E [ B H α , 0 , T 2 ] .

引理3.7得证。

Step4:对 X ε X ¯ 进行估计。

应用基本不等式,由式(8)和(12)可以得到

E [ X ε X ¯ α , λ ( τ ) 2 ] 2 E [ X ε X ^ ε α , λ ( τ ) 2 ] + 2 E [ X ^ ε X ¯ α , λ ( τ ) 2 ] .

根据定理3.6和定理3.7,有

E [ X ε X ¯ α , ( τ ) 2 ] C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) + C ( δ + ε δ 1 ) + C R 1 E [ B H α , 0 , T 2 ] .

综上所述,

E [ X ε X ¯ α , ( τ ) 2 ] e λ T E [ X ε X ¯ α , λ ( τ ) 2 ] C δ 2 ( 1 + δ 1 2 α ( 1 + ε 1 e β δ ε ) ) + C ( δ + ε δ 1 ) + C R 1 E [ B H α , 0 , T 2 ] .

δ = ε ( ln ε ) 1 / 2 ,当 R 时,有

lim ε 0 E [ X ε X ¯ α , ( τ ) 2 ] = 0.

综上所述,即证明了定理3.1。

4. 结束语

在本篇论文中,我们研究了由分数布朗运动驱动的非自治双尺度随机时滞微分方程的平均原理,其中该方程的慢分量由Hurst参数 H ( 1 / 2 , 1 ) 的分数布朗运动驱动,快分量由布朗运动驱动。首先,我们使用广义Stieltjes积分对快慢过程进行估计,得到方程(2.1)的解 ( X ε , Y ε ) 的矩估计。其次,我们利用Khasminskii技巧,构造了辅助过程 ( X ^ t ε , Y ^ t ε ) ,估计了 X t ε X ^ t ε Y t ε Y ^ t ε 这两项。然后,基于快方程的指数遍历性,我们引入停时序列来控制分数布朗运动的大小进而估计辅助过程 X ^ t ε 与平均方程解 X ¯ t 之间的误差。最后,证明了慢分量 X t ε 在均方意义下收敛于平均方程的解 X ¯ t

基金项目

长安大学中央高校基本科研业务费专项资金资助(CHD300102122113),陕西省自然科学基础研究计划项目(2023-JC-QN-0009)。

文章引用

贺 鑫. 双尺度随机时滞微分方程的平均原理
The Averaging Principle of Two-Scale Stochastic Delay Differential Equation[J]. 应用数学进展, 2024, 13(02): 788-805. https://doi.org/10.12677/AAM.2024.132077

参考文献

  1. 1. Larter, R., Steinmetz, C.G. and Aguda, B.D. (1988) Fast-Slow Variable Analysis of the Transition to Mixed-Mode Os-cillations and Chaos in the Peroxidase Reaction. Journal of Chemical Physics, 89, 6506-6514. https://doi.org/10.1063/1.455370

  2. 2. Chow, P.L. (1973) Thermoelastic Wave Propagation in a Random Medium and Some Related Problems. International Journal of Engineering Science, 11, 953-971. https://doi.org/10.1016/0020-7225(73)90010-4

  3. 3. Krauskopf, B., Osinga, H.M. and Galan-Vioque, J. (2007) Numerical Continuation Methods for Dynamical Systems: Path Following and Boundary Value Problems. Springer, Dordrecht. https://doi.org/10.1007/978-1-4020-6356-5

  4. 4. Dubbeldam, J.L. and Krauskopf, B. (1999) Self-Pulsations of Lasers with Saturable Absorber: Dynamics and Bifurcations. Optics Communications, 159, 325-338. https://doi.org/10.1016/S0030-4018(98)00568-9

  5. 5. Bogoliubov, N.N. and Mitropolsky, Y.A. (1961) Asymptotic Methods in the Theory of Non-Linear Oscillations. Gordon and Breach Science Publishers, New York.

  6. 6. Gikhman, I.I. (1947) A Method of Constructing Random Process. Proceedings of the USSR Academy of Sciences, 58, 961-964.

  7. 7. Volosov, V.M. (1962) Averaging in Systems of Ordinary Differential Equations. Uspekhi Matematich-eskikh Nauk, 108, 3-126. https://doi.org/10.1070/RM1962v017n06ABEH001130

  8. 8. Besjes, J.G. (1968) On the Asymptotic Methods for Non-Linear Differential Equations. Journal of Mecanique, 8, 357-372.

  9. 9. Khasminskii, R.Z. (1968) On the Principle of Averaging the Itô Stochastic Differential Equations. Kybernetika, 4, 260-279.

  10. 10. Freidlin, M.I. and Wentzell, A.D. (1998) Random Perturbation of Dynamical Systems. 2nd Edition, Springer-Verlag, New York. https://doi.org/10.1007/978-1-4612-0611-8

  11. 11. Freidlin, M.I. and Wentzell, A.D. (2006) Long-Time Behavior of Weakly Coupled Oscillators. Journal of Statistical Physics, 123, 1311-1337. https://doi.org/10.1007/s10955-006-9133-8

  12. 12. Golec, J. and Ladde, G. (1990) Averaging Principle and Systems of Singularly Perturbed Stochastic Differential Equations. Journal of Statistical Mechanics, 31, 1116-1123. https://doi.org/10.1063/1.528792

  13. 13. Givon, D., Kevrekidis, I.G. and Kupferman, R. (2006) Strong Convergence of Projective Integration Schemes for Singularly Perturbed Stochastic Differential Systems. Communications in Mathe-matical Sciences, 4, 707-729. https://doi.org/10.4310/CMS.2006.v4.n4.a2

  14. 14. Cerrai, S. and Freidlin, M. (2009) Averaging Principle for A Class of Stochastic Reaction-Diffusion Equations. Probability Theory and Related Fields, 144, 137-177. https://doi.org/10.1007/s00440-008-0144-z

  15. 15. Fu, H.B. and Duan, J.Q. (2011) An Averaging Principle for Two Time-Scales Stochastic Partial Differential Equations. Stochastics and Dynamics, 11, 353-367. https://doi.org/10.1142/S0219493711003346

  16. 16. Fu, H.B. and Liu, J.C. (2011) Strong Convergence in Stochastic Averaging for Two Time-Scales Stochastic Partial Differential Equations. The Journal of Mathematical Analysis and Ap-plications, 384, 70-86. https://doi.org/10.1016/j.jmaa.2011.02.076

  17. 17. Xu, J. (2017) Lp-Strong Convergence of the Averaging Principle for Slow-Fast SPDEs with Jumps. The Journal of Mathematical Analysis and Applications, 445, 342-373. https://doi.org/10.1016/j.jmaa.2016.07.058

  18. 18. Xu, Y., Pei, B. and Guo, R. (2015) Stochastic Averaging for Slow-Fast Dynamical Systems with Fractional Brownian Motion. Discrete and Continuous Dynamical Systems-Series B, 20, 2257-2267. https://doi.org/10.3934/dcdsb.2015.20.2257

  19. 19. 郭蓉. 分数阶时滞微分方程的随机平均法[J]. 山西大学学报(自然科学版), 2019, 42(2): 300-306.

  20. 20. Mao, W., You, S., Wu, X. and Mao, X. (2015) On the Averaging Principle for Stochastic Delay Differential Equations with Jumps. Advances in Difference Equations, 2015, Article No. 70. https://doi.org/10.1186/s13662-015-0411-0

  21. 21. Wang, P.G. and Xu, Y. (2020) Averaging Method for Neutral Stochastic Delay Differential Equations Driven by Fractional Brownian Motion. Journal of Function Spaces, 2020, Arti-cle ID: 5212690. https://doi.org/10.1155/2020/5212690

  22. 22. Bao, J., Song, Q., Yin, G. and Yuan, C. (2017) Ergodicity and Strong Limit Results for Two-Time-Scale Functional Stochastic Differential Equations. Stochastic Analysis and Applications, 35, 1030-1046. https://doi.org/10.1080/07362994.2017.1349613

  23. 23. Han, M., Xu, Y., Pei, B. and Wu, J.L. (2022) Two-Time-Scale Stochastic Differential Delay Equations Driven by Multiplicative Fractional Brownian Noise: Averaging Principle. Journal of Mathematical Analysis and Applications, 510, Article ID: 126004. https://doi.org/10.1016/j.jmaa.2022.126004

  24. 24. Zähle, M. (1998) Integration with Respect to Fractal Functions and Stochastic Calculus. I. Probability Theory and Related Fields, 111, 333-374. https://doi.org/10.1007/s004400050171

  25. 25. Nualart, D. and Rășcanu, A. (2002) Differential Equations Driven by Fractional Brownian Motion. Collectanea Mathematica, 53, 55-81.

  26. 26. Atienza, M.J., Lu, K. and Schmalfuß, B. (2010) Random Dynamical Systems for Stochastic Partial Differential Equations Driven by a Fractional Brownian Mo-tion. Discrete and Continuous Dynamical Systems-Series B, 14, 473-493. https://doi.org/10.3934/dcdsb.2010.14.473

  27. 27. Shevchenko, G. (2013) Mixed Stochastic Delay Differential Equa-tions. Theory of Probability and Mathematical Statistics, 89, 181-195. https://doi.org/10.1090/S0094-9000-2015-00944-3

  28. 28. Oksendal, B. (2003) Stochastic Differential Equations. 6th Edition, Springer, Berlin.

期刊菜单