第8章 常微分方程初值问题的数值解法.pptVIP

* 第8章 常微分方程初值问题的数值解法 ? 考虑一阶常微分方程的初值问题 : ? 要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 ?把[a, b]n等分,节点间距 为步长，通常采用等距节点，即取 hi = h (常数)。 §1 几种简单的数值解法 ? 欧拉(Euler)方法 (1) 用向前差商近似导数 记为 向前欧拉方法的计算公式： Euler方法的几何意义 初值问题的解曲线 y=y(x) :过点 P0 (x0, y0)，从 x0??出发，以 f(x0 , y0)为斜率作一直线，与直线 x=x1?交点于 P1 (x1, y1) ，显然有 ??????????????????? ，再从 P1??出发，以 f(x1 , y1) 为斜率作直线,与直线 x=x2?交点于 P2 (x2, y2) ，其余类推，这样得到解曲线的一条近似曲线，它就是折线 P0 P1 P2?... 由于欧拉方法是用一条折线近似地代替解 y(x)，所以欧拉方法也叫折线法 例1： 用Euler方法解方程 取 h =0.1，在区间 [0, 0.5]上计算 解：由于 由Euler计算公式得： 计算结果见下表 1.056100 0.4 1.090490 0.5 1.0