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

def gm11(x,n):app x1 = x.cumsum()    z1 = (x1[:len(x1) - 1] + x1[1:])/2.0#紧邻均值    z1 = z1.reshape((len(z1),1))    B = np.append(-z1,np.ones_like(z1),axis=1)    Y = x[1:].reshape((len(x) - 1,1))    [[a

>>阅读原文<<