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