approx_fprime - ValueError：使用序列设置数组元素

Question

当我执行approx_fprime时，我收到错误（标题）

import numpy as np
import scipy.optimize as op
def J(x):
    return x*x

xk = np.arange(1,10).reshape(3,3).astype('float32')
print(op.approx_fprime(xk.ravel(), J, 0.01))

我做错了什么？

---

我必须改变函数，它应该返回一个标量而不是一个向量。

def J(x): 
    return np.sum(x**2)

Answer 1

你使approx_fprime的假设无效：

标量函数梯度的有限差分逼近。

参数：

...

f：可赎回

其功能是确定梯度（偏导数）。应该将xk作为第一个参数，f的其他参数可以在* args中提供。应该返回一个标量，函数的值在xk。

在这种情况下，您将返回一个大小为9的向量。这不是标量函数应返回的内容。原因很简单：

x = np.array([1,2,3])
print((x*x).shape)
# (3,)