作业帮matlab求解x1’=a*x1-b*x1*x2 X2’=c*x1*x2-d*x2 其中a=2,b=0.01,c=0.001,d=0.7.x1(0)=300,x2(0)=100function dx=preyer(t,x,flag,a,b,c,d);dx(1,1)=a*x(1)-b*x(1)*x(2);dx(2,1)=c*x(1)*x(2)-d*x(2);tspan=[0,14];x0=[300,100];a=2;b=0.01;c=0.001;d=0.7;[t,x]=od
来源：学生作业帮助网 编辑：作业帮 时间：2024/11/25 07:43:37
![matlab求解x1’=a*x1-b*x1*x2 X2’=c*x1*x2-d*x2 其中a=2,b=0.01,c=0.001,d=0.7.x1(0)=300,x2(0)=100function dx=preyer(t,x,flag,a,b,c,d);dx(1,1)=a*x(1)-b*x(1)*x(2);dx(2,1)=c*x(1)*x(2)-d*x(2);tspan=[0,14];x0=[300,100];a=2;b=0.01;c=0.001;d=0.7;[t,x]=od](/uploads/image/z/11489245-61-5.jpg?t=matlab%E6%B1%82%E8%A7%A3x1%E2%80%99%3Da%2Ax1-b%2Ax1%2Ax2+X2%E2%80%99%3Dc%2Ax1%2Ax2-d%2Ax2+%E5%85%B6%E4%B8%ADa%3D2%2Cb%3D0.01%2Cc%3D0.001%2Cd%3D0.7.x1%280%29%3D300%2Cx2%280%29%3D100function+dx%3Dpreyer%28t%2Cx%2Cflag%2Ca%2Cb%2Cc%2Cd%29%3Bdx%281%2C1%29%3Da%2Ax%281%29-b%2Ax%281%29%2Ax%282%29%3Bdx%282%2C1%29%3Dc%2Ax%281%29%2Ax%282%29-d%2Ax%282%29%3Btspan%3D%5B0%2C14%5D%3Bx0%3D%5B300%2C100%5D%3Ba%3D2%3Bb%3D0.01%3Bc%3D0.001%3Bd%3D0.7%3B%5Bt%2Cx%5D%3Dod)
matlab求解x1’=a*x1-b*x1*x2 X2’=c*x1*x2-d*x2 其中a=2,b=0.01,c=0.001,d=0.7.x1(0)=300,x2(0)=100function dx=preyer(t,x,flag,a,b,c,d);dx(1,1)=a*x(1)-b*x(1)*x(2);dx(2,1)=c*x(1)*x(2)-d*x(2);tspan=[0,14];x0=[300,100];a=2;b=0.01;c=0.001;d=0.7;[t,x]=od
matlab求解x1’=a*x1-b*x1*x2 X2’=c*x1*x2-d*x2 其中a=2,b=0.01,c=0.001,d=0.7.x1(0)=300,x2(0)=100
function dx=preyer(t,x,flag,a,b,c,d);
dx(1,1)=a*x(1)-b*x(1)*x(2);
dx(2,1)=c*x(1)*x(2)-d*x(2);
tspan=[0,14];
x0=[300,100];
a=2;
b=0.01;
c=0.001;
d=0.7;
[t,x]=ode45('preyer',tspan,x0,a,b,c,d);
figure;
plot(t,x(:,1),'k-');
hold on;
plot(t,x(:,2),'k');
set(gca,'Fontsize',12);
xlabel('\itt','Fontsize',16);
L=legend('{\itx}_1','{\itx}_2',0);
set(L,'Fontname','Times New Roman');
set(gcf,'color','w',[121 233 1121 406]);
这个程序哪里不对,最好把改好后的程序和图像一起,
matlab求解x1’=a*x1-b*x1*x2 X2’=c*x1*x2-d*x2 其中a=2,b=0.01,c=0.001,d=0.7.x1(0)=300,x2(0)=100function dx=preyer(t,x,flag,a,b,c,d);dx(1,1)=a*x(1)-b*x(1)*x(2);dx(2,1)=c*x(1)*x(2)-d*x(2);tspan=[0,14];x0=[300,100];a=2;b=0.01;c=0.001;d=0.7;[t,x]=od
1.新建function函数
function dx=preyer(t,x);
a=2;b=0.01;c=0.001;d=0.7;
dx=zeros(2,1);
dx(1)=a*x(1)-b*x(1)*x(2);
dx(2)=c*x(1)*x(2)-d*x(2);
2.
[t,x]=ode45(@preyer,[0 14],[300 100]);
figure('position',[121 233 1121 406])
plot(t,x)
set(gcf,'color','w')
set(gca,'fontsize',12)
xlabel('\itt','fontsize',16)
ylabel('\itx','fontsize',16)
L=legend('{\itx}_1','{\itx}_2',0);
legend boxoff
set(L,'fontname','Times New Roman')
看看是不是你想要的结果