微分方程习题.docxVIP

精品文档 精品文档 PAGE PAGE9 精品文档 PAGE 微分方程和差分方程作业题 专业：土规1101班姓名：刘迈克学号：2011306200521微分方程模型作业： 1．用matlab求解微分方程组 dx y 0 x dt dy y 0 x dt （1）求在初始条件x|t01,y|t0 0下的特解，并画出解函数y f(x)的图形． （2）分别用ode23、ode45求此微分方程组初值问题的数值解 (近似解)，求解 区间为t[0,2]．利用画图来比较两种求解器之间的差异． 解： 程序： [x,y]=dsolve(Dx+x+y=0,Dy+x-y=0,... x(0)=1,y(0)=0,t) ezplot(x,y,[0,5]); 1） 1 e 2t 1 2e2t 2 x 2e2t 4 4e2t 2 y 2 2e2t 4 4e2t 2） 先编写函数文件 verderpol.m functionxprime=verderpol(t,x) xprime=[-x(1)-x(2);-x(1)+x(2)]; 再编写脚本文件 vdpl.m，在命令窗口直接运行该文件 clear; y0=[1;0]; [t,x]=ode45 （23）(verderpol,[1,40],y0); plot(t,x(:,1),or-); ode45求解器微分方程组初值问题的数值解 (近似解) Ode23求解器微分方程组