logistic规律可以用以下式子表示
N(t+1)=k*N(t)*(N-N(t))
其中N(t)为t时刻种群的数量,N为环境最大容纳的的种群数
k的取值通常与N有关,大约在1/N和4/N之间
题目如下
假设一些参数后的代码如下
T=100;
T1=1:T;
N1=1000;
a1=3/N1;
Q=800;
N2=700;
a2=-0.2/N2;
q=720;
x0=200;
x=zeros(T,1);
x(1)=x0;
temp=0;%检验的标志
for t=2:T
if temp==0
if x(t-1)<Q
x(t)=a1*x(t-1)*(N1-x(t-1));
else
x(t)=a2*x(t-1)*(N2-x(t-1));
temp=1;
end
else
if x(t-1)>q
x(t)=a2*x(t-1)*(N2-x(t-1));
else
x(t)=a1*x(t-1)*(N1-x(t-1));
temp=0;
end
end
end
plot(T1,x);
效果如下