3d绘图表面2x ^ 2 + 3y ^ 2 + z ^ 2 = 6

Question

我试图在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)会导致负6-2*x.^2-3*y.*2，但我不知道如何解决这个问题，因为我无法修剪{{}的任何一部分1}}或x。任何人都可以帮忙吗？感谢

Answer 1

这取决于你想要对z的非实际值做什么。 您可以做的一件事是将所有这些值设置为零或NaN（根据@ hbaderts＆＃39;评论）：

z = sqrt(6-2*x.^2-3*y.^2);
z( imag(z)~=0 ) = NaN;

还有一件事：您的代码可能存在问题，因为z是一个长度为1000的向量，您希望它是1000x1000矩阵。您应该使用meshgrid()和x上的y来获取二维矩阵：

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'更好的可视化建议。）

在八度音程上运行上面的代码给出了这个图：