第10章-常微分方程数值解 2.pptVIP

　　若将(4.2)展开式写成 则 称为局部截断误差主项. 以上定义对隐式单步法(4.1)也是适用的.例如，对隐式欧拉法(3.3)其局部截断误差为 这里p=1是1阶方法，局部截断误差主项为 　　同样对梯形法(3.5)有 所以梯形方法(3.5)是二阶的. 其局部截断误差主项为 * 分别将 y(xn+1), y(xn-1) 在 xn 点处用泰勒公式展开： 二步欧拉法的局部截断误差： 二步欧拉法： 2 阶精度 * 各种欧拉法的比较 方法 精度 评述 显式欧拉法 1 最简单，精度低 隐式欧拉法 1 不便计算，稳定性好 二步欧拉法 2 需要两步初值，且第 2 个初值只能由其它方法给出，可能对后面的递推精度有影响 梯形公式 欧 拉 法 2 精度有所提高，但为隐式，需要迭代求解，计算量大 * 提高计算精度的方法 增加展开项数 增加y(xn+1)函数的泰勒展开式的阶数，即逐次取其高次项，这样就可以减少截断误差。其误差系列为 O(h2), O(h3), O(h4), … 缩小步长值h 由于差分间隔减小，即h值变得很小，可使计算精度提高。在理论上步长值h越小越好，而且h趋近于0最好。但是，h越小会使计算次数增加，浪费计算时间，而且还会使舍入误差增大 * 10.4.2单步法的收敛性 显式单步法可统一写成： 增量函数，仅依赖于函数 f，且仅仅是 xn, yn, h 的函数 求 y = y(x) 求 y(xn)，xn = x0 + nh 离散化 求 yn ? y(xn) 某种数值方法 h ? 0 时，近似解 是否收敛到精确解 ，它应当 是一个固定节点，因 此 h ? 0 时应同时附 带 n ? ? * 对于 p 阶的常微分方程数值算法，当 h ? 0， n ? ? 时，是否 yn+1 ? y(xn+1)？ p 阶算法的局部截断误差为： 显然： 局部截断误差的前提假设是： 局部截断误差 ? 0 并不能保证算法收敛 * 定义：若求解某初值问题的单步数值法，对于固定的 当 h ? 0 且 n ? ? 时，它的近似 解趋向于精确解 y(xn)，即： 则称该单步法是收敛的。 定义：称 y(xn) - yn 为单步法的近似解 yn 的整体截断 误差。 单步法收敛 h ? 0 且 n ? ? 时，yn 的整体截断误差 ? 0 * 收敛性定理 若某单步法满足以上条件，则该方法是收敛的 则该单步法的整体截断误差为： 　　若单步法 具有 p 阶精度，且增量函数 关于 y 满足： Lipschitz 条件： 初值 y0 是准确的 计 算 方 法 课 件 1、常微分方程与解 为n阶常微分方程。 如果函数 在区间[a,b]内n阶可导，称方程 满足方程的函数 称为微分方程的解。 则 如 为任意常数） 一般称为方程的通解。 为方程的解。 如果 则有 为方程满足定解条件的解。 第10章 常微分方程的数值解 10.1引 言 科学研究和工程技术中的问题往往归结为求某个常微分方程的定解问题. 常微分方程的理论指出，除少数简单情况能获得初值问题的初等解（用初等函数表示的解）外，绝大多数情况下是求不出初等解的.有些初值问题即便有初等解，也往往由于形式过于复杂而不便处理。 常微分方程的数值解法常用来求近似解，由于它提供的算法能通过计算机便捷地实现，因此近年来得到迅速的发展和广泛的应用。 10.2 初值问题解法的基本概念 科学技术中常常需要求解常微分方程的定解问题. 这类问题最简单的形式，是本章将着重考察的一阶方程的初值问题 我们知道，只有f(x, y)适当光滑—譬如关于y满足利普希茨(Lipschitz)条件 理论上就可以保证初值问题的解y＝f(x)存在并且唯一.我们以下的讨论，都在满足上述条件下进行。 虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，实际问题中归结出来的微分方程主要靠数值解法. 所谓数值解法, 就是寻求解y(x)在一系列离散节点 上的近似值 y1,y2,?,yn,yn+1,?. 相邻两个节点的间距hn=xn+1-xn称为步长. 今后如不特别说明，总是假定 hi=h(i=1,2,?)为定数, 这时节点为xn=x0+nh(i=0,1,2,?) (等距节点). 常微分方程数值解是一组离散的函数值数