谢尔平斯基地毯维数代码

谢尔平斯基地毯是一种分形图形，其维数为1.585。下面是一个简单的 Python 代码实现：

import numpy as np
import matplotlib.pyplot as plt

def sierpinski_carpet(n):
    if n == 0:
        return np.array([[1]])
    else:
        prev_carpet = sierpinski_carpet(n-1)
        prev_size = prev_carpet.shape[0]
        new_size = 3 * prev_size
        new_carpet = np.zeros((new_size, new_size))
        for i in range(3):
            for j in range(3):
                if i == 1 and j == 1:
                    new_carpet[prev_size:2*prev_size, prev_size:2*prev_size] = prev_carpet
                else:
                    new_carpet[i*prev_size:(i+1)*prev_size, j*prev_size:(j+1)*prev_size] = prev_carpet
        return new_carpet

carpet = sierpinski_carpet(4)
plt.imshow(carpet, cmap='gray')
plt.axis('off')
plt.show()

这个代码使用递归的方式生成谢尔平斯基地毯，最终输出一个 numpy 数组，可以使用 matplotlib 库将其可视化。