Logistic Map
x k + 1 = β ∗ x k ∗ ( 1 − x k ) x_{k+1} = \beta*x_k*(1-x_k) xk+1=β∗xk∗(1−xk)
选择不同的 β \beta β可以获得混沌的逻辑斯特图。对应的MATLAB程序如
clear all;close all;clc
x0 = 0.5;
xvals = [];
for beta = 0:0.01:4
xold = x0;
% transient
for i=1:2000
xnew = beta*xold*(1-xold);
xold = xnew;
end
xss = xnew;
for i = 1:1000
xnew = beta*xold*(1-xold);
xold = xnew;
xvals(1,length(xvals)+1) = beta;
xvals(2,length(xvals)) = xnew;
if(abs(xnew-xss)<0.001)
break
end
end
end
plot(xvals(1,:),xvals(2,:),'.r','LineWidth',.5,'MarkerSize',1.5);
计算获得的图为: