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

《数学实验》报告 实验名称 常微分方程求解 ??? ???? ? 学??号?????? ???? ?? 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]*