【Matlab】产生混沌序列并绘制气泡图

%% 产生混沌序列
x0 = 0.2983;
len = 1000;
squence = zeros(1,len);
xn = x0;
for i=1:len
    xn = 4*xn*(1-xn);
	squence(i) = xn;
end

%% 绘制混沌气泡图
y = linspace(1,len,len);
x = squence;
sz = linspace(50,200,length(x));
c = linspace(0,1,length(x));
scatter(x,y,sz,c,'filled');

function []=hd()
    %% 产生混沌序列
    % 3.569946~4
    t = 4;
    x1 = 0.9952;
    len1 = 2000;
    squence1 = zeros(1,len1);
    xn = x1;
    for i=1:len1
        xn = t*xn*(1-xn);
    	squence1(i) = xn;
    end
    
    x2 = 0.1123;
    len2 = 2000;
    squence2 = zeros(1,len2);
    xn = x2;
    for i=1:len2
        xn = t*xn*(1-xn);
    	squence2(i) = xn;
    end
    
    %% 绘制混沌气泡图
    x = squence2;
    y = squence1;
    %sz = linspace(50,200,length(x));
    %c = linspace(0,1,length(x));
    %scatter(x,y,sz,c,'filled');
    plotmatrix([x',y']);
end

matlab plotmatrix函数

function [] = ck()
%% -------------------虫口模型-分岔图-----------
% 虫口模型 x(n+1) = lambda * x(n) * (1 - x(n))
    lambda = 0:5e-3:4;
    N = 300;
    f = zeros(1,N);
    for j=1:length(lambda)
        x = 0.4;     %初值
        for i = 1:N
            x = lambda(j) * x * (1 - x);
            f(i) = x;
        end
        plot(repmat(lambda(j),1,length(f)),f,'k.');
        hold on;
    end

    xlabel('\lambda');
    ylabel('x');
    title('Logistic虫口模型-虫口随着λ的变化');
end

虫口模型---昆虫变化的简单数学模型。（昆虫一次产卵的数量极大，子代出生后，其亲代的数量可以忽略不计）认为昆虫的种群是由单一的世代构成，在历代之间没有重叠。因此，昆虫种群的增长（虫口数）是分步进行的，例如，每个3年、5年或7年出现一次成虫期的蝉就是这样。