分数阶微分方程
第三讲 分数阶微分方程基本理论
1、 分数阶微分方程的出现背景及研究现状
1、出现背景
分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。
整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有较完善的理论。但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到以下问题:
(1) 需要构造非线性方程,并引入一些人为的
经验
班主任工作经验交流宣传工作经验交流材料优秀班主任经验交流小学课改经验典型材料房地产总经理管理经验
参数和与实际不符的假
设条件;
(2) 因材料或外界条件的微小改变就需要构造新的模型;
(3) 这些非线性模型无论是理论求解还是数值求解都非常繁琐。
基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。
2、研究现状
在近三个世纪里,对分数阶微积分理论的研究主要在数学的纯理论领域里进行,似乎它只对数学家们有用。然而在近几十年来,分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。分数阶微积分理论也受到越来越多的国内外学者的广泛关注,特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。随着分数阶微分方程在越来越多的科学领域里出现,无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。然而由于分数阶微分是拟微分算子,它的保记忆性(非局部性)对现实问题进行了优美刻画的同时,也给我们的分析和计算造成很大困难。
在理论研究方面,几乎所有结果全都假定了满足李氏条件,而且证明方法也和经典微积分方程一样,换句话说,这些工作基本上可以说只是经典微积分方程理论的一个延拓。对分数阶微分方程的定性分析很少有系统性的结果,大多只是给出了一些非常特殊的方程的求解,且常用的求解方法都是具有局限性的。
在数值求解方面,现有分数阶方程数值算法还很不成熟,主要表现为:
(1)在数值计算中一些挑战性难题仍未得到彻底解决,如长时间历程的计算和大空间域的计算等;
(2)成熟的数值算法比较少,现在研究较多的算法主要集中在有限差分方法与有限单元法;
(3)未形成成熟的数值计算软件,严重滞后于应用的需要。
鉴于此,发展新数值算法,特别是在保证计算可靠性和精度的前提下,提高计算效率,解决分数阶微分方程计算量和存储量过大的难点问题,发展相应的计算力学应用软件成为迫切需要关注的课题。
2、 预备知识
1、 分数阶微积分经典定义回顾
作为分数阶微积分方程的基础,本
书
关于书的成语关于读书的排比句社区图书漂流公约怎么写关于读书的小报汉书pdf
在第二章中对分数阶微积分的定义及性质做了系统的介绍,为了接下来讨论的需要,我们首先对其进行一个简要的回顾。
(1)分数阶微积分的主要思想
如上图所示,分数阶微积分的主要思想是推广经典的整数阶微积分,从而将微积分的概念延拓到整个实数轴,甚至是整个复平面。但由于延拓的方法多种多样,因而根据不同的需求人们给出了分数阶微积分的不同定义方式。然而这些定义方式不仅只能针对某些特定条件下的函数给出,而且只能满足人们的某些特定需求,迄今为止,人们仍然没能给出分数阶微积分的一个统一的定义, 这对分数阶微积分的研究与应用造成了一定的困难。
(2)几种经典的分数阶微积分定义
下面我们试图从理论依据、定义域、表达式和优缺点几个方面给出常见的四种分数阶微积分定义的比较图。
从上图我们看到,在分数阶微积分的发展过程中,人们根据不同的需求,从不同角度给出了分数阶微积分的定义,但这些定义无论从对象上还是从表达式上都无法实现统一,它们之间的关系大致可以用下图来表示。
注:
条件1:
在
上逐段连续,且在任何有限子区间上可积;
条件2:
在
上具有
阶连续导数;
条件3:
,
;
条件4:
,
。
由上图我们可以看到,对于不同的分数阶微积分定义方式有着不同定义域,即便是在公共区域内,不同的定义方式之间也无法实现完全的统一,这对分数阶微积分的应用和研究造成了一定的困难,因此人们迫切期望着分数阶微积分的一种哪怕是形式上的统一定义方式。
2、 M-R序列分数阶微分的定义
为了满足实际需要,下面我们试图从形式上对分数阶微积分给出一种统一的表达式。
分数阶微积分的主要思想是推广经典的累次微积分,所有推广方法的共同目标是以非整数参数
取代经典微积分符号中的整数参数
,即:
实际上,任意的
阶微分都可以看成是一列一阶微分的叠加:
(1)
由此,我们可以给出一种在很多实际应用中十分重要的分数阶微积分的推广方式。首先,我们假设已有一种合适的推广方式来将一阶微分推广为
(
)阶微分,即
是可实现的。那么类似地可得到(1)的推广式为:
(2)
这种推广方式最初是由
和
提出来的,其中
采用的是
分数阶微分定义,他们称之为序列分数阶微分。序列分数阶微分的其他形式可以通过将
替换为
分数阶微分、
分数阶微分或其他任意形式分数阶微分来得到。
进一步,如果我们将(2)中的分数阶微分
替换为不同阶数的分数阶微分可得到序列分数阶微分更一般的表达式:
(3)
根据问题的需要,
可以是
分数阶微分、
分数阶微分、
分数阶微分或其他任意形式的分数阶微分,从这一点看来,我们可以说序列分数阶微分从形式上给出了分数阶微积分在时域上的一个统一表达式,
分数阶微分、
分数阶微分和
分数阶微分都只是序列分数阶微分的一种特殊情况。故而,下面我们在对分数阶微积分方程进行理论分析的时候可以仅仅针对序列分数阶微积分来给出结论。
3、M-R序列分数阶微分的Laplace变换
下面我们考虑如下形式的序列分数阶微分的Laplace变换。
(4)
(5)
(6)
在R-L分数阶微分定义下有:
(7)
重复利用上式
次可得:
(8)
注:虽然上述序列分数阶微分的Laplace变换是在R-L分数阶微分定义下进行证明的,但是该结论对其他几种分数阶微积分也是成立的。
4、泛函理论基础
定理1(Schauder不动点定理)
设
是
空间
的有界闭子集,如果
是连续映射,那么
在
中存在不动点,即使得
的点存在。
定义1(Lipschitz条件)
设
是距离空间,
是从
到
的映射,如果存在常数
,使得对所有的
,
则称
满足
条件,
成为
的
常数。
特别的,如果
,则
称为压缩映射。
定理2( Banach压缩映像原理)
设
是距离空间,
是压缩映射,则
在
中恰有一个不动点。设这个不动点为
,则对任何初始点
,逐次迭代点列
,
收敛于
,且关于收敛速度有如下估计式:
其中,
是
的
常数。
3、 解的存在唯一性理论
近年来,分数阶微分方程已经在国内外引起极大的研究兴趣,尤其是关于其解的性质的研究,诸如存在性及唯一性等,其中大多数的研究方法是通过把分数阶初值问题转换成等价的分数阶积分方程,然后运用不动点定理来得到分数阶初值问题解的存在唯一性结果。已有研究结果主要有以下限制:
(1) 函数的定义区间为有限区间
;
(2) 函数在定义域上需满足
条件;
因此,目前人们在这方面所做的工作都是希望设法在放宽上述两个限制条件后给出分数阶微积分方程的解的存在唯一性定理。
下面我们对分数阶微分方程初值问题的现有理论结果作一个简单的介绍,相应的结论都是针对定义在有限区间
上的M-R序列分数阶微分形式,在满足
条件下给出的,当然,由前面的介绍可知,这些结论也可直接推广到其他分数阶微分形式。
1、 线性分数阶微分方程解的存在唯一性定理
考虑如下形式的初值问题:
且
,即
(11)
第一步:假设
,考虑由此得到的退化问题解的存在唯一性。
定理1 如果
,则方程
(12)
有满足初值条件(10)的唯一解
。
定理的证明过程如下:
步骤一 通过Laplace变换证明解的存在性;
下面我们设法构造一个待求解问题解,对式(12)做Laplace变换可得:
(13)
其中,
、
分别是
、
的Laplace变换。利用初值条件(10)可得:
(14)
对上式做Laplace逆变换可得:
(15)
步骤二 由分数阶微分的线性性和Laplace变换的性质证明唯一性。
假设有存在两个满足上述初值问题的解
、
令
,有分数阶微分方程的线性性可得:
(16)
从而有
(17)
由Laplace变换的性质可知:
在
上几乎处处成立。
故原方程的解在
上唯一。
注:上述证明过程中用到的Laplace变换法是一种常用的分数阶微分方程求解方法,该方法步骤简单,适用范围较广,在实际中有着重要应用,后面将对其进行详细介绍。
第二步:运用第一步的结论证明原初值问题解的存在唯一性。
定理2 如果
且
是
上的连续函数,则初值问题(9)—(10)有唯一解
。
定理的证明过程如下:
步骤一 化微分方程为积分方程
假设原方程有解
并记
,那么运用定理1可得:
(18)
将上式代入到原微分方程表达式(9)可得:
(19)
其中
(20)
(21)
步骤二 证明变换后的积分方程有唯一解
用不动点定理易证结论成立。
步骤三 说明原微分方程有唯一解
由定理1易得。
2、 一般形式的分数阶微分方程的存在唯一性定理
考虑如下形式的微分方程:
(22)
,
(23)
其中,
的定义域为平面
上的一个子区域
,且存在
上的子区域
满足:
,
(24)
定理3 设
为
上的连续实值函数,且在
上关于
满足
条件,即
(25)
从而
,对任意
且
那么,方程(22)—(23)在区域
有唯一的连续解。
定理的证明过程如下:
步骤一 化微分方程为等价积分方程;
对方程(22)按
,
,
,
逐次进行分部积分可得:
(26)
步骤二 证明上述等价积分方程解的存在性;
构造函数序列
,
,
,
如下:
(27)
(28)
首先,我们可以证明对任意的
及任意的
有
。
(29)
进一步,我们可由数学归纳法证明,对任意的
有下式成立:
(30)
证明过程如下:
在式(29)中令
可得:
(31)
假设当
时,式(30)成立,即下式成立:
(32)
那么,当
时有:
(33)
从而由归纳法可知,对任意的
,式(30)成立。
进而,有
的收敛性可知,函数序列
收敛。
令
,易证
是等价积分方程(26)的解,也即是原微分方程(22)的解。
步骤三 证明上述等价积分方程解的唯一性;
假设
也是等价积分方程(26)的解,令
,有:
(34)
由
的连续性可知,存在常数
,使得对任意的
,
。
利用式(34)可得:
(35)
将该估计过程重复
次可得:
(36)
又
故
,也即
。
注:有上面的介绍可知,整个线性分数阶微积分方程解的存在唯一性理论的证明过程都是建立在不动点理论的基础上的,使得我们必须将讨论范围限制在有限区间内的满足Lipschitz条件的函数上,如何打破这个限制是一个值得思考的问题。在某些情况下,定理3可直接作为分数阶微分方程的求解方法,通常称之为存在唯一性解法。
由上面的介绍,我们可将分数阶微积分方程存在唯一性理论及其所面临的问题描述如下:
3、 分数阶微分方程初值问题解的依赖性
下面我们来考察初值条件的微小变化将对方程的解造成怎样的影响,为此,我们在初值条件中引入一个微小的改变量。
,
(37)
其中
为任意常数。
定理4 设
是初值问题(22)—(23)的解,
是初值问题(22)、(37)的解,那么对任意的
有:
(38)
其中
为
函数。
证明:
步骤一 用定理3的方式构造两组函数序列
,
,
,
和
,
,
,
使得
,
。
步骤二 由数学归纳法容易证明
(39)
步骤三 对上式两端取极限
可得:
(40)
4、 Laplace变换求解法
随着分数阶微分方程在
工程
路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理
应用中出现得越来越频繁,给出分数阶微分方程的有效而简便的求解方法便显得越来越重要,然而现有的求解方法都有着各种各样的缺陷。下面我们介绍一种基于Laplace变换的分数阶微分方程求解方法,该方法简单、直观,适用于常系数线性分数阶微分方程的求解。
1、 Laplace变换求解法
(1) Laplace变换求解法的主要步骤
步骤一:对原微分方程做Laplace变换,化微分方程为代数方程;
步骤二:求解该代数方程,得到原问题在变换域上的解;
步骤三:对该变化域上的解做Laplace逆变换得到原问题的时域解。
(2) Laplace变换求解法的应用
下面我们通过两个例子来说明Laplace变换法的应用方法。
例1 我们考虑用Laplace变换法对如下的非齐次标准分数阶微分方程的初值问题进行求解。
,
(41)
,
(42)
其中,
解:对方程(41)两端做Laplace变换,并利用初值条件(42)可得:
从而
(43)
对式(31)做Laplace逆变换可得原微分方程的解为:
(44)
注:某些文献中也给出了该问题用迭代法进行求解的过程,虽然两种解法的结果相同,但显然Laplace求解法更为直观、简便。
例2 下面我们考虑用Laplace变换法对序列分数阶微分方程的初值问题进行求解。
(45)
,
(46)
解:对方程(45)两端做Laplace变换,并利用初值条件(46)可得:
从而
(47)
对式(35)做Laplace逆变换可得原微分方程的解为:
(48)
其中,
。
注:对比上面两个初值问题容易看到他们在形式上非常相似,唯一的差别体现在一个是基于经典分数阶微积分定义的标准分数阶微分方程,一个是基于序列分数阶微积分定义的序列分数阶微分方程,从而在初值地给法不一样。但我们发现它们的解在表达式上也非常地相近,对比结果如下:
通过上面的对比可以发现,对应的标准分数阶微分方程和序列分数阶微分方程的解有一个共同点,即它们具有同样的
函数,下面我们就
展开讨论。
2、 Green函数
考虑如下的初值问题:
,
,
(49)
其中
(1) 定义
若函数
满足如下条件,则称其为方程(37)的
函数:
1)
对任意的
;
2)
,
(
是Kronecker delta函数);
3)
,
。
(2) 性质
1)
是方程(49)的解;
2)对常系数分数阶微分方程有:
;
3)对
的适当微分可得到一组齐次方程(
)的线性无关解。
下面我们利用
函数的定义来证明上述三个性质。
证明1):计算
如下:
(50)
将上述等式所表示的
,
,
,
累加起来有:
(51)
证明2):由Laplace求解法可得。
证明3):第一步,取
,令
,由交换律的条件及定义3)有:
(52)
故
是原齐次方程的解。
第二步,由结合律可知:
(53)
有上面的分析可知,对线性常系数分数阶微分方程的求解就转换为寻找该方程的分数阶
函数,有该
函数便可直接写出该方程的解,解的具体表达式如下:
(54)
下面我们给出寻找一般的线性常系数分数阶微分方程的
函数的方法。
(3) 求法
对于常系数分数阶微分方程而言,我们可以通过对其做
变换来得到相应的
函数,下面我们通过一个例子来对其进行说明。
例3 考虑如下的一阶常系数微分方程:
(55)
解:由于初值对
函数无影响,因此我们对上式做
变换并忽略初值部分可得:
(56)
从而
(56)
故
(57)
对其做Laplace逆变换可得:
(58)
注:对如下的
阶分数阶微分方程:
(59)
假设
可得:
(60)
(61)
五 小结
有关分数阶微分方程的理论分析部分我们主要介绍了两方面的内容,一是分数阶微分方程解的性质,一是分数阶微分方程的求解方法。由于对分数阶微分方程的研究还不够成熟,因此对其所做出的理论分析还处于探索阶段。已有成果多半是对经典微积分方程理论的简单推广,且只能覆盖部分特殊形式的分数阶微分方程。现有的很多工作都是试图寻找新的理论方法,以打破现有的限制条件,力求构建一套完善的分数阶微分方程理论。
_1349036716.unknown
_1349036780.unknown
_1349036813.unknown
_1349036829.unknown
_1349036837.unknown
_1349036845.unknown
_1349036849.unknown
_1349036851.unknown
_1349036853.unknown
_1349036854.unknown
_1349036855.unknown
_1349036852.unknown
_1349036850.unknown
_1349036847.unknown
_1349036848.unknown
_1349036846.unknown
_1349036841.unknown
_1349036843.unknown
_1349036844.unknown
_1349036842.unknown
_1349036839.unknown
_1349036840.unknown
_1349036838.unknown
_1349036833.unknown
_1349036835.unknown
_1349036836.unknown
_1349036834.unknown
_1349036831.unknown
_1349036832.unknown
_1349036830.bin
_1349036821.unknown
_1349036825.bin
_1349036827.unknown
_1349036828.unknown
_1349036826.unknown
_1349036823.unknown
_1349036824.unknown
_1349036822.unknown
_1349036817.unknown
_1349036819.unknown
_1349036820.unknown
_1349036818.unknown
_1349036815.unknown
_1349036816.unknown
_1349036814.unknown
_1349036796.unknown
_1349036805.unknown
_1349036809.unknown
_1349036811.unknown
_1349036812.unknown
_1349036810.unknown
_1349036807.unknown
_1349036808.unknown
_1349036806.unknown
_1349036801.bin
_1349036803.bin
_1349036804.unknown
_1349036802.unknown
_1349036799.unknown
_1349036800.unknown
_1349036798.unknown
_1349036788.unknown
_1349036792.unknown
_1349036794.unknown
_1349036795.unknown
_1349036793.unknown
_1349036790.unknown
_1349036791.unknown
_1349036789.unknown
_1349036784.unknown
_1349036786.unknown
_1349036787.unknown
_1349036785.unknown
_1349036782.unknown
_1349036783.unknown
_1349036781.unknown
_1349036748.unknown
_1349036764.unknown
_1349036772.unknown
_1349036776.unknown
_1349036778.unknown
_1349036779.unknown
_1349036777.unknown
_1349036774.unknown
_1349036775.unknown
_1349036773.unknown
_1349036768.unknown
_1349036770.unknown
_1349036771.unknown
_1349036769.unknown
_1349036766.unknown
_1349036767.unknown
_1349036765.unknown
_1349036756.unknown
_1349036760.unknown
_1349036762.bin
_1349036763.unknown
_1349036761.unknown
_1349036758.unknown
_1349036759.unknown
_1349036757.unknown
_1349036752.unknown
_1349036754.unknown
_1349036755.unknown
_1349036753.unknown
_1349036750.unknown
_1349036751.unknown
_1349036749.unknown
_1349036732.unknown
_1349036740.unknown
_1349036744.unknown
_1349036746.unknown
_1349036747.unknown
_1349036745.unknown
_1349036742.unknown
_1349036743.unknown
_1349036741.unknown
_1349036736.unknown
_1349036738.unknown
_1349036739.unknown
_1349036737.unknown
_1349036734.unknown
_1349036735.unknown
_1349036733.unknown
_1349036724.unknown
_1349036728.unknown
_1349036730.unknown
_1349036731.unknown
_1349036729.unknown
_1349036726.unknown
_1349036727.unknown
_1349036725.unknown
_1349036720.unknown
_1349036722.unknown
_1349036723.unknown
_1349036721.unknown
_1349036718.unknown
_1349036719.unknown
_1349036717.unknown
_1349036652.unknown
_1349036684.unknown
_1349036700.unknown
_1349036708.unknown
_1349036712.unknown
_1349036714.unknown
_1349036715.unknown
_1349036713.unknown
_1349036710.unknown
_1349036711.unknown
_1349036709.unknown
_1349036704.unknown
_1349036706.unknown
_1349036707.unknown
_1349036705.unknown
_1349036702.unknown
_1349036703.unknown
_1349036701.unknown
_1349036692.unknown
_1349036696.unknown
_1349036698.unknown
_1349036699.unknown
_1349036697.unknown
_1349036694.unknown
_1349036695.unknown
_1349036693.unknown
_1349036688.unknown
_1349036690.unknown
_1349036691.unknown
_1349036689.unknown
_1349036686.unknown
_1349036687.unknown
_1349036685.unknown
_1349036668.unknown
_1349036676.unknown
_1349036680.unknown
_1349036682.unknown
_1349036683.unknown
_1349036681.unknown
_1349036678.unknown
_1349036679.unknown
_1349036677.unknown
_1349036672.unknown
_1349036674.unknown
_1349036675.unknown
_1349036673.unknown
_1349036670.unknown
_1349036671.unknown
_1349036669.unknown
_1349036660.unknown
_1349036664.unknown
_1349036666.unknown
_1349036667.unknown
_1349036665.unknown
_1349036662.unknown
_1349036663.unknown
_1349036661.unknown
_1349036656.unknown
_1349036658.unknown
_1349036659.unknown
_1349036657.unknown
_1349036654.unknown
_1349036655.unknown
_1349036653.unknown
_1349036619.unknown
_1349036635.unknown
_1349036643.unknown
_1349036647.unknown
_1349036649.unknown
_1349036650.unknown
_1349036648.unknown
_1349036645.unknown
_1349036646.unknown
_1349036644.unknown
_1349036639.unknown
_1349036641.unknown
_1349036642.unknown
_1349036640.unknown
_1349036637.unknown
_1349036638.unknown
_1349036636.unknown
_1349036627.unknown
_1349036631.unknown
_1349036633.unknown
_1349036634.unknown
_1349036632.unknown
_1349036629.unknown
_1349036630.unknown
_1349036628.unknown
_1349036623.unknown
_1349036625.unknown
_1349036626.unknown
_1349036624.unknown
_1349036621.unknown
_1349036622.unknown
_1349036620.unknown
_1349036603.unknown
_1349036611.unknown
_1349036615.unknown
_1349036617.unknown
_1349036618.unknown
_1349036616.unknown
_1349036613.unknown
_1349036614.unknown
_1349036612.unknown
_1349036607.unknown
_1349036609.unknown
_1349036610.unknown
_1349036608.unknown
_1349036605.unknown
_1349036606.unknown
_1349036604.unknown
_1349036595.unknown
_1349036599.unknown
_1349036601.unknown
_1349036602.unknown
_1349036600.bin
_1349036597.unknown
_1349036598.unknown
_1349036596.unknown
_1349036591.unknown
_1349036593.unknown
_1349036594.unknown
_1349036592.unknown
_1349036589.unknown
_1349036590.unknown
_1349036588.bin