plotcylinder matlab,Matlab在任意两点之间绘制三维圆柱

Matlab在任意两点之间绘制三维圆柱

Matlab在任意两点之间绘制三维圆柱

此函数可能存在一些不足，请多多指教！

function plotcylinder(u1,u2,color_a,r)

L=norm(u1-u2);

ROD=u2-u1;

[X,Y,Z]=cylinder(r,100);

x1=X*0;y1=Y*0;z1=Z*0;

Z=L*Z-L/2;

ROD_midpoint=(u1+u2)/2;

x=ROD_midpoint(1);

y=ROD_midpoint(2);

z=ROD_midpoint(3);

a=[1 0 0];b=[0 1 0];c=[0 0 1];

if(ROD(2)==0||ROD(1)==0)

if(ROD(2)==0) % 在XZ平面

angel=acos(dot(ROD,c)/norm(ROD)/norm(c));

if(ROD(1)<0)

angel=-angel; %%

end

A2=[cos(angel) 0 sin(angel);0 1 0; -sin(angel) 0 cos(angel)]; % 绕Y轴旋转

for i=1:length(X(1,:))

u=[X(1,i) Y(1,i) Z(1,i)]';

u1=A2*u;

x1(1,i)=u1(1);y1(1,i)=u1(2);z1(1,i)=u1(3);

u=[X(1,i) Y(1,i) Z(2,i)]';

u1=A2*u;

x1(2,i)=u1(1);y1(2,i)=u1(2);z1(2,i)=u1(3);

end

end

if(ROD(1)==0) % 在YZ平面

angel=acos(dot(ROD,c)/norm(ROD)/norm(c));

angel=2*pi-angel; % 此处maybe要增加类似上面的判定是否变为负数

A1=[1 0 0;0 cos(angel) -sin(angel);0 sin(angel) cos(angel)]; % 绕X轴旋转

for i=1:length(X(1,:))

u=[X(1,i) Y(1,i) Z(1,i)]';

u1=A1*u;

x1(1,i)=u1(1);y1(1,i)=u1(2);z1(1,i)=u1(3);

u=[X(1,i) Y(1,i) Z(2,i)]';

u1=A1*u;

x1(2,i)=u1(1);y1(2,i)=u1(2);z1(2,i)=u1(3);

end

end

else

% 先绕Z轴旋转，在绕Y轴旋转

angel=acos(dot(ROD,c)/norm(ROD)/norm(c));

A2=[cos(angel) 0 sin(angel);0 1 0; -sin(angel) 0 cos(angel)];

angel=acos(dot(ROD,a)/norm(ROD)/norm(a));

if(ROD(2)<0)

angel=2*pi-angel;

end

A3=[cos(angel) -sin(angel) 0;sin(angel) cos(angel) 0;0 0 1];

A=A3*A2;

for i=1:length(X(1,:))

u=[X(1,i) Y(1,i) Z(1,i)]';

u1=A*u;

x1(1,i)=u1(1);y1(1,i)=u1(2);z1(1,i)=u1(3);

u=[X(1,i) Y(1,i) Z(2,i)]';

u1=A*u;

x1(2,i)=u1(1);y1(2,i)=u1(2);z1(2,i)=u1(3);

end

end

fill3(x1(1,:)+x,y1(1,:)+y,z1(1,:)+z,color_a,'EdgeColor','none')

hold on

fill3(x1(2,:)+x,y1(2,:)+y,z1(2,:)+z,color_a,'EdgeColor','none')

hold on

surf(x1+x,y1+y,z1+z,'facecolor',color_a,'edgecolor','none')

axis equal

xlabel('x')

ylabel('y')

zlabel('z')

set(gcf,'color','w')

主函数如下：

figure(1)

set(gcf,'color','w')

a=[2 4 6];

b=[10 10 8];

plotcylinder(a,b,'b',0.1)

xlabel('x')

ylabel('y')

zlabel('z')

axis equal

camlight;

lighting flat;

效果图

Matlab在任意两点之间绘制三维圆柱相关教程