常微分方程的求解 实验六.doc

《数学实验》报告 实验名称 常微分方程求解 ??? ???? ? 学??号?????? ???? ?? 2013年5月【实验目的】 学习在MATLAB中如何求解微分方程的方法； 基本的微分求解，学过的基础知识求解方程； 运用的解法即数值解法解微分方程 注意不同方法下求得微分方程的优缺点。 求解微分方程为 用数值方法求解下列方程，用不同颜色和线将画在同一个图形窗口里 初始时间：：条件： 三、 1. y=dsolve(Dy=x*sinx/cosy,x) 2. 定义的程序： function xdot=exf(t,x) xdot=[0 1;1 -t]*x+[0;1]*(1-2*t); ： clc clf [t,x]=ode23(exf,[0,pi],[0.1;0.2]); y=x(:,1) dy=x(:,2) plot(t,y,r,t,dy,b) 1. 2.先将高阶微分方程转化为一阶微分方程。令即原方程为： 写成矩阵形式为： 放入函数，命令如下：[t,x]=ode23(exf,[0,pi],[0.1;0.2])[t,x]中求出的按列排列，故用ode23即为y。 表达式的右端写成exf.m文件： function xdot=exf(t,x) xdot=[0 1;1 -t]*x+[0;1]*(1-2*t); 如下： clc clf [t,x]=ode23(exf,[0,pi],[0.1;0.2]); y=x(:,1) dy=x(:,2) plot(t,y,r,t,dy,b) 结果如图： 对于高阶常微分方程，需要先将它转化为一阶常微分方程组，即状态方程所以求解常微分方程需要先用学过的知识对方程进行处理； 运用数值解法求解微分方程时需要建立模型的函数文件，然后运用函数调用的知识中调用子程序解方程的解； 方程有数值积分解也有解析解两者各自优缺点，不同的求解方法不同的场合，具体方程可以将两者互相对比画出对比的； 常微分方程的步骤改写为矩阵形式函数文件微分方程图形显示。 6