發(fā)布時(shí)間：2020-04-14 10:24

【摘要】：考慮如下時(shí)空分?jǐn)?shù)階變系數(shù)擴(kuò)散方程其中,Ω=(0,1),D是一階導(dǎo)數(shù)算子,0Dx1-β和xD11-β分別表示1-β(0β1)階左、右Riemman-Liouville分?jǐn)?shù)階導(dǎo)數(shù)算子,0CDtαu表示α(0α1)階Caputo分?jǐn)?shù)階導(dǎo)數(shù)算子,當(dāng)α = 1時(shí),0CDtαu=ut為時(shí)間的一階導(dǎo)數(shù);K(x)為具有正的上下界的擴(kuò)散系數(shù),即存在Kmin、Kmaz0,0Kmin ≤ K(x)≤ Kmax+∞,x ∈(0,1);f ∈ L2(Ω)為源項(xiàng)或匯項(xiàng),0≤θ ≤ 1表示粒子向前相對(duì)于向后轉(zhuǎn)移概率的權(quán)重.對(duì)上述變系數(shù)時(shí)空分?jǐn)?shù)階擴(kuò)散方程直接采用Galerkin方法所得到的雙線形式不能保證其強(qiáng)制性[15],因此,不能直接利用有限元方法求解.為克服此困難,我們通過(guò)引入通量p=-K(x)Du作為中間變量,將方程(0.0.1)化為變系數(shù)整數(shù)階問(wèn)題和系數(shù)為1的雙邊分?jǐn)?shù)階問(wèn)題,建立混合型變分形式,據(jù)此構(gòu)造了混合型有限元方法.進(jìn)一步,我們對(duì)α階時(shí)間導(dǎo)數(shù)進(jìn)行離散獲得全離散混合型有限元格式,即當(dāng)α = 1時(shí),我們利用向后Euler方法對(duì)時(shí)間導(dǎo)數(shù)進(jìn)行離散,對(duì)空間分?jǐn)?shù)階導(dǎo)數(shù)分別采用混合型有限元和混合型間斷有限元方法離散;當(dāng)0α1時(shí),對(duì)Caputo分?jǐn)?shù)階時(shí)間導(dǎo)數(shù)應(yīng)用L1方法離散,對(duì)空間分?jǐn)?shù)階導(dǎo)數(shù)采用混合型有限元方法離散,從而,構(gòu)造出了數(shù)值模擬方程(0.0.1)的全離散混合型有限元格式.進(jìn)而巧妙地利用反對(duì)稱矩陣分和塊矩陣的代數(shù)性質(zhì),證明了由混合型有限元格式導(dǎo)出的有限元方程組的系數(shù)矩陣可逆和有限元解存在唯一.在有限元方程組的求解過(guò)程中,我們發(fā)現(xiàn)由分?jǐn)?shù)階算子的非局部性而產(chǎn)生的系數(shù)矩陣是由四個(gè)分塊矩陣構(gòu)成的非對(duì)稱非稀疏矩陣.傳統(tǒng)的Gauss消去方法求解需要O(N3)的計(jì)算量和O(N2)的存儲(chǔ)量(N為剖分網(wǎng)格的節(jié)點(diǎn)數(shù)),為提高計(jì)算效率,引入快速傅里葉變換(FFT),局部處理Toeplitz-like結(jié)構(gòu)的分塊矩陣,降低局部計(jì)算量為O(NlogN)和存儲(chǔ)量O(N),進(jìn)而保證系數(shù)矩陣的計(jì)算量為O(N2)、存儲(chǔ)量為O(N).注意到求解有限元方程組解的計(jì)算效率不僅與計(jì)算量和存儲(chǔ)量有關(guān),還與迭代次數(shù)有關(guān).而系數(shù)矩陣為具有壞條件數(shù)的非稀疏矩陣,其迭代次數(shù)會(huì)隨著剖分的加密顯著增長(zhǎng).為降低迭代次數(shù),我們改進(jìn)傳統(tǒng)的基于Toeplitz分塊矩陣和循環(huán)矩陣表達(dá)的預(yù)處理子,結(jié)合穩(wěn)定的雙共軛梯算法(BiCG-STAB),構(gòu)造了預(yù)條件快速穩(wěn)定的雙共軛梯度算法(PFBiCG).新的PFBiCG算法在保證每次迭代計(jì)算量為O(N2)、存儲(chǔ)量為O(N)的前提下,較傳統(tǒng)的Guass消去法和穩(wěn)定的雙共軛梯度算法(BiCG-STAB),明顯降低了迭代次數(shù),縮短了計(jì)算時(shí)間.本文中的數(shù)值實(shí)驗(yàn)提高了計(jì)算效率,也驗(yàn)證了實(shí)驗(yàn)結(jié)果與理論預(yù)測(cè)的一致性.
【學(xué)位授予單位】：山東師范大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2018
【分類號(hào)】：O241.82

【相似文獻(xiàn)】

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

1 高哲;廖曉鐘;;一種線性分?jǐn)?shù)階系統(tǒng)穩(wěn)定性的頻域判別準(zhǔn)則[J];自動(dòng)化學(xué)報(bào);2011年11期

2 李周紅;蔣涓;;基于積極控制的兩個(gè)不同分?jǐn)?shù)階混沌系統(tǒng)的反同步[J];玉溪師范學(xué)院學(xué)報(bào);2018年04期

3 楊宸;;基于對(duì)角占優(yōu)準(zhǔn)則的分?jǐn)?shù)階系統(tǒng)同步控制[J];科技展望;2014年13期

4 陸瑤;;基于調(diào)制函數(shù)法的分?jǐn)?shù)階系統(tǒng)參數(shù)辨識(shí)[J];科學(xué)技術(shù)創(chuàng)新;2018年33期

5 楊朋飛;劉冬梅;高寧;白珍龍;;分?jǐn)?shù)階控制系統(tǒng)的穩(wěn)定性理論研究[J];儀器儀表用戶;2018年05期

6 毛北行;;兩類分?jǐn)?shù)階系統(tǒng)的觀測(cè)器同步[J];吉林大學(xué)學(xué)報(bào)(理學(xué)版);2017年01期

7 劉煥霞;趙鑫;林崇;馬瑞蘭;;廣義分?jǐn)?shù)階混沌系統(tǒng)的魯棒同步研究[J];青島大學(xué)學(xué)報(bào)(工程技術(shù)版);2018年02期

8 賈雅琦;;時(shí)滯分?jǐn)?shù)階混沌系統(tǒng)的完全同步[J];計(jì)算機(jī)產(chǎn)品與流通;2018年07期

9 劉瑞娟;康美玲;徐欣;林海燕;;一類分?jǐn)?shù)階不確定重復(fù)控制系統(tǒng)的穩(wěn)定性分析[J];廈門理工學(xué)院學(xué)報(bào);2018年05期

10 李慶賓;李亮;;一類不確定分?jǐn)?shù)階混沌系統(tǒng)的參數(shù)辨識(shí)[J];數(shù)學(xué)的實(shí)踐與認(rèn)識(shí);2018年08期

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

1 蔣靜菲;曹登慶;陳華濤;;一類可變分?jǐn)?shù)階系統(tǒng)的穩(wěn)定性與控制[A];第十五屆全國(guó)非線性振動(dòng)暨第十二屆全國(guó)非線性動(dòng)力學(xué)和運(yùn)動(dòng)穩(wěn)定性學(xué)術(shù)會(huì)議摘要集[C];2015年

2 申永軍;楊紹普;溫少芳;牛江川;;分?jǐn)?shù)階系統(tǒng)的動(dòng)力學(xué)與控制[A];第十屆動(dòng)力學(xué)與控制學(xué)術(shù)會(huì)議摘要集[C];2016年

3 薛定宇;白鷺;;分?jǐn)?shù)階系統(tǒng)的仿真方法(英文)[A];系統(tǒng)仿真技術(shù)及其應(yīng)用學(xué)術(shù)論文集（第15卷）[C];2014年

4 王在華;;分?jǐn)?shù)階系統(tǒng)的實(shí)驗(yàn)建模、穩(wěn)定性分析與數(shù)值求解[A];第六屆全國(guó)動(dòng)力學(xué)與控制青年學(xué)者學(xué)術(shù)研討會(huì)論文摘要集[C];2012年

5 申永軍;楊紹普;邢海軍;;單自由度分?jǐn)?shù)階系統(tǒng)的動(dòng)力學(xué)分析[A];第九屆全國(guó)動(dòng)力學(xué)與控制學(xué)術(shù)會(huì)議會(huì)議手冊(cè)[C];2012年

6 胡軼超;李東海;;分?jǐn)?shù)階系統(tǒng)的PID控制器整定[A];中國(guó)自動(dòng)化學(xué)會(huì)控制理論專業(yè)委員會(huì)D卷[C];2011年

7 許建強(qiáng);;參數(shù)不確定分?jǐn)?shù)階統(tǒng)一混沌系統(tǒng)的自適應(yīng)同步[A];中國(guó)自動(dòng)化學(xué)會(huì)控制理論專業(yè)委員會(huì)C卷[C];2011年

8 劉瑞娟;聶卓峗;潘秋萍;;基于等價(jià)輸入干擾方法的分?jǐn)?shù)階系統(tǒng)擾動(dòng)抑制[A];第36屆中國(guó)控制會(huì)議論文集（G）[C];2017年

9 馬周健;王軍;;分?jǐn)?shù)階超混沌系統(tǒng)基于廣義同步的混沌遮掩保密通信[A];第19屆中國(guó)系統(tǒng)仿真技術(shù)及其應(yīng)用學(xué)術(shù)年會(huì)論文集（19th CCSSTA 2018）[C];2018年

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

1 王飛;分?jǐn)?shù)階網(wǎng)絡(luò)動(dòng)態(tài)分析與控制方法研究[D];江南大學(xué);2017年

2 杜s，

本文編號(hào)：2627197