浅谈常微分方程的数值解法及其应用文献综述Word文件下载.docx
《浅谈常微分方程的数值解法及其应用文献综述Word文件下载.docx》由会员分享,可在线阅读,更多相关《浅谈常微分方程的数值解法及其应用文献综述Word文件下载.docx(10页珍藏版)》请在冰豆网上搜索。
.
(1)
为常微分方程.其中出现的最高阶导数的阶数,叫做常微分方程的阶.例如
,
,是一阶常微分方程.是二阶常微分方程.设定义于
区间上,有直到阶的导数,将它代入
(1),使
(1)变成关于的恒等式,即
.
就称=为
(1)的一个定义于上的解,并称为该解的定义区间.[4]
如果一个微分方程中出现多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程.
2.2微分方程产生的历史背景
微分方程差不多是和微积分同时先后产生的,苏格兰数学家耐普尔创立对数的时候,就讨论过微分方程的近似解.牛顿在建立微积分的同时,对简单的微分方程用级数来求解。
后来瑞士数学家雅各布•贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。
微分方程的形成与发展是和力学、天文学、物理学,以及其他科学技术的发展密切相关的.数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对微分方程的发展产生了深刻的影响,当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具.[5]
牛顿研究天体力学和机械力学的时候,利用了微分方程这个工具,从理论上得到了行星运动规律.后来,法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置.这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量.
微分方程的理论逐步完善的时候,利用它就可以精确地表述事物变化所遵循的基本规律,只要列出相应的微分方程,有了解方程的方法.微分方程也就成了最有生命力的数学分支.总之,力学、天文学、几何学等领域的许多问题都导致微分方程.在当代,甚至许多社会科学的问题亦导致微分方程,如人口发展模型、交通流模型等.因而微分方程的研究是与人类社会密切相关的.[6]
2.3微分方程发展现状及其基本功能
在数学学科内部的许多分支中,微分方程是常用的重要工具之一,微分方程进一步发展的需要,有推动着其它数学分支的发展;
相反,微分方程每一步进展都离不开其他数学分支的支援.数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对微分方程的发展产生了深刻的影响.当前计算机的发展更是为微分方程的应用及理论研究提供了非常有力的工具.时至今日,可以说微分方程在所有自然科学领域和众多社会科学领域都有着广泛的应用,如自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等.只要能够列出相应的微分方程,有了解方程的方法,利用它就可以精确地表述事物变化所遵循的基本规律.从微积分理论形成以来,人们一直用微分方程来描述、解释或预见各种自然现象,不断的取得了显著的成效.[7]
2.4常微分方程的数值求解方法
2.4.1Euler法
Euler法是最简单的数值方法,为求解良态初值问题,的区间。
实际上,下面的过程不是要找到满足该初值问题的可微函数,而是要生成点集,并且将这些点作为近似解,即。
如何构造“近似满足微方程”的“点集”呢?
首先为这些点选择横坐标,为方便起见,将区间划分为个等距子区间,并选择网络点
,k=0,1,……,其中
(1)
值称为步长。
然后近似解
在上,
(2)
设,和连续,;
;
利用泰勒定理将在处展开,对每个值,存在一个和之间的值,使得
(3)
将和代人等式(3),得到的表示:
(4)
如果步长h足够小,则可以忽略2次项(包含的项),得到
(5)
这就是欧拉近似。
重复该过程,就能得到近似解曲线的一个点序列。
欧拉方法的一般步骤是
,其中k=0,1,……,M-1[8](6)
2.4.2泰勒级数法
泰勒级数法有着广泛的应用,并且是比较求解初值问题的各种不同数值方法的标准,它可设计为任意指定的精度。
下面首先将泰勒定理用新的公式表示,使之适合于求解微分方程。
定理9.5(泰勒定理)设,且在不动点处有N次泰勒级数展开:
(1)
其中,
(2)
表示函数关t的()次全导数。
求导公式可以递归地计算:
并且一般有
其中为导数算子
区间上的初值问题的近似数值解可由各子区间上的公式
(1)来推导。
次泰勒方法的一般步骤为
其中在各步有。
次泰勒方法的最终全局误差是阶的,因此可选择所需大小的,使得误差足够小。
如果是固定,则理论上可以推导出步长,使之满足任意的最终全局误差。
然而在实际运算中,通常用和计算两个近似结果集,然后比较其结果[9]。
2.4.3龙格—库塔方法
泰勒方法的优点是最终全局误差的阶为,并且可以通过选择较大的N来得到较小的误差。
然而泰勒方法的缺点是,需要先确定N,并且要计算高阶导数,它们可能十分复杂。
每个龙格一库塔(Runge-Kutta)方法都由一个合适的泰勒方法推导而来,使得其最终全局误差为。
一种折中方法是每步进行若干次函数求值,从而省去高阶导数计算。
这种方法可构造任意N阶精度的近似公式。
最常用的是N=4的龙格一库塔方法,它适用于一般的应用,因为它非常精确、稳定,且易于编程。
许多专家声称,没有必要使用更高阶的方法,因为提高的精度与增加的计算量相抵消。
如果需要更高的精度,则应该使用更小的步长或某种自适应方法。
4阶龙格一库塔方法(RK4)可模拟N=4的泰勒方法的精度。
它基于如下对,的计算:
其中,,和形如
通过与N=4阶的泰勒级数方法的系数匹配,使得局部误差为,龙格和库塔得出了如下方程组:
该方程组有11个方程和13个未知量,必须补充两个条件才可以求解。
最有用的选择是
,(4)
其余变量的解为
将式(4)和(5)中的值代入式
(2)和式
(1),得到标准的阶龙格—库塔方法,其描述如下。
自初始点开始,利用
(6)
生成近似值序列,其中[10]
(7)
2.4.4预报—校正方法
欧拉方法、休恩方法、泰勒方法以及龙格一库塔方法都称为单步长方法,因为它们只利用前一个点的信息来计算下一个点,即计算时只使用了初始点。
一般地,只有用来。
当计算出若干个点之后,就可以利用几个已计算出的点来计算下一个点。
以亚当斯一巴什福斯4步法的推导为例,计算需要,,和。
该方法不是自启动的,要生成点,必须先给出其4个初始点,,,(可用前面各节中的方法完成)。
多步法的一个优点是,可以确定它的局部截断误差(localtruncationerror,简称L.T.E.),并可以包含一个校正项,用于在每一步计算中改善解的精确度。
该方法还可以确定步长是否小到能得到的精确值,同时又大到能够免除不必要的和费时的计算。
使用预报子和校正子的组合在每一步只需要进行两次函数求值[11]。
2.4.4.1亚当斯一巴什福斯一莫尔顿方法
亚当斯一巴什福斯一莫尔顿方法(Adams—Bashforth—Moulton)是由基本微积分定理推导出的多步法:
预报子使用基于点和的的拉格朗日多项式逼近值,并在区间上对式
(1)积分,这个过程产生亚当斯一巴什福斯预报子:
(2)
校正子的推导类似。
这时可以实用刚刚计算出的值。
基于点,和新的点构造的一个新的拉格朗日多项式逼近,然后在区间上对该多项式积分,即可得到亚当斯一莫尔顿校正子[12]:
2.4.4.2米尔恩—辛普森方法
米尔恩—辛普森方法是预报子基于区间上的对的积分:
预报子使用基于和的拉格朗日多项式逼近,在区间上对它积分,得到米尔恩预报子:
此时值已知,基于点,和新点构造的新的拉格朗日多项式,然后在区间上对该多项式积分,结果为大家所熟悉的辛普森公式[13]:
2.5常微分方程的应用
2.5.1微分方程在传染病模型中的应用
随着卫生设施的改善、医疗水平的提高以及人类文明的不断发展,诸如霍乱、天花、等曾经肆虐全球的传染性疾病已经得到有效的控制。
但是一些新的、不断变异着的传染病毒却悄悄向人类袭来。
20世纪80年代四分险恶的艾滋病毒开始肆虐全球,至今仍在蔓延;
2003年春来历不明的非典病毒突袭人间,给人们的生命财产带来极大的危害。
长期以来,建立传染病的数学模型来描述传染病的传播过程,分析受感染人数的变化规律,探索制止传染病蔓延的手段等,一直是各国有关专家和官员关注的问题。
这里我们介绍的是按照一般传播机理建立的数学模型如下。
,
,
(1)
,
对上述常微分方程组进行分析求解,就可以了解不同时刻传染病的动力学特征(比如:
传染病病情的发展趋势,即各类人的人口数量的分布情况)[14].
2.5.2常微分方程在捕鱼业的应用
设某水域现有鱼量,由于受资源限制所能容纳的最大鱼量,高自然增长率,捕捞增长率,按人口的逻辑模型建立微分方程。
要保持鱼量平衡,设平衡点为,解得
设考虑在的泰勒展式
当>0时与同号为不稳定平衡点
当<0时与异号为稳定平衡点
<0即>
设由于<
曲线与有交点,因在原点切线为
解得,易知当时,取得最大捕捞量
,
最大捕捞量为[15]
三、总结部分
微分方程是在解决实际问题的过程中产生的,微分方程的研究又促进实际问题的解决,同时也促进其他学科的发展.微分方程在物理、工程、力学、天文学、生物学、医学、经济学等诸多领域都有重要作用.目前,微分方程的实际背景广、应用性强的特点已受到广泛关注.由于同一类型的微分方程往往可以用来描述许多性质上颇不相同的自然现象,对一些重要的微分方程开展研究,