2022年分数阶微分方程课件 .pdf

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

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

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

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

4、非常特殊的方程的求解，且常用的求解方法都是具有局限性的。在数值求解方面，现有分数阶方程数值算法还很不成熟，主要表现为：（1）在数值计算中一些挑战性难题仍未得到彻底解决，如长时间历程的计算和大空间域的计算等；（2）成熟的数值算法比较少，现在研究较多的算法主要集中在有限差分方法与有限单元法；（3）未形成成熟的数值计算软件，严重滞后于应用的需要。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 15 页优秀学习资料欢迎下载鉴于此，发展新数值算法， 特别是在保证计算可靠性和精度的前提下，提高计算效率，解决分数阶微分方程计算量和存储量过大的难点问题

5、，发展相应的计算力学应用软件成为迫切需要关注的课题。二、 预备知识1、 分数阶微积分经典定义回顾作为分数阶微积分方程的基础， 本书在第二章中对分数阶微积分的定义及性质做了系统的介绍，为了接下来讨论的需要， 我们首先对其进行一个简要的回顾。（1）分数阶微积分的主要思想如上图所示， 分数阶微积分的主要思想是推广经典的整数阶微积分，从而将微积分的概念延拓到整个实数轴，甚至是整个复平面。 但由于延拓的方法多种多样，因而根据不同的需求人们给出了分数阶微积分的不同定义方式。然而这些定义方式不仅只能针对某些特定条件下的函数给出，而且只能满足人们的某些特定需求，迄今为止，人们仍然没能给出分数阶微积分的一个统一

6、的定义, 这对分数阶微积分的研究与应用造成了一定的困难。（2）几种经典的分数阶微积分定义下面我们试图从理论依据、 定义域、表达式和优缺点几个方面给出常见的四种分数阶微积分定义的比较图。依据1：整数阶微分的差分定义( )001( )lim( 1)(),kkrkhrkftf trhkNrhG-L分数阶微积分定义：001( )lim( 1)()tanhGpratphrpDf tf trhrh1( )1 ( )()( ),0()tCpn pnataDf ttufudu pnp11 ( )()( ),0( )tCppataDf ttuf u du pp依据2：整数阶积分的柯西公式()11( )()( )

7、( )tnnafttuf u dunR-L分数阶微积分定义：()11 ( )()( )()tRn pn pataDf ttuf udunp() ( ) ( ),1,nRpRnpatatdDf tDf tnpnnNdtCaputo 分数阶微积分定义：先积分后微分先微分后积分缺点：定义域较窄计算复杂优点：定义域较宽缺点：Laplace 变换较复杂优点：Laplace 较简洁缺点：定义域太窄分数阶微积分定义定义域：( )( )f tf ta,b在上逐段连续且在任意有限子区间上可积定义域：( )( )f tf ta,bp+1在具有阶连续导数定义域：( )( )f tf ta,bp+1在具有阶连续导数基

8、于广义函数的分数阶微积分定义：( )( )()patpD f tf tt1,0( )( )0,0pptttpt定义域：( )( )0;f tf tta优点：有利于工程中对系统的描述从上图我们看到， 在分数阶微积分的发展过程中，人们根据不同的需求， 从不同角度给出了分数阶微积分的定义，但这些定义无论从对象上还是从表达式上都无法实现统一，它们之间的关系大致可以用下图来表示。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 15 页优秀学习资料欢迎下载注：条件 1：( )f t在 , a b上逐段连续，且在任何有限子区间上可积；条件 2：( )

9、f t在 , a b上具有1p阶连续导数；条件 3：()( )0kfa，1,2,kp；条件 4：( )0f t，ta。由上图我们可以看到，对于不同的分数阶微积分定义方式有着不同定义域，即便是在公共区域内， 不同的定义方式之间也无法实现完全的统一，这对分数阶微积分的应用和研究造成了一定的困难，因此人们迫切期望着分数阶微积分的一种哪怕是形式上的统一定义方式。2、M-R序列分数阶微分的定义为了满足实际需要， 下面我们试图从形式上对分数阶微积分给出一种统一的表达式。分数阶微积分的主要思想是推广经典的累次微积分，所有推广方法的共同目标是以非整数参数p 取代经典微积分符号中的整数参数n，即：nnddtpp

10、ddt实际上，任意的n阶微分都可以看成是一列一阶微分的叠加：( )( )nnnd f td ddf tdtdt dtdt（1）由此，我们可以给出一种在很多实际应用中十分重要的分数阶微积分的推广方式。首先，我们假设已有一种合适的推广方式来将一阶微分推广为（01）阶微分，即dDdt是可实现的。那么类似地可得到（1）的推广式为：( )( )nnDf tD DDf t（2）这种推广方式最初是由. .K S Miller 和.B Ross提出来的，其中D 采用的是精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 15 页优秀学习资料欢迎下载RL 分

11、数阶微分定义，他们称之为序列分数阶微分。序列分数阶微分的其他形式可以通过将 D 替换为 GL 分数阶微分、Caputo分数阶微分或其他任意形式分数阶微分来得到。进一步，如果我们将（ 2）中的分数阶微分 D 替换为不同阶数的分数阶微分可得到序列分数阶微分更一般的表达式：12( )( )nDf tDDDf 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( ); ( )( )tttLDf tss F sDf t（7）重复利用上式m次可得：110000( ); ( )( )mmmmkmkmtttkLDf tssF ssDf t（8）注：虽然上述序列分数阶微分的Lapla

13、ce 变换是在 R-L 分数阶微分定义下进行证明的，但是该结论对其他几种分数阶微积分也是成立的。4、泛函理论基础定理 1（Schauder 不动点定理）设U 是 Banach空间 X 的有界闭子集，如果:T UU 是连续映射，那么 T 在U 中存在不动点，即使得Txx 的点存在。定义 1（Lipschitz条件）设,X d 是距离空间， T 是从 X 到 X 的映射，如果存在常数0q，使得对所有的,x yX，精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 15 页优秀学习资料欢迎下载(,)( ,)d Tx Tyqd x y则称 T 满足

14、Lispschitz条件， q成为 T 的Lispschitz常数。特别的，如果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、 线性分数阶微分

16、方程解的存在唯一性定理考虑如下形式的初值问题：1001100( )( )( )( )( )( ) , (0)(9)( ),1,2,(10)njnkntjtnjttkDy tptDy tpt y tf ttTDy tbkn且1( )(0,)f tLT，即0( )Tf t dt（11）第一步：假设( )0 , (1,2, )kptkn，考虑由此得到的退化问题解的存在唯一性。定理 1 如果1( )(0,)f tLT，则方程0( )( )ntDy tf t（12）有满足初值条件（ 10）的唯一解1( )(0,)y tLT。定理的证明过程如下：步骤一 通过 Laplace 变换证明解的存在性；精选学习资

17、料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 15 页优秀学习资料欢迎下载下面我们设法构造一个待求解问题解，对式（12）做 Laplace 变换可得：11000( )( )( )nnnknknttks Y ssDy tF s（13）其中，( )Y s、( )F s分别是( )y t、( )f t的 Laplace 变换。利用初值条件（ 10）可得：10( )( )nnknn kkY ssF sbs（14）对上式做 Laplace 逆变换可得：11101( )()( )()()nitniiniby ttfdt（15）步骤二 由分数阶微分的线性性和

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

19、定理 2 如果1( )(0,)f tLT且( )(1,2, )jp tjn是0,T上的连续函数， 则初值问题（9）（ 10）有唯一解1( )(0,)y tLT。定理的证明过程如下：步骤一 化微分方程为积分方程假设原方程有解( )y t并记0( )( )ntDy tt，那么运用定理 1 可得：11101( )()( )()()intniinity ttt dtb（18）精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 15 页优秀学习资料欢迎下载将上式代入到原微分方程表达式（9）可得：0( )( , ) ( )( )ttK tdtg t（1

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

21、( ,)R h K满足：0th，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逐次进行分部积分可得：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7

22、页，共 15 页优秀学习资料欢迎下载11101( )()( ,( ,( )()()intniiinby tttffyd（26）步骤二 证明上述等价积分方程解的存在性；构造函数序列0( )yt，1( )y t，2( )yt，如下：101( )()iniiibytt（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)nkkkknMAtytytk（32）那么，当1mk时有：110100(1)( )( )()()( )()1()(1)()(1)(1)(1)(1)(1(1)nnnnnnnntkkkkntkknnkktnkknnnnkknAytyttyydMAtdkMADtkMAktkkMA tk（33）精选学习资料 - - - - - - - - - 名师归纳总结 -

24、 - - - - - -第 8 页，共 15 页优秀学习资料欢迎下载从而由归纳法可知，对任意的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 。利用式（ 34

25、）可得：( ),(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、论及其所面临的问题描述如下：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 15 页优秀学习资料欢迎下载线性分数阶微积分的方程存在唯一性理论非线性分数阶微积分的方程存在唯一性理论分数阶微积分方程的存在唯一性理论理论基础：Laplace变换法、不动点原理理论基础：等价积分方程、不动点原理限制条件：（1）有限区间(2) 即绝对可积1( )(0,)f tLT( )f t0,T限制条件：（1）有限区间（2）条件0, TLipschitz？打破限制条件？修改现有理论基础3、 分数阶微分方程初值问题解的依赖性下面我们来考察初值条件的微小变化将对方程

27、的解造成怎样的影响，为此，我们在初值条件中引入一个微小的改变量。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( )yt，和0( )y t，1( )y t，2( )yt，使得( )lim( )mmy tyt ，( )lim( )mmy tyt

28、。步骤二 由数学归纳法容易证明110( )( )()nikknmmmiikniA tytyttk（39）步骤三 对上式两端取极限 m可得：1101,1( )( )()()niinnikknmiikniniiA ty ty ttktEAt（40）精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 15 页优秀学习资料欢迎下载四、Laplace 变换求解法随着分数阶微分方程在工程应用中出现得越来越频繁，给出分数阶微分方程的有效而简便的求解方法便显得越来越重要，然而现有的求解方法都有着各种各样的缺陷。下面我们介绍一种基于Laplace 变换的分

29、数阶微分方程求解方法，该方法简单、直观，适用于常系数线性分数阶微分方程的求解。1、 Laplace 变换求解法（1） Laplace 变换求解法的主要步骤步骤一：对原微分方程做Laplace 变换，化微分方程为代数方程；步骤二：求解该代数方程，得到原问题在变换域上的解；步骤三：对该变化域上的解做Laplace 逆变换得到原问题的时域解。（2） Laplace 变换求解法的应用下面我们通过两个例子来说明Laplace 变换法的应用方法。例 1 我们考虑用 Laplace 变换法对如下的非齐次标准分数阶微分方程的初值问题进行求解。0( )( )( )tD y ty th t，(0)t（41）00(

30、 )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 求解法更为直观、简便。例 2 下面我们考虑用Laplace 变换法对序列分数阶微分方程的初值

31、问题进行求解。2100( )( )( )ttDDy ty th t（45）2110001( )tttDDy tb，11002( )ttDy tb（46）解：对方程（ 45）两端做 Laplace 变换，并利用初值条件（ 46）可得：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 15 页优秀学习资料欢迎下载12221()( )( )sY sH ss bb从而22121( )( )H ss bbY ss（47）对式（ 35）做 Laplace 逆变换可得原微分方程的解为：111112,1,0( )()()()( () ) ( )ty t

32、b tEtbtEttEthd（48）其中，12。注：对比上面两个初值问题容易看到他们在形式上非常相似，唯一的差别体现在一个是基于经典分数阶微积分定义的标准分数阶微分方程，一个是基于序列分数阶微积分定义的序列分数阶微分方程，从而在初值地给法不一样。但我们发现它们的解在表达式上也非常地相近，对比结果如下：1,1,10()()()( () ) ( )tnkkkkytbtEttEthd令1,( )()G ttEt,11( )()( )()nkkkkytbtEtG th t令1,( )()G ttEt111112,1,0( )()()()( () ) ( )ty tbtEtbtEttEthd11112,

33、1,( )()()( )( )y tbtEtbtEtG th t通过上面的对比可以发现，对应的标准分数阶微分方程和序列分数阶微分方程的解有一个共同点，即它们具有同样的Green函数，下面我们就 Green展开讨论。2、 Green函数考虑如下的初值问题：0( )( )tL y tf t，100( )0kttDy t，(1,2, )kn（49）其中10001( )( )( )( )( ) ( )nnknttktnkL y tDy tptDy tp t y t（1） 定义若函数( , )G t满足如下条件，则称其为方程（37）的 Green函数：1）( , )0tL G t对任意的(0, ) t；

34、2）1,lim( , )ktk ntDG t，0,1,kn（,k n是 Kronecker delta函数） ；3）,0lim ( , )0ktttDG t，0,1,kn。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页，共 15 页优秀学习资料欢迎下载（2） 性质1）0( )( , )( )ty tG tfd是方程（ 49）的解；2）对常系数分数阶微分方程有：( , )()G tG t；3）对( , )G t的适当微分可得到一组齐次方程（( )0f t）的线性无关解。下面我们利用 Green函数的定义来证明上述三个性质。证明 1） ：计算

35、0( )ktDy t如下：111110000010100( )( )( , )( )( , )( )lim( , )( )( , )( )lim( , )( )( , )( )kkkkkkkkkkkkktttttttttttttttttttDy tDDy tDDG tfdDDG tfdDDG tfDDG tfdDG tfDG tfdkn0( , )( )( )kttDG tfdf tkn,k n（ 定义2）（50）将上述等式所表示的10( )tDy t，20( )tDy t，0( )ntDy t累加起来有 : 00( )( , ) ( )( )( )tttL y tLG tfdf tf t0（

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

37、d（54）精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页，共 15 页优秀学习资料欢迎下载100( )kkttbDy t0( )( )nkkttDG t下面我们给出寻找一般的线性常系数分数阶微分方程的Green函数的方法。（3） 求法对于常系数分数阶微分方程而言，我们可以通过对其做Laplace变换来得到相应的 Green函数，下面我们通过一个例子来对其进行说明。例 3 考虑如下的一阶常系数微分方程：0( )( )ta Dy tf t（55）解：由于初值对 Green函数无影响，因此我们对上式做Laplace变换并忽略初值部分可得：( )

38、( )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）01202021101211001200,()()120()1()()1( 1)( )(;,)!()()nnnnnnnjjjinnnnnnnjjjmnnmkkkmnkkmknkiinmnmknGtm kkkamataaEta（

39、61）精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 14 页，共 15 页优秀学习资料欢迎下载五 小结有关分数阶微分方程的理论分析部分我们主要介绍了两方面的内容，一是分数阶微分方程解的性质， 一是分数阶微分方程的求解方法。由于对分数阶微分方程的研究还不够成熟， 因此对其所做出的理论分析还处于探索阶段。已有成果多半是对经典微积分方程理论的简单推广，且只能覆盖部分特殊形式的分数阶微分方程。现有的很多工作都是试图寻找新的理论方法，以打破现有的限制条件，力求构建一套完善的分数阶微分方程理论。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页，共 15 页