黃　勝，穆　攀，張　睿，袁建國(guó)

(重慶郵電大學(xué) 光通信與網(wǎng)絡(luò)重點(diǎn)實(shí)驗(yàn)室，重慶 400065)

?

基于大衍數(shù)列的規(guī)則QC-LDPC碼構(gòu)造方法

黃勝，穆攀，張睿，袁建國(guó)

(重慶郵電大學(xué) 光通信與網(wǎng)絡(luò)重點(diǎn)實(shí)驗(yàn)室，重慶 400065)

結(jié)合楊輝三角的結(jié)構(gòu)形式，基于大衍數(shù)列提出了一種列重為3或4的規(guī)則準(zhǔn)循環(huán)低密度奇偶校驗(yàn)(QC-LDPC)碼的新構(gòu)造方法。該方法構(gòu)造的校驗(yàn)矩陣圍長(zhǎng)至少為6，碼長(zhǎng)可靈活變化，并且可節(jié)省存儲(chǔ)空間。仿真結(jié)果表明: 在相同的仿真參數(shù)下，當(dāng)誤碼率(BER)為10-6時(shí)，所構(gòu)造的列重為3的QC-LDPC(1260,620)碼的凈編碼增益(NCG)比二次函數(shù)碼改善了1 dB左右；列重為4的QC-LDPC(6056, 3028)碼相對(duì)于WMC-OCS、QC-OCS碼分別有0.1 dB和0.2 dB的NCG提升。

規(guī)則QC-LDPC碼; 楊輝三角; 大衍數(shù)列

低密度奇偶校驗(yàn)碼(LDPC)是一類性能接近香農(nóng)限、糾錯(cuò)能力強(qiáng)的線性分組碼，自從它被重新發(fā)現(xiàn)以來(lái)，LDPC碼的設(shè)計(jì)、構(gòu)造、分析、編譯碼器的實(shí)現(xiàn)以及在數(shù)字通信系統(tǒng)中的應(yīng)用都已經(jīng)成為了人們研究的焦點(diǎn)，成為T(mén)urbo碼的有力競(jìng)爭(zhēng)者。LDPC碼的構(gòu)造方法根據(jù)校驗(yàn)矩陣的結(jié)構(gòu)特點(diǎn)可分為兩種：隨機(jī)構(gòu)造和代數(shù)構(gòu)造。隨機(jī)構(gòu)造法有Gallager構(gòu)造[1]、Mackay構(gòu)造[2]、Hu等人為盡可能增大圍長(zhǎng)提出的經(jīng)典PEG算法[3]等。這些基于計(jì)算機(jī)搜索構(gòu)造出的碼字參數(shù)選擇靈活，性能優(yōu)異。但是，隨機(jī)構(gòu)造所產(chǎn)生的校驗(yàn)矩陣沒(méi)有固定的結(jié)構(gòu)，導(dǎo)致編碼復(fù)雜度高、存儲(chǔ)困難、硬件難以實(shí)現(xiàn)。為了克服隨機(jī)構(gòu)造的這些缺點(diǎn)，代數(shù)的思想被納入到LDPC碼的構(gòu)造中，人們提出了基于有限域[4]、有限幾何[5]、組合數(shù)學(xué)[6]等構(gòu)造法，這些方法構(gòu)造的碼字具有比隨機(jī)碼更低的編碼復(fù)雜度，尤其是具有循環(huán)或準(zhǔn)循環(huán)特性的LDPC碼，已被廣泛應(yīng)用于通信系統(tǒng)中。

近些年來(lái)，人們提出了許多基于組合數(shù)學(xué)構(gòu)造LDPC碼的設(shè)計(jì)方案，比較典型的是BIBD構(gòu)造法[7]，此外，還有很多利用差集、數(shù)列等構(gòu)造LDPC碼的方法[8-9]。根據(jù)大衍數(shù)列固定項(xiàng)差對(duì)應(yīng)的值單調(diào)遞增的性質(zhì)，朱磊基提出了一種基于大衍數(shù)列構(gòu)造LDPC碼的方法[10]，方法簡(jiǎn)單，所構(gòu)造的校驗(yàn)矩陣具有準(zhǔn)循環(huán)結(jié)構(gòu)，長(zhǎng)碼長(zhǎng)的QC-LDPC碼具有優(yōu)異的性能，但此方法只限于構(gòu)造優(yōu)異的長(zhǎng)碼字。相對(duì)于長(zhǎng)碼而言，中短碼的編譯復(fù)雜度低，被廣泛地應(yīng)用于高速率無(wú)線通信及數(shù)據(jù)存儲(chǔ)通信等領(lǐng)域，因此對(duì)中短碼長(zhǎng)的研究具有極好的理論和實(shí)踐意義。本文結(jié)合楊輝三角的構(gòu)造形式，基于大衍數(shù)列提出了一種適宜于構(gòu)造中短碼長(zhǎng)的列重為3或4的規(guī)則QC-LDPC碼的新方法(簡(jiǎn)稱為DY-QC-LDPC碼)。

1　融合大衍數(shù)列和楊輝三角的QC-LDPC碼構(gòu)造

行重為L(zhǎng)、列重為J、碼長(zhǎng)N=P×L的(J,L)QC-LDPC的校驗(yàn)矩陣可以表示為

(1)

式中：IPj,l為循環(huán)置換矩陣；Pj,l表示P×P的單位矩陣右移移位次數(shù)；Pj,l取值范圍為{0,1,…,P-1}或者∞，Pj,l=0時(shí)表示IPj,l為單位矩陣，Pj,l=∞時(shí)表示IPj,l為全零矩陣。

將校驗(yàn)矩陣H中的循環(huán)置換矩陣用循環(huán)移位次數(shù)Pj,l表示，得到參數(shù)移位次數(shù)Pj,l構(gòu)成的基矩陣P，為

(2)

當(dāng)基礎(chǔ)矩陣確定后，QC-LDPC碼的校驗(yàn)矩陣也就唯一確定了。校驗(yàn)矩陣H中長(zhǎng)度為2k的環(huán)可以由基矩陣P中長(zhǎng)為2k的序列(Pj1,l1,Pj1,l2,Pj2,l1,…，Pjk,lk,Pjk,l1)表示，對(duì)于由循環(huán)置換矩陣擴(kuò)展而成的QC-LDPC碼，F(xiàn)ossorier證明了長(zhǎng)為2k的環(huán)存在的充要條件[11]為

(3)

式中：jk=j0，jb+1≠jb，lb+1≠lb，P為維數(shù)。不滿足式(3)則校驗(yàn)矩陣中就不存在長(zhǎng)為2k的環(huán)。

1.1楊輝三角

楊輝三角的構(gòu)造簡(jiǎn)單，是大家所熟知的一種數(shù)表，如圖1所示。

將楊輝三角沿著逆時(shí)針?lè)较蛐D(zhuǎn)45°，可以得到

111121133114641..................

圖1楊輝三角

(4)

矩陣h具有這樣的特點(diǎn)：除第一行與第一列元素以外，其他位置的元素值均為位于本行、上一列的元素與本列、上一行的元素之和。若將這種結(jié)構(gòu)運(yùn)用到校驗(yàn)矩陣中，只需要存儲(chǔ)校驗(yàn)矩陣的第一行和第一列元素，其他元素可以通過(guò)簡(jiǎn)單的加法計(jì)算得到，這樣能夠節(jié)省大量的存儲(chǔ)空間。若將h矩陣定義為校驗(yàn)矩陣的基矩陣，h中的每個(gè)元素都代表著循環(huán)移位的次數(shù)，由fossorier定理可得擴(kuò)展后的校驗(yàn)矩陣中存在大量的四環(huán)，這會(huì)導(dǎo)致譯碼器不能快速收斂甚至不能收斂。所以，如何將楊輝三角這種簡(jiǎn)單的構(gòu)造運(yùn)用到LDPC碼的中，并且使得校驗(yàn)矩陣中不存在四環(huán)是需要考慮的一個(gè)重要問(wèn)題。

1.2DY-QC-LDPC碼構(gòu)造

對(duì)于一個(gè)數(shù)列f(n)，若n取正偶數(shù)，f(n)=(n×n)/2，若n取正奇數(shù)，f(n)=(n×n-1)/2，滿足這樣條件的數(shù)列稱為大衍數(shù)列。結(jié)合楊輝三角的構(gòu)造形式，基于大衍數(shù)列構(gòu)造列重為3或4的規(guī)則QC-LDPC碼的步驟如下：

1)構(gòu)造校驗(yàn)矩陣的基矩陣：將基矩陣的第一行置0，第一列除首元素之外依次用大衍數(shù)列中的奇數(shù)項(xiàng)排列，第二行除首元素之外依次用大衍數(shù)列中的偶數(shù)項(xiàng)排列，其余位置上的元素采用楊輝三角形式構(gòu)造。

2)由Fossorier定理檢測(cè)發(fā)現(xiàn)，基矩陣的第二行的前兩個(gè)元素與第三行的前兩個(gè)元素構(gòu)成了四環(huán)的存在，這里將第二行的首元素0改為1，這樣可以避免四環(huán)的出現(xiàn)。

3)基矩陣擴(kuò)展成校驗(yàn)矩陣，具體方法是：將基矩陣中的0元素用單位矩陣替換，將非零元素用相應(yīng)的移位循環(huán)矩陣替換。

由步驟1)基矩陣的構(gòu)造方法可以看出：編碼器只需要存儲(chǔ)第一、二行與第一列的元素，其他位置元素可以通過(guò)簡(jiǎn)單的加法計(jì)算得到，這樣有利于減少存儲(chǔ)空間。基矩陣中的各位置上的元素值如式(5)，在基矩陣上進(jìn)行擴(kuò)展，可以構(gòu)造出性能優(yōu)異的中短碼長(zhǎng)碼字。

(5)

式(5)中，對(duì)Pj,l=Pj-1,l+Pj,l-1歸納總結(jié)得

(6)

P3,l=l[f(3)+f(2)]+f(5)+(l-lk+1)f(2lk),

2≤lk≤l

(7)

根據(jù)構(gòu)造規(guī)則，式(6)易得。式(7)成立的證明如下：

1)當(dāng)l=1時(shí)，左式=P3,l=l[f(3)+f(2)]，與實(shí)際相符。

綜上，式(7)得證。

2　環(huán)分析

新方法構(gòu)造的基矩陣的行、列上的元素滿足這樣的兩個(gè)特點(diǎn)：1)除第一行以外，每行的元素依次呈遞增關(guān)系；2)每列元素從上到下依次呈遞增關(guān)系。校驗(yàn)矩陣不含四環(huán)，滿足式(8)，其中j0≠j1，l0≠l1。

Pj0,l0-Pj1,l0+Pj1,l1-Pj0,l1≠0 modP

(8)

1)j0位于第0行、j1位于其他行

Pj0,l0-Pj1,l0+Pj1,l1-Pj0,l1=Pj1,l1-Pj1,l0

(9)

2)j0位于第一行、j1位于第二行

P1,l0-P2,l0+P2,l1-P1,l1=f(2l0)-[f(3)+

(10)

3)j0位于第一行、j1位于第三行

P1,l0-P3,l0+P3,l1-P1,l1=f(2l0)-f(2l1)+{l1[f(3)+

(11)

4)j0位于第二行、j1位于第三行

(12)

綜上可知：本文構(gòu)造的(4,L)QC-LDPC碼的校驗(yàn)矩陣中不含四環(huán)，(3,L)QC-LDPC碼中不含四環(huán)的證明與上類似。

3　仿真分析

本節(jié)給出DY-QC-LDPC碼在相同參數(shù)條件下同其他QC-LDPC碼的性能對(duì)比，選取1/2碼率、BP譯碼算法、BPSK調(diào)制方式、最大迭代次數(shù)40次，在AWGN信道下進(jìn)行仿真實(shí)驗(yàn)。首先給出列重為3的DY-QC-LDPC碼的性能，以(3,6)DY-QC-LDPC碼為例，該碼的基矩陣為

(13)

將基矩陣中的0元素用維數(shù)為210×210的單位矩陣替換，其他位置用對(duì)應(yīng)的循環(huán)置換矩陣替換。圖2給出了相同的碼長(zhǎng)、碼率以及仿真環(huán)境下，DY-QC-LDPC(1260,630)碼與大衍數(shù)列碼[10]和兩類二次函數(shù)碼(QF-LDPC)[12]的性能對(duì)比。由圖可以明顯地看出大衍數(shù)列不適合構(gòu)造短碼長(zhǎng)的碼字，DY-QC-LDPC碼性能優(yōu)于二次函數(shù)碼。在誤碼率為10-6時(shí)，DY-QC-LDPC碼比QF-LDPC的兩種碼字性能提高約1 dB左右，具有更好的糾錯(cuò)性能。

圖2　列重為3的DY-QC-LDPC(1 260, 630)碼和其他碼的性能比較

式(10)給出了(4,8)DY-QC-LDPC碼的基礎(chǔ)矩陣，將基礎(chǔ)矩陣中的元素用相應(yīng)的757×757的循環(huán)移位矩陣替換，可以得到碼長(zhǎng)為6 056的DY-QC-LDPC碼。圖3給出了DY-QC-LDPC(6 056, 3 028)碼和相同參數(shù)(碼長(zhǎng)、碼率以及仿真環(huán)境)下的其他碼字性能對(duì)比。在誤碼率為10-6時(shí)，DY-QC-LDPC碼相對(duì)于文獻(xiàn)[13]中所提出的WMC-OCS和QC-OCS碼分別有0.1dB、0.2dB的凈編碼增益，且沒(méi)有明顯的錯(cuò)誤平層。同時(shí)由大衍數(shù)列碼的仿真曲線可以看出文獻(xiàn)[10]中的構(gòu)造方法也不適合于中碼長(zhǎng)碼字的構(gòu)造。

(14)

圖3　列重為4的DY-QC-LDPC(6 056, 3 028)碼和其他碼的性能比較

4　總結(jié)

結(jié)合楊輝三角結(jié)構(gòu)，本文基于大衍數(shù)列構(gòu)造出了一種圍長(zhǎng)至少為6、列重為3或4的規(guī)則DY-QC-LDPC碼。準(zhǔn)循環(huán)特性使得它易于硬件實(shí)現(xiàn)，并且由于基矩陣的構(gòu)造借鑒了楊輝三角的結(jié)構(gòu)形式，所以編碼器只需要存儲(chǔ)0元素以及大衍數(shù)列的項(xiàng)，其他元素可以根據(jù)其位置通過(guò)簡(jiǎn)單的加法計(jì)算得到，這樣可節(jié)省存儲(chǔ)空間。在相同的參數(shù)條件下，與大衍數(shù)列碼、QF-LDPC、WMC-OCS和QC-OCS碼相比，仿真結(jié)果表明中短長(zhǎng)度的DY-QC-LDPC碼具有優(yōu)異的糾錯(cuò)性能。

[1]GALLAGERG.Low-densityparity-checkcodes[J].IEEEtransactionsoninformationtheory, 1962, 8(1)：21-28.

[2]MACKAYDJC.Gooderror-correctingcodesbasedonverysparsematrices[J].IEEEtransactionsoninformationtheory, 1999, 45(2)：399-431.

[3]HUXY,ELEFTHERIOUE,ARNOLDD.M.Regularandirregularprogressiveedge-growthtannergraphs[J].IEEEtransactionsoninformationtheory, 2005, 51(1)：386-398.

[4]HUANGS,TIANFF,JIAXT.AnovelconstructionmethodofLDPCcodesoverfinitefieldintheopticalcommunicationsystem[J].Actaphotonicasinica, 2014,43(6)：1-4.

[5]HUANGQ,DIAOQJ,LINS.CyclicandQuasi-CyclicLDPCcodesonconstrainedparity-checkmatricesandtheirtrappingsets[J].IEEEtransactionsoninformationtheory, 2012, 58(5):2648-2671.

[6]YUANJG,LICY,HUANGS,etal.AnovelconstructionmethodofQC-LDPCcodesforopticalcommunicationsbasedontheBIBDandcirculantdecomposition[J].Journalofoptoelectronics·laser, 2013, 24(9)：1698-1701.

[7]LAN L, TAI Y Y, LIN S. New constructions of quasi-cyclic LDPC codes based on special classes of BIBD's for the AWGN and binary erasure channels[J]. IEEE transactions on communications, 2008, 56(1)：39-48.

[8]ESMAEILI M. 4-Cycle free LDPC codes based on difference sets[J]. IEEE transactions on communications, 2012, 60(12)：3579-3586.

[9]ZHANG L J, LI B. Constructions of QC LDPC codes based on integer sequences[J]. Science China, 2014, 57(6)：1-14.

[10]ZHU L J, WANG H, SHI Y S. Method for constructing QC-LDPC codes using the dayan sequence[J]. Journal of Xidian University, 2012, 39(3)：144-160.

[11]FOSSORIER M P C. Quasi-cyclic low-density parity-check codes from circulant permutation matrices[J]. IEEE transactions on information theory, 2004, 50(8)：1788-

1793.

[12]ZHANG G H, SUN R. Deterministic construction of girth-eight (3, L) QC-LDPC codes from quadratic function[J]. Electronics letters, 2013,49(9)：600-602.

[13]HUANG J F, HUANG C M, YANG C C. Construction of one-coincidence sequence quasi-cyclic LDPC codes of large girth[J]. IEEE transactions on information theory, 2012, 58(3)：1825-1836.

黃勝(1974— )，教授，主要研究方向?yàn)樾诺谰幋a和下一代網(wǎng)絡(luò)；

穆攀(1990— )，女，碩士生，主研信道編碼；

張睿(1991— )，女，碩士生，主研信道編碼；

袁建國(guó)(1968— )，教授，主要研究方向?yàn)楣馔ㄐ偶夹g(shù)與光電子技術(shù)。

責(zé)任編輯：薛京

Method of construction of regular QC-LDPC codes based on Dayan sequence

HUANG Sheng, MU Pan, ZHANG Rui, YUAN Jianguo

(KeyLaboratoryofOpticalCommunicationandNetwork,ChongqingUniversityofPostsandTelecommunications,Chongqing400065，China)

According to structure of Yanghui Triangle and based on the Dayan sequence, a novel construction method of regular quasi-cyclic low-density parity-check (QC-LDPC) codes whose column weight is 3 or 4 is put forward. Girth of parity-check matrix is at least 6. Length of the code can vary flexibly. And storage spaces can be saved. Under the same parameters, at the bit error rate(BER) of 10-6, simulation results show QC-LDPC(1260, 620)code with column weight 3 has a net coding gain(NCG) nearly 1 dB more than quadratic function code and the NCG of QC-LDPC(6056, 3028) code with column weight 4 is 0.1 dB and 0.2 dB more than WMC-OCS and QC-OCS code respectively.

regular QC-LDPC code; Yanghui triangle; Dayan sequence

TN911

A

10.16280/j.videoe.2016.09.015

國(guó)家自然科學(xué)基金項(xiàng)目(61171158)；重慶市自然科學(xué)基金項(xiàng)目(cstc2013jcyjA40052;cstc2012jjA40060);重慶市教委科學(xué)技術(shù)研究項(xiàng)目(KJ130515)

2016-01-08

文獻(xiàn)引用格式：黃勝，穆攀，張睿，等. 基于大衍數(shù)列的規(guī)則QC-LDPC碼構(gòu)造方法[J]. 電視技術(shù)，2016,40(9)：77-80.

HUANG S, MU P, ZHANG R，et al. Method of construction of regular QC-LDPC codes based on Dayan sequence [J]. Video engineering，2016,40(9)：77-80.