常微分实习报告常微分方程数值求解问题教程分析.docVIP

（数学类课程） 课程实习报告 课程名称： 常微分方程课程实习 实习题目： 常微分方程数值求解问题的实习 姓 名： 系： 专 业： 年 级： 学 号： 指导教师： 职 称： 年 月 日  课程实习报告结果评定 评语： 成绩： 指导教师签字： 评定日期： 目 录 1. 实习的目的和任务 1 2. 实习要求 1 3. 实习地点 1 4. 主要仪器设备 1 5. 实习内容 1 5.1 用不同格式对同一个初值问题的数值求解及其分析……………………..1 5.1.1求精确解 1 5.1.2用欧拉法求解 6 5.1.3用改进欧拉法求解 9 5.1.4用4级4阶龙格—库塔法求解 12 5.1.5 问题讨论与分析 …………………………………………………………………… 15 5.2 一个算法不同不长求解同一个初值问题及其分析…………………………………..18 6. 结束语 29 参考文献 29 常微分方程课程实习 实习的目的和任务 目的：通过课程实习能够应用MATLAB软来计算微分方程（组）的数值解；了解常微分方程数值解。 任务：通过具体的问题，利用MATLAB软件来计算问题的结果，分析问题的结论。 实习要求 能够从案例的自然语言描述中，抽象出其中的数学模型；能够熟练应用所学的数值解计算方法；能够熟练使用MATLAB软件；对常微分方程数值解有所认识，包括对不同算法有所认识和对步长有所认识。 实习地点 数学实验室、学生宿舍 主要仪器设备 计算机、 Microsoft Windows 7 Matlab 7.0 实习内容 5.1 用欧拉方法，改进欧拉方法，4阶龙格—库塔方法分别求下面微分方程的初值dy/dx=y*cos(x+2) y(-2)=1 x∈[-2，0] 5.1.1求精确解 ①变量分离方程情形:形如的方程,这里分别是的连续函数.如果,我们可将方程改写成,这样,变量就”分离”开来了,两边同时积分即可:为任意常数. ②常数变易法:一阶线性微分方程,其中在考虑区 间上是的连续函数.可先解出方程的解,这是属于变量分 离方程情形,可解得:,这里是任意常数.然后将变 易为的待定函数,令,将其代入原方程可得:所以可解得, 这里是任意常数.将代入可得原方程的通解: 为任意常数. ③恰当微分方程情形:形如的一阶微分方程,这里 假设在某矩形域内是的连续函数,且具有连续的一阶偏导数. 若,则为恰当微分方程.判断为恰当微分方程后,则可用如下解法： 设是原方程的解，则，所以设， 则，所以，由此,由此可解得，所以原方程的通解为为任意常数。 首先可以求得其精确解为：y=exp(sin(x+2)) x=-2:0.1:2; y= exp(sin(x+2)) plot(x,y,r.-); Data=[x,y] y = Columns 1 through 4 1.00000000000000 1.10498683033169 1.21977855600062 1.34382524373165 Columns 5 through 8 1.47612194644573 1.61514629644208 1.75881884576699 1.90449653438673 Columns 9 through 12 2.04900865016427 2.18874191260461 2.31977682471585 2.43807150515633 Columns 13 through 16 2.53968253238078 2.62100592628670 2.67901644757271 2.71148101768216 Columns 17 through 20 2.71712300843128 2.69571859920382 2.64811384739078 2.57616043684702 Columns 21 through 24 2.48257772801500 2.37075712617031 2.24453058057755 2.10792744704554 Columns 25 through 28 1.96494288855677 1.81933699108106 1.67447782737116 1.53323499677732 Columns 29 through 32 1.39792381994500 1.2702952188