一种自适应QC-LDPC码的确定性构造方法与流程

一种自适应qc-ldpc码的确定性构造方法
技术领域
1.本发明涉及一种自适应qc-ldpc码的确定性构造方法，属于卫星信道编码技术领域。


背景技术：

2.随着卫星通信业务的持续扩大以及对卫星组网需求的不断提升，卫星转发器也朝着向具备再生式处理功能的方向发展。相比于透明转发器，处理转发器虽然在频谱利用率、链路性能、通信容量、组网灵活性等方面具有优越性，但由于增大了星上处理设备的复杂度，也增加了对卫星载荷和电源功率的要求，因此限制了编码方式的灵活选择。
3.码的结构决定了ldpc码的性能。借助于计算机搜索，人们已经提出一些大围长ldpc码的构造方法，虽然纠错性能较好，但其校验矩阵中非零元素的分布没有规律，导致其节点信息无法用映射关系或者递推公式表征，因此在编解码时需要存储校验矩阵和生成矩阵的所有列向量，从而在vlsi上的实现变得困难。以dvb-s2标准为例，它对长码和短码分别规定了11种和10种码率，由此可见，其qc-ldpc码的信息码向量长度、校验码向量长度都已具体规定而不能是自适应选择。当码率需要自适应调整时，会根据已经存储的21个校验矩阵列向量地址表进行编码选择，因此会带来计算复杂度的增加并占用大量的存储空间，这对于处理转发器在实现自适应抗衰落的应用中显然是不利的。
4.因此，研究信道编码时码参数的动态调整策略，以适应不同衰减强度下的卫星信道环境，同时尽可能地减小多样化参数带来的大规模集成电路存储问题，是实现自适应编码优化设计的关键。


技术实现要素：

5.本发明解决的技术问题为：针对qc-ldpc码中准循环基矩阵的移位系数确定问题，提出基于代数构造的确定方法。该方法构造的校验矩阵围长为8，码参数可自适应选取，移位系数由简单的数学表达式确定，编码复杂度与码长呈线性关系，节省了编解码存储空间。
6.本发明采用的技术方案为：
7.一种自适应qc-ldpc码的确定性构造方法，步骤如下：
8.(1)将稀疏校验矩阵分为校验部分和信息部分，其中校验部分采用双对角线结构；
9.(2)利用信息部分与单位矩阵生成等效矩阵，提取等效矩阵的循环移位系数形成基矩阵；
10.(3)按照递推公式对基矩阵进行赋值；
11.(4)求解基矩阵中各元素的最大值，作为循环移位矩阵尺寸的下界；
12.(5)根据卫星信道状态对编码参数进行自适应调整；
13.(6)利用构造好的校验矩阵直接对信息进行编码。
14.所述步骤(1)中的稀疏校验矩阵h分为校验部分ha和信息部分hb：
[0015][0016]
其中，m、n、p均为正整数，p
m,n
为循环移位系数，i(p
m,n
)表示一个p
×
p 的循环移位矩阵；ha为双对角线结构，大小为mp
×
mp，对角线位置由单位阵i 和零矩阵0构成；hb大小为mp
×
np。
[0017]
所述步骤(2)中的等效矩阵由信息部分hb与单位矩阵生成：
[0018][0019]
等效矩阵q具备与校验矩阵h相一致的围长特性，由此将构造过程等效为对q的构造。
[0020]
提取等效矩阵q的循环移位系数形成基矩阵p：
[0021][0022]
其中，p
1,0
＝p
2,0
＝...＝p
m,0
＝0。
[0023]
所述步骤(3)中对p进行赋值的递推公式为：
[0024][0025]
其中，d为非负整数，k表示行指数因子，，d
i,j
＝p
i,j-p
i,j-1
，1≤i≤m， 1≤j≤n。由此构造的等效矩阵q的围长至少为8。
[0026]
所述步骤(4)中基矩阵p的最大元素值为p
m,n
，而循环移位矩阵的维数p 应大于循环移位系数，故p的取值下界为：
[0027]
inf p＝p
m,n
1
[0028]
所述步骤(5)中的编码参数为m、n、d、p，可根据卫星信道状态进行自适应调整，从
而按照码长和码率要求确定性的得到校验矩阵h。
[0029]
所述步骤(6)中构造好的校验矩阵h是非奇异的，因此可回避生成矩阵 g而直接由校验矩阵进行编码。
[0030]
本发明与现有技术相比具有如下有益效果：
[0031]
(1)本发明设计的准循环基矩阵由数学公式确定地得到，无需依靠计算机搜索或存储大量矩阵形式就能得到任意参数的校验矩阵；
[0032]
(2)本发明可根据卫星信道状态或系统性能需求，仅调整初始参量便可灵活选取不同的编码方案，同时保持一致的大围长特性；
[0033]
(3)本发明可回避复杂的生成矩阵而直接由校验矩阵参与编码，编码器的实际存储量需求非常小；
[0034]
(4)本发明构造的循环移位矩阵维数具有连续取值的下界，与经典码型相比构造方式更加灵活。
附图说明
[0035]
图1为本发明方法的处理流程框图；
[0036]
图2为校验矩阵中6环存在的六种形式；
[0037]
图3为本发明方法设计的qc-ldpc码与peg码的性能比较。
具体实施方式
[0038]
下面结合附图和实施例对本发明作进一步介绍。
[0039]
如图1所示，本发明针对具备深度覆盖及高速灵活的卫星通信业务需求，旨在通过优化参数调整策略以进一步完善自适应编码关键技术。作为一种确定性方法，无需计算机搜索或存储大量不同形式的校验矩阵，就可由数学表达式直接计算得到循环移位系数；编码时只需存储少量初始值和生成式，节省了存储空间；任意参数下都能保证大围长特性，有效降低硬件资源的开销。
[0040]
具体步骤如下：
[0041]
(1)将稀疏校验矩阵h分为校验部分ha和信息部分hb：
[0042][0043]
其中，m、n、p均为正整数，p
m,n
为循环移位系数，i(p
m,n
)表示一个p
×
p 的循环移位矩阵；ha为双对角线结构，大小为mp
×
mp，对角线位置由单位阵i 和零矩阵0构成；hb大小为mp
×
np。
[0044]
(2)为简化构造过程，利用信息部分hb与单位矩阵设计等效矩阵q：
[0045][0046]
易知q与h具备相一致的围长特性，因此提取q的循环移位系数形成基矩阵p：
[0047][0048]
其中，p
1,0
＝p
2,0
＝...＝p
m,0
＝0。由此将对校验矩阵h的构造转化为对基矩阵p的构造。
[0049]
(3)对基矩阵p进行代数化构造的递推公式为：
[0050][0051]
其中，d为非负整数，k表示行指数因子，，d
i,j
＝p
i,j-p
i,j-1
，1≤i≤m， 1≤j≤n。
[0052]
下面对p的围长特性进行分析：
[0053]
不失一般性，令1≤r＜s＜t≤m，0≤i＜j＜l≤n。因此行指数r、s、t存在以下4种关系：
[0054]
(i)4k 1≤r＜s＜t≤4k 4,k≥0；(ii)4k 1≤r＜s≤4k 4＜t,k≥0；
[0055]
(iii)r＜4k 1≤s＜t≤4k 4,k≥1；(iv)r＜4k 1≤s≤4k 4＜t,k≥1。
[0056]
由递推公式易知p中任一行元素均构成单调递增数列，因此显然不存在4 环及图2(a)～(d)所示6环。故首先考虑关系(iv)和图2(e)的情形，此时有：
[0057][0058]
由于式中各项不满足同时取得最值的情形，因此经讨论可得：
[0059]
(p
t,l-p
t,j
) (p
r,j-p
r,i
)-(p
s,l-p
s,i
)
[0060]
≥n
·
(d
4k 4,1-d) d 1 n
·
d-n
·d4k 4,1
[0061]
＝d 1＞0
[0062]
即p中不存在图2(e)所示的6环。
[0063]
同理可得，当行指数满足关系(i)、(ii)、(iii)，以及图2(f)情形时，由递推公式定
义的基矩阵p的围长至少为8。
[0064]
(4)可以发现，基矩阵p的最大元素值为p
m,n
，而循环移位矩阵的维数p 应大于循环移位系数，故p具有连续取值的稳定下界：
[0065]
inf p＝p
m,n
1
[0066]
(5)信道编码参数包括码长n和码率r：
[0067][0068]
当码长一定时，维数p越大，则基矩阵尺寸越小；同时p的取值与d密切相关。因此当卫星信道状态发生变化时，m、n、d、p等变量可根据码长和码率要求进行自适应调整，从而确定性的得到校验矩阵h。
[0069]
(6)由于校验部分ha满秩，因此h是非奇异的，即可回避生成矩阵g而直接由校验矩阵进行编码。采用lu分解法，令信息比特向量为s＝[s
1 s2ꢀ…ꢀ
sn]，校验比特向量为p＝[p
1 p2ꢀ…ꢀ
pm]，编码器输出行向量为c，则有:
[0070]
c＝[p s]
[0071]
根据校验等式h
·ct
＝0可得：
[0072]
ha·
p
t
hb·st
＝0
[0073]
由于运算在gf(2)中进行，所以有：
[0074][0075]
将上式展开并由第1式可得：
[0076][0077]
将上式回代到方程组的第2式可得：
[0078][0079]
依次类推可得到各个pi，编码完毕。
[0080]
编码复杂度主要关注编码过程的运算量、运算复杂度和编码所需存储的参数。运算量即乘法和加法次数，运算复杂度即运算量与码长的变化关系。 qc-ldpc码的各个子矩阵都是稀疏矩阵，因此按照稀疏矩阵的运算方式可大大减小运算量。现计算pi(1≤i≤m)的运算量为：
[0081]
乘法次数：n-m
·
p
[0082]
加法次数：n-(m-i 2)
·
p
[0083]
可以明显地看到，计算各个校验分向量的运算复杂度为o(n)，即运算复杂度与码
长呈线性关系。同时由于采用代数构造法，校验矩阵中的元素由简单的数学运算得到，只需要存储一组初始参数和几个数学表达式即可，编码器的实际存储量需求非常小。
[0084]
实施例：
[0085]
利用本发明方法构造2种(4,8)准循环基矩阵，设码率r＝1/2，其它参数分别为：码长n＝512，维数p＝64，d＝p
1,1
＝8，p
2,1
＝9，p
3,1
＝14，p
4,1
＝ 15；码长n＝1200，维数p＝150，d＝p
1,1
＝30，p
2,1
＝31，p
3,1
＝36，p
4,1
＝37。在awgn信道下进行仿真，译码采用bp算法，最大迭代次数为30，调制方式为bpsk，选取同码长码率的peg码进行性能比较，仿真结果如图 3所示。
[0086]
按照本实施例设置条件，n＝512时二者的译码性能差异不大；当n＝ 1200、ber为10-5
时，与peg码相比所构造码字信噪比增益约为0.2db。此外，peg方法需要对度数的分布进行优化处理，其算法复杂度也高于所发明方法。
[0087]
以上所述，仅为本发明的较佳实施例及所运用技术原理，本领域技术人员会理解。本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
[0088]
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。