numpy.polynomial.laguerre.lagfit¶

numpy.polynomial.laguerre.lagfit(x, y, deg, rcond=None, full=False, w=None)[source]¶

Laguerre系列的最小二乘拟合数据。

返回作为对在点x给出的数据值y的最小二乘拟合的度deg的Laguerre系数的系数。如果y是1-D，则返回的系数也将是1-D。如果y是2-D多重拟合，对于y的每一列进行一次，并且所得到的系数存储在2-D返回的相应列中。拟合的多项式是形式

其中n是deg。

参数：	x：array_like，shape（M，） M个采样点`（x [i]， y [i]）`的x坐标。 y：array_like，shape（M，）或（M，K） y坐标。通过传递每列包含一个数据集的2D阵列，可以一次拟合共享相同x坐标的样本点的若干数据集。 deg：int或1-D array_like 拟合多项式的度（s）。如果deg是单个整数，则包括deg项的所有项包括在拟合中。对于Numpy版本> = 1.11，可以使用指定要包括的术语的度数的整数列表。 rcond：float，可选相对条件编号。相对于最大奇异值小于该值的奇异值将被忽略。默认值为len（x）* eps，其中eps是float类型的相对精度，在大多数情况下约为2e-16。 full：bool，可选开关确定返回值的性质。当它为False（默认值）时，只返回系数，当来自奇异值分解的True诊断信息也返回时。 w：array_like，shape（M，），可选重量。如果不是无，则通过w [i]加权每个点`(x[i],y[i])`对拟合的贡献。理想地，选择权重使得乘积`w[i]*y[i]`的误差都具有相同的方差。默认值为“无”。
返回：	coef：ndarray，shape（M，）或（M，K） Laguerre系数从低到高排序。如果y是2-D，则y的列k中的数据的系数在列k中。 [residuals，rank，singular_values，rcond]：list 只有full = True时，才会返回这些值 resid – sum of squared residuals of the least squares fit rank – the numerical rank of the scaled Vandermonde matrix sv – singular values of the scaled Vandermonde matrix rcond – value of rcond. 有关详细信息，请参阅linalg.lstsq。
警告：	RankWarning 最小二乘法拟合中的系数矩阵的秩是不足的。只有在满 = False时，才会发出警告。警告可以通过关闭 >>> import warnings >>> warnings.simplefilter('ignore', RankWarning)

参数：

x：array_like，shape（M，）

M个采样点（x [i]， y [i]）的x坐标。

y：array_like，shape（M，）或（M，K）

y坐标。通过传递每列包含一个数据集的2D阵列，可以一次拟合共享相同x坐标的样本点的若干数据集。

deg：int或1-D array_like

拟合多项式的度（s）。如果deg是单个整数，则包括deg项的所有项包括在拟合中。对于Numpy版本> = 1.11，可以使用指定要包括的术语的度数的整数列表。

rcond：float，可选

相对条件编号。相对于最大奇异值小于该值的奇异值将被忽略。默认值为len（x）* eps，其中eps是float类型的相对精度，在大多数情况下约为2e-16。

full：bool，可选

开关确定返回值的性质。当它为False（默认值）时，只返回系数，当来自奇异值分解的True诊断信息也返回时。

w：array_like，shape（M，），可选

重量。如果不是无，则通过w [i]加权每个点(x[i],y[i])对拟合的贡献。理想地，选择权重使得乘积w[i]*y[i]的误差都具有相同的方差。默认值为“无”。

coef：ndarray，shape（M，）或（M，K）

Laguerre系数从低到高排序。如果y是2-D，则y的列k中的数据的系数在列k中。

[residuals，rank，singular_values，rcond]：list

只有full = True时，才会返回这些值

resid – sum of squared residuals of the least squares fit rank – the numerical rank of the scaled Vandermonde matrix sv – singular values of the scaled Vandermonde matrix rcond – value of rcond.

有关详细信息，请参阅linalg.lstsq。

警告：

RankWarning

最小二乘法拟合中的系数矩阵的秩是不足的。只有在满 = False时，才会发出警告。警告可以通过关闭
>>> import warnings
>>> warnings.simplefilter('ignore', RankWarning)

也可以看看

chebfit，legfit，polyfit，hermfit，hermefit

lagval: 评估Laguerre系列。
lagvander: Laguerre系列的假Vandermonde矩阵。
lagweight: Laguerre重量函数。
linalg.lstsq: 从矩阵计算最小二乘拟合。
scipy.interpolate.UnivariateSpline: 计算样条拟合。

笔记

解决方案是Laguerre系列p的系数，其使加权平方误差的和最小化

其中 $w_j$ 是权重。这个问题通过设置为（通常）过度确定的矩阵方程来解决

其中V是x的加权伪Vandermonde矩阵，c是要求解的系数，w权重和y是观察值。然后使用V的奇异值分解来求解该方程。

如果V的一些奇异值如此小以至于被忽略，则将发出RankWarning。这意味着可能不良地确定系数值。使用较低的顺序通常会摆脱警告。rcond参数也可以设置为小于其默认值的值，但是所得到的拟合可能是假的并且具有来自舍入误差的较大贡献。

使用Laguerre系列的拟合可能是最有用的，当数据可以近似为sqrt（w（x）） * p >，其中w（x）是Laguerre权重。在这种情况下，应当与数据值y[i]/sqrt(w(x[i])一起使用权重sqrt(w(x[i])。权重函数可用作lagweight。

参考文献

[R63]

维基百科，“曲线拟合”，http://en.wikipedia.org/wiki/Curve_fitting

例子

>>> from numpy.polynomial.laguerre import lagfit, lagval
>>> x = np.linspace(0, 10)
>>> err = np.random.randn(len(x))/10
>>> y = lagval(x, [1, 2, 3]) + err
>>> lagfit(x, y, 2)
array([ 0.96971004,  2.00193749,  3.00288744])