-
水下大功率LED光通信系统中的光信号为单极性正脉冲,码字序列中会包含直流分量,在传输过程中,由于被调制的信号中存在不均匀的连"0"或连"1"伪随机码,导致直流分量的幅值产生一定的变化。当码字序列经过具有高通特性的交流耦合网络时,其中包含的直流分量会被滤除,导致输出脉冲出现基线浮动。为了更好地分析大功率LED调制信号特性,需要建立基线漂移的理论模型。调制信号可表示为[15]:
$ \boldsymbol{x}(t)=\sum\limits_{q} \boldsymbol{A}_{q} \boldsymbol{p}(t-q T) $
(1) 式中,Aq表示输入的信息序列数据(q为输入信息序列数据的总个数),t为时间,T是发送数据的间隔,p(t)是发射系统的脉冲信号。
本文中使用的系统接收端交流耦合网络为一个1阶高通滤波器,其频率响应函数Z(jω)可表示为:
$ Z(\mathrm{j} \omega)=\frac{\mathrm{j} \omega}{\mathrm{j} \omega+1 / \tau} $
(2) 式中,j为虚数单位,ω表示频率,τ为时间常数,决定了高通滤波器的截止频率。
计算得到1阶高通滤波器的脉冲信号v(t)为:
$ \boldsymbol{v}(t)=\left\{\begin{array}{l} \exp \left(\frac{-t}{\tau}\right),(0<t<T) \\ \exp \left(\frac{-t}{\tau}\right)\left[1-\exp \left(\frac{-T}{\tau}\right)\right],(t \geqslant T) \end{array}\right. $
(3) 将(3)式代入(1)式,得到发送信号经过交流耦合网络之后的信号vs(t)为:
$ \boldsymbol{v}_{\mathrm{s}}(t)=\sum\limits_{q} \boldsymbol{A}_{q} \boldsymbol{v}(t-q T) $
(4) 计算采样后的发送序列经过交流耦合网络之后的信号vs(uT)为:
$ \begin{gathered} \boldsymbol{v}_{\mathrm{s}}(u T)=\boldsymbol{a}_{u}+ \left\{\boldsymbol{b}_{u-1}-[1-\exp (-T / \tau)] \times\right. \\ \left.\left(\boldsymbol{a}_{u-1}+\boldsymbol{b}_{u-1}\right)\right\} \end{gathered} $
(5) 式中,信号由两部分组成,第一部分表示希望接收到的信号au,第二部分表示经过交流耦合网络后所带来的基线漂移部分,其中bu-1表示上一时刻的基线漂移量,u表示某一时刻的采样点。根据(5)式可以利用MATLAB对接收端交流耦合后的时域信号进行仿真。
图 1是接收端交流耦合后的时域信号仿真图。其中实线部分表示原始发送信号,虚线、点线、点划线部分分别表示经过不同截止频率下的交流耦合网络后的接收信号。通过分析可知,交流耦合网络的截止频率越高,接收信号的基线漂移现象越严重。
-
由于光束扩散和水下粒子等因素对光信号的吸收、散射作用[16],对光信号的传输有很大影响,以下是水下无线光信道传输模型[17]:
$ P_{\mathrm{r}}=P_{\mathrm{t}} \exp (-c d)\left(\frac{a_{\mathrm{r}}}{a_{\mathrm{t}}+d \theta}\right)^{2}+z $
(6) 式中,d为传输距离;at, ar分别为发射机的发射孔径半径和接收天线的孔径半径;θ为光束发散半角;Pr和Pt分别为接收信号光功率和发射信号光功率;此处的噪声z近似为高斯白噪声;c(λ)为衰减系数,由吸收系数a(λ)和散射系数b(λ)组成。表 1中给出了不同的海水类型中衰减系数的相应取值[18]。
Table 1. Typical attenuation coefficients for different water types
water types a(λ)/m-1 b(λ)/m-1 c(λ)/m-1 pure water 0.053 0.003 0.056 clear water 0.114 0.037 0.151 costal water 0.179 0.219 0.398 作者所在课题组通过水下信道外场实验发现,在清澈水质下测得的接收光功率符合水下信道传输模型[19](见(6)式),因此, 本文中参量将根据课题组外场实验而进行设置。
-
在实际应用中,随着码字长度的增加,译码出现错误的可能性会接近于0,但相应的复杂性和计算量也会增加。而水下通信需要快速的处理数据,不需要过于复杂的算法。为了在复杂度和性能之间取得平衡,本文中设计了级联交织码方案,采用LDPC码作为内码,RS码作为外码,并加入改进的块交织技术,提高码字方案处理错误的能力。
本文中设计的方案如图 2所示。图中,m为原始信息序列,m′为RS编码后的校验位比特序列;N为进行级联交织后的交织序列,同时N也作为LDPC编码前的信息序列,N′为LDPC编码后的校验位比特序列。
-
设计级联交织码方案时,首先进行RS编码。为了使内码的信息长度与外码相匹配,选择使用缩短RS码。RS(n, k)是一种纠正集中错误能力较强的线性分组码,其中n为总码长,k为信息位长度,进行RS编码首先定义本原元为α,生成多项式为:
$ \begin{gathered} \boldsymbol{g}(\boldsymbol{x})=(\boldsymbol{x}-\alpha)\left(\boldsymbol{x}-\alpha^{2}\right) \cdots\left(\boldsymbol{x}-\alpha^{2 e}\right)= \\ \prod\limits_{h=1}^{2 e}\left(\boldsymbol{x}-\alpha^{h}\right) \end{gathered} $
(7) 式中,e为RS编解码可纠错的个数,h的取值范围为1, 2, …, 2e。
信息多项式可以表示为(原始信息序列m=(m0, m1, …, mk-1)):
$ \boldsymbol{M}(\boldsymbol{x})=m_{0}+m_{1} \boldsymbol{x}+\cdots+m_{k-2} \boldsymbol{x}^{k-2}+m_{k-1} \boldsymbol{x}^{k-1} $
(8) 校验多项式如下式所示:
$ \boldsymbol{r}(\boldsymbol{x})=\boldsymbol{M}(\boldsymbol{x}) \cdot \boldsymbol{x}^{n-k}(\bmod \boldsymbol{g}(\boldsymbol{x})) $
(9) 式中,r(x)为n-k次校验多项式; M(x)乘以xn-k表示在原始的信息码字之后添加n-k个监督位, 则编码后的码多项式为:
$ \boldsymbol{c}(\boldsymbol{x})=\boldsymbol{x}^{n-k} \boldsymbol{M}(\boldsymbol{x})+\boldsymbol{r}(\boldsymbol{x}) $
(10) 对于RS译码,可分为以下三部分:(1)利用接收到的码字R(x)/g(x),得到伴随多项式;(2)由(1)式结果进一步得到错误图样Ê(x);(3)根据R(x)-Ê(x)计算出最大可能发送的码字Ĉ(x),完成RS码纠错。
-
在级联码方案设计过程中,要对RS码编码后的码字进行随机化处理,而且实际的水下信道中可能存在着突发干扰,会导致系统产生连串的突发错误,所以在级联码编码方案中引入交织技术,增强码字的纠错能力。
交织是将原始消息序列按某种规则进行打乱后传输,按相应规则在接收端还原后再进行纠错的一种技术。设置交织器是为了使码字序列之间的相关性变小[20-21],对于二元系统,假定每个信息位取0或1的概率为1/2,定义长度为l的两个序列X和Y的相关系数为:
$ \gamma_{\boldsymbol{XY}}=\sum\limits_{L=0}^{l}\left(2 \boldsymbol{X}_{L}-1\right)\left(2 \boldsymbol{Y}_{L}-1\right) $
(11) 式中,L的取值范围为0, 1, …, l。由(11)式可知,对于码字序列来说,交织前后相应位置上的不动点个数、同一位置的元素是否取同一值都与序列之间的相关性有关,其中同一位置的元素是否取同一值对相关系数的影响是一致的,所以研究交织前后相应位置上的不动点个数。本文中采用块交织技术,将要传输的信息序列存储在块交织矩阵中,按特定读写规则来处理数据,不同规则对应的输出序列不同。
图 3a中展示出了一种读写规则下的A型块交织器,将原始信息序列按从左至右的顺序,以行序写到交织器矩阵中,之后按由上往下的顺序,以列序从交织器矩阵中读出。为进一步降低信息序列的相关性,本文中改进此块交织器的读写规则,按从左至右的行序写入,但在读出时首先读取次列数据,之后依次读取次列之后的数据,最后读取首列数据,将其定义为B型块交织器,经过A、B两种块交织器后输出的信息序列如图 3b所示。
Figure 3. A-type block interleaver and interleaver output information sequence under different read and write rules
可以看出,对比两种规则下的交织器输出序列与原始消息序列,A型交织器的输出序列在交织前后的不动点个数为5,改进后B型交织器的输出序列在交织前后的不动点个数为0,由此可见, 改进后的交织器可明显降低序列之间的相关性,因此本文中使用改进后的B型块交织器技术。
-
LDPC码是一种(n, k)线性分组码,其中n为总码长,k为信息位长度,LDPC码可由校验矩阵H来唯一确定。本文中采用随机法来生成H矩阵,使用的是规则LDPC码。在进行LDPC编码时,首先得到校验矩阵H和生成矩阵G,两者之间的关系如下所示:
$ \boldsymbol{H}=\left[\boldsymbol{I}_{(n-k) \times(n-k)} \mid-\left(\boldsymbol{P}_{k \times(n-k)}\right)^{\mathrm{T}}\right] $
(12) $ \boldsymbol{G}=\left[\boldsymbol{P}_{k \times(n-k)} \mid \boldsymbol{I}_{k \times k}\right] $
(13) 式中,P是通过高斯消元后得到的矩阵,I为单位矩阵。
根据生成矩阵和原始信息序列,即可得到编码后的码字c:
$ \boldsymbol{c}=\boldsymbol{m} \times \boldsymbol{G} $
(14) 然后进行LDPC码译码,考虑到级联码会存在门限效应,所以采用迭代译码。本文中使用了置信度传播(belief propagation,BP)算法,其中基于OOK调制的BP译码算法初始概率密度函数如下式所示[20]:
$ \begin{array}{c} P\left(i_{f}\right)= \\ \left\{\begin{array}{l} \frac{\exp (2 W-1)}{1+\exp \left(\frac{2 W-1}{2 \sigma^{2}}\right)}=\frac{1}{1+\exp \left(\frac{1-2 W}{2 \sigma^{2}}\right)},(f=1) \\ 1-P\left(i_{1}\right)=1-\frac{1}{1+\exp \left(\frac{1-2 W}{2 \sigma^{2}}\right)},(f=0) \end{array}\right. \end{array} $
(15) 式中,P(if)表示解码信号i解码为1或0的概率,W表示接收符号,σ2描述了噪声的强度。得到初始概率密度函数后,完成节点更新和最终判决,得到译码结果。
-
水下LED光通信系统框图如图 4所示。系统分为发射系统、海水信道、接收系统三部分,发射系统完成信号的产生、编码以及调制并进行传送,经过海水信道,在接收系统进行接收和信号处理。
本文中各项参量的设置如表 2所示。
Table 2. System parameter settings
name parameter name parameter LED wavelength 450nm receiver sensitivity <1μW transmit power 5W receiver dadius 100mm divergence half angle of light source 0.087rad detector size 10mm×10mm transmitter radius 50mm collecting efficiency of receiving antenna 76.46% -
为研究级联交织码方案中不同码字参量对系统误比特率性能的影响,通过MATLAB搭建水下LED光通信系统模型,考虑水下无线光信道传输中光信号的衰减以及基线漂移特性,对比系统BER性能最终确定出级联交织码方案的优化参量,其中水质衰减系数取0.056m-1(表示纯海水),通信距离选为10m,系统其它参量按实验系统参量设置。
在设计码字方案时,级联RS码会使级联交织码的整体码率降低,导致加在编码后码字上的噪声能量增加,当RS码的码率较低时,可能会使码字纠错性能下降;其次在级联码译码过程中,先由LDPC码处理错误,若存在LDPC码处理不了的错误码字,RS译码器可进一步对其进行改正,所以在设计级联交织码中的RS码字方案时,可以适当减少相应的校验位冗余。本文中选取RS(240,158), RS(240,208), RS(240,224)分别与LDPC(2880,1920)级联,交织深度为4,得到不同RS码率下级联交织码系统BER曲线,如图 5所示。由结果可知,RS(240,208)+LDPC(2880,1920)级联交织码方案下的水下LED光通信系统性能较好。由此可见,在单独使用RS码纠错时,码率越低纠错能力越强,但使用级联交织码方案时可适当增加RS码率,减少相应的校验位冗余。
RS码组确定后内码LDPC的信息位个数也随之确定,下面研究级联码中LDPC码的码长对系统性能的影响。选用RS(240,208)作为外码,交织深度为4,得到不同LDPC码码长下级联交织码的系统BER曲线,如图 6所示。
由对比结果可知RS(240,208)+LDPC(2880,1920)级联交织码方案下的系统性能较好。与使用RS(240,208)+LDPC(2560,1920), RS(240,208)+LDPC(2340,1920),RS(240,208)+LDPC(2136,1920)方案的系统相比分别具有的增益为1dB, 1.7dB, 2dB。因此在设计级联码方案时可以增加内码LDPC码的校验位个数。但需注意噪声会对码字校验位造成影响,使水下LED光通信系统纠错性能变差,同时码长较长时还会导致较大时延,所以要合理设置LDPC码长使系统达到最佳性能。
在RS(240,208)+LDPC(2880,1920)级联码方案中使用B型块交织技术,研究不同的交织深度对系统BER性能的影响。图 7是不同交织深度参量下的性能对比。从结果可以看出,加入交织技术后系统BER曲线收敛速度更快,且随着交织深度增大,系统性能越好。
通过分析,最终确定采用RS(240,208)+LDPC(2880,1920)级联码方案进行模拟实验。但在实际应用中随着交织深度增大,会导致数据处理的时间变长,因此本文中选取10作为交织深度。
-
用MATLAB产生随机信息序列并进行级联交织编码,下载到现场可编程逻辑门阵列(field-programmable gate array,FPGA)上,使用FPGA对发送序列添加一定信噪比的高斯白噪声,通过光发送系统实现光信号的OOK调制。使用衰减系数为0.151m-1的水质模拟清澈海水信道,光信号通过该种水质信道之后,在接收端进行数据采集,之后利用MATLAB将采集到的信号进行离线处理。离线处理包括信号的同步、译码等过程。实验中统计1s内传输数据中错误的信息个数,每组数据量为993600bits。
在接收端接收到的未编码信号波形如图 8所示。可观察到经过直流滤波后连"0"、连"1"信号引起的基线漂移现象(圆圈标识出的部分)以及水下信道存在的各类干扰对信号的衰减,若直接进行信号判决将会产生较大的误码。本文中将采集到的基线漂移信号在MATLAB上进行了离线处理,采用未编码、LDPC(2880,1920),RS(240,208),RS(240,208)+LDPC(2880,1920)等码字方案,固定通信距离为10m,得到的系统BER曲线如图 9所示。
由图 9可知,系统误比特率为3.8×10-3时,与未编码系统相比,级联交织码系统具有大约3.8dB的编码增益,可有效提升水下LED光通信系统的BER性能。与单独加入LDPC码和RS码的系统相比,级联交织码系统的增益分别为1.2dB和2dB。
通过固定系统发射端,移动接收端来控制通信距离,采用RS(240,208)+LDPC(2880,1920),RS(240,208),LDPC(2880,1920)码字方案,对发送序列添加的高斯白噪声信噪比为10dB,研究不同方案下的水下LED光通信系统的通信距离与系统BER性能的关系,如图 10所示。
由图 10可知,当通信距离为25m时,采用RS(240,208)+LDPC(2880,1920)方案的系统BER为1.04×10-4;采用LDPC(2880,1920)码字方案系统BER为7.01×10-4;采用RS(240,208)码字方案的系统BER为8.67×10-3。当通信距离增加到30m时,使用LDPC码和RS码的系统BER均高于3.8×10-3,已无法满足可靠通信要求,而LDPC-RS级联交织码系统BER仍低于3.8×10-3。对比可知, 当系统传输相同距离时,与LDPC码或RS码系统相比,级联交织码系统性能更优。
通过选取1Mbit/s,2Mbit/s,4Mbit/s,8Mbit/s几种不同的传输速率,研究了系统传输速率与系统误比特率之间的关系,信噪比为9dB,采用RS(240,208)+LDPC(2880,1920)码字方案,如图 11所示。
由图 11可知,当系统传输速率为1Mbit/s时,使用级联交织码方案的系统误比特率为2.41×10-5;当系统传输速率为2Mbit/s时,使用级联交织码方案的系统误比特率为2.59×10-5;当系统传输速率为4Mbit/s时,使用级联交织码方案的系统误比特率为2.82×10-5;当系统传输速率为8Mbit/s时,使用级联交织码方案的系统误比特率为3.03×10-5;对比可知, 在同样的信噪比环境下,随着系统传输速率的增加,系统误比特率在逐渐增大。这是因为当信噪比不变时,随着传输速率的增加,每比特对噪声的容忍度会急剧降低,出现误码的概率会增加,导致误比特率变差。另外在系统带宽一定的条件下,随着传输速率的增加,基线漂移现象加重,因此也会导致系统误比特率增加,但是通过纠错编码可提高误比特率性能。
基于级联交织码的水下LED光通信误码性能分析
Analysis of error code performance of underwater LED optical communication based on cascaded interleaving code
-
摘要: 为了提升系统误比特率,减小基线漂移以及海水信道的吸收散射等特性对光信号产生的影响,采用了基于水下发光二极管(LED)光通信系统的低密度奇偶校验码(LDPC)-里所(RS)级联交织码方案,在模拟水下LED光通信实验系统的情况下,分析码字方案中RS码、LDPC码以及交织参量对系统误比特率性能的影响,得到了级联交织码方案的优化参量,并进行了实验模拟验证。结果表明,优化后的级联交织码系统与未编码系统、RS码系统、LDPC码系统相比分别可获得3.8dB,2dB,1.2dB的增益,可有效提高系统的误比特率性能。该研究为提高水下无线光通信系统的可靠性提供了参考。Abstract: In order to improve the bit error rate of the system, reduce the baseline drift and the absorption scattering of the seawater channel and other characteristics on the optical signal, a low density parity check (LDPC)-Reed Solomon(RS) cascaded interleaving code scheme based on the underwater light-emitting diode(LED) optical communication system was adopted. Under the simulated underwater LED optical communication experimental system, the effects of RS code, LDPC code, and interleaved parameters on the bit-error rate performance of the system were analyzed in the code scheme, and the optimized parameters of the cascaded interleaving code scheme were obtained and verified by the experimental simulation. The results show that the gain of 3.8dB, 2dB, and 1.2dB can be obtained with the optimized cascaded interleaving code system, compared with the uncoded system, RS code system, and LDPC code system, respectively, which can effectively improve the bit error rate performance of the system. This study provides a reference for improving the reliability of the underwater wireless optical communication system.
-
Key words:
- optical communication /
- baseline drift /
- cascaded code /
- interleave
-
Table 1. Typical attenuation coefficients for different water types
water types a(λ)/m-1 b(λ)/m-1 c(λ)/m-1 pure water 0.053 0.003 0.056 clear water 0.114 0.037 0.151 costal water 0.179 0.219 0.398 Table 2. System parameter settings
name parameter name parameter LED wavelength 450nm receiver sensitivity <1μW transmit power 5W receiver dadius 100mm divergence half angle of light source 0.087rad detector size 10mm×10mm transmitter radius 50mm collecting efficiency of receiving antenna 76.46% -
[1] TAO Y Sh, ZHANG Ch. Research on underwater LED communication based on FBMC technology[J]. Electronic Measurement Technology, 2020, 43(10): 153-157(in Chinese). [2] CHI N, WANG Ch F, LI W P, et al. Research progress of underwater visible light cummunication technology based on blue/green LED[J]. Journal of Fudan University(Natural Science Edition), 2019, 58(5): 537-548(in Chinese). [3] CHI N, CHEN H. Progress and prospect of high-speed visible light communication[J]. Opto-Electronic Engineering, 2020, 47(3): 6-17(in Chinese). [4] QIN L, ZHANG Y Q, LI B S, et al. LED visible light communication systems based on MIMO technology[J]. Laser Technology, 2019, 43(4): 103-109(in Chinese). [5] LIN T, HUANG Zh T, JI Y F. Hybrid run length limited code and pre-emphasis technique to reduce wander and jitter on on-off keying nonreturn-to-zero visible light communication systems[J]. Optical Engineering, 2016, 55(11): 1105031. [6] ZENG Zh Q, SHU F, ZHANG H H, et al. A survey of underwater optical wireless communications[J]. IEEE Communications Surveys & Tutorials, 2017, 19(1): 204-238. [7] FORNEY D, GEORG D. Concatenated codes[M]. Cambridge, MA, USA: MIT Press, 1966: 113-124. [8] DONG Y, WANG J M, FANG X Y, et al. Simulation research of RS-Turbo code in indoor OFDM based visible light communication system[J]. Electronic Science & Technology, 2018, 31(10): 11-14(in Ch-inese). [9] MA Ch B, TANG Ch P, AO J, et al. Design of laser communication system with variable code length serial concatenated code[J]. Optical Communication Technology, 2019, 43(5): 40-43(in Chinese). [10] CAO Y, ZHANG X, PENG X F, et al. Cascade scheme based on multiple-output in spatial optical communication[J]. Acta Optical Sinica, 2018, 38(1): 0106003 (in Chinese). doi: 10.3788/AOS201838.0106003 [11] COSSU G, CORSINI R, KHALID A M, et al. Experimental demonstration of high speed underwater visible light communications[C]//2013 2nd International Workshop on Optical Wireless Communications. New York, USA: IEEE, 2014: 11-15. [12] WANG W P, ZHENG B. The simulation design of LED-based closerange underwater optical communication system[C]//International Computer Conference Wavelet Active Media Technology and Information Processing. New York, USA: IEEE, 2013: 283-285. [13] WANG J, ZENG F J, LIANG Q Q, et al. Research on LDPC coding for free space optical communication with modified Bessel-Gaussian beams based on K-distribution[C]// Matec Web of Conferences. Paris, France: EDP Science, 2018: 1-9. [14] HU S Q, MI L, ZHOU T H, et al. 35.88 attenuation lengths and 332 bits/photon underwater optical wireless communication based on photon-counting receiver with 256-PPM[J]. Optics Express, 2018, 26(17): 21685-21699. doi: 10.1364/OE.26.021685 [15] CHEN L, MIAO D. Study on the baseline excursion in optical fiber network of missile control system[J]. Optical Fiber & Electric Cable, 2006, 13(3): 90-105(in Chinese). [16] YANG Y, HE F T, GUO Q P, et al. Analysis of underwater wireless optical communication system performance. [J]. Applied Optics, 2019, 58(36): 9808-9814. doi: 10.1364/AO.58.009808 [17] YANG Y, ZHANG J N, YIN Y F. Performance analysis of LDPC code in seawater channel model[J]. Optical Communication Technology, 2018, 42(6): 29-32(in Chinese). [18] YE D M, LIU Y, WANG Y T, et al. Underwater LED optical communication technology of real-time error detection[J]. Infrared and Laser Engineering, 2019, 48(9): 151-157(in Chinese). [19] WANG H X, HE F T, ZHAN F, et al. Design of underwater LED blue light communication system based on PPM[J]. Optical Communication Technology, 2018, 42(7): 46-50(in Chinese). [20] LIU D H, TANG Ch J. Design of interleaver in Turbo code system[J]. Wireless Communication Technology, 2000, 9(1): 41-44(in Chinese). [21] JIA M J, LU A J. Design and simulation of a new digital interleaving[J]. Modern Information Technology, 2019, 3(5): 53-55(in Chinese).