计算机数学基础—微分方程数值解(02—14).docVIP

《计算机数学基础(2)》辅导 ―第14章 常微分方程的数值解法(2002级(春)用) 中央电大 冯 泰 第14章 常微分方程的数值解法 一、重点内容 ( 欧拉公式： 局部截断误差是O(h2). ( 改进欧拉公式： 梯形公式： y(xk+1)(yk+1=yk+[f(xk,yk)+f(xk+1,yk+1)](k=0,1,2,…,n－1) 预报－校正公式： 即 平均的形式公式： 改进欧拉法的局部截断误差是O(h3)。 ( 龙格－库塔法 二阶龙格－库塔法的局部截断误差是O(h3) 三阶龙格－库塔法的局部截断误差是O(h4) 四阶龙格(库塔法公式： 其中 (1=f(xk,yk)；(2=f(xn+h，yk+h(1)； (3=f(xk+h，yn+h(2)；(4=f(xk+h，yk+h(3) 四阶龙格－库塔法的局部截断误差是O(h5). 二、实例 例1 用欧拉法解初值问题，取步长h=0.2.计算过程保留4位小数. 解h=0.2, f(x)=－y－xy2.首先建立欧拉迭代公式 当k=0，x1=0.2时，已知x0=0,y0=1，有 y(0.2)(y1=0.2×1(4－0×1)＝0.800 0 当k＝1，x2=0.4时，已知x1=0.2, y1=0.8，有 y(0.4)(y2=0.2×0.8×(4－0.2×0.8)＝0.614 4 当k=2，x3=0.6时，已知x2=0.4,y2=0.614 4，有 y(0.6)(y3=0.2×0.614 4×(4－0.4×0.4613)＝0.800 0 例2 用欧拉预报－校正公式求解初值问题，取步长h=0.2,计算 y(0.2),y(0.4)的近似值，计算过程保留5位小数.l 解 步长h=0.2, 此时f(x,y)=－y－y2sinx. 欧拉预报－校正公式为： 有迭代公式： 当k=0，x0=1, y0=1时，x1=1.2，有 当k=1，x1=1.2, y1=0.71549时，x2=1.4，有 =0.52608 例3 用改进的欧拉法平均公式，取步长h=0.1，求解初值问题 计算过程保留4位小数. 解 首先建立迭代格式： (7分) 当k=0时，x0=0,y0=1,x1=0.1,有 (11分) 当k=1时，x1=0.1, y1=1.11, x2=0.2,有 所求y(0.1)(1.11; y(0.2)(1.242 1 例4 写出用四阶龙格－库塔法求解初值问题的计算公式，取步长h=0.2计算y(0.4)的近似值.计算过程保留4位小数. 解 此处f(x,y)=8－3y, 四阶龙格－库塔法公式为 其中 (1=f(xk,yk)；(2=f(xn+h，yk+h(1)；(3=f(xk+h，yn+h(2)；(4=f(xk+h，yk+h(3) 本例计算公式为： 其中 (1=8－3 yk；(2=5.6－2.1 yk；(3=6.32－2.37yk; (4=4.208＋1.578yk 当x0=0,y0==2, 例5 选择填空题： 1. 取步长h=0.1, 用欧拉法求解初值问题的计算公式是 答案： 解答：欧拉法的公式 此处，迭代公式为 2. 改进欧拉法的平均形式公式是( ) .(A) (B) (C) (D) 答案：(D) 解答：见改进欧拉法平均形式公式. 3. 解微分方程初值问题的方法，( )的局部截断误差为O(h3). (A) 欧拉法 (B)改进欧拉法 (C)三阶龙格－库塔法 (D) 四阶龙格－库塔法 答案：(B) 解答：改进欧拉法的局部截断误差是二阶精度，O(h3)。 4. 解常微分方程初值问题的改进欧拉法预报――校正公式是 预报值：，校正值：yk+1= ． 答案：yk+hf(xk+1,(yk+1) 解答：见预报－校正公式。 5设初值问题 把区间[0,1]10等分，用欧拉法解该初值问题的公式为 . 三、练习题 1.求解初值问题欧拉法的局部截断误差是( ); 改进欧拉法的局部截断误差是( ); 四