问题描述：
用matlab求解常系数微分方程组并画出图像
哪位大神帮我用matlab解以下微分方程组并作出图像
麻烦给出语句
dx1/dt=x1*(1-x1/150000-0.5*x2/30000)
dx2/dt=x2*(-1-x2/30000+0.5*x1/150000-2*x3/2200)
dx3/dt=x3*(-1-x3/2200+2*x2/30000)
x1(0)=121000
x2(0)=27000
x3(0)=2000
哪位大神帮我用matlab解以下微分方程组并作出图像
麻烦给出语句
dx1/dt=x1*(1-x1/150000-0.5*x2/30000)
dx2/dt=x2*(-1-x2/30000+0.5*x1/150000-2*x3/2200)
dx3/dt=x3*(-1-x3/2200+2*x2/30000)
x1(0)=121000
x2(0)=27000
x3(0)=2000
问题解答：
我来补答
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 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结果
展开全文阅读