我有一组数据，x 和 y 值，具有 n 个拟合参数。我有一个已知的函数，并且以某种方式涉及 n 个拟合参数。有没有办法对函数进行编码，使其接受任意数量的拟合参数？
我正在尝试编写 levenberg-marquardt 算法以使参数适合数据。我如何使用 numpy 或 scipy 计算关于拟合参数的导数的 hessian 和 jacobian 矩阵？
我试过类似的东西
def func(X,P):
return P[0]+P[1]*X+P[2]*X**2
但是我不知道如何让 scipy 计算每个 P[i] 的梯度