python - 获取 scipy 的 LSQUnivariateSpline 或 make_lsq_spline 的正态矩阵或方差-协方差矩阵（在给定权重的情况下）

Question

我正在使用 scipy.interpolate 中的 LSQUnivariateSpline 或 make_lsq_spline 来获取给定的权重。有没有一种简单的方法不仅可以获得结果样条曲线，还可以获得样条曲线系数的方差 - 协方差矩阵？有没有一种简单的方法来获取在这些函数中生成并用于 LSQ 问题的法线矩阵？

Answer 1

没有支持的方式，没有。

您当然可以从以下源代码中获取您需要的内容make_interp_spline：

https://github.com/scipy/scipy/blob/v1.5.2/scipy/interpolate/_bsplines.py#L854-L1009

（正常的 eq 矩阵是用 构造的_bspl._norm_eq_lsq）。但请注意，内部结构如有更改，恕不另行通知。

特别是，正规方程的使用可能会消失。

对于 LSQUnivariateSpline，您需要阅读 FITPACK 的 Fortran 源代码。（有点）好消息是它在未来的 scipy 版本中不太可能发生变化。