CN108768593B - 一种基于dmt调制的qc-ldpc码编、解码方法与系统 - Google Patents

一种基于dmt调制的qc-ldpc码编、解码方法与系统 Download PDF

Info

Publication number
CN108768593B
CN108768593B CN201810604042.4A CN201810604042A CN108768593B CN 108768593 B CN108768593 B CN 108768593B CN 201810604042 A CN201810604042 A CN 201810604042A CN 108768593 B CN108768593 B CN 108768593B
Authority
CN
China
Prior art keywords
check
node
variable
decoding
nodes
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201810604042.4A
Other languages
English (en)
Other versions
CN108768593A (zh
Inventor
王瑾
梁晴晴
曾福江
张亚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences
Original Assignee
China University of Geosciences
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences filed Critical China University of Geosciences
Priority to CN201810604042.4A priority Critical patent/CN108768593B/zh
Publication of CN108768593A publication Critical patent/CN108768593A/zh
Application granted granted Critical
Publication of CN108768593B publication Critical patent/CN108768593B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0056—Systems characterized by the type of code used
    • H04L1/0057—Block codes
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/11—Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
    • H04B10/114—Indoor or close-range type systems
    • H04B10/116—Visible light communication
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0047—Decoding adapted to other signal detection operation
    • H04L1/005—Iterative decoding, including iteration between signal detection and decoding operation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0047—Decoding adapted to other signal detection operation
    • H04L1/005—Iterative decoding, including iteration between signal detection and decoding operation
    • H04L1/0051—Stopping criteria
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0056—Systems characterized by the type of code used
    • H04L1/0061—Error detection codes
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03171—Arrangements involving maximum a posteriori probability [MAP] detection
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03178—Arrangements involving sequence estimation techniques
    • H04L25/03312—Arrangements specific to the provision of output signals
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2602—Signal structure
    • H04L27/2605—Symbol extensions, e.g. Zero Tail, Unique Word [UW]
    • H04L27/2607—Cyclic extensions
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2626—Arrangements specific to the transmitter only
    • H04L27/2627—Modulators
    • H04L27/2628—Inverse Fourier transform modulators, e.g. inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Error Detection And Correction (AREA)

Abstract

本发明公开了一种基于DMT调制的QC‑LDPC码编、解码方法与系统，在进行编码过程中，将信源产生的随机信息比特与生成矩阵相乘得到编码后的码字；其中，生成矩阵对应的校验矩阵为非正则QC‑LDPC码，解码时，先进行初始化，然后进行迭代处理直至最终满足终止条件，并输出译码结果。利用本发明的QC‑LDPC码进行编码，解码简单，且降低了VLC系统的ISI，进而提高系统的频谱利用率和数据传输速率，同时提高了系统传输信息的可靠性以及精确性。

Description

一种基于DMT调制的QC-LDPC码编、解码方法与系统
技术领域
本发明涉及信道编解码领域，更具体地说，涉及一种基于DMT调制的QC-LDPC码编、解码方法与系统。
背景技术
可见光通信(Visible Light Communication，VLC)作为一项近几年兴起的新型无线接入技术，以其功耗低、频带宽、保密性强和无电磁污染等优点，得到了世界范围内的广泛关注。在传统的VLC系统中，如图1所示，首先发射机将电信号转换为光信号通过LED1-4发送出去，随后接收端光电二极管(photodiode，PD)将接收光功率转换为电流。在VLC系统中，LED可以同时进行数据传输以及照明。
在可见光通信系统中，DMT调制可以有效地利用调制带宽实现较高的频谱利用率，同时对由多径传播引起的ISI具有适应性。基于快速傅里叶变换(fast Fouriertransform，FFT)算法，可以通过DMT来实现多载波调制解调技术。与正交频分复用(orthogonal frequency-division multiplexing，OFDM)相反，DMT调制器经过逆FFT(inverse FFT，IFFT)后，其输出端的信号是一个实值信号。因此必须寻找一种低错误平层，高速率通信和性能优异的VLC系统信道编码方案。在该背景下，一些简单的分组码比如Reed-Solomon(RS)码或者具有较强纠错能力的LDPC码和turbo码就可以应用到VLC系统中。LDPC码是一类可以用稀疏校验矩阵H(Parity Check Matrix H)或二分图(BipartiteGraph)来描述的线性分组纠错码。LDPC码的最小汉明距离会随着码长的增加而线性增加。在进行后验概率迭代译码时，随着码字长度的增加，系统的BER会随之降低；而且采用迭代译码算法时LDPC码的性能非常接近香农极限。因此可以利用随机结构或者代数方法来构造LDPC码，但是当码长较短时，随机构造的LDPC码，其性能往往不是很理想。而且在这种情况下比较容易出现较高的误码平层，换句话说就是随机构造的码会使得系统在译码时比较困难。
发明内容
本发明要解决的技术问题在于，针对现有技术中随机构造的LDPC码会使得系统在译码时比较困难的技术缺陷，提供一种基于DMT调制的QC-LDPC码编解码方案。
根据本发明的第一方面，本发明解决其技术问题所采用的技术方案是：构造一种基于DMT调制的QC-LDPC码编码方法，在进行编码过程中，将信源产生的随机信息比特与生成矩阵相乘得到编码后的码字；其中，所述生成矩阵对应的校验矩阵为：
Figure GDA0002375049550000021
式中，H表示校验矩阵，I是大小为L×L的单位矩阵，P是经过循环移位之后的单位矩阵，P的上标表示循环移位的次数，该校验矩阵H的基矩阵的大小为m×n，并且L是质数满足L＞n＞m。
根据本发明的又一方面，本发明为解决其技术问题，还提供了一种基于DMT调制的QC-LDPC码解码方法，用于对上述的基于DMT调制的QC-LDPC码编码方法编码后的信息进行解码，包括如下步骤：
S1、初始化：计算得出信道传递给各个变量节点i的初始概率似然比消息L(Pi)；根据下述公式(1)初始化变量节点i传递给其相邻的校验节点j的初始消息L(0)(qij)：
L(0)(qij)＝L(Pi) (1)；
式中，i＝1、2、…、N，j＝1、2、…、M，N＝L*n，M＝L*m；
S2、迭代处理；其中，在进行第任意迭代次数l时包括如下依次进行的步骤：
S21、校验节点信息处理：
对所有的校验节点j和与其相邻的变量节点i∈R(j)在第l次迭代时，根据公式(2)计算出变量节点传递给校验节点的消息得L(l)(rji)：
Figure GDA0002375049550000022
L(l-1)(qi’j)表示对所有的变量节点i’和与其相邻的校验节点j∈C(i)在第l-1次迭代时，校验节点传递给变量节点的消息，qi’j代表的是从变量节点i’传递给校验节点j的外部概率信息；
S22、变量节点信息处理：
对所有的变量节点i和与其相邻的校验节点j∈C(i)在第l次迭代时，根据公式(3)计算出校验节点传递给变量节点的消息L(l)(qij)：
Figure GDA0002375049550000031
rj’i代表的是从校验节点j’传递给变量节点i的外部概率信息；
S23、译码判决：根据公式(4)计算得出所有变量节点的硬判决信息L(l)(qi)：
Figure GDA0002375049550000032
然后根据硬判决信息L(l)(qi)得出译码结果
Figure GDA0002375049550000033
若L(l)(qi)＞0，则
Figure GDA0002375049550000034
否则
Figure GDA0002375049550000035
S24、迭代停止判断：若满足
Figure GDA0002375049550000036
或者已经达到预先设定的最大迭代次数，那么就停止译码并以最后一次译码的译码结果作为最终的译码结果进行输出，否则返回步骤S21继续迭代；式中，
Figure GDA0002375049550000037
为本次的译码结果的转置；
其中，所述QC-LDPC码为二进制LDPC码，rji(b)代表的是从校验节点j传递给变量节点i的外部概率信息，qij(b)代表的是从变量节点i传递给校验节点j的外部概率信息，C(i)代表的是所有与变量节点i相连的校验节点的集合，R(j)代表的是所有与校验节点j相连的变量节点的集合，C(i)\j代表的是除校验节点j外所有与变量节点i相连的校验节点的集合，R(j)\i代表的是除变量节点i外所有与校验节点j相连的变量节点的集合，b表示二进制码0和1。
进一步地，在本发明的基于DMT调制的QC-LDPC码解码方法中，步骤S24中，译码结果c等于[c1 c2 … cN]，c1至cN为第1至N个译码输出。
根据本发明的再一方面，一种基于DMT调制的QC-LDPC码编码系统，在进行编码过程中，将信源产生的随机信息比特与生成矩阵相乘得到编码后的码字；其中，所述生成矩阵对应的校验矩阵为：
Figure GDA0002375049550000038
式中，H表示校验矩阵，I是大小为L×L的单位矩阵，P是经过循环移位之后的单位矩阵，P的上标表示循环移位的次数，该校验矩阵H的基矩阵的大小为m×n，并且L是质数满足L＞n＞m。
本发明为解决其技术问题，还提供一种基于DMT调制的QC-LDPC码解码系统，用于对上述的基于DMT调制的QC-LDPC码编码系统编码后的信息解码，包括如下模块：
初始化模块，用于进行初始化：计算得出信道传递给各个变量节点i的初始概率似然比消息L(Pi)；根据下述公式(1)初始化变量节点i传递给其相邻的校验节点j的初始消息L(0)(qij)：
L(0)(qij)＝L(Pi) (1)；
式中，i＝1、2、…、N，j＝1、2、…、M，N＝L*n，M＝L*m；
迭代处理模块，用于进行迭代处理；其中，在进行第任意迭代次数l时采用如下子模块依次进行处理：
校验节点信息处理子模块，用于对所有的校验节点j和与其相邻的变量节点i∈R(j)在第l次迭代时，根据公式(2)计算出变量节点传递给校验节点的消息得L(l)(rji)：
Figure GDA0002375049550000041
L(l-1)(qi’j)表示对所有的变量节点i’和与其相邻的校验节点j∈C(i)在第l-1次迭代时，校验节点传递给变量节点的消息，qi’j代表的是从变量节点i’传递给校验节点j的外部概率信息；
变量节点信息处理子模块，用于对所有的变量节点i和与其相邻的校验节点j∈C(i)在第l次迭代时，根据公式(3)计算出校验节点传递给变量节点的消息L(l)(qij)：
Figure GDA0002375049550000042
译码判决子模块，用于根据公式(4)计算得出所有变量节点的硬判决信息L(l)(qi)：
Figure GDA0002375049550000043
然后根据硬判决信息L(l)(qi)得出译码结果
Figure GDA0002375049550000044
若L(l)(qi)＞0，则
Figure GDA0002375049550000045
否则
Figure GDA0002375049550000046
迭代停止判断子模块，用于若满足
Figure GDA0002375049550000051
或者已经达到预先设定的最大迭代次数，那么就停止译码并以最后一次译码的译码结果作为最终的译码结果进行输出，否则返回步骤S21继续迭代；式中，
Figure GDA0002375049550000052
为本次的译码结果的转置；
其中，所述QC-LDPC码为二进制LDPC码，rji(b)代表的是从校验节点j传递给变量节点i的外部概率信息，qij(b)代表的是从变量节点i传递给校验节点j的外部概率信息，C(i)代表的是所有与变量节点i相连的校验节点的集合，R(j)代表的是所有与校验节点j相连的变量节点的集合，C(i)\j代表的是除校验节点j外所有与变量节点i相连的校验节点的集合，R(j)\i代表的是除变量节点i外所有与校验节点j相连的变量节点的集合，b表示二进制码0和1。
进一步地，在本发明的基于DMT调制的QC-LDPC码解码系统中，迭代停止判断子模块中，译码结果c等于[c1 c2 … cN]，c1至cN为第1至N个译码输出。
实施本发明的一种基于DMT调制的QC-LDPC码编、解码方法与系统，具有以下有益效果：利用本发明的QC-LDPC码进行编码，解码简单，且降低了VLC系统的ISI，进而提高系统的频谱利用率和数据传输速率，同时提高了系统传输信息的可靠性以及精确性。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是传统的可见光通信系统模型；
图2是基于DMT调制的QC-LDPC码解码方法流程图；
图3是基于DMT调制的VLC系统的结构框图；
图4是VLC信道和AWGN信道情况下BER的对比图；
图5是不同长度循环前缀情况下DMT-VLC系统BER的对比图；
图6是不同码率情况下DMT-VLC系统BER的对比图；
图7是不同迭代次数情况下DMT-VLC系统BER的对比图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
在本发明中提出基于VLC系统的DMT调制技术用来降低ISI，同时提出将QC-LDPC码应用于DMT-VLC系统中用于有效编码。
在本发明中，主要利用QC-LDPC编码的DMT调制系统，通过QC-LDPC码对要传输的信息进行编码，同时利用基于对数似然比的置信传播(logarithmic likelihood ratiobelief propagation，LLR-BP)译码算法对输出端的信号进行译码。并且分析两种不同信道下QC-LDPC编码的DMT调制系统的性能。
VLC系统架构介绍如下。
A.信道分析
一个传统的光无线系统比如常见的基于LED的室内VLC系统(如图1所示)，可以用下式表示
Figure GDA0002375049550000061
其中y(t)和x(t)分别表示接收信号和发射信号的幅度，R表示接收端的光电转换系数，h(t)表示可见光信道的脉冲响应，n(t)表示加性高斯白噪声(additive whiteGaussian noise，AWGN)，假设该噪声是独立同分布的。
由于循环前缀(cyclic prefix，CP)的存在，DMT调制对于由多径效应引入的散射信号具有强鲁棒性。反射光功率可以忽略不计，只考虑直射链路。直射信道的光信道衰落根据LED的朗伯辐射模式可以用下式来模拟：
Figure GDA0002375049550000062
其中，朗伯模型的阶数m可以用下式表示m＝-ln2/ln(cosΦ1/2)，其中Φ1/2为LED的半功率角。A表示光电探测器的物理面积，d表示发射端与接收端的距离，φ表示发射光与发射光法线方向的夹角，必表示入射光与接收器平面之间的夹角，则在接收端，每个子载波在时域内的输出信号的幅度可以用下式表示：
Figure GDA0002375049550000063
B.码的构造
目前，可见光通信系统的研究主要集中在调制技术、信道均衡技术以及快速传输速率等方面。随着可见光通信关键技术和芯片的快速研发，未来的可见光通信将朝着超宽带，超高速，大覆盖范围和高可靠性的方向发展，包括移动大数据和高清影像等对误码性能要求高的多媒体业务。因此必须寻找一种低错误平层，高速率通信和性能优异的VLC系统信道编码方案。在该背景下，一些简单的分组码比如Reed-Solomon(RS)码或者具有较强纠错能力的LDPC码和turbo码就可以应用到VLC系统中。LDPC码是一类可以用稀疏校验矩阵H(Parity Check Matrix H)或二分图(Bipartite Graph)来描述的线性分组纠错码。LDPC码的最小汉明距离会随着码长的增加而线性增加。在进行后验概率迭代译码时，随着码字长度的增加，系统的BER会随之降低；而且采用迭代译码算法时LDPC码的性能非常接近香农极限。因此可以利用随机结构或者代数方法来构造LDPC码，但是当码长较短时，随机构造的LDPC码，其性能往往不是很理想。而且在这种情况下比较容易出现较高的误码平层。换句话说就是随机构造的码会使得系统在译码时比较困难。为了解决这个问题，本发明提出一种QC-LDPC码对系统的传输信息进行有效编码，其校验矩阵如下修正数组阵列所示：
Figure GDA0002375049550000071
这是一种非正则QC-LDPC码，I是单位矩阵，其大小为L×L(×之前为行数，之后为列数，其他位置与此相同)，P是经过循环移位之后的单位矩阵。该校验矩阵基矩阵的大小为m×n，则校验矩阵H的大小为Lm×Ln。其中L＝i，m＝j，n＝k，并且i是质数满足i＞k＞j。很明显校验矩阵H(i，j，k)是一个满秩矩阵，因为它是一个上三角矩阵，对角线上的元素是非零元素。因此它的码率可以表示为
Figure GDA0002375049550000072
因为校验矩阵H(i，j，k)是一种上三角的形式，所以它能够进行有效地编码。并且它能够很容易的校验因为它所对应的二分图中不存在4环。
QC-LDPC码在进行编码的过程中，将信源产生的随机信息比特与生成矩阵G相乘，得到编码后的码字，即
c＝s·G (5)
其中，s表示随机信息比特，c表示编码后的码字，G表示生成矩阵。任何一个生成矩阵都能够被置换简化为：
G＝[I0|Q0] (6)
这种固定形式的矩阵被称之为系统的生成矩阵。其中I0为单位矩阵。通常都是对LDPC码的校验矩阵进行LU分解或者高斯消元法来得到生成矩阵。对校验矩阵进行置换，从而使其转换成：H1＝[P1|I1]，其中对P1进行转置可以得到Q0，即Q0＝P’1。
根据系统传输信息的表达式，置信传播(beliefpropagation，BP)译码算法可以分为基于概率的BP译码算法和基于对数似然比的BP译码算法。在本发明中，使用LLR-BP译码算法对系统接收端的信息进行译码，LLR-BP译码算法节点间传递的信息是对数似然比值，将大量运算变为加法运算，在很大程度上减少了运算时间，加快了译码速度。参考图2，其译码步骤具体如下：
本发明采用二进制LDPC码，首先定义以下变量：rji(b)(b＝0，1)代表的是从校验节点j传递给变量节点i的外部概率信息，qij(b)代表的是从变量节点i传递给校验节点j的外部概率信息，C(i)代表的是所有与变量节点i相连的校验节点的集合，R(j)代表的是所有与校验节点j相连的变量节点的集合，C(i)\j代表的是除校验节点j外所有与变量节点i相连的校验节点的集合，R(j)\i代表的是除变量节点i外所有与校验节点j相连的变量节点的集合。
S1、初始化
首先要计算得出信道传递给各个变量节点i的初始概率似然比消息L(Pi)，i＝1、2、…、n；然后对每一变量节点i和与其相邻的校验节点j∈C(i)，根据式(7)初始化变量节点传递给校验节点的初始消息L(0)(qij)：
L(0)(qij)＝L(Pi) (7)
S2、迭代处理；其中，在进行第l次迭代时包括如下依次进行的步骤：
S21：校验节点信息处理
也可称作横向步骤，对所有的校验节点j和与其相邻的变量节点i∈R(j)在第l次迭代时，根据式(8)或(9)计算出变量节点传递给校验节点的消息L(l)(rji)：
Figure GDA0002375049550000091
或
Figure GDA0002375049550000092
L(l-1)(qi’j)表示对所有的变量节点i’和与其相邻的校验节点j∈C(i)在第l-1次迭代时，校验节点传递给变量节点的消息，qi’j代表的是从变量节点i’传递给校验节点j的外部概率信息；
S22：变量节点信息处理
也可称作纵向步骤，对所有的变量节点i和与其相邻的校验节点j∈C(i)在第l次迭代时，根据式(10)计算出校验节点传递给变量节点的消息L(l)(qij)：
Figure GDA0002375049550000093
rj’i代表的是从校验节点j’传递给变量节点i的外部概率信息；S23：译码判决
计算得出所有变量节点的硬判决信息为
Figure GDA0002375049550000094
若L(l)(qi)＞0，则
Figure GDA0002375049550000095
否则
Figure GDA0002375049550000096
S24：迭代停止判决
若满足
Figure GDA0002375049550000097
或者已经达到预先设定的最大迭代次数，那么就停止译码并以最后一次译码的译码结果作为最终的译码结果进行输出，否则返回S21继续迭代。其中
Figure GDA0002375049550000098
等于[c1 c2 … cN]，c1至cN为第1至N个译码输出。
C.DMT调制分析
如上所述，DMT调制是一种特殊的多载波调制技术，能够把一个高速串行数据流划分为多个并行低速数据流并且调制到不同频率的子载波上，而这些子载波之间要满足相互正交，这样就可以保证每个子载波在信道中传送数据时不会相互干扰，从而降低多径衰落的影响同时有效的抑制由多径效应引起的ISI问题，进而提高系统的频谱利用率和数据传输速率。
基于DMT调制的VLC系统其结构框图如图3所示。输入的高速串行数据流经过串并转换变成Q组并行低速二进制数据流。每组比特流经过正交幅度调制(quadratureamplitude modulation，QAM)映射到星座图上。经过QAM星座映射后的数据为复值信号，可以表示为Cq＝Aq+jBq，其中q＝0，1，…，Q-1。然后，这些信号在DMT发射端经过快速傅里叶变换(fast Fourier transform，FFT)调制到相应的载波上。基带DMT调制的传输序列包含有Q个子载波，为了获得一个实值信号，则需要2Q个IFFT点。对于IFFT模块的2Q个输入信号，前面一半信号赋值为Cq，后面一半赋值为Cq的共轭，满足下面的公式：
C2Q-q-1＝Cq * (12)
其中q＝0，1，…Q-1。通常情况下，将第一个信号以及第Q个信号值设为0，即C0＝CQ＝0，那么最终的DMT序列中不包含任何直流信号。则2Q个IFFT点的输出端的信号为实值信号，这些实值信号可以表示为：
Figure GDA0002375049550000101
其中g＝0，1，…2Q-1。
在每2Q个IFFT点的分组中加入长度为P的循环前缀用来降低ISI的影响。在实际情况中，循环前缀一般是复制信号U(g)的后面一部分然后插入到信号U(g)的前端。那么此时信号可以表示为
Figure GDA0002375049550000102
其中g＝0，1，…2Q+P-1。
在IFFT模块输出端通过并串转换可以得到每个长度为(Q+P)的单位序列所对应的DMT时域信号的离散时间采样值。随后，输出的串行数据经过数模转换传输到信道中。在接收过程中经过光敏二极管，放大器，滤波器后的信号可以通过模数转换恢复到二进制数据。假设不考虑模数转换之前数模转换器的脉冲响应以及滤波器响应，接收端接收到的DMT帧信号(模数转换之前)可以表示为：
Figure GDA0002375049550000103
经过串并转换和删除循环前缀后的信号传输到FFT模块。理想情况下，此时FFT模块输出端的信号应该是从发射端发送到IFFT模块的原始信号。另外，利用信道的均衡处理达到消除或减小由发射机和接收机中有源或无源设备的参数变化引起的码间干扰的目的。最终信号进行解映射，接着通过并串转换，信号完成在DMT-VLC系统中的传输。
实验结果
在这一部分，主要研究循环前缀的长度，迭代次数以及码率对DMT-VLC系统性能的影响。本发明利用在不同信噪比情况下接收信号的误码率作为接收信号质量的判决标准。同时，利用LLR-BP译码算法，研究QC-LDPC码在两种不同信道下的性能。
根据在上述构造的QC-LDPC码，其中构造出了大小为Lm×Ln的校验矩阵，其中L＝31，m＝50，n＝70。接收端的信号可以表示为：
1.当发射端的信号为0时，接收信号的条件概率密度函数可以表示为：
Figure GDA0002375049550000111
2.当发射端的信号为1时，接收信号的条件概率密度函数可以表示为：
Figure GDA0002375049550000112
其中H0和H1分别表示发射信号为0和1。σ2表示对数振幅方差。I表示光强。因为选择的是二进制QC-LDPC码，所以采用LLR-BP算法对经过QC-LDPC编码的光信号进行译码。
在两种不同信道情况下，即VLC信道和AWGN信道，对QC-LDPC码的编码性能进行仿真。实验结果如4所示，可以看出，在没有进行QC-LDPC编码时，VLC信道的性能明显比AWGN信道的性能要差。这是因为相较于AWGN信道，在VLC信道中存在更多的背景噪声，ISI以及多径干扰。因此，在VLC信道中有必要采用QC-LDPC编码来提高VLC系统的性能。正如预期的那样，从图4中可以看出，相比于没有采用QC-LDPC编码的系统，利用QC-LDPC编码的VLC系统能够取得较高的编码增益。
从图5可以看出，当循环前缀的长度设为10时，系统的可靠性明显提高。但是当循环前缀的长度增加到30时，系统的BER反而增加了。根据仿真结果以及数值分析，可以得出在VLC系统中，DMT调制模块中循环前缀的设置对系统抗干扰能力的提高占有至关重要的地位。但是，循环前缀过长会导致系统硬件复杂度提高，进而影响系统的数据传输速率，同时会占用较多的硬件资源。因此，对于不同的信道环境，需要选择一个合适长度的循环前缀，从而确保系统处于最优性能的状态，同时避免由DMT调制所引起的硬件复杂度过高的问题。
作为LDPC码的重要参数，码率的选择对于LDPC码的编码性能的影响较大。图6为四种相同码长但不同码率的LDPC码在DMT-VLC系统下的性能表现。根据仿真结果可以看出，随着码率的增大，QC-LDPC码的性能有所降低，这是因为当码率增大时，码中的校验位逐渐减少，码的纠错能力也随之减弱。因此选择一个合适的码率才能得到最佳编码性能。对于本发明所提出的DMT-VLC系统中，取得最优性能时的码率应该设为0.4。
迭代次数是影响译码质量的重要因素。在DMT-VLC系统中，对不同迭代次数的译码算法进行仿真，检验迭代次数对译码性能的影响，仿真结果如图7所示。从图7可以看出，随着迭代次数的增加，系统传输数据的精确度越高。但是当迭代次数增加到某一数值时，系统的性能就不再提高。考虑到译码迭代次数越大，系统的计算复杂度就越高，不能进行有效译码，因此在本发明中将译码迭代次数设置为50。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。

Claims (5)

1.一种基于DMT调制的QC-LDPC码编码方法，其特征在于，在进行编码过程中，将信源产生的随机信息比特与生成矩阵相乘得到编码后的码字；其中，所述生成矩阵对应的校验矩阵为：
Figure FDA0002457377200000011
式中，H表示校验矩阵，I是大小为L×L的单位矩阵，P是经过循环移位之后的单位矩阵，P的上标表示循环移位的次数，该校验矩阵H的基矩阵的大小为m×n，并且L是质数满足L＞n＞m。
2.一种基于DMT调制的QC-LDPC码解码方法，用于对权利要求1所述的基于DMT调制的QC-LDPC码编码方法编码后的信息解码，其特征在于，包括如下步骤：
S1、初始化：计算得出信道传递给各个变量节点i的初始概率似然比消息L(Pi)；根据下述公式(1)初始化变量节点i传递给其相邻的校验节点j的初始消息L(0)(qij)：
L(0)(qij)＝L(Pi) (1)；
式中，i＝1、2、…、N，j＝1、2、…、M，N＝L*n，M＝L*m；
S2、迭代处理；其中，在进行第任意迭代次数l时包括如下依次进行的步骤：
S21、校验节点信息处理：
对所有的校验节点j和与其相邻的变量节点i∈R(j)在第l次迭代时，根据公式(2)计算出变量节点传递给校验节点的消息得L(l)(rji)：
Figure FDA0002457377200000012
L(l-1)(qi′j)表示对所有的变量节点i’和与其相邻的校验节点j∈C(i)在l-1次迭代时，校验节点传递给变量节点的消息，qi′j代表的是从变量节i’传递给校验节点j的外部概率信息；
S22、变量节点信息处理：
对所有的变量节点i和与其相邻的校验节点j∈C(i)在第l次迭代时，根据公式(3)计算出校验节点传递给变量节点的消息L(l)(qij)：
Figure FDA0002457377200000021
rj′i代表的是从校验节点j’传递给变量节点i的外部概率信息；
S23、译码判决：根据公式(4)计算得出所有变量节点的硬判决信息L(1)(qi)：
Figure FDA0002457377200000022
然后根据硬判决信息L(l)(qi)得出译码结果
Figure FDA0002457377200000023
若L(l)(qi)＞0，则
Figure FDA0002457377200000024
否则
Figure FDA0002457377200000025
S24、迭代停止判断：若满足
Figure FDA0002457377200000026
或者已经达到预先设定的最大迭代次数，那么就停止译码并以最后一次译码的译码结果作为最终的译码结果进行输出，否则返回步骤S21继续迭代；式中，
Figure FDA0002457377200000027
为本次的译码结果的转置；
其中，所述QC-LDPC码为二进制LDPC码，rji(b)代表的是从校验节点j传递给变量节点i的外部概率信息，qij(b)代表的是从变量节点i传递给校验节点j的外部概率信息，C(i)代表的是所有与变量节点i相连的校验节点的集合，R(j)代表的是所有与校验节点j相连的变量节点的集合，C(i)\j代表的是除校验节点j外所有与变量节点i相连的校验节点的集合，R(j)\i代表的是除变量节点i外所有与校验节点j相连的变量节点的集合，b表示二进制码0和1。
3.根据权利要求2所述的基于DMT调制的QC-LDPC码解码方法，其特征在于，步骤S24中，译码结果
Figure FDA0002457377200000028
等于[c1 c2…cN]，c1至cN为第1至N个译码输出。
4.一种基于DMT调制的QC-LDPC码解码系统，用于对基于DMT调制的QC-LDPC码编码系统编码后的信息解码，所述基于DMT调制的QC-LDPC码编码系统在进行编码过程中，将信源产生的随机信息比特与生成矩阵相乘得到编码后的码字；其中，所述生成矩阵对应的校验矩阵为：
Figure FDA0002457377200000031
式中，H表示校验矩阵，I是大小为L×L的单位矩阵，P是经过循环移位之后的单位矩阵，P的上标表示循环移位的次数，该校验矩阵H的基矩阵的大小为m×n，并且L是质数满足L＞n＞m；
其特征在于，所述基于DMT调制的QC-LDPC码解码系统包括如下模块：
初始化模块，用于进行初始化：计算得出信道传递给各个变量节点i的初始概率似然比消息L(Pi)；根据下述公式(1)初始化变量节点i传递给其相邻的校验节点j的初始消息L(0)(qij)：
L(0)(qij)＝L(Pi) (1)；
式中，i＝1、2、…、N，j＝1、2、…、M，N＝L*n，M＝L*m；
迭代处理模块，用于进行迭代处理；其中，在进行第任意迭代次数l时采用如下子模块依次进行处理：
校验节点信息处理子模块，用于对所有的校验节点j和与其相邻的变量节点i∈R(j)在第l次迭代时，根据公式(2)计算出变量节点传递给校验节点的消息得L(l)(rji)：
Figure FDA0002457377200000032
L(l-1)(qi′j)表示对所有的变量节点i’和与其相邻的校验节点j∈C(i)在l-1次迭代时，校验节点传递给变量节点的消息，qi′j代表的是从变量节i’传递给校验节点j的外部概率信息；
变量节点信息处理子模块，用于对所有的变量节点i和与其相邻的校验节点j∈C(i)在第l次迭代时，根据公式(3)计算出校验节点传递给变量节点的消息L(l)(qij)：
Figure FDA0002457377200000033
rj′i代表的是从校验节点j’传递给变量节点i的外部概率信息；
译码判决子模块，用于根据公式(4)计算得出所有变量节点的硬判决信息L(l)(qi)：
Figure FDA0002457377200000041
然后根据硬判决信息L(l)(qi)得出译码结果
Figure FDA0002457377200000042
若L(l)(qi)＞0，则
Figure FDA0002457377200000043
否则
Figure FDA0002457377200000044
迭代停止判断子模块，用于若满足
Figure FDA0002457377200000045
或者已经达到预先设定的最大迭代次数，那么就停止译码并以最后一次译码的译码结果作为最终的译码结果进行输出，否则返回步骤S21继续迭代；式中，
Figure FDA0002457377200000046
为本次的译码结果的转置；
其中，所述QC-LDPC码为二进制LDPC码，rji(b)代表的是从校验节点j传递给变量节点i的外部概率信息，qij(b)代表的是从变量节点i传递给校验节点j的外部概率信息，C(i)代表的是所有与变量节点i相连的校验节点的集合，R(j)代表的是所有与校验节点j相连的变量节点的集合，C(i)\j代表的是除校验节点j外所有与变量节点i相连的校验节点的集合，R(j)\i代表的是除变量节点i外所有与校验节点j相连的变量节点的集合，b表示二进制码0和1。
5.根据权利要求4所述的基于DMT调制的QC-LDPC码解码系统，其特征在于，迭代停止判断子模块中，译码结果
Figure FDA0002457377200000047
等于[c1 c2…cN]，c1至cN为第1至N个译码输出。
CN201810604042.4A 2018-06-12 2018-06-12 一种基于dmt调制的qc-ldpc码编、解码方法与系统 Expired - Fee Related CN108768593B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810604042.4A CN108768593B (zh) 2018-06-12 2018-06-12 一种基于dmt调制的qc-ldpc码编、解码方法与系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810604042.4A CN108768593B (zh) 2018-06-12 2018-06-12 一种基于dmt调制的qc-ldpc码编、解码方法与系统

Publications (2)

Publication Number Publication Date
CN108768593A CN108768593A (zh) 2018-11-06
CN108768593B true CN108768593B (zh) 2020-07-07

Family

ID=64021143

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810604042.4A Expired - Fee Related CN108768593B (zh) 2018-06-12 2018-06-12 一种基于dmt调制的qc-ldpc码编、解码方法与系统

Country Status (1)

Country Link
CN (1) CN108768593B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112187402B (zh) * 2019-07-05 2024-05-17 北京京东振世信息技术有限公司 一种数据处理的方法、装置和存储介质
CN115396065B (zh) * 2022-10-26 2023-04-28 南京邮电大学 一种稀疏随机线性网络编码的低时延解码方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1715588A1 (en) * 2005-04-22 2006-10-25 Broadcom Corporation Algebraic construction of LDPC (Low Density Parity Check) codes with corresponding parity check matrix having CSI (Cyclic Shifted Identity) sub-matrices
CN1941635A (zh) * 2005-09-30 2007-04-04 电子科技大学 一种非正则低密度奇偶校验码的构造方法
CN101159435A (zh) * 2007-11-14 2008-04-09 中国人民解放军理工大学 基于移位矩阵分级扩展的低密度校验码校验矩阵构造方法
US7617441B2 (en) * 2005-07-18 2009-11-10 Broadcom Corporation Efficient construction of LDPC (Low Density Parity Check) codes with corresponding parity check matrix having CSI (Cyclic Shifted Identity) sub-matrices
CN103152056A (zh) * 2013-01-30 2013-06-12 北京大学 一种基于原模图的准循环ldpc码构造方法及装置
CN103957015A (zh) * 2014-05-12 2014-07-30 福州大学 用于ldpc码解码的非均匀量化编码方法及其在解码器的应用
CN103973314A (zh) * 2013-01-24 2014-08-06 中国科学院声学研究所 一种基于ldpc的信号编解码方法、及接收端和发送端
CN106571832A (zh) * 2016-11-04 2017-04-19 华南理工大学 一种多进制ldpc码级联神经网络解码方法及装置
CN107612558A (zh) * 2017-10-18 2018-01-19 重庆邮电大学 一种基于斐波那契‑卢卡斯序列的大围长qc‑ldpc码构造方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7162684B2 (en) * 2003-01-27 2007-01-09 Texas Instruments Incorporated Efficient encoder for low-density-parity-check codes

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1715588A1 (en) * 2005-04-22 2006-10-25 Broadcom Corporation Algebraic construction of LDPC (Low Density Parity Check) codes with corresponding parity check matrix having CSI (Cyclic Shifted Identity) sub-matrices
US7617441B2 (en) * 2005-07-18 2009-11-10 Broadcom Corporation Efficient construction of LDPC (Low Density Parity Check) codes with corresponding parity check matrix having CSI (Cyclic Shifted Identity) sub-matrices
CN1941635A (zh) * 2005-09-30 2007-04-04 电子科技大学 一种非正则低密度奇偶校验码的构造方法
CN101159435A (zh) * 2007-11-14 2008-04-09 中国人民解放军理工大学 基于移位矩阵分级扩展的低密度校验码校验矩阵构造方法
CN103973314A (zh) * 2013-01-24 2014-08-06 中国科学院声学研究所 一种基于ldpc的信号编解码方法、及接收端和发送端
CN103152056A (zh) * 2013-01-30 2013-06-12 北京大学 一种基于原模图的准循环ldpc码构造方法及装置
CN103957015A (zh) * 2014-05-12 2014-07-30 福州大学 用于ldpc码解码的非均匀量化编码方法及其在解码器的应用
CN106571832A (zh) * 2016-11-04 2017-04-19 华南理工大学 一种多进制ldpc码级联神经网络解码方法及装置
CN107612558A (zh) * 2017-10-18 2018-01-19 重庆邮电大学 一种基于斐波那契‑卢卡斯序列的大围长qc‑ldpc码构造方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Researchon LDPC Coding for Free Space Optical Communication with Modified Bessel-Gaussian BeamsBased on K-distribution;Jin Wanga,FujiangZeng,Qingqing Liang and Xiaofeng Zhang;《MATEC Web of Conferences 2018》;20180131;全文 *
光通信系统中LDPC码构造方法的研究进展;袁建国，汪哲，孙雪敏，郭乔，胡潇月，吴英冬;《激光杂志》;20170531;第2017年第38卷(第11期);全文 *
基于改进2-DGRS码的QC-LDPC码高效构造;赵明; 张晓林;《通信学报》;20150225;第36卷(第2期);全文 *

Also Published As

Publication number Publication date
CN108768593A (zh) 2018-11-06

Similar Documents

Publication Publication Date Title
Eleftheriou et al. Low-density parity-check codes for digital subscriber lines
KR101103472B1 (ko) 무선 통신 시스템들에 대한 멀티레벨 쉐이핑을 위한 시스템, 방법 및 컴퓨터-판독가능 매체
US20190199573A1 (en) Trellis Shaping on LDPC-Coded Modulation DMT Systems
CN108768593B (zh) 一种基于dmt调制的qc-ldpc码编、解码方法与系统
CN108696470B (zh) 一种可见光通信的dmt调制方法
Jiang et al. Direct current bias optimization of the LDPC coded DCO-OFDM systems
CN111092663B (zh) 一种基于比特加权分布的光正交频分复用系统和通信方法
Tikka et al. Error Correction Using Hybrid Coding Algorithm for BER Reduction in MIMO-OFDM
Chen et al. Design and performance of the polar coded modulation for high mobility communications
Zhang et al. Complementary M‐ary orthogonal spreading OFDM architecture for HF communication link
Aljubouri et al. Performance Evaluation of Coded OFDM Using 4QPSK and 16QAM
El-Gohary et al. Utilization of raptor codes for OFDM-system performance enhancing
Alkhasraji et al. Performance investigation of coded IM/DD optical OFDM for short range underwater wireless channel
Selvi et al. Review of performance of ldpc codes for various ofdm systems
Arafa et al. Performance of combined constellation shaping and bit interleaved coded modulation with iterative decoding (bicm-id)
Band et al. BER performance of WIMAX system using wavelet packet modulation technique
CN116094893B (zh) 一种基于无码率码的海洋浮标ofdm机会通信方法
Yildiz et al. Comparison of communication performance of block coding methods in different channels using DWT
Lian et al. A Research on Modulation and Detection in OFDM-IM System
Zhang et al. Performance comparison of WOFDM with different coding schemes
Abd Al-Kareem LOW DENSITY PARITY CHECK (LDPC) CODES FOR A PROPOSED SLANTLET TRANSFORM OFDM SYSTEM IN A RAYLEIGH FADING CHANNELS WITH PERFECT AND PILOT CHANNEL ESTIMATION FOR M_ARY PSK MODULATION
Joshi et al. Performance Analysis and PAPR Reduction of Coded OFDM (with RSCC and Turbo coding) System using Modified SLM, PTS and DHT Precoding in Fading Environments
CN101783893B (zh) 一种多媒体广播无线信号抗噪声传输方法
Chandio et al. Peak to average power ratio (PAPR) reduction technique in orthogonal frequency division multiplexing (OFDM) using block coding
Song et al. Coded modulation with APSK constellations for power line communication

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200707

CF01 Termination of patent right due to non-payment of annual fee