function p=prediction(x,m,t,xx,nn)
n=length(x);
N=n-(m-1)*t;
for i=1:N %將時(shí)間序列嵌入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 %預(yù)測(cè)向量點(diǎn)
num1=k+(m-1)*t %最近鄰近點(diǎn)
YC=x(k+(m-1)*t+1) %預(yù)測(cè)值
Zhen=xx(nn+1)真值
GuJierror=sqrt(s)
error=YC-Zhen
%估計(jì)誤差
%實(shí)際預(yù)測(cè)誤差