数值计算方法后退欧拉法,龙格库塔法(三阶,四阶方法).docVIP

数值计算方法实验课测验作业——微分方程求解的后退欧拉法,龙格库塔法（三阶，四阶方法） 日期： 20一、 实验目的 1. 学习matlab的使用方法。 2. 掌握常微分方程的几种数值解法：，，3. 比较各方法的数值解及误差，了解各方法的优缺点。 二、 实验题目 给定的初值问题 y 按 (1)，步长h=0., h=0.1; (2) 龙格库塔法三阶方法，步长h=; (3) 龙格库塔法四阶方法，步长h=0.1; 求在节点……10)处的数值解及误差比较各方法的优缺点。三、 实验原理 1.对于： 利用Yn+1 = Yn + 1/2*K1 + 1/2*K2① n = 1, 2, 3…… K1 = hf（Xn, Yn） ② K2 = hf（Xn + h, Yn + K1） ③ 三式可以完成计算 需要将微分方程表达式和精度计算表达式作为两个函数保存在m文件里并在程序中调用： ①微分方程（wf） function z=wf(x,y) z=-y+x+2 end ②精确解计算(j_q) function z=ouleij_q(x) z=exp(-x)+x end 2.对于龙哥库塔阶： 利用Yn+1 = Yn + /6(K1 + 4K2 + K3 ) ①