matlab求解微分方程组dx/dt=-0.008u-0.018v-0.001wdy/dt=-0.07u-0.09v-0.007wdz/dt=-1.11u-1.07v-0.03wdu/dt=-0.03x-0.06y-0.1zdv/dt=-0.1x-0.47y-0.01zdw/dt=-1.7x-1.6y-0.03z初值条件:t=0时,x=17,y=195,z=1500,u=4,v=62,w=200最好能贴上代码,

Question

matlab求解微分方程组
dx/dt=-0.008u-0.018v-0.001w
dy/dt=-0.07u-0.09v-0.007w
dz/dt=-1.11u-1.07v-0.03w
du/dt=-0.03x-0.06y-0.1z
dv/dt=-0.1x-0.47y-0.01z
dw/dt=-1.7x-1.6y-0.03z
初值条件:t=0时,x=17,y=195,z=1500,u=4,v=62,w=200
最好能贴上代码,

Answer 1

function hh
X0=[17,195,1500,4,62,200];
TSPAN=[0 4];
[t,X] = ode45(@myfun,TSPAN,X0)
plot(t,X)
legend('x','y','z','u','v','w','Location','NorthWest')
xlabel('t'),ylabel('X')
function dX=myfun(t,X)
x=X(1);y=X(2);z=X(3);u=X(4);v=X(5);w=X(6);
dxdt=-0.008*u-0.018*v-0.001*w;
dydt=-0.07*u-0.09*v-0.007*w ;
dzdt=-1.11*u-1.07*v-0.03*w ;
dudt=-0.03*x-0.06*y-0.1*z ;
dvdt=-0.1*x-0.47*y-0.01*z;
dwdt=-1.7*x-1.6*y-0.03*z;
dX=[dxdt;dydt;dzdt;dudt;dvdt;dwdt];
图形见：
http://hi.baidu.com/chemical%5Fengineering/album/item/a2f0be0988baeb016a60fb8a.html