logistic.m

混沌系统需要的各种仿真程序,希望对大家(特别是初学者)有用.

function  Logistic

% Logistic映射程序
n=40;
u=4;
X(1)=0.2;                      % X取值范围是[0,1]
for i=1:n-1
    X(i+1)=u*X(i)*(1-X(i));  % Logistic映射程序公式
end;
x=mean(X)            % x为数组的平均值

for i=1:n            % 量化数组
    if X(i)>=x
        X(i)=0
    else  
        X(i)=1
        
    
    end
   a(i)=X(i)          % 定义数组y（i）为a量化后的新数组
end
 
% 量化后的数组调用计算步数的程序

leth=length(a)
 num=zeros(leth,1); %存放各个exhausting history中的comment的长度
 
if n>=500           % 数目较小的可直接运行
    for i=1:leth
        num(i)=0
    end;               % 分配内存，加快运行速度
end;

 t=2;

n=length(a);        % n为矩阵a的长度
e=0;                % e为步长，先清零
p=1;                % p为指针
if a(2)==a(1)
    e=e+1;
    p=p+1;
    num(1)=1;        % 前两个数的比较，相同则是表示断点在1;
else
    e=e+1;
    p=p+2;
    num(1)=2;
end;                  %  前两个数的比较，不同则是表示断点在2;

while p<n
    ps=p;    % ps为指针
    el=0;    % 一次小循环内每次比较两数组值后所得exhausting history长度；穷尽步长
    emax=0;  % 寻找一个exhausting history中所有el的最大值；最长步长的字数
    pmax=0;  
    for i=1:ps-1
         j=i;
         while(a(j)==a(p)) % 作比较,逻辑比较
             j=j+1;        
             p=p+1;         
             el=el+1; 
             if p==n       % 逻辑比较
                 p=p-1;
                 break;
             end;
         end;
                          
         if p>=pmax
             pmax=p;
         end;
         el=0;    
         p=ps;
     end;
     
     num(t)=pmax+1-ps;
     t=t+1;
    p=pmax+1;
   e=e+1;
end;

if p==n
    e=e+1;
    num(t)=1
end;

disp('*****');
disp(e);
% disp(num);