MATLAB元胞自动机

元胞自动机执行时需要计算每一个格子周围格子的情况，而这里需要并行执行，不能利用矩阵循环，否则的话，后面的格子参考的是新值，而不是原来的旧的值，这是不行的。
那么如何实现并行，直觉来看应该将矩阵存起来，然后用旧的计算新的，其实不用。可以直接利用矩阵计算：如A(i,j)，可设i和j都是向量，然后就算操作的是整个矩阵

x = 2:n-1;
y = 2:n-1;
sum(x,y) = cells(x,y-1) + cells(x,y+1) + cells(x-1, y) + cells(x+1,y)...
            + cells(x-1,y-1) + cells(x-1,y+1) + cells(x+1,y-1) + cells(x+1,y+1);

代码来源：https://blog.csdn.net/qq_40527086/article/details/86798384