MATLAB---CAD绘制单位球面

function DrawUnitSphere()
%采用球面参数方程绘制单位球面

M = 51;%在u参数的区域[0,2*pi]上取51个点
N = 26;%在v参数的区域[-pi/2,pi/2]上取26个点

stepu = 2*pi/(M-1);%u参数方向的步长,也就是两个邻近参数之间的间隔
stepv = pi/(N-1);%v参数方向的步长,也就是两个邻近参数之间的间隔

for i=1:N%这是i方向,就是v向,也就是参数域的纵向
    v = -pi/2 + (i-1)*stepv;
    for j=1:M%这是j方向,就是u向,也就是参数域的横向
        u = (j-1)*stepu;
        x(i,j) = cos(u)*cos(v);
        y(i,j) = sin(u)*cos(v);
        z(i,j) = sin(v);
    end
end
mesh(x,y,z,'edgecolor','k','facecolor','w')
axis equal