分数阶微分方程-课件
- 格式:doc
- 大小:1.41 MB
- 文档页数:15
分数阶微分方程解法1、分数阶微积分介绍分数阶微积分是传统微积分的推广和扩展,在这门学科中,函数的导数和积分的阶数可以为分数,有时也可以是负数或复数。
与传统微积分相比,分数阶微积分的应用更加广泛,可以通过它来解释和研究各种复杂的自然现象,例如金融市场的非平稳性、地震的时序性等。
2、分数阶微分方程简介分数阶微分方程是指微分方程的微分阶数为分数,例如阶为1.5或2.7的微分方程。
在实际应用中,分数阶微分方程被广泛地用于描述自然现象的动态行为,例如分形、非线性动力学、力学、电动力学和生物学等。
3、分数阶微分方程的解法分数阶微分方程的解法是与传统微分方程不同的。
下面介绍两种常用的解法。
3.1、分式变换法分式变换法是最常用的解分数阶微分方程的方法之一。
它的基本思想是将分数阶微分方程转化为一些常见的函数或微分方程。
例如,我们考虑一个分数阶微分方程:D^βy(t)=f(t),其中D^β表示分数阶导数运算符,y(t)是未知函数,f(t)是已知函数。
现在我们把分数阶微分方程改写为下面的形式:y(t)=1/Γ(α)(d/dt)^α∫_0^t f(u)(t-u)^{α-1}du其中α和β之间的关系可以用以下公式表示:α=β-n这里,n是一个正整数,它满足0<n<=β。
通过这个公式,分数阶微分方程就被转化为常数阶微分方程和分式变换的形式。
3.2、拉普拉斯变换法拉普拉斯变换法也是解分数阶微分方程的有效方法。
它的基本思想是将分数阶微分方程转化为常数阶微分方程,然后通过拉普拉斯变换及其逆变换来得到方程的解。
例如,我们考虑一个分数阶微分方程:D^βy(t)+ay(t)=f(t),其中a是常数,D^β表示分数阶导数运算符,y(t)是未知函数,f(t)是已知函数。
现在我们把分数阶微分方程改写为下面的形式:L{D^βy(t)}+aL{y(t)}=L{f(t)}其中L表示拉普拉斯变换,而L{D^βy(t)}和L{f(t)}分别是D^βy(t)和f(t)的拉普拉斯变换。
分数阶微分方程第三讲分数阶微分方程基本理论一、分数阶微分方程的出现背景及研究现状1、出现背景分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。
整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有较完善的理论。
但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到以下问题:(1)需要构造非线性方程,并引入一些人为的经验参数和与实际不符的假设条件;(2)因材料或外界条件的微小改变就需要构造新的模型;(3)这些非线性模型无论是理论求解还是数值求解都非常繁琐。
基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。
分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。
2、研究现状在近三个世纪里,对分数阶微积分理论的研究主要在数学的纯理论领域里进行,似乎它只对数学家们有用。
然而在近几十年来,分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。
分数阶微积分理论也受到越来越多的国内外学者的广泛关注,特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。
随着分数阶微分方程在越来越多的科学领域里出现,无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。
然而由于分数阶微分是拟微分算子,它的保记忆性(非局部性)对现实问题进行了优美刻画的同时,也给我们的分析和计算造成很大困难。
在理论研究方面,几乎所有结果全都假定了满足李氏条件,而且证明方法也和经典微积分方程一样,换句话说,这些工作基本上可以说只是经典微积分方程理论的一个延拓。
capulo 分数阶微分
Caputo分数阶微分是数学中的一个概念,它用于描述一个函数在某个点或某个范围内的阶数导数。
与传统的整数阶导数不同,分数阶导数可以更好地描述一些复杂的物理或工程问题。
Caputo分数阶微分是在Caputo分数阶积分的基础上定义的。
它的一般形式为:
D^α_a^x f(x) = 1 / Γ(n - α) * ∫_a^x (x - t) ^ {n - α - 1} * f(t) * dt
其中,D^α是Caputo分数阶导数,f(x)是要分析的函数,n是使得f(t)具有意义的整数,a和x是积分上下限。
Caputo分数阶微分在描述一些具有记忆性和遗传性的问题时非常有用,例如在描述材料力学、控制理论、信号处理等领域的问题时。
在这些情况下,整数阶导数无法准确描述系统的动态特性,而分数阶导数可以更好地描述这些特性。
使用Caputo分数阶微分可以更好地模拟系统的非线性行为和记忆效应。
例如,在描述一个材料的蠕变行为时,Caputo分数阶微分可以更好地描述材料在不同应力水平下的行为。
然而,Caputo分数阶微分也有其局限性。
由于其定义涉及到积分,因此计算上可能会比整数阶导数更加复杂和困难。
此外,对于一些特殊的问题,可能需要特殊的方法来处理分数阶导数。
综上所述,Caputo分数阶微分是一种有用的数学工具,可以用于描述一些具有记忆性和遗传性的问题。
虽然计算上可能比整数阶导数更加复杂,但在一些特定的问题中,它可以提供更好的描述和预测能力。
分数阶微分方程第三讲分数阶微分方程基本理论一、分数阶微分方程的出现背景及研究现状1、出现背景分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。
整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有较完善的理论。
但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到以下问题:(1)需要构造非线性方程,并引入一些人为的经验参数和与实际不符的假设条件;(2)因材料或外界条件的微小改变就需要构造新的模型;(3)这些非线性模型无论是理论求解还是数值求解都非常繁琐。
基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。
分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。
2、研究现状在近三个世纪里,对分数阶微积分理论的研究主要在数学的纯理论领域里进行,似乎它只对数学家们有用。
然而在近几十年来,分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。
分数阶微积分理论也受到越来越多的国内外学者的广泛关注,特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。
随着分数阶微分方程在越来越多的科学领域里出现,无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。
然而由于分数阶微分是拟微分算子,它的保记忆性(非局部性)对现实问题进行了优美刻画的同时,也给我们的分析和计算造成很大困难。
在理论研究方面,几乎所有结果全都假定了满足李氏条件,而且证明方法也和经典微积分方程一样,换句话说,这些工作基本上可以说只是经典微积分方程理论的一个延拓。
(四)分数阶微积分我们重点考察R −L 型分数阶微积分的性质,简记RL 0D βt =D βt ,若⽆特殊说明。
a). 线性性D βt [f (t )+g (t )]=D βt f (t )+D βt g (t )D βt λf (t )=λD βt f (t )证明:直接带⼊定义验算即可.设m =[β]+1RL 0D βtf (t )=1Γ(m −β)d m dt m ∫t0(t −τ)m −β−1f (τ)d τb). 积分的叠加性D −αt D −βt f (t )=D −α−βt f (t ) (α,β>0)证明:对整数阶积分结论是显然的,对于分数阶R-L 积分仍然具有叠加性。
由定义知0D −βt f (t )=1Γ(β)∫t0(t −x )β−1f (x ):=g (t )那么D −αt g (t )=1Γ(α)∫t0(t −τ)α−1g (τ)d τ=1Γ(α)Γ(β)∫t 0(t −τ)α−1d τ∫τ0(τ−x )β−1f (x )dx=1Γ(α)Γ(β)∫t 0f (x )dx ∫tx (t −τ)α−1(τ−x )β−1d τ(交换积分次序)=1Γ(α)Γ(β)∫t 0f (x )dx ∫10(t −x )α+β−1(1−ξ)α−1ξβ−1d ξ (Let ξ=τ−xt −x )=B (α,β)Γ(α)Γ(β)∫t0(t −x )α+β−1f (x )dx=1Γ(α+β)∫t0(t −x )α+β−1f (x )dx=0D −α−βt f (t )由此我们也得到了积分满⾜交换性,即D −αt D −βt f (t )=D −α−βt f (t )=D −βt D −αt f (t ) (α,β>0)c). 上式考虑了积分叠加的情形,对于连续函数f (t )考虑混合运算“先积分再微分”.(还记得R-L 定义思路D β=D m D −(m −β))0D αt 0D −βt f (t )=0D α−βt f (t ) (α>0,β>0)证明:先探讨⼀种特殊的情形D λD −λf (t )=f (t ) (λ>0)当λ为整数时结论显然成⽴。
含caputo分数阶导数的分数阶微分方程分数阶微分方程是指其解中包含阶级高于一的分数和跨越多个时间段的定积分形式。
该方程研究混合不同时段内的过程变化,特别是近几年在复杂科学问题研究方面得到了广泛的应用,比如金融理论等。
Caputo分数阶微分方程是一种特殊的分数阶微分方程,它以概念全面解决了分数阶微分方程的研究。
Caputo分数阶微分方程的研究始于20世纪七八十年末意大利,由M. Caputo首先提出。
在Caputo的原始文献中,他推导了一类分数阶微分方程,主要研究过程当中时间和空间变化的关系。
Caputo分数阶微分方程有如下定义:假设函数f(t)是在区间(a,b)上可微,则在该区间上满足条件$${D_a^nf(t)=\int_a^b{K(t,s)f^{(n)}(s)}ds}$$的函数称为阶级为n的Caputo分数阶微分方程的解,K(t,s)即为积分算子。
其中n为分数,当n<1,Caputo分数阶微分方程就是一类特殊的非定常方程,可以用于模拟效应非定常变化的过程,变化函数f(t)的值取决于导数的序号。
Caputo分数阶微分方程主要用于研究各种变分问题,重点关注函数解决定积分参数变化中无视时间关系影响的问题。
Caputo分数阶微分方程应用于工程领域也是一个激动人心的问题。
比如利用Caputo分数阶微分方程求解微分方程模型,可以有效模拟工程实际的过程。
在工程和物理领域,Caputo分数阶微分方程应用最广泛,尤其是其算子在时间参数函数f(t)中表现出来的使用,及其形成模型包含秩序和变化关系。
当具体系统存在改变或不同状态时,Caputo分数阶微分方程可以用于描述和模拟这一系列模型。
从此可见,Caputo分数阶微分方程在研究复杂科学问题中有着广泛的应用,它求解出的模型有助于系统解析和模拟,开拓了解决科学问题的思路和方法,在各个领域得到了贴切的应用和深入的发展。
分数阶微分方程第三讲分数阶微分方程基本理论一、分数阶微分方程的出现背景及研究现状1、出现背景分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。
整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有较完善的理论。
但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到以下问题:(1)需要构造非线性方程,并引入一些人为的经验参数和与实际不符的假设条件;(2)因材料或外界条件的微小改变就需要构造新的模型;(3)这些非线性模型无论是理论求解还是数值求解都非常繁琐。
基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。
分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。
2、研究现状在近三个世纪里,对分数阶微积分理论的研究主要在数学的纯理论领域里进行,似乎它只对数学家们有用。
然而在近几十年来,分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。
分数阶微积分理论也受到越来越多的国内外学者的广泛关注,特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。
随着分数阶微分方程在越来越多的科学领域里出现,无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。
然而由于分数阶微分是拟微分算子,它的保记忆性(非局部性)对现实问题进行了优美刻画的同时,也给我们的分析和计算造成很大困难。
在理论研究方面,几乎所有结果全都假定了满足李氏条件,而且证明方法也和经典微积分方程一样,换句话说,这些工作基本上可以说只是经典微积分方程理论的一个延拓。
对分数阶微分方程的定性分析很少有系统性的结果,大多只是给出了一些非常特殊的方程的求解,且常用的求解方法都是具有局限性的。
在数值求解方面,现有分数阶方程数值算法还很不成熟,主要表现为:(1)在数值计算中一些挑战性难题仍未得到彻底解决,如长时间历程的计算和大空间域的计算等;(2)成熟的数值算法比较少,现在研究较多的算法主要集中在有限差分方法与有限单元法;(3)未形成成熟的数值计算软件,严重滞后于应用的需要。
鉴于此,发展新数值算法,特别是在保证计算可靠性和精度的前提下,提高计算效率,解决分数阶微分方程计算量和存储量过大的难点问题,发展相应的计算力学应用软件成为迫切需要关注的课题。
二、预备知识1、分数阶微积分经典定义回顾作为分数阶微积分方程的基础,本书在第二章中对分数阶微积分的定义及性质做了系统的介绍,为了接下来讨论的需要,我们首先对其进行一个简要的回顾。
(1)分数阶微积分的主要思想如上图所示,分数阶微积分的主要思想是推广经典的整数阶微积分,从而将微积分的概念延拓到整个实数轴,甚至是整个复平面。
但由于延拓的方法多种多样,因而根据不同的需求人们给出了分数阶微积分的不同定义方式。
然而这些定义方式不仅只能针对某些特定条件下的函数给出,而且只能满足人们的某些特定需求,迄今为止,人们仍然没能给出分数阶微积分的一个统一的定义, 这对分数阶微积分的研究与应用造成了一定的困难。
(2)几种经典的分数阶微积分定义下面我们试图从理论依据、定义域、表达式和优缺点几个方面给出常见的四种分数阶微积分定义的比较图。
依据1:整数阶微分的差分定义()1()lim(1)(),kk rkhrkf t f t rh k Nrh→=⎛⎫=--∈⎪⎝⎭∑G-L分数阶微积分定义:[]1[()]lim(1)()t anhG p ra t phrpD f t f t rhrh-=→=⎛⎫=--⎪⎝⎭∑1()1[()]()(),0()tC p n p na t aD f t t u f u du pn p--=->Γ-⎰11[()]()(),0()tC p pa t aD f t t u f u du pp--=->Γ⎰依据2:整数阶积分的柯西公式()11()()()()tn naf t t u f u dun--=-Γ⎰R-L分数阶微积分定义:()11[()]()()()tR n p n pa t aD f t t u f u dun p----=-Γ-⎰()[()][()],1,nR p R n pa t a tdD f t D f t n p n n Ndt--⎛⎫⎡⎤=-≤<∈⎪⎣⎦⎝⎭Caputo分数阶微积分定义:先积分后微分先微分后积分缺点:定义域较窄计算复杂优点:定义域较宽缺点:Laplace变换较复杂优点:Laplace较简洁缺点:定义域太窄分数阶微积分定义定义域:(){()}f tf t[a,b]在上逐段连续且在任意有限子区间上可积定义域:(){()}f tf t[a,b][p]+1在具有阶连续导数定义域:(){()}f tf t[a,b][p]+1在具有阶连续导数基于广义函数的分数阶微积分定义:()()()pa t pD f t f t t-=*Φ1,0()()0,0ppttt pt-⎧>⎪Φ=Γ⎨⎪≤⎩定义域:{()()0;}f t f t t a=<优点:有利于工程中对系统的描述从上图我们看到,在分数阶微积分的发展过程中,人们根据不同的需求,从不同角度给出了分数阶微积分的定义,但这些定义无论从对象上还是从表达式上都无法实现统一,它们之间的关系大致可以用下图来表示。
注:条件1:()f t 在[,]a b 上逐段连续,且在任何有限子区间上可积; 条件2:()f t 在[,]a b 上具有[]1p +阶连续导数; 条件3:()()0k f a =,1,2,,[]k p =; 条件4:()0f t =,t a <。
由上图我们可以看到,对于不同的分数阶微积分定义方式有着不同定义域,即便是在公共区域内,不同的定义方式之间也无法实现完全的统一,这对分数阶微积分的应用和研究造成了一定的困难,因此人们迫切期望着分数阶微积分的一种哪怕是形式上的统一定义方式。
2、 M-R 序列分数阶微分的定义为了满足实际需要,下面我们试图从形式上对分数阶微积分给出一种统一的表达式。
分数阶微积分的主要思想是推广经典的累次微积分,所有推广方法的共同目标是以非整数参数p 取代经典微积分符号中的整数参数n ,即:n n d dt ppddt实际上,任意的n 阶微分都可以看成是一列一阶微分的叠加:()()n n nd f t d d df t dt dt dt dt= (1) 由此,我们可以给出一种在很多实际应用中十分重要的分数阶微积分的推广方式。
首先,我们假设已有一种合适的推广方式来将一阶微分推广为α(01α≤≤)阶微分,即dD dtα→是可实现的。
那么类似地可得到(1)的推广式为:()()n nD f t D D D f t αααα= (2)这种推广方式最初是由..K S Miller 和.B Ross 提出来的,其中D α采用的是R L -分数阶微分定义,他们称之为序列分数阶微分。
序列分数阶微分的其他形式可以通过将D α替换为G L -分数阶微分、Caputo 分数阶微分或其他任意形式分数阶微分来得到。
进一步,如果我们将(2)中的分数阶微分D α替换为不同阶数的分数阶微分可得到序列分数阶微分更一般的表达式: 12()()n D f t D D D f t αααα= (3) 12n αααα=+++根据问题的需要,D α可以是R L -分数阶微分、G L -分数阶微分、Caputo分数阶微分或其他任意形式的分数阶微分,从这一点看来,我们可以说序列分数阶微分从形式上给出了分数阶微积分在时域上的一个统一表达式,R L -分数阶微分、G L -分数阶微分和Caputo 分数阶微分都只是序列分数阶微分的一种特殊情况。
故而,下面我们在对分数阶微积分方程进行理论分析的时候可以仅仅针对序列分数阶微积分来给出结论。
3、M-R 序列分数阶微分的Laplace 变换下面我们考虑如下形式的序列分数阶微分的Laplace 变换。
11m m m a t a t a t at D D D D σααα-= (4)1111m m m at a t a t at D D D D σααα---= (5)101,(1,2,,)mm jj j j m σαα==<≤=∑ (6)在R-L 分数阶微分定义下有:1000{();}()[()]t t t L D f t s s F s D f t ααα-==- (7)重复利用上式m 次可得:110000{();}()[()]mmm m k m k m t t t k L D f t s s F s s D f t σσσσσ-----===-∑ (8)注:虽然上述序列分数阶微分的Laplace 变换是在R-L 分数阶微分定义下进行证明的,但是该结论对其他几种分数阶微积分也是成立的。
4、泛函理论基础定理1(Schauder 不动点定理)设U 是Banach 空间X 的有界闭子集,如果:T U U →是连续映射,那么T 在U 中存在不动点,即使得Tx x =的点存在。
定义1(Lipschitz 条件)设,X d 是距离空间,T 是从X 到X 的映射,如果存在常数0q >,使得对所有的,x y X ∈,(,)(,)d Tx Ty qd x y ≤则称T 满足Lispschitz 条件,q 成为T 的Lispschitz 常数。
特别的,如果1q <,则T 称为压缩映射。
定理2( Banach 压缩映像原理)设,X d 是距离空间,:T X X →是压缩映射,则T 在X 中恰有一个不动点。
设这个不动点为x ,则对任何初始点0x X ∈,逐次迭代点列1n n x Tx +=,1,2,n =收敛于x ,且关于收敛速度有如下估计式:100(,)(1)(,)n n d x x q q d Tx x -≤- 其中,q 是T 的Lipschitz 常数。
三、 解的存在唯一性理论近年来,分数阶微分方程已经在国内外引起极大的研究兴趣,尤其是关于其解的性质的研究,诸如存在性及唯一性等,其中大多数的研究方法是通过把分数阶初值问题转换成等价的分数阶积分方程,然后运用不动点定理来得到分数阶初值问题解的存在唯一性结果。
已有研究结果主要有以下限制:(1) 函数的定义区间为有限区间[,]a b ; (2) 函数在定义域上需满足Lipschitz 条件; 因此,目前人们在这方面所做的工作都是希望设法在放宽上述两个限制条件后给出分数阶微积分方程的解的存在唯一性定理。