近年來,隨著分?jǐn)?shù)階微積分在數(shù)學(xué)建模中的大量應(yīng)用,分?jǐn)?shù)階偏微分方程越來越受到學(xué)術(shù)界的關(guān)注,并已經(jīng)開始在電磁學(xué)、多孔介質(zhì)力學(xué)、經(jīng)濟(jì)金融、環(huán)境科學(xué)、控制理論和高分子材料等諸多科學(xué)領(lǐng)域的研究中扮演重要角色。作為一類新興的數(shù)學(xué)模型,此類方程能夠?yàn)閹в洃浌δ�、自相似性質(zhì)和遺傳特征的復(fù)雜動(dòng)力學(xué)行為提供更為深刻和準(zhǔn)確的物理闡述,故具備傳統(tǒng)偏微分方程不可比擬的優(yōu)越性。由于解析技術(shù)的局限性,開展有關(guān)的數(shù)值方法研究具有重要的理論意義和實(shí)際應(yīng)用價(jià)值。本文研究時(shí)間和空間分?jǐn)?shù)階對(duì)流擴(kuò)散方程的數(shù)值解法,涉及有限差分和配置法、運(yùn)算矩陣法和微分求積法,主要內(nèi)容包括以下四個(gè)部分:首先,討論常系數(shù)時(shí)間分?jǐn)?shù)階擴(kuò)散方程的一種高精度數(shù)值方法。采用一類高階差分格式離散Caputo分?jǐn)?shù)階導(dǎo)數(shù),在空間上應(yīng)用指數(shù)樣條插值并給予擬一致配點(diǎn),構(gòu)造一種有限差分-指數(shù)B-樣條配置法。分析一階數(shù)值格式的唯一可解性和無條件穩(wěn)定性,其中,穩(wěn)定性通過分?jǐn)?shù)階von Neumann分析法證明。數(shù)值結(jié)果驗(yàn)證了理論分析結(jié)果。其次,研究變系數(shù)時(shí)間分?jǐn)?shù)階對(duì)流擴(kuò)散方程的運(yùn)算矩陣法。采用Chebyshev cardinal函數(shù)作為基函數(shù),推導(dǎo)Riemann-Liouville分?jǐn)?shù)階積分的一個(gè)運(yùn)算矩陣,給出兩種計(jì)算矩陣元素的方案,并分析它們的優(yōu)缺點(diǎn)。進(jìn)一步地,求得一、二階空間導(dǎo)數(shù)的矩陣近似,利用分?jǐn)?shù)階微積分的性質(zhì)將控制方程轉(zhuǎn)化成含分?jǐn)?shù)階積分的等價(jià)形式,以此提出一種高效Chebyshev cardinal運(yùn)算矩陣法。數(shù)值結(jié)果及與現(xiàn)有算法的比較表明該方法具備較高的計(jì)算精度,并且適用于長時(shí)間歷程分?jǐn)?shù)階問題的數(shù)值模擬。再次,建立一、二維變系數(shù)空間分?jǐn)?shù)階對(duì)流擴(kuò)散方程的樣條微分求積法。試函數(shù)選為三次B-樣條,介紹Riemann-Liouville分?jǐn)?shù)階導(dǎo)數(shù)的微分求積公式,給出確定加權(quán)系數(shù)的方法。同時(shí),為了快速計(jì)算分?jǐn)?shù)階導(dǎo)數(shù)在各個(gè)離散節(jié)點(diǎn)上的函數(shù)值,對(duì)分?jǐn)?shù)階導(dǎo)數(shù)的積分部分應(yīng)用分部積分公式并反復(fù)進(jìn)行遞歸,推出三次B-樣條的Riemann-Liouville分?jǐn)?shù)階導(dǎo)數(shù)的顯式表達(dá)式。導(dǎo)出的常微分方程組采用加權(quán)平均差分格式離散。該方法繼承了傳統(tǒng)微分求積法計(jì)算量小、精度高和易于編程的特點(diǎn),在相同的離散參數(shù)和誤差量級(jí)下,樣條微分求積法的CPU時(shí)耗遠(yuǎn)低于有限元法。最后,研究不規(guī)則區(qū)域上帶分?jǐn)?shù)階方向?qū)?shù)的二維變系數(shù)空間分?jǐn)?shù)階擴(kuò)散方程的徑向基微分求積法。對(duì)全域離散節(jié)點(diǎn)上的函數(shù)值進(jìn)行加權(quán)線性組合,引入Caputo分?jǐn)?shù)階方向?qū)?shù)的微分求積公式。試函數(shù)分別選為Multiquadric、Inverse Multiquadric和Gaussian徑向基函數(shù),給出三種確定加權(quán)系數(shù)的方法。以之為基礎(chǔ),發(fā)展一種適用于任意區(qū)域上的空間分?jǐn)?shù)階擴(kuò)散方程的Crank-Nicolson徑向基微分求積法,并給出算法流程。數(shù)值算例包含有矩形、梯形、圓形和L-型區(qū)域問題,數(shù)值結(jié)果說明了方法的靈活性和對(duì)不規(guī)則區(qū)域問題的適應(yīng)性。由于對(duì)維數(shù)變化不敏感,所提方法可以進(jìn)一步推廣至三維任意區(qū)域空間分?jǐn)?shù)階問題,并且不會(huì)引起計(jì)算量的大幅度增加。
【學(xué)位單位】：西北工業(yè)大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位年份】：2018
【中圖分類】：O241.8
【部分圖文】：

3-1 算例 1 中 α = 0.9, N = 5000和 M = 20時(shí)算法在 t = 1時(shí)刻的誤差隨張力參數(shù)的表 3-1 算例 1 中 α = 0.9, p = 5.16和 N = 5000時(shí)算法在 t = 1時(shí)刻的數(shù)值結(jié)果2e (τ , M)Cov. order e ( , M)∞τCov. or0 8.6951e-04 - 1.3969e-03 -0 2.2384e-04 1.9578 3.6042e-04 1.9540 5.7089e-05 1.9712 9.2262e-05 1.9650 1.5162e-05 1.9127 2.4710e-05 1.900 由于0( ) ( )CtD E t = E tα α αα αλ λ λ ，給定2κ =4 /π ，初邊值條件 ( x ) = sin(π x/ 2),1g (t ) = 0,2g (t ) = E ( t)αα,次的右端項(xiàng)，易驗(yàn)證在區(qū)間 Λ =[0,1]上方程的精確解是u ( x , t ) = E ( t ) sin( x/ 2)ααπ , E ( )α為 Mittag-Leffler 函數(shù)。選擇 α =0.3和張力參數(shù) p = 1.52，表 3-200時(shí)樣條配置法在 t = 1時(shí)刻時(shí)間方向的數(shù)值結(jié)果；表 3-3 則列出了 N =500

不同時(shí)刻下p=0.01,N=500*t和M=100時(shí)精確解和數(shù)值解的比較

(c) p = 2.53(d) p =3.35圖 3-3 不同的張力參數(shù)下 α = 0.6, N = 2500和 M = 50時(shí)空間網(wǎng)格上的絕對(duì)誤差分布表 3-5 算例 4 中 α = 0.6和 p = 2.53時(shí)兩種算法在 t = 1時(shí)刻的絕對(duì)誤差比較

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 邵克勇;郭浩軒;韓峰;王婷婷;;異結(jié)構(gòu)分?jǐn)?shù)階混沌系統(tǒng)的柔性變結(jié)構(gòu)同步控制[J];揚(yáng)州大學(xué)學(xué)報(bào)(自然科學(xué)版);2019年04期

2 擺玉龍;楊陽;魏強(qiáng);段濟(jì)開;范滿紅;;分?jǐn)?shù)階混沌系統(tǒng)的同步研究及電路實(shí)現(xiàn)[J];西北師范大學(xué)學(xué)報(bào)(自然科學(xué)版);2019年06期

3 盧寧;司輝;鄭永愛;;基于狀態(tài)觀測器的分?jǐn)?shù)階混沌系統(tǒng)的同步[J];電子設(shè)計(jì)工程;2019年22期

4 王嫻;李東;;分?jǐn)?shù)階混沌系統(tǒng)的間歇控制同步[J];重慶工商大學(xué)學(xué)報(bào)(自然科學(xué)版);2018年04期

5 王詩靖;王廷江;;分?jǐn)?shù)階RL_α-C_β并聯(lián)諧振電路的等效電路[J];物理通報(bào);2017年01期

6 王廷江;;分?jǐn)?shù)階RL_α-C_β并聯(lián)諧振頻率簡易表達(dá)式[J];物理通報(bào);2017年04期

7 王婷;張麗娟;達(dá)佳麗;;淺析分?jǐn)?shù)階微分方程三點(diǎn)共振邊值問題正解的存在[J];課程教育研究;2017年32期

8 李明;陳旭;鄭永愛;;一類分?jǐn)?shù)階混沌系統(tǒng)的自適應(yīng)滑模同步[J];揚(yáng)州大學(xué)學(xué)報(bào)(自然科學(xué)版);2016年03期

9 李特;袁建寶;吳瑩;;一類不確定分?jǐn)?shù)階混沌系統(tǒng)同步的自適應(yīng)滑�？刂品椒╗J];動(dòng)力學(xué)與控制學(xué)報(bào);2017年02期

10 賴滿豐;金玲玉;;分?jǐn)?shù)階Klein-Gordon-Schr?dinger方程弱解的存在性[J];佛山科學(xué)技術(shù)學(xué)院學(xué)報(bào)(自然科學(xué)版);2017年03期

相關(guān)會(huì)議論文 前10條

1 羅紹凱;李琳;許艷麗;賀金滿;張笑天;彭旺;戴蕓;;分?jǐn)?shù)階動(dòng)力學(xué)基本理論與方法的研究進(jìn)展[A];第十六屆全國交叉科學(xué)會(huì)議論文摘要集[C];2016年

2 羅紹凱;;分?jǐn)?shù)階動(dòng)力學(xué)基本理論與方法的研究進(jìn)展[A];第十二屆全國分析力學(xué)學(xué)術(shù)會(huì)議摘要集[C];2016年

3 孟瑞繁;殷德順;;分?jǐn)?shù)階本構(gòu)模型描述流變現(xiàn)象時(shí)間效應(yīng)的物理意義[A];中國力學(xué)大會(huì)-2015論文摘要集[C];2015年

4 李西成;;經(jīng)皮吸收的分?jǐn)?shù)階藥物動(dòng)力學(xué)模型[A];中國力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年

5 顧葆華;單梁;李軍;王執(zhí)銓;;一種新分?jǐn)?shù)階混沌系統(tǒng)及其復(fù)合快速同步控制[A];2009年中國智能自動(dòng)化會(huì)議論文集（第七分冊(cè)）[南京理工大學(xué)學(xué)報(bào)（增刊）][C];2009年

6 謝勇;;分?jǐn)?shù)階模型神經(jīng)元的動(dòng)力學(xué)行為及其同步[A];第四屆全國動(dòng)力學(xué)與控制青年學(xué)者研討會(huì)論文摘要集[C];2010年

7 張碩;于永光;王亞;;帶有時(shí)滯和隨機(jī)擾動(dòng)的不確定分?jǐn)?shù)階混沌系統(tǒng)準(zhǔn)同步[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

8 李常品;;分?jǐn)?shù)階動(dòng)力學(xué)的若干關(guān)鍵問題及研究進(jìn)展[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

9 王淑英;趙建峰;常迎香;李險(xiǎn)峰;;基于假分?jǐn)?shù)階統(tǒng)一混沌系統(tǒng)同步的圖像加密[A];第十五屆全國非線性振動(dòng)暨第十二屆全國非線性動(dòng)力學(xué)和運(yùn)動(dòng)穩(wěn)定性學(xué)術(shù)會(huì)議摘要集[C];2015年

10 劉杰;董鵬真;尚鋼;;分?jǐn)?shù)階非線性系統(tǒng)動(dòng)力學(xué)分析中數(shù)值算法可靠性及其誘導(dǎo)的復(fù)雜現(xiàn)象[A];中國力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年

相關(guān)博士學(xué)位論文 前10條

1 高獻(xiàn)成;兩維隨機(jī)旋轉(zhuǎn)歐拉流和一些隨機(jī)偏微分方程的穩(wěn)定性的研究[D];南京師范大學(xué);2019年

2 沈天龍;隨機(jī)分?jǐn)?shù)階偏微分方程的動(dòng)力學(xué)[D];國防科技大學(xué);2017年

3 朱曉鋼;分?jǐn)?shù)階對(duì)流擴(kuò)散方程的幾種數(shù)值方法研究[D];西北工業(yè)大學(xué);2018年

4 林世敏;若干分?jǐn)?shù)階微分方程的理論和數(shù)值方法研究[D];廈門大學(xué);2017年

5 呂春婉;分?jǐn)?shù)階微分方程的若干高階數(shù)值方法研究[D];廈門大學(xué);2017年

6 王楠;幾類空間分?jǐn)?shù)階偏微分方程的高效算法研究[D];華中科技大學(xué);2019年

7 費(fèi)明發(fā);幾類分?jǐn)?shù)階偏微分方程的譜方法研究[D];華中科技大學(xué);2019年

8 程廷治;非線性高階微分方程和分?jǐn)?shù)階方程解的定性分析[D];上海交通大學(xué);2018年

9 龍騰;Caputo型分?jǐn)?shù)階微分方程高精度數(shù)值方法研究[D];哈爾濱工業(yè)大學(xué);2018年

10 王康樂;幾類分?jǐn)?shù)階微分方程的近似解析解[D];西安電子科技大學(xué);2017年

相關(guān)碩士學(xué)位論文 前10條

1 翟如夢(mèng);自適應(yīng)分?jǐn)?shù)階的偏微分方程圖像去噪方法研究[D];重慶郵電大學(xué);2019年

2 朱姍姍;關(guān)于分?jǐn)?shù)階趨化-流體模型解的全局存在性及漸近行為的研究[D];揚(yáng)州大學(xué);2019年

3 龍騰;具Hilfer導(dǎo)數(shù)的分?jǐn)?shù)階微分方程邊值問題的研究[D];湘潭大學(xué);2019年

4 胡冬冬;分?jǐn)?shù)階非線性反應(yīng)擴(kuò)散方程的高精度算法[D];湘潭大學(xué);2019年

5 黃小輝;分?jǐn)?shù)階延遲積分及偏微分代數(shù)方程的波形松弛方法[D];湘潭大學(xué);2019年

6 蔣慧靈;時(shí)—空分?jǐn)?shù)階擴(kuò)散方程的擬緊差分算法[D];湘潭大學(xué);2019年

7 郭韜;一類帶分?jǐn)?shù)階量子修正的Zakharov系統(tǒng)的Fourier譜方法[D];湘潭大學(xué);2019年

8 呂龍波;非線性分?jǐn)?shù)階拉普拉斯方程的波節(jié)解[D];江西師范大學(xué);2019年

9 劉妍月;多維空間分?jǐn)?shù)階Schr(?)dinger方程的Fourier擬譜方法[D];湘潭大學(xué);2019年

10 謝柳柳;分?jǐn)?shù)階p-Laplace方程解的性質(zhì)的研究[D];南京航空航天大學(xué);2019年

本文編號(hào)：2807679