馬 麗,嚴(yán)良清,韓新方

(海南師范大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,海南?？?71158)

1 引言

中立延遲隨機(jī)微分方程歷年來(lái)在生物、工程、金融等各個(gè)領(lǐng)域引起了學(xué)者的廣泛關(guān)注.文[1]系統(tǒng)地介紹了不帶跳的隨機(jī)泛函微分方程的基本理論及其在金融、隨機(jī)游戲、人口問(wèn)題中的應(yīng)用;文[2]給出了帶有Lévy 跳隨機(jī)泛函微分方程解的存在唯一性;文[3]得到了一類(lèi)帶Lévy 跳的中立隨機(jī)泛函微分方程解的存在唯一性;文[4]研究了帶有特殊跳(泊松跳)的中立隨機(jī)延遲微分方程的數(shù)值逼近;文[5]得到了Lévy 噪聲擾動(dòng)的混合隨機(jī)微分方程的Euler 近似解;文[6]和文[7]研究了帶Markov 狀態(tài)轉(zhuǎn)換的跳擴(kuò)散方程的數(shù)值解.

設(shè)(?,F,P)是完備概率空間,(Ft)t≥0是其上一個(gè)滿足通常條件的適應(yīng)流. 設(shè){=(t),t ≥0}是一個(gè)關(guān)于(Ft)t≥0適應(yīng)的穩(wěn)定的Rn值泊松點(diǎn)過(guò)程. 設(shè)B(Rn?{0})為Rn?{0}上的波萊爾σ-代數(shù),對(duì)A ∈B(Rn?{0}),定義與聯(lián)系的泊松計(jì)數(shù)測(cè)度N(t,A)=N((0,t]×A)如下

則存在一個(gè)σ 有限測(cè)度π 使得

這里的測(cè)度π 稱(chēng)為L(zhǎng)évy 測(cè)度.由Doob-Meyer 分解定理,存在關(guān)于(Ft)t≥0適應(yīng)的唯一的鞅和唯一的增過(guò)程使得

設(shè)|·|表示歐式空間Rd中的范數(shù),τ 為一個(gè)正的固定的延遲,C([?τ,0],Rd)為[?τ,0]到Rd上的連續(xù)函數(shù)類(lèi),其上的范數(shù)為設(shè)ξ(t)為關(guān)于F0可測(cè)的C([?τ,0],Rd)隨機(jī)變量且滿足Eξp<∞,其中p 為大于等于2 的任意正整數(shù). 設(shè)W(t)是(?,F,P)上關(guān)于流(Ft)t≥0適應(yīng)的標(biāo)準(zhǔn)的r 維布朗運(yùn)動(dòng)且與跳N 獨(dú)立. 設(shè)Z ∈B(Rn?{0})且π(Z)<∞,設(shè)0

在系數(shù)滿足局部Lipschitz 條件和線性增長(zhǎng)條件,中立項(xiàng)D(t,x(t ?τ))關(guān)于第二個(gè)分量為壓縮映射的條件下,類(lèi)似于文[2],我們可得方程(1.1)存在唯一解.由于解沒(méi)有顯示表達(dá),因此有必要研究其數(shù)值解.如果數(shù)值解逼近于真實(shí)解,我們可以用數(shù)值解來(lái)估計(jì)真實(shí)解.

本文內(nèi)容安排如下:第二節(jié)介紹了方程(1.1)的Euler 的數(shù)值算法,并給出主要結(jié)果即定理1;第三節(jié)給出定理1 的證明.本文推廣了文[4]的結(jié)果,考慮的中立項(xiàng)是時(shí)間和狀態(tài)的二元函數(shù),在逼近的時(shí)候?qū)χ辛㈨?xiàng)需要加一定的條件才可以放縮.此外,有中立項(xiàng)時(shí)需要把它看成一個(gè)整體,進(jìn)而用伊藤公式,再由基本不等式及壓縮映射最終得到數(shù)值解穩(wěn)定于真實(shí)解.方程(1.1)中f,g,h 也依賴于時(shí)間,因此需要三個(gè)函數(shù)關(guān)于時(shí)間t 是局部Lipschitz 的.

2 EM 近似解及主要結(jié)果

下面給出(1.1)式的EM 逼近解.

其中?Wk=Wtk+1?Wtk.假設(shè)則EM 逼近解y(t)的連續(xù)形式如下

本文對(duì)系數(shù)做如下假設(shè).

(H1) 對(duì)任意的正整數(shù)m,存在正整數(shù)使得對(duì)任意的t1,t2∈[0,+∞),任意的且|x|≤m,|y|≤m,|有

(H2)對(duì)任意的p ≥2,存在正數(shù)k1使得對(duì)任意t ∈[0,+∞),x,y ∈Rn,有

(H3)存在正數(shù)k2∈(0,1),使得對(duì)任意t1,t2∈[0,+∞),x,y ∈Rn,有

設(shè)T ∈[0,+∞)為任一常數(shù),本文主要結(jié)果如下.

定理1在(H1)–(H3)條件下,方程(1.1)的Euler 數(shù)值解收斂到真實(shí)解.即

3 定理1 的證明

在證明定理1 之前,需要一些重要的引理.

引理1在(H1)–(H3)條件下,對(duì)任意p ≥2,存在一個(gè)獨(dú)立于?的常數(shù)M>0,使得

證不失一般性,假定x(t)是有界的,否則的話,對(duì)每個(gè)整數(shù)n,定義停時(shí)τn=inf{t ∈[0,T]:|x(t)|≥n},考慮停止過(guò)程x(t ∨τn)即可.由基本不等式、假設(shè)(H3)及不等式可得

因此對(duì)任意的t1∈[0,T],有

顯然有

由BDG 不等式及假設(shè)(H2)可得

其中Cp為與p 有關(guān)的正的常數(shù).對(duì)于跳部分,由假設(shè)(H2)及文[3]引理3.2,得

其中Dp為正的常數(shù).注意到對(duì)任意的t1∈[0,T],有

將(3.2)–(3.4)式代入(3.5)式得

因此

從而定理得證.

下面先建立兩個(gè)停時(shí),

引理2在假設(shè)(H1)–(H3)下,有

這里C1是獨(dú)立于?的正的常數(shù).

證對(duì)任意的t1∈[0,T],

證明方法與引理1 的方法相同,這里其證明省略.

推論3在假設(shè)(H1)–(H3)下有

這里的M2是一個(gè)正的常數(shù)且獨(dú)立于?.

證對(duì)任意的t ∈[0,T ∧ρd],存在k 使得t ∈[tk,tk+1),注意到

因此

又

因此

由基本不等式及假設(shè)(H2)可得

由推論3,文[2]引理3.2 及文[9]中Lyapunov 不等式得

取M3=3k1T[1+r+2π(Z)](1+2M2)即可,其中r 為布朗運(yùn)動(dòng)的維數(shù).

定理1 的證明假設(shè)e(t)=x(t)?y(t),易知

因此對(duì)任意的δ>0,有

根據(jù)引理1,有

由基本不等式得

將(3.13),(3.14)式代入(3.12)式得

根據(jù)x(t)和y(t)的定義,有

類(lèi)似于引理4 中的證明可知

不妨令M4=3k1[1+2M2](?+m+2π(Z)),則

由假設(shè)(H3)及(3.16)式知

由假設(shè)(H1)及H?lder 不等式可得

由假設(shè)(H1)及BDG 不等式可得

由假設(shè)(H1)及文[2]引理3.2 可得

因此

設(shè)

將(3.15)、(3.18)式代入(3.11)式得