《常微分方程的求解与定性分析》由会员分享，可在线阅读，更多相关《常微分方程的求解与定性分析（17页珍藏版）》请在金锄头文库上搜索。
1、学生实验报告（4）学生姓名陈宋学号008888同组人：无实验项目Matlab程序设计与作图必修 选修 演示性实验 验证性实验 操作性实验 综合性实验实验地点H110实验仪器台号 D02 指导教师谢建实验日期及节次2011-04-18第八九十节一、实验综述. 归纳和学习求解常微分方程(组)的基本原理和方法； 2. 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； 3. 熟悉MATLAB软件关于微分方程求解的各种命令； 4. 通过范例学习建立微分方程方面的数学模型以及求解全过程； 通过该实验的学习，使学生掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方
2、程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是十分必要的。二、实验过程（实验步骤、记录、数据、分析）1开启MATLAB软件平台，开启MATLAB编辑窗口； 2根据问题，建立的线性规划模型，并编写求解规划模型的M文件； 3保存文件并运行； 4观察运行结果(数值或图形)，并不断地改变参数设置观察运行结果； 5根据观察到的结果和体会，写出实验报告。三、实验要求与任务 根据实验内容和步骤，完成以下实验，要求写出实验报告（实验目的问题数学模型
3、算法与编程计算结果分析、检验和结论） 1求微分方程的解析解，并画出它们的图形。 y = y + 2 x, y (0) = 1, 0 x eps*x(i)x(i+1)=x(i)-myfun6(x(i)*(x(i)-x(i-1)/(myfun6(x(i)-myfun6(x(i-1);b=x(i+1)-x(i);i=i+1;if(in) error(n is full);endenddisp(i-2);y=x(i);运行结果如下：（2）建立模型如下：代入初始条件可以得出x=-10*200cos*exp(-0.1t)+10*200cos建立myfun7函数如下：function y=myfun7(x)
4、y=200*sin(x)*(-10*log(1-360/2000/cos(x)-4.9*(-10*log(1-360/2000/cos(x).2)-160建立newton2函数如下：function y=newton2(x0,x1,n,tol)x(1)=x0;x(2)=x1;b=1;i=2;while(abs(b)eps*x(i)x(i+1)=x(i)-myfun7(x(i)*(x(i)-x(i-1)/(myfun7(x(i)-myfun7(x(i-1);b=x(i+1)-x(i);i=i+1;if(in) error(n is full);endenddisp(i-2);y=x(i);运行结果如下：错后进行调制、再进行运行得到结果，结束。心得体会：通过自己的一次次锻炼，总算对matlab的线性规划有所了解了，虽然操作中有很多错误，尤其是模型没有建好所带来的错误，在一次次的修改调试中不断提升自己的能力，希望自己进一步加强学习Matlab。四、指导教师评语及成绩：成绩： 指导教师签名： 批阅日期： 友情提示：方案范本是经验性极强的领域，本范文无法思考和涵盖全面，供参考!最好找专业人士起草或审核后使用。
