陳 琛, 馮曉莉, 陳漢章

(西安電子科技大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院, 西安 710126)

0 引 言

眾所周知,反問題的應(yīng)用非常廣泛,相應(yīng)的研究成果也越來越多[1-5]．然而,絕大多數(shù)反問題都是不適定的,尤其是穩(wěn)定性往往不成立,這也是研究反問題的困難和意義所在．不僅如此,實(shí)際應(yīng)用中會(huì)有很多隨機(jī)因素的影響,比如測量數(shù)據(jù)帶有噪聲,數(shù)學(xué)模型的建立含有模型誤差等．這些因素對(duì)于不適定反問題的研究不容小覷,否則求得的數(shù)值結(jié)果會(huì)出現(xiàn)災(zāi)難性的問題．一種自然而又巧妙的處理方法是將這些因素視為隨機(jī)的,相應(yīng)的問題就看作隨機(jī)問題．此外,隨機(jī)技術(shù)的引入可以很好地耦合不同度量尺度之間的干擾[6-9]．綜上可知,關(guān)于隨機(jī)反問題的研究既非常重要又很有實(shí)際意義．

本文關(guān)心的是如下一類隨機(jī)微分方程：

(1)

其中a(x),h(t)已知,且a(x)?0,h(t)≥h0>0,D=(0,1),BH(t)為概率空間(Ω,F,{Ft}t≥0,)上的分式Brown運(yùn)動(dòng),H∈(0,1)為Hurst參數(shù),為彩噪聲．為了簡便,后文將u(x,t,ω)簡記為u(x,t)．本文主要探討：

(P1) 正問題 已知φ(x),ψ(x),求溫和解u(x,t);

(P2) 反問題 給定T時(shí)刻的樣本數(shù)據(jù)u(x,T,ω), 反演問題(1)中隨機(jī)源項(xiàng)的部分統(tǒng)計(jì)量φ(x)和ψ2(x)．

對(duì)于正問題(P1),由于隨機(jī)源項(xiàng)F(x,t)的低正則性[10],該方程不是幾乎處處存在的,所以方程的適定性與經(jīng)典意義下截然不同,相應(yīng)的解在什么樣的空間存在是值得研究的．而關(guān)于反問題(P2),F(x,t)為隨機(jī)場,如何由u(x,t,ω)的統(tǒng)計(jì)量來反演F(x,t)的統(tǒng)計(jì)量φ(x)與ψ2(x)是很困難且有意義的,而其不適定性又將使得該反問題的求解難上加難．

由于以上的眾多原由,近年來反隨機(jī)源問題倍受學(xué)者們的重視[11-12],關(guān)于分式Brown運(yùn)動(dòng)驅(qū)動(dòng)下的隨機(jī)微分方程的反源問題,例如,隨機(jī)的時(shí)空分?jǐn)?shù)階方程[13],隨機(jī)波動(dòng)方程[14],一個(gè)隨機(jī)的時(shí)間分?jǐn)?shù)階擴(kuò)散方程[15],一維隨機(jī)對(duì)流擴(kuò)散方程[16]等已有一些研究結(jié)果．這些工作均先是利用Laplace 變換或Fourier變換或廣義Fourier 級(jí)數(shù)展開得到相應(yīng)的溫和解,然后再進(jìn)一步討論．本文所要討論的問題(1),是一種更為廣泛的隨機(jī)問題,比如隨機(jī)熱方程：

(2)

其中L為一致橢圓算子．借助于廣義Fourier級(jí)數(shù)展開,問題(2)可轉(zhuǎn)化為

(3)

1 理 論 基 礎(chǔ)

定義1[15,17]設(shè)0

則稱BH(t)為分式Brown運(yùn)動(dòng),其中H是Hurst參數(shù),E[·]表示期望．特別地,當(dāng)H=1/2時(shí),即為標(biāo)準(zhǔn)Brown運(yùn)動(dòng)B(t)．

對(duì)于分式Brown運(yùn)動(dòng)BH(t),其中H∈(0,1),相應(yīng)的隨機(jī)積分有如下公式[15]：

2) 當(dāng)H∈(0,1/2)時(shí),

(4)

其中

3) 當(dāng)H=1/2時(shí),

(5)

4) 當(dāng)H∈(1/2,1)時(shí),

(6)

其中αH=H(2H-1)．

2 正 問 題

不難得到問題(1)的溫和解如下．

定義2 一個(gè)隨機(jī)過程u∈L2(D)稱為問題(1)的溫和解,如果

(7)

幾乎處處成立．

下面進(jìn)一步說明在一定條件下,溫和解u(x,t)是適定的．為了推導(dǎo)的方便,先給出如下引理．

引理1[15]當(dāng)H∈(0,1/2)時(shí),對(duì)?t>τ有

由式(7)知,

2I1(t)+2I2(t)．

(8)

由h(t)≥h0>0及積分中值定理可知,存在某個(gè)ξ∈(0,t),有

其中

若令M(x)max{1,e-a(x)T},則

關(guān)于I2(t),有

情形1 當(dāng)H∈(0,1/2)時(shí),由式(4)可知

J1(t)+J2(t)+J3(t),

(9)

這里“ab”表示a≤Cb,其中C為某個(gè)正常數(shù)．顯然有

(10)

(11)

由引理1可知

J2(t)τ1-2H(t4H-2+τ4H-2)dτ·M2(x)t2HM2(x)．

關(guān)于J3(t),由于

1) 當(dāng)|a(x)|≤1時(shí),

|e-a(x)(t-u)-e-a(x)(t-τ)|≤|a(x)|M(x)|u-τ|≤M(x)|u-τ|;

2) 當(dāng)|a(x)|>1時(shí),

根據(jù)?KH(u,τ)/?u的表達(dá)式,易得

(12)

將式(10)—(12)代入式(9),可得

此時(shí)有

情形2 當(dāng)H=1/2時(shí),由式(5)可知

所以

情形3 當(dāng)H∈(1/2,1)時(shí),由式(6)可知

(13)

因此

綜上可知如下定理成立．

定理1 如果函數(shù)f(x),g(x)∈L2(D),則有估計(jì)式

注1 由定理1可知：1) 若a(x)有界,則‖M‖L∞(D)也有界;2) 若a(x)≥0,則‖M‖L∞(D)亦有界;3) 若a(x)→-∞,則‖M‖L∞(D)無界．因此第1)、2)兩種情況下問題(1)的解是穩(wěn)定的,但第3)種情況下問題(1)的解不穩(wěn)定．

3 反 問 題

本節(jié)我們考慮通過u(x,T,ω)的統(tǒng)計(jì)量來反演問題(1)隨機(jī)源項(xiàng)中的統(tǒng)計(jì)量φ(x)和ψ2(x)．若h(t)≡1時(shí),φ(x)和ψ2(x)就是隨機(jī)源項(xiàng)的期望和方差．由式(7)可知

(14)

(15)

其中var(·)表示方差．

進(jìn)一步可得

(16)

下面分析由式(14)和(15)來反演φ(x)與ψ2(x)的唯一性與穩(wěn)定性如何．

3.1 唯一性

由h(τ)≥h0>0易知,對(duì)于固定的x,總有

情形1 當(dāng)H∈(0,1/2)時(shí),

對(duì)于固定的x,當(dāng)a(x)≥0時(shí),e-a(x)(T-s)關(guān)于s是遞增的,所以

e-a(x)(T-u)-e-a(x)(T-τ)≥0．

又由KH(u,τ)的表達(dá)式可知KH(T,τ)與?KH(u,τ)/?u均大于0且連續(xù),所以一定存在一個(gè)c1(x)>0,有

當(dāng)a(x)<0時(shí),e-a(x)(T-s)關(guān)于s是遞減的,所以(d/ds)e-a(x)(T-s)<0,此時(shí)

根據(jù)分部積分,可知

由e-a(x)(T-s)關(guān)于s是遞減的,且KH(u,τ)>0,可知

(17)

情形2 當(dāng)H=1/2時(shí),由It等距公式(5)可知：當(dāng)a(x)≥0時(shí),有

當(dāng)a(x)<0時(shí),有

(18)

情形3 當(dāng)H∈(1/2,1)時(shí),有

當(dāng)a(x)≥0時(shí),有

注意到式(13),有

當(dāng)a(x)≤0,類似地有

(19)

綜上分析可知,如下定理成立．

定理2 若φ(x),ψ(x)∈L2(D)且‖ψ‖L2(D)≠0,h(x)≥h0>0,h(x)∈L∞(0,T),則源項(xiàng)中φ(x)與ψ2(x)能由u(x,T,ω)的期望和方差唯一確定．

3.2 穩(wěn)定性分析

反源問題在醫(yī)學(xué)、地質(zhì)勘探等方面都有非常廣泛的應(yīng)用,然而反源問題往往是不穩(wěn)定的．本小節(jié)將分析反演φ(x)與ψ2(x)時(shí)的穩(wěn)定性．先分析a(x)>0時(shí)的穩(wěn)定性．

一方面,我們有

情形1 若a(x)>1/T2,由式(9)可知

對(duì)于H∈(0,1/2),借助于不等式

(20)

有

根據(jù)式(20)易得

并且

因此

此外,

J21(T)+J22(T)+J23(T)．

根據(jù)引理1,有

利用證明引理1相同的方法,我們得到

同樣地,利用引理1,有

所以,對(duì)于J2(T),我們有

此外,

因此

對(duì)于H=1/2,利用It等距公式(5)并注意到a(x)>1/T2可得

對(duì)于H∈(1/2,1),利用式(6)并注意到估計(jì)式(20),將其化為4個(gè)部分,每個(gè)部分分別借用估計(jì)式(20)進(jìn)行計(jì)算,所以有

其中

由于對(duì)稱性

因此

對(duì)于H=1/2,利用It等距公式(5)易得

對(duì)于H∈(1/2,1),利用公式(6)并進(jìn)行簡單的計(jì)算可得

綜合本小節(jié)的推導(dǎo)以及3.1小節(jié)中的式(17)—(19),我們得到以下定理．

定理3 對(duì)于任意的a(x),有估計(jì)式

對(duì)于a(x)≥1/T2,有

對(duì)于a(x)<0或0

注2 由定理3可知,當(dāng)a(x)<0或0

4 結(jié) 論

本文討論了一類隨機(jī)微分方程的反源問題,根據(jù)方程的溫和解,討論了正問題的適定性．利用溫和解的統(tǒng)計(jì)性質(zhì),根據(jù)T時(shí)刻的數(shù)據(jù)反演方程的源項(xiàng),證明了源項(xiàng)反演的唯一性,并分析了關(guān)于源項(xiàng)φ(x)和ψ2(x)反演的穩(wěn)定性情況．后面我們將進(jìn)一步考慮相應(yīng)的數(shù)值結(jié)果．