我试图在octave/matlab中绘制上面的曲面,但遇到了这个问题。我的代码如下:
x = linspace(-sqrt(3),sqrt(3),1000);
y = linspace(-sqrt(2),sqrt(2),1000);
z = sqrt(6-2*x.^2-3*y.^2);
surf(x,y,z)
我得到了错误:
error: mesh: X, Y, Z arguments must be real.
我理解这是因为一些(x,y)
s会导致负6-2*x.^2-3*y.*2
,但我不知道如何解决这个问题,因为我不能裁剪x
或y
的任何部分。有人能帮上忙吗?谢谢
发布于 2015-05-18 13:58:36
这取决于您想如何处理z
的非实数值。您可以做的一件事是将所有这些值设置为0或NaN
(根据@hbaderts的评论):
z = sqrt(6-2*x.^2-3*y.^2);
z( imag(z)~=0 ) = NaN;
不过还有一件事:您的代码可能会有问题,因为z
是一个length-1000
向量,而您希望它是一个1000x1000
矩阵。您应该在x
和y
上使用meshgrid()
来获取任意位置的二维矩阵:
x = linspace(-sqrt(3),sqrt(3),1000);
y = linspace(-sqrt(2),sqrt(2),1000);
[xx,yy] = meshgrid(x,y);
z = sqrt(6-2*xx.^2-3*yy.^2);
z( imag(z)~=0 ) = NaN;
surf(xx,yy,z,'edgecolor','none')
(感谢@LuisMendo为更好的可视化而提出的'edgecolor','none'
建议。)
在octave上运行上面的代码会得到这样的结果:
https://stackoverflow.com/questions/30295960
复制相似问题