利用元胞自动机建立沙漏模型的matlab代码

以下是一个简单的利用元胞自动机建立沙漏模型的 Matlab 代码：

% 设定模型参数
N = 51;         % 网格尺寸
T = 100;        % 迭代次数
p = 0.5;        % 概率参数
s = floor(N/2); % 沙漏大小

% 初始化沙漏模型
grid = zeros(N);
for i = 1:N
    for j = 1:N
        if (i+j >= s+1 && i+j <= 3*s+1 && i-j >= s-1 && j-i >= 1-s)
            grid(i,j) = 1;
        end
    end
end

% 迭代更新沙漏模型
for t = 1:T
    new_grid = zeros(N);
    for i = 2:N-1
        for j = 2:N-1
            if (grid(i,j) == 1)
                if (rand < p)
                    new_grid(i,j) = 1;
                else
                    new_grid(i,j) = 0;
                end
            else
                if (grid(i-1,j)+grid(i+1,j)+grid(i,j-1)+grid(i,j+1) >= 2)
                    new_grid(i,j) = 1;
                else
                    new_grid(i,j) = 0;
                end
            end
        end
    end
    grid = new_grid;
end

% 输出沙漏模型
imagesc(grid);
colormap gray;
axis equal;
axis off;

运行代码后，会生成一个大小为 51x51 的沙漏模型，如下图所示：


其中黑色表示沙漏内部，白色表示沙漏外部。可以通过修改模型参数来改变沙漏大小、迭代次数和概率参数等。