帶有Caputo-Fabrizio分?jǐn)?shù)階導(dǎo)數(shù)的Fokker-Planck方程的算法研究
發(fā)布時(shí)間:2021-08-11 02:57
分?jǐn)?shù)階微積分作為整數(shù)階微積分推廣而出現(xiàn),一方面其方程在一定程度上推動(dòng)了分?jǐn)?shù)階微積分理論的發(fā)展,另一方面引起學(xué)者們的重視并被廣泛應(yīng)用于各大領(lǐng)域.比如海洋動(dòng)力學(xué)、流體動(dòng)力學(xué)、超導(dǎo)和經(jīng)濟(jì)金融等領(lǐng)域.雖然,分?jǐn)?shù)階微分方程能準(zhǔn)確地解釋一些數(shù)學(xué)物理領(lǐng)域內(nèi)的非線性問(wèn)題,但是一般情況下,分?jǐn)?shù)階微分方程的解析解較難得到,從而數(shù)值解的獲取理論意義和應(yīng)用價(jià)值顯得尤為重要.本文對(duì)帶有Caputo-Fabrizio分?jǐn)?shù)階導(dǎo)數(shù)的Fokker-Planck方程的兩種數(shù)值解法進(jìn)行探究.首先考慮用Ritz近似解決帶有Caputo-Fabrizio分?jǐn)?shù)階導(dǎo)數(shù)的線性Fokker-Planck方程,并使用多項(xiàng)式函數(shù)得到代數(shù)方程組,然后求解了線性代數(shù)方程組,最后通過(guò)算例說(shuō)明了該方法的有效性和適用性.其次,基于再生核理論,對(duì)Caputo-Fabrizio分?jǐn)?shù)階導(dǎo)數(shù)的Fokker-Planck方程進(jìn)行數(shù)值求解,構(gòu)造二維再生核空間,將問(wèn)題轉(zhuǎn)化為求解非線性方程組,這種求解方法避免使用Gram-Schmidt正交化過(guò)程.并且提高計(jì)算的速度,最后通過(guò)數(shù)值算例說(shuō)明了方法良好有效.
【文章來(lái)源】:哈爾濱師范大學(xué)黑龍江省
【文章頁(yè)數(shù)】:41 頁(yè)
【學(xué)位級(jí)別】:碩士
【部分圖文】:
圖3.2:數(shù)值解與解析解比較??-[42]
—4?8.08258x10一5?4.30474x10—5?2.81919x10—5??0.7?1.64276xl0—3?8.77138x10—4?5.9744x10—4?4.52876x10—4??0.8?3.27658x10—3?1,70738x10—3?1.15397x10—3?8.71425xl0—4??0.9?4.33614x10—3?2.25792x10—3?L52567x?10—3?1.15196x10—3??。海?=?0.5,?a?=?0.5,0.7,0.9,1數(shù)值結(jié)果如圖4-1所示??2????_?|?|?—0—?a?=?0.5? ̄A ̄?a?=?0.7??■?—〇—?a?=?0.9?—D—?a?=?1.0??丁?—丨1?■??I?I?I?I?I??0?—??i?r ̄n?i?i?i?i?i?i?i?i?i?i?i?i?i?i?i?i?i?|??0?0.2?0.4?0.6?0.8?1??圖4.1:?rr?=?0.5時(shí)的數(shù)值解結(jié)果比較??取a?=?1時(shí)數(shù)值解與解析解的三維曲面如圖4-2所示??-22-??
第4掌再生核法求解非淺性Caputo-Fabrizio分?jǐn)?shù)階導(dǎo)數(shù)的RAke&Planck方程??_??1^)??□?數(shù)值解?a?解析解??圖4.2:。?=?1.0時(shí)的數(shù)值解與解析解三維曲面比較??4.4.2算例二??取乂(M,=?f?—?f,B(M,?q?R考慮如下算例??^?9?_?f>u{x^)\?+?-§^U2{^,?t),?〇?e?(0,?l],x,t?e?[0,1]?(4_24)??U{x^?0)?=?x2??當(dāng)a?=?1時(shí),我們可以得到精確解??由于方程在a不是整數(shù)時(shí),解析解不易求出,所以為了比較誤差,如表4-2,我??們?nèi)。?=?H,=M,=||,=M時(shí)的數(shù)值解來(lái)與〇■?=?i時(shí)的解析解來(lái)比較誤墓??—?23?—??
【參考文獻(xiàn)】:
期刊論文
[1]關(guān)于邊值問(wèn)題的一種新的再生核數(shù)值算法[J]. 朱慧,林迎珍. 數(shù)學(xué)的實(shí)踐與認(rèn)識(shí). 2015(23)
[2]三步五階迭代方法解非線性方程組[J]. 張旭,檀結(jié)慶. 計(jì)算數(shù)學(xué). 2013(03)
[3]W22(D)空間函數(shù)的近似再生[J]. 吳勃英. 哈爾濱師范大學(xué)自然科學(xué)學(xué)報(bào). 1995(03)
[4]W22[a,b]空間中的最佳Hermite插值算子[J]. 吳勃英,崔明根,鄧中興. 計(jì)算物理. 1988(02)
[5]Solutions to the Definite Solution Problem of Differential Equations in Space W2l[J]. 崔明根,鄧中興,張少謙. 數(shù)學(xué)進(jìn)展. 1988(03)
[6]W21空間中的最佳插值逼近算子[J]. 崔明根,鄧中興. 計(jì)算數(shù)學(xué). 1986(02)
本文編號(hào):3335304
【文章來(lái)源】:哈爾濱師范大學(xué)黑龍江省
【文章頁(yè)數(shù)】:41 頁(yè)
【學(xué)位級(jí)別】:碩士
【部分圖文】:
圖3.2:數(shù)值解與解析解比較??-[42]
—4?8.08258x10一5?4.30474x10—5?2.81919x10—5??0.7?1.64276xl0—3?8.77138x10—4?5.9744x10—4?4.52876x10—4??0.8?3.27658x10—3?1,70738x10—3?1.15397x10—3?8.71425xl0—4??0.9?4.33614x10—3?2.25792x10—3?L52567x?10—3?1.15196x10—3??。海?=?0.5,?a?=?0.5,0.7,0.9,1數(shù)值結(jié)果如圖4-1所示??2????_?|?|?—0—?a?=?0.5? ̄A ̄?a?=?0.7??■?—〇—?a?=?0.9?—D—?a?=?1.0??丁?—丨1?■??I?I?I?I?I??0?—??i?r ̄n?i?i?i?i?i?i?i?i?i?i?i?i?i?i?i?i?i?|??0?0.2?0.4?0.6?0.8?1??圖4.1:?rr?=?0.5時(shí)的數(shù)值解結(jié)果比較??取a?=?1時(shí)數(shù)值解與解析解的三維曲面如圖4-2所示??-22-??
第4掌再生核法求解非淺性Caputo-Fabrizio分?jǐn)?shù)階導(dǎo)數(shù)的RAke&Planck方程??_??1^)??□?數(shù)值解?a?解析解??圖4.2:。?=?1.0時(shí)的數(shù)值解與解析解三維曲面比較??4.4.2算例二??取乂(M,=?f?—?f,B(M,?q?R考慮如下算例??^?9?_?f>u{x^)\?+?-§^U2{^,?t),?〇?e?(0,?l],x,t?e?[0,1]?(4_24)??U{x^?0)?=?x2??當(dāng)a?=?1時(shí),我們可以得到精確解??由于方程在a不是整數(shù)時(shí),解析解不易求出,所以為了比較誤差,如表4-2,我??們?nèi)。?=?H,=M,=||,=M時(shí)的數(shù)值解來(lái)與〇■?=?i時(shí)的解析解來(lái)比較誤墓??—?23?—??
【參考文獻(xiàn)】:
期刊論文
[1]關(guān)于邊值問(wèn)題的一種新的再生核數(shù)值算法[J]. 朱慧,林迎珍. 數(shù)學(xué)的實(shí)踐與認(rèn)識(shí). 2015(23)
[2]三步五階迭代方法解非線性方程組[J]. 張旭,檀結(jié)慶. 計(jì)算數(shù)學(xué). 2013(03)
[3]W22(D)空間函數(shù)的近似再生[J]. 吳勃英. 哈爾濱師范大學(xué)自然科學(xué)學(xué)報(bào). 1995(03)
[4]W22[a,b]空間中的最佳Hermite插值算子[J]. 吳勃英,崔明根,鄧中興. 計(jì)算物理. 1988(02)
[5]Solutions to the Definite Solution Problem of Differential Equations in Space W2l[J]. 崔明根,鄧中興,張少謙. 數(shù)學(xué)進(jìn)展. 1988(03)
[6]W21空間中的最佳插值逼近算子[J]. 崔明根,鄧中興. 計(jì)算數(shù)學(xué). 1986(02)
本文編號(hào):3335304
本文鏈接:http://www.sikaile.net/kejilunwen/yysx/3335304.html
最近更新
教材專著