种群竞争模型

一、模型背景
当两个种群为争夺同一食物来源和生存空间相互竞争时，常见的结局是，竞争力弱的灭绝，竞争力强的达到环境容许的最大容量。

使用种群竞争模型可以描述两个种群相互竞争的过程，分析产生各种结局的条件。

二、数学理论
1.甲乙两个种群，它们独自生存时数量变化均服从Logistic规律

2.两种群在一起生存时，乙对甲增长的阻滞作，用与乙的数量成正比；甲对乙有同样作用乙

其中x(t)，y(t)分别为甲乙两种群的数量，r1r2为它们的固有增长率，n1n2为它们的最大容量。s1s1的含义是对于供养甲的资源来说，单位数量的乙（相对n2n2）的消耗为单位数量甲（相对n1n1）消耗的s1s1倍，s2s2同理。

三、源程序
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=[r1x(1)(1-x(1)/n1-s1x(2)/n2);r2x(2)(1-s2x(1)/n1-x(2)/n2)];

% r：自然增长率 n：物种容量 s：生存能力（当s1s2都小于1且相差不大可近似认为两物种共存） % 以上都是可变量（该模型也可套用于企业产品销售）

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 %作相轨线

% 召唤图2按任意键，图2中横坐标为种群1数量，纵坐标为种群2数量 % 另外，为使图像标准，可以在图中插入箭头及注释。 % x：物种初值

原文链接：https://blog.csdn.net/weixin_40108753/article/details/81282157