《一种Caputo型时间分数阶波动方程的差分方法.》.pdfVIP

2014年 7月 四川师范大学学报(自然科学版) July，2014 第37卷 第4期 JournalofSichuanNormalUniversity(NaturalScience) Vo1．37．No．4 一 种 Caputo型时间分数阶波动方程的差分方法 余跃玉 (四川文理学院数学与财经学院，四川达州 635000) 摘要：将带阻尼项的波动方程中的阻尼项和对时间的二阶导数，用Caputo分数阶导数替换，从而得到 一 个带Caputo分数阶阻尼项的分数阶波动方程．对该方程，建立了一种差分格式，证明了此格式差分解的 存在唯一性 ，分析了差分解的收敛性和稳定性，并用数值试验验证了格式的有效性． 关键词：波动方程；Caputo分数阶导数；有限差分格式；收敛性 ；稳定性 中图分类号：O241．4 文献标志码：A 文章编号：1001—8395(2014)04—0524—05 doi：10．3969／j．issn．1001—8395．2014．04．016 1 引言及预备知识 +A = + ， 随着科技的发展，计算能力的提高，分数阶微 用分数阶Caputo导数：D (，t)去替换二阶时间导 分方程不仅越来越多的应用于光学、热学、材料学、 数构造了一种差分格式，并对格式进行研究．本文 力学、信号处理和辨别、控制学、生物学、金融及其 将在此基础上，对其阻尼项的一阶导数也用Caputo 其它领域 J，而且能更好地拟合某些 自然物理过 导数：D (，t)进行替换，再用有限差分的方 程和动态系统过程，特别是分数阶时间导数模型能 法 ¨-19]对带分数阶阻尼项的时间分数阶波动方程 用较少的参数精确地拟合实验数据，因此分数阶微 进行数值求解． 分方程引起了广大学者的关注．然而，分数阶微分 考虑如下带分数阶阻尼项的时间分数阶波动 方程的解析解往往形式复杂，很难用简单的函数表 方程初边值问题 示而且有些非线性的方程的解析解是不易求得的， DM(，)+AD (，t)= 于是研究分数阶微分方程的数值解引起了广大学 者的关注l6 j． +z／s( )， 分数 阶微分算子主要有 3种形式的定义： 0≤ ≤L， 0tT， (2) Grumwald—Letnikov定义、Riemann—Liouville定义 (，0)= )， 和Caputo定义，其中分数阶Caputo导数：D “(，) ：g()，o≤ ≤L， (3) 定义如下 1 ，t Z／(0，t)= “(L，t)=0， 0≤t≤T， (4) nC (，) ：南 J。(一)一一× 其中，1OL≤2，0 ≤1，A0， 是常数，f(0)= — ，(L)=0 )和g()是性质充分好的实函数． d ， n一1 n， (1) 如果 =2，／3=1，则方程 (2)称为电报方程 ，它 其中，，为Gamma函数． 描述了电子在电报电缆中的传输过程，也可