专利名称：一种优化ldpc码校验矩阵列重分布的方法
技术领域：
本发明涉及分组信道编码技术领域，尤其涉及一种优化LDPC(Low Density Parity Check code，低密度奇偶校验码)校验矩阵列重分布的方法。
背景技术：
现代编码技术的终极目标是以逼近香农极限Shannon limit的有效功耗实现可靠 通信，LDPC码的诞生和发展使人们更加接近这一目标。LDPC码最早在1962年由Gallager 提出，当时没有受到编码界的重视，直到1996年英国的Mackay教授等人通过仿真证明了 LDPC码的优秀性能，才轰动了编码界，成为自信息论提出以来最重大的研究进展之一。研究 表明，基于非规则双向图的LDPC码的性能优于Turbo码，具有更低的线性译码复杂度，没有 错误平层，因此受到广泛关注。LDPC码是一类线性分组码，由它的校验矩阵来定义，设码长为N，信息位为K，校验 位M = N-K，码率为R = K/N，则校验矩阵是一个MXN的矩阵。校验矩阵中很少一部分元素 非零，其他大部分元素都是零，所以LDPC码的校验矩阵是稀疏矩阵。矩阵中每一列非零元 素的数目称为该列的重。LDPC码的校验矩阵可以用一个与其对应的双向图来表示，例如，一个码长为10， 码率为0. 5的LDPC码的校验矩阵如图1所示，该校验矩阵对应的双向图如图2所示。双 向图的下边有N个节点，每个节点表示码字的信息位，这N个节点称为信息节点{\，j = 1，. . .，N}，分别对应于校验矩阵的N列；双向图的上边有M个节点，每个节点表示码字的一 个校验集，这M个节点称为校验节点{Zi，i = 1，...，M}，分别对应于校验矩阵的各行，代表 着校验方程。与校验矩阵中“1”相对应的左右两节点之间存在连接边，将这条边两端的节 点称为相邻节点，每个节点相连的边的数目称为该节点的度数。若LDPC码对应的双向图为规则双向图，则此LDPC码称为规则LDPC码，规则双向 图是指所有校验节点度数相等，信息节点度数也相等，规则LDPC码的校验矩阵中每行或列 中的“ 1 ”的个数是相等的；若LDPC码对应的双向图为非规则双向图，则此LDPC码称为非规 则LDPC码，非规则LDPC码校验矩阵中每行或列中的“1”的个数不相等。模拟实验说明适 当构造的不规则码的性能优于规则码的性能。这一点可以从构成LDPC码的双向图中直观 的得到解释信息节点的度数越大，它从相关联的校验节点得到的信息越多，便越能准确地 判断它的正确值，而对于校验节点，情况则相反，校验节点的度数越小，它能反馈给其邻接 信息节点的信息便越有价值。显然，不规则码比规则码能够更好的平衡这两种需求。PEG (Progressive Edge-Growth，渐进边增长)方法是一种基于随机搜索方式构造 LDPC码的重要方法。它以校验矩阵对应的双向图中环的长度最大化为目标，PEG方法具体 是从左到右逐列以当前的节点度数为循环次数，依次将连接到信息节点的边设定为与校验 节点之间的边，通过迭代搜索逐列构造LDPC码的校验矩阵，是一种参数配置非常灵活的方 法。研究表明只要列重为某一数值的列所占总列数的百分比即列重分布确定了，不管其具体分布方式如何，也不管随机搜索中随机数的生成方式如何，PEG方法构造出来的 LDPC校验矩阵，虽然具体形式各有不同，但是纠错性能都是非常接近的。因此，一个LDPC码 的校验矩阵的列重分布在很大程度上决定了 LDPC码的性能。
发明内容
本发明要解决的技术问题是，提供一种优化LDPC校验矩阵列重分布的方法，进一 步提高LDPC校验矩阵的性能，降低通信时的有效功耗。本发明采用的技术方案是，所述优化LDPC校验矩阵列重分布的方法，包括如下步 骤步骤一、根据不同的码率选择校验矩阵列重的取值集合，基于所述列重的取值集 合随机生成合法的列重分布，估计每一个列重分布下LDPC码的性能；步骤二、将性能最差的LDPC码校验矩阵对应的列重分布记为X，再将性能倒数第 二差的LDPC码校验矩阵对应的列重分布记为r，计算除了 χ之外剩下的列重分布的均值，并 将其记为1 ；步骤三、构造新的列重分布ζ = (1+α)γ-αχ,0 < α < 1 ；步骤四、将新的列重分布ζ合法化，然后估计其对应LDPC码的性能，若ζ对应LDPC 码的性能好于r对应LDPC码的性能，则将χ替换为ζ，跳转步骤五，否则令ζ = (z+y)/2，重 复执行步骤四；步骤五、判断所有列重分布之间的差距是否小于某一个阈值，若是，则算法结束， 否则重复步骤二。设步骤一中所述校验矩阵列重的取值集合中包含的元素个数为L，基于所述列重 的取值集合随机生成合法的列重分布的具体过程如下步骤1、任意设定一个合法的第一列重分布向量，所述第一列重分布向量中分量的 个数与所述列重的取值集合中元素的个数相等；步骤2、采用随机算法生成第二列重分布向量到第2L-2列重分布向量中的前L-I 个分量；步骤3、顺序检查所有列重分布向量是否合法，对不合法的列重分布向量进行调 整；步骤4、由分量和为1的约束计算出每个列重分布向量的最后一个分量。步骤一和步骤四中所述估计的方法是密度进化方法或者蒙特卡洛仿真。步骤四中所述将新的列重分布ζ合法化包括将所述新的列重分布ζ向量中的各 分量的数值调整到0和1之间，并保证各分量之和等于1。步骤五中所述列重分布之间的差距是列重分布向量模之间的差值，所述阈值根据 实际应用需求确定。采用上述技术方案，本发明至少具有下列优点本发明所述优化LDPC校验矩阵列重分布的方法，根据不同的码率选择校验矩阵 列重的取值集合，基于所述列重的取值集合随机生成合法的列重分布，可以在随机编码的 前提下，优化校验矩阵的列重分布。同时仿真实验表明，在相同的码率下，与用其他方式构 造的相同码长的LDPC码误码率曲线相比，采用本发明技术方案构造的校验矩阵的LDPC码的性能有明显的优势，能在进一步降低功耗的情况下实现可靠的通信。
图1为码长为10，码率为0. 5的LDPC码的校验矩阵结构示意图；图2为码长为10，码率0. 5的LDPC码的双向图；图3为现有技术中中短码长的LDPC码在BPSK调制AWGN信道下的误码率仿真曲 线.
一入 ，图4为采用本发明所述方法构造的校验矩阵后终端中短码长的LDPC码在BPSK调 制AWGN信道下的误码率仿真曲线；图5为本发明所述优化LDPC校验矩阵列重分布的方法流程图。
具体实施例方式为更进一步阐述本发明为达成预定目的所采取的技术手段及功效，以下结合附图 及较佳实施例，对本发明提出的所述优化LDPC校验矩阵列重分布的方法，详细说明如后。本发明的具体实施例，一种优化LDPC校验矩阵列重分布的方法，如图5所示，包括 如下步骤步骤一、根据不同的码率选择校验矩阵列重的取值集合，设集合中元素的个数即 集合的大小为1。由于在不同的码率下，相同列重数值所在的列占整个校验矩阵总列数的比 例是不同的，本领域技术人员可以根据实际中不同码率下关心的列重数值选择一个校验矩 阵列重的取值集合，也可以根据不同的码率下列重数值所在的列占整个校验矩阵总列数的 比例的经验值等因素选择一个校验矩阵列重的取值集合。步骤二、基于所述列重的取值集合随机生成21-2个合法的列重分布。一个合法的 列重分布指的是一个长度为列重数值取值集合大小1的向量，可以称为列重分布向量，其 各个分量为列重取值集合中各列重所在的列的个数占总列数的比例，因此各个分量取值均 在0，1之间，且所有分量的和为1。令这些向量为Λ j (j = 1，2，. . .，21-2)，Λ丨是Λj的第i 个分量，i = 1，2，. . .，1。具体地，所述随机生成21-2个合法的列重分布过程如下步骤a、首先任意设定一个满足合法分布要求的第一列重分布向量A1 ；步骤b、构造余下的21-3个列重分布向量，具体地，在构造第二列重分布向量A2 时，将w，l]范围内的随机数依次赋值给前1-1个分量Λ (ζ· = 1,…，/-1)，按照同样的方法构 造第三列重分布向量A3........第21-2列重分布向量Λ21_2 ；步骤C、按照j的取值从小到大的顺序检查$Λ； <1的约束条件是否得到满足，若 在j取值的范围内,》Κ< 1的约束条件都能得到满足，则所述随机生成21-2个合法的列重 分布过程结束，若在j的前η个取值下，<1分布约束条件得到满足，l^n^ 21-3，而在 j的第n+1个取值时不满足，则进行下面的替换操作Λη+1 = (An+1+t)/2，其中 =，
5重复本步骤c，直到对所有的j的取值，ΣΥ<:l分布约束条件均能得到满足；步骤d、计算21-2个的列重分布的最后一个分量，即进行如下计算 步骤三、估计每一个列重分布下LDPC码的性能，所述估计的方法可以是密度进化 方法或者蒙特卡洛仿真方法。密度进化方法是在实际译码算法上抽象地、理想化地迭代过 程，它的应用可以看作是一种数值仿真实验。密度进化可以在给定LDPC矩阵列重分布的 前提下，通过迭代演算，得到码长足够长且LDPC矩阵无环的理想情况下该列重分布对应的 LDPC码的性能，可以通过误码率曲线表现出来。在本技术领域中，通常采用蒙特卡洛仿真方法对中短码长LDPC码的性能进行估 计，中短码长为2000 5000比特，而采用密度进化方法对长码LDPC码的性能进行估计，长 码长为5000比特以上，因此根据不同的LDPC码选择适当的性能估计方法，可以使得性能评 估的结果针对性更强。因为上述两种估计的方法是本领域技术人员在估计LDPC码性能时 常用的技术手段，所以此处不详述。步骤四、将性能最差的LDPC码校验矩阵对应的列重分布记为Aw，再将性能倒数 第二差的LDPC码校验矩阵对应的列重分布记为if，计算除Aw之外剩下的列重分布的均值 步骤五、构造新的列重分布Ar= (1+α) ΛΕ-α AwjO < α <1，优选的，α = 0. 5 ；步骤六、将新的列重分布Λ1合法化，然后按照步骤三中的估计方法分别估计” 和if对应LDPC码的性能，若Λ 对应LDPC码的性能好于Zf对应LDPC码的性能，则将Λ1替 换为Λ\跳转步骤七，否则令” =(Λ4Λκ)/2，重复执行步骤六；所述将新的列重分布Λr合法化的过程如下步骤A、对Λ1中的各分量的数值进行调整，使0<《<l,i = 1,2,... ,I-I0具体 地，若Λ〗<0 ,则进行赋值操作Λ〗=Λ〗若Λ〗>1,则进行赋值操作Λ〗=乂-5’ δ通
常设定为10_5。步骤B、检查约束条件
是否满足，如果满足，则所述将新的列重分布Λ1合
法化的过程结束；如果不满足，则进行赋值操作
7^重复步骤A。步骤七、判断所有列重分布之间的差距是否小于预设的阈值，若是，则算法结束， 否则重复步骤四。所述列重分布之间的差距可以用列重分布向量模的差值来衡量，所述阈 值根据实际应用需求确定。采用本发明所述优化LDPC校验矩阵列重分布的方法确定的8种码率下中短码长 LDPC码校验矩阵的理想列重分布如下，其中Yi表示列重为i的列占总列数的比例 经过实验仿真，在8种码率下，现有技术中LDPC码在BPSK调制AWGN信道下的误 码率如图3所示，采用本发明所述方法优化LDPC码校验矩阵后，LDPC码在BPSK调制AWGN信道下的误码率如图4所示，解码均采用BP (Belief Propagation，置信概率传播)算法，图 中横坐标为传递1比特信息所用的能量与噪声的比值，纵坐标为误码率，显然，图4中误码 率曲线即有效功耗更逼近香浓极限Shannon limit，提高了 LDPC码的性能，在进一步降低 功耗的情况下实现可靠的通信。 通过具体实施方式
的说明，应当可对本发明为达成预定目的所采取的技术手段及 功效得以更加深入且具体的了解，然而所附图示仅是提供参考与说明之用，并非用来对本 发明加以限制。
权利要求
一种优化LDPC校验矩阵列重分布的方法，其特征在于包括如下步骤步骤一、根据不同的码率选择校验矩阵列重的取值集合，基于所述列重的取值集合随机生成合法的列重分布，估计每一个列重分布下LDPC码的性能；步骤二、将性能最差的LDPC码校验矩阵对应的列重分布记为x，再将性能倒数第二差的LDPC码校验矩阵对应的列重分布记为r，计算除了x之外剩下的列重分布的均值，并将其记为y；步骤三、构造新的列重分布z＝(1+α)y αx，0＜α＜1；步骤四、将新的列重分布z合法化，然后估计其对应LDPC码的性能，若z对应LDPC码的性能好于r对应LDPC码的性能，则将x替换为z，跳转步骤五，否则令z＝(z+y)/2，重复执行步骤四；步骤五、判断所有列重分布之间的差距是否小于某一个阈值，若是，则算法结束，否则重复步骤二。
2.根据权利要求1所述优化LDPC校验矩阵列重分布的方法，其特征在于设步骤一中所 述校验矩阵列重的取值集合中包含的元素个数为L，基于所述列重的取值集合随机生成合 法的列重分布的具体过程如下步骤1、任意设定一个合法的第一列重分布向量，所述第一列重分布向量中分量的个数 与所述列重的取值集合中元素的个数相等；步骤2、采用随机算法生成第二列重分布向量到第2L-2列重分布向量中的前L-I个分量；步骤3、顺序检查所有列重分布向量是否合法，对不合法的列重分布向量进行调整； 步骤4、由分量和为1的约束计算出每个列重分布向量的最后一个分量。
3.根据权利要求2所述优化LDPC校验矩阵列重分布的方法，其特征在于步骤一和步骤 四中所述估计的方法是密度进化方法或者蒙特卡洛仿真。
4.根据权利要求1或2或3所述优化LDPC校验矩阵列重分布的方法，其特征在于步骤 四中所述将新的列重分布ζ合法化包括将所述新的列重分布ζ向量中的各分量的数值调整到0和1之间，并保证各分量之和 等于1。
5.根据权利要求4所述优化LDPC校验矩阵列重分布的方法，其特征在于步骤五中所述 列重分布之间的差距是列重分布向量模之间的差值，所述阈值根据实际应用需求确定。
全文摘要
本发明公开了一种优化LDPC校验矩阵列重分布的方法，根据不同的码率选择校验矩阵列重的取值集合，基于所述列重的取值集合随机生成合法的列重分布，可以在随机编码的前提下，优化校验矩阵的列重分布。在相同的码率下，与用其他方式构造的相同码长的LDPC码误码率曲线相比，采用本发明技术方案构造的校验矩阵的LDPC码的性能有明显的优势，能在进一步降低功耗的情况下实现可靠的通信。
文档编号H03M13/11GK101908893SQ20091008518
公开日2010年12月8日 申请日期2009年6月3日 优先权日2009年6月3日
发明者吴智勇, 徐淑正, 金庆文, 高鹏 申请人:国家广播电影电视总局广播科学研究院