维特比译码介绍liuchao.ppt

The Viterbi Algorithm 刘超 liuchao@hziee.edu.cn 杭州电子科技大学通信学院网络通信教研室 教学内容： 卷积码的简要介绍 维特比译码的基本原理 维特比译码的基本过程 教学目标 掌握维特比译码的基本原理 熟悉用栅格描述维特比译码的过程 卷积码编码器 卷积码编码器结构框图 [例1] (2,1,2)码的状态向量为S=(σ2σ1)，共有4种状态S0=(0,0),S1=(0,1),S2=(1,0),S3=(1,1)，如图所示。 该码的状态转移方程和输出方程分别为 σ1’=U σ2’=σ1 V1=U +σ1+σ2 V2=U +σ2 卷积码的状态转移图 卷积码的栅格图(篱笆图) 状态图不能反映出状态转移与时间的关系 栅格图/篱笆图：将开放型的状态转移图按时间顺序级联形成一个栅格图。 编码路径：状态序列σ在栅格图中形成的一条有向路径。 当有向路径始于全“0”状态S0，又终于S0时，表明此时编码器又回到全“0”状态， 红实线表示U=0时输入产生的转移分支; 黄虚线表示U=1时输入产生的转移分支; 转移分支上数字表示输出的编码比特V1和V2。 最大似然译码/最小距离译码 待编码的信息序列M：M=[M0, M1,…, ML－1]； 编码器输入序列的总长度：k(L+m)； 编码器输出的码序列C：C=[C0, C1,…,CL－1]，其中每个子码Ci含有n个比特； 经离散无记忆信道(DMC)传输后，译码器接收的序列 R：R=[R0, R1,…,RL－1]； 对于DMC信道： 码序列 C 的路径度量 M(R/C)：计算第 l 时刻到达状态 i 的最大似然路径的相似度—log p(R/C)； 子码 Ci 度量M(Ri/Ci) ：计算第 l 时刻接收子码 Ri 相对于各码字的相似度— log p(Ri/Ci)，也称为分支度量。 最大似然译码/最小距离译码 译码器接收到 R 序列后，按最大似然法则力图寻找编码器在篱笆图上原来走过的路径，也就是寻找具有最大度量的路径； 对BSC信道，就是寻找与 R 有最小汉明距离的路径，即计算和寻找 min[d(R, Cj)]，j=1,2,…,2Lk 。 注：二进制对称信道BSC（Binary Symmetry Channel） 最大似然译码/最小距离译码 最大似然译码方法只是提供了一个译码准则，实现起来尚有一定困难。因为它是考虑了长度为 (L+m)n 的接收序列来译码的，这样的序列可能有 2Lk 条； 若实际接收序列中，L=50，k=2，则可能的路径有 2100 条。译码器每接收一个序列 R，就要计算 1030 个似然函数才能做出译码判决。若 kL 再大一些，译码器按最大似然译码准则译码将是很困难的。 维特比译码工作原理 维特比提出了一种算法：译码器不是在篱笆图上一次就计算和比较 2Lk 条路径，而是接收一段，就计算、比较一段，从而在每个状态时，选择进入该状态的最可能的分支。 维特比译码的基本思想：将接收序列 R 与篱笆图上的路径逐分支地比较，比较的长度一般取 (5~6)mn，然后留下与 R 距离最小的路径，称为幸存路径，而去掉其余可能的路径，并将这些幸存路径逐分支地延长并存储起来。 幸存路径的数目等于状态数：2km 以 (2,1,2) 卷积码为例说明维特比译码的一般过程： 设发送序列 C 为全0； 接收序列 R=[10,00,01,00,00,00,00,…] 假设译码器的初始状态为全0； 第0个时刻：接收序列的第0个分支 R0=10 进入译码器。从 S0 状态有两个分支，它们是 00 和 11，R0与这两个分支比较，比较的结果和到达的状态如表1 所示： 每个状态/节点都有两个存储器： 路径存储器：存储该状态的部分路径； 路径值存储器：存储达到该状态的部分路径值 (累加距离)。 第一个时刻：进入译码器的接收码组 R1=00 和此时刻出发的四条分支比较，比较结果和达到状态如表2所示： 从第一个时刻到第二个时刻：共有四条路径，到达S0, S1, S2和S3。在第二个时刻以前译码器不做任何选择和判决。 每个状态的路径存储器存储下此时刻的幸存路径：0000,0011,1110,1101； 每个状态的路径值存储器存储了此时刻到达该状态的幸存路径累加值 (累加距离)。 维特比译码的基本原理 从第二个时刻起：第二个接收码组 R2=01 进入译码器，从篱笆图上可见，从第二个时刻到第三个时刻，进入每个状态的分支有两个（或者说在第三个时刻，进入每个状态的路径有两条）。译码器将接收码组 R2 与进入每个状态的两个分支进行比较和判决，选择一个累加距离（部分路径值）最小的路径作为进入该状态的幸存路径。这样的