微分方程习题.doc

微分方程和差分方程作业题 专业：土规1101班 姓名：刘迈克 学号：2011306200521 微分方程模型作业： 1．用matlab求解微分方程组 （1）求在初始条件下的特解，并画出解函数的图形． （2）分别用 ode23、ode45 求此微分方程组初值问题的数值解(近似解)，求解区间为．利用画图来比较两种求解器之间的差异． 解： 程序： [x,y]=dsolve(Dx+x+y=0,Dy+x-y=0, ... x(0)=1, y(0)=0, t) ezplot(x,y,[0,5]); （1） （2） 先编写函数文件 verderpol.m function xprime=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求解器微分方程组初值问题的数值解(近似解) 两种求解器之间的差异： 由图像可知，Ode45求解器的图像中点数比较多，更加精确。 2．设初始时容器里盛放着含净盐10千克的盐水100升，现对其以每分钟3升的速率注入清水，容