分数阶扩散-波动方程数值求解

摘要

针对Caputo分数阶导数意义下的时间分数阶扩散-波动方程进行数值研究.利用Caputo分数阶导数与Grunwald-Letnikov分数阶导数的关系对时间分数阶导数进行时间离散化处理,再利用二阶中心差商离散方程中的二阶空间导数,并结合边值条件的离散化,把离散化方程的求解转化为一个线性方程组的求解.利用Mat-lab编程实现了离散化方程组的求解,并绘制了不同参数下的数值解曲面.