我想绘制以下二次形式:
f(x) = 1/2 x^T * A * x -b^T *x
其中A是2x2矩阵。 x1轴:[ - 6,4] x2轴:[ - 4,6]
我做了以下事情:
x = arange(-6,4,0.1)
y = arange(-4,6,0.1)
因此x和y都能获得100分。现在我想评估所有这些点对的二次形式,即我得到100 * 100 = 10000点。
通过以下代码完成:
res= []
for a in x:
for b in y:
val = f(array([a,b]))
res.append(val)
现在我想制作一种曲面图。任何人都可以帮我弄清楚如何做到这一点?我试过了:
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
surf = ax.plot_surface(x, y, res)
我认为以下方法应该有所帮助,但我真的不知道具体如何:
griddata-从非均匀间隔网格插值到其他网格。
任何想法?