stationary_measure_Vivian.py 文件源码

def measure(y):
    x=np.linspace(1,183,183)
    y_ex=[]
    y_ex=np.array(y_ex)

    pred=Pipeline([('poly',PolynomialFeatures(10)),
                   ('linear',LinearRegression(fit_intercept=False))])
    pred.fit(x[:,np.newaxis],y) 
    y_ex=pred.predict(x[:,np.newaxis])

    t=comp(y_ex,y)
        return t