prediction.m

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

function p=prediction(x,m,t,xx,nn)
n=length(x);
N=n-(m-1)*t;
for i=1:N %将时间序列嵌入m维相空间
      for j=l:m
          Y(i,j)=x(i+(j-1)*t);
      end
end
 s=1000;k=1;
for i=1:N-1
      d=dadd(y(i,:)-y(N,:),m);
          if(d<s)
              s=d;
             k=i;
          end 
end 
%Li=Lyapunov(y,m,N)
Li=0.6694;
e=exp(Li);
s=e*(y(k,1)-y(N,1))^2;
for i=2:m
    s=(e-1)*(y(k,i)-y(N,i))^2>
end
n+1 %预测向量点
num1=k+(m-1)*t %最近邻近点
YC=x(k+(m-1)*t+1) %预测值
Zhen=xx(nn+1)真值
GuJierror=sqrt(s)
error=YC-Zhen
%估计误差
%实际预测误差