基于外信息符号变化的低密度校验码译码方法转让专利
申请号 : CN200910051181.X
文献号 : CN101552613B
文献日 : 2011-08-31
发明人 : 宫良 , 归琳 , 熊箭 , 刘勃 , 苗瑞琴
申请人 : 上海交通大学
摘要 :
本发明涉及一种通信系统中信道纠错译码技术领域的基于外信息符号变化的低密度奇偶校验码译码方法。对信号的硬判决向量判断是否满足校验方程，若满足，则译码成功，否则开始迭代译码。对初始软信息和上次迭代译码时与其相连的所有校验节点的外信息求和，得到第l次迭代译码时信息节点传递给校验节点的信息。对此信息和上次迭中传递给校验节点的信息进行符号比较，若相同，则将此信息传递给校验节点；若不同，则对此信息乘α，传递给校验节点。校验节点利用此信息得到新软信息。对其进行硬判断，若满足校验方程，则译码成功，否则进入下一次迭代。本发明缓解了现有基于软信息的迭代译码方法过程中错误符号数波动而导致译码器不收敛的问题。
权利要求 :
1.一种基于外信息符号变化的低密度校验码译码方法，其特征在于，包括如下步骤：①.对接收信号向量进行硬判决，得到码元硬判决向量，判断硬判决后的接收向量是否满足校验方程，如果满足则宣布译码成功；否则，计算软信息向量，开始迭代译码，置迭代译码次数l＝1，最大迭代次数为L；将软信息分量赋给对应的信息节点作为初始软信息；
②.第l次迭代译码时，对每个信息节点利用初始的软信息和上次迭代译码时来自与其相连的所有校验节点的外信息，计算信息，然后，对此信息和上次迭代中信息节点传递给校验节点的信息进行符号比较，如果相同，则将此信息传递给校验节点，如果不同，将此信息乘以α，α∈(0，1]，传递给校验节点；
③.每个校验节点利用来自与其相关联的信息节点的信息，计算新外信息并传回给对应的信息节点，每个信息节点利用初始的软信息和来自与其相连的校验节点的新外信息，计算该信息节点的新软信息，然后，对新软信息向量进行硬判决，判断硬判决向量是否满足校验方程，若满足，则译码成功；否则，继续下一步译码；
④.对迭代次数进行判断，如果迭代次数l＝L，则将硬判决向量作为最终的译码结果进行输出；否则，迭代译码次数l＝l+1，跳至步骤②，继续译码。
2.根据权利要求1所述的基于外信息符号变化的低密度校验码译码方法，其特征是，所述的校验方程，其表示如下：T
mod(uH，2)＝0
其中：u为硬判决向量，H为一个m行n列二元矩阵。
3.根据权利要求1所述的基于外信息符号变化的低密度校验码译码方法，其特征是，所述的软信息向量s＝[s1，…，sn]，其公式如下：2
其中：x＝[x1，x2，…，xn]为经过高斯白噪声信道后码元采样向量，σ 为高斯白噪声方差。
4.根据权利要求1所述的基于外信息符号变化的低密度校验码译码方法，其特征是，所述的外信息 ，其公式如下：i′表示与校验节点cj相连的信息节点序号集合中除去序号i后余下的信息节点，N(j)\i是与校验节点cj相邻的信息节点序号集合中除去序号i的剩余集合； 表示第l次迭代时，在来自与校验节点cj相连的信息节点集合中，除去信息节点vi之后，来自其它信息节点的信息； 其中e为自然对数的底。
5.根据权利要求1所述的基于外信息符号变化的低密度校验码译码方法，其特征是，所述的信息，其公式如下：其中： 表示第l次迭代译码时，信息节点vi传递给校验节点cj的信息；si表示初始软信息，j′表示与信息节点vi相连的校验节点序号集合中除去序号j的其余校验节点序号；N(i)\j是与信息节点vi相邻的校验节点集合中除去节点j的剩余集合； 表示上一次迭代译码时，来自校验节点cj′的外信息；N(i)为与信息节点vi相邻的校验节点序号集合。
6.根据权利要求1所述的基于外信息符号变化的低密度校验码译码方法，其特征是，所述的新软信息，其公式如下：其中 表示第l次迭代译码时的软信息，si表示初始软信息，j为序号；N(i)为与信息节点vi相邻的校验节点序号集合； 表示第l次迭代译码时，来自校验节点cj的外信息。
说明书 :
基于外信息符号变化的低密度校验码译码方法
技术领域
[0001] 本发明涉及的是一种通信技术领域的方法，具体是一种基于外信息符号变化的低密度校验码译码方法。
背景技术
[0002] 通信系统中，各种噪声的干扰会使得接收机所接收到的信息符号出现错误，而信道纠错编码可以用来对抗这种干扰。信道纠错编码是在发送的信息中加入冗余信息，从而接收机可以利用被传输的信息与冗余信息之间的特定关系，消除错误，恢复被传输的信息。
[0003] 迭代译码是将译码过程变为一个往复循环的过程，通过一次次的迭代循环，逐步消除传输信息中的错误；相比于传统的一次译码方法，迭代译码具有更强的纠错能力。低密度奇偶校验码(Low-Density Parity-Check Codes)是一种线性分组码，校验矩阵是低密度矩阵(即矩阵中1的个数远大于0的个数)，如果低密度校验矩阵中每行的1个数恒定，并且每列中1的个数也恒定，这样的码称为规则码；若行列中的1个数不恒定，称为非规则码。低密度奇偶校验码的译码采用基于置信传播的迭代译码方法。
[0004] 低密度奇偶校验码的迭代译码根据译码时所利用的信息和对信息的处理方式不同可以分为：比特翻转算法(Bit-Flipping Algorithm)，和积算法(Sum-Product Algorithm，SPA)，最小和算法(Min-Sum Algorithm，MSA)以及修正的最小和方法(Modified Min-Sum Algorithm，MMSA)等。其中，比特翻转算法是基于硬判决信息的迭代译码方法，后三种是基于软信息的译码方法。采用基于软信息的译码方法获得的译码性能要优于基于硬判决的译码方法(通常会优于后者2-3dB)。但是，基于软信息的译码方法在迭代译码过程中会出现译码不收敛的现象，即随着迭代次数的增加，错误符号数不呈现下降趋势，而是出现波动，最终无法收敛到0。
[0005] 经对现有技术文献的检索发现，S.Gounai等在《IEEE VehicularTechnology Conference论文 集》(pp.1467-1471，2006)上 发表 了“ModifiedBelief Propagation Decoding Algorithm for Low-Density Parity Check CodeBased on Oscillation”(“基于错误抖动的置信传播改进算法”)。该技术针对基于软信息译码不收敛的现象，在每次迭代译码时，比较本次信息节点传递给校验节点的信息的符号与上次的信息符号之间的变化，如果相同则直接传递给校验节点，如果不同，要将本次迭代信息与上次迭代的信息合并再传输，该技术的缺点是需要保存上次迭代时信息节点传递给校验节点的信息。
发明内容
[0006] 本发明的目的在于克服现有技术的不足，提供一种基于外信息符号变化的低密度奇偶校验码译码方法。本发明通过检测相继的两次迭代译码过程中，信息节点发送给校验节点的外信息符号编码，抑制符号改变的外信息数值，从而缓解现有的基于软信息的迭代译码方法过程中错误符号数波动而导致译码器不收敛的问题。
[0007] 本发明通过以下技术方案实现，包括如下步骤：
[0008] 1.对接收信号向量x＝[x1，x2，L，xn]进行硬判决，得到判决向量u＝[u1，u2，L，un]，判断硬判决后的接收向量是否满足校验方程，即检验硬判决向量的各个分量与校验矩阵的每一列各个分量相乘后再模2求和的结果是否为都零，其数学表达式如下：
[0009] mod(uHT，2)＝0
[0010] 如果上式满足；则宣布译码成功；否则，根据接收信号向量计算软信息向量s＝[s1，L，sn]，表达式为：
[0011]
[0012] 其中：H为一个m行n列二元矩阵，σ2为高斯白噪声方差；开始迭代译码。
[0013] 2.迭代译码初始化，置迭代译码次数l＝1，最大迭代次数为L；将软信息分量si赋给对应的信息节点vi作为初始软信息，i＝1，L，n；置信息节点到校验节点信息符号：
[0014] λi→j＝sgn(si)，j∈N(i)，i＝1，L，n
[0015] 置校验节点传递给信息节点的外信息：
[0016]
[0017] 其中，λi→j表示信息节vi点传递给校验节点cj的信息的符号，qj→i0为校验节点传递给信息节点的初始外信息，vi表示第i个信息节点，N(i)为与信息节点vi相邻的校验节点序号集合；N(j)为与校验节点cj相邻的信息节点序号集合。
[0018] 3.第l次迭代译码时，对每个信息节点vi利用初始的软信息si和上次迭代译码l-1 l时来自与其相连的所有校验节点cj的外信息qj→i ，j∈N(i)，计算信息ri→j，[0019]
[0020] 然后，将sgn(ri→jl)的值与λi→j的值进行比较，如果相同，则将ri→jl传递给校验节点cj，如果不同，置 令信息 传递给校验节点cj。l
[0021] 其中：ri→j 表示第l次迭代译码时，信息节点vi传递给校验节点cj的信息；j′表示与信息节点vi相连的校验节点序号集合中除去序号j的其余校验节点序号；N(i)\j是l-1与信息节点vi相邻的校验节点集合中除去节点j的剩余集合；qj′→i 表示上一次迭代译码时，来自校验节点cj′的外信息；其中α∈(0，1]。
l
l
[0022] 4.每个校验节点cj利用来自与其相关联的信息节点vi的信息ri→j，i∈N(j)，l计算新外信息qj→i 并传回给对应的信息节点vi；其中：
l
l
[0023] N(j)为与校验节点cj相邻的信息节点序号集合；ri→j 表示第l次迭代译码时信l息节点vi传递给校验节点cj的信息；qj→i 表示第l次迭代译码时，校验节点cj传递给信l
息节点vi的外信息，qj→i 的计算可以采用现有的基于软信息的译码方法中的计算方法。
息节点vi的外信息，qj→i 的计算可以采用现有的基于软信息的译码方法中的计算方法。
[0024] 5.每个信息节点vi利用初始的软信息si和来自与其相连的校验节点的外信息l lqj→i，计算该信息节点的新软信息si，
[0025]
[0026] 然后，对各个sil进行硬判决，得到判决向量ul，若硬判决后的信息向量满足校验方程，则译码成功；否则，继续下一步译码；其中：sil表示第l次迭代译码时的信息节点的软信息；ul表示第1次迭代译码时的硬判决向量。
[0027] 6.对迭代次数进行判断，如果迭代次数l＝L，则将ul作为最终的译码结果进行输出；否则，迭代译码次数l＝l+1，跳至3，继续译码。
[0028] 本发明不需要存储信息节点发送给校验节点的信息值，只需要存储每次迭代信息节点发送给校验节点的符号，从而大大降低了存储的复杂度；在信息节点发送给校验节点的外信息符号发生变化时，对新信息乘以预先设定的抑制因子α，从而抑止相邻的两次迭代译码过程中的符号突变造成的错误符号蔓延，进而缓解迭代译码过程中错误符号数的波动现象，从而提高译码成功的概率；与现有的基于软判决的低密度校验码译码方法相比，译码器输出的误比特率下降50％。同时通过适当选择抑制因子，对信息的乘积运算在实际硬件实现时，可以直接采取寄存器的移位操作，从而进一步使得计算复杂度降低。
附图说明
[0029] 图1为校验矩阵H的二分图表示。
[0030] 图2为信息节点处外信息的更新与传递。
[0031] 图3为校验节点出外信息的更新与传递。
[0032] 图4为一个给定的接收信号向量样本采用本发明方法译码和SPA方法译码，每次迭代译码后的误比特率变化统计曲线。
[0033] 图5为1300个接收信号样本采用本发明方法译码和SPA方法译码方式下不同迭代循环中的平均误比特数统计。
[0034] 图6本发明译码方式和SPA译码方式下的误比特率和误帧率性能曲线。
具体实施方式
[0035] 下面结合附图对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
[0036] 本实施例涉及一种基于外信息符号变化的低密度校验码译码方法，包括如下步骤：
[0037] 1.对接收信号向量x＝[x1，x2，L，xn]进行硬判决，硬判决规则为：
[0038]
[0039] 得到判决向量u＝[u1，u2，L un]，判断硬判决后的接收向量是否满足校验方程，即计算校验向量p＝mod(uHT，2)，如果结果为全0向量，即满足则宣布译码成功；否则，计算作为译码的初始信息赋给信息节点，并置迭代译码次数l＝1，最大迭代次数为L；其中：
[0040] x＝[x1，x2，L，xn]为经过高斯白噪声信道后码元采样向量p＝mod(uHT，2)＝[p1，T Tp2，...，pm]；mod(uH，2)表示对uH 取模2；H为m×n阶低密度校验矩阵；s＝[s1，s2，L，sn]
2
为初始软信息向量；σ 为高斯白噪声的方差；
2
为初始软信息向量；σ 为高斯白噪声的方差；
[0041] 如图1所示vi(i＝1，2，3，...，n)节点(圆圈表示)代表信息节点，分别对应码字的各个符号，即校验方程的列；cj(j＝1，2，...，m)节点(方框表示)，代表校验节点，分别对应每个校验方程的行。
[0042] 2.迭代译码初始化，置迭代译码次数l＝1，最大迭代次数为L；将软信息分量si赋给对应的信息节点vi作为初始软信息，i＝1，L，n；置信息节点到校验节点信息符号，[0043] λi→j＝sgn(si)，j∈N(i)，i＝1，L，n；
[0044] 置校验节点传递给信息节点的外信息
[0045] i∈N(j)，j＝1，L，m；
[0046] 其中，λi→j表示信息节vi点传递给校验节点cj的信息的符号，sgn(x)为取符号操作，即sgn(x)＝1，若x≥0；sgn(x)＝-1，若x＜0；vi表示第i个信息节点，cj表示第j个校验节点，；N(i)为与信息节点vi相邻的校验节点序号集合；N(j)为与校验节点cj相邻的信息节点序号集合。
[0047] 3.在第l次迭代译码时，对每个信息节点vi利用初始的软信息si和来自与其相l-1 l连的校验节点的外信息qj→i ，j∈N(i)，计算信息ri→j，
[0048]
[0049] 然后，将sgn(ri→jl)与λi→j进行比较，如果相同，则将新的外信息传递给校验节l点cj，；如果不同， 则并对信息ri→j 乘以α，传递给校验节点cj；
[0050] 其中：ri→jl表示第l次迭代译码时，信息节点vi传递给校验节点cj的信息；j′表示与信息节点vi相连的校验节点序号集合中除去序号j的其余校验节点序号；N(i)\j是l-1与信息节点vi相邻的校验节点集合中除去节点j的剩余集合；qj′→i 表示上一次迭代译码时，来自校验节点cj′的外信息；其中α∈(0，1]。
[0051] 如图2所示：图(a)为在第l次迭代译码循环中，信息节点vi利用上次迭代时，来自与其相关联的|N(i)|个校验节点 的外信息对这|N(i)|个校验节点 jk∈N(i)分别计算本次迭代译码的
信息ri→jl，其中， 表示与信息节点vi相连的第k个校验节点，其序号为jk。
信息ri→jl，其中， 表示与信息节点vi相连的第k个校验节点，其序号为jk。
[0052] 如图2(b)所示：首先，信息节点vi计算外信息 然后将sgn(ri→jl)与λi→j进行符号比较。如果同号，则ri→jl保持不变传递给节点 否则，令计算 将ri→jl传递给校验节点 其中：α为预先设定的抑制因
子α，取值区间为(0，1]，本实施例中α＝0.125。
子α，取值区间为(0，1]，本实施例中α＝0.125。
[0053] 4.每个校验节点cj利用来自与其相连的信息节点的信息ri→jl，i∈N(j)，计算新的外信息qj→il并传给各个对应的信息节点vi，其中：
[0054] N(j)为与校验节点cj相邻的信息节点序号集合；qj→il表示第l次迭代译码时，校验节点cj传递给信息节点vi的外信息，计算公式如下：
[0055]
[0056] i′表示与校验节点cj相连的信息节点序号集合中除去序号i后余下的信息节l点，N(j)\i是与校验节点cj相邻的信息节点序号集合中除去序号i的剩余集合；ri′→j 表示第l次迭代时，在来自与校验节点cj相连的信息节点集合中，除去信息节点vi之后，来自其它信息节点的信息； 其中e为自然对数的底；|x|对于数字表示取绝对值；对于集合表示取集合的元素个数。
[0057] 如图3所示：校验节点cj接收来自与其相关联的|N(j)|个信息节点的信息：利用这些信息分别为各个信息节点 ik∈N(j)计算新的外信息
并将其传回相对应的信息节点 其中 表示与第j个校验
节点相关联的第k个信息节点，其对应的信息节点序号为ik。
并将其传回相对应的信息节点 其中 表示与第j个校验
节点相关联的第k个信息节点，其对应的信息节点序号为ik。
[0058] 5.每个信息节点vi利用初始的软信息si和来自与其关联的所有校验节点的外信l息qj→i 计算关于该信息节点的新软信息
[0059]
[0060] 对sl进行硬判决得到判决向量ul；若硬判决后的信息向量满足校验方程，则译码l成功；否则，继续下一步译码；其中，si 表示第l次迭代译码时的软信息，s＝[s1，L sn]表l
示第l次迭代译码时的软信息向量，u 表示第l次迭代译码时的硬判决向量；
示第l次迭代译码时的软信息向量，u 表示第l次迭代译码时的硬判决向量；
[0061] 6.对迭代次数进行判断，如果迭代次数l＝L，则将ul作为最终的译码结果进行输出；否则，迭代译码次数l＝l+1，跳至3继续译码。
[0062] 图4为BPSK调试方式下，选取IEEE802.16e标准中长度为576的低密度奇偶校验码，经过高斯白噪声信道后的一个输出样本，分别使用SPA译码方法(使用实心圆点表示)及基于外信息符号变化译码方法(使用实心方块表示)仿真时，50次译码迭代过程中的错误比特数变化曲线。图中的横坐标为迭代次数，纵坐标为误比特数。可以看到，SPA方法在各个迭代译码循环中的错误比特数呈现波动，并且在最大迭代次数50次达到后，无法收敛到0；而采用基于外信息符号变化的译码方法时，由于变号的外信息幅度被抑制，译码器在第22次译码循环时错误比特数收敛到0，从而成功译码。
[0063] 图5为BPSK调试方式下，选取IEEE802.16e标准中长度为576的低密度校验码，经过高斯白噪声信道后的1300个输出样本，分别使用SPA译码方法(使用实心圆点表示)及基于外信息符号变化译码方法(使用实心方块表示)仿真时，50次译码迭代过程中的平均错误比特数变化曲线。图中的横坐标为迭代次数，纵坐标为平均误比特数。可以看到，SPA方法在各个迭代译码循环中的平均错误比特数呈现波动，并在最大迭代次数50达到后，仍未呈现收敛趋势；而采用基于外信息符号变化的译码方法后，平均错误比特数随着迭代次数的增加波动趋势减小，并呈现稳定的下降趋势。
[0064] 图6为BPSK调试方式下，高斯白噪声信道中，对IEEE802.16e标准中长度为576的低密度校验码，最大迭代译码次数50次，分别采用SPA译码方法(使用圆形标志)及基于外信息符号变化译码方法(使用方块符号标志)仿真时，误比特率(使用实心符号标志)的变化曲线和误帧率(使用空心符号表示)变化曲线。图中的横坐标表示接收信号的信噪比，纵坐标表示译码后的错误比率，即误比特率和误帧率。可以看到，采用基于外信息符号变化的译码方法无论是误帧率还是误比特率曲线都低于SPA方法获得的曲线；并且相同信噪比下前者的误比特率和误帧率均为后者的50％左右。因此，基于外信息符号变化的译码方法，降低了译码的误比特和误帧率，即提高了译码成功率。