numpy.polynomial.polynomial.Polynomial.fit#

方法

classmethod polynomial.polynomial.Polynomial.fit(x, y, deg, domain=None, rcond=None, full=False, w=None, window=None, symbol='x')[source]#

数据的最小二乘拟合。

返回一个级数实例,该实例是对在 x 处采样的 y 数据的最小二乘拟合。可以指定返回实例的域,这通常会产生更好的拟合,并减少病态条件的机会。

参数:
x类数组,形状 (M,)

M 个样本点 (x[i], y[i]) 的 x 坐标。

y类数组,形状 (M,)

M 个样本点 (x[i], y[i]) 的 y 坐标。

degint 或 1-D 类数组

拟合多项式的度。如果 deg 是单个整数,则拟合中包含直至并包括 deg 次的所有项。对于 NumPy 版本 >= 1.11.0,可以使用整数列表来指定要包含的项的度数。

domain{None, [beg, end], []},可选

用于返回级数的域。如果为 None,则选择覆盖点 x 的最小域。如果为 [],则使用类域。默认值在 NumPy 1.4 中是类域,在后续版本中是 None[] 选项在 numpy 1.5.0 中添加。

rcond浮点数,可选

拟合的相对条件数。小于此值(相对于最大奇异值)的奇异值将被忽略。默认值为 len(x)*eps,其中 eps 是浮点类型的相对精度,在大多数情况下约为 2e-16。

full布尔值,可选

决定返回值的性质的开关。当为 False(默认值)时,只返回系数;当为 True 时,还会返回来自奇异值分解的诊断信息。

w类数组,形状 (M,),可选

权重。如果不是 None,权重 w[i] 适用于 x[i] 处未平方的残差 y[i] - y_hat[i]。理想情况下,应选择权重,使乘积 w[i]*y[i] 的误差都具有相同的方差。当使用逆方差加权时,使用 w[i] = 1/sigma(y[i])。默认值为 None。

window{[beg, end]},可选

用于返回级数的窗口。默认值是默认类域

symbol字符串,可选

表示自变量的符号。默认为 'x'。

返回:
new_series级数

表示数据最小二乘拟合的级数,具有调用中指定的域和窗口。如果对未缩放和未平移的基多项式的系数感兴趣,请执行 new_series.convert().coef

[resid, rank, sv, rcond]列表

这些值仅在 full == True 时返回

  • resid – 最小二乘拟合的残差平方和

  • rank – 缩放后的范德蒙矩阵的数值秩

  • sv – 缩放后的范德蒙矩阵的奇异值

  • rcond – rcond 的值。

更多详细信息,请参阅 linalg.lstsq