α-稳定过程驱动的非线性随机微分方程解的存在唯一性
Existence and Uniqueness Solution of Nonlinear SDE Driven by α-Stable Process
摘要: 本文的主要目的是研究α-稳定过程驱动的非线性随机微分方程解的存在唯一性。首先我们给出了漂移系数和扩散系数都是高度非线性的随机微分方程解存在唯一的假设条件,其次我们证明了该解的存在唯一性,最后,我们给出了在有限时间内解的收敛性证明。
Abstract: The main purpose of this paper is to investigate the existence and uniqueness solution of nonlinear SDE driven by α-stable process. Firstly, we give the assumptions to obtain the existence and uniqueness solution of such model which both drift coefficient and diffusion coefficient are both highly nonlinear. Then we prove the existence and uniqueness of the solution. Finally, we prove the convergence of the solution in finite time.
文章引用:陆芸芳, 童金英. α-稳定过程驱动的非线性随机微分方程解的存在唯一性[J]. 理论数学, 2019, 9(1): 54-61. https://doi.org/10.12677/PM.2019.91008

1. 引言

近些年来,由布朗运动驱动的随机微分方程已经在很多文章中被广泛讨论了,例如Mao [1] ,Mao [2] 等。考虑标量随机微分方程

d x ( t ) = ( a x ( t ) b x 3 ( t ) ) d t + c x ( t ) d B ( t ) (1.1)

其中 B ( t ) 是一个标准的布朗运动, a , b , c 是三个正常数,我们称这样的随机微分方程为Ginzburg-Landan方程。这一类的随机微分方程有很多好的性质,例如: L r 强收敛性,解的二阶矩渐近稳定性(见Guo [3] )。

最近,关于 d x ( t ) = f ( ( x ( t ) ) d t + g ( x ( t ) ) ) d B ( t ) (highly nonlinear SDEs)的数值解得到了广泛的讨论。其中 f g 都是 R R 上的非线性函数。然而在现实世界中,许多问题都需要我们用α-稳定过驱动的随机微分方程搭建的数学模型来解决,例如:Mandelbrot [4] 指出1890年至1937年月度羊毛价格变化的分布遵循α = 1.7的α-稳定分布。从Embrechts [5] 我们可知SDE驱动因子为0.018时,对于 d x ( t ) = 0.018 x ( t ) d Z ( t ) ,初始值为 x 0 = 1 x ( t ) 的取值可能随着时间而取负值,也就是说这个SDE的数值解随时会爆炸,然而现实生活中的许多问题都需要我们用收敛的数值解去解决,但是,在讨论数值解之前,我们必须确保SDE的解是存在的,所以我们有足够的理由去讨论以下SDE:

d x ( t ) = ( a x ( t ) b x 3 ( t ) ) d t + c x ( t ) d Z ( t ) (1.2)

其中 Z ( t ) 是一个α-稳定过程,且 0 < α < 2 R * 上的Lévy测度 ν 定义为:

ν α ( d z ) = c 1 | z | α + 1 1 ( 0 , + ) ( z ) d z + c 2 | z | α + 1 1 ( , 0 ) ( z ) d z

N ( d s , d z ) [ 0 , + ) × R * 上的泊松测度, N ˜ ( d s , d z ) 为补偿泊松测度,写做:

N ˜ ( d s , d z ) = N ( d s , d z ) ν ( d z ) d s

α-稳定过程 ( Z ( t ) ) t 0 定义如下:

{ Z ( t ) = 0 t R * z N ( d s , d z ) , α ( 0 , 1 ) Z ( t ) = 0 t R * z N ˜ ( d s , d z ) , α ( 1 , 2 ) ,

过程 ( Z ( t ) ) t 0 含有参数 c 1 , c 2 ,当 c 1 = c 2 时,我们称该过程为对称α稳定过程。本文仅对对称α稳定过程进行讨论,且规定 α ( 1 , 2 ) ,也就是说

ν ( d z ) = C α | z | α + 1 d z ,

其中 C α = α 2 α - 1 Γ ( ( 1 + α ) / 2 ) π 1 / 2 Γ ( 1 α / 2 ) 以及 Γ ( ) 是一个Gamma函数,定义为:

Γ ( s ) = 0 + t s 1 e t d t , s ( 0 , + ) .

更为一般的,我们讨论漂移算子和积分算子都是非线性的随机微分方程,即:

d x ( t ) = f ( x ( t ) ) d t + g ( x ( t ) ) d Z ( t ) (1.3)

其中 Z ( t ) 是一个定义在 α ( 1 , 2 ) t 0 上的对称α-稳定过程,规定初始值为 x ( 0 ) = x 0 R ,其中

f : R R , g : R R .

假设 f g 能够分解为

f ( x ) = F 1 ( x ) + F ( x ) , g ( x ) = G 1 ( x ) + G ( x ) ,

其中函数 F 1 , F : R R G 1 , G : R R .

由α-稳定过程驱动的SDE在近几年已经被很多学者探讨(参见Applebaum [6] ),然而,目前为止还没有人给出SDE(1.3)的数值解,在此之前,讨论SDE (1.3)解的存在唯一性就成了我们需要解决的问题。在解决这个问题之前,我们先来对文中出现的符号和需要用到的预备知识进行说明。

2. 符号说明和预备知识

2.1. 符号说明

首先,我们队文本中出现的随机过程进行符号解释。 ( Ω , F , P ) 是一个完备概率空间,其中 F 是域流,满足 F t 右连续并且 F 0 包含了所有 P 空集。 E 表示 P 对应的数学期望。 B ( t ) 表示定义在 ( Ω , F , P ) 上的标准布朗运动。 Z ( t ) 表示定义在 ( Ω , F , P ) 上的α-稳定过程。如果 G 是一个集合,它的示性函数表示为当 x G 时, 1 G ( x ) = 1 否则为0。我们用 R R * 表示实数集以及不包含0的实数集,定义

N ˜ ( d t , d z ) = N ( d s , d z ) ν ( d z ) d t 以及 x ( t ) = lim s t x ( s ) .

2.2. α-稳定过程的无穷小生成元

我们介绍由谱正稳定过程驱动的带马尔科夫调制穷小生成元。

只有正跳的稳定过程叫做谱正稳定过程,谱正稳定过程的Lévy测度 ν

ν ( d z ) = { C α 1 ( z > 0 ) z α + 1 d z , z 0 , 0 , z = 0.

对于由 α ( 1 , 2 ) 的普正稳定过程驱动的模型

d x ( t ) = f ( x ( t ) ) d t + g ( x ( t ) ) d Z t , t [ 0 , T ] ,

其中 f ( x ( t ) ) g ( x ( t ) ) 均为 R 上的函数。

对任意 V ( x ) R + ,随机过程 x ( t ) 的无穷小生成元 L V 满足:

L V ( x ) = f ( x ) V ( x ) x + R * [ V ( x + g ( x ) z ) V ( x ) V ( x ) x g ( x ) z 1 { | z | 1 } ] ν (dz)

3. 解的存在并且唯一满足的假设条件

为了方便下文的证明,我们提出以下几个基本假设。

3.1. 假设A1

假设 F 1 , F , G 1 , G 满足以下条件:

| F 1 ( x ) F 1 ( y ) | | G 1 ( x ) G 1 ( y ) | L 1 | x y | (3.1)

| F ( x ) F ( y ) | | G ( x ) G ( y ) | L 1 ( 1 + | x | γ + | y | γ ) | x y | (3.2)

从(3.1)可得线性增长条件,即存在一个 K 1 > 0 使得

| F 1 ( x ) | | G 1 ( x ) | K 1 ( 1 + | x | ) (3.3)

对所有的 x R 都成立。

3.2. 假设A2

假设 g ( x ) 满足条件:对 α ( 1 , 2 ) , z ( 0 , 1 ] , p ¯ ( 0 , 1 ] 以及 p ( 0 , p ¯ ) ,存在一对常数 L 2 0 θ 0 ,令 ω ( x ) = min ( | | x | | g ( x ) z | | ) ,则

( C α 2 ( 2 α ) ( | g ( x ) | ω ( x ) ) 2 p C α α p ) L 2 ( x 2 + θ 2 ) p 2 2 , x ( 0 , 1 ) . (3.4)

3.3. 假设A3

假设 f ( x ) g ( x ) 满足条件:存在常数 θ > 0 , p ( 0 , p ¯ ) , M 1 sgn ( g ( x ) ) L 2 L 3 > 0 使得

p x f ( x ) + M 1 | g ( x ) | p L 3 ( x 2 + θ 2 ) (3.5)

对所有 x R 成立。

4. 解存在唯一性定理及解的收敛性

这小结分为两部分,第一部分给出SDE (1.3)解存在唯一的证明,第二部分给出解收敛的证明。

注:从此刻起, C 代表一个取值依赖于 T , p , L 1 , L 2 , L 3 , X 0 的常数。

4.1. 解的存在唯一性

定理B1:若假设A1~A3成立,则SDE (1.3)有一个唯一的全局解 x ( t )

证明:令 k 0 > 0 充分大且 x 0 < k 0 。对各个整数规定 k > k 0 ,定义停时

τ k = inf { t [ 0 , τ e ) : x ( t ) ( 1 k , k ) }

因为SDE (1.3)的系数满足局部Lipschitz条件,所以,对于给定的初始值 x 0 R ,在 t [ 0 , τ e ) 上存在一个唯一的局部解 x ( t ) ,其中 τ e 是爆炸的时刻。(参见Bao等 [7] ,Bao和Yuan [8] ,Bass [9] )。令 k 0 > 0 充分大且 x 0 < k 0 。对于各个整数 k > k 0 ,设 τ = lim k τ k ,其中 τ < τ e a.s.若我们得到 τ = a.s.则 x ( t ) R a.s.对所有的 t 0 。对于 θ > 0 , 0 < p < 1 ,令

V θ ( x ) = ( x 2 + θ 2 ) p 2 , x R (3.6)

T > 0 ,对任意 0 t T ,由 I t o ^ 公式(参见Applebaum [6] )可得

V θ ( x ( t τ k ) ) = V θ ( x 0 ) + 0 t τ k L V θ ( x ( s ) ) d s + M 1 ( t τ k ) (3.7)

其中 L V θ : R R 定义为:

L V θ = p x ( x 2 + θ 2 ) p 2 2 f ( x ) + 0 1 [ V θ ( x + g ( x ) z ) V θ ( x ) V θ ( x ) g ( x ) z 1 { 0 z < 1 } ] C α | z | 1 + α d z = : A ( x ) + B ( x ) + C ( x ) (3.8)

其中 M 1 ( t ) 是一个局部鞅,定义如下:

M 1 ( t τ k ) = 0 t τ k R * ( ( V θ ( x ( s ) ) + g ( x ( s ) ) ) V θ ( x ( s ) ) ) N ˜ ( d s , d z )

A ( x ) = p x ( x 2 + θ 2 ) p 2 2 f ( x ) , B ( x ) = 0 1 [ V θ ( x + g ( x ) z ) V θ ( x ) V θ ( x ) g ( x ) z ] C α | z | 1 + α d z , C ( x ) = 0 + [ V θ ( x + g ( x ) z ) V θ ( x ) ] C α | z | 1 + α d z

g ( x ) z = y ,首先,我们计算 B (x)

由Taylor公式可知,存在一个 ξ 取值介于 x x + y 之间,使得

B ( x ) = 0 1 [ V θ ( x + g ( x ) z ) V θ ( x ) V θ ( x ) g ( x ) z ] C α | z | 1 + α d z sgn ( g ( x ) ) | g ( x ) | α 0 g ( x ) [ V θ ( x + y ) V θ ( x ) V θ ( x ) y ] C α | y | 1 + α d y = sgn ( g ( x ) ) | g ( x ) | α 0 g ( x ) 1 2 V θ ' ' ( ξ ) C α | y | 1 + α d y (3.9)

其中

V θ ( ξ ) = p ( ξ 2 + θ 2 ) p 2 2 ( ( p 1 ) ξ 2 + θ 2 ) (3.10)

由基本不等式可知

ξ 2 θ 2 ( p 1 ) ξ 2 + θ 2 ξ 2 + θ 2

从这个式子我们立即得到

| ( p 1 ) ξ 2 + θ 2 | | ξ 2 + θ 2 | (3.11)

| x | | y | ,即 | g ( x ) z | | x | ,将式(3.11)带入式(3.10)我们得到

V θ ( ξ ) p ( ξ 2 + θ 2 ) p 2 1 p ξ p 2

结合式(3.9)可得

B ( x ) sgn ( g ( x ) ) | g ( x ) | α 0 g ( x ) 1 2 y 2 p | | x | | g ( x ) z | | p 2 C α | y | 1 + α d y sgn ( g ( x ) ) C α p | ω ( x ) | p 2 2 ( 2 α ) | g ( x ) | 2 (3.12)

| g ( x ) z | = | x | ,从式(3.9)我们可得

B ( x ) = sgn ( g ( x ) ) | g ( x ) | α 0 g ( x ) 1 2 y 2 V θ ' ' C α | y | 1 + α d y sgn ( g ( x ) ) | g ( x ) | α 0 g ( x ) 1 2 y 2 p ( ξ 2 + θ 2 ) p 2 2 C α | y | 1 + α d y sgn ( g ( x ) ) | g ( x ) | α 0 g ( x ) 1 2 y 2 p ( | | x | | g ( x ) z | | + θ 2 ) p 2 2 C α | y | 1 + α d y C α p θ p 2 2 ( 2 α ) | g ( x ) | 2 (3.13)

从我们选择的p值不难看出

lim θ C α p θ p 2 2 ( 2 α ) | g ( x ) | 2 = 0 (3.14)

即,对一个充分大的 θ B ( x ) | g ( x ) z | = | x | 条件下是收敛的。

接下来,我们处理 C ( x )

f ( θ ) = V θ ( x + y ) V θ ( x ) ,则

f ( θ ) = p θ [ ( ( x + y ) 2 + θ 2 ) p 2 2 ( x 2 + θ 2 ) p 2 2 ]

| x + y | > x , f ( θ ) < 0 f ( θ ) < 0 ,可得 f ( θ ) < f (0)

| x + y | < x , f ( θ ) < 0 f ( θ ) > 0 ,可得 f ( θ ) = | f ( θ ) | < f (0)

那么,我们的到

| f ( θ ) | = | V θ ( x + y ) V θ ( x ) | < | f ( 0 ) | = | ( x + y ) p x p |

结合我们对 p 的选择,由基本不等式可得

| V θ ( x + y ) V θ ( x ) | < | y | p (3.15)

C ( x ) = 1 + [ V θ ( x + g ( x ) z ) V θ ( x ) ] C α | z | 1 + α d z sgn ( g ( x ) ) | g ( x ) | α g ( x ) + | V θ ( x + y ) V θ ( x ) | C α | y | 1 + α d y sgn ( g ( x ) ) | g ( x ) | α g ( x ) + | y | p C α | y | 1 + α d y = sgn ( g ( x ) ) | g ( x ) | α C α α p | g ( x ) | p (3.16)

x ( 1 , 1 ) ,将式(3.12)和式(3.16)加入式(3.8)我们可以得到

L V θ ( x ) p x ( x 2 + θ 2 ) p 2 2 f ( x ) + sgn ( g ( x ) ) C α p | ω ( x ) | p 2 2 ( 2 α ) | g ( x ) | 2 + sgn ( g ( x ) ) C α α p | g ( x ) | p = p x ( x 2 + θ 2 ) p 2 2 f ( x ) + sgn ( g ( x ) ) C α 2 ( 2 α ) | g ( x ) | p ( | g ( x ) | | ω ( x ) | ) 2 p + sgn ( g ( x ) ) C α α p | g ( x ) | p = p x ( x 2 + θ 2 ) p 2 2 f ( x ) + sgn ( g ( x ) ) C α | g ( x ) | p ( 1 2 ( 2 α ) ( | g ( x ) | | ω ( x ) | ) 2 p + 1 α p ) (3.17)

结合假设式(3.2),可得

L V θ ( x ) p x ( x 2 + θ 2 ) p 2 2 f ( x ) + L 2 sgn ( g ( x ) ) C α | g ( x ) | p ( x 2 + θ 2 ) p 2 2 ( x 2 + θ 2 ) p 2 2 ( p x f ( x ) + M 1 | g ( x ) | p ) ( x 2 + θ 2 ) p 2 2 L 3 ( x 2 + θ 2 ) C V θ ( x ) (3.18)

x = 0 V θ ( x ) = V θ ( 0 ) = θ p ,则, V θ ( x ) = 0

L V θ ( x ) L 2 sgn ( g ( 0 ) ) C α | g ( 0 ) | p ( θ 2 ) p 2 2 C

I t o ^ 公式,依据 x ( t τ k ) R 0 t T

E V θ ( x ( T τ k ) ) = V θ ( x 0 ) + E 0 t τ k L V θ ( x ( s ) ) d s (3.19)

则,由式(3.19)我们可得

E V θ ( x ( T τ k ) ) V θ ( x 0 ) + C E ( T τ k ) V θ ( x 0 ) + C T (3.20)

注意到对于所有的 σ { τ k T } ,对于一个充分大的数 k ,满足 x ( τ k , σ ) k 或者 x ( τ k , σ ) 1 k

由于 V θ ( x ) x ( , 0 ) 是非增的并且在 x ( 0 , + ) 是非减的,所以 V θ ( x ( τ k , σ ) ) 不小于 ( k 2 + θ 2 ) p 2 ( 1 / k 2 + θ 2 ) p 2 ,即

V θ ( x ( τ k , σ ) ) ( k 2 + θ 2 ) p 2 ( 1 / k 2 + θ 2 ) p 2 (3.21)

从式(3.2.23)和式(3.2.24)可得

( k 2 + θ 2 ) p 2 ( 1 / k 2 + θ 2 ) p 2 P ( τ k T ) E ( V θ ( τ k , σ ) 1 { τ k T } ) V θ ( x 0 ) + C T

k 可得

P ( τ T ) = 0

T 的任意性可知

P ( τ = ) = 1

综上,SDE (1.3)在 t 0 上存在一个唯一的全局解 x ( t ) R

4.2. 解的收敛性

若SDE (1.3)的解存在且唯一,即,定理B1成立下,对任意的 p ¯ ( 0 , 1 ] , p ( 0 , p ¯ ) ,有

sup 0 t T E | x ( t ) | p < C , T > 0. (4.1)

证明:由 I t o ^ 公式,对所有的 x ( 0 ) = x 0 R 0 t T 满足

E V θ ( x ( t ) ) = V θ ( x 0 ) + E 0 t L V θ ( x ( s ) ) d s (4.2)

由式(3.18)可得

E V θ ( x ( t ) ) = V θ ( x 0 ) + E 0 t L V θ ( x ( s ) ) d s V θ ( x 0 ) + E C V θ ( x ( s ) ) d s

由Gronwall不等式

E | x ( t ) | p E V θ ( x ( t ) ) V θ ( x 0 ) exp ( C t ) V θ ( x 0 ) exp ( C T ) C

证毕。

参考文献

[1] Mao, X. (2015) The Truncated Euler-Maruyama Method for Stochastic Differential Equations. Journal of Computational and Applied Mathematics, 290, 370-384.
https://doi.org/10.1016/j.cam.2015.06.002
[2] Mao, X. (2016) Convergence Rates of the Truncated Euler-Maruyama Method for Stochastic Differential Equations. Journal of Com-putational and Applied Mathematics, 296, 362-375.
https://doi.org/10.1016/j.cam.2015.09.035
[3] Guo, Q., Liu, W. and Mao, X. (2017) The Partially Truncated Euler-Maruyama Method and Its Stability and Boundedness. Applied Numerical Mathematics, 115, 235-251.
[4] Mandelbrot, B.B. (1963) The Variation of Certain Speculative Prices. The Journal of Business, 36, 394-419.
https://doi.org/10.1086/294632
[5] Embrechts, P. (2012) Stochastic Modelling and Applied Probability. Springer, Germany.
[6] Applebaum, D. (2004) Lévy Process and Stochastic Calculus. Cambridge University Press, Cambridge.
https://doi.org/10.1017/CBO9780511755323
[7] Bao, J. and Yuan, C. (2011) Comparison Theorem for Sto-chastic Differential Delay Equations with Jumps. Springer, Germany.
[8] Bao, J. and Yuan, C. (2012) Stochastic Population Dynamics Driven by Lévy Noise. Journal of Mathematical Analysis and Applications, 391, 365-375.
https://doi.org/10.1016/j.jmaa.2012.02.043
[9] Bass, R. (1997) Diffusions and Elliptic Operators. Springer, New York.