fun.m

阻滞增长模型

% 本程序为差分形式的阻滞增长模型
function f=fun(n)%n为迭代次数
x=0:0.01:1;
b(1)=2.5;
hold on;
x(1)=0.2000;%选定初值为0.2
for j=1:100
    b(j+1)=b(j)+0.015;
    for i=1:n
       x(i+1)=b(j)*x(i)*(1-x(i));
       if i>50   %画出50次以后的点
            plot(b(j),x(i));
            legend('收敛，分岔和混沌');
       end
    end  
end
hold off;