混沌时间序列的最大李雅普诺夫指数问题

原来是小飞

使用程序运行出了最大李雅普诺夫指数的图片，但不知道如何从图片中看出最大李雅普诺夫指数，特此请教帮助，图片和代码如下：

1. clear all;
   
2. data=load('b.txt');
   
3. N=length(data);
   
4. m=2;
   
5. tau=6;
   
6. P=6;
   
7. delt_t=1;
   
8. Y=reconstitution(data,N,m,tau );%reconstitute state space
   
9. M=N-(m-1)*tau;%M is the number of embedded points in m-dimensional space
   
10. for j=1:M
    
11.     d_max=1e+100;
    
12.     for jj=1:M                                              %寻找相空间中每个点的最近距离点，并记下
    
13.         d_s=0;                                              %该点下标
    
14.         if abs(j-jj)>P                                      %限制短暂分离
    
15.             for i=1:m
    
16.                 d_s=d_s+(Y(i,j)-Y(i,jj))*(Y(i,j)-Y(i,jj));
    
17.                 d_min=d_max;
    
18.                 if d_s<d_min
    
19.                    d_min=d_s;
    
20.                    idx_j=jj;
    
21.                end
    
22.             end
    
23.         end
    
24.     end
    
25. %     index(j)=idx_j;
    
26.     max_i=min((M-j),(M-idx_j));%计算点j的最大演化时间步长i
    
27.     for k=1:max_i              %计算点j与其最近邻点在i个离散步后的距离
    
28.         d_j_i=0;
    
29.         for kk=1:m
    
30.             d_j_i=d_j_i+(Y(kk,j+k)-Y(kk,idx_j+k))*(Y(kk,j+k)-Y(kk,idx_j+k));
    
31.             d(k,j)=d_j_i;
    
32.         end
    
33.     end
    
34. end
    
35. 
    
36. %对每个演化时间步长i，求所有的j的lnd(i,j)平均
    
37. [l_i,l_j]=size(d);
    
38. for i=1:l_i
    
39.     q=0;
    
40.     y_s=0;
    
41.     for j=1:l_j
    
42.         if d(i,j)~=0
    
43.             q=q+1;
    
44.             y_s=y_s+log(d(i,j));
    
45.         end
    
46.     end
    
47.     y(i)=y_s/(q*delt_t);
    
48. end
    
49. x=1:length(y);
    
50. pp=polyfit(x,y,1);
    
51. lambda_1=pp(1);
    
52. yp=polyval(pp,x);
    
53. plot(x,y,'-o',x,yp,'--');grid on

复制代码

图片如下：

huadianlizhan

为啥不行那/