请扮演matlab元胞自动机代码高手

好的，我会尽力帮助您，以下是一个简单的 MATLAB 元胞自动机代码：

% 设置元胞自动机参数
sizeX = 100; % 网格大小
sizeY = 100;
numSteps = 100; % 迭代次数
density = 0.5; % 初始状态密度

% 初始化网格和随机状态
grid = zeros(sizeX, sizeY);
grid(rand(sizeX, sizeY) < density) = 1;

% 迭代更新网格状态
for i = 1:numSteps
    % 复制上一时刻的网格状态
    newGrid = grid;
    % 遍历所有元胞
    for x = 2:sizeX-1
        for y = 2:sizeY-1
            % 统计周围元胞的数量
            neighbors = sum(sum(grid(x-1:x+1,y-1:y+1))) - grid(x,y);
            % 根据规则更新元胞状态
            if grid(x,y) == 1 && (neighbors < 2 || neighbors > 3)
                newGrid(x,y) = 0;
            elseif grid(x,y) == 0 && neighbors == 3
                newGrid(x,y) = 1;
            end
        end
    end
    % 将新状态赋值给网格
    grid = newGrid;
    % 绘制网格图像
    imagesc(grid);
    colormap(gray);
    drawnow;
end

这个代码实现了一个简单的生命游戏（Conway's Game of Life）元胞自动机。它使用一个二维网格表示元胞状态，每个元胞的状态只有两种可能：活着（1）或死亡（0）。在每一次迭代中，根据周围元胞的状态进行更新，从而模拟元胞自动机的演化过程。最后，代码使用 `imagesc` 函数将网格状态绘制成图像，并在每次迭代后更新。