本发明涉及光通信领域，更具体地说，涉及一种可见光通信的DMT调制方法与系统。
背景技术：
可见光通信(Visible Light Communication，VLC)作为一项近几年兴起的新型无线接入技术，以其功耗低、频带宽、保密性强和无电磁污染等优点，得到了世界范围内的广泛关注。在传统的VLC系统中，如图1所示，首先发射机将电信号转换为光信号通过LED1-4发送出去，随后接收端光电二极管(photodiode,PD)将接收光功率转换为电流。在VLC系统中，LED可以同时进行数据传输以及照明。
在可见光通信系统中，DMT调制可以有效地利用调制带宽实现较高的频谱利用率，同时对由多径传播引起的ISI具有适应性。基于快速傅里叶变换(fast Fourier transform,FFT)算法，可以通过DMT来实现多载波调制解调技术。本发明旨在FFT的基础上利用QC-LDPC码对VLC系统中传输的信息进行编码，目的在于降低VLC系统的ISI，进而提高系统的频谱利用率和数据传输速率，同时提高系统传输信息的可靠性以及精确性。
技术实现要素：
本发明要解决的技术问题在于，提出在FFT的基础上利用QC-LDPC码对VLC系统中传输的信息进行编码，从而形成的可见光通信的DMT调制方法与系统的技术方案，目的在于降低VLC系统的ISI，进而提高系统的频谱利用率和数据传输速率，同时提高系统传输信息的可靠性以及精确性。
根据本发明的第一方面，本发明解决其技术问题所采用的技术方案是：构造一种可见光通信的DMT调制方法，包含如下步骤：
(1)、将发送信息进行非正则LDPC编码后形成的串行数据流进行串并转换形成多组并行的二进制数据流；
(2)、每组二进制数据流分别经过QAM调制映射到星座图上；
(3)、将映射后形成的信号进行IFFT调制；
(4)、对经过IFFT调制的信号中每2Q个IFFT输出点形成的一个分组添加一个循环前缀；
(5)、将添加循环前缀后的信号进行并串转换，形成串行的数据流；
其中，每组二进制数据流经过QAM调制映射到星座图上所形成的信号作为IFFT调制的输入信号的前面一半信号，后面一半信号满足：
C2Q-q-1＝Cq*,
式中，q＝0、1、…、Q-1，C0至CQ-1为所述前面一半信号，*表示共轭，CQ至C2Q-1为所述后面一半信号，Q为步骤(1)中并行的二进制数据流的组数。
进一步地，在本发明的可见光通信的DMT调制方法中，C0以及CQ被设置为0。
进一步地，在本发明的可见光通信的DMT调制方法中，步骤(5)中形成的串行的数据流用于经过数模转换后，转换为模拟信号，然后通过可见光发送出去。
进一步地，在本发明的可见光通信的DMT调制方法中，所述LDPC编码在进行编码过程中，将信源产生的随机信息比特与生成矩阵相乘得到编码后的码字；其中，所述生成矩阵对应的校验矩阵为：
式中，H表示校验矩阵，I是大小为L×L的单位矩阵，P是经过循环移位之后的单位矩阵，P的上标表示循环移位的次数，该校验矩阵H的基矩阵的大小为m×n，并且L是质数满足L>n>m。
进一步地，在本发明的可见光通信的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)：
S22、变量节点信息处理：
对所有的变量节点i和与其相邻的校验节点j∈C(i)在第l次迭代时，根据公式(3)计算出校验节点传递给变量节点的消息L(l)(qij)：
S23、译码判决：根据公式(4)计算得出所有变量节点的硬判决信息L(l)(qi)：
然后根据硬判决信息L(l)(qi)得出译码结果若L(l)(qi)>0，则否则
S24、迭代停止判断：若满足或者已经达到预先设定的最大迭代次数，那么就停止译码并以最后一次译码的译码结果作为最终的译码结果进行输出，否则返回步骤S21继续迭代；式中，为本次的译码结果的转置；
其中，所述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调制方法中，步骤S24中，译码结果等于[c1c2…cN]。
进一步地，在本发明的可见光通信的DMT调制方法中，步骤(3)中经过IFFT调制后，每组2Q个IFFT点对应的输出端的信号为实值信号，这些实值信号为：
其中，U(g)为实值信号，g＝0,1,…2Q-1。
进一步地，在本发明的可见光通信的DMT调制方法中，步骤(4)中，所述实值信号经过添加一个循环前缀后的信号为：
根据本发明的另一方面，本发明为解决其技术，还提供了一种可见光通信的DMT调制系统，采用上述任一项的可见光通信的DMT调制方法对可见光通信进行DMT调制。
实施本发明的见光通信的DMT调制方法与系统，具有以下有益效果：降低了VLC系统的ISI，进而提高了系统的频谱利用率和数据传输速率，同时提高了系统传输信息的可靠性以及精确性。进一步地利用本发明的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 ratio belief propagation,LLR-BP)译码算法对输出端的信号进行译码。并且分析两种不同信道下QC-LDPC编码的DMT调制系统的性能。
VLC系统架构介绍如下。
A.信道分析
一个传统的光无线系统比如常见的基于LED的室内VLC系统(如图1所示)，可以用下式表示
其中y(t)和x(t)分别表示接收信号和发射信号的幅度，R表示接收端的光电转换系数，h(t)表示可见光信道的脉冲响应，n(t)表示加性高斯白噪声(additive white Gaussian noise,AWGN),假设该噪声是独立同分布的。
由于循环前缀(cyclic prefix,CP)的存在，DMT调制对于由多径效应引入的散射信号具有强鲁棒性。反射光功率可以忽略不计，只考虑直射链路。直射信道的光信道衰落根据LED的朗伯辐射模式可以用下式来模拟：
其中，朗伯模型的阶数m可以用下式表示m＝-ln2/ln(cosΦ1/2)，其中Φ1/2为LED的半功率角。A表示光电探测器的物理面积，d表示发射端与接收端的距离，φ表示发射光与发射光法线方向的夹角，ψ表示入射光与接收器平面之间的夹角，则在接收端，每个子载波在时域内的输出信号的幅度可以用下式表示：
B.码的构造
目前，可见光通信系统的研究主要集中在调制技术、信道均衡技术以及快速传输速率等方面。随着可见光通信关键技术和芯片的快速研发，未来的可见光通信将朝着超宽带，超高速，大覆盖范围和高可靠性的方向发展，包括移动大数据和高清影像等对误码性能要求高的多媒体业务。因此必须寻找一种低错误平层，高速率通信和性能优异的VLC系统信道编码方案。在该背景下，一些简单的分组码比如Reed-Solomon(RS)码或者具有较强纠错能力的LDPC码和turbo码就可以应用到VLC系统中。LDPC码是一类可以用稀疏校验矩阵H(Parity Check Matrix H)或二分图(Bipartite Graph)来描述的线性分组纠错码。LDPC码的最小汉明距离会随着码长的增加而线性增加。在进行后验概率迭代译码时，随着码字长度的增加，系统的BER会随之降低；而且采用迭代译码算法时LDPC码的性能非常接近香农极限。因此可以利用随机结构或者代数方法来构造LDPC码，但是当码长较短时，随机构造的LDPC码，其性能往往不是很理想。而且在这种情况下比较容易出现较高的误码平层。换句话说就是随机构造的码会使得系统在译码时比较困难。为了解决这个问题，本发明提出一种QC-LDPC码对系统的传输信息进行有效编码，其校验矩阵如下修正数组阵列所示：
这是一种非正则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)是一个满秩矩阵，因为它是一个上三角矩阵，对角线上的元素是非零元素。因此它的码率可以表示为
因为校验矩阵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＝P1’。
根据系统传输信息的表达式，置信传播(belief propagation,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)：
或
S22：变量节点信息处理
也可称作纵向步骤，对所有的变量节点i和与其相邻的校验节点j∈C(i)在第l次迭代时，根据式(10)计算出校验节点传递给变量节点的消息L(l)(qij)：
S23：译码判决
计算得出所有变量节点的硬判决信息为
若L(l)(qi)>0，则否则
S24：迭代停止判决
若满足或者已经达到预先设定的最大迭代次数，那么就停止译码并以最后一次译码的译码结果作为最终的译码结果进行输出，否则返回S21继续迭代。其中等于[c1c2…cN]。
C.DMT调制分析
如上所述，DMT调制是一种特殊的多载波调制技术，能够把一个高速串行数据流划分为多个并行低速数据流并且调制到不同频率的子载波上，而这些子载波之间要满足相互正交，这样就可以保证每个子载波在信道中传送数据时不会相互干扰，从而降低多径衰落的影响同时有效的抑制由多径效应引起的ISI问题，进而提高系统的频谱利用率和数据传输速率。
基于DMT调制的VLC系统其结构框图如图3所示。输入的高速串行数据流经过串并转换变成Q组并行低速二进制数据流。每组比特流经过正交幅度调制(quadrature amplitude 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点的输出端的信号为实值信号，这些实值信号可以表示为：
其中g＝0,1,…2Q-1。
在每2Q个IFFT点的分组中加入长度为P的循环前缀用来降低ISI的影响。在实际情况中，循环前缀一般是复制信号U(g)的后面一部分然后插入到信号U(g)的前端。那么此时信号可以表示为
其中g＝0,1,…2Q+P-1。
在IFFT模块输出端通过并串转换可以得到每个长度为(Q+P)的单位序列所对应的DMT时域信号的离散时间采样值。随后，输出的串行数据经过数模转换传输到信道中。在接收过程中经过光敏二极管，放大器，滤波器后的信号可以通过模数转换恢复到二进制数据。假设不考虑模数转换之前数模转换器的脉冲响应以及滤波器响应，接收端接收到的DMT帧信号(模数转换之前)可以表示为：
经过串并转换和删除循环前缀后的信号传输到FFT模块。理想情况下，此时FFT模块输出端的信号应该是从发射端发送到IFFT模块的原始信号。另外，利用信道的均衡处理达到消除或减小由发射机和接收机中有源或无源设备的参数变化引起的码间干扰的目的。最终信号进行解映射，接着通过并串转换，信号完成在DMT-VLC系统中的传输。
实验结果
在这一部分，主要研究循环前缀的长度，迭代次数以及码率对DMT-VLC系统性能的影响。本发明利用在不同信噪比情况下接收信号的误码率作为接收信号质量的判决标准。同时，利用LLR-BP译码算法，研究QC-LDPC码在两种不同信道下的性能。
根据在上述构造的QC-LDPC码，其中构造出了大小为Lm×Ln的校验矩阵，其中L＝31，m＝50，n＝70。接收端的信号可以表示为：
1.当发射端的信号为0时，接收信号的条件概率密度函数可以表示为：
2.当发射端的信号为1时，接收信号的条件概率密度函数可以表示为：
其中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。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。