close all;clear all;
%% 椭球体、椭圆锥体、单叶双曲面、双叶双曲面
[X,Y] = meshgrid(-5:0.1:5,-5:0.1:5);
figure()
Z = 2*sqrt(25 - X.^2 - Y.^2).*(X.^2 + Y.^2<=25);
waterfall(X,Y,Z);
hold on;
Z = 2*-sqrt(25 - X.^2 - Y.^2).*(X.^2 + Y.^2<=25);
waterfall(X,Y,Z)
hold off;
hidden off;
xlabel('x');
ylabel('y');
zlabel('z');
title('椭球体:x^2+y^2+z^2/2^2=25');
figure()
Z = sqrt(X.^2 + Y.^2).*(X.^2 + Y.^2<=25) + 5*(X.^2 + Y.^2>25);
Z2 = -sqrt(X.^2 + Y.^2).*(X.^2 + Y.^2<=25) -5*(X.^2 + Y.^2>25);
waterfall(X,Y,Z);
hold on;
waterfall(X,Y,Z2);
hold off;
colormap('winter');
hidden off;
xlabel('x');
ylabel('y');
zlabel('z');
title('椭圆锥体:z^2=x^2+y^2,z<=25');
figure()
Z = sqrt(X.^2 + Y.^2-4).*(X.^2 + Y.^2>=4).*(X.^2 + Y.^2<=16) + 4*(X.^2 + Y