非线性空间分数阶Fisher方程的数值解法-厦门大学学报（自然科学版）.docVIP

非线性空间分数阶Fisher方程的数值解法 陈雪娟*，陈景华 （集美大学理学院，福建 厦门 361021） 摘要：分数阶微分方程是经典的整数阶微分方程的推广，这类方程在物理、工程和金融等领域得到了广泛的运用. 本文考虑非线性空间分数阶Fisher方程的数值解，提出一种基于二次多项式样条函数的数值解法，并证明该方法具有无条件稳定性和收敛性. 为了验证所构造格式的有效性，引入分数阶行方法 (FMOL) 与之进行比较. 最后通过一个数值算例说明本文的理论分析是正确的，所构造的离散格式是有效的. 关键词：分数阶扩散方程；Caputo分数阶导数；二次多项式样条函数；行方法 中图分类号：O 241.82 文献标识码：A 分数阶微分方程(时间、空间或时间空间分数阶)是传统整数阶微分方程的推广，即用分数阶导数代替相应的整数阶导数. 由于分数阶微积分具有记忆和遗传特性，目前已被广泛地应用于模拟工程、物理、化学和其它科学领域的许多现象[1-4]. 众所周知，要得到偏微分方程的解析解是很困难的，对于分数阶偏微分方程而言更是如此，比较典型和常用的方法是借助于各种积分变换和特殊函数的方法[5-7]. 因此，研究分数阶方程的数值模拟有着重要的理论与应用价值,近年来许多学者在这个领域里做了大量的研究工作[8-11]. 本文考虑非线性空间分数阶Fisher方程的数值解问题，Fisher方程主要用来模拟物种增长和扩散问题. 在已有的研究成果中，主要的数值解法是差分法、有限元法和谱方法等，但是采用二次多项式样条函数进行数值逼近的研究文献却较缺乏[12,13]. 本文提出一种基于二次样条函数的数值解法，并分析所构造迭代格式的稳定性和收敛性. 1 一些记号和简单的结论 非线性空间分数阶Fisher方程如下： (1) 收稿日期：2015-07-05 录用日期：2015-10-25 基金项目：福建省课题（）资助 *通信作者： 初始条件和边界条件分别为： (2) (3) 其中，和为充分光滑的已知函数，非线性源项是关于满足 Lipschitz 条件. 阶的Caputo分数阶导数定义如下： (4) 阶的Riemann-Liouville分数阶导数定义如下： (5) 以上两种分数阶导数的定义之间具有如下关系式[14]： . (6) 在实际求解微分方程初值问题的过程中，Caputo导数比Riemann-Liouville导数的应用更为广泛而且更具有物理背景[15]. 本文采用Caputo分数阶导数定义，即 记，，其中，分别是空间步长和时间步长(和是正整数). 设点集，表示的近似解. 本文首先利用二次多项式样条函数提出一种逼近Fisher方程数值解的迭代格式，并证明该格式具有无条件稳定性和收敛性. 然后给出非线性空间分数阶Fisher方程的分数阶行方法(FMOL)，用于验证所构造格式的有效性. 最后通过一个数值算例说明本文理论分析的的正确性和可行性. 2 基于多项式样条函数的数值解法 2.1 二次多项式样条函数 我们考虑如下二次多项式样条函数: (7) 为了确定函数的系数表达式，首先定义 (8) (9) (10) 由等式(7) (8)和(9)，可得 (11) 其中 又由等式(10)和Caputo分数阶导数的定义，有 (12) 因此，系数的表达式为 (13) 2.2 数值解法的迭代格式 为了使二次函数在处满足连续性条件： ,, 由等式(7)和(13)，可以推出 (14) 其中，, 引理1 设迭代格式(14)的截断误差为，则截断误差满足: (15) 证明 由于 , (16) 将等式(16)的右边在点处进行Taylor展开，得到 . (17) 证毕. 所以，迭代格式(14)的截断误差阶为，我们得到 (18) 利用向后差分法 , (19) 和等式 , (20) 由方程(