混沌系统的分岔图画法

大家好，请教一下，混沌系统的分岔图怎么画？有没有人做混沌保密通信的?想结识一下，课题组就我一人做这个。

你说的是逻辑斯蒂映射吧

clear
scale = 10000; % determines the level of rounding
maxpoints = 200; % determines maximum values to plot
N = 3000; % number of "r" values to simulate
a = 2.0; % starting value of "r"
b = 4; % final value of "r"... anything higher diverges.
rs = linspace(a,b,N); % vector of "r" values
M = 500; % number of iterations of logistics equation
% Loop through the "r" values
for j = 1:length(rs)
    
    r=rs(j); % get current "r"
    x=zeros(M,1); % allocate memory
    x(1) = 0.5; % initial condition (can be anything from 0 to 1)
    
    for i = 2:M, % iterate
        x(i) = r*x(i-1)*(1-x(i-1));
    end
    % only save those unique, semi-stable values
    out{j} = unique(round(scale*x(end-maxpoints:end)));
end
% Rearrange cell array into a large n-by-2 vector for plotting
data = [];
for k = 1:length(rs)
    n = length(out{k});
    data = [data;  rs(k)*ones(n,1),out{k}];
end
% Plot the data
figure(97);clf
h=plot(data(:,1),data(:,2)/scale,'k.');
set(h,'markersize',1)
axis tight
set(gca,'units','normalized','position',[0 0 1 1])
set(gcf,'color','white')
axis off

分岔图产生的条件其实就是方程出现不止一个解（2、4、8个等等），也就是一个x对应多个y。而混沌则可以简单的理解为x对应无数多个y。
要学习画分插图、乃至混沌图，可以先从最简单的一维系统入手，在经典系统上增删项，调整迭代次数看呈现的效果
（做过一年左右的混沌构造研究，《复杂网络算法与应用》这本书对入门挺有用的）