[论文解读] Probabilistic initial value problem for cellular automaton rule 172

Henryk Fukś|arXiv (Cornell University)|Jul 7, 2010

Cellular Automata and Applications参考文献 8被引用 6

一句话总结

本文提出一种组合方法，通过分析长度为3的块的原像树的拓扑结构，求解元胞自动机规则172的概率初值问题。该方法推导出时间n时块概率的精确表达式，涉及斐波那契数与卢卡斯数，表明1的密度表现出指数衰减的阻尼振荡，且其闭式解涉及这些数列。

ABSTRACT

We consider the problem of computing a response curve for binary cellular automata -- that is, the curve describing the dependence of the density of ones after many iterations of the rule on the initial density of ones. We demonstrate how this problem could be approached using rule 130 as an example. For this rule, preimage sets of finite strings exhibit recognizable patterns, and it is therefore possible to compute both cardinalities of preimages of certain finite strings and probabilities of occurrence of these strings in a configuration obtained by iterating a random initial configuration $n$ times. Response curves can be rigorously calculated in both one- and two-dimensional versions of CA rule 130. We also discuss a special case of totally disordered initial configurations, that is, random configurations where the density of ones and zeros are equal to 1/2.

研究动机与目标

• 在初始测度为p=1/2的随机伯努利测度下，求解元胞自动机规则172的概率初值问题。
• 推导出规则172经过n次迭代后，所有长度为3的块的概率的精确表达式。
• 通过分析长度为3的有限字符串的原像树的拓扑结构，组合地枚举n步原像。
• 证明1的密度在配置中随时间演化时表现出阻尼振荡，这是由于斐波那契数与卢卡斯数项的相互作用所致。
• 展示该方法在具有规则原像树结构的其他初等元胞自动机规则中的适用性。

提出的方法

• 该方法为规则172下所有长度为3的块构建原像树，通过分析其层级集来统计n步原像的数量。
• 采用经典的组合枚举技术，计算每个3位块在时间n时原像集的基数。
• 块的概率由原像计数结果按4^n（n步后可能的3位块总数）进行缩放后得出。
• 该方法利用斐波那契数与卢卡斯数的闭式表达式，以精确代数形式表达概率与密度。
• 该方法依赖于原像树的规则结构，通过与斐波那契和卢卡斯数列相关的递推关系实现精确枚举。
• 该方法可推广至任意初始测度μ_p，通过马尔可夫链进行概率计算，尽管此部分工作留待未来完成。

实验结果

研究问题

• RQ1 从p=1/2的随机初始配置出发，经规则172迭代n次后，每个3位块的精确出现概率是多少？
• RQ2 在规则172下，1的密度如何随时间演化，其表现出何种动力学特征？
• RQ3 是否可以利用初等元胞自动机原像树的结构，推导出块概率的精确解析表达式？
• RQ4 为何斐波那契数与卢卡斯数会出现在规则172块概率的精确公式中？
• RQ5 该方法在多大程度上可推广至具有规则原像树拓扑结构的其他元胞自动机规则？

主要发现

• 规则172在n次迭代后1的密度为5/8 − 2^{−n−2}F_{n+3} − 2^{−n−4}F_{n+2}，表现出在指数衰减之上的阻尼振荡。
• 概率P_n(001) = P_n(100) = 1/8对所有n保持不变，表明其在动力学下具有不变性。
• 概率P_n(011) = P_n(110) = 2^{−n−4}L_{n+2}涉及卢卡斯数L_{n+2}，将数论序列与CA动力学联系起来。
• 概率P_n(101) = 2^{−n−4}F_{n+2}与P_n(111) = 2^{−n−3}F_{n+3}明确依赖于斐波那契数。
• 块000的n步原像集基数为5·4^n − 2^{n+1}F_{n+3} − 2^{n−1}F_{n+2}，通过组合枚举推导得出。
• 该方法成功推广至任意初始测度μ_p，使用马尔可夫链技术，尽管此扩展在当前工作中未完全计算。