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

Posted 2022-10-05 slandarer

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

绘制效果

---

完整代码

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

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')获取的比较合适的数值，以下展示不同视角的图像：

---