huise.rar_GM_GREY (1_1)_grey forecast_灰色预测_灰色预测1.1

function GM=huise(data,N) data1=1e+4*[9.55 11.95 15.28 18.03 21.52 26.72]; data=[32.5 34.5 35.4 40.2 42.9 42.1 45.2 43.2 43.2 43.7 45.5 43.4 44.1 46.6 47.1]; data=[5 4 2 6 4 4 3 9 4 1 2 7 7 2 5]; data=[6004.21 6680.34 7199.95 7697.82 8582.73] T=length(data); N=1; X0=data; for i=2:T X1(1)=X0(1); X1(i)=X1(i-1)+X0(i); %用AGO生成一阶累加生成模块 end for i=1:T-1 M(i)=-(0.5*(X1(i)+X1(i+1))); end B=zeros(T-1,2); %构造累加矩阵B for i=1:T-1 for j=1:2 if j<2 B(i,j)=M(i); elseif j>1 B(i,j)=1; end end end for i=2:T %构造常数项向量Y Y(i-1)=X0(i); end HCS=inv(B'*B)*B'*Y' %用最小二乘法求灰参数HCS H=HCS'; XR1(1)=X0(1); for i=2:T+N %预测出累加序列 XR1(i)=(X0(1)-H(2)/H(1))*exp(-1*H(1)*(i-1))+H(2)/H(1); end K(1)=XR1(1); for i=2:T+N %还原计算出预测值 K(i)=XR1(i)-XR1(i-1); end GM=K; for i=1:T E(i)=X0(i)-GM(i); end for i=1:T Q(i)=E(i)/X0(i); end