数值分析实验报告之常微分方程数值解汇总.docx

数学与计算科学学院 数学与计算科学学院 实 验 报 告 实验项目名称 常微分方程数值解 所属课程名称 数值方法 B 实 验 类 型 验证 实 验 日 期 2013.11.11 班 级 学 号 姓 名 成 绩 PAGE PAGE 10 一、实验概述： 【实验目的】 掌握求解常微分方程的欧拉法； 掌握求解常微分方程的预估校正法； 掌握求解常微分方程的经典的四阶龙格库塔法； 能用 C 语言或 MATLAB 将上述三种算法用程序运行出来； 将算法实例化，并得出三种算法的相关关系，如收敛性、精度等； 附带书中例题的源程序见附录 1。 【实验原理】 欧拉格式 显式欧拉格式： yn 1 yn hf ( xn , yn ) 2 2 局部截断误差： y( xn 1 ) yn 1 h y ( ) h y ( x ) n o( h2 ) 2 2 隐式欧拉格式： yn 1 yn hf ( xn 1 , yn 1) 局部截断误差： y( x ) y h y ( x ) o(h2 ) 2n 1 n 1 n 2 2 预估校正法 预估： yn 1 yn hf (xn , yn ) h 校正： y y [ f ( x , y ) f ( x , y )] n 1 n 2 n n n 1 n 1 统一格式： y y h [ f (x , y ) f (x h, y hf ( x , y ))] n 1 n n n n n n n 2 yp yn hf ( xn , yn ), p c平均化格式： y p c yn hf (xn 1 , yp ), yn 1 1 ( y y ). 2 四阶龙格库塔方法的格式（经典格式） y y h (K 2 K  2 K K ), n 1 n 6 1 2 3 4 K 1 f (xn , yn ), K f ( x h , y h K ), 2 n n 1 2 2 K 3 f h h , y2(xn n K , y 2 2 2 K 4 f ( xn h, yn hK 3 ). 【实验环境】 硬件环境： HP Microsoft 76481-640-8834005-23929 HP Corporation Intel(R) Core(TM) I5-2400 CPU @ 3.10GHz 3.09GHz,3.16GB 的内存 软件环境： Microsoft Windows XP Professional 版 本 2002 Service Pack 3 二、实验内容： 【实验方案】 方案一： 用欧拉法，预估校正法，经典的四阶龙格库塔方法求解下列 ODE 问题： 例题：在区间【 0 ，1】上以 h=0.1 用欧拉法，预估校正法 ,经典的四阶龙格库塔法求解微分方程 dy/dx=-y+x+1 ，初值 y（0）=1 ；其精确解为 y=x+exp(-x), 且将计算结果与精确解进行比较，对三个算法的收敛性的进行分析比较。 方案二： 用欧拉法，预估校正法 , 经典的四阶龙格库塔方法 求解初值问题 dy/dx= xe x y ，初值 y（0）=1 ； 将计算结果与精确解为 1 ( x2 2 2)e  x 比较在区 间[0,1] 上分别取步长 h=0.1; 0.05 时进行计算。对三个算法的收敛性进行分析比较， 【实验过程】（实验步骤、记录、数据、分析） 注：以下图形是通过 Excel 表格处理数据得出，并未通过 MATLAB 编程序所得！ dy y x 1 1、 dx y(0) 1 由题可知精确解为： y x e x ，当 x=0 时， y(x)=0 。 h=0.1 表 1 h=0.1 时三个方法与精确值的真值表 步长 Euler 法 预估校正法 经典四阶库 精确值 0.1 1.010000 1.005000 1.004838 1.249080 0.2 1.029000 1.019025 1.018731 1.055455 0.3 1.056100 1.041218 1.040818 1.091217 0.4 1.090490 1.070802 1.070320 1.131803 0.5 1.131441 1.107076 1.106531 1.176851 0.6 1.178297 1.149404 1.148812 1.226025 0.7 1.230467 1.197211 1.196586 1.279016 0.8 1.287421 1.249975 1.249329 1.335536 0.9 1.348678 1.307228 1.306570 1.395322 1.0 1.413811 1.368541 1.367880 1.458127 图 1 h