MATLAB动态程序：数学模型课程里的人口logistic离散模型

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

忙活一晚上，终于给弄出来了，可以好好去睡一觉了，看着结果莫大的快乐！