元胞自动机matlab 代码

以下是一个简单的 MATLAB 元胞自动机示例代码，其中使用了一维元胞自动机规则30：

% 元胞自动机规则
rule = [0 0 0 1 1 1 1 0];

% 初始化元胞自动机
n = 100; % 元胞数
t = 50; % 迭代次数
cellular_automaton = zeros(t, n); % 存储元胞自动机状态的矩阵
cellular_automaton(1, round(n/2)) = 1; % 初始化中心元胞为1

% 迭代更新元胞自动机状态
for i = 2:t
    for j = 2:n-1
        % 计算元胞周围三个元胞的状态
        left_cell = cellular_automaton(i-1, j-1);
        center_cell = cellular_automaton(i-1, j);
        right_cell = cellular_automaton(i-1, j+1);
        
        % 计算当前元胞的状态
        if left_cell == 1 && center_cell == 1 && right_cell == 1
            cellular_automaton(i, j) = rule(1);
        elseif left_cell == 1 && center_cell == 1 && right_cell == 0
            cellular_automaton(i, j) = rule(2);
        elseif left_cell == 1 && center_cell == 0 && right_cell == 1
            cellular_automaton(i, j) = rule(3);
        elseif left_cell == 1 && center_cell == 0 && right_cell == 0
            cellular_automaton(i, j) = rule(4);
        elseif left_cell == 0 && center_cell == 1 && right_cell == 1
            cellular_automaton(i, j) = rule(5);
        elseif left_cell == 0 && center_cell == 1 && right_cell == 0
            cellular_automaton(i, j) = rule(6);
        elseif left_cell == 0 && center_cell == 0 && right_cell == 1
            cellular_automaton(i, j) = rule(7);
        elseif left_cell == 0 && center_cell == 0 && right_cell == 0
            cellular_automaton(i, j) = rule(8);
        end
    end
end

% 可视化元胞自动机
imagesc(cellular_automaton);
colormap(gray);

此代码首先定义了一个规则向量，它表示元胞自动机规则30。然后，初始化元胞自动机状态矩阵并设置初始状态为一个中心元胞为1的一维向量。接下来，使用嵌套循环进行迭代更新元胞自动机状态。最后，使用`imagesc`函数可视化元胞自动机状态矩阵。