PDSCH物理层处理过程（四）：信道编码（基于LDPC码）

下文为对3GPP TS38.212 Clause5.3部分的学习。

LDPC编码

针对每一个码块进行编码，给定一个码块，其比特序列输入为$c_0,c_1,c_2,...,c_{K-1}$，码块分割部分见PDSCH物理层处理过程第三部分。信道编码后的输出可以表示为$d_0,d_1,d_2,...,d_{N-1}$，其中对于LDPC基图1，$N=66Z_c$；对于LDPC基图2，$N=50Z_c$。$Z_c$的确定见PDSCH物理层处理过程第三部分。

编码流程

1. 根据$Z_c$确定索引$i_{LS}$(下表为 TS 38.212中表格5.3.2-1)；
   举例：若根据码块分割部分计算得到 $Z_c=384$，则根据表格得到 $i_{LS}$为1。
   注：表格中每一行的值分别为2，3，5，7，9，11，13的倍数。
   

2. 信息比特映射

   for k = 2Zc : K-1
   	if c[k] != <NULL>
   		d[k-2Zc] = c[k];
   	else
   		c[k] = 0;
   		d[k-2Zc] = <NULL>;
   	end
   end
   

3. 生成校验比特。校验比特长度为$N+2Z_c-K$，可以表示为$\mathbf{w}=[w_0,w_1,w_2,...,w_{N+2Z_c-K-1}{]}^T$，使得$$\mathbf{H}\times \left[\begin{array}{c}\mathbf{c}\\ \mathbf{w}\end{array}\right]=0,$$其中$\mathbf{c}=[c_0,c_1,...,c_{K-1}{]}^T$；$0$为全为0的列向量。编码在GF(2)中进行。
   注：GF(2)为素数$p=2$的有限域，只有0、1二元以及异或、与运算

   对于LDPC基图1，基矩阵 ${\mathbf{H}}_{BG}$ 的大小为 $46\times 68$ (行数*列数)；对于LDPC基图2，基矩阵 ${\mathbf{H}}_{BG}$ 的大小为 $42\times 52$ 。${\mathbf{H}}_{BG}$ 中的元素在TS38.212的表格5.3.2-2（LDPC基图1）及表格5.3.2-3（LDPC基图2）中给出，即表格中包括的行列索引处的元素值为1，其余没有包括在表格内的行列索引处的元素值为0。

   将${\mathbf{H}}_{BG}$中的每个元素按以下方式替换为 $Z_c\times Z_c$ 矩阵从而得到矩阵$H$：

   • ${\mathbf{H}}_{BG}$ 中元素值为0的位置替换成大小为 $Z_c\times Z_c$元素值全为0的矩阵；

   • ${\mathbf{H}}_{BG}$ 中元素值为1的位置替换成大小为 $Z_c\times Z_c$ 的循环置换矩阵 $\mathbf{I}(P_{i,j})$，其中 $i$ 和 $j$ 为 ${\mathbf{H}}_{BG}$ 矩阵中元素的行列索引值，$\mathbf{I}(P_{i,j})$通过大小为 $Z_c\times Z_c$ 的单位矩阵向右循环移位 $P_{i,j}$ 得到，其中 $P_{i,j}=mod(V_{i,j},Z_c)$，$V_{i,j}$ 的值在TS38.212的表格5.3.2-2（LDPC基图1）及表格5.3.2-3（LDPC基图2）中给出，根据对应的LDPC基图及 $i_{LS}$ 找到 $V_{i,j}$ 的值。

     下图为LDPC基图1的部分内容。
     举例： ${\mathbf{H}}_{BG}$ 中行索引为0，列索引为0的元素值为1，参照编码步骤1例子中 $i_{LS}=1$得到 $V_{i,j}$ 的值为307，且$Z_c=384$，则$P_{i,j}=mod(V_{i,j},Z_c)=307$，因此将该位置的元素替换为单位矩阵向右循环移位$P_{i,j}=307$的大小为 $Z_c\times Z_c$的矩阵；${\mathbf{H}}_{BG}$ 中行索引为0，列索引为4没有包含在表格中，表示该位置元素值为0，因此将该位置的元素替换为大小为 $Z_c\times Z_c$的全为0的矩阵。
     

4. 校验比特映射

   for k = K : N+2Zc-1
   	d[k-2Zc] = w[k-K];
   end
   

LDPC编码实现

参考原文：https://zhuanlan.zhihu.com/p/455310565
5G标准下有两种分别对应于LDPC基图1和LDPC基图2的基矩阵，

其中，$I$ 区域表示该部分整体是一个单位阵，且对角线上元素1对应的子矩阵的循环移位系数全为0，$0$ 区域表示该区域在 $H$ 矩阵中对应的部分全是零，两种基矩阵下 $\mathbf{B}$ 区域均是 $4\times 4$矩阵且对应的循环移位系数分布情况只有 $4$ 种：

注：编码实现中矩阵中的元素，其中 1 表示单位矩阵向左移一位， 0 表示不偏移， -1 表示 0 矩阵。

由于生成校验位时需要满足$\mathbf{H}\times \left[\begin{array}{c}\mathbf{c}\\ \mathbf{w}\end{array}\right]=0,$，且及由于LDPC基图1及LDPC基图2均具有特殊结构，即矩阵$\mathbf{H}$可以表示为$$\mathbf{H}=\left[\begin{array}{ccc}A& B& 0\\ C& D& I\end{array}\right],$$令$\mathbf{w}=\left[\begin{array}{c}{{\mathbf{p}}_1}^T\\ {{\mathbf{p}}_2}^T\end{array}\right]$，其中${\mathbf{p}}_1$长度为4，${\mathbf{p}}_1$长度为$N+2Z_c-K-4$，可以得到$$\left[\begin{array}{ccc}A& B& 0\\ C& D& I\end{array}\right]\times \left[\begin{array}{c}{\mathbf{c}}^T\\ {{\mathbf{p}}_1}^T\\ {{\mathbf{p}}_2}^T\end{array}\right]=0,$$$$⟹\{\begin{array}{c}A{\mathbf{c}}^T+B{{\mathbf{p}}_1}^T=0\\ C{\mathbf{c}}^T+D{{\mathbf{p}}_1}^T+{{\mathbf{p}}_2}^T=0\end{array}$$ $$⟹\{\begin{array}{c}{{\mathbf{p}}_1}^T=B^{-1}A{\mathbf{c}}^T\\ {{\mathbf{p}}_2}^T=C{\mathbf{c}}^T+D{{\mathbf{p}}_1}^T\end{array}$$
直接通过上述式子计算得到校验位$\mathbf{p}1$较难，其中涉及到矩阵求逆，这对软件或者硬件的消耗很大，为了尽量减小这种消耗，5G标准在$H$矩阵的结构上做了特殊处理，基图中矩阵$B$较小且内容简单，可以通过式子$A{\mathbf{c}}^T+B{\mathbf{p}1}^T=0$计算（由于其中均为模二计算，则$A{\mathbf{c}}^T=B{{\mathbf{p}}_1}^T$），通过该方法可以降低编码硬件实现的复杂度，实现快速编码。

以LDPC基图1为例，计算${{\mathbf{p}}_1}^T$方法如下(下图中的$\mathbf{s}$为上述式子中的$\mathbf{c}$)：