“单个 trellis 上这一步是精确 MAP 边缘化;两个 trellis 经 interleaver 连接后形成有圈图,轮流运行 SISO 相当于按特定 schedule 使用置信传播,整…”
形式陈述 ​
固定二元线性码的Tanner 图,并设码字通过二元输入的离散无记忆信道独立发送。对输出
作为信道对数似然比。置信传播在有向边上维护两类消息。变量到校验的更新为
偶校验到变量的更新为
一次完整迭代交替完成两类更新。变量的暂定后验为
据其符号作硬判决,并可在所有校验满足时提前停止。每条外发消息都排除收件邻居提供的旧消息;否则同一证据会在一步内被重复计算。
在有限无圈 factor graph 上,sum-product 给出的边缘分布是精确的。随机 LDPC 图在块长趋大且迭代轮数固定时,局部计算邻域以高概率无圈,因此可由树递推分析。对一张有限有圈图运行的 loopy BP 则是近似算法:更新式仍有定义,但不附带一般的收敛、唯一不动点或 MAP 最优保证。
直觉
变量节点把自己的信道证据与“其他校验怎样看我”相加,再把尚未包含目标校验的信息发过去。校验节点则回答:“若其余变量的软判断可信,为使奇偶和为零,这个变量应偏向哪一值?”正负号传播奇偶关系,绝对值表达置信强度。弱证据经过一致的多条路径可累积,冲突证据会互相抵消。
这与连续消除译码的固定串行条件化不同。BP 保留软消息并可并行或分层反复更新,早期没有必须永久接受的单个硬判决;SC 则按预定次序把先前估计代入后续 bit-channel。两者都利用图上的概率分解,但信息状态和调度语义并不相同。
例子与边界
考虑一个无圈的三元偶校验
若
可得两个双曲正切分别是
若
在 BEC 上,消息只需取“已知零、已知一、擦除”三种状态,BP 退化为 peeling:某校验只剩一个未知变量时即可恢复它。这个简化不能原样搬到 BSC 或软输出信道。有限图上的四圈会使消息很快携带自己的旧证据;trapping set 可能让所有 LLR 看似稳定却仍有错误。数值实现还要处理
推论与应用
若图度有界,一轮 BP 访问每条边常数次,成本为
BP 还是 Turbo 迭代译码、隐马尔可夫模型与一般 factor graph 推断的共同语言;在每个场景中,必须分别说明图是否有圈、局部因子是什么、消息是否归一化以及停止条件,而不能仅以“传递置信度”替代算法定义。
参考资料
- Frank R. Kschischang, Brendan J. Frey, and Hans-Andrea Loeliger, “Factor Graphs and the Sum-Product Algorithm,” IEEE Transactions on Information Theory 47(2), 2001, 498–519.
- Tom Richardson and Rüdiger Urbanke, Modern Coding Theory, Cambridge University Press, 2008, Chs. 3–4.
- Niclas Wiberg, Codes and Decoding on General Graphs, Linköping University, 1996.