x=[];%時(shí)間序列的輸入
n=length(x);
t=3;%t(time dely)
de=0.5;%error threshold
DD(I)=0; m1=0;
for m=2:8
    N=n-(m-1) * t;
      for i=1:N
            for I=I:m
              y(i,l)=x(i+(l-I)*t);
              end
      end
    D=GP(y,m,N);%調(diào)用子函數(shù)
    if (abs(DD(m-1)-D)>de)
          m1=m; yy=y;
      end
    DD(m)=D;
end
lyapunov=lyapunov(yy,m1,N)%調(diào)用子函數(shù)
plot(DD,'b:.');
plot(mI,DD(ml),'r0');