郭嘉瑋,王同科
( 天津師范大學(xué)數(shù)學(xué)科學(xué)學(xué)院,天津300387)
考慮如下形式的第二類Fredholm積分方程
其中,0是參數(shù),f(x)、k(x,y)為已知函數(shù),u(x)為待求函數(shù),k(x,y)稱為積分方程(1.1)的核函數(shù).假設(shè)核函數(shù)k(x,y)與右端項(xiàng)f(x)具有滿足需要的最低階光滑性使得積分方程(1.1)存在惟一的連續(xù)解[1].
退化核法[2?7]是一種求解積分方程的簡單方法.通常而言,退化核的構(gòu)造是基于Taylor公式和Lagrange插值導(dǎo)出的.但如果在求解區(qū)間上,核函數(shù)有一個(gè)或多個(gè)奇點(diǎn),由于函數(shù)在奇點(diǎn)處不存在通常意義上的Taylor公式,且多項(xiàng)式插值逼近精度低,所以,此種情形的退化核方法計(jì)算精度顯著下降.文[8]針對(duì)單點(diǎn)奇異核函數(shù)構(gòu)造了一種分?jǐn)?shù)階退化核方法.本文將針對(duì)在區(qū)間兩個(gè)端點(diǎn)處非充分光滑的核函數(shù),設(shè)計(jì)基于分段混合線性插值的退化核方法.
假定在積分方程(1.1)中,核函數(shù)k(x,y)在積分區(qū)間的兩個(gè)端點(diǎn)均奇異,可以使用核函數(shù)在奇點(diǎn)處的分?jǐn)?shù)階Taylor級(jí)數(shù)來近似.與整數(shù)階Taylor級(jí)數(shù)性質(zhì)類似,分?jǐn)?shù)階Taylor級(jí)數(shù)[9?10]也僅在奇點(diǎn)處有較高的精度,在遠(yuǎn)離奇點(diǎn)的時(shí)候精度逐漸降低,故本文利用分段混合插值構(gòu)造近似退化核,在包含奇點(diǎn)的小區(qū)間上使用分?jǐn)?shù)階Taylor公式,在其它區(qū)間上使用標(biāo)準(zhǔn)的分段線性插值來逼近核函數(shù),由此得到一種分段混合插值退化核方法.數(shù)值結(jié)果表明對(duì)于兩端奇異的積分方程,傳統(tǒng)方法如Nystrm方法[1,7]收斂階下降,本文所提出的插值方法仍保持二階逼近精度.
考慮積分方程(1.1),假設(shè)其核函數(shù)k(x,y)在y=a和y=b處奇異且成立分?jǐn)?shù)階Taylor展開式
將積分區(qū)間[a,b]劃分為n個(gè)子區(qū)間,第i個(gè)子區(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),取有限項(xiàng)如下
由于函數(shù)的分?jǐn)?shù)階Taylor展開式僅在奇點(diǎn)附近具有較高的逼近程度,因此在其它區(qū)間利用插值方法來逼近.為了保證逼近函數(shù)的整體連續(xù)性,需要對(duì)(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)確定,計(jì)算可得
在下文中為方便推導(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)兩邊同時(shí)乘以(x ?a)αυ,υ=1,2,··· ,ma,并在區(qū)間[a0,a1]上對(duì)x積分,得
其中Bi,j表示Bi的第j個(gè)分量.
將(2.15)寫成向量形式,即
其中
其次,將方程(2.13)兩邊同時(shí)乘以lj,γ(x),γ=1,2,并在區(qū)間[aj?1,aj],j=2,3,··· ,n?1上對(duì)x積分,得
將其寫成向量形式
最后,將方程(2.13)兩邊同時(shí)乘以(b ?x)βω,ω=1,2,··· ,mb,并在區(qū)間[an?1,an]上對(duì)x積分,得
將其寫成向量形式
其中
聯(lián)立(2.16)、(2.18)和(2.20),并令
則積分方程(1.1)離散為一個(gè)線性代數(shù)方程組
需要指出的是,為使(2.21)的解存在唯一,需要假定λ不是系數(shù)矩陣A的特征值.
記Ln(x,y)為上節(jié)構(gòu)造的混合插值退化核,離散后的方程(2.13)寫成算子方程的形式
其中I為單位算子.為分析算法的收斂性,下面給出一些引理.
引理1[1]記
根據(jù)該引理,為說明算法的收斂性,需要證明當(dāng)h →0時(shí),ρn→0.經(jīng)分析,ρn由三部分構(gòu)成,分別是線性插值的誤差以及函數(shù)兩個(gè)奇點(diǎn)所在區(qū)間利用分?jǐn)?shù)階Taylor展開而產(chǎn)生的誤差.考慮后兩種誤差,對(duì)于分?jǐn)?shù)階插值公式(2.10)和(2.11),它們的插值余項(xiàng)分別為[11]
其中
由此可得
進(jìn)一步,由(2.26)以及線性插值誤差估計(jì)式,得
其中
定理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時(shí),由分?jǐn)?shù)階混合線性插值退化核方法所得到的近似解un(x) 一致收斂于精確解u(x),而且收斂階為min(αma+1,βmb+1,2).
需要指出的是,當(dāng)函數(shù)接近奇異時(shí),普通插值的精度將顯著下降,所以,在實(shí)際應(yīng)用本文算法時(shí),步長h1和hn要適當(dāng)增大,相應(yīng)地,核函數(shù)k(x,y) 在x=a和x=b的Taylor 展開式也應(yīng)該多取幾項(xiàng),以保證算法在包含奇點(diǎn)的小區(qū)間上的誤差小于總體計(jì)算誤差.
例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展開式中
此時(shí)
當(dāng)h1<2時(shí),收斂.同理當(dāng)hn<2時(shí),也收斂,則由定理1可得,當(dāng)h →0 時(shí),所得到的近似解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時(shí),u(x)?un(x)的圖像,此時(shí)un(x)=0.999970534x2.
計(jì)算最大誤差Eh=∥u(x)?un(x)∥∞及數(shù)值收斂階結(jié)果如表1所示.
表1 例1計(jì)算結(jié)果
圖1 例1當(dāng)h=0.0125時(shí)的誤差圖形
另外,在表1中Condh指系數(shù)矩陣A的2-范數(shù)條件數(shù),不同的h所對(duì)應(yīng)的Condh雖然隨著h的減小而增加,但增長相對(duì)緩慢,并沒有出現(xiàn)數(shù)量級(jí)的變化,說明本文所提出的方法是良性的,表1中的結(jié)果顯示數(shù)值收斂階近似為2,與理論收斂階一致.需要指出的是,本例中的核函數(shù)k(x,y)在積分區(qū)間端點(diǎn)處函數(shù)值為∞,不能直接采用線性插值的方法來構(gòu)造退化核.
例2用混合線性插值法求解如下的第二類Fredholm積分方程
其中
解在這個(gè)例子中,根據(jù)(2.1)和(2.2),核函數(shù)k(x,y)在y=0和y=1處的分?jǐn)?shù)階Taylor展開式為
其中Γ(x)表示gamma函數(shù),此時(shí)
當(dāng)h1<1時(shí),級(jí)數(shù)收斂.同樣地,當(dāng)hn<1時(shí)也收斂.由定理1可得,當(dāng)h →0時(shí),近似解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).計(jì)算最大誤差Eh及數(shù)值收斂階Oh,結(jié)果如表2所示.
表2 例2混合插值方法計(jì)算結(jié)果
圖2 例2當(dāng)h=0.0125時(shí)的誤差圖形
由表2可以看出,在例2中系數(shù)矩陣條件數(shù)變化很小,說明算法始終是良性的,誤差的收斂階近似等于2,與理論分析相吻合,說明本文所提出的算法是正確且有效的.
圖2是當(dāng)h=0.0125時(shí)的誤差圖形.由該圖形可以看到,誤差u(x)?un(x)始終穩(wěn)定在10?5數(shù)量級(jí),說明混合插值算法對(duì)端點(diǎn)奇異核函數(shù)的處理是成功的.
此外,本例中的積分方程可以在全區(qū)間上直接利用插值構(gòu)造退化核或采用Nystrm方法來求解.首先,對(duì)本例中的核函數(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ū)間上的分段線性插值計(jì)算結(jié)果
表4 例2Nystrm方法計(jì)算結(jié)果
表4 例2Nystrm方法計(jì)算結(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í),這兩種方法的計(jì)算精度很差.此算例說明對(duì)于非光滑核函數(shù),本文構(gòu)造的分?jǐn)?shù)階混合插值退化核方法是成功的.