求解一類多項四階時間分數(shù)階擴散波方程的有限差分方法
發(fā)布時間:2021-02-17 06:11
分數(shù)階微分方程是由整數(shù)階微分方程推廣得到的.由于全局相關(guān)性,它能更好地刻畫各種模型的物理過程,因此,分數(shù)階微分方程的理論和數(shù)值方法是目前的熱點研究課題之一.本文主要研究求解一類多項四階時間分數(shù)階擴散波方程的有限差分方法.首先,研究求解一類多項四階時間分數(shù)階慢擴散方程的有限差分方法.應用L1公式逼近時間分數(shù)階導數(shù),用降階法處理空間四階導數(shù)項,再借助離散能量方法證明所得差分格式是無條件穩(wěn)定的且在無窮范數(shù)下其收斂階為O(τ2-α2+h2),其中τ和h分別為時間方向和空間方向步長,α2為時間分數(shù)階導數(shù)的最大階數(shù).最后用數(shù)值實驗驗證所提出差分格式的精度和有效性.其次,討論求解此類多項四階時間分數(shù)階慢擴散方程的高階數(shù)值算法.先用降階法,將原方程等價轉(zhuǎn)化為一個低階方程組,再對相應的離散方程兩邊作用一個平均算子,應用L1公式逼近時間分數(shù)階導數(shù),空間導數(shù)采用緊逼近,建立高階差分格式.借助離散能量方法,分析所得差分格式在L∞范數(shù)下是無條件穩(wěn)定且收斂階為O(t2-a2...
【文章來源】:南京郵電大學江蘇省
【文章頁數(shù)】:57 頁
【學位級別】:碩士
【部分圖文】:
方程(2.78)-(2.81)的精確解(左)、格式(2.23)-(2.25)(中)和(2.54)-(2.56)(右)求得的數(shù)值解
圖 2.1 方程(2.74)-(2.77)的精確解(左)、格式(2.23)-(2.25)(中)和(2.54)-(2.56) (右)求得的數(shù)值解圖 2.2 方程(2.78)-(2.81)的精確解(左)、格式(2.23)-(2.25)(中)和(2.54)-(2.56) (右)求得的數(shù)值解表 2.1 (實驗一)當1 2 0.3, 0.5時的最大誤差和空間方向收斂階h格式(2.23)-(2.25) τ 1/10000 格式(2.54)-(2.56) τ 1/200000 2 ,E h
該問題的精確解為 3u x , t t 1 sin x.圖 3.2 分別給出四階時間分數(shù)階波方程(3.70)-(3.73)的精確解,應用差分格式(3.24)-(3.26)以及差分格式(3.54)-(3.56)計算所得的數(shù)值解 1 2 1.2, 1.7, M 100, N 100.由圖可見,兩個差分格式數(shù)值解均與精確解非常吻合,進而表明兩個差分格式都是有效性的.定義 0, maxn nn NE h U u .在固定時間步長 充分小時,表 3.1 和表 3.2 分別給出實驗三和實驗四在不同空間步長下差分格式(3.24)-(3.26)和差分格式(3.54)-(3.56)的計算結(jié)果.可以看出,差分格式(3.24)-(3.26)空間方向的收斂階為二階,而差分格式(3.54)-(3.56)空間方向可達到四階收斂.表3.3和表3.4分別測試兩個差分格式時間方向的數(shù)值收斂階.當1 2 1.2, 1.7時,時間方向的收斂階均約為23 1.3;當1 2 1.4, 1.9時,時間方向的收斂階均約為23 1.1,與理論結(jié)果吻合.
【參考文獻】:
期刊論文
[1]一維連續(xù)函數(shù)的Riemann-Liouville分數(shù)階微積分[J]. 梁永順,蘇維宜. 中國科學:數(shù)學. 2016(04)
[2]高維分數(shù)階cable方程隱式差分逼近[J]. 馬亮亮,劉冬兵. 遼寧工程技術(shù)大學學報(自然科學版). 2014(04)
[3]數(shù)字圖像處理中的偏微分方程方法綜述[J]. 丁暢,尹清波,魯明羽. 計算機科學. 2013(S2)
[4]兩類分數(shù)階對流-擴散方程的有限差分方法[J]. 張紅玉,崔明榮. 山東大學學報(理學版). 2012(06)
[5]二維分數(shù)階變系數(shù)對流-擴散方程的數(shù)值解[J]. 劉青霞,劉發(fā)旺. 高等學校計算數(shù)學學報. 2011(01)
[6]用物理學熱傳導理論剖析二階齊次線性偏微分方程的本質(zhì)[J]. 董麗華. 長春師范學院學報. 2007(04)
博士論文
[1]分數(shù)階滑?刂评碚摷捌鋺醚芯縖D]. 鄧立為.哈爾濱工業(yè)大學 2014
[2]幾類分數(shù)階微分方程的數(shù)值方法研究[D]. 肖靜宇.哈爾濱工業(yè)大學 2013
[3]基于偏微分方程的圖像分割和腫瘤生物醫(yī)學建模研究[D]. 張魯筠.山東大學 2012
碩士論文
[1]兩類Caputo型分數(shù)階差分方程的初值問題[D]. 鄭金娣.安徽大學 2017
[2]基于分數(shù)階非牛頓流體本構(gòu)關(guān)系的雙孔介質(zhì)波動方程[D]. 熊繁升.中國地質(zhì)大學(北京) 2016
[3]時間分數(shù)階偏微分方程高階數(shù)值解法[D]. 紀翠翠.東南大學 2015
[4]熱量傳遞的分數(shù)階微分方程模型與數(shù)值模擬[D]. 李丹.華北理工大學 2015
[5]加權(quán)Bergman空間上Toeplitz算子的零積問題和Mellin變換[D]. 韓丹丹.河北師范大學 2014
本文編號:3037554
【文章來源】:南京郵電大學江蘇省
【文章頁數(shù)】:57 頁
【學位級別】:碩士
【部分圖文】:
方程(2.78)-(2.81)的精確解(左)、格式(2.23)-(2.25)(中)和(2.54)-(2.56)(右)求得的數(shù)值解
圖 2.1 方程(2.74)-(2.77)的精確解(左)、格式(2.23)-(2.25)(中)和(2.54)-(2.56) (右)求得的數(shù)值解圖 2.2 方程(2.78)-(2.81)的精確解(左)、格式(2.23)-(2.25)(中)和(2.54)-(2.56) (右)求得的數(shù)值解表 2.1 (實驗一)當1 2 0.3, 0.5時的最大誤差和空間方向收斂階h格式(2.23)-(2.25) τ 1/10000 格式(2.54)-(2.56) τ 1/200000 2 ,E h
該問題的精確解為 3u x , t t 1 sin x.圖 3.2 分別給出四階時間分數(shù)階波方程(3.70)-(3.73)的精確解,應用差分格式(3.24)-(3.26)以及差分格式(3.54)-(3.56)計算所得的數(shù)值解 1 2 1.2, 1.7, M 100, N 100.由圖可見,兩個差分格式數(shù)值解均與精確解非常吻合,進而表明兩個差分格式都是有效性的.定義 0, maxn nn NE h U u .在固定時間步長 充分小時,表 3.1 和表 3.2 分別給出實驗三和實驗四在不同空間步長下差分格式(3.24)-(3.26)和差分格式(3.54)-(3.56)的計算結(jié)果.可以看出,差分格式(3.24)-(3.26)空間方向的收斂階為二階,而差分格式(3.54)-(3.56)空間方向可達到四階收斂.表3.3和表3.4分別測試兩個差分格式時間方向的數(shù)值收斂階.當1 2 1.2, 1.7時,時間方向的收斂階均約為23 1.3;當1 2 1.4, 1.9時,時間方向的收斂階均約為23 1.1,與理論結(jié)果吻合.
【參考文獻】:
期刊論文
[1]一維連續(xù)函數(shù)的Riemann-Liouville分數(shù)階微積分[J]. 梁永順,蘇維宜. 中國科學:數(shù)學. 2016(04)
[2]高維分數(shù)階cable方程隱式差分逼近[J]. 馬亮亮,劉冬兵. 遼寧工程技術(shù)大學學報(自然科學版). 2014(04)
[3]數(shù)字圖像處理中的偏微分方程方法綜述[J]. 丁暢,尹清波,魯明羽. 計算機科學. 2013(S2)
[4]兩類分數(shù)階對流-擴散方程的有限差分方法[J]. 張紅玉,崔明榮. 山東大學學報(理學版). 2012(06)
[5]二維分數(shù)階變系數(shù)對流-擴散方程的數(shù)值解[J]. 劉青霞,劉發(fā)旺. 高等學校計算數(shù)學學報. 2011(01)
[6]用物理學熱傳導理論剖析二階齊次線性偏微分方程的本質(zhì)[J]. 董麗華. 長春師范學院學報. 2007(04)
博士論文
[1]分數(shù)階滑?刂评碚摷捌鋺醚芯縖D]. 鄧立為.哈爾濱工業(yè)大學 2014
[2]幾類分數(shù)階微分方程的數(shù)值方法研究[D]. 肖靜宇.哈爾濱工業(yè)大學 2013
[3]基于偏微分方程的圖像分割和腫瘤生物醫(yī)學建模研究[D]. 張魯筠.山東大學 2012
碩士論文
[1]兩類Caputo型分數(shù)階差分方程的初值問題[D]. 鄭金娣.安徽大學 2017
[2]基于分數(shù)階非牛頓流體本構(gòu)關(guān)系的雙孔介質(zhì)波動方程[D]. 熊繁升.中國地質(zhì)大學(北京) 2016
[3]時間分數(shù)階偏微分方程高階數(shù)值解法[D]. 紀翠翠.東南大學 2015
[4]熱量傳遞的分數(shù)階微分方程模型與數(shù)值模擬[D]. 李丹.華北理工大學 2015
[5]加權(quán)Bergman空間上Toeplitz算子的零積問題和Mellin變換[D]. 韓丹丹.河北師范大學 2014
本文編號:3037554
本文鏈接:http://www.sikaile.net/kejilunwen/yysx/3037554.html
最近更新
教材專著