5G全流程详解4：LDPC码---译码

weixin_43831432

已于 2025-03-26 16:29:54 修改

阅读量623

点赞数 3

分类专栏：无线知识文章标签： 5G

于 2024-11-08 16:21:31 首次发布

本文链接：https://blog.csdn.net/weixin_43831432/article/details/143594153

版权

无线知识专栏收录该内容

8 篇文章

订阅专栏

原理

举个例子来说明，假设校验矩阵是

$A=\begin{Bmatrix} 1 &1 &1 &0 &0 &1 &1 &0 &0 &1 \\ 1 &0 &1 &0 &1 &1 &0 &1 &1 &0 \\ 0 &0 &1 &1 &1 &0 &1 &0 &1 &1 \\ 0 &1 &0 &1 &1 &1 &0 &1 &0 &1 \\ 1 &1 &0 &1 &0 &0 &1 &1 &1 &0 \end{Bmatrix}$

把接收到的码字表示成一个向量：

$r=\begin{Bmatrix} r_{1} &r_{2} &r_{3} &r_{4} &r_{5} &r_{6} &r_{7} &r_{8} &r_{9} &r_{10} \end{Bmatrix}$

根据接收到的码字 $r$ ，可以计算出伴随式syndrome

$s=r\times A^T$

可以展开成5个校验方程

$s_1=r_1+r_2+r_3+r_6+r_7+r_{10}$

$s_2=r_1+r_3+r_5+r_6+r_8+r_9$

$s_3=r_3+r_4+r_5+r_7+r_9+r_{10}$

$s_4=r_2+r_4+r_5+r_6+r_8+r_{10}$

$s_5=r_1+r_2+r_4+r_7+r_8+r_{9}$

举一个具体例子，发送一个码字，接收到的码字有一个错误，用bit-flipping算法做一个译码。

发送码字为：

$c=\begin{Bmatrix} 0 &0 &0 &1 &0 &1 &0 &1 &0 &1 \end{Bmatrix}$

接收到的码字有错误：

$c=\begin{Bmatrix} 0 &0 &0 &1 &1 &1 &0 &1 &0 &1 \end{Bmatrix}$

可以看到 $r_5$ 是错误的，用比特翻转算法，尝试做一个译码：

$s_1=r_1+r_2+r_3+r_6+r_7+r_{10}=0 \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; z_1$

$s_2=r_1+r_3+r_5+r_6+r_8+r_9=1\, \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \, \: \: \: \: z_2$

$s_3=r_3+r_4+r_5+r_7+r_9+r_{10}=1\: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \, \: \: \: \: z_3$

$s_4=r_2+r_4+r_5+r_6+r_8+r_{10}=1\: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \, \: \: \: \: z_4$

$s_5=r_1+r_2+r_4+r_7+r_8+r_{9}=0\: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \, \: \: \: \: z_5$

伴随式不是都为0，查看不等于0的方程里 $r_n$ 出现个数，可以发现只有 $r_5$ 是都出现了，一共出现了3次。

然后我们可以尝试将 $r_5$ 翻转成0，再一次计算伴随式，我们会发现伴随式全部为0，译码结束。

为什么 $r_5$ 我们明明接收错误了，还能够被纠正呢？因为 $r_5$ 的信息不光在 $r_5$ 中，还通过LDPC编码分布在了 $r_n(n!=5)$ 里面，所以我们才可能通过其他的 $r_n$ 得到正确的 $r_5$ 。

我们用 $\lambda$ 来表示对数似然比，很容易可以得到

$\lambda(c_2|r)=log\frac{p(c_2=1|r)}{p(c_2=0|r)}=\frac{2}{\alpha^2 }r_2+\sum_{m}^{}\lambda(z_{m,2}|r_i,i\neq 2)\: \: \: \: \: \: \: \: \: (1)$

在上面公式里 $\frac{2}{\alpha^2 }r_2$ 就是 $r_2$ 本身的一个似然比，而 $\sum_{m}^{}\lambda(Z_{m,2}|r_i,i\neq 2)$ 则是除 $r_2$ 之外的其他接收码字给出的关于 $r_2$ 的信息量。

$Z_{m,2}$ 是指校验方程 $Z_{m}$ 中，除掉比特2之外其他比特给出的关于比特2的信息。

$c_2$ 参与校验的方程有 $z_1$ ， $z_4$ ， $z_5$ 这三个，所以

$\sum_{m}^{}\lambda(Z_{m,2}|r_i,i\neq 2)=\sum_{m\in{\{1,4,5}\}}^{}\lambda(Z_{m,2}|r_i,i\neq 2)=\lambda(z_{1,2}|\{r_i,i\neq 2\})+\lambda(z_{4,2}|\{r_i,i\neq 2\})+\lambda(z_{5,2}|\{r_i,i\neq 2\})$

因为第一个校验方程比特有1，2，3，6，7，10，去掉比特2还有1，3，6，7，10，所以

$\lambda(z_{1,2}|\{r_i,i\neq 2\}) \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: =\lambda(c_1+c_3+c_6+c_7+c_{10}|r_1,r_3,r_4,r_5,r_6,r_7,r_8,r_9,r_{10})=log\frac{p(c_1+c_3+c_6+c_7+c_{10}=1|r_1,r_3,r_4,r_5,r_6,r_7,r_8,r_9,r_{10})}{p(c_1+c_3+c_6+c_7+c_{10}=0)|r_1,r_3,r_4,r_5,r_6,r_7,r_8,r_9,r_{10})}$

下面引入一条定理，具体推导后续更新

$tanh(-\frac{1}{2}\lambda(x_1\bigoplus x_2\bigoplus x_3 \cdots \bigoplus x_n))=tanh(-\frac{\lambda(x_1)}{2})tanh(-\frac{\lambda(x_2)}{2})tanh(-\frac{\lambda(x_3)}{2})\cdots tanh(-\frac{\lambda(x_n)}{2})$

那么有

$\lambda(x_1\bigoplus x_2\bigoplus x_3 \cdots \bigoplus x_n)=-2tanh^{-1}( \prod_{i}^{}tanh(-\frac{\lambda(x_n)}{2}))$

将上面定理应用于 $\lambda(z_{m,n}|{r_i,i\neq n})$ 的计算中，可以得到

$\lambda(z_{m,n}|\{r_i,i\neq n\})=\lambda((\sum_{j\in N_{m,n} }^{}C_j)|\{r_i,i\neq n\})=-2tanh^{-1}(\prod_{j\in N_{m,n}}^{}tanh(-\frac{\lambda(c_j|\{r_i,i\neq n\})}{2})$

上式中 $N_{m,n}$ 表示第 $m$ 个校验方程中除去 $n$ 的所有比特的集合。

令

$\eta_{m,n}= -2tanh^{-1}(\prod_{j\in N_{m,n}}^{}tanh(-\frac{\lambda(c_j|\{r_i,i\neq n\})}{2})\: \: \: \: \: \: \: \: \: \: \: \: (2)$

将上式代入(1)得

$\lambda(c_2|r)=log\frac{p(c_2=1|r)}{p(c_2=0|r)}=\frac{2}{\alpha^2 }r_2+\sum_{m}^{}\eta_{m,n}\: \: \: \: \: \: \: \: \: \: \: (3)$

公式(2)(3)一起构成了迭代循环，我们就可以多次迭代，直至校验方程全部满足。

Matlab代码

function x = ldpc_decode(A,r,Nloop)

[M,N] = size(A);
clear Nl Ml
Nl = cell(M);
Ml = cell(N);
for m=1:M
    Nl{m} = [];
end
for n=1:N
    Ml{n} = [];
end

for m=1:M
    for n=1:N
        if(A(m,n))
            Nl{m} = [Nl{m} n];
            Ml{n} = [Ml{n} m];
        end
    end
end

eta = zeros(M,N);
lasteta = zeros(M,N);
lambda = r;

for loop = 1:Nloop
    fprintf(1,'loop=%d\n',loop);
    for m = 1:M
        for n=Nl{m}
            pr = 1; 
            for np=Nl{m}
                if np == n
                    continue;
                end
                pr = pr*tanh((-lambda(np)+lasteta(m,np))/2);
            end
            eta(m,n) = -2*atanh(pr);
        end
    end
    lasteta = eta;

    for n=1:N
        lambda(n) = r(n);
        for m = Ml{n}
            lambda(n) = lambda(n) + eta(m,n);
        end
    end

    x = lambda >= 0;
    z1 = mod(A*x',2)';

    if(all(z1==0))
        break
    end
end

%X = [1 0 1; 0 1 0; 1 1 0];
%Y = gf(X, 1);
%Z = inv(Y);
%U=X*Z;