[求助]请教大家一个简单的元胞自动机程序

[feirfei]

请教大家一个很简单的元胞自动机程序
元胞自动机程序状态演化规则 ①如果一个元胞的状态是0,则当它有状态为1的邻元时,该元胞的状态以概率a×b变为1,并且当它的临域中状态为1的邻元越多,其状态转变为1的概率越大;②如果一个元胞的状态是1,则该元胞的状态保持不变。初始传播者位于中心单元格25×25处

取50×50的网格，临域为Moore型，经过50个仿真时钟,知识传播结果的模拟图见附件




请教大家 此元胞自动机的matlab程序怎么写。 多谢大家了 感激不尽