我经常在博士论文中使用scipy.optimize.leastsq()
,但是我不知道如何从leastsq()
返回的数据中得到雅可比矩阵的估计值我需要知道在最小化中使用的雅可比矩阵的估计,并与有限差分近似进行比较。
有没有人有一个公式如何得到它?
当你检查如何计算协方差矩阵时,这可能有点棘手。
最佳答案
我认为答案是没有人能从optimize.leatsq中恢复Jacobian。
最新威胁http://mail.scipy.org/pipermail/scipy-user/2011-August/030320.html
从leatsq返回的协方差矩阵的一些例子中,我认为在许多情况下,使用默认设置的精度不是很高。
要检查它有多好,可以将返回的协方差矩阵与雅可比矩阵的外积进行比较。