常微分方程初值问题的Runge-Kutta解法开题报告.docVIP

开题报告 常微分方程初值问题的Runge-Kutta解法 一、选题的背景、意义 常微分方程是指只有一个自变量的微分方程。而且常微分方程在很多学科领域内有着重要的作用，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等等，这些问题都可以化为求微分方程的解，或者化为研究解的性质的问题。。现代科学、技术、工程中的大量数学模型都可以用常微分方程的初值问题来描述，所以研究常微分方程初值问题的一些性质、解法很有必要。 1.1 历史背景 常微分方程在微积分概念出现后即已出现。从莱布尼兹专门研究用变量变换解决一阶微分方程的求解问题的“求通解”时代，到1841年刘维尔的里卡蒂方程不存在一般初等解和柯西的初值问题的提出，常微分方程转向“求定解”时代。再到19世纪末天体力学的太阳系稳定性问题眼界需要常微分方程解的大范围性态，从而发展到了“求所有解问题”。最后到了20世纪六七十年代，常微分方程因为计算机的发展进入“求特殊解”阶段，发现了具有新性质的特殊的解和方程，如混沌（解）、奇异吸引子及孤立子等等。 常微分方程初值问题的数值求解是求近似解的一种经典方法。其中常微分初值问题的数值求解问题的方法又包括单步法和线性多步法。单步法主要有欧拉法和Runge-Kutta法。Runge-Kutta法是最重要也最便于应用的单步法。 1.2 国内外研究现状 常微分方程在科学和工程中常用于建立数学模型，通常它们没有解析解，而需要数值方法来近似求解。在科学的计算机化进程中，常微分方程也在不断地精进中，更在不同的领域得到了前所未有的发展和应用。求解常微分方程的初值问题的数值方法有单步法和多步法。其中单步法中的Runge-Kutta方法是目前工程上求常微分方程数值解的基本方法。 现在的运用非常广泛，国内外都对其研究非常多，以仅仅是运用在一些简单的方程和实际问题中，现在科学家都希望将其能够运用在更广泛和更复杂的情况下，且能构造满足各种不同目的需要。，例如：它计算函数值的次数较多。尤其是在级数高于4级时，这一缺点特别明显：计算函数值的工作量增加较快，但是精度即阶数提高较慢。还有各阶Runge-Kutta方法的绝对稳定区域不够大，对于坏条件的初值问题不尽适用等。但是它作为一类便于应用的单步方法，其实用价值还是很高的，特别是在现如今计算机发展飞快的时代。一定会有更先进的方法来使得Runge-Kutta方法得到优化，从而得到更好也更为实用的数值方法。 1.4计算前提 对于一阶的常微分方程初值问题 ， （1） 的数值方法，其中是和的已知函数，为给定的初值。为使上面的解存在、唯一且连续依赖初值，则必须关于满足Lipschitz条件：即存在常数L，使得对于任意，均成立不等式 单步显式公式的一般形式为 ， （2） 称为增量函数。一般来说微分方程问题（1）的精确解不满足式（2），即一般地有 定义1 称 为单步显式公式（2）在处的截断误差。 一个求解公式的局部截断误差刻画了其逼近微分方程的准确程度。根据上述定义可以直接求的各式的局部截断误差。例如欧拉公式的局部截断误差为 定义2 如果一个求解公式的局部截断误差为，则称该求解公式是阶的，或称具有阶精度。 二、研究的基本内容与拟解决的主要问题 2.1 基本内容 （1）欧拉法 单步法中最简单的数值解法是欧拉法。欧拉法的基本思想是由推导出下一个节点，其推导公式为 ，. 其局部截断误差为。欧拉法的几何意义是用一条过的折线来近似替代过的积分曲线。 虽然这个方法的计算很简单，但精确度较低，所以有了精确度较高的梯形方法。推导公式为，，它的局部截断误差为。此为隐型格式，迭代求解所需计算量大。 结合上述两种方法的优点，从而得到改进欧拉方法。其计算公式为， 。 （2）线性多步法 常用的多步法主要有Adams法和Milne法，其中Adams法又包括显式Adams法和隐式Adams法。显式Adams法的线性步公式为 其中 ， 。 且隐式Adams公式为 其中