2023年分数阶微分方程-课件.pdf

《2023年分数阶微分方程-课件.pdf》由会员分享，可在线阅读，更多相关《2023年分数阶微分方程-课件.pdf（15页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精品资料 欢迎下载 分数阶微分方程 第三讲 分数阶微分方程基本理论 一、分数阶微分方程的出现背景及研究现状 1、出现背景 分数阶微积分是关于任意阶微分和积分的理论，它与整数阶微积分是统一的，是整数阶微积分的推广。整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受，很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题，其无论在理论分析还是数值求解方面都已有较完善的理论。但当人们进入到复杂系统和复杂现象的研究时，经典整数阶微积分方程对这些系统的描述将遇到以下问题：（1）需要构造非线性方程，并引入一些人为的经验参数和与实际不符的假 设条件；（2）因材料或外界条件的微小改变

2、就需要构造新的模型；（3）这些非线性模型无论是理论求解还是数值求解都非常繁琐。基于以上原因，人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程，其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势，因而成为复杂力学与物理过程数学建模的重要工具之一。2、研究现状 在近三个世纪里，对分数阶微积分理论的研究主要在数学的纯理论领域里进行，似乎它只对数学家们有用。然而在近几十年来，分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问

3、题。分数阶微积分理论也受到越来越多的国内外学者的广泛关注，特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。随着分数阶微分方程在越来越多的科学领域里出现，无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。然而由于分数阶微分是拟微分算子，它的保记忆性（非局部性）对现实问题进行了优美刻画的同时，也给我们的分析和计算造成很大困难。在理论研究方面，几乎所有结果全都假定了满足李氏条件，而且证明方法也和经典微积分方程一样，换句话说，这些工作基本上可以说只是经典微积分方程理论的一个延拓。对分数阶微分方程的定性分析很少有系统性的结果，大多只是给出了一些非常特殊的方程的求解，且

4、常用的求解方法都是具有局限性的。在数值求解方面，现有分数阶方程数值算法还很不成熟，主要表现为：（1）在数值计算中一些挑战性难题仍未得到彻底解决，如长时间历程的计算和大空间域的计算等；（2）成熟的数值算法比较少，现在研究较多的算法主要集中在有限差分方法与有限单元法；（3）未形成成熟的数值计算软件，严重滞后于应用的需要。精品资料 欢迎下载 鉴于此，发展新数值算法，特别是在保证计算可靠性和精度的前提下，提高计算效率，解决分数阶微分方程计算量和存储量过大的难点问题，发展相应的计算力学应用软件成为迫切需要关注的课题。二、预备知识 1、分数阶微积分经典定义回顾 作为分数阶微积分方程的基础，本书在第二章中对

5、分数阶微积分的定义及性质做了系统的介绍，为了接下来讨论的需要，我们首先对其进行一个简要的回顾。（1）分数阶微积分的主要思想 如上图所示，分数阶微积分的主要思想是推广经典的整数阶微积分，从而将微积分的概念延拓到整个实数轴，甚至是整个复平面。但由于延拓的方法多种多样，因而根据不同的需求人们给出了分数阶微积分的不同定义方式。然而这些定义方式不仅只能针对某些特定条件下的函数给出，而且只能满足人们的某些特定需求，迄今为止，人们仍然没能给出分数阶微积分的一个统一的定义,这对分数阶微积分的研究与应用造成了一定的困难。（2）几种经典的分数阶微积分定义 下面我们试图从理论依据、定义域、表达式和优缺点几个方面给出

6、常见的四种分数阶微积分定义的比较图。依据 1：整数阶微分的差分定义()001()lim(1)(),kkrkhrkftf trh kNrhG-L 分数阶微积分定义：001()lim(1)()t anhGpratphrpDf tf trhrh 1()1()()(),0()tCpn pnataDf tt ufu du pnp 11()()(),0()tCppataDf tt uf u du pp依据 2：整数阶积分的柯西公式()11()()()()tnnaftt uf u dunR-L 分数阶微积分定义：()11()()()()tRn pn pataDf tt uf u dunp ()()(),1,

7、nRpRn patatdDf tDf tnp n n Ndt Caputo分数阶微积分定义：先积分后微分先微分后积分缺点：定义域较窄计算复杂优点：定义域较宽缺点：Laplace 变换较复杂优点：Laplace 较简洁缺点：定义域太窄分数阶微积分定义定义域：()()f tf ta,b在上逐段连续且在任意有限子区间上可积定义域：()()f tf ta,bp+1在具有阶连续导数定义域：()()f tf ta,bp+1在具有阶连续导数基于广义函数的分数阶微积分定义：()()()patpD f tf tt1,0()()0,0pptttpt 定义域：()()0;f t f tta优点：有利于工程中对系统的

8、描述 从上图我们看到，在分数阶微积分的发展过程中，人们根据不同的需求，从不同角度给出了分数阶微积分的定义，但这些定义无论从对象上还是从表达式上都无法实现统一，它们之间的关系大致可以用下图来表示。作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 注：条件 1：()f t在,a b上逐段连续，且在任何有限子区间上可积；条件 2：()f t在,a b上具有 1p 阶连续导数；条件 3：()()0kfa，1,2,kp；条件 4：()0f t，ta。

9、由上图我们可以看到，对于不同的分数阶微积分定义方式有着不同定义域，即便是在公共区域内，不同的定义方式之间也无法实现完全的统一，这对分数阶微积分的应用和研究造成了一定的困难，因此人们迫切期望着分数阶微积分的一种哪怕是形式上的统一定义方式。2、M-R 序列分数阶微分的定义 为了满足实际需要，下面我们试图从形式上对分数阶微积分给出一种统一的表达式。分数阶微积分的主要思想是推广经典的累次微积分，所有推广方法的共同目标是以非整数参数p取代经典微积分符号中的整数参数n，即：nnddt ppddt 实际上，任意的n阶微分都可以看成是一列一阶微分的叠加：()()nnnd f td ddf tdtdt dtdt

10、 （1）由此，我们可以给出一种在很多实际应用中十分重要的分数阶微积分的推广方式。首先，我们假设已有一种合适的推广方式来将一阶微分推广为（01）阶微分，即dDdt是可实现的。那么类似地可得到（1）的推广式为：()()nnDf tD DDf t （2）这种推广方式最初是由.K S Miller和.B Ross提出来的，其中D采用的是作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 RL分数阶微分定义，他们称之为序列分数阶微分。序列分数阶微分的其

11、他形式可以通过将D替换为GL分数阶微分、Caputo分数阶微分或其他任意形式分数阶微分来得到。进一步，如果我们将（2）中的分数阶微分D替换为不同阶数的分数阶微分可得到序列分数阶微分更一般的表达式：12()()nD f tD DDf t （3）12n 根据问题的需要，D可以是RL分数阶微分、GL分数阶微分、Caputo分数阶微分或其他任意形式的分数阶微分，从这一点看来，我们可以说序列分数阶微分从形式上给出了分数阶微积分在时域上的一个统一表达式，RL分数阶微分、GL分数阶微分和Caputo分数阶微分都只是序列分数阶微分的一种特殊情况。故而，下面我们在对分数阶微积分方程进行理论分析的时候可以仅仅针对

12、序列分数阶微积分来给出结论。3、M-R 序列分数阶微分的 Laplace 变换 下面我们考虑如下形式的序列分数阶微分的 Laplace 变换。11mmmatatatatDDDD （4）1111mmmatatatatDDDD （5）101,(1,2,)mmjjjjm （6）在 R-L分数阶微分定义下有：1000();()()tttLD f tss F sDf t （7）重复利用上式m次可得：110000();()()mmmm km kmtttkLDf tssF ssDf t （8）注：虽然上述序列分数阶微分的 Laplace 变换是在 R-L分数阶微分定义下进行证明的，但是该结论对其他几种分数阶

13、微积分也是成立的。4、泛函理论基础 定理 1（Schauder 不动点定理）设U是Banach空间X的有界闭子集，如果:T UU是连续映射，那么T在U中存在不动点，即使得Txx的点存在。定义 1（Lipschitz条件）设,X d是距离空间，T是从X到X的映射，如果存在常数0q，使得对所有的,x yX，作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载(,)(,)d Tx Tyqd x y 则称T满足Lispschitz条件，q成为T的Lis

14、pschitz常数。特别的，如果1q，则T称为压缩映射。定理 2（Banach 压缩映像原理）设,X d是距离空间，:TXX是压缩映射，则T在X中恰有一个不动点。设这个不动点为x，则对任何初始点0 xX，逐次迭代点列1nnxTx，1,2,n 收敛于x，且关于收敛速度有如下估计式：100(,)(1)(,)nnd xxqqd Txx 其中，q是T的Lipschitz常数。三、解的存在唯一性理论 近年来，分数阶微分方程已经在国内外引起极大的研究兴趣，尤其是关于其解的性质的研究，诸如存在性及唯一性等，其中大多数的研究方法是通过把分数阶初值问题转换成等价的分数阶积分方程,然后运用不动点定理来得到分数阶初

15、值问题解的存在唯一性结果。已有研究结果主要有以下限制：（1）函数的定义区间为有限区间,a b；（2）函数在定义域上需满足Lipschitz条件；因此，目前人们在这方面所做的工作都是希望设法在放宽上述两个限制条件后给出分数阶微积分方程的解的存在唯一性定理。下面我们对分数阶微分方程初值问题的现有理论结果作一个简单的介绍，相应的结论都是针对定义在有限区间0,T上的 M-R序列分数阶微分形式，在满足Lipschitz条件下给出的，当然，由前面的介绍可知，这些结论也可直接推广到其他分数阶微分形式。1、线性分数阶微分方程解的存在唯一性定理 考虑如下形式的初值问题：1001100()()()()()(),(

16、0)(9)(),1,2,(10)njnkntjtnjttkDy tp tDy tp t y tf ttTDy tbkn 且1()(0,)f tLT，即 0()Tf t dt （11）第一步：假设()0,(1,2,)kp tkn，考虑由此得到的退化问题解的存在唯一性。定理 1 如果1()(0,)f tLT，则方程 0()()ntDy tf t （12）有满足初值条件（10）的唯一解1()(0,)y tLT。定理的证明过程如下：步骤一 通过 Laplace 变换证明解的存在性；作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需

17、要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 下面我们设法构造一个待求解问题解，对式（12）做 Laplace 变换可得：11000()()()nnn kn knttks Y ssDy tF s （13）其中，()Y s、()F s分别是()y t、()f t的 Laplace 变换。利用初值条件（10）可得：10()()nn knn kkY ssF sbs （14）对上式做 Laplace 逆变换可得：11101()()()()()nitniiniby ttfdt （15）步骤二 由分数阶微分的线性性和 Laplace 变换的性质证明唯一性。假设有存在两

18、个满足上述初值问题的解1()y t、2()y t 令12()()()z ty ty t，有分数阶微分方程的线性性可得：0100()0(),1,2,nktttkDz tDy tbkn （16）从而有()0Z s （17）由 Laplace 变换的性质可知：()0z t 在(0,)T上几乎处处成立。故原方程的解在1(0,)LT上唯一。注：上述证明过程中用到的 Laplace变换法是一种常用的分数阶微分方程求解方法，该方法步骤简单，适用范围较广，在实际中有着重要应用，后面将对其进行详细介绍。第二步：运用第一步的结论证明原初值问题解的存在唯一性。定理 2 如果1()(0,)f tLT且()(1,2,)

19、jp tjn是0,T上的连续函数，则初值问题（9）（10）有唯一解1()(0,)y tLT。定理的证明过程如下：步骤一 化微分方程为积分方程 假设原方程有解()y t并记0()()ntDy tt，那么运用定理 1 可得：11101()()()()()intniinity ttt dtb （18）作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 将上式代入到原微分方程表达式（9）可得：0()(,)()()ttK tdtg t （19）其中 1

20、11()1()()(,)()()()nnknn tn kknnkttK tppt （20）111111()()()()()()iiknnnnin kiiki kiikttg tf tp tbptb （21）步骤二 证明变换后的积分方程有唯一解 用不动点定理易证结论成立。步骤三 说明原微分方程有唯一解 由定理 1 易得。2、一般形式的分数阶微分方程的存在唯一性定理 考虑如下形式的微分方程：0()(,)ntDy tf t y （22）100()kttkDy tb，1,2,kn （23）其中，(,)f t y的定义域为平面(,)t y上的一个子区域G，且存在G上的子区域(,)R h K满足：0th，

21、1111()()iniiitty tbK （24）定理 3 设(,)f t y为G上的连续实值函数，且在G上关于y满足Lipschitz条件，即 1212(,)(,)f t yf t yA yy （25）从而(,)f t yM，对任意(,)t yG 且 11(1)nnMhK 那么，方程（22）（23）在区域(,)R h k有唯一的连续解。定理的证明过程如下：步骤一 化微分方程为等价积分方程；对方程（22）按nD,1nD，1D逐次进行分部积分可得：作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因

22、人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 11101()()(,(,()()()intniiinby tttffyd （26）步骤二 证明上述等价积分方程解的存在性；构造函数序列0()y t，1()y t，2()y t，如下：101()()iniiiby tt （27）111101()()(,()()()intnimmiiibytttfyd （28）首先，我们可以证明对任意的0th 及任意的m有()(,)mytR h k。1111111111011()()(,()()()(1)(1)innntnmimiinnntttytbtfydMtMhK （29）进一步，我们可由数学归纳法证明

23、，对任意的m有下式成立：11()()(1)nmmmmnMAtytytm （30）证明过程如下：在式（29）中令1m 可得：10()()(1)nnMty ty t （31）假设当mk时，式（30）成立，即下式成立:11()()(1)nkkkknMAty tytk （32）那么，当1mk 时有：110100(1)()()()()()()1()(1)()(1)(1)(1)(1)(1(1)nnnnnnnntkkkkntkknnkktnkknnnnkknAytyttyydMAtdkMADtkMAktkkMA tk （33）作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微

24、积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 从而由归纳法可知，对任意的m，式（30）成立。进而，有1(1)nmmnMAtm 的收敛性可知，函数序列()myt收敛。令()lim()mmy tyt，易证()y t是等价积分方程（26）的解，也即是原微分方程（22）的解。步骤三 证明上述等价积分方程解的唯一性；假设()y t也是等价积分方程（26）的解，令()()()z ty ty t，有：101()()(,()()ntnz ttfzd （34）由()z t的连续性可知，存在常数B，使得对任意的0th，()z tB。利用

25、式（34）可得：(),(0)(1)nnABtz tth （35）将该估计过程重复j次可得：(),1,2,(1)njjnA Btz tj （36）又 0(1)limnjjjnA Bt 故()0,(0)z tth，也即()(),(0)y ty tth。注：有上面的介绍可知，整个线性分数阶微积分方程解的存在唯一性理论的证明过程都是建立在不动点理论的基础上的，使得我们必须将讨论范围限制在有限区间内的满足 Lipschitz条件的函数上，如何打破这个限制是一个值得思考的问题。在某些情况下，定理 3 可直接作为分数阶微分方程的求解方法，通常称之为存在唯一性解法。由上面的介绍，我们可将分数阶微积分方程存在唯

26、一性理论及其所面临的问题描述如下：作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 线性分数阶微积分的方程存在唯一性理论非线性分数阶微积分的方程存在唯一性理论分数阶微积分方程的存在唯一性理论理论基础：Laplace变换法、不动点原理理论基础：等价积分方程、不动点原理限制条件：（1）有限区间(2)即 绝对可积 1()(0,)f tLT()f t 0,T限制条件：（1）有限区间（2）条件 0,T Lipschitz？打破限制条件？修改现有理论基

27、础 3、分数阶微分方程初值问题解的依赖性 下面我们来考察初值条件的微小变化将对方程的解造成怎样的影响，为此，我们在初值条件中引入一个微小的改变量。100()kttkkDy tb，1,2,kn （37）其中(1,2,)kkn为任意常数。定理 4 设()y t是初值问题（22）（23）的解，()y t是初值问题（22）、（37）的解，那么对任意的(0,th有：1,1()()()inininiy ty ttEAt （38）其中,niE 为ML函数。证明：步骤一 用定理3的方式构造两组函数序列0()y t，1()y t，2()y t，和 0()y t，1()y t，2()y t，使得()lim()mm

28、y tyt，()lim()mmy tyt。步骤二 由数学归纳法容易证明 110()()()nikknmmmiikniA tytyttk （39）步骤三 对上式两端取极限m 可得：1101,1()()()()niinnikknmiikniniiA ty ty ttktEAt （40）作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 四、Laplace 变换求解法 随着分数阶微分方程在工程应用中出现得越来越频繁，给出分数阶微分方程的有效而简便的

29、求解方法便显得越来越重要，然而现有的求解方法都有着各种各样的缺陷。下面我们介绍一种基于 Laplace 变换的分数阶微分方程求解方法，该方法简单、直观，适用于常系数线性分数阶微分方程的求解。1、Laplace 变换求解法（1）Laplace 变换求解法的主要步骤 步骤一：对原微分方程做 Laplace 变换，化微分方程为代数方程；步骤二：求解该代数方程，得到原问题在变换域上的解；步骤三：对该变化域上的解做 Laplace 逆变换得到原问题的时域解。（2）Laplace 变换求解法的应用 下面我们通过两个例子来说明 Laplace 变换法的应用方法。例 1 我们考虑用 Laplace 变换法对如

30、下的非齐次标准分数阶微分方程的初值问题进行求解。0()()()tD y ty th t，(0)t （41）00()kttkDy tb，(1,2,)kn （42）其中，1nn 解：对方程（41）两端做 Laplace 变换，并利用初值条件（42）可得：11()()()nkkks Y sY sH sb s 从而 11()()knkkH ssY sbss （43）对式（31）做 Laplace 逆变换可得原微分方程的解为：1,1,10()()()()()tnkkkky tb tEttEthd （44）注：某些文献中也给出了该问题用迭代法进行求解的过程，虽然两种解法的结果相同，但显然 Laplace

31、求解法更为直观、简便。例 2 下面我们考虑用 Laplace 变换法对序列分数阶微分方程的初值问题进行求解。2100()()()ttDDy ty th t （45）2110001()tttDDy tb，11002()ttDy tb （46）解：对方程（45）两端做 Laplace 变换，并利用初值条件（46）可得：作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 12221()()()sY sH ss bb 从而 22121()()H ss

32、 bbY ss （47）对式（35）做 Laplace 逆变换可得原微分方程的解为：111112,1,0()()()()()()ty tb tEtbtEttEthd （48）其中，12 。注：对比上面两个初值问题容易看到他们在形式上非常相似，唯一的差别体现在一个是基于经典分数阶微积分定义的标准分数阶微分方程，一个是基于序列分数阶微积分定义的序列分数阶微分方程，从而在初值地给法不一样。但我们发现它们的解在表达式上也非常地相近，对比结果如下：1,1,10()()()()()tnkkkky tbtEttEthd 令 1,()()G ttEt ,11()()()()nkkkky tbtEtG th t

33、 令 1,()()G ttEt 111112,1,0()()()()()()ty tbtEtbtEttEthd 11112,1,()()()()()y tb tEtbtEtG th t 通过上面的对比可以发现，对应的标准分数阶微分方程和序列分数阶微分方程的解有一个共同点，即它们具有同样的Green函数，下面我们就Green展开讨论。2、Green 函数 考虑如下的初值问题：0()()tL y tf t，100()0kttDy t，(1,2,)kn （49）其中 10001()()()()()()nn knttktnkL y tDy tp tDy tp t y t（1）定义 若函数(,)G t满

34、足如下条件，则称其为方程（37）的Green函数：1）(,)0tLG t对任意的(0,)t；2）1,lim(,)ktk ntDG t，0,1,kn（,k n是 Kronecker delta函数）；3）,0lim(,)0ktttD G t，0,1,kn。作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载（2）性质 1）0()(,)()ty tG tfd 是方程（49）的解；2）对常系数分数阶微分方程有：(,)()G tG t；3）对(,)G

35、t的适当微分可得到一组齐次方程（()0f t）的线性无关解。下面我们利用Green函数的定义来证明上述三个性质。证明 1）：计算0()ktDy t如下：111110000010100()()(,)()(,)()lim(,)()(,)()lim(,)()(,)()kkkkkkkkkkkkktttttttttttttttttttDy tDDy tDDG tfdDDG tfdDDG tfDDG tfdDG tfD G tfdkn 0(,)()()kttD G tfdf tkn ,k n（定义2）（50）将上述等式所表示的10()tDy t，20()tDy t，0()ntDy t累加起来有:00()(

36、,)()()()tttL y tLG tfdf tf t 0（定义1）（51）证明 2）：由 Laplace 求解法可得。证明 3）：第一步，取0n ，令0()()tytD G t，由交换律的条件及定义 3）有：00000()()()0tttttL ytLD G tDLG t （52）故()yt是原齐次方程的解。第二步，由结合律可知：1100000100()()()1nnntttttttDytDD G tDG t （53）有上面的分析可知，对线性常系数分数阶微分方程的求解就转换为寻找该方程的分数阶Green函数，有该Green函数便可直接写出该方程的解，解的具体表达式如下：10()()()()

37、tnkkky tbtG tfd （54）作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 100()kkttbDy t 0()()nkkttDG t 下面我们给出寻找一般的线性常系数分数阶微分方程的Green函数的方法。（3）求法 对于常系数分数阶微分方程而言，我们可以通过对其做Laplace变换来得到相应的Green函数，下面我们通过一个例子来对其进行说明。例 3 考虑如下的一阶常系数微分方程：0()()ta D y tf t （55）解

38、：由于初值对Green函数无影响，因此我们对上式做Laplace变换并忽略初值部分可得：()()as Y sF s （56）从而 1()()Y sF sas （56）故 1()g sas （57）对其做 Laplace 逆变换可得：11()()tG saa （58）注：对如下的n阶分数阶微分方程：101110()()()()()nnnna Dy taDy ta Dy ta Dy tf t （59）假设110nn可得：01()knnkkgsa s （60）01202 021101211001200,()()120()1()()1(1)()(;,)!()()nnnnnnnjjjinnnnnnnjj

39、jmnnmkkkmnkkmknkiinmnmknG tm k kkamataaEta （61）作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和精品资料 欢迎下载 五 小结 有关分数阶微分方程的理论分析部分我们主要介绍了两方面的内容，一是分数阶微分方程解的性质，一是分数阶微分方程的求解方法。由于对分数阶微分方程的研究还不够成熟，因此对其所做出的理论分析还处于探索阶段。已有成果多半是对经典微积分方程理论的简单推广，且只能覆盖部分特殊形式的分数阶微分方程。现有的很多工作都是试图寻找新的理论方法，以打破现有的限制条件，力求构建一套完善的分数阶微分方程理论。作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受很的研究时经典整数阶微积分方程对这些系统的描述将遇到以下问题需要求解都非常繁琐基于以上原因人们迫切期待着有一种可用的数学工具和