%%
圆环面
R=6;r=2;
symsuv;
ezmesh((R+r*cos(u))*cos(v),(R+r*cos(u))*sin(v),r*sin(u));
axisequal;
%%
圆盘
R=6;r=2;
theta=linspace(0,2*pi,90);ph=linspace(r,R,30);
[t,p]=meshgrid(theta,ph);
r=t*0;
[x,y,z]=pol2cart(t,p,r);
mesh(x,y,z);
axisequal;
%%a.
画椭球
a=1;b=2;c=3;
figure(1),title('
椭球
')
ellipsoid(0,0,0,a,b,c)
%
shadinginterp
boxon;
gridoff
colormaphsv
%%b.
单叶双曲面
[x,y]=meshgrid(-2:0.2:2,-2:0.2:2);
z=sqrt(c.^2*(x.^2/(a.^2)+y.^2/(b.^2)-1));
%
figure(2)
subplot(121),title('
单叶双曲面
')
surf(x,y,real(z))
holdon;
surf(x,y,-real(z))
%