天堂国产午夜亚洲专区-少妇人妻综合久久蜜臀-国产成人户外露出视频在线-国产91传媒一区二区三区

當(dāng)前位置:主頁 > 科技論文 > 數(shù)學(xué)論文 >

時(shí)—空分?jǐn)?shù)階擴(kuò)散方程的快速算法以及MT-TSCR-FDE的快速數(shù)值解法

發(fā)布時(shí)間:2018-07-26 20:47
【摘要】:如今,分?jǐn)?shù)階微積分已成為流行在社會(huì)科學(xué)與工程的重要工具。特別是時(shí)空分?jǐn)?shù)階擴(kuò)散方程正越來越多地應(yīng)用于研究許多領(lǐng)域的反常擴(kuò)散現(xiàn)象。由于分?jǐn)?shù)階導(dǎo)數(shù)的非局部性,其數(shù)值求解會(huì)生成滿系數(shù)矩陣,并且此矩陣的求解需要消耗大量的計(jì)算成本及存儲(chǔ)量。所以我們展開研究快速數(shù)值方法來解決這一問題。第一步:我們給出時(shí)-空分?jǐn)?shù)階雙邊擴(kuò)散方程的一般形式:對(duì)于時(shí)間分?jǐn)?shù)階,我們采用Caputo分?jǐn)?shù)階導(dǎo)數(shù);對(duì)于空間分?jǐn)?shù)階,采用左和右的Ricmann-Liouville空間分?jǐn)?shù)階導(dǎo)數(shù),并利用修正的Grunwald-Letnikov近似。給出相應(yīng)的有限差分格式及矩陣格式。第二步:分析有限差分格式的滿系數(shù)矩陣,它可以分解成Toeplitz矩陣與向量乘積之和。根據(jù)Toeplitz矩陣與循環(huán)矩陣的關(guān)系,以及循環(huán)矩陣的性質(zhì),提出用傅里葉變換法求解矩陣向量的乘積,開發(fā)一個(gè)基于快速傅里葉變換(FFT)的快速解決方案。第三步:基于FFT的快速解決方案,我們開發(fā)了兩種快速數(shù)值方法:一種是O(N log N)的最小剩余共軛梯度平方快速迭代方法,一種是O(N log2 N)的快速有限差分法。與常規(guī)有限差分法相比能夠極大地減少計(jì)算成本和存儲(chǔ)空間,同時(shí)保持相同的精度。第四步:對(duì)于多項(xiàng)分?jǐn)?shù)階時(shí)-空Caputo-Resz方程(MT-TSCR-FDE), P(Dt)u(x,t)=p(x)Rxβ+q(x)Rxγ-h(x)u(x,t)+F(x,t)首先利用預(yù)估-校正法對(duì)此多項(xiàng)分?jǐn)?shù)階方程進(jìn)行數(shù)值逼近,然后應(yīng)用上面提出的方法和步驟,分析研究數(shù)值求解的快速方法,以極大地減少計(jì)算成本和存儲(chǔ)空間。第五步:數(shù)值實(shí)驗(yàn)。分別針對(duì)時(shí)-空Caputo-Riesz分?jǐn)?shù)階擴(kuò)散方程,一個(gè)有解析解的時(shí)-空分?jǐn)?shù)階擴(kuò)散方程,以及最后的MT-TSCR-FDE給出實(shí)例數(shù)值模擬,給出相應(yīng)的誤差分析,以及CPU時(shí)間分析,通過優(yōu)良的數(shù)值結(jié)果證明本文提出的幾個(gè)快速方法的有效性。
[Abstract]:Today, fractional calculus has become an important tool in social science and engineering. In particular, fractional diffusion equations are increasingly used to study anomalous diffusion phenomena in many fields. Because of the nonlocality of fractional derivative, the numerical solution of the matrix will generate a full coefficient matrix, and the solution of this matrix needs a lot of computation cost and storage. So we study the fast numerical method to solve this problem. The first step: we give the general form of the time-space fractional two-sided diffusion equation: for the time fractional order, we use the Caputo fractional derivative; for the space fractional order, we use the left and right Ricmann-Liouville space fractional derivative. The modified Grunwald-Letnikov approximation is used. The corresponding finite difference scheme and matrix scheme are given. The second step is to analyze the full coefficient matrix of finite difference scheme, which can be decomposed into the sum of the product of Toeplitz matrix and vector. According to the relationship between Toeplitz matrix and cyclic matrix, and the properties of cyclic matrix, a fast solution based on fast Fourier transform (FFT) is developed by using Fourier transform method to solve the product of matrix vector. Step 3: based on the fast solution of FFT, we develop two fast numerical methods: one is the fast iterative method of minimum residual conjugate gradient square of O (N log N) and the other is the fast finite difference method of O (N log2 N). Compared with the conventional finite difference method, the computational cost and storage space can be greatly reduced, while maintaining the same accuracy. The fourth step: for the multiterm fractional space-time Caputo-Resz equation (MT-TSCR-FDE), P (Dt) u (XT) p (x) Rx 尾 q (x) Rx 緯 -h (x) u (XT) F (XT, we first use the predictor-correction method to approximate the multivariate fractional order equation, and then apply the methods and steps mentioned above to analyze and study the fast method of numerical solution. To greatly reduce computational costs and storage space. Step 5: numerical experiment. For time-space Caputo-Riesz fractional order diffusion equation, a time-space fractional diffusion equation with analytic solution, and the final MT-TSCR-FDE numerical simulation, the corresponding error analysis and CPU time analysis are given. The validity of several fast methods proposed in this paper is proved by good numerical results.
【學(xué)位授予單位】:山東大學(xué)
【學(xué)位級(jí)別】:碩士
【學(xué)位授予年份】:2015
【分類號(hào)】:O241.82

【相似文獻(xiàn)】

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

1 王德金;鄭永愛;;分?jǐn)?shù)階混沌系統(tǒng)的延遲同步[J];動(dòng)力學(xué)與控制學(xué)報(bào);2010年04期

2 楊晨航,劉發(fā)旺;分?jǐn)?shù)階Relaxation-Oscillation方程的一種分?jǐn)?shù)階預(yù)估-校正方法[J];廈門大學(xué)學(xué)報(bào)(自然科學(xué)版);2005年06期

3 王發(fā)強(qiáng);劉崇新;;分?jǐn)?shù)階臨界混沌系統(tǒng)及電路實(shí)驗(yàn)的研究[J];物理學(xué)報(bào);2006年08期

4 夏源;吳吉春;;分?jǐn)?shù)階對(duì)流——彌散方程的數(shù)值求解[J];南京大學(xué)學(xué)報(bào)(自然科學(xué)版);2007年04期

5 張隆閣;;一類參數(shù)不確定混沌系統(tǒng)的分?jǐn)?shù)階自適應(yīng)同步[J];中國科技信息;2009年15期

6 陳世平;劉發(fā)旺;;一維分?jǐn)?shù)階滲透方程的數(shù)值模擬[J];高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào);2010年04期

7 辛寶貴;陳通;劉艷芹;;一類分?jǐn)?shù)階混沌金融系統(tǒng)的復(fù)雜性演化研究[J];物理學(xué)報(bào);2011年04期

8 黃睿暉;;分?jǐn)?shù)階微方程的迭代方法研究[J];長春理工大學(xué)學(xué)報(bào);2011年06期

9 蔣曉蕓,徐明瑜;分形介質(zhì)分?jǐn)?shù)階反常守恒擴(kuò)散模型及其解析解[J];山東大學(xué)學(xué)報(bào)(理學(xué)版);2003年05期

10 陳玉霞;高金峰;;一個(gè)新的分?jǐn)?shù)階混沌系統(tǒng)[J];鄭州大學(xué)學(xué)報(bào)(理學(xué)版);2009年04期

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

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

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

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

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

5 李常品;;分?jǐn)?shù)階動(dòng)力學(xué)簡介[A];第三屆海峽兩岸動(dòng)力學(xué)、振動(dòng)與控制學(xué)術(shù)會(huì)議論文摘要集[C];2013年

6 蔣曉蕓;徐明瑜;;時(shí)間依靠分?jǐn)?shù)階Schr銉dinger方程中的可動(dòng)邊界問題[A];中國力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年

7 王花;;分?jǐn)?shù)階混沌系統(tǒng)的同步在圖像加密中的應(yīng)用[A];第二屆全國隨機(jī)動(dòng)力學(xué)學(xué)術(shù)會(huì)議摘要集與會(huì)議議程[C];2013年

8 王在華;;分?jǐn)?shù)階動(dòng)力系統(tǒng)的若干問題[A];第三屆全國動(dòng)力學(xué)與控制青年學(xué)者研討會(huì)論文摘要集[C];2009年

9 張碩;于永光;王莎;;帶有時(shí)滯和隨機(jī)擾動(dòng)的分?jǐn)?shù)階混沌系統(tǒng)同步[A];第十四屆全國非線性振動(dòng)暨第十一屆全國非線性動(dòng)力學(xué)和運(yùn)動(dòng)穩(wěn)定性學(xué)術(shù)會(huì)議摘要集與會(huì)議議程[C];2013年

10 李西成;;一個(gè)具有糊狀區(qū)的分?jǐn)?shù)階可動(dòng)邊界問題的相似解研究[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

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

1 陳善鎮(zhèn);兩類空間分?jǐn)?shù)階偏微分方程模型有限差分逼近的若干研究[D];山東大學(xué);2015年

2 任永強(qiáng);油藏與二氧化碳埋存問題的數(shù)值模擬與不確定性量化分析以及分?jǐn)?shù)階微分方程的數(shù)值方法[D];山東大學(xué);2015年

3 蔣敏;分?jǐn)?shù)階微分方程理論分析與應(yīng)用問題的研究[D];電子科技大學(xué);2015年

4 卜紅霞;基于分?jǐn)?shù)階傅里葉域稀疏表征的CS-SAR成像理論與算法研究[D];北京理工大學(xué);2015年

5 王少偉;分?jǐn)?shù)階微積分理論在粘彈性流體力學(xué)及量子力學(xué)中的某些應(yīng)用[D];山東大學(xué);2007年

6 郭霄怡;分?jǐn)?shù)階微積分在量子力學(xué)和非牛頓流體力學(xué)研究中的某些應(yīng)用[D];山東大學(xué);2007年

7 陳景華;空間和時(shí)間分?jǐn)?shù)階偏微分方程[D];廈門大學(xué);2007年

8 趙瀟;分?jǐn)?shù)階粘彈性本構(gòu)參數(shù)識(shí)別[D];大連理工大學(xué);2012年

9 趙振剛;三類分?jǐn)?shù)階偏微分方程的有限元計(jì)算[D];上海大學(xué);2011年

10 郭鵬;隨機(jī)中的分?jǐn)?shù)階微分方程的數(shù)值計(jì)算[D];上海大學(xué);2012年

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

1 黃志穎;非線性時(shí)間分?jǐn)?shù)階微分方程的數(shù)值解法[D];華南理工大學(xué);2015年

2 趙九龍;基于分?jǐn)?shù)階微積分的三維圖像去噪增強(qiáng)算法研究[D];寧夏大學(xué);2015年

3 楚彩虹;單載波分?jǐn)?shù)階傅里葉域均衡系統(tǒng)及關(guān)鍵技術(shù)研究[D];鄭州大學(xué);2015年

4 全曉靜;非線性分?jǐn)?shù)階積分方程的Adomian解法[D];寧夏大學(xué);2015年

5 黃潔;非線性分?jǐn)?shù)階Volterra積分微分方程的小波數(shù)值解法[D];寧夏大學(xué);2015年

6 莊嶠;復(fù)合介質(zhì)中時(shí)間分?jǐn)?shù)階熱傳導(dǎo)正逆問題及其應(yīng)用研究[D];山東大學(xué);2015年

7 高素娟;分?jǐn)?shù)階延遲偏微分方程的緊致有限差分方法[D];山東大學(xué);2015年

8 趙珊珊;時(shí)—空分?jǐn)?shù)階擴(kuò)散方程的快速算法以及MT-TSCR-FDE的快速數(shù)值解法[D];山東大學(xué);2015年

9 王珍;分?jǐn)?shù)階奇異邊值問題的研究[D];山東師范大學(xué);2015年

10 馮靜;一類分?jǐn)?shù)階奇異脈沖邊值問題正解的存在性研究[D];山東師范大學(xué);2015年

,

本文編號(hào):2147243

資料下載
論文發(fā)表

本文鏈接:http://sikaile.net/kejilunwen/yysx/2147243.html


Copyright(c)文論論文網(wǎng)All Rights Reserved | 網(wǎng)站地圖 |

版權(quán)申明:資料由用戶b0557***提供,本站僅收錄摘要或目錄,作者需要?jiǎng)h除請(qǐng)E-mail郵箱bigeng88@qq.com
国产亚洲不卡一区二区| 亚洲中文字幕在线观看四区| 午夜激情视频一区二区| 九九热视频免费在线视频| 天堂网中文字幕在线观看| 午夜精品福利视频观看| 麻豆国产精品一区二区| 国产精品人妻熟女毛片av久| 国产丝袜美女诱惑一区二区| 欧美久久一区二区精品| 老司机精品线观看86| 国产精品超碰在线观看| 好吊妞在线免费观看视频| 福利一区二区视频在线| 精品国产亚洲一区二区三区| 久久99亚洲小姐精品综合| 草草视频福利在线观看| 白白操白白在线免费观看| 欧美加勒比一区二区三区| 夫妻性生活黄色录像视频| 美日韩一区二区精品系列| 午夜精品在线视频一区| 一区二区三区亚洲天堂| 国产美女精品午夜福利视频| 美女被草的视频在线观看| 免费特黄一级一区二区三区| 国产精品色热综合在线| 日本东京热加勒比一区二区| 色播五月激情五月婷婷| 草草草草在线观看视频| 久久经典一区二区三区| 99久久国产精品亚洲| 日韩人妻中文字幕精品| 亚洲欧美日韩另类第一页| 人妻熟女欲求不满一区二区| 欧美丝袜诱惑一区二区| 精品视频一区二区三区不卡| 亚洲女同一区二区另类| 福利视频一区二区三区| 久热久热精品视频在线观看| 亚洲一区二区三区福利视频|