main_gp.m

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

      end
   D=GP(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(m1,DD(m1),'r0');