如何从leastsq()返回的数据中获取雅可比矩阵？

Question

我在我的博士论文中经常使用 scipy.optimize.leastsq()，但是我不知道如何从 leastsq() 返回的数据中获取雅可比矩阵的估计值。我需要知道在最小化过程中使用的雅可比矩阵的估计值，以便与在最小值处的有限差分近似进行比较。

有没有人知道怎么计算这个？

这可能有点棘手，比如说在 leastsq() 内部是如何计算协方差矩阵的。

Answer 1

我觉得答案是，没人能从 optimize.leastsq 中恢复出雅可比矩阵。

关于这个问题的最新讨论可以在这里找到：http://mail.scipy.org/pipermail/scipy-user/2011-August/030320.html

从 leastsq 返回的一些协方差矩阵的例子来看，我觉得在默认设置下，精度并不是很高的情况有很多。

要检查结果的好坏，你可以把返回的协方差矩阵和你的雅可比矩阵的外积进行比较。