import matplotlib.pyplot as plt
from matplotlib import cm
from numpy import nan, linspace, meshgrid
x1=linspace(0,2,50)
x2=linspace(0,2,50)
x1, x2 = meshgrid(x1, x2)
f=((x1+1.5)**2+5*(x2-1.7)**2)*((x1-1.4)**2+0.6*(x2-0.5)**2)
f[-x1<=0] = nan
f[-x2<=0] = nan
f[3*x1-x1*x2+4*x2-7<=0] = nan
f[2*x1+x2-3<=0] = nan
f[3*x1-4*x2**2-4*x2<=0] = nan
fig, ax = plt.subplots(subplot_kw={"projection": "3d"})
ax.plot_surface(x1, x2, f, cmap=cm.jet,linewidth=0, antialiased=False,label="Kurva $f(x_1,x_2)$",alpha=1)
ax.set_xlabel('$x_1$')
ax.set_ylabel('$x_2$')
ax.set_zlabel('$f(x_1,x_2)$')
plt.show()

我想用给定的5个约束来绘制f=((x1+1.5)**2+5*(x2-1.7)**2)*((x1-1.4)**2+0.6*(x2-0.5)**2):
-x1<=0-x2<=03*x1-x1*x2+4*x2-7<=02*x1+x2-3<=03*x1-4*x2**2-4*x2<=0当我运行上面的代码时,图中没有显示任何内容。我犯了什么错?怎么修呢?
发布于 2022-01-22 12:45:46
您将x1和x2定义为:
x1=linspace(0,2,50)
x2=linspace(0,2,50)因此,它们的每个值都在0和2之间。然后用以下内容编辑它们:
x1, x2 = meshgrid(x1, x2)此时,x1和x2不再是数组,而是矩阵。无论如何,和以前一样,它们中的每个元素仍然在0和2之间。
如果应用此筛选器:
f[-x1<=0] = nan然后f的每个元素都变成nan,因为x1的所有元素都是正的,所以表达式-x1<=0对x1的每个元素计算为False。同样的逻辑也适用于过滤器:
f[-x2<=0] = nan因此,实际上,f的每个元素都是nan,而图是空的。
如果尝试移除上述两个筛选器,则会得到:
f[3*x1-x1*x2+4*x2-7<=0] = nan
f[2*x1+x2-3<=0] = nan
f[3*x1-4*x2**2-4*x2<=0] = nan

https://stackoverflow.com/questions/70811099
复制相似问题