Matlab实验报告五(微分方程求解Euler折线法)

                                                                                                                                                                                                                                                                                                               

数学与信息科学系实验报告

实验名称      微分方程求解      

所属课程      数学软件与实验    

实验类型      综合型实验        

          信息与计算科学    

                             

                             

                             

指导教师                         

 

第二篇:实验五 欧拉法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

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

四、实验题目

  用欧拉法和改进欧拉法求解初值问题,试取步长计算的近似值。

五、实验结果                      

  

六、总结

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

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

…….

相关推荐