gmnew1.m

来自「这是一个有预测功能的程序」· M 代码 · 共 38 行

function []=gmnew1(x,T,N)
%对数列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));
C
P
a
b
x
Xm0
Predict