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/18 00:33:47
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')
看看是不是你想要的结果