main_gp.asv

关于GP算法,不过好像有些问题.大家一起解决吧

x=[];%时间序列的输入
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);%调用子函数
    if (abs(DD(m-1)-D)>de)
          m1=m; yy=y;
      end
    DD(m)=D;
end
lyapunov=lyapunov(yy,m1,N)%调用子函数
plot(DD,'b:.');
plot(mI,DD(ml),'r0');