GM(1,1)算法代码（python）

def gm11(x,n):

1. x1 = x.cumsum()  
2.  z1 = (x1[:len(x1) - 1] + x1[1:])/2.0#紧邻均值  
3.  z1 = z1.reshape((len(z1),1))  
4.  B = np.append(-z1,np.ones_like(z1),axis=1)  
5.  Y = x[1:].reshape((len(x) - 1,1))  
6.  [[a],[b]] = np.dot(np.dot(np.linalg.inv(np.dot(B.T, B)), B.T), Y)#计算参数  
7.  return (x[0]-b/a)*np.exp(-a*(n-1))-(x[0]-b/a)*np.exp(-a*(n-2))