Matlab画三维图像给xy设定范围，定义域可是任意形状

用matlab画曲面图时，比如画一个旋转抛物面,如果如下写程序：

x=-1:0.02:1;
y=x;
[x1,y1]=meshgrid(x,y);%生成画图用的xy网格数据
z=x1.^2+y1.^2;
%C=0.1*z;%设置图像颜色的矩阵，跟z同形状
mesh(x1,y1,z)%画三维网格图
figure
surf(x1,y1,z)%画三维曲面图

图像就会这样：

 

 跟我们想象的旋转抛物面不怎么一样，边界是方的。一种方法是用极坐标生成网格数据，这种方法的局限性较大，生成的xy轴定义域网格数据都是圆形的（能否生成其他形状定义域没有好好研究），下面的方法基本是任意形状的定义域都可以。先上程序和图：

x=-1:0.001:1;
y=x;
[x1,y1]=meshgrid(x,y);%生成画图用的xy网格数据
z=x1.^2+y1.^2;
[m,n]=size(x1);
for i=1:m
    for j=1:m
        if x1(i,j)^2+y1(i,j)^2>1
            x1(i,j)=nan;%将不想要的点横坐标或纵坐标值改为nan
        end
    end
end
figure
mesh(x1,y1,z)%画三维网格图
figure
surf(x1,y1,z)%画三维曲面图
%shading flat;%去掉绘图时网格线的黑边，防止点太密时曲面呈黑色
shading interp%绘图时网格线颜色自动插值，让图像像颜色更自然

因为去掉了surf画图的黑色网格线，所以上面这段程序用surf和mesh画图结果一样。 

为防止边缘出现毛边，该段程序取点时取得比较密集。

原理通过程序大家应该看清楚了，就是把不想要的点横坐标或纵坐标改成nan。比如本例中不想要的点，就将其横坐标都改为了nan，这样三维坐标缺少一维，图像自然就不显示了。

再来个马鞍面：

x=0:0.01:1;
y=x;
[x1,y1]=meshgrid(x,y);
[m,n]=size(x1);
for i=1:m
    for j=1:n
        if abs(x1(i,j))+abs(y1(i,j))>1
            x1(i,j)=nan;
            y1(i,j)=nan;
        end
    end
end
z=x1.*y1;
mesh(x1,y1,z)