用matlab求解常系数微分方程组并画出图像dx/dt=220-0.25x-0.000015xy²dy/dt=0.000015xy²-0.64y+0.15zdz/dt=0.64y-0.15z-0.7z初始值为x=68,y=109,z=82跪求好人,求解和画图,时间t定义为0-60,需要代码,跪谢了
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/16 22:03:01
用matlab求解常系数微分方程组并画出图像dx/dt=220-0.25x-0.000015xy²dy/dt=0.000015xy²-0.64y+0.15zdz/dt=0.64y-0.15z-0.7z初始值为x=68,y=109,z=82跪求好人,求解和画图,时间t定义为0-60,需要代码,跪谢了
用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,需要代码,跪谢了
用matlab求解常系数微分方程组并画出图像dx/dt=220-0.25x-0.000015xy²dy/dt=0.000015xy²-0.64y+0.15zdz/dt=0.64y-0.15z-0.7z初始值为x=68,y=109,z=82跪求好人,求解和画图,时间t定义为0-60,需要代码,跪谢了
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结果
以上回答你满意么?