对于给定的权重,我使用的是scipy.interpolate中的LSQUnivariateSpline或make_lsq_spline。有没有一种简单的方法不仅可以得到得到的样条曲线,还可以得到样条曲线系数的方差-协方差矩阵?有没有一种简单的方法来获得在这些函数中生成的用于LSQ问题的正规矩阵?
没有支持的方式,没有
你当然可以从make_interp_spline
:的源代码中获得你需要的东西
https://github.com/scipy/scipy/blob/v1.5.2/scipy/interpolate/_bsplines.py#L854-2009年1月
(正规的eq矩阵是用_bspl._norm_eq_lsq
构造的(。但请注意,内部结构可能会发生变化,恕不另行通知。
特别是,正态方程的使用可能会消失。
对于LSQUnivariateSpline,您需要阅读FITPACK的Fortran源代码。(有点(好消息是,它在未来的科幻版本中不太可能改变。