Least squares fit of Laguerre series to data. Return the coefficients of a Laguerre series of degree `deg` that is the least squares fit to the data values `y` given at points `x`. If `y` is 1-D the returned coefficients will also be 1-D. If `y` is 2-D multiple fits are done, o
(x, y, deg, rcond=None, full=False, w=None)
| 1270 | |
| 1271 | |
| 1272 | def lagfit(x, y, deg, rcond=None, full=False, w=None): |
| 1273 | """ |
| 1274 | Least squares fit of Laguerre series to data. |
| 1275 | |
| 1276 | Return the coefficients of a Laguerre series of degree `deg` that is the |
| 1277 | least squares fit to the data values `y` given at points `x`. If `y` is |
| 1278 | 1-D the returned coefficients will also be 1-D. If `y` is 2-D multiple |
| 1279 | fits are done, one for each column of `y`, and the resulting |
| 1280 | coefficients are stored in the corresponding columns of a 2-D return. |
| 1281 | The fitted polynomial(s) are in the form |
| 1282 | |
| 1283 | .. math:: p(x) = c_0 + c_1 * L_1(x) + ... + c_n * L_n(x), |
| 1284 | |
| 1285 | where ``n`` is `deg`. |
| 1286 | |
| 1287 | Parameters |
| 1288 | ---------- |
| 1289 | x : array_like, shape (M,) |
| 1290 | x-coordinates of the M sample points ``(x[i], y[i])``. |
| 1291 | y : array_like, shape (M,) or (M, K) |
| 1292 | y-coordinates of the sample points. Several data sets of sample |
| 1293 | points sharing the same x-coordinates can be fitted at once by |
| 1294 | passing in a 2D-array that contains one dataset per column. |
| 1295 | deg : int or 1-D array_like |
| 1296 | Degree(s) of the fitting polynomials. If `deg` is a single integer |
| 1297 | all terms up to and including the `deg`'th term are included in the |
| 1298 | fit. For NumPy versions >= 1.11.0 a list of integers specifying the |
| 1299 | degrees of the terms to include may be used instead. |
| 1300 | rcond : float, optional |
| 1301 | Relative condition number of the fit. Singular values smaller than |
| 1302 | this relative to the largest singular value will be ignored. The |
| 1303 | default value is len(x)*eps, where eps is the relative precision of |
| 1304 | the float type, about 2e-16 in most cases. |
| 1305 | full : bool, optional |
| 1306 | Switch determining nature of return value. When it is False (the |
| 1307 | default) just the coefficients are returned, when True diagnostic |
| 1308 | information from the singular value decomposition is also returned. |
| 1309 | w : array_like, shape (`M`,), optional |
| 1310 | Weights. If not None, the weight ``w[i]`` applies to the unsquared |
| 1311 | residual ``y[i] - y_hat[i]`` at ``x[i]``. Ideally the weights are |
| 1312 | chosen so that the errors of the products ``w[i]*y[i]`` all have the |
| 1313 | same variance. When using inverse-variance weighting, use |
| 1314 | ``w[i] = 1/sigma(y[i])``. The default value is None. |
| 1315 | |
| 1316 | Returns |
| 1317 | ------- |
| 1318 | coef : ndarray, shape (M,) or (M, K) |
| 1319 | Laguerre coefficients ordered from low to high. If `y` was 2-D, |
| 1320 | the coefficients for the data in column *k* of `y` are in column |
| 1321 | *k*. |
| 1322 | |
| 1323 | [residuals, rank, singular_values, rcond] : list |
| 1324 | These values are only returned if ``full == True`` |
| 1325 | |
| 1326 | - residuals -- sum of squared residuals of the least squares fit |
| 1327 | - rank -- the numerical rank of the scaled Vandermonde matrix |
| 1328 | - singular_values -- singular values of the scaled Vandermonde matrix |
| 1329 | - rcond -- value of `rcond`. |