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( )()( )( )t
7、nnafttuf 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:( )f t在 ,
9、 a b上具有1p阶连续导数;条件 3:()( )0kfa,1,2,kp;条件 4:( )0f t,ta。由上图我们可以看到,对于不同的分数阶微积分定义方式有着不同定义域,即便是在公共区域内, 不同的定义方式之间也无法实现完全的统一,这对分数阶微积分的应用和研究造成了一定的困难,因此人们迫切期望着分数阶微积分的一种哪怕是形式上的统一定义方式。2、M-R序列分数阶微分的定义为了满足实际需要, 下面我们试图从形式上对分数阶微积分给出一种统一的表达式。分数阶微积分的主要思想是推广经典的累次微积分,所有推广方法的共同目标是以非整数参数p 取代经典微积分符号中的整数参数n,即:nnddtppddt实际上
10、,任意的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)注:虽然上述序列分数阶微分的Laplace 变换是在
13、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 满足Lispschitz
14、条件, 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)步骤二 由分数阶微分的线性性和Laplace 变换的性
18、质证明唯一性。假设有存在两个满足上述初值问题的解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 变换法是一种常用的分数阶微分方程求解方法,该方法步骤简单,适用范围较广,在实际中有着重要应用,后面将对其进行详细介绍。第二步:运用第一步的结论证明原初值问题解的存在唯一性。定理 2 如果1( )(
19、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(19)其中111( )1()(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年分数阶微分方程-课件 2022 年分 微分方程 课件
限制150内