10级常微分实习.docVIP

福建农林大学计算机与信息学院 （数学类课程） 课程实习报告 课程名称： 常微分方程课程实习 实习题目： 常微分方程数值求解问题的实习 姓 名： 邱育华 系： 信息与计算科学 专 业： 信息与计算科学 年 级： 2010 学 号： 102260003049 指导教师： 陈永雪 职 称： 讲师 2011 年 11 月 19 日  福建农林大学计算机与信息学院数学类 课程实习报告结果评定 评语： 成绩： 指导教师签字： 评定日期： 目 录 1、实习的目的和任务…………………………………………4 2、实习要求……………………………………………………4 3、实习地点……………………………………………………4 4、实验主要仪器………………………………………………4 5、实验主要内容………………………………………………4 5.1、用不同格式对一个初值问题的数值求解极其分析……………4 5.1.1、求精确解………………………………………………………4 5.1.2、用欧拉法求解…………………………………………………5 5.1.3、用4级4阶龙格库塔法求解…………………………………8 5.1.4、问题的讨论与分析……………………………………………10 5.2、一个算法不同步长求解同一个初值问题及分析………………16 6、数值求解Lorenz方程………………………………………18 7、结束语……………………………………………………………………19 参考文献 常微分方程课程实习 实习的目的和任务 目的：通过课程实习能够应用MATLAB软来计算微分方程（组）的数值解；了解常微分方程数值解。 任务：通过具体的问题，利用MATLAB软件来计算问题的结果，分析问题的结论。 实习要求 能够从案例的自然语言描述中，抽象出其中的数学模型；能够熟练应用所学的数值解计算方法；能够熟练使用MATLAB软件；对常微分方程数值解有所认识，包括对不同算法有所认识和对步长有所认识。 实习地点 数学实验室、学生宿舍 主要仪器设备 计算机、 Microsoft Windows 7 Matlab（R2009a） 实习内容 5.1 先求精确解，然后用欧拉方法、4阶龙格—库塔方法分别求下面微分方程的初值 dy/dx=(y^2-x-2)/(4(x+1)) y(0)=2 x∈[0，2] .1求精确解 首先可以求得其精确解为：y=sqrt(x+1)+1 5.1.1 程序代码： (1) x=0:0.1:2； y1= sqrt(x+1)+1; plot(x,y,b*-); y = Columns 1 through 10 2.0000 2.0488 2.0954 2.1402 2.1832 2.2247 2.2649 2.3038 2.3416 2.3784 Columns 11 through 20 2.4142 2.4491 2.4832 2.5166 2.5492 2.5811 2.6125 2.6432 2.6733 2.7029 Column 21 2.7321 即对应的[x,y]为： 0 2 0.1000 2.0488 0.2000 2.0954 0.3000 2.1402 0.4000 2.1832 0.5000 2.2247 0.6000 2.2649 0.7000 2.3038 0.8000 2.3416 0.9000 2.3784 1.0000 2.4142 1.1000 2.4491 1.2000 2.4832 1.3000 2.5166 1.4000 2.5492 1.5000 2.5811 1.6000 2.6125 1.7000 2.6432 1.8000 2.6733 1.9000 2.7029 2.0000 2.7321 对应的图像为： 5.1.2 用欧拉法求解 程序如下： 建立函数文件yu.m function [x,y]=yu(fun,x,y0,h) x=x(1):h:x(2); y(1)=y0; for n=1:length(x)-1