在Julia中使用Plot或PyPlot绘制3D曲面
在Julia中使用Plot或PyPlot绘制3D曲面
提问于 2020-12-29 02:31:06
我想绘制一个两个变量的函数(代码中的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)));using Plots
plot(KX,KY,e_pos, st=:surface,cmap="inferno")如果我以这种方式使用绘图，有时我会得到一个没有曲面的空3D平面。我做错了什么？我认为这可能与我为kx和ky所做的网格有关，但我不确定。
编辑:我还得到了以下错误：
回答 1
Stack Overflow用户
发布于 2020-12-30 09:30:03
我在我的代码中修改了一些东西。
首先，我将变量保留为范围。其次，我只是简单地计算了我需要的函数，而没有将变量映射到它们上。代码如下：
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)")我得到了我想要的！
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/65482345
复制相关文章
相似问题
腾讯云开发者
