论文部分内容阅读
本文的主题是研究分数阶微分方程的数值方法及其稳定性和收敛性。基于不同的数值方法,论文的主体分为三部分:求解分数阶微分方程的有限差分方法,有限元方法和多项式配置方法。论文的绪论介绍了课题的背景及研究意义,给出了分数阶导数的基本定义,回顾了分数阶微分方程数值方法的历史及研究现状,概括了论文的主要内容。论文的第二章研究了一种分数阶偏微分方程的有限差分方法,该方法首先由Diethelm在1997年提出,并成功用于求解简单的单变量左侧分数阶导数问题。在本文中,将该方法推广到右侧分数阶导数问题,并进一步推广到求解分数阶偏微分方程的初边值问题上,给出了具体的离散格式。进而利用代数矩阵方法对此离散格式进行了稳定性分析,给出了关于网比的稳定性条件以及收敛阶估计。最后,数值试验验证了方法的稳定性条件及收敛阶。论文主体的第二部分研究了求解不同分数阶偏微分方程的有限元方法,由三,四,五章组成。在第三章中研究了一类时间空间Riesz分数阶对流-弥散方程的有限元方法。时间空间Riesz分数阶对流-弥散方程的特点是具有时间方向的左侧Caputo型分数阶导数和空间方向的Riesz分数阶导数,本章研究了此类方程的初边值问题的数值方法。首先基于实数阶Sobolev空间的定义给出了方程的变分形式,利用Lax-Milgram定理证明了弱解的存在唯一性。然后,时间方向上采用有限差分法进行离散,空间方向上采用有限元方法进行离散,讨论了方法的收敛阶,随后证明了全离散格式的无条件稳定性。最后用数值算例验证了方法的理论结果。第四章研究了一类时间复合型分数阶偏微分方程的有限元方法,所研究的方程模型特点是在时间方向上具有复合型的左侧Caputo型分数阶导数。针对该类方程类似上一章的方法,给出了此类方程的弱形式,并证明了弱解的存在唯一性。进而应用Diethelm的分数阶离散化方法,给出了时间方向上的复合型分数阶导数的有限差分离散格式并证明了方法的收敛阶。在空间方向上应用有限元方法进行离散,并得到收敛阶。进而结合时间方向和空间方向的离散化格式,给出了方法无条件稳定的证明。最后用数值算例验证了方法的收敛性,并和几种已有的分数阶数值方法进行了比较,验证了方法的有效性。第五章研究了一类时间空间复合型Riesz分数阶对流-扩散方程的有限元方法。此类方程的模型具有时间方向上的左侧Caputo分数阶复合型微分算子以及空间方向上的Riemann-Liouville型Riesz分数阶微分算子。理论方面,首先给出了此类方程的弱形式,并证明了弱解的存在唯一性。数值方法方面,分别就时间和空间方向进行了离散化,其中时间方向利用差分法进行离散给出了收敛阶,空间方向用有限元方法进行离散给出了方法的收敛阶,随后证明了全离散格式的无条件稳定性,最后给出了具体的数值算例验证了方法的有效性。本文主体的第三部分,第六章讨论了具有弱奇核的分数阶积分微分方程的配置方法。首先,证明了具有弱奇核的分数阶积分微分方程等价于一类具有弱奇性的线性Volterra积分方程。其次,利用分片多项式配置方法求解该Volterra积分方程。然后证明了当精确解充分光滑时,用一致网格构造的配置方法能够达到最优的收敛阶;当真解不够光滑时,可以通过变换的方法使得对变换后的方程应用配置法也能得到最优的收敛阶。最后在数值实验中,考虑了梯形配置方法,Simpson配置方法,Newton3/8配置方法以及Cotes配置方法的收敛性。