2021-06-20

第五小组单叶双曲面

clc
clear  all
hold on
title("第五小组实验：单叶双曲面")
axis([-15 15 -15 15 -15 15])
axis on
grid on
view([0.4178 0.2247 0.8803])
a=3;
b=4;
c=5;
x=-a-5:a+5;
y=-b-5:b+5;
[X,Y]=meshgrid(x,y);
Z=zeros(size(X))+10;
for i=1:80
h=mesh(X,Y,Z);
k=sqrt(1+(Z(1)).^2/c.^2);
Thetal=linspace(0,2*pi,50);
X1=a*k.*cos(Thetal);
Y1=b*k.*sin(Thetal);
Z1=zeros(length(X1))+Z(1);
g=plot3(X1,Y1,Z1,'b');
Z=Z-1/4;
F=getframe;
delete(h);
end
movie(F,2,30)
hold on