半球面：[fia,theta]=meshgrid([linspace(0,pi,100),pi]);x=sin(theta).*cos(fia);y=sin(theta).*sin(fia);z=cos(theta);surf(x,y,z)shading interp椭球面：theta=0:0.1:2*pi;phi=theta';x=cos(phi)*cos(theta);y=co...