sierpinski地毯 matlab代码

下面是一个用Matlab实现Sierpinski地毯的代码示例：

function sierpinski_carpet(n)
    % 创建一个空白的图像
    img = ones(3, 3);

    % 递归绘制Sierpinski地毯
    img = draw_carpet(img, n);

    % 显示图像
    imshow(img);
end

function img = draw_carpet(img, n)
    if n == 0
        return;
    else
        % 计算每个方块的大小
        [h, w] = size(img);
        block_h = h / 3;
        block_w = w / 3;

        % 绘制中心方块
        img(block_h+1:2*block_h, block_w+1:2*block_w) = 0;

        % 递归绘制其他8个方块
        for i = 1:3
            for j = 1:3
                if i ~= 2 || j ~= 2
                    sub_img = img((i-1)*block_h+1:i*block_h, (j-1)*block_w+1:j*block_w);
                    img((i-1)*block_h+1:i*block_h, (j-1)*block_w+1:j*block_w) = draw_carpet(sub_img, n-1);
                end
            end
        end

        % 返回绘制完成的图像
        img = img;
    end
end

使用该代码，您可以调用`sierpinski_carpet(n)`函数来生成Sierpinski地毯，其中`n`是递