使用Plots或PyPlot在Julia中打印三维曲面

我想画一个两变量函数（代码中的e_pos和e_neg）。这里，t和a是常数，我已经给出了1的值

绘制此函数的代码如下所示：

t = 1
a = 1

kx = ky = range(3.14/a, step=0.1, 3.14/a)

# Doing a meshgrid for values of k
KX, KY = kx'.*ones(size(kx)[1]), ky'.*ones(size(ky)[1])


e_pos = +t.*sqrt.((3 .+ (4).*cos.((3)*KX*a/2).*cos.(sqrt(3).*KY.*a/2) .+ (2).*cos.(sqrt(3).*KY.*a)));
e_neg = -t.*sqrt.((3 .+ (4).*cos.((3)*KX*a/2).*cos.(sqrt(3).*KY.*a/2) .+ (2).*cos.(sqrt(3).*KY.*a)));

如果以这种方式使用打印，有时会得到一个没有曲面的空3D平面。我做错了什么？我想这可能与我为kx和ky做的网格有关，但我不确定

编辑：我还得到以下错误：

---

我修改了代码中的一些内容

首先，我将变量保留为范围。其次，我只是计算了我需要的函数，而没有将变量映射到它们上。代码如下：

t = 2.8
a = 1

kx = range(-pi/a,stop = pi/a, length=100)
ky = range(-pi/a,stop = pi/a, length=100)

#e_pos = +t*np.sqrt(3 + 4*np.cos(3*KX*a/2)*np.cos(np.sqrt(3)*KY*a/2) + 2*np.cos(np.sqrt(3)*KY*a))

e_pos(kx,ky) = t*sqrt(3+4cos(3*kx*a/2)*cos(sqrt(3)*ky*a/2) + 2*cos(sqrt(3)*ky*a))
e_neg(kx,ky) = -t*sqrt(3+4cos(3*kx*a/2)*cos(sqrt(3)*ky*a/2) + 2*cos(sqrt(3)*ky*a))

# Sort of broadcasting?
e_posfunc = e_pos.(kx,ky);
e_negfunc = e_neg.(kx,ky);

对于绘图，我只使用GR后端：

using Plots
gr()

plot(kx,ky,e_pos,st=:surface)
plot!(kx,ky,e_neg,st=:surface, xlabel="kx", ylabel="ky",zlabel="E(k)")

我得到了我想要的！

using Plots
gr()

plot(kx,ky,e_pos,st=:surface)
plot!(kx,ky,e_neg,st=:surface, xlabel="kx", ylabel="ky",zlabel="E(k)")