clc;clear;
x=1; %初始人口
k=20; %人口承载量
j=0; %计数用
r1=0; %作为画图的自变量
for r=1.8:0.01:3 %r为1.8到3,步长为0.01
j=j+1;%计数
n=ones(1,100); %赋初始值
for t=1:49 %按照离散型的logistic回归公式迭代
n(t+1)=(1+r)*n(t)-r/k.*n(t).^2;
end
for t=50:99 %画出50到99所有的点
n(t+1)=(1+r)*n(t)-r/k.*n(t).^2;
x(j)=n(t)/k; %映射到0到1.5上
r1(j)=r;
plot(r1,x,'b*');%画图
axis([1.8 3 0 1.5]);%设置坐标轴范围
grid on;%网格线
hold on; %一定要有,否则x只有最后一点
end
M(j)=getframe;%动态图
end
忙活一晚上,终于给弄出来了,可以好好去睡一觉了,看着结果莫大的快乐!
网友评论