coneball.m

来自「初学者可以参考一下 作为MATLAB的入门教程与实践」· M 代码 · 共 52 行

clear                          %清空工作空间
R=300;H=400;                   %圆锥半径和高
r=200;                         %圆球半径
h=100;                         %圆球垂直偏移量
a=150;                         %圆球横向偏移量
th=(0:2:362)*pi/180;           %圆周角
th1=(50:5:290)*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;Y1=y1;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));

%画圆球轮廓在圆锥底面的投影
x4=r*cos(th1)+a;
y4=r*sin(th1);
z4=0*ones(size(th1));

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,.7],'EdgeColor','none','FaceLighting','phong')
hold on
plot3(x4',y4',z4','-w')
%hold on
%plot3(x,y,z,'w-')
hold off,axis equal
view(-50,30)
camlight(90,-20)
camlight(-90,-20)
camlight
axis off
pause
view([0,0,1]),
pause,
view([0,-1,0]),
pause,
view([0,0,-1]),
pause
view(70,-10)