微分方程和差分方程作业题专业:土规1101班姓名:刘迈克学号:2011306200521微分方程模型作业:1.用matlab求解微分方程组00dxxydtdyxydt(1)求在初始条件0|,1|00ttyx下的特解,并画出解函数()yfx的图形.(2)分别用ode23、ode45求此微分方程组初值问题的数值解(近似解),求解区间为[0,2]t.利用画图来比较两种求解器之间的差异.解:程序:[x,y]=dsolve('Dx+x+y=0','Dy+x-y=0',
'x(0)=1','y(0)=0','t')ezplot(x,y,[0,5]);(1)xtttteeee222242422121424222tteey(2)先编写函数文件verderpol
mfunctionxprime=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升的速率注入清水,容器内装有搅拌器能将溶液迅时搅拌均匀,并同时以每分钟2升的速率放出盐水,求1小时后容器里的盐水中还含有多少净盐
解:分析:设t时刻(单位为分钟)净盐量为x(t)(单位:kg/l),考虑时间区间为[t,t+△