有限元或者其他数值计算程序中经常会涉及到对求解空间的剖分（均匀剖分或者非均匀剖分），在数据可视化过程中，需要绘制剖分单元，反映分布结果。下面给出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');结果显示如下：
