分数阶微分方程数值解的一种逼近方法讲解.pdfVIP

分数阶微分方程数值解的一种逼近方法 By:Pankaj Kumar, Om Prakash Agrawal 摘要 本文提出了一类分数阶微分方程（FDEs）的数值解方案.在这种方法中,FDEs 被Caputo 型分数阶导数所表现. Caputo 型分数阶导数的属性可以让一个分数阶 微分方程减弱为一个 Volterra 型积分方程. 这样做了之后,许多研究 Volterra 型积分方程的数值方法也可以应用于寻找FDEs 的数值解. 本文总时间被划分为 一组小区间,在两个连续区间中，用二次多项式逼近未知函数. 这些近似被替换 成转化的Volterra 型积分方程由此获得一组方程. 这些方程的解提供了FDE 的 解. 这种方法被应用于解决两种类型的 FDEs,线性和非线性. 用这里给出的方 法得到的解能与解析解和其他方法的数值解较好的吻合 . 同时结果说明这种数 值方法是稳定的. 1.引言 本文讨论分数阶微分方程的数值解. 分数阶导数和分数阶积分近年来收到 了广泛的关注. 在许多实际应用中，分数阶导数和分数阶积分为考虑的系统提供 了更加精确地模型. 比如，分数阶导数已经被成功地运用到模拟许多粘性材料的 依赖频率的阻尼行为.1980 年之前，Bagley 和Torvi 提出了这个领域已经被研 究的工作的一个回顾，并且说明了半阶导数模型可以非常好地描述阻尼材料的频 率以来. 另一些学者说明了分数阶导数和分数阶积分在电化学过程，电解质极 化，有色噪声，粘性材料和混沌领域的应用. Mainardi，Rossikhin 和Shitikova 提出了分数阶导数和分数阶积分在一般固体力学，特定粘弹性阻尼模型中的应用 的调查. Magin 提出了分数阶微积分在生物工程的三个关键部分的回顾. 分数阶 导数和分数阶积分在其他领域的应用以及相关的数学工具和技巧还可以在许多 其他文献上找到. 系统模型中分数阶导数的引进大多会导致分数阶微分方程的出现 . 对某些 特定的分数阶微分方程在通常系统条件下的解，已经有几种方法被找到. 这些方 法包括，拉普拉斯变换，傅里叶变换，模态综合法和特征向量展开法，数值法以 及基于Laguerre 积分公式的方法. 然而，这些方法中大多数不能被应用到非线 性分数阶微分方程. 更进一步的，正如Diethelm 等人指出的，这些方法很多只 能应用到特定类型的分数阶微分方程，并且人们并不知道他们能否被推广 . 并 且，在很多作者的研究成果中，并没有出现系统性的收敛性分析. 最近，对于能被应用到线性和非线性分数阶微分方程的数值稳定数值逼近技 巧，人们的兴趣愈发浓厚. 这些方法技巧大多利用了分数阶微分方程可以被减弱 为Volterra 型积分方程的特性. 因此，Volterra 型积分方程的数值解法也可以 应用到分数阶微分方程的解当中. Diethelm 等人提出了分数阶微分方程数值解 的一种PECE 方法，其中P,C,E 分别代表预测，校正和估计. 这样一来很多学者 又推广了应用于常微分方程和分数阶微分方程的 Adams –Bashfort –Moulton 型预测-校正格式. 这种方法的提出也是利用分数阶微分方程可以被转化为 Volterra 型积分方程的特点. 这些作者同时提出了误差分析和用Richardson 外 推法改善数值精度的延伸. Ford 和Simpson 提出了一种阶数大于 1 的分数阶微 分方程的数值解法. 在该公式中，阶大于1 的分数阶微分方程被减弱为阶小于1 的分数阶微分方程，然后用相应的数值解法解由此导出的系统. 在所有这些方法 当中，节点之间的未知函数用线性函数逼近. Kumar and Agrawal提出了阶数大 于 1 的分数阶微分方程的数值解法. 这种方法要求就 y(t)和它的导数在时间节 点上连续. 本文基于古典分数阶微分方程可以转化为Volterra 型积分方程的特点也提 出了一种数值方法来逼近分数阶微分方程的解. 特别地，我们用二次逼近函数来 建立这种算法，结果说明这种方法可以被应用到寻求分数阶微分方程的数值解. 我们还通过两个例子，线性和非线性问题的解决，说明了这种方法的高效和准确， 并且这种数值方法是稳定的. 2.数值算法 关于分数阶导数的定义已经出现有好几种，它们包括Riemann –Liouville, Grun-wald –Letnikov, Wey