分数阶微分方程-课件

分数阶微分方程第三讲分数阶微分方程基本理论一、分数阶微分方程的出现背景及研究现状1、出现背景分数阶微积分是关于任意阶微分和积分的理论，它与整数阶微积分是统一的，是整数阶微积分的推广。整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受，很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题，其无论在理论分析还是数值求解方面都已有较完善的理论。但当人们进入到复杂系统和复杂现象的研究时，经典整数阶微积分方程对这些系统的描述将遇到以下问题：（1）需要构造非线性方程，并引入一些人为的经验参数和与实际不符的假设条件；（2）因材料或外界条件的微小改变就需要构造新的模型；（3）这些非线性模型无论是理论求解还是数值求解都非常繁琐。基于以上原因，人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程，其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势，因而成为复杂力学与物理过程数学建模的重要工具之一。2、研究现状在近三个世纪里，对分数阶微积分理论的研究主要在数学的纯理论领域里进行，似乎它只对数学家们有用。然而在近几十年来，分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。分数阶微积分理论也受到越来越多的国内外学者的广泛关注，特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。随着分数阶微分方程在越来越多的科学领域里出现，无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。然而由于分数阶微分是拟微分算子，它的保记忆性（非局部性）对现实问题进行了优美刻画的同时，也给我们的分析和计算造成很大困难。在理论研究方面，几乎所有结果全都假定了满足李氏条件，而且证明方法也和经典微积分方程一样，换句话说，这些工作基本上可以说只是经典微积分方程理论的一个延拓。对分数阶微分方程的定性分析很少有系统性的结果，大多只是给出了一些非常特殊的方程的求解，且常用的求解方法都是具有局限性的。在数值求解方面，现有分数阶方程数值算法还很不成熟，主要表现为：（1）在数值计算中一些挑战性难题仍未得到彻底解决，如长时间历程的计算和大空间域的计算等；（2）成熟的数值算法比较少，现在研究较多的算法主要集中在有限差分方法与有限单元法；（3）未形成成熟的数值计算软件，严重滞后于应用的需要。鉴于此，发展新数值算法，特别是在保证计算可靠性和精度的前提下，提高计算效率，解决分数阶微分方程计算量和存储量过大的难点问题，发展相应的计算力学应用软件成为迫切需要关注的课题。二、预备知识1、分数阶微积分经典定义回顾作为分数阶微积分方程的基础，本书在第二章中对分数阶微积分的定义及性质做了系统的介绍，为了接下来讨论的需要，我们首先对其进行一个简要的回顾。（1）分数阶微积分的主要思想如上图所示，分数阶微积分的主要思想是推广经典的整数阶微积分，从而将微积分的概念延拓到整个实数轴，甚至是整个复平面。但由于延拓的方法多种多样，因而根据不同的需求人们给出了分数阶微积分的不同定义方式。然而这些定义方式不仅只能针对某些特定条件下的函数给出，而且只能满足人们的某些特定需求，迄今为止，人们仍然没能给出分数阶微积分的一个统一的定义,这对分数阶微积分的研究与应用造成了一定的困难。（2）几种经典的分数阶微积分定义下面我们试图从理论依据、定义域、表达式和优缺点几个方面给出常见的四种分数阶微积分定义的比较图。依据1：整数阶微分的差分定义()001()lim(1)(),kkrkhrkftftrhkNrhG-L分数阶微积分定义：[]001[()]lim(1)()tanhGpratphrpDftftrhrh1()1[()]()(),0()tCpnpnataDfttufudupnp11[()]()(),0()tCppataDfttufudupp依据2：整数阶积分的柯西公式()11()()()()tnnafttufudunR-L分数阶微积分定义：()11[()]()()()tRnpnpataDfttufudunp()[()][()],1,nRpRnpatatdDftDftnpnnNdtCaputo分数阶微积分定义：先积分后微分先微分后积分缺点：定义域较窄计算复杂优点：定义域较宽缺点：Laplace变换较复杂优点：Laplace较简洁缺点：定义域太窄分数阶微积分定义定义域：(){()}ftft[a,b]在上逐段连续且在任意有限子区间上可积定义域：(){()}ftft[a,b][p]+1在具有阶连续导数定义域：(){()}ftft[a,b][p]+1在具有阶连续导数基于广义函数的分数阶微积分定义：()()()patpDftftt1,0()()0,0pptttpt定义域：{()()0;}ftftta优点：有利于工程中对系统的描述从上图我们看到，在分数阶微积分的发展过程中，人们根据不同的需求，从不同角度给出了分数阶微积分的定义，但这些定义无论从对象上还是从表达式上都无法实现统一，它们之间的关系大致可以用下图来表示。注：条件1：()ft在[,]ab上逐段连续，且在任何有限子区间上可积；条件2：()ft在[,]ab上具有[]1p阶连续导数；条件3：()()0kfa，1,2,,[]kp；条件4：()0ft，ta。由上图我们可以看到，对于不同的分数阶微积分定义方式有着不同定义域，即便是在公共区域内，不同的定义方式之间也无法实现完全的统一，这对分数阶微积分的应用和研究造成了一定的困难，因此人们迫切期望着分数阶微积分的一种哪怕是形式上的统一定义方式。2、M-R序列分数阶微分的定义为了满足实际需要，下面我们试图从形式上对分数阶微积分给出一种统一的表达式。分数阶微积分的主要思想是推广经典的累次微积分，所有推广方法的共同目标是以非整数参数p取代经典微积分符号中的整数参数n，即：nnddtppddt实际上，任意的n阶微分都可以看成是一列一阶微分的叠加：()()nnndftdddftdtdtdtdt（1）由此，我们可以给出一种在很多实际应用中十分重要的分数阶微积分的推广方式。首先，我们假设已有一种合适的推广方式来将一阶微分推广为（01）阶微分，即dDdt是可实现的。那么类似地可得到（1）的推广式为：()()nnDftDDDft（2）这种推广方式最初是由..KSMiller和.BRoss提出来的，其中D采用的是RL分数阶微分定义，他们称之为序列分数阶微分。序列分数阶微分的其他形式可以通过将D替换为GL分数阶微分、Caputo分数阶微分或其他任意形式分数阶微分来得到。进一步，如果我们将（2）中的分数阶微分D替换为不同阶数的分数阶微分可得到序列分数阶微分更一般的表达式：12()()nDftDDDft（3）12n根据问题的需要，D可以是RL分数阶微分、GL分数阶微分、Caputo分数阶微分或其他任意形式的分数阶微分，从这一点看来，我们可以说序列分数阶微分从形式上给出了分数阶微积分在时域上的一个统一表达式，RL分数阶微分、GL分数阶微分和Caputo分数阶微分都只是序列分数阶微分的一种特殊情况。故而，下面我们在对分数阶微积分方程进行理论分析的时候可以仅仅针对序列分数阶微积分来给出结论。3、M-R序列分数阶微分的Laplace变换下面我们考虑如下形式的序列分数阶微分的Laplace变换。11mmmatatatatDDDD（4）1111mmmatatatatDDDD（5）101,(1,2,,)mmjjjjm（6）在R-L分数阶微分定义下有：1000{();}()[()]tttLDftssFsDft（7）重复利用上式m次可得：110000{();}()[()]mmmmkmkmtttkLDftssFssDft（8）注：虽然上述序列分数阶微分的Laplace变换是在R-L分数阶微分定义下进行证明的，但是该结论对其他几种分数阶微积分也是成立的。4、泛函理论基础定理1（Schauder不动点定理）设U是Banach空间X的有界闭子集，如果:TUU是连续映射，那么T在U中存在不动点，即使得Txx的点存在。定义1（Lipschitz条件）设,Xd是距离空间，T是从X到X的映射，如果存在常数0q，使得对所有的,xyX，(,)(,)dTxTyqdxy则称T满足Lispschitz条件，q成为T的Lispschitz常数。特别的，如果1q，则T称为压缩映射。定理2（Banach压缩映像原理）设,Xd是距离空间，:TXX是压缩映射，则T在X中恰有一个不动点。设这个不动点为x，则对任何初始点0xX，逐次迭代点列1nnxTx，1,2,n收敛于x，且关于收敛速度有如下估计式：100(,)(1)(,)nndxxqqdTxx其中，q是T的Lipschitz常数。三、解的存在唯一性理论近年来，分数阶微分方程已经在国内外引起极大的研究兴趣，尤其是关于其解的性质的研究，诸如存在性及唯一性等，其中大多数的研究方法是通过把分数阶初值问题转换成等价的分数阶积分方程,然后运用不动点定理来得到分数阶初值问题解的存在唯一性结果。已有研究结果主要有以下限制：（1）函数的定义区间为有限区间[,]ab；（2）函数在定义域上需满足Lipschitz条件；因此，目前人们在这方面所做的工作都是希望设法在放宽上述两个限制条件后给出分数阶微积分方程的解的存在唯一性定理。下面我们对分数阶微分方程初值问题的现有理论结果作一个简单的介绍，相应的结论都是针对定义在有限区间[0,]T上的M-R序列分数阶微分形式，在满足Lipschitz条件下给出的，当然，由前面的介绍可知，这些结论也可直接推广到其他分数阶微分形式。1、线性分数阶微分方程解的存在唯一性定理考虑如下形式的初值问题：1001100()()()()()(),(0)(9)[()],1,2,,(10)njnkntjtnjttkDytptDytptytfttTDytbkn且1()(0,)ftLT，即0()Tftdt（11）第一步：假设()0,(1,2,,)kptkn，考虑由此得到的退化问题解的存在唯一性。定理1如果1()(0,)ftLT，则方程0()()ntDytft（12）有满足初值条件（10）的唯一解1()(0,)ytLT。定理的证明过程如下：步骤一通过Laplace变换证明解的存在性；下面我们设法构造一个待求解问题解，对式（12）做Laplace变换可得：11000()[()]()nnnknknttksYssDytFs（13）其中，()Ys、()Fs分别是()yt、()ft的Laplace变换。利用初值条件（10）可得：10()()nnknnkkYssFsbs（14）对上式做Laplace逆变换可得：11101()()()()()nitniinibyttfdt（15）步骤二由分数阶微分的线性性和Laplace变换的性质证明唯一性。假设有存在两个满足上述初值问题的解1()yt、2()yt令12()()()ztytyt，有分数阶微分方程的线性性可得：0100()0[()],1,2,,nktttkDztDytbkn（16）从而有()0Zs（17）由Laplace变换的性质可知：()0zt在(0,)T上几乎处处成立。故原方程的解在1(0,)LT上唯一。注：上述证明过程中用到的Laplace变换法是一种常用的分数阶微分方程求解方法，该