在同一图形上组合绘图和冲浪（matlab）

Question

我想在同一张图上用Matlab在表面图上绘制一些磁盘。 我的问题是计算由小的圆形障碍物（磁盘）散射的波。 该波在网格（X，Y）上计算，并存储为矩阵U。可以使用以下命令对其进行绘制：surf（X，Y，U，'EdgeColor'，'None'，'facecolor'， '的interp'）; 这样可以得到一张漂亮的图片 （我无法在stackoverflow上上传图片）

现在，我要在同一张图片上添加磁盘。 注意，“ O”是磁盘中心的矩阵，“ a”是磁盘半径的向量，而N_scat是磁盘数。

hold on
surf(X,Y, U, 'EdgeColor','None', 'facecolor', 'interp');
theta = [1:361]*2*pi/360;
nb = length(theta);
for p = 1:N_scat
    plot(O(1,p).*ones(1,nb) + a(p).*cos(theta), O(2,p).*ones(1,nb) + a(p).*sin(theta), 'k', 'LineWidth', 1)
end
hold off
view(2);
colorbar;

但是，这部分给出了预期的结果（请参见此处的图片 ）。 确实，磁盘已显示，但是如您所见，大多数磁盘只是部分绘制的...

有人有主意吗？

对于以下情况也会发生相同的结果：-网格而不是冲浪-冲浪中没有选项（经典的平面表示和边缘着色）-磁盘的线宽更大或更小

在此先感谢您，祝您度过愉快的一天！

---

Answer 1

问题是因为圆是在z = 0处绘制的，所以它们在某些点处在表面下

您可以在绘图调用中指定线条的zdata，我建议在曲面的最大高度处绘制所有圆（更高的圆可能会影响颜色图）

在循环外部定义zdata，例如

myzdata=max(U(:))*ones(1,nb);

然后在循环内将其设置为该行的zdata

for p = 1:N_scat
    plot(... what you had before ...,'zdata',myzdata)
end