#LDPC解码算法
##和积算法
SPA(Sum-Product Algorithm,和积算法) or BPA(Belief Progration Algorithm,置信传播算法)
-
重复码的LLR
c ∈ { 0 , 1 } c\in\{0,1\} c∈{ 0,1}重复编码后经过AWGN信道的接收信号为: r ⃗ = c + n ⃗ \vec{r}=c+\vec{n} r=c+n 则在独立等概条件下,MAP接收机计算出的LLR:
L ( c ∣ r ⃗ ) = log Pr ( c = 0 ∣ r ⃗ ) Pr ( c = 1 ∣ r ⃗ ) = log ∏ l Pr ( r l ∣ c = 0 ) ∏ l Pr ( r l ∣ c = 1 ) = ∑ l log Pr ( r l ∣ c = 0 ) Pr ( r l ∣ c = 1 ) = ∑ l L ( r l ∣ c ) L(c|\vec{r})=\log{\dfrac{\Pr(c=0|\vec{r})}{\Pr(c=1|\vec{r})}}=\log{\dfrac{\prod_l \Pr(r_l|c=0)}{\prod_l \Pr(r_l|c=1)}}=\sum_l\log\dfrac{\Pr(r_l|c=0)}{\Pr(r_l|c=1)}=\sum_lL(r_l|c) L(c∣r)=logPr(c=1∣r)Pr(c=0∣r)=log∏lPr(rl∣c=1)∏lPr(rl∣c=0)=l∑logPr(rl∣c=1)Pr(rl∣c=0)=l∑L(rl∣c)
结论:重复码的MAP解码器等于所有信道输出 r l r_l rl的LLR的总和。 -
SPC码的LLR
随机变量 a ⃗ \vec{a} a中第 l l l个元素为1的概率是 p 1 ( l ) p_1^{(l)} p1(l),则 a ⃗ \vec{a} a中有偶数个1的概率是 0.5 + 0.5 ∏ l ( 1 − 2 p 1 ( l ) ) 0.5+0.5\prod_l(1-2p_1^{(l)}) 0.5+0.5l∏(1−2p1(l))有奇数个1的概率是 0.5 − 0.5 ∏ l ( 1 − 2 p 1 ( l ) ) 0.5-0.5\prod_l(1-2p_1^{(l)}) 0.5−0.5l∏(1−2p1(l))
因为SPC码的形式是 c d = c 0 + c 1 + … + c d − 1 c_d=c_0+c_1+\ldots+c_{d-1} cd=c0+c1+…+cd−1,码字中一定有偶数个1,所以 c 0 = 0 c_0=0 c0=0等价于剩下的元素中有偶数个1: Pr ( c 0 = 0 ∣ r ⃗ ) = Pr ( c 1 , … ,