Question

用matlab求解常系数微分方程组并画出图像

dx/dt=220-0.25x-0.000015xy²
dy/dt=0.000015xy²-0.64y+0.15z
dz/dt=0.64y-0.15z-0.7z
初始值为x=68,y=109,z=82
跪求好人，求解和画图，时间t定义为0-60，需要代码，跪谢了

Answer 1

1定义函数：
function y=fun(t,x)
y=zeros(3,1);
x1=x(1);x2=x(2);x3=x(3);
y(1)=x1*(1-x1/150000-0.5*x2/30000);
y(2)=x2*(-1-x2/30000+0.5*x1/150000-2*x3/2200);
y(3)=x3*(-1-x3/2200+2*x2/30000);
2解方程并画图
[t,x]=ode23(@fun,[0,1],[121000,27000,2000]);
hold on
plot(t,x(:,1))
plot(t,x(:,2))
plot(t,x(:,3))
3结果

以上回答你满意么？

追问

我输入进去怎么说不对啊，大神求帮助啊
这个我也找到了，这个定义函数容易出错吗，是不是定义一个function，输入里面的，完了再command window输入t,x那个后面的，我一点都不懂啊，求帮助