论文部分内容阅读
时间分数阶微分方程是微分方程研究领域的一个重要分支,在各类扩散系统中有着广泛的应用。但是由于时间分数阶微分方程具有分数阶导数,不易求解,从现有的研究方法看,求解此类微分方程的主要方法是将其离散化进而求解数值解。由于B样条函数是对称单峰值函数,具有光滑性好、紧支集等特点,同时保留了分段插值多项式计算简单、稳定性好、收敛性有保证且易于在计算机上实现的优点。并且B样条配点法构造简单,数值精度高,易于处理复杂的边界问题,目前已成为求解偏微分方程的重要数值方法之一。所以本文提出使用三次B样条方法求解时间分数阶微分方程数值解以及参数反演,研究其求解此类问题的收敛性和稳定性。为了研究三次B样条方法能否较好地求解时间分数阶微分方程数值解,本文将利用Caputo分数阶微分定义,将时间分数阶微分转变为积分形式,同时对于时间维度的一阶偏导采用一阶差分格式,离散化时间分数阶微分方程。实验结果表明,三次B样条方法能够有效快速地求解时间分数阶微分方程,并且能够将误差控制在10-5数量级。随后,本文根据三次B样条插值方法的误差,通过理论证明来验证三次B样条方法求解时间分数阶微分方程的有效性。通过证明我们得到三次B样条方法求解时间分数阶微分方程的收敛阶为(Δt2-α+Δx2),进一步证明三次B样条方法求解时间分数阶微分方程的有效性。源项反演也是微分方程数值解的一个重要部分,在实际问题中通常需要通过观测值来推测源项、参数或者初边值,由于时间分数阶微分方程的特殊性,当前的分数阶微分与之前所有时间的函数值相关,所以很难使用整数阶微分方程的方法求解,因此本文提出了一种新的泛函,对方程进行平方计算并在时间空间上积分产生一个泛函,极小化该泛函所得到的源项即所求源项。但是仅仅对泛函偏微分等于零求源项具有不适定性,因此本文在这里引入正则化方法,然后通过最小二乘法的思想求解参数方程。实验结果表明,改进的三次B样条源项反演算法虽然因为三次B样条函数的分段特点增加了源项中的参数个数,但是算法不需要通过迭代,避免了复杂的计算,拟合效果整体上要优于多项式拟合算法。