非线性分数阶微分方程的hp,型Legendre,谱配置法
时间：2023-01-18 13:55:08 来源：雅意学习网 本文已影响 人 
李 珊，安 筱，孙桂磊
（上海理工大学 理学院，上海 200093）
分数阶微分方程最重要的是非局部性质,它能有效地避免整数阶导数的局部性,所以，分数阶导数可以更加精确地描述对历史有依赖的问题。这种非局部性质使得分数阶微分方程在许多工程和科学领域，特别是在描述具有记忆过程、遗传特性[1-3]和异质材料等方面越来越有吸引力。分数阶微分方程已被广泛应用于各种复杂的黏性流体流动模型，如异常扩散过程[4]、分数动力学[5]、粘弹性材料[6]、输送管道边界层效应[7]、氧气通过毛细血管输送到组织[8]等。同时，利用分数阶导数描述问题时，使用少量的参数就能获得与实际相吻合的结果，所以，与整数阶微分方程相比较，在处理复杂系统时，使用分数阶微分方程建模将更加简单、参数意义更加清楚。一般来说，对分数阶微分方程的求解是非常困难的。近几年，学者们对分数阶积分微分方程的基本理论及其数值求解进行了大量的研究。因为，一方面分数阶微分方程的求解可以转化为等价的积分方程进行求解；
另一方面为使得数学模型更加符合客观实际，对分数阶积分微分方程的求解已经不可避免，所以，对分数阶积分微分方程的理论及其数值方法的研究越来越受到学者们的关注。本文提出了一种求解以下非线性分数阶初值问题的有效数值方法。
其中，α∈(0,1),f∈C(0,T),并且 Dα是根据下面定义的 α阶Caputo 导数。
Caputo 导数 Dα不同于传统的微分算子，它是一种全局算子。目前分数阶导数主要分为3 种类型：Riesz 意义下的分数阶导数、Riemann-Liouville分数阶导数和Caputo 型分数阶导数。1832 年Liouville 采用级数的形式给出了分数阶导数，1853年Riemann 采用定积分的形式给出了另一种分数阶微分，Grünwald 和Krug 将Riemann 和Liouville的结论进行了统一，得到了Riemann-Liouville 分数阶导数的定义。1967 年Caputo 提出了Caputo 型分数阶导数。Caputo 型分数阶导数因其广泛的应用背景和相对简便的运算被广泛地应用于各个领域。目前在分数阶微分方程理论的研究中，关于其数值解进行了大量的研究。在Caputo 型分数阶定义下的初值条件具有明确的物理意义，因而在工程物理学领域该方程具有重要的应用价值。但由于解的复杂性，不能直接应用于实际工程中，所以，求解Caputo 型分数阶微分方程的数值解受到广泛关注。
谱方法作为一种全球性的高精度数值方法，在求解微分方程过程中扮演着重要的角色，非常适用于求解分数阶微分方程。然而，现有的谱方法主要是一步法。这种方法要求解具有高度正则性。然而，通常对于分数阶问题，解不可能是任意光滑的，即使源项非常光滑。因此，一步法缺乏有效处理局部弱奇异的能力。最近，Guo 等[9]提出了一种求解单分数阶边值问题的hp 型Chebyshev谱配置方法，并从理论上和数值上得到了hp 收敛性。在他们工作的基础上，本文提出了一种求解非线性分数阶初值问题的hp 型Legendre 谱配置法，该方法灵活地采用时间步长和局部逼近阶数，可以实现任意精度的数值解。并且，本文所提出的数值格式对于数值实现和误差分析来说更加简单易行。
设Ih是区间I=[0,T]上的一个网格,Ih:={tn:0=t0＜t1＜···＜tN=T}，并且设hn=tn-tn-1，In=(tn-1,tn]。通常情况下，对于一个给定区间 Λ和某一权函数 χ(x)，定义
设Lk(x),x∈(-1,1)是标准的阶数为k的Legendre多项式。Legendre 多项式集合构成一个完备的L2(-1,1)正交系统，即
其中，δk,j是Kronecker 符号。
定义k阶变换的Legendre 多项式
Ln,k(t),k≥0的集合构成一个完备的L2(In)正交系统，即
根据标准Legendre-Gauss 求积的性质，可得
正如文献[10]中提到的，由于弱奇点的存在(e.g.,ν∈(0,1))，考虑到加权插值求积公式的权值取决于核中的弱奇异因子(t-s)ν是很自然的,对于任意的 φ(s)∈,引入加权插值求积公式
方程(1)可以转换为以下等价的Volterra 积分方程[12]：
本文的主要目的是基于等价式(13) 提出并分析一种有效的求解方程(1)的 hp型Legendre-Gauss谱配置法。
设un(t)是第n个区间上式(13)的解，即
un(t):=u(t),t∈In,1≤n≤N
由式(13)可得，对于任意t∈In，
求解式(14)的hp 型谱配置法的目的是寻找Un(t)∈，使得
其中，Uk(t)∈是区间Ik上uk(t)的数值解。用U表示全局数值解,即
U(t)=Uk(t),t∈Ik,1≤k≤N
现在描述数值实现,并给出了式(15)的算法。设
由式(6)和式(9)可得，对于 1≤k≤n-1,
其中，
由式(6)，(8)和(11)可得
其中，
由式(15)～(20)可得
根据 Legendre 多项式的正交性，比较式(11)的展开系数得到，对于 0≤p≤Mn，
可采用迭代法计算系统式（22）的展开系数（例如，Newton-Raphson 迭代）。
现提供一系列数值实验来验证所提出算法的有效性。首先，定义在网格点上离散的L2误差和L∞误差如下：
例1考虑问题
问题(23)的真解是u(t)=tα。
在图1（a）和1（b）的第1 个图中，分别列出了问题（23）在相同的M（多项式阶数）以及变化的h（网格长度）下的数值误差（h型）和，其中，选取M=1，2，4，6。由此可以观察到h 型的快速收敛速度。
在图1(a)和1(b)的第2 个图中，分别列出了问题(23)在相同的h（网格长度）以及变化的M（多项式阶数）下的数值误差，其中，选取h=1，，由此可以观察到数值误差的收敛速度随着M的减小按代数衰减。
