gmnewpred.m

这是一个有预测功能的程序

function [Predict]=gmnewpred(x,T,N)
n=length(x);
m=N+n;
xT_1=gmT(T-1,x);
xT=gmT(T,x);
q=getab(xT_1,xT);
a=q(1);
b=q(2);
for i=1:m
  xm1(i)=(x(1)-b/a)*exp(-a*(i-1))+b/a;
end
  xm0=xm1;
Xm0=gm_T(T,xm0);
Err=Xm0(2:n)-x(2:n);
C=std(Err)/std(x);
%RelErr=abs((Err)mean(Err)*ones(size(Err)));
RelErr=abs(Err-mean(Err)*ones(size(Err)));
[Nrow, Ncol]=size(RelErr);
count=0;
for j=1:Ncol
    if RelErr(1:j)<0.6745*std(x)
        count=count+1;
    else count=count;
    end
end
P=count/Ncol;
Predict=Xm0((length(Xm0)-N+1):length(Xm0));