一、种群竞争模型

二、分析

两个种群：
自然增长率r1 = 1、r2 = 1
初始数量 n1 = 100、n2 = 100
资源竞争力 s1 = 0.5、s2 = 2

（1）未改变初值

（2）改变自然增长率r

（3）改变该环境种群最大容量

（4）改变两个种群初始数量

（5）改变资源竞争力

三、MATLAB执行代码

fun.m:function dx=fun(t,x,r1,r2,n1,n2,s1,s2)
r1=1;
r2=1;
n1=100;
n2=100;
s1=0.5;
s2=2;
dx=[r1*x(1)*(1-x(1)/n1-s1*x(2)/n2);r2*x(2)*(1-s2*x(1)/n1-x(2)/n2)];p3.m:h=0.1;%所取时间点间隔
ts=[0:h:30];%时间区间
x0=[10,10];%初始条件
opt=odeset('reltol',1e-6,'abstol',1e-9);%相对误差1e-6，绝对误差1e-9
[t,x]=ode45(@fun,ts,x0,opt);%使用5级4阶龙格—库塔公式计算
plot(t,x(:,1),'r',t,x(:,2),'b','LineWidth',2),grid;
pause;
plot(x(:,1),x(:,2),'LineWidth',2),grid  %作相轨线