实验五 欧拉法Matlab实验报告

北京理工大学珠海学院实验报告

ZHUHAI CAMPAUS OF BEIJING INSTITUTE OF TECHNOLOGY

班级2012电气2班学号120109021010姓名陈冲指导教师张凯成绩

实验题目（实验五）欧拉法实验地点及时间JD501 2014/1/2（6-7节）

一、实验目的

1.掌握用程序语言来编辑函数。

2.学会用MATLAB编写Euler.m以及TranEuler.m函数。

二、实验环境

Matlab软件

三、实验内容

1、以书中第124页题目11为例编辑程序来实现计算结果。

2、使用MATLAB进行编写：

第一步：编写Euler.m函数，代码如下

编写TranEuler.m函数，代码如下

第二步：利用上述函数编辑命令：（可见实验结果中的截图）

在此之前先建立一个名为f.m的M文件，代码如下

function z=f(x);

z=8-3y;

再编辑代码：

得到了欧拉法的结果：y（0.4）=2.47838030901267

编辑另一段命令：

得到改进欧拉法的结果：y（0.4）=2.46543714659780

在此基础上，我还编辑龙格库达的命令窗口代码，如下：

四、实验题目

用欧拉法和改进欧拉法求解初值问题'83,(0)2

y

=-=，试取步长0.2

y y y

h=计算(0.4)的近似值。

五、实验结果

六、总结

通过这次实验我掌握了将得到的解进一步精确，而且要学会比较这几种方法的精确性，显然，四阶龙格库达比改进欧拉发精确，改进欧拉发比欧拉法精确。

实验难度不大，要比较n的取值不同，产生的影响不同。

…….