灰色预测代码及检验

x0=[113.73 114.98 121.98 130.88 138.68 145.64];%2013~2018年单独二胎政策下专家预测的新生儿数（万人）

x1=[x0(1),x0(1)+x0(2),x0(1)+x0(2)+x0(3),x0(1)+x0(2)+x0(3)+x0(4),x0(1)+x0(2)+x0(3 )+x0(4)+x0(5),x0(1)+x0(2)+x0(3)+x0(4)+x0(5)+x0(6)];

y = (x0(2:6))';

B=[-1/2*[x1(1)+x1(2) ;x1(2)+x1(3);x1(3)+x1(4) ;x1(4)+x1(5);x1(5)+x1(6)],[1;1;1;1;1] ] ;

(B'*B)^-1;

U = (B'*B)^-1*B'*y;

a=U(1);

u=U(2);

for k= 0:8;

X1(k+1)=(x1(1)-u/a)*exp(-a*k)+u/a;

end;

x2020=X1(8)-X1(7) %二胎政策下的2020年人口预计值

X0=[X1(1) X1(2)-X1(1) X1(3)-X1(2) X1(4)-X1(3) X1(5)-X1(4) X1(6)-X1(5)]

E = x0(2:6) - X0(2:6)

e = E./x0(2:6)

Ax0 = mean(x0)

Vx0=var(x0)

AE = mean(E)

VE= var(E)

由上表可见，数据检验表明最大误差为0.73％，均小于1%，拟合精度较高。后验差C=0.0035=0.35%，预测精度好。