MATLAB | 趣味编程——三维彭罗斯三角形可视化

光整复杂画图和算法大家看着也累，这里时不时发点有趣的小程序，这期是绘制一个三维的彭罗斯三角形:

绘制效果

完整代码

就是画几个立方体拼在一起，其中一个立方体削掉一块，啥也不说了直接上代码，反正就几行：

function penroseTriangle % @author : slandarer
vertex=[0 0 0;0 0 1;0 1 0;0 1 1;1 0 0;1 0 1;1 1 0;1 1 1]; facet=[1 2 4 3;1 2 6 5;1 3 7 5;2 4 8 6;3 4 8 7;5 6 8 7]; faceColor=[.9,.9,.9]; for i=1:5,patch('Vertices',vertex+[i,3,0],'Faces',facet,'FaceColor',faceColor);end for i=1:4,patch('Vertices',vertex+[1,3,i],'Faces',facet,'FaceColor',faceColor);end for i=1:2,patch('Vertices',vertex+[1,3-i,4],'Faces',facet,'FaceColor',faceColor);end
vertex=[0 1 0;0 0 1;0 1 0;0 1 1;1 1 0;1 0 1;1 1 0;1 1 1]; patch('Vertices',vertex+[1,0,4],'Faces',facet,'FaceColor',faceColor) ax=gca;hold on;axis off;light ax.XLim=[0,6]; ax.YLim=[0,6]; ax.ZLim=[0,6]; ax.DataAspectRatio=[1 1 1]; set(gcf,'Color',[1,1,1]) view(-45.1099,35.3223) end

至于那个视角-45.1099,35.3223数值是咋设定的，其实就是鼠标转动调整到合适视角后get(gca,'View')获取的比较合适的数值，以下展示不同视角的图像：