有意思的是,Sierpinski三角形和杨辉三角的01表示是对应的。
基本单位是三角形,对于原三角形,删去各边中点构成的小三角形,不断重复下去。
function draw_Sierpinski_Triangle3D(P1, P2, P3, P4, N, options)
% Sierpinski_Triangle 的基本单位是三棱锥
% P1, P2,P3, P4是三棱锥的四顶点(必须是1*3的向量),N是递归层数
% options:三棱柱表面颜色
default_options = [0 1 0];
if nargin == 5
options = default_options;
end
Color = options;
if N == 1
plane = [[P1; P2; P3] [P1; P2; P4] [P1; P3; P4] [P2; P3; P4]];
index = [1 4 7 10];
fill3(plane(:,index), plane(:,index+1), plane(:,index+2), Color);
else
draw_Sierpinski_Triangle3D(P1, (P1+P2)/2, (P1+P3)/2, (P1+P4)/2, N-1);
draw_Sierpinski_Triangle3D(P2, (P2+P1)/2, (P2+P3)/2, (P2+P4)/2, N-1);
draw_Sierpinski_Triangle3D(P3, (P3+P1)/2, (P3+P2)/2, (P3+P4)/2, N-1);
draw_Sierpinski_Triangle3D(P4, (P4+P1)/2, (P4+P2)/2, (P4+P3)/2, N-1);
end
扩展到三维也是很简单的,然而颜色的渐变不知道怎么处理,图像不好看。
function draw_Sierpinski_Triangle2D(P1, P2, P3, N, options)
% Sierpinski_Triangle 的基本单位是三角形
% P1, P2,P3是三角形的三顶点(必须是1*2的向量),N是递归层数
% options(1):线段宽度
% options(2:4):线段颜色
default_options = [3 0 1 0];
if nargin == 4
options = default_options;
end
LineWidth = options(1);
Color = options(2:4);
if N == 1
plot([P1(1) P2(1) P3(1) P1(1)], [P1(2) P2(2) P3(2) P1(2)], ...
'Color', Color, 'LineWidth', LineWidth, 'EraseMode', 'none');
else
draw_Sierpinski_Triangle2D(P1, (P1+P2)/2, (P1+P3)/2, N-1, options);
draw_Sierpinski_Triangle2D(P2, (P2+P1)/2, (P2+P3)/2, N-1, options);
draw_Sierpinski_Triangle2D(P3, (P3+P1)/2, (P3+P2)/2, N-1, options);
end