分数阶微分方程的高精度算法.docVIP

A Dissertation Submitted in Partial Ful?llment of the Requirements for the Degree of Doctor of Philosophy in Science High accuracy algorithms for fractional differential equations Ph.D. Candidate : Ma Xiaohua Major  : Computer software and theory Supervisor  : Prof.  Huang Cheng- ming Huazhong University of Science and Technology Wuhan 430074, P. R. China May, 2013 独创性声明 本人声明所呈交的学位论文是我个人在导师的指导下进行的研究工作及取得的 研究成果。尽我所知，除文中已标明引用的内容外，本论文不包含任何其他人或集 体已经发表或撰写过的研究成果。对本文的研究做出贡献的个人和集体，均已在文 中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。 学位论文作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，即：学校有权 保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。 本人授权华中科技大学可以将本学位论文的全部或部分内容编入有关数据库进行检 索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 本论文属于  保密 ，在 不保密 。  年解密后适用本授权书。 （请在以上方框内打“ ”） 学位论文作者签名：  指导教师签名： 日期： 年 月 日 日期： 年 月 日√ √ 摘  要 分数阶微分方程可用于模拟物理、力学、生物学、工程、金融、水文学、分数 阶控制器等领域中的许多现象。由于其深厚的物理背景，分数微分方程近年来已经 变成一个非常重要的热门课题。一般情况下，大多数分数阶微分方程的解析解是难 以获得的。所以，有必要开展其数值方法研究。在本文中，我们考虑几类分数阶常 微分方程和分数阶偏微分方程的高阶近似算法。 第一章，我们主要回顾分数阶微积分和分数阶微分方程数值分析的发展历史和 研究背景，并给出本论文的工作概要。同时为后文方便，我们给出分数阶导数的一 些定义和基本性质。 第二章我们研究一类非线性分数阶积分微分方程的数值解法。我们考虑混合配 置方法，即在第一个子区间使用非多项式配置，在剩余的子区间使用分级网格分段 多项式配置，并详细给出了该方法收敛阶的理论分析。 第三章我们针对一般的线性分数阶积分微分方程构造了一个Jacobi配置方法。 首先，我们使用一些函数和变量变换把方程转化成一个定义在标准区间[?1, 1]上 的Volterra积分方程。然后用Jacobi-Gauss点作为配置节点且用Jacobi-Gauss求积公式 近似方程中的积分，得到离散格式，且我们研究了该方法在L∞意义下的收敛阶。 第四章我们考虑Lane-Emden型分数阶常微分方程奇异初值问题。我们针对这 个方程构造了一个Jacobi配置方法。这个方法首先把这个Lane-Emden型的分数阶 常微分方程的求解转化成Volterra积分方程的求解。然后用Jacobi-Gauss点作为配置 节点且用Jacobi-Gauss求积公式近似方程中的积分。最后分别给出Jacobi配置方法 在L∞和L2范数下的误差分析。 第五章，我们讨论非线性加热下分数阶广义二阶流体Stokes第一类问题的数值 解法。我们在时间方向采用有限差分，在空间方向采用Legendre谱方法，获得高阶 有限差分谱方法，并给出方法的稳定性和L2范数下的收敛性分析。 关键词： 分数阶微积分; 配置方法; 谱方法; 有限差分方法; 收敛阶; 稳定性. I Abstract fractional differential equations have been proved to be a valuable tool in the model- ing of many phenomena，especially in physics, mechanics, biology, engineering, ?nance, hydrology, and fractional-order controllers．In recent years, due to its profound physical background, the subject of the fractional differential equations i