coneball2.m

MATLAB的基本操作与使用

clear                          %清空工作空间
R=300;H=400;                   %圆锥半径和高
r=200;                         %圆球半径
h=100;                         %圆球垂直偏移量
a=150;                         %圆球横向偏移量
th=(0:2:362)*pi/180;           %圆周角
%画圆球
z1=(-r:2:r)';                  %垂向坐标（一维）
x1=sqrt(r^2-z1.^2)*cos(th);    %横向坐标
y1=sqrt(r^2-z1.^2)*sin(th);    %纵向坐标
z1=z1*ones(size(th));          %垂向坐标（二维）
x1=a+x1;z1=z1+h;               %平移
%画圆锥
r2=[0,R]';                     %母线
x2=r2*cos(th);
y2=r2*sin(th);
z2=[H;0]*ones(size(th));
%画圆锥底
r3=[0,R]';
x3=r3*cos(th);
y3=r3*sin(th);
z3=[0;0]*ones(size(th));

%画圆球轮廓在圆锥底面的投影
th1=(65:5:295)*pi/180;         %圆球轮廓在底面的投影转角
x4=r*cos(th1)+a;
y4=r*sin(th1);
z4=0*ones(size(th1));
subplot(221)
surf(x1',y1',z1','FaceColor',[.7,.1,.1]','EdgeColor','none','FaceLighting','phong')
camlight
axis equal,
axis off
subplot(222)
surf(x2',y2',z2','FaceColor',[.1,.7,.1],'EdgeColor','none','FaceLighting','phong')
camlight
axis equal,
axis off
subplot(223)
surf(x1',y1',z1','FaceColor',[.7,.1,.1]','EdgeColor','none','FaceLighting','phong')
hold on
surf(x2',y2',z2','FaceColor',[.1,.7,.1],'EdgeColor','none','FaceLighting','phong')
hold on
surf(x3',y3',z3','FaceColor',[1,1,0],'EdgeColor','none','FaceLighting','phong')
hold on,axis equal,
camlight (45,45,'infinite')
camlight (135,-25,'infinite')
axis off
view(70,-10)

subplot(224)
surf(x1',y1',z1','FaceColor',[.7,.1,.1]','EdgeColor','none','FaceLighting','phong')
hold on
surf(x2',y2',z2','FaceColor',[.1,.7,.1],'EdgeColor','none','FaceLighting','phong')
hold on
surf(x3',y3',z3','FaceColor',[1,1,0],'EdgeColor','none','FaceLighting','phong')
hold on,axis equal,
camlight (135,-25,'infinite')
axis off
plot3(x4',y4',z4','-b')
hold off,
view([0,0,-1]),