Matlab中的这个3D图有什么问题

Question

我正在尝试在 Matlab 中创建 3D 绘图，但我几乎没有经验。我真的很想绘制这些方程式所描述的数字：

x = cos(u)*(2 + cos(v))
y = sin(u)*(2 + cos(v))
z = x^2 - y^2 + 2*x*y*tan^2(v)
-pi <= u <= pi, -pi <= v <= pi

哪个情节本来应该是这样的

_{（来源：paulbourke.net）}

我已经尽力了并想出了：

% volume data
u = linspace(-pi,pi,40);
v = linspace(-pi,pi,40);
[u, v] = meshgrid(u, v); 
X =cos(u).*(2 + cos(v));
Y =sin(u).*(2 + cos(v));
Z =(X.^2)-(Y.^2)+(X.*Y.*2.*((tan(v)).^2));
surf(X,Y,Z,'FaceColor','none','EdgeColor','interp')
 
%shaded surface
colormap(jet)
camlight right

我得到以下信息：

怎么了？

Answer 1

您的图表实际上没有任何问题。您之所以看到这样的图表是因为z 轴太大。实际定义其大部分形状的z 值比该图中看到的最高值要小几个数量级……大约数千。原因是当u 和v 的值是+/- pi/2 时，tan 是未定义的，这就是为什么当你倾向于这些角度时高度如此之高。

如果z 值太大，我的一个建议是剪裁它们。在z 值上设置某种阈值并裁剪它们。 100 表示高正值，-100 表示大负值。此外，调整视角以获得更好的视图，我还建议增加点数，以便更好地插入阴影......类似于 1000。要获得更光滑的效果，请尝试将 FaceColor 属性更改为 @ 987654331@ 而不是none。

这样：

u = linspace(-pi,pi,1000); %// Define 1000 points per dimension
v = linspace(-pi,pi,1000);
[u, v] = meshgrid(u, v); 
X =cos(u).*(2 + cos(v));
Y =sin(u).*(2 + cos(v));
Z =(X.^2)-(Y.^2)+(2*X.*Y.*(tan(v).^2));
Z(Z <= -100) = -100; %// Enforce threshold
Z(Z >= 100) = 100;
surf(X,Y,Z,'FaceColor','interp','EdgeColor','interp')

%shaded surface
colormap(jet)
view(30,50); %// Change viewing angle
camlight right

我明白了：