分数阶微分方程课件.docx

《分数阶微分方程课件.docx》由会员分享，可在线阅读，更多相关《分数阶微分方程课件.docx（34页珍藏版）》请在三一办公上搜索。

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

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

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

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

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

6、微积分定义依据1：整数阶微分的差分定义依据2：整数阶积分的柯西公式f(k)(t)=lim1kk(-1)rf(t-rh),kNh0hkr=0rf(-n)(t)=1t(t-u)n-1f(u)duG(n)a定义域：f(t)p+1阶连续导数f(t)在a,b具有先积分后微分定义域：定义域：f(t)f(t)在a,b上逐段连续且在任意有限子区间上可积f(t)p+1阶连续导数f(t)在a,b具有先微分后积分定义域： f(t)f(t)=0;taG-L分数阶微积分定义：GaDtpf(t)=lim1h0hpn=r=0t-ahR-L分数阶微积分定义：dR-(n-p)Rpf(t)aDtf(t)=aDt,n-1p0aDt

7、f(t)=G(n-p)a1tC-p(t-u)p-1f(u)du,p0aDtf(t)=G(p)a基于广义函数的分数阶微积分定义：ap(-1)rf(t-rh)r%pf(t)=f(t)*F(t)Dt-ptp-1,t0Fp(t)=G(p)0,t0缺点：定义域较窄计算复杂优点：定义域较宽缺点：Laplace变换较复杂优点：Laplace较简洁缺点：定义域太窄优点：有利于工程中对系统的描述从上图我们看到，在分数阶微积分的发展过程中，人们根据不同的需求，从不同角度给出了分数阶微积分的定义，但这些定义无论从对象上还是从表达式上都无法实现统一，它们之间的关系大致可以用下图来表示。 注： 条件1：f(t)在a,b

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

9、p取代经典微积分符号中的整数参数n，即： dn dpdtndtp实际上，任意的n阶微分都可以看成是一列一阶微分的叠加： dnf(t)ddd=Lf(t) dtndtdtdt14243n由此，我们可以给出一种在很多实际应用中十分重要的分数阶微积分的推广方式。首先，我们假设已有一种合适的推广方式来将一阶微分推广为a阶微分，即dDa是可实现的。那么类似地可得到的推广式为： dtnaaa Dnaf(t)=1 DDL43D(f ) t 42这种推广方式最初是由K.S.Miller和B.Ross提出来的，其中Da采用的是R-L分数阶微分定义，他们称之为序列分数阶微分。序列分数阶微分的其他形式可以通过将Da替

10、换为G-L分数阶微分、Caputo分数阶微分或其他任意形式分数阶微分来得到。 进一步，如果我们将中的分数阶微分Da替换为不同阶数的分数阶微分可得到序列分数阶微分更一般的表达式： Daf(t)=Da1Da2LDanf(t) a=a1+a2+L+an 根据问题的需要，Da可以是R-L分数阶微分、G-L分数阶微分、Caputo分数阶微分或其他任意形式的分数阶微分，从这一点看来，我们可以说序列分数阶微分从形式上给出了分数阶微积分在时域上的一个统一表达式，R-L分数阶微分、G-L分数阶微分和Caputo分数阶微分都只是序列分数阶微分的一种特殊情况。故而，下面我们在对分数阶微积分方程进行理论分析的时候可以

11、仅仅针对序列分数阶微积分来给出结论。 3、M-R序列分数阶微分的Laplace变换 下面我们考虑如下形式的序列分数阶微分的Laplace变换。 saaa aDtm=aDtmaDtm-1LaDt1 sm-1am-1am-1a1 D=DDLDatatat atm sm=aj=1j00，使得对所有的x,yX， d(Tx,Ty)qd(x,y) 则称T满足Lispschitz条件，q成为T的Lispschitz常数。 特别的，如果q1，则T称为压缩映射。 定理2 T:XX是压缩映射，设X,d是距离空间，则T在X中恰有一个不动点。设这个不动点为x，则对任何初始点x0X，逐次迭代点列xn+1=Txn，n=1

12、,2,L收敛于x，且关于收敛速度有如下估计式： d(xn,x)qn(1-q)-1d(Tx0,x0) 其中，q是T的Lipschitz常数。 三、 解的存在唯一性理论 近年来，分数阶微分方程已经在国内外引起极大的研究兴趣，尤其是关于其解的性质的研究，诸如存在性及唯一性等，其中大多数的研究方法是通过把分数阶初值问题转换成等价的分数阶积分方程,然后运用不动点定理来得到分数阶初值问题解的存在唯一性结果。已有研究结果主要有以下限制： 函数的定义区间为有限区间a,b； 函数在定义域上需满足Lipschitz条件； 因此，目前人们在这方面所做的工作都是希望设法在放宽上述两个限制条件后给出分数阶微积分方程的解

13、的存在唯一性定理。 下面我们对分数阶微分方程初值问题的现有理论结果作一个简单的介绍，相应的结论都是针对定义在有限区间0,T上的M-R序列分数阶微分形式，在满足Lipschitz条件下给出的，当然，由前面的介绍可知，这些结论也可直接推广到其他分数阶微分形式。 1、 线性分数阶微分方程解的存在唯一性定理 考虑如下形式的初值问题： n-jy(t)+pn(t)y(t)=f(t),(0tT)0Dty(t)+pj(t)0Dtsnn-1j=1s(9)(10)0Dtsk-1y(t)t=0=bk,k=1,2,L,nT且f(t)L1(0,T)，即 0f(t)dt 第一步：假设pk(t)0,(k=1,2,L,n)，

14、考虑由此得到的退化问题解的存在唯一性。 定理1 如果f(t)L1(0,T)，则方程 0Dtsny(t)=f(t) 有满足初值条件的唯一解y(t)L1(0,T)。 定理的证明过程如下： 步骤一 通过Laplace变换证明解的存在性； 下面我们设法构造一个待求解问题解，对式做Laplace变换可得： sY(s)-ssn-sn-k0Dtsn-k-1y(t)t=0=F(s) snk=0n-1其中，Y(s)、F(s)分别是y(t)、f(t)的Laplace变换。利用初值条件可得： Y(s)=s-sn F(s)+bn-ks-sn-k k=0n-1对上式做Laplace逆变换可得： n1bisi-1sn-1

15、 y(t)=(t-t)f(t)dt+t G(sn)G(s)i=1i0t步骤二 由分数阶微分的线性性和Laplace变换的性质证明唯一性。 假设有存在两个满足上述初值问题的解y1(t)、y2(t) 令z(t)=y1(t)-y2(t) ，有分数阶微分方程的线性性可得： s0Dtnz(t)=0sk-1 0Dty(t)t=0=bk,k=1,2,L,n从而有 Z(s)=0由Laplace变换的性质可知：z(t)=0在(0,T)上几乎处处成立。 故原方程的解在L1(0,T)上唯一。 注：上述证明过程中用到的Laplace变换法是一种常用的分数阶微分方程求解方法，该方法步骤简单，适用范围较广，在实际中有着重

16、要应用，后面将对其进行详细介绍。 第二步：运用第一步的结论证明原初值问题解的存在唯一性。 定理2 如果f(t)L1(0,T)且pj(t)(j=1,2,L,n)是0,T上的连续函数，则初值问题有唯一解y(t)L1(0,T)。 定理的证明过程如下： 步骤一 化微分方程为积分方程 假设原方程有解y(t)并记0Dtsny(t)=j(t)，那么运用定理1可得： tsi-1n1ty(t)=(t-t)sn-1j(t)dt+bi G(sn)0G(s)i=1i将上式代入到原微分方程表达式可得： j(t)+K(t,t)j(t)dt=g(t) t0 其中 ntsi-1n-1tsi-sk-1g(t)=f(t)-pn(

17、t)bi-pn-k(t)bi G(s)G(s-s)i=1k=1i=k+1iikn(t-t)sn-1n-1(t-t)sn-sk-1K(t,t)=pn(t)+pn-k(t) G(sn)G(s-s)k=1nk步骤二 证明变换后的积分方程有唯一解 用不动点定理易证结论成立。 步骤三 说明原微分方程有唯一解 由定理1易得。 2、 一般形式的分数阶微分方程的存在唯一性定理 考虑如下形式的微分方程： s 0Dtny(t)=f(t,y) sk-1Dy(t)t=0=bk， k=1,2,L,n 0t 其中，f(t,y)的定义域为平面(t,y)上的一个子区域G，且存在G上的子区域R(h,K)满足： 0th，t1-s

18、1tsi-s1y(t)-biK G(si)i=1n定理3 设f(t,y)为G上的连续实值函数，且在G上关于y满足Lipschitz条件，即 f(t,y1)-f(t,y2)Ay1-y2 从而 Mhsn-s1+1 f(t,y)M，对任意(t,y)G 且 KG(1+sn)那么，方程在区域R(h,k)有唯一的连续解。 定理的证明过程如下： 步骤一 化微分方程为等价积分方程； 对方程按Dan,Dan-1，L，Da1逐次进行分部积分可得： b1 y(t)=itsi-1+(t-t)sn-1f(t,f(t,y(t)dt G(sn)0i=1G(si)步骤二 证明上述等价积分方程解的存在性； nt构造函数序列y0

19、(t)，y1(t)，y2(t)，L 如下： y0(t)=i=1nbisi-1t G(si)b1sn-1 ym(t)=itsi-1+ (t-t)f(t,ym-1(t)dt nti=1G(si)G(si)0首先，我们可以证明对任意的0th及任意的m有ym(t)R(h,k)。 nttytsi-s1t1-s1-s11sm(t)-bi=(t-t)n-1f(ti-1G(si)G(sn),ym-1(t)dt0sMtn-s1+1sn-s1+1G(1+s)MhG(1+sKnn)进一步，我们可由数学归纳法证明，对任意的m有下式成立： yMAm-1tmsnm(t)-ym-1(t)G(1+ms n)证明过程如下： 在

20、式中令m=1可得： yyMtsn1(t)-0(t)G(1+s n)假设当m=k时，式成立，即下式成立: yMAk-1tksnk(t)-yk-1(t)G(1+ks n)那么，当m=k+1时有： Ayk+1(t)-yk(t)(t-t)snyk(t)-yk-1(t)dtG(sn)0MAk1(t-t)sn-1tksndtG(1+ksn)G(sn)0MAk-snksnt0DtG(1+ksn)MAkG(1+ksn)tksn+snG(1+ksn)G(1+ksn+sn)MAkt(k+1)snG(1+(k+1)sn)tt从而由归纳法可知，对任意的m，式成立。 MAm-1tmsn进而，有的收敛性可知，函数序列ym

21、(t)收敛。 G(1+msn)令y(t)=limym(t)，易证y(t)是等价积分方程的解，也即是原微分方程m的解。 步骤三 证明上述等价积分方程解的唯一性； %(t)也是等价积分方程的解，令z(t)=y(t)-y%(t)，有： 假设y1sn-1 z(t)=(t-t)f(t,z(t)dt G(sn)0由z(t)的连续性可知，存在常数B，使得对任意的0th，z(t)0) Dtay(t)-ly(t)=h(t)， 0Dta-ky(t)t=0=bk，(k=1,2,L,n) 其中， n-1an 解：对方程两端做Laplace变换，并利用初值条件可得： sY(s)-lY(s)=H(s)+bksk-1 ak

22、=1n从而 nH(s)sk-1Y(s)=a+bkas-lk=1s-l对式做Laplace逆变换可得原微分方程的解为： y(t)=bktk=1nta-k Ea,a-k+1(lt)+(t-t)a-1Ea,a(l(t-t)a)h(t)dt 0a注：某些文献中也给出了该问题用迭代法进行求解的过程，虽然两种解法的结果相同，但显然Laplace求解法更为直观、简便。 例2 下面我们考虑用Laplace变换法对序列分数阶微分方程的初值问题进行求解。 0 Dta2(0Dta1y(t)-ly(t)=h(t) 0Dta2-1(0Dta1y(t)t=0=b1，0Dta1-1y(t)t=0=b2 解：对方程两端做La

23、place变换，并利用初值条件可得： (sa1+a2-l)Y(s)=H(s)+sa2b2-b1 从而 H(s)+sa2b2+b1Y(s)= sa2+a1-l对式做Laplace逆变换可得原微分方程的解为： ty(t)=b2ta1-1 Ea,a1(lta)+b1ta-1Ea,a(lta)+(t-t)a-1Ea,a(l(t-t)a)h(t)dt 0其中，a=a1+a2。 注：对比上面两个初值问题容易看到他们在形式上非常相似，唯一的差别体现在一个是基于经典分数阶微积分定义的标准分数阶微分方程，一个是基于序列分数阶微积分定义的序列分数阶微分方程，从而在初值地给法不一样。但我们发现它们的解在表达式上也非

24、常地相近，对比结果如下： y(t)=bkta-kEa,a-k+1(lta)+(t-t)a-1Ea,a(l(t-t)a)h(t)dtk=10ntty(t)=b2tEa,a1(lt)+b1tEa,a(lt)+(t-t)a-1Ea,a(l(t-t)a)h(t)dta1-1aa-1a0令 G(t)=ta-1Ea,a(lta) a1-1 y(t)=bE,al(at)+2ta1令 G(t)=at-1aE,a(lat)1y(t)=bkta-kEa,a-k+1(lta)+G(t)*h(t)k=1na-1bta,aEl(a+t)(G*)t(h)t通过上面的对比可以发现，对应的标准分数阶微分方程和序列分数阶微分方

25、程的解有一个共同点，即它们具有同样的Green函数，下面我们就Green展开讨论。 2、 Green函数 考虑如下的初值问题： 0 (k=1,2,L,n) Lty(t)=f(t)，0Dtsk-1y(t)t=0=0，其中 Lty(t)0Dty(t)+pk(t)0Dtsn-ky(t)+pn(t)y(t) snk=1n-10 定义 若函数G(t,t)满足如下条件，则称其为方程的Green函数： 1）tLGt(t,t)=0对任意的t(0,t)； sk-1DG(t,t)=dk,n，k=0,1,L,n2）lim(； tt-tt3）lim+(tDtskG(t,t)=0，k=0,1,L,n。 t,t0tt 性

26、质 1）y(t)=G(t,t)f(t)dt是方程的解； 0t2）对常系数分数阶微分方程有：G(t,t)G(t-t)； 3）对G(t,t)的适当微分可得到一组齐次方程的线性无关解。 下面我们利用Green函数的定义来证明上述三个性质。 证明1）：计算0Dtsky(t)如下： 0Dty(t)=0Dt(0Dtt0tskaksk-1y(t)=0Dtakt0tDtk-1G(t,t)f(t)dtttssk-1kk-1=tDaG(t,t)f(t)dt+lim-tDtak-1(tDsG(t,t)f(t)t(tDtt=0tDtk(tDtask-1G(t,t)f(t)dt+lim-tDtttsk-1G(t,t)f

27、(t) dk,）n tDskG(t,t)f(t)dtt0t=ttDtskG(t,t)f(t)dt+f(t)0knk=n将上述等式所表示的0Dts1y(t)，0Dts2y(t)，L，0Dtsny(t)累加起来有: 0Lty(t)=tLGt(t,t)f(t)dt+f(t)=f(t)0t0）证明2）：由Laplace求解法可得。 证明3）：第一步，取0lbn-1Lb1b0可得： gn(s)=1nasbkk k=056） 57）60）E(m)(bn-bn-1)m+bn+n-2j=0(bn-1-bj)kj(-an-1bn-bn-1t)an五 小结 有关分数阶微分方程的理论分析部分我们主要介绍了两方面的内容，一是分数阶微分方程解的性质，一是分数阶微分方程的求解方法。由于对分数阶微分方程的研究还不够成熟，因此对其所做出的理论分析还处于探索阶段。已有成果多半是对经典微积分方程理论的简单推广，且只能覆盖部分特殊形式的分数阶微分方程。现有的很多工作都是试图寻找新的理论方法，以打破现有的限制条件，力求构建一套完善的分数阶微分方程理论。