chaos.m

来自「Chaos (logistic function) test」· M 代码 · 共 13 行

% Chaos test
% logistic function plot
i = disp('please enter the number of iteration : ');
cx = [];
cx(1) = rand(1);
if (cx(1) ~= 0.25) & (cx(1) ~= 0.5) & (cx(1) ~= 0.75)
    for k = 1: 30
        cx(k+1) = 4*cx(k)*(1-cx(k));
    end
else
    break
end
plot (cx)