种群相互竞争模型

种群竞争模型、截图加代码

种群相互竞争模型

调用函数：

function xdot=shier(t,x)

r1=2.5;r2=1.8;q1=0.5;q2=1.6;n1=1.6;n2=1;

xdot=[r1*x(1).*(1-x(1)./n1-q1.*(x(2)./n2));r2*x(2).*(1-q2*x(1)./n1-x(2)./n2)]; 计算：

s=0:0.1:8;

x0=[0.1,0.1];

[t,x]=ode45('shier',ts,x0);[t,x],

axis([0 8, 0 2]),

hold on

plot(t,x),grid,gtext('x1(t)'),gtext('x2(t)'),

pause,

plot(x(:,1),x(:,2)),

grid,

结果截图：

调用函数：

function xdot=shier(t,x)

r1=2.5;r2=1.8;q1=0.5;q2=1.6;n1=1.6;n2=1;

种群竞争模型、截图加代码

xdot=[r1*x(1).*(1-x(1)./n1-q1.*(x(2)./n2));r2*x(2).*(1-q2*x(1)./n1-x(2)./n2)];

计算：

s=0:0.1:8;

x0=[1,2];

[t,x]=ode45('shier',ts,x0);[t,x],

axis([0 8, 0 2]),

hold on

plot(t,x),grid,gtext('x1(t)'),gtext('x2(t)'),

pause,

plot(x(:,1),x(:,2)),

grid,

结果截图：

调用函数：

function xdot=shier(t,x)

r1=2.5;r2=1.8;q1=0.5;q2=1.6;n1=1.6;n2=1;

xdot=[r1*x(1).*(1-x(1)./n1-q1.*(x(2)./n2));r2*x(2).*(1-q2*x(1)./n1-x(2)./n2)];

计算：

ts=0:0.1:8;

x0=[0.1,0.1];

axis([0 2, 0 2]),

种群竞争模型、截图加代码

hold on

[t,x]=ode45('shier',ts,x0); hold on

plot(x(:,1),x(:,2)),

hold on

x0=[1,2];

[t,x]=ode45('shier',ts,x0); hold on

plot(x(:,1),x(:,2)),

hold on

x1=0:0.1:1;

y1=-x1+1;

x2=0:0.1:2;

y2=-1.25*x2+2;

plot(x1,y1,'--',x2,y2,'--') grid

结果截图：