[½»Á÷] ¡¾ÇóÖú¡¿ÎªºÎÉú²»³öÈýάÇòÌå ÒÑÓÐ1È˲ÎÓë

%%PSL@CSU
%%QQ:547423688
%%Email:anyuezhiji@qq.com

d=6;%°ë¾¶
C=[0,0,0];%ÇòÐÄ×ø±ê

figure
hold on
fl = @(x,y) -sqrt(d^2-(x-C(1)).^2 - (y-C(2)).^2)+C(3);
ezsurf(fl)
fh = @(x,y) sqrt(d^2-(x-C(1)).^2 - (y-C(2)).^2)+C(3);
ezsurf(fh)
axis([-d d -d d -d d])
View(36,6)
title(['(x-',num2str(C(1)),')^2 + (y-',num2str(C(2)),')^2+(z-',num2str(C(3)),')^2=',num2str(d),'^2'])

figure
[x,y,z] = sphere(36);
X = C(1)+d*x;
Y = C(2)+d*y;
Z = C(3)+d*z;
surf(X,Y,Z ,'FaceColor','yellow')
axis equal
title(['(x-',num2str(C(1)),')^2 + (y-',num2str(C(2)),')^2+(z-',num2str(C(3)),')^2=',num2str(d),'^2'])

figure
syms a b c theta phi;
a=d*cos(theta)*cos(phi)+C(1);
b=d*sin(theta)*cos(phi)+C(2);
c=d*sin(phi)+C(3);
ezsurf(a,b,c,[0,2*pi],[0,2*pi])
title(['(x-',num2str(C(1)),')^2 + (y-',num2str(C(2)),')^2+(z-',num2str(C(3)),')^2=',num2str(d),'^2'])
axis equal

figure
theta=(0:pi/18:2*pi);
phi=(0:pi/18:2*pi);
[theta,phi]=meshgrid(theta,phi)
a=d*cos(theta).*cos(phi)+C(1);
b=d*sin(theta).*cos(phi)+C(2);
c=d*sin(phi)+C(3);
surf(a, b, c ,'FaceColor','yellow','FaceAlpha',0.5)
axis equal