一类非线性R-L分数阶积分微分方程的数值解法.pdfVIP

河北科技师范学院学报第29卷第2期，2015年6月 ofHebeiNormM of Vd．29No．2Jun．2015 JoumM UnivemityScience＆Techndogy DOI：10．3969／J．1SSN．1672-7983．2015．02．010 一类非线性R—L分数阶积分微分方程的数值解法 张盼盼1’2，任正杰2 (1宁夏大学数学计算机学院，宁夏银川，750021；2甘肃省永昌县第一高级中学) 摘要：利用Adomian多项式将分数阶积分微分方程中的积分项离散化，进而得到原方程解的级数表达形式， 数值算例验证了该分解方法的有效性。 关键词：Adomian多项式；分数阶积分微分方程；数值解 中图分类号：0241。8 文献标志码：A 文章编号：1672-7983(2015)02-0047-05 考虑分数阶微积分，但由于缺乏应用等众多原因的支撑，使得分数阶微积分的发展相对滞后。近年来， 越来越多的分数阶微积分方程广泛出现在各大工程领域。例如，PID控制理论、粘弹性材料及混沌现象 等反常问题¨’2J，其中有些问题建模后得到的方程大多数都是非线性分数阶积分微分方程。由于分数 阶微积分复杂的定义，增加了分数阶方程的求解难度，分数阶积分微分方程的求解对众多学者来说更是 挑战。目前关于分数阶积分微分方程数值解的文献已有一些，主要是用小波的方法求解Fredholm或 Voherra型的分数阶积分微分方程¨“J。LI 微分方程。然而，关于非线性分数阶积分微分方程数值解的文献甚少。 利用Adomian分解法求解代数方程具有强大的优越性，其解为级数形式、收敛快、计算方便且容易 在计算机上运算。它的核心思想是：把方程拆成尼个部分，把解拆成忍个项，同时非线性项用一种特殊 的多项式A。进行替代，然后由低阶解分量逐渐向高阶解分量逐一解出，从而得到方程的近似解析解，也 可以得到精确解。该方法无需进行任何变换，很好的保持了原方程的物理性质。在实际生活中不仅会 遇到线性问题和确定问题，而且经常会遇到非线性问题甚至是随机问题。对于这些非线性问题，通常的 办法是将其线性化近似处理或采用扰动技术。然而，这些方法在求解时无形中改变了原问题，得到的解 往往不能满足实际需要。Adomian分解法不会存在这样的问题，计算量很小，而且很快就能收敛到真 解。所以，Adomian分解法是求解非线性方程的有力工具旧’9J，笔者利用Adomian多项式将R．L分数阶 积分微分方程离散化，得到原方程解的级数表达式，并使用数值算例证明该算法的有效性。 1预备知识 定义‘11设火戈)∈L[a，b]，ao。则称 agf(￡)2高J。(t—r)”狄r)dr (1) 为Riemann-Liouville分数阶积分，其中t“口，b]，F(a)为Gamma函数。 性质…2’ (I)若a，卢o，八f)E c，则有：。片。JF(t)=。r+袄￡)=。刀。gf(t)。 (2)设八￡)=^秽一l，则有：。刀￡”=下黼∥4。 2主要结果 设 D4菇(t)=A J K(f，r九厂(省(7-))dr+g(z) (2) 收稿日期：2015-04-02；修改稿收到日期：2015-05．24 万方数据 河北科技师范学院学报 29卷 程(2)两端作用算子C，则有 (3) 础)=击小t刊川咖)dr+击小t刊川，：附∽八小))dsd丁 由Adomian分解法知，设 八t)=静。(菇。(￡)，菇。(￡)，…，戈。(￡)) (4) 其中