第二類兩端奇異Fredholm積分方程的分?jǐn)?shù)階線性插值方法

郭嘉瑋,王同科

( 天津師范大學(xué)數(shù)學(xué)科學(xué)學(xué)院,天津300387)

1.引言

考慮如下形式的第二類Fredholm積分方程

其中,0是參數(shù),f(x)、k(x,y)為已知函數(shù),u(x)為待求函數(shù),k(x,y)稱為積分方程(1.1)的核函數(shù).假設(shè)核函數(shù)k(x,y)與右端項f(x)具有滿足需要的最低階光滑性使得積分方程(1.1)存在惟一的連續(xù)解[1].

退化核法[2?7]是一種求解積分方程的簡單方法.通常而言,退化核的構(gòu)造是基于Taylor公式和Lagrange插值導(dǎo)出的.但如果在求解區(qū)間上,核函數(shù)有一個或多個奇點(diǎn),由于函數(shù)在奇點(diǎn)處不存在通常意義上的Taylor公式,且多項式插值逼近精度低,所以,此種情形的退化核方法計算精度顯著下降.文[8]針對單點(diǎn)奇異核函數(shù)構(gòu)造了一種分?jǐn)?shù)階退化核方法.本文將針對在區(qū)間兩個端點(diǎn)處非充分光滑的核函數(shù),設(shè)計基于分段混合線性插值的退化核方法.

假定在積分方程(1.1)中,核函數(shù)k(x,y)在積分區(qū)間的兩個端點(diǎn)均奇異,可以使用核函數(shù)在奇點(diǎn)處的分?jǐn)?shù)階Taylor級數(shù)來近似.與整數(shù)階Taylor級數(shù)性質(zhì)類似,分?jǐn)?shù)階Taylor級數(shù)[9?10]也僅在奇點(diǎn)處有較高的精度,在遠(yuǎn)離奇點(diǎn)的時候精度逐漸降低,故本文利用分段混合插值構(gòu)造近似退化核,在包含奇點(diǎn)的小區(qū)間上使用分?jǐn)?shù)階Taylor公式,在其它區(qū)間上使用標(biāo)準(zhǔn)的分段線性插值來逼近核函數(shù),由此得到一種分段混合插值退化核方法.數(shù)值結(jié)果表明對于兩端奇異的積分方程,傳統(tǒng)方法如Nystrm方法[1,7]收斂階下降,本文所提出的插值方法仍保持二階逼近精度.

2.混合線性插值方法

考慮積分方程(1.1),假設(shè)其核函數(shù)k(x,y)在y=a和y=b處奇異且成立分?jǐn)?shù)階Taylor展開式

將積分區(qū)間[a,b]劃分為n個子區(qū)間,第i個子區(qū)間的長度為hi,且記a0=a,ai=ai?1+hi,i=1,2,··· ,n.由積分區(qū)間的可加性,積分方程(1.1)則為

下面構(gòu)造k(x,y)基于分?jǐn)?shù)階Tylor展開和分段線性插值的離散退化核.

在區(qū)間[a0,a1]和[an?1,an]上,核函數(shù)k(x,y)分別在區(qū)間的左右端點(diǎn)處具有分?jǐn)?shù)階Taylor展開式(2.1)和(2.2),取有限項如下

由于函數(shù)的分?jǐn)?shù)階Taylor展開式僅在奇點(diǎn)附近具有較高的逼近程度,因此在其它區(qū)間利用插值方法來逼近.為了保證逼近函數(shù)的整體連續(xù)性,需要對(2.4)和(2.5)進(jìn)行修正.令

其中κ1(x),κ2(x)為待定函數(shù),分別由插值條件k1(x,a1)=k(x,a1),kn(x,an?1)=k(x,an?1)確定,計算可得

在下文中為方便推導(dǎo),仍用ξma(x)、ηmb(x)代替κ1(x)、κ2(x),即將(2.6)和(2.7)分別記為

在區(qū)間[ai?1,ai],i=2,3··· ,n ?1上,做核函數(shù)k(x,y)關(guān)于自變量y的線性插值,有

其中

用上面構(gòu)造的分段函數(shù)ki(x,y),i=1,2···n近似代替原方程中的核函數(shù)k(x,y),并記得到的近似解為un(x),則un(x)滿足

為了方便推導(dǎo),將上式改寫為向量形式,即

其中

顯然Bi,i=1,2,··· ,n為未知向量,下面給出確定這組未知向量的方法.

首先,方程(2.13)兩邊同時乘以(x ?a)αυ,υ=1,2,··· ,ma,并在區(qū)間[a0,a1]上對x積分,得

其中Bi,j表示Bi的第j個分量.

將(2.15)寫成向量形式,即

其中

其次,將方程(2.13)兩邊同時乘以lj,γ(x),γ=1,2,并在區(qū)間[aj?1,aj],j=2,3,··· ,n?1上對x積分,得

將其寫成向量形式

最后,將方程(2.13)兩邊同時乘以(b ?x)βω,ω=1,2,··· ,mb,并在區(qū)間[an?1,an]上對x積分,得

將其寫成向量形式

其中

聯(lián)立(2.16)、(2.18)和(2.20),并令

則積分方程(1.1)離散為一個線性代數(shù)方程組

需要指出的是,為使(2.21)的解存在唯一,需要假定λ不是系數(shù)矩陣A的特征值.

記Ln(x,y)為上節(jié)構(gòu)造的混合插值退化核,離散后的方程(2.13)寫成算子方程的形式

其中I為單位算子.為分析算法的收斂性,下面給出一些引理.

引理1[1]記

根據(jù)該引理,為說明算法的收斂性,需要證明當(dāng)h →0時,ρn→0.經(jīng)分析,ρn由三部分構(gòu)成,分別是線性插值的誤差以及函數(shù)兩個奇點(diǎn)所在區(qū)間利用分?jǐn)?shù)階Taylor展開而產(chǎn)生的誤差.考慮后兩種誤差,對于分?jǐn)?shù)階插值公式(2.10)和(2.11),它們的插值余項分別為[11]

其中

由此可得

進(jìn)一步,由(2.26)以及線性插值誤差估計式,得

其中

定理1當(dāng)核函數(shù)k(x,y)關(guān)于x在[a,b]上連續(xù),關(guān)于y在(a,b) 上二階可導(dǎo),且在y=a和y=b點(diǎn)存在分?jǐn)?shù)階Taylor展開式,則當(dāng)h →0時,由分?jǐn)?shù)階混合線性插值退化核方法所得到的近似解un(x) 一致收斂于精確解u(x),而且收斂階為min(αma+1,βmb+1,2).

需要指出的是,當(dāng)函數(shù)接近奇異時,普通插值的精度將顯著下降,所以,在實際應(yīng)用本文算法時,步長h1和hn要適當(dāng)增大,相應(yīng)地,核函數(shù)k(x,y) 在x=a和x=b的Taylor 展開式也應(yīng)該多取幾項,以保證算法在包含奇點(diǎn)的小區(qū)間上的誤差小于總體計算誤差.

3.數(shù)值算例

例1用混合線性插值法求解如下的第二類Fredholm積分方程

其中

解該題中核函數(shù)k(x,y)本身就是退化核,方程可化為

求得c=6π,再將c帶回原方程可得精確解為u(x)=(6πx2?5πx2)/π=x2.

下面采用混合線性插值法求解該問題.根據(jù)(2.1)和(2.2),本例中核函數(shù)k(x,y)在y=0和y=4處的分?jǐn)?shù)階Taylor展開式中

此時

當(dāng)h1<2時,收斂.同理當(dāng)hn<2時,也收斂,則由定理1可得,當(dāng)h →0 時,所得到的近似解un(x)一致收斂于精確解u(x).

由于核函數(shù)k(x,y)越靠近積分區(qū)間的端點(diǎn),函數(shù)的奇性越強(qiáng),因此區(qū)間端點(diǎn)所在的區(qū)間長度不宜過小.將積分區(qū)間端點(diǎn)所在的小區(qū)間長度固定為0.2,其余部分分別按照步長h為0.2、0.1、0.05、0.025、0.0125均勻劃分,并按照第二節(jié)給出的方法,取分別得到近似解un(x).圖1是當(dāng)h=0.0125時,u(x)?un(x)的圖像,此時un(x)=0.999970534x2.

計算最大誤差Eh=∥u(x)?un(x)∥∞及數(shù)值收斂階結(jié)果如表1所示.

表1 例1計算結(jié)果

圖1 例1當(dāng)h=0.0125時的誤差圖形

另外,在表1中Condh指系數(shù)矩陣A的2-范數(shù)條件數(shù),不同的h所對應(yīng)的Condh雖然隨著h的減小而增加,但增長相對緩慢,并沒有出現(xiàn)數(shù)量級的變化,說明本文所提出的方法是良性的,表1中的結(jié)果顯示數(shù)值收斂階近似為2,與理論收斂階一致.需要指出的是,本例中的核函數(shù)k(x,y)在積分區(qū)間端點(diǎn)處函數(shù)值為∞,不能直接采用線性插值的方法來構(gòu)造退化核.

例2用混合線性插值法求解如下的第二類Fredholm積分方程

其中

解在這個例子中,根據(jù)(2.1)和(2.2),核函數(shù)k(x,y)在y=0和y=1處的分?jǐn)?shù)階Taylor展開式為

其中Γ(x)表示gamma函數(shù),此時

當(dāng)h1<1時,級數(shù)收斂.同樣地,當(dāng)hn<1時也收斂.由定理1可得,當(dāng)h →0時,近似解un(x)一致收斂于精確解u(x).

與例1的求解思路類似,將積分區(qū)間端點(diǎn)所在的小區(qū)間長度固定為0.2,其余部分分別按照步長h為0.2、0.1、0.05、0.025、0.0125均勻劃分,并按照第二節(jié)方法取分別得到近似解un(x).計算最大誤差Eh及數(shù)值收斂階Oh,結(jié)果如表2所示.

表2 例2混合插值方法計算結(jié)果

圖2 例2當(dāng)h=0.0125時的誤差圖形

由表2可以看出,在例2中系數(shù)矩陣條件數(shù)變化很小,說明算法始終是良性的,誤差的收斂階近似等于2,與理論分析相吻合,說明本文所提出的算法是正確且有效的.

圖2是當(dāng)h=0.0125時的誤差圖形.由該圖形可以看到,誤差u(x)?un(x)始終穩(wěn)定在10?5數(shù)量級,說明混合插值算法對端點(diǎn)奇異核函數(shù)的處理是成功的.

此外,本例中的積分方程可以在全區(qū)間上直接利用插值構(gòu)造退化核或采用Nystrm方法來求解.首先,對本例中的核函數(shù)k(x,y) 在積分區(qū)間[0,1] 上直接進(jìn)行分段線性插值,取步長h分別為0.2、0.1、0.05、0.025、0.0125可得近似解、最大誤差以及數(shù)值收斂階,如表3所示.

表3 例2全區(qū)間上的分段線性插值計算結(jié)果

表4 例2Nystrm方法計算結(jié)果

表4 例2Nystrm方法計算結(jié)果

h Eh Oh Condh 0.2 9.13371×10?2 1.73324 0.1 4.07416×10?2 1.16470 1.83738 0.05 1.86295×10?2 1.12891 1.88729 0.025 9.48273×10?3 0.97421 1.90905 0.0125 5.76217×10?3 0.71869 1.91836

由表3和表4可以看出,傳統(tǒng)的基于插值構(gòu)造的退化核方法以及Nystrm方法的收斂階均依賴核函數(shù)以及解的光滑性,當(dāng)核函數(shù)與解存在奇點(diǎn)時,這兩種方法的計算精度很差.此算例說明對于非光滑核函數(shù),本文構(gòu)造的分?jǐn)?shù)階混合插值退化核方法是成功的.