谢尔宾斯基地毯matlab

谢尔宾斯基地毯是一种分形图形，可以使用Matlab进行绘制。以下是两种不同的绘制方法：

方法一：

function ditan(n)
a=[0 1 1 0]; 
b=[0 0 1 1]; 
fill(a,b,'b')%将正方形填充蓝色 
for k=1:n%迭代n次 
    for m=1:3^(k-1)%按行 
        for p=1:3^(k-1)%按列 
            fill(1/3^k+1/3^(k-1)*(p-1)+a/3^k,1/3^k+1/3^k+1/3^(k-1)*(m-1)+b/3^k,'w') 
        end 
    end 
end 
hold off 
axis off%不要坐标轴 
axis equal%各坐标轴同比例 
set(findobj(gcf,'type','patch'),'edgecolor','none')%将线条隐去。

方法二：

a=[0 1 1 0]; 
b=[0 0 1 1]; 
n = input('请输入迭代次数：'); 
for k=1:n 
    a=[a(1:end-1) (a(2:end)+a(1:end-1))/2 (a(2:end)+a(1:end-1))/2 a(end)]; 
    b=[b(1:end-1) (b(2:end)+b(1:end-1))/2 b(2:end) b(end)+(b(end)-b(end-1))/2]; 
end 
for k=1:8^n 
    fill(a(4*k-3:4*k), b(4*k-3:4*k), 'b'); 
end 
axis off%不要坐标轴 
axis equal%各坐标轴同比例 
set(findobj(gcf,'type','patch'),'edgecolor','none')%将线条隐去。