4.待定系数矩阵[a;b]=inv(B'*B)*B'*Y,其中inv表示求逆，B'表示转置
5.得到X0（k+1）=（1-exp（a））*（X0（1）-b/a）*exp（-a*k） 取k=7，8，……，19可得预测值
以上只是计算过程，具体原理可以查阅相关论文。另外，可以考虑将X0（k）先做差分，对差分序列进行灰色预测后再还原。
125.73657 141.55648 159.36681 179.41801 201.99201 227.40621 411.28518
256.01799 463.03222 288.22963 521.28995 324.49407 365.32123