function lyapunov(x).txt

来自「求解logistic映射的李雅普诺夫指数」· 文本 代码 · 共 38 行

function Lyapunov(x)
if length(x)<400
   helpdlg('数据量需要大于400，请重新选定时间','提示信息');pause(3),
end   
dm=26;tao=5;set_error=0.001;
for i=0:(dm-1)
    for j=1:(length(x)/1.2)
       k=j+i*tao;
           Y(i+1,j)=x(k);
    end
end 

l=1; i=1; p=1;                      %  i为Y的一列，j为Y的一列 
while i<length(Y)
     for j=p:length(Y)               %求两点之间距离
        if j>i
        e_normal=max(abs(Y(:,i)-Y(:,j)));
        end
     end    
    flag=e_normal-set_error;
            if flag>0
              d(l)=e_normal;l=l+1;
              %i=i+1;
                 for j=1:length(Y)             %找离点最小距离的点  
                    l_normal(j)=max(abs(Y(:,i)-Y(:,j)));
                 end
              p=find((l_normal-min(l_normal))==0);            %找到最近点Y(:,j)
           else 
              j=j+1;i=i+1;                    %范数小于set_error，Y(:,i)和Y(:,j)继续向前一步  
           end
end
 
LE1=0;
for i=1:l
    f1=d(i);f2=d(1);
 LE1=log10(f1/f2)+LE1;
end
LE1=1/(l*0.014)*LE1,