用Julia和PyPlot绘制球体

Question

最近，我试图用pyplot/julia绘制一个球体，不幸的是，它比我想象的要难。可能积分生成有问题，但我不明白为什么我的实现不起作用。尽管使用原始python代码一切都很好。

我已经尝试将matplotlib曲面图doc中的demo2改编为MWE：

using PyPlot
u = linspace(0,2*π,100);
v = linspace(0,π,100);

x = cos(u).*sin(v);
y = sin(u).*sin(v);
z = cos(v);

surf(x,y,z)

那么，我的Julia实现到底出了什么问题？

Answer 1

X、Y和Z应该是矩阵，而不是向量--否则你只能在球面上画一条曲线，而不是曲面本身。

using PyPlot
n = 100
u = linspace(0,2*π,n);
v = linspace(0,π,n);

x = cos(u) * sin(v)';
y = sin(u) * sin(v)';
z = ones(n) * cos(v)';

# The rstride and cstride arguments default to 10
surf(x,y,z, rstride=4, cstride=4)

最初绘制的曲线对应于那些矩阵的对角线。

plot( diag(x), diag(y), diag(z), color="yellow", linewidth=3 )