数值分析实验报告之常微分方程数值解.pdfVIP

数学与计算科学学院 实 验 报 告 实验项目名称 常微分程数值解 所属课程名称 数值法 B 实 验 类 型 验证 实 验 日 期 2013.11.11 班 级 学 号 姓 名 成 绩 专业资料 . 一、实验概述： 【实验目的】 1.掌握求解常微分程的欧拉法； 2.掌握求解常微分程的预估校正法； 3.掌握求解常微分程的经典的四阶龙格库塔法； 4.能用 C 语言或 MATLAB 将上述三种算法用程序运行出来； 5.将算法实例化，并得出三种算法的相关关系，如收敛性、精度等； 6.附带书中例题的源程序见附录 1。 【实验原理】 1．欧拉格式 （1 ）显式欧拉格式： y y hf (x , y ) n 1 n n n 2 2 h h 2 局部截断误差： y( x ) y y ( ) y (x ) o (h ) n 1 n 1 n 2 2 （2 ）隐式欧拉格式： yn 1 yn hf (xn 1 , y n 1) 2 h 2 局部截断误差： y(x ) y y (x ) o(h ) n 1 n 1 n 2 2 ．预估校正法 预估： y y hf (x , y ) n 1 n n n 资料 . . h 校正： y y [ f (x , y ) f ( x , y )] n 1 n n n n 1 n 1 2 h 统一格式： y y [ f (x , y ) f (x h, y hf (x , y ))] n 1 n n n n n n n 2 y