《Matlab在数学建模中的应用》笔记3-灰色预测（GM（1，1））

灰色预测

今晚看了GM（1，1），发觉还不是很懂原理，貌似是：

• 先把原始数据累加，进而构建一阶微分；
• 再求均值生成新矩阵B和常数项矩阵Yn；
• 再利用最小二乘法（没给过程，直接公式求解，记忆中线代里有）求出所谓的灰参数a1；
• 再代回原微分方程，求出新序列x1
• 最后，对新序列的相邻数据求差，“还原”为拟合（预测）x0的X0序列

• 最后的最后，检验（这一步不是很懂，待以后弄懂更新）

  把原书代码修改了一下，加新注释以及检验部分（检验部分应该理解有误，到时记得再改吧）

clear;clc
syms a u;
a1=[a,u]';
X0=[89677,99215,109655,120333,135823,159878,182321,209407</