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

        ?

        全直線上四階方程的Laguerre-Laguerre復(fù)合譜逼近

        2021-04-15 01:20:10蔡耀雄莊清渠
        關(guān)鍵詞:代數(shù)方程四階對角

        蔡耀雄, 莊清渠

        (華僑大學(xué) 數(shù)學(xué)科學(xué)學(xué)院, 福建 泉州362021)

        近三十年來,關(guān)于求解四階微分方程的譜方法的研究有很大進展.Shen[1]采用Legendre-Galerkin譜方法求解二階和四階微分方程;Kwan等[2]考慮可分離橢圓型問題的并行譜元計算;Shen等[3]采用Legendre Petrov-Galerkin逼近法數(shù)值求解四階方程;文獻[4-5]分別研究一維和二維四階方程的譜元計算;莊清渠等[6]研究四階常微分方程的Birkhoff配點法;文獻[7-8]分別研究兩類四階積分微分方程的譜逼近;Chen[9]研究四階Cahn-Hilliard方程的Legendre-Galerkin譜逼近.上述文獻研究的是有界區(qū)域上的四階方程,對于無界區(qū)域上的四階方程也有一些研究工作.葉小華[10]研究一維半直線區(qū)域四階方程的Legendre-Laguerre復(fù)合譜方法;Zhuang等[11]研究一維半無界區(qū)域四階方程的Legendre-Laguerre耦合譜元計算;李敏等[12]研究半無界條狀區(qū)域四階方程的Laguerre-Legendre混合譜逼近;李珊等[13]研究半直線區(qū)域上四階橢圓型方程的有理Legendre函數(shù)全對角化譜方法;Yu等[14]研究全直線區(qū)域上的對角化Legendre有理譜方法.本文用Laguerre-Laguerre復(fù)合譜方法求解全直線上的四階方程,并與一類Hermite譜方法進行對比,主要考慮方程的數(shù)值計算.

        1 問題及復(fù)合逼近形式

        記I:=(-∞,∞),考慮如下的四階問題,即

        (1)

        d(u,v)=(f,v), ?v∈V,

        (2)

        式(2)中:d(u,v):=λ2(u,v)+λ1(ux,vx)+(uxx,vxx),?u,v∈V.

        對問題(2),用Laguerre-Laguerre復(fù)合譜方法進行求解.首先,將(-∞,∞)剖分成I1=(-∞,0),I2=(0,∞)兩部分.然后,在兩個區(qū)間上分別采用Laguerre譜方法進行逼近.

        記uI1:=u|I1,uI2:=u|I2,N=(N1,N2) ,令PN(Ω) 為Ω上次數(shù)不超過N的全體多項式組成的空間,并記

        (3)

        (4)

        則問題(1)的Laguerre-Laguerre復(fù)合逼近形式如下:找uN∈VN,使得

        d(uN,vN)=(f,vN), ?vN∈VN.

        (5)

        2 計算實施

        Laguerre多項式Lm的基本性質(zhì)[15]如下,即

        (6)

        (7)

        (8)

        (9)

        (10)

        (11)

        類似于文獻[11,16]求解高階方程的過程,Laguerre-Laguerre復(fù)合逼近問題(5)所對應(yīng)的線性系統(tǒng)可通過如下3個步驟進行求解.

        (12)

        (13)

        (14)

        3) 結(jié)合.由式(13),(14)可知,對任意的vN∈VN,有

        (15)

        因此,問題(5)的解為

        (16)

        問題(5)分解成了兩個相對獨立的子問題,子問題(12)和(13)分別由九對角問題和五對角問題組成,因而容易進行求解.

        (λ2C+λ1B+A)U2=F2.

        (17)

        由Laguerre函數(shù)的正交性可得

        由此可知,矩陣都是對稱的,且最多是五對角陣,因而式(17)可以有效地求解.子問題(13)在I1上的代數(shù)方程組與其在I2上的代數(shù)方程組類似.

        子問題(12)對應(yīng)的代數(shù)方程組的系數(shù)矩陣與子問題(13)相同.子問題(12)對應(yīng)的代數(shù)方程組的右端項可通過對φ1(x),φ2(x)進行求導(dǎo),有

        又因為

        從而在區(qū)間I1上,有

        以及

        而在區(qū)間I2上,有

        以及

        因此,可得子問題(12)對應(yīng)代數(shù)方程組的右端項.

        對于子問題(14),由于

        圖1 最大誤差隨的變化情況Fig.1 Maximum errors change with

        3 數(shù)值算例

        算例1在問題(1)中,固定λ1=λ2=1,取精確解為

        u(x)=sin(x)e-σx2.

        其次,當(dāng)σ=0.01和σ=0.02時,分別利用Laguerre-Laguerre復(fù)合譜方法及Hermite譜方法進行計算,得到的最大誤差,如表1所示.由表1可知:當(dāng)σ比較小時,用Laguerre-Laguerre復(fù)合譜方法進行逼近的誤差比用Hermite譜方法進行逼近的誤差要小得多.因此,Laguerre-Laguerre復(fù)合譜方法逼近更具優(yōu)越性.

        表1 兩種方法在u(x)=sin(x)e-σx2時的最大誤差Tab.1 Maximum errors by two ways for u(x)=sin(x)e-σx2

        表2 兩種方法在時的最大誤差Tab.2 Maximum errors by two ways for

        4 結(jié)束語

        將全直線區(qū)域剖分為兩部分,進而構(gòu)造了求解全直線上四階方程的Laguerre-Laguerre復(fù)合譜方法.數(shù)值結(jié)果表明,復(fù)合方法具有譜收斂性.同時,通過與純Hermite譜方法的比較可以看出,復(fù)合譜方法對求解具有衰減緩慢解析解的問題具有優(yōu)越性.

        猜你喜歡
        代數(shù)方程四階對角
        四階p-廣義Benney-Luke方程的初值問題
        擬對角擴張Cuntz半群的某些性質(zhì)
        基于置換思想的代數(shù)方程求解理論探析
        未知量符號x的歷史穿越
        拉格朗日代數(shù)方程求解中的置換思想
        矩陣代數(shù)方程在城市燃?xì)夤芫W(wǎng)水力計算中的應(yīng)用研究
        上海煤氣(2016年1期)2016-05-09 07:12:37
        帶參數(shù)的四階邊值問題正解的存在性
        四階累積量譜線增強方法的改進仿真研究
        基于四階累積量和簡化粒子群的盲分離算法
        非奇異塊α1對角占優(yōu)矩陣新的實用簡捷判據(jù)
        亚洲男人堂色偷偷一区 | 国产精品无码一区二区在线观一 | 国产亚洲青春草在线视频| 一区二区三区亚洲视频| 国产精品白浆在线观看免费| 野花在线无码视频在线播放 | 亚洲国产成人资源在线桃色| 日本人妻高清免费v片| 成人欧美一区二区三区在线| 国产福利酱国产一区二区| 思思久久96热在精品不卡| 久久久亚洲一区二区三区| 国产日产欧产精品精品蜜芽| 九九99久久精品国产| 中文字幕久久久人妻无码| 视频一区二区三区国产| 欧美成人猛片aaaaaaa| 亚洲av无码国产剧情| 国产成人无精品久久久| 久久精见国产亚洲av高清热| 久久精品国产亚洲av香蕉| a在线观看免费网站大全| 无码伊人久久大香线蕉| 人妖国产视频一区二区| 无码精品人妻一区二区三区av| 欧美亚洲日本在线| 国产成人精品自拍在线观看| 欧美最猛性xxxx| 99精品国产高清一区二区麻豆| 日韩精品永久免费播放平台| 在线观看国产一区二区av| 2020无码专区人妻系列日韩| 伊人久久网国产伊人| 国产丝袜美腿诱惑在线观看| 亚洲国产av自拍一区| 亚洲精品97久久中文字幕无码| 国产国拍亚洲精品永久69| 亚洲最大一区二区在线观看| 极品少妇一区二区三区四区| 国产偷国产偷高清精品| 国产成人av三级三级三级在线|