常微分方程实验汇报.docVIP

《常微分方程》实验报告 班 级 学 号 姓名 课程编号任课教师 评分 特别注意： 上交作业只能为1个压缩文件. 请把实验报告（.doc，word文件）和源文件（.m，Matlab程序文件）放在一个文件夹下后压缩，文件夹及压缩文件命名规则：学号+姓名。例如：201202000526刘大壮.rar。 每一M文件只能用于求解一个小题，m文件命名格式：ex1_1.m,ex2.m. 实验报告当场做完后进行提交，作为最终成绩的评判材料。 实验目的及要求 实验要求学生熟练掌握Matlab在常微分方程中的应用，主要涉及以下四个方面： 1、掌握Matlab中用于求常微分方程解析解的方法； 2、了解微分方程数值解的求解方法，向量场、等高线等图形的做法； 3、熟悉矩阵、行列式与代数方程组的求解方法； 4、熟悉拉普拉斯变换方法求解初值问题的过程。 二、用Matlab完成以下作业 1. 求下列微分方程（组）的解析解，并验证解是否正确： （1） （2） （3） ２．求解微分方程 先求解析解，再求数值解，并作图比较。 3. 求解方程 首先看是否存在解析解，若不存在，求其数值解并作图.不妨取. 4．求下面方程组的奇点，并通过变换将奇点变为原点，进一步判断奇点的类型及稳定性。 5. 做V函数的等高线图。 6. 做方程组 的轨线图及向量场。 7. 模拟laplace变换方法求解方程组初值问题： 其中。 三、问题及建议 实事求是的分析你在实验过程或常微分方程课程学习中遇到的问题，给出你的建议。 实验报告 附录：Matlab在常微分方程中常用命令 一． 解析解： 1． dsolve(equ1,equ2,…)：Matlab 求微分方程的解析解．equ1、equ2、…为方程（或条件）．写方程（或条件）时用 Dy 表示y 关于自变量的一阶导数，用用 D2y 表示 y 关于自变量的二阶导数，依此类推． 2． simplify(s)：对表达式 s 使用 maple 的化简规则进行化简． 二．数值解： [T,Y] = solver(odefun,tspan,y0) 求微分方程的数值解．说明： 1． 其中的 solver为命令 ode45、ode23、ode113、ode15s、ode23s、ode23t、ode23tb 之一． 2. odefun 是显式常微分方程： 3. 在积分区间 tspan=上，从到，用初始条件求解． 4．要获得问题在其他指定时间点上的解，则令 tspan= （要求是单调的）． 5． 因为没有一种算法可以有效地解决所有的 ODE 问题，为此，Matlab 提供了多种求解器 Solver，对于不同的ODE 问题，采用不同的Solver． 求解器 Solver ODE类型 特点 说明 ode45 非刚性 单步算法；4、5阶Runge-Kutta方程；累计截断误差达 大部分场合的首选算法 ode23 非刚性 单步算法；2、3阶Runge-Kutta方程；累计截断误差达 使用于精度较低的情形 ode113 非刚性 多步法；Adams算法；高低精度均可到 计算时间比 ode45 短 ode23t 适度刚性 采用梯形算法 适度刚性情形 ode15s 刚性 多步法；Gears反向数值微分；精度中等 若 ode45 失效时，可尝试使用 ode23s 刚性 单步法；2阶 Rosebrock 算法；低精度 当精度较低时，计算时间比 ode15s 短 ode23tb 刚性 梯形算法；低精度 当精度较低时，计算时间比 ode15s 短 6． 要特别的是：ode23、ode45 是极其常用的用来求解非刚性的标准形式的一阶常微分方程(组)的初值问题的解的 Matlab 的常用程序，其中： ode23 采用龙格-库塔2 阶算法，用3 阶公式作误差估计来调节步长，具有低等的精度． ode45 则采用龙格-库塔4 阶算法，用5 阶公式作误差估计来调节步长，具有中等的精度． 三．相关作图函数： 1. 向量场： （1） quiver(x,y,u,v)　该函数使用箭头来直观的显示矢量场表示通过在(x, y)指定的位置绘制小箭头来表示以该点为起点的向量(u,v) （2） quiver3(x,y,z,u,v,w) 该函数在(x,y,z)处显示(u,v,w)向量。该函数用于绘制等高线 （2）contour3(x,y,z,[a,b]) 绘制z在[a,b]范围内的等高线立体图。 四．辅助计算（矩阵、行列式与代数方程组解）： 1. 指数函数： （1） exp(A)　 （2） exp2(A)　 （3） exp3(A)