matlab空间中任意一点上画圆柱

%function [Cylinder EndPlate1 EndPlate2] = cylinder3(X1,X2,r,n,cyl_color,closed,lines) % for test X1=[10 10 10];%输入三维坐标 X2=[10 10 30];%根据显示方向调整坐标值，本例选择以Z坐标为显示方向，即原始坐标值加上高度 r=0.25;%显示圆柱的半径 n=20;%圆柱面由n个矩形组成 cyl_color='b';%显示颜色 closed=1;%选择封闭实体；1为选择空心实体显示 lines=0;%显示曲面，不显示棱边；1为显示棱边 length_cyl=norm(X2-X1); t=linspace(0,2*pi,n)'; x2=r*cos(t); x3=r*sin(t); x1=[0 length_cyl]; xx1=repmat(x1,length(x2),1); xx2=repmat(x2,1,2); xx3=repmat(x3,1,2); if closed==1 hold on EndPlate1=fill3(xx1(:,1),xx2(:,1),xx3(:,1),'r'); EndPlate2=fill3(xx1(:,2),xx2(:,2),xx3(:,2),'r'); end Cylinder=mesh(xx1,xx2,xx3); unit_Vx=[1 0 0]; angle_X1X2=acos( dot( unit_Vx,(X2-X1) )/( norm(unit_Vx)*norm(X2-X1)) )*180/pi; axis_rot=cross([1 0 0],(X2-X1) ); if angle_X1X2~=0 rotate(Cylinder,axis_rot,angle_X1X2,[0 0 0]) if closed==1 rotate(EndPlate1,axis_rot,angle_X1X2,[0 0 0]) rotate(EndPlate2,axis_rot,angle_X1X2,[0 0 0]) end end if closed==1 set(EndPlate1,'XData',get(EndPlate1,'XData')+X1(1)) set(EndPlate1,'YData',get(EndPlate1,'YData')+X1(2)) set(EndPlate1,'ZData',get(EndPlate1,'ZData')+X1(3)) set(EndPlate2,'XData',get(EndPlate2,'XData')+X1(1)) set(EndPlate2,'YData',get(EndPlate2,'YData')+X1(2)) set(EndPlate2,'ZData',get(EndPlate2,'ZData')+X1(3)) end set(Cylinder,'XData',get(Cylinder,'XData')+X1(1)) set(Cylinder,'YData',get(Cylinder,'YData')+X1(2)) set(Cylinder,'ZData',get(Cylinder,'ZData')+X1(3)) set(Cylinder,'FaceColor',cyl_color) if closed==1 set([EndPlate1 EndPlate2],'FaceColor',cyl_color) else EndPlate1=[]; EndPlate2=[]; end if lines==0 set(Cylinder,'EdgeAlpha',0) end axis([0,50,0,50,0,50]);%控制坐标轴范围，分别为X,Y,Z最大和最小值