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