亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

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

        2019-06-27 09:58:58郭嘉瑋王同科
        應(yīng)用數(shù)學(xué) 2019年3期
        關(guān)鍵詞:展開式奇點(diǎn)端點(diǎn)

        郭嘉瑋,王同科

        ( 天津師范大學(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)與右端項(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]收斂階下降,本文所提出的插值方法仍保持二階逼近精度.

        2.混合線性插值方法

        考慮積分方程(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ì)算誤差.

        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展開式中

        此時(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ù)階混合插值退化核方法是成功的.

        猜你喜歡
        展開式奇點(diǎn)端點(diǎn)
        例說二項(xiàng)展開式的通項(xiàng)公式應(yīng)用
        非特征端點(diǎn)條件下PM函數(shù)的迭代根
        校中有笑
        校中有笑
        校中有笑
        泰勒展開式在函數(shù)中的應(yīng)用
        奇點(diǎn)迷光(上)
        軍事文摘(2020年14期)2020-12-17 06:27:46
        不等式求解過程中端點(diǎn)的確定
        參數(shù)型Marcinkiewicz積分算子及其交換子的加權(quán)端點(diǎn)估計(jì)
        基丁能雖匹配延拓法LMD端點(diǎn)效應(yīng)處理
        少妇无套裸按摩呻吟无呜 | 一本大道久久东京热无码av| 亚洲一区二区三区成人在线| 亚洲av毛片在线网站| 大地资源网高清在线播放| 伊人久久综合精品无码av专区 | 日本不卡一区二区三区在线| 亚洲av日韩一卡二卡| 国产日产欧洲系列| 亚洲精品老司机在线观看| 亚洲 美腿 欧美 偷拍| 草逼视频污的网站免费| 欧美xxxx做受欧美| 欧美色aⅴ欧美综合色| 亚洲青青草视频在线播放| 国产精品人伦一区二区三| 中国老熟女重囗味hdxx| 狠狠色狠狠色综合| 一区二区日本影院在线观看| 蜜桃视频在线免费观看| 成人久久久久久久久久久| 手机看片国产日韩| 日韩人妻一区二区中文字幕| 中文亚洲av片不卡在线观看| 台湾佬娱乐中文22vvvv| 国产福利美女小视频| 丝袜美腿国产一区二区| 国产乱人激情h在线观看| 国产91成人精品亚洲精品| 精品少妇白浆一二三区| 久久久99精品免费视频| a级毛片无码免费真人| 日本专区一区二区三区| 丝袜美腿精品福利在线视频| 亚洲日韩精品一区二区三区| 精品国产成人亚洲午夜福利| 亚洲国产精一区二区三区性色| 一区二区三区人妻少妇| 久久久精品人妻一区二区三区蜜桃 | 免费va国产高清大片在线| 亚洲区小说区图片区|