微分方程求解的后退欧拉法、龙格库塔法.docVIP

数值计算方法实验上机测验 微分方程求解的后退欧拉法、龙格库塔法（三阶、四阶） 日期： 2011-06-16 一、测验目的 1. 学习matlab的使用方法。 2. 掌握常微分方程的几种数值解法：后退欧拉法、龙格库塔法（三阶、四阶）。 3. 比较各方法的数值解及误差，了解各方法的优缺点。 二、实验题目 给定的初值问题 及其精确解 按 (1)后退欧拉法，步长h=0.05, h=0.1; (2)三阶龙格—库塔法，步长h=0.05, h=0.1; (3)四阶标准龙格—库塔法，步长h=0.05，h=0.1; 求在节点处的数值解及误差比较各方法的优缺点。 三、实验原理 1.对于后退欧拉法： 利用 进行迭代求解可以完成计算 需要将微分方程表达式和精度计算表达式作为两个函数保存在m文件里并在程序中调用： ①微分方程（wangdaren_wf） function z=wangdaren(x,y) z=-2*x*(y*exp(x^2)-1)/exp(x^2) end ②精确解计算(wangdaren _jq) function z= wangdaren _jq(x) z=x*x/exp(x^2) end 龙格—库塔法基本思想： 用在几个不同点的数值加权平均斜率，使截断误差的阶数尽可能高。即取不同点的斜率的加权平均作为平均斜率，以便提高阶数。 2.对于三阶龙格—库塔法： 利用 可以完成计算 3.对于四阶龙格—库塔法： 利用 可以完成计算 四、实验内容 由上述实验原理叙述的后退欧拉法，三阶龙格—库塔法，四阶龙格—库塔法几种常微分方程数值解法分别对已给定的初值问题进行求解，比较各方法的数值解及误差，了解各方法的优缺点。 五、实验结果 1. 对于后退欧拉法： ①若h=0.05 h=0.05; y=0; x=0; for i=1:20 k1=h*wangdaren_wf(x,y); k2=h*wangdaren_wf(x+h,y+k1); y=y+0.5*k1+0.5*k2 x=x+h; z=wangdaren_jq(x) t=abs(y-z); A=[x y z t] End 结果：A = 1.0000 0.3671 0.3679 0.0007 自变量(x) 近似值(y) 精确值(z) 误差绝对值(t) 0.0500 0.0025 0.0025 0.0000 0.1000 0.0099 0.0099 0.0000 0.1500 0.0220 0.0220 0.0000 0.2000 0.0384 0.0384 0.0000 0.2500 0.0586 0.0587 0.0001 0.3000 0.0822 0.0823 0.0001 0.3500 0.1082 0.1084 0.0001 0.4000 0.1362 0.1363 0.0002 0.4500 0.1652 0.1654 0.0002 0.5000 0.1944 0.1947 0.0003 0.5500 0.2232 0.2235 0.0003 0.6000 0.2508 0.2512 0.0004 0.6500 0.2765 0.2769 0.0005 0.7000 0.2997 0.3002 0.0005 0.7500 0.3199 0.3205 0.0006 0.8000 0.3368 0.3375 0.0006 0.8500 0.3501 0.3508 0.0007 0.9000 0.3596 0.3603 0.0007 0.9500 0.3653 0.3660 0.0007 1.0000 0.3671 0.3679 0.0007 ②若h=0.1 h=0.1; y=0; x=0; for i=1:10 k1=h*wangdaren_wf(x,y); k2=h*wangdaren_wf(x+h,y+k1); y=y+0.5*k1+0.5*k2 x=x+h; z=wangdaren_jq(x); t=abs(y-z); A=[x y z t] end 结果为A= 1.0000 0.3648 0.3679 0.0031 自变量 近似值 精确值 误差 0.1000 0.0099 0.0099 0.0000 0.2000 0.0383 0.0384 0.0001 0.3000 0.0819 0.0823 0.0003 0.4000 0.1357 0.1363 0.0006 0.5000 0.1937 0.1947 0.0010 0.6000 0.2496 0.2512 0.0015 0.7000 0.2982 0.300