利用Octave做混沌实验

clc;clear;
y=@(k,x)k*x-k*x^2;
for k=0:0.04:4
    x0=0.1;
    for i=1:100
        x0=y(k,x0);
        if i>50
        plot(k,x0,'.b')
        hold on;
         end
    end
end

 

当0<a<1时，系统只有一个稳定点，且为零点。当1<a<3.2时，该方程迭代结果出现一个分支，趋于一个稳定点，且方程稳定。3<a<3.6时，方程的迭代出现两个分支，出现周期倍增。当3.6<a<4时，Feigenbaum图进入混沌状态，不存在极限点。有图像可知，an各分支点估计值为3，3.4，3.6，3.7，3.85，为4，2，1.75，预测该值收敛。