用matlab求解常系数微分方程组并画出图像
来源:学生作业帮 编辑:百度作业网作业帮 分类:数学作业 时间:2024/07/30 15:11:31
用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
![用matlab求解常系数微分方程组并画出图像](/uploads/image/z/704156-68-6.jpg?t=%E7%94%A8matlab%E6%B1%82%E8%A7%A3%E5%B8%B8%E7%B3%BB%E6%95%B0%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E5%B9%B6%E7%94%BB%E5%87%BA%E5%9B%BE%E5%83%8F)
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结果
![](http://img.wesiedu.com/upload/f/e7/fe7fb19ee7e28e56acb3f7017c530ed3.jpg)
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结果
![](http://img.wesiedu.com/upload/f/e7/fe7fb19ee7e28e56acb3f7017c530ed3.jpg)