Logistic Map

Logistic Map的数学表示形式为：

clear,clc
% when a = 2.8
x0 = 0.2;
n = 1000;
a = 2.8;
x(1) = a * x0 * (1 - x0);
for i = 2:100
    x(i) = a * x(i-1) * (1 - x(i-1));
end
x = [x0,x];
t = 1:1:101;
subplot(2,2,1)
plot(t,x)
txt = {'a = 2.8','x0 = 0.2'};
text(4,1,txt)
set(gca,'xlim',[min(t),max(t)],'ylim',[0,1.2])
xlabel('n');
ylabel('x(n)');
clear x x0

% when a = 3.2
x0 = 0.2;
n = 1000;
a = 3.2;
x(1) = a * x0 * (1 - x0);
for i = 2:100
    x(i) = a * x(i-1) * (1 - x(i-1));
end
x = [x0,x];
t = 1:1:101;
subplot(2,2,2)
plot(t,x)
txt = {'a = 3.2','x0 = 0.2'};
text(4,1,txt)
set(gca,'xlim',[min(t),max(t)],'ylim',[0,1.2])
xlabel('n');
ylabel('x(n)');
clear x x0

% when a = 3.5 
x0 = 0.2;
n = 1000;
a = 3.5;
x(1) = a * x0 * (1 - x0);
for i = 2:100
    x(i) = a * x(i-1) * (1 - x(i-1));
end
x = [x0,x];
t = 1:1:101;
subplot(2,2,3)
plot(t,x)
txt = {'a = 3.5','x0 = 0.2'};
text(4,1,txt)
set(gca,'xlim',[min(t),max(t)],'ylim',[0,1.2])
xlabel('n');
ylabel('x(n)');
clear x x0

% when a = 4
x0 = 0.2;
n = 1000;
a = 4;
x(1) = a * x0 * (1 - x0);
for i = 2:100
    x(i) = a * x(i-1) * (1 - x(i-1));
end
x = [x0,x];
t = 1:1:101;
subplot(2,2,4)
plot(t,x)
txt = {'a = 4','x0 = 0.2'};
text(4,1.05,txt)
set(gca,'xlim',[min(t),max(t)],'ylim',[0,1.2])
xlabel('n');
ylabel('x(n)');
clear x x0

结果