三维立方体剖分绘图（matlab）

有限元或者其他数值计算程序中经常会涉及到对求解空间的剖分（均匀剖分或者非均匀剖分），在数据可视化过程中，需要绘制剖分单元，反映分布结果。下面给出matlab程序。

function [x,y,z] = plotcube( p1, cube_x, cube_y, cube_z )
%% p1: 立方体左下角顶点坐标
%% cube_x, cube_y, cube_z 是立方体三个方向边长
v = zeros(8,3);
v(1,:)=p1;
v(2,:)=p1+[0 0 cube_z];
v(3,:)=p1+[0 cube_y 0];
v(4,:)=p1+[0 cube_y cube_z];
v(5,:)=p1+[cube_x 0 0];
v(6,:)=p1+[cube_x 0 cube_z];
v(7,:)=p1+[cube_x cube_y 0];
v(8,:)=p1+[cube_x cube_y cube_z];

x=[v(1,1) v(1,1) v(6,1) v(8,1) v(2,1) v(1,1);
   v(2,1) v(2,1) v(5,1) v(7,1) v(4,1) v(3,1);
   v(4,1) v(6,1) v(7,1) v(3,1) v(8,1) v(7,1);
   v(3,1) v(5,1) v(8,1) v(4,1) v(6,1) v(5,1)];

y=[v(1,2) v(1,2) v(6,2) v(8,2) v(2,2) v(1,2);
   v(2,2) v(2,2) v(5,2) v(7,2) v(4,2) v(3,2);
   v(4,2) v(6,2) v(7,2) v(3,2) v(8,2) v(7,2);
   v(3,2) v(5,2) v(8,2) v(4,2) v(6,2) v(5,2)];

z=[v(1,3) v(1,3) v(6,3) v(8,3) v(2,3) v(1,3);
   v(2,3) v(2,3) v(5,3) v(7,3) v(4,3) v(3,3);
   v(4,3) v(6,3) v(7,3) v(3,3) v(8,3) v(7,3);
   v(3,3) v(5,3) v(8,3) v(4,3) v(6,3) v(5,3)];

end

测试程序：

p1=[0 0 0];%%%%%%%左下顶点
cube_x=1;
cube_y=1;
cube_z=1;

p2=[1 1 1];
p3=[2 2 2];
[x1,y1,z1] = plotcube( p1, cube_x, cube_y, cube_z );
[x2,y2,z2] = plotcube( p2, cube_x, cube_y, cube_z );
[x3,y3,z3] = plotcube( p3, cube_x, cube_y, cube_z );
 x=[x1 x2 x3];
 y=[y1 y2 y3];
 z=[z1 z2 z3];

patch(x,y,z,[1;1;1;1;1;1;0.5;0.5;0.5;0.5;0.5;0.5;0;0;0;0;0;0]);
view(3);
grid on;
axis([-5 5 -5 5 -5 5]);
xlabel('X');
ylabel('Y');
zlabel('Z');

结果显示如下：