数值分析数值计算小论文.docxVIP

Word可编辑 Runge-Kutta法的历史开展与应用 摘要Runge-Kutta法是极其重要的常微分方程数值解法，本文仅就其起源及开展脉络加以简要研究。对Runge、Heun以及Kutta等人的奉献做出适当评述，指出Runge-Kutta 方法起源于Euler折线法。同时对Runge-Kutta法的应用做简要研究。 关键词 Euler折线法 标准四阶Runge-Kutta法 应用 一、开展历史 Euler折线法 在微分方程研究之初，瑞士数学家(1707.4—1783.9)做出了开创性的工作。他和其他一些数学家在解决力学、物理学问题的过程中创立了微分方程这门学科。在常微分方程方面，Euler 在1743年发表的论文中，用代换给出了任意阶常系数线性微分方程的古典解法，最早引入了“通解〞和“特解〞的概念。 1768年，Euler在其有关月球运行理论的著作中，创立了广泛用于求初值问题 的数值解的方法，次年又把它推广到二阶方程。欧拉的想法如下：我们选择，然后在情况下用解函数的切线 代替解函数。这样对于点 就得到 。 在重复如上的程序再次计算新的方向就会得到所谓的递推公式： 这就是Euler 方法。通过连接所有这些切线得到的函数被称为Euler 折线。如果我们令， 这些折线就会越来越接近解函数。 Euler 折线法是最早出现的，虽然它亦是常微分方程初值问题的最简单的数值解法， 但它的一些特性和研究方法对于更复杂的方法却具有普遍意义。几十年后，法国数学家A .L .Cauchy(—1857.5)在历史上首次研究了常微分方程的局部性态。对于给定的初值问题()和()，在连续可微的假设下， 他用类似于欧拉折线的方法构造逼近解， 利用微分中值定理估计逼近解之间差的上界，严格证明了以为中心的一个小区间上逼近解收敛， 其极限函数即为所提问题的解。同时Cauchy 指出，这种方法也适用于常微分方程组，所以欧拉方法有时又称Euler-Cauchy 折线法。 2.2 Runge-Kutta方法 德国数学家(1856—1927)是数值方法开展史上具有里程碑作用的人物。1895年，他在Hanover发表了关于微分方程数值解法的经典论文?常微分方程数值解法?， 此文成为常微分方程Runge-Kutta方法的发端。此后，Runge结合教学活动积极投身于开展一般的数值分析特别是各种实际应用中的Runge-Kutta方法〔严格来说，此方法在Kutta做出工作后才能称作Runge-Kutta方法〕。 Runge伟大创造的思想是什么呢？他的灵感来自于初值问题〔1.1〕和〔1.2〕与积分问题 〔〕 之间的比照，显然，等式〔1,3〕右侧数值积分的精确度决定的精确度，Runge发现， Euler方法采用的是左矩形公式 ， 即用高为宽为的矩形代替数值积分， 而这个公式的精确度并不高。因此他说：最好通过插入上述Euler步骤的结果来代替未出现的值， 把精度更高的中点法那么和梯形法那么拓展到微分方程。 ， 。 其中和分别表示用中点法和梯形法算得的数值积分。与他的后继者一样，Runge 用Taylor 展开式说明上述两方法的局部误差是，方法的阶为2。不过他的梦想却是使用具有更高精度的Simpson法那么。但是众所周知，的微小变化往往易产生假象，令人误以为可以获得更高的阶。Taylor 级数展开说明，如果依赖于，事实上这个表达形式只是2 阶的。接着Runge发现，通过重复使用Euler步骤对梯形法那么做些许调整，会使形式成为3 阶方法。Runge还把他的方法及方法的展开式拓展到微分方程组。 1900 年，Runge的同胞K .Heun评论说，Runge获得的上述方法是归纳性的而且是令人费解的，他提出采用“更具一般性〞的Gauss方法。于是一般的Gauss求积公式 ， 被扩充为 〔〕 把〔〕式的右端进行二元Taylor展开后与的Taylor 展开式的对应项的系数比较，适中选取参数使方法具有尽可能高的精度。 Runge的另一个同胞(1867—1944)，1894到1909年在Munich做助教，在那里受到Runge文章的吸引并在Heun论文的鼓励下发表了他的研究结果。他认为，为什么不让已经求得的导数值全部进入到新的求值点的计算中呢？基于这样的想法，Heun格式就被Kutta代替为如下格式： (其中称为级)这个格式在满足所需的阶条件上能够允许更多的自由度。在古典的Runge-Kutta方法中，对系数的选择极大地取决于由这些系数构成的方法是否方便进行桌上计算，而所谓的古典方法是指在前计算机时代得到的方法。而这些方法对于在自动计算机上使用那么未必是最适合的方法