浅谈常微分方程的数值解法及其应用[文献综述]

毕业论文文献综述信息与计算科学浅谈常微分方程的数值解法及其应用一、前言部分微分方程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解.牛顿在建立微积分的同时，对简单的微分方程用级数来求解.后来瑞士数学家雅各布•贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论.微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法.微分方程也就成了最有生命力的数学分支.总之，力学、天文学、几何学等领域的许多问题都导致微分方程.在当代，甚至许多社会科学的问题亦导致微分方程，如人口发展模型、交通流模型等.因而微分方程的研究是与人类社会密切相关的. [1]“常微分方程”是理学院数学系所有专业学生的重要专业基础课之一，也是工科、经济等专业必学内容之一.其重要性在于它是各种精确自然科学、社会科学中表述基本定律和各种问题的根本工具之一，换句话说，只要根据实际背景，列出了相应的微分方程，并且能（数值地或定性地）求出这种方程的解，人们就可以预见到，在已知条件下这种或那种“运动”过程将怎样进行，或者为了实现人们所希望的某种“运动”应该怎样设计必要的装置和条件等等.例如，我们要设计人造卫星轨道，首先，根据力学原理，建立卫星运动的微分方程，列出初始条件，然后求出解，即卫星运行轨道.随着物理科学所研究的现象在广度和深度两方面的扩展，微分方程的应用范围更广泛. [2]从数学自身的角度看，微分方程的求解促使数学在函数论、变分法、级数展开、常微分方程、代数、微分几何等各方面进行发展.从这个角度说，微分方程变成了数学的中心. [3]总之，微分方程从它诞生起即日益成为人类认识并进而改造自然、社会的有力工具，成为数学科学联系实际的主要途径之一.文章就常微分的数值解法以及应用展开简单的论述。

二、主体部分2.1微分方程概念介绍2.1.1 微分方程概况由一元函数得到的方程.即：称含有自变量，未知函数及其导数的关系式22(,,,,...,)0n n dy d y d y F x y dx dx dx=. （1） 为常微分方程.其中出现的最高阶导数的阶数，叫做常微分方程的阶.例如 dy dx=x ，dy y dx = ，是一阶常微分方程. 22sin 0d g dt pθθ+=是二阶常微分方程.设)(x y ϕ=定义于 区间J 上，有直到n 阶的导数，将它代入（1），使（1）变成关于x 的恒等式，即()()(,(),,...,)0,n n d x d x F x x x J dx dxϕϕϕ=∈. 就称y ＝()x ϕ为（1）的一个定义于J 上的解，并称J 为该解的定义区间. [4]如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程.2.2微分方程产生的历史背景微分方程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解.牛顿在建立微积分的同时，对简单的微分方程用级数来求解。

后来瑞士数学家雅各布•贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。

微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的.数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具. [5]牛顿研究天体力学和机械力学的时候，利用了微分方程这个工具，从理论上得到了行星运动规律.后来，法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置.这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量.微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法.微分方程也就成了最有生命力的数学分支.总之，力学、天文学、几何学等领域的许多问题都导致微分方程.在当代，甚至许多社会科学的问题亦导致微分方程，如人口发展模型、交通流模型等.因而微分方程的研究是与人类社会密切相关的. [6]2.3 微分方程发展现状及其基本功能在数学学科内部的许多分支中，微分方程是常用的重要工具之一，微分方程进一步发展的需要，有推动着其它数学分支的发展；相反，微分方程每一步进展都离不开其他数学分支的支援.数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对微分方程的发展产生了深刻的影响.当前计算机的发展更是为微分方程的应用及理论研究提供了非常有力的工具.时至今日，可以说微分方程在所有自然科学领域和众多社会科学领域都有着广泛的应用，如自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等.只要能够列出相应的微分方程，有了解方程的方法，利用它就可以精确地表述事物变化所遵循的基本规律.从微积分理论形成以来，人们一直用微分方程来描述、解释或预见各种自然现象，不断的取得了显著的成效.[7] 2.4常微分方程的数值求解方法2.4.1 Euler 法Euler 法是最简单的数值方法，[,]a b 为求解良态初值问题'(,)y f t y =，0()y a y =的区间。

实际上，下面的过程不是要找到满足该初值问题的可微函数，而是要生成点集{(,)}k k t y ，并且将这些点作为近似解，即()k k y t y ≈。

如何构造“近似满足微方程”的“点集”呢？首先为这些点选择横坐标，为方便起见，将区间[,]a b 划分为M 个等距子区间，并选择网络点k t a kh =+， k=0，1，……，M 其中b a h M-=（1） 值h 称为步长。

然后近似解'(,)y f t y = 在0[,]M t t 上， 00()y t y = （2） 设()y t ，'()y t 和''()y t 连续，；；利用泰勒定理将()y t 在0t t =处展开，对每个值t ，存在一个0t 和t 之间的值1c ，使得 ''2'10000()()()()()()2y c t t y t y t y t t t -=+-+ （3） 将'00()(,())y t f t y t =和10h t t =-代人等式（3），得到1()y t 的表示：''211000()()()(,())2y c h y t y t hf t y t =++ （4） 如果步长 h 足够小，则可以忽略 2 次项（包含2h 的项），得到1000(,)y y hf t y =+ （5） 这就是欧拉近似。

重复该过程，就能得到近似解曲线()y y t =的一个点序列。

欧拉方法的一般步骤是 1k k t t h +=+ ， 1(,)k k k k y y hf t y +=+ 其中 k = 0，1，……，M-1[8]（6）2.4.2 泰勒级数法泰勒级数法有着广泛的应用，并且是比较求解初值问题的各种不同数值方法的标准，它可设计为任意指定的精度。

下面首先将泰勒定理用新的公式表示，使之适合于求解微分方程。

定理9.5（泰勒定理）设1()N y t C+∈ 0[,]t b ，且()y t 在不动点0[,]k t t t b =∈处有N 次泰勒级数展开：1()()(,())()N k k N k k y t h y t hT t y t O h ++=++ （1） 其中，()11()(,())!j Nj k N k k j y t T t y t h j -==∑ （2） ()1()(,())j j y t f t y t -=表示函数f 关t 的（1j -）次全导数。

求导公式可以递归地计算： ''''(3)''''22(4)''2'''''''''3232()()()2()2()33()33()(33)(2)t y t y tt yt y yy tt yt yy y t y ttt ytt yyt ty y yy yyy ttt ytt yyt yyy y tt yt yy y t fy t f f y f f fy t f f y f y f y f f f f f f f f f y f f y f y f y f y f y y f y f f f f f f f f f f f f f ==+=+=+++=++++=++++++=++++++23()()()t y yt yy y t y f f f f f f f f f f +++++ （3）并且一般有()(1)()(,())N N yt P f t y t -= （4） 其中P 为导数算子()P f t y∂∂=+∂∂区间0[,]M t t 上的初值问题'()(,)y t f t y =的近似数值解可由各子区间1[,]k k t t +上的公式（1）来推导。

N 次泰勒方法的一般步骤为 323211...2!3!!N N k k d h d h d h y y d h N +=+++++ （5） 其中在各步0,1,......,1k M =-有()(),1,2,......,j j k d y t j N ==。

N 次泰勒方法的最终全局误差是1()N O h +阶的，因此可选择所需大小的N ，使得误差足够小。

如果N 是固定，则理论上可以推导出步长h ，使之满足任意的最终全局误差。

然而在实际运算中，通常用h 和/2h 计算两个近似结果集，然后比较其结果[9]。

2.4.3 龙格—库塔方法泰勒方法的优点是最终全局误差的阶为()NO h ，并且可以通过选择较大的 N 来得到较小的误差。

然而泰勒方法的缺点是，需要先确定 N ，并且要计算高阶导数，它们可能十分复杂。

每个龙格一库塔（Runge-Kutta ）方法都由一个合适的泰勒方法推导而来，使得其最终全局误差为()N O h 。

一种折中方法是每步进行若干次函数求值，从而省去高阶导数计算。

这种方法可构造任意 N 阶精度的近似公式。

最常用的是N= 4 的龙格一库塔方法，它适用于一般的应用，因为它非常精确、稳定，且易于编程。

许多专家声称，没有必要使用更高阶的方法，因为提高的精度与增加的计算量相抵消。

如果需要更高的精度，则应该使用更小的步长或某种自适应方法。