常微分方程数值解法48822.pptVIP

常微分方程数值解法 7.0 基本概念 1. 一阶常微分方程的初值问题 (7.0-1) 注：若f 在D = {a ? x ? b , | y |+?}内连续，且满足Lip条件： ?L ?0，使 |f (x，y1) – f (x，y2)| ? L|y1 – y2| (7.0-2) 则(7.0-1)的连续可微解y(x)在[a，b]上唯一存在。 2. 初值问题的数值解 称(7.0-1)的解y(x)在节点xi处的近似值 yi ? y(xi) 为其数值解，方法称为数值方法。 注： ① 考虑等距节点: xi = a + ih，h = (b – a)/n. ② 从初始条件y(a) = y0出发，依次逐个计算y1，y2，…，yn的值，称为步进法。 两种：单步法、多步法。 ③ 二阶常微分方程y(x) = f (x，y，y)可设为一阶常微分方程组的初值问题： 引进新的未知函数z(x) = y，则 其初始条件为： 称为一阶微分方程组的初值问题，方法类似。 7.1 初值问题数值解法的构造及其精度 7.1.1 构造方法 对于(7.0-1)可借助Taylor展开（导数法）、差商法、积分法实现离散化来构造求积公式： 1. 用差商来表示： 得差分方程： yi+1 = yi + hf (xi，yi). 即为Euler公式。 若记 ? yi+1 = yi + hf (xi+1，yi+1). 称为向后Euler法。 注：① Euler法为显式，向后Euler法为隐式——须解出yi+1. 2. 梯形公式： ? ? (7.1-4) 称(7.1-4)为梯形公式??隐式公式。 显化：预估值： 校正值： (7.1-5) 称(7.1-5)为改进的Euler公式 4. 几何意义 Euler法??折线法 向后Euler法 改进Euler法??平均斜率折线法 7.1.2 截断误差与代数精度 定义7.1-1 ① 称 ?i = y(xi) – yi 为数值解yi的（整体）截断误差。 ② 若yk = y(xk)，k = 0，1，2，…，i – 1. 由求解公式得数值解 ，则称 为yi的局部截断误差。 注：局部截断误差ei是指单步计算产生的误差，而（整体）截断误差?i则考虑到每步误差对下一步的影响。 解析解 7.2.2 Runge?Kutta方法 1. 基本思想 因为 = y(xi) + hf (?，y(?)) = y(xi) + hK? 其中K? = f (?，y(?))称为y(x)在[xi，xi+1]上的平均斜率。 若取 K1 = f (xi，y(xi)) ——Euler公式 取 K2 = f (xi+1，y(xi+1)) —— 向后Euler公式 一阶精度 取 —— 梯形公式 二阶精度 猜想：若能多预测几个点的斜率，再取其加权平均作为K?，可望得到较高精度的数值解，从而避免求f 的高阶导数。 2. R?K公式 (7.2-4) 其中Kj为y = y(x)在xi + ajh (0 ? aj ? 1)处的斜率预测值。 aj，bjs，cj为特定常数。 3. 常数的确定 确定的原则是使精度尽可能高。 以二阶为例： (7.2-5) (希望y(xi+1) – yi+1 = O(hp)的阶数p尽可能高) 首先： 另一方面： 将K2在(xi，yi)处展开。 K2 = f (xi，yi) + a2hf x(xi，yi) + b21hK1 f y(xi，yi) + O(h2). 代入(7.2-5)得： yi+1 = yi + hc1 f (xi，yi) + hc2 f (xi，yi) + h2c2[a2 f