-
在水下可见光通信系统中,光传输的两大主要影响因素分别为吸收和散射[2]。图 1是一个边长为Δr、体积为ΔV的水立方体。当波长为λ、功率为Pt(λ)的入射光束通过水立方体时,功率Pa(λ)被水立方吸收,功率Ps(λ)被水立方体散射,剩余的功率Pr(λ)不受任何影响地穿过水立方[19]。
由能量守恒定律可得:
$ P_{{\mathrm{t}}}(\lambda)=P_{{\mathrm{a}}}(\lambda)+P_{{\mathrm{s}}}(\lambda)+P_{{\mathrm{r}}}(\lambda) $
(1) 令吸收度为A,散射度为B, 分别定义为:
$ \left\{\begin{array}{l} A(\lambda)=\frac{P_{{\mathrm{a}}}(\lambda)}{P_{{\mathrm{t}}}(\lambda)} \\ B(\lambda)=\frac{P_{{\mathrm{s}}}(\lambda)}{P_{{\mathrm{t}}}(\lambda)} \end{array}\right. $
(2) 当水立方体的边长Δr趋近无穷小时,吸收系数和散射系数分别为:
$ \left\{\begin{array}{l} a(\lambda)=\lim\limits _{\Delta r \rightarrow 0} \frac{\Delta A(\lambda)}{\Delta r}=\frac{{\mathrm{d}} A(\lambda)}{{\mathrm{d}} r} \\ b(\lambda)=\lim\limits _{\Delta r \rightarrow 0} \frac{\Delta B(\lambda)}{\Delta r}=\frac{{\mathrm{d}} B(\lambda)}{{\mathrm{d}} r} \end{array}\right. $
(3) 水下总衰减系数为:
$ c(\lambda)=a(\lambda)+b(\lambda) $
(4) 此时,水下可见光传输的损耗系数[3]可表示为:
$ L(\lambda, d)=\exp [-c(\lambda) d] $
(5) 式中,d表示为水下可视通信距离。
图 2为水下信道模型。联合考虑发射设备孔径at、接收设备孔径ar和光源发散角θ等因素时,水下可见光信道传输公式[20]可表示为:
$ \begin{aligned} P_{{\mathrm{r}}}(\lambda)=& P_{{\mathrm{t}}}(\lambda) \eta_{{\mathrm{t}}} \eta_{{\mathrm{r}}}\left[\frac{a_{{\mathrm{t}}}^{2}}{\left(d \tan \theta+a_{{\mathrm{r}}}\right)^{2}}\right] \times \\ & \exp \left[-c\left(\lambda, \rho_{1}, \rho_{{\mathrm{p}}}\right) d\right] \end{aligned} $
(6) 式中,ηt为发射设备传输效率,ηr为接收设备传输效率,ρl为叶绿素的质量浓度,ρp为悬浮粒子质量浓度。
-
图 3为本文中所采用的DCO-OFDM水下可见光通信系统模型。在系统发送端,首先将串行的原始数据进行正交振幅调制(orthogonal amplitude modulation, QAM),得到离散频域信号{Xk}k=1N/2-1,接着对调制后的信号进行串并(serial/parallel, S/P)变换。图中, D/A(digital/analog)表示数模转换,PD(photodetector)表示光电探测器,CP(cycle prelix)表示循环流程。
不同于其它通信介质下的OFDM系统,DCO-OFDM水下可见光通信系统的时域信号需要符合实数特性要求。为了满足这个条件,需对Xk进行Hermitian共轭对称变换[9]。即:
$ X_{k}=\left\{\begin{array}{l} X_{k}, (k=1, 2, \cdots, N / 2-1) \\ X_{N-k}{ }^{*}, (k=N / 2+1, \cdots, N-1) \\ 0, (k=0, N / 2) \end{array}\right. $
(7) 式中,上标*表示共轭对称变换。变换后的信号为X=[X0, X1, …, XN-1]T,再通过IFFT运算后得到N点离散时域信号x=[x0, x1, …, xN-1]T,即:
$ \begin{array}{c} x_{n}=\frac{1}{\sqrt{N}} \sum\limits_{n=0}^{N-1} X_{n} \exp \left(\frac{{\mathrm{j}} 2 {\rm{ \mathsf{ π}}} k n}{N}\right), \\ (n=0, 1, 2, \cdots, N-1, k=1, 2, \cdots, N / 2-1) \end{array} $
(8) 对并串变换后的时域信号x添加循环前缀,以减小码间干扰的影响; 再经过数模转换将离散时域信号转换为连续时域信号x(t); 最后,需要对时域信号添加一个合适的直流偏置xbias使其转换为非负实数信号xDCO(t),从而可以驱动LED将所要传输的电信号以光的形式发送出去, 即:
$ x_{{\mathrm{DCO}}}(t)=\left\{\begin{array}{l} x(t)+x_{\rm {bias }}, \left(x(t) \geqslant-x_{\rm {bias }}\right) \\ 0, \left(x(t)<-x_{\rm {bias }}\right) \end{array}\right. $
(9) 式中,t是时间; xbias=$\mu \sqrt{E\left[x(t)^{2}\right]}$,其中, μ为比例常数, E表示对x(t)的数学期望,即平均功率。
经水下信道的传输到达接收端后,依次进行光电信号转换、模数转换、去循环前缀以及串并转换后得到离散时域信号y。然后,对信号y进行快速傅里叶变换(fast Fourier transform, FFT),以获得离散频域信号Y=[Y0, Y1, …, YN-1]。最后经反Hermitian共轭对称、并串转换和QAM解调,还原出原始数据。
-
PAPR定义为(单位为dB):
$ R_{{\mathrm{PAPR}}}=10 \lg \frac{\max \left|x_{n}\right|^{2}}{E\left[\left|x_{n}\right|^{2}\right]} $
(10) 式中,max|xn|2为信号峰值功率,E[|xn|2]为信号平均功率。
为了更直观地描述系统PAPR的性能好坏,通常使用互补累计分布函数(complementary cumulative distribution function, CCDF)来表示其PAPR的概率分布情况[21]。CCDF定义如下:
$ F_{{\mathrm{CCDF}}}=\left\{\begin{array}{l} P, \left(R_{{\mathrm{PAPR}}}>R_{0}\right) \\ 1-P, \left(R_{{\mathrm{PAPR}}} \le R_{0}\right) \end{array}\right. $
(11) $ P=1-\left[1-\exp \left(-R_{0}\right)\right]^{N} $
(12) 式中,R0为预设的PAPR阈值,P为累积分布函数定义。
-
设α1, α2, …, αj为j个实数或者复数,pi(x)(i=0, 1, …, j-1)是i次多项式,且满足递推关系式:
$ \left\{\begin{array}{l} p_{0}(x)=1 \\ p_{1}(x)=\theta_{1}\left(x-\beta_{1}\right) p_{0}(x) \\ p_{i}(x)=\theta_{i}\left(x-\beta_{i}\right) p_{i-1}(x)-\gamma_{i} p_{i-2}(x), (i \geqslant 2) \end{array}\right. $
(13) 式中,β, θ, γ表示可以满足VLM矩阵表达式的任意常量值。θi≠0(i=1, 2, …, j-1),则由pi(x)构成的j阶方阵p称为VLM[22]:
$ {\boldsymbol{p}}=\left[\begin{array}{cccc} p_{0}\left(\alpha_{1}\right) & p_{0}\left(\alpha_{2}\right) & \cdots & p_{0}\left(\alpha_{j}\right) \\ p_{1}\left(\alpha_{1}\right) & p_{1}\left(\alpha_{2}\right) & \cdots & p_{1}\left(\alpha_{j}\right) \\ \ldots & \ldots & \cdots & \cdots \\ p_{j-1}\left(\alpha_{1}\right) & p_{j-1}\left(\alpha_{2}\right) & \cdots & p_{j-1}\left(\alpha_{j}\right) \end{array}\right] $
(14) 在VLM中取切比雪夫多项式为pi(x),即:
$ p_{i}(x)=\cos (i \arccos x) $
(15) 满足:
$ \left\{\begin{array}{l} p_{0}(x)=1 \\ p_{1}(x)=x \\ p_{i}(x)=2 x p_{i-1}(x)-p_{i-2}(x), (i \geqslant 2) \end{array}\right. $
(16) 根据以上的条件可以得到VLM[23],如下式所示:
$ \left\{\begin{array}{l} {\boldsymbol{p}}_{i, j}^{(1)}=\sqrt{\frac{2}{N+1}}\left\{\cos \left[\frac{(i-1)(j-1) \pi}{N-1}\right]\right\} \\ {\boldsymbol{p}}_{i, j}^{(2)}=\sqrt{\frac{2}{N+1}} \cos \left[\frac{(i-1)(j-1 / 2) \pi}{N}\right] \end{array}\right. $
(17) 式中,下标i, j表示矩阵的行列,上标(1)和(2)表示两种不同的矩阵。
-
图 4为VLM预编码与改进的自适应缩放联合PAPR抑制算法系统模型。图中,IVLM表示逆范德蒙类矩阵(inverse VLM)。
假设有效子载波个数为N/2-1,通过对原始数据进行QAM调制和串并变换的操作得到离散频域信号:Xm=[Xm, 1, Xm, 2, …, Xm, N/2-1]T。由(17)式生成一个(N/2-1)×(N/2-1)的VLM预编码方阵p:
$ {\boldsymbol{p}}=\left[\begin{array}{ccc} p_{1 \times 1} & \cdots & p_{1 \times(N / 2-1)} \\ \vdots & \ddots & \vdots \\ p_{(N / 2-1) \times 1} & \cdots & p_{(N / 2-1) \times(N / 2-1)} \end{array}\right] $
(18) 将生成的方阵p与离散频域信号相乘,即Xm′= pXm=[X1, X2, …, XN/2-1]T,便可得到VLM变换后的信号。
再根据(7)式的方法对Xm′进行Hermitian共轭对称变换得到X=[X0, X1, …, XN-1]T。对X进行IFFT运算和并串变换输出一个离散时域信号x=[x0, x1, …, xN-1],同时在x的最前面添加循环前缀。
由于LED的电压线性工作范围十分有限,因此需要对信号x进行自适应缩放,以此满足这个范围[24]。参考文献[24]中所提出的自适应缩放方法与预编码方法相结合不能很好地抑制PAPR,本文作者在此基础上做出了改进,使其更好地与预编码方法相结合,达到最佳的PAPR抑制效果。假设时域信号的幅值范围为[xmin, xmax],LED的电压范围为[ul, uh],自适应缩放因子S由下式得到:
$ S=\frac{x_{\max }-x_{\min }}{u_{{\mathrm{h}}}-u_{1}}, (0<S<1) $
(19) 则经过缩放后的信号为xs=Sx,并将其进行数模转换为模拟信号x(t)。
最后在[ul, uh]之间选择一个最优的直流偏置xbias添加到模拟信号x(t)上,得到xDCO(t)=x(t)+xbias。其中:
$ x_{\rm {bias }}=\frac{u_{{\mathrm{h}}}+u_{1}}{2} $
(20) 同时为了让发送信号在[ul, uh]范围之内,对其进行限幅操作,限幅后得到发送信号xDCO′(t),由下式所示:
$ x_{{\mathrm{DCO}}}^{\prime}(t)=\left\{\begin{array}{l} u_{{\mathrm{h}}}, \left(x_{{\mathrm{DCO}}}(t)>u_{{\mathrm{h}}}\right) \\ x_{{\mathrm{DCO}}}(t), \left(u_{1} \le x_{{\mathrm{DCO}}}(t) \le u_{{\mathrm{h}}}\right) \\ u_{1}, \left(x_{{\mathrm{DCO}}}(t)<u_{1}\right) \end{array}\right. $
(21) 将待发送的电信号xDCO′(t)转换为光信号,通过LED发送出去。经由水下信道的传输,到达信号接收端执行发射端的逆操作。
-
DCO-OFDM水下可见光通信系统的PAPR与离散频域信号Xm的非周期自相关函数ρ(k)有关,它可以表示为:
$ \begin{array}{l} \rho(k)=\sum\limits_{i=1}^{N / 2-1} X_{{\mathrm{m}}, i+k} X_{{\mathrm{m}}, i}{ }^{*}, \\ (k=0, 1, \cdots, N / 2-2) \end{array} $
(22) 其与PAPR的关系可以表示为:
$ R_{{\mathrm{PAPR}}} \leqslant 1+\frac{2}{N / 2-1} \sum\limits_{k=1}^{N / 2-2}|\rho(k)| $
(23) 从(22)式和(23)式可以看出,具有较低自相关的输入信号产生的PAPR较低。VLM是由正交多项式构成,具有较好的去相关性,经过VLM预编码的输入信号的自相关性会变小,从而可以很好地抑制DCO-OFDM水下可见光通信系统的PAPR。
-
从缩放因子的比值计算(19)式可以看出,其分母是固定不变的LED上下门限差值,分子是动态可变的离散时域信号峰峰值,其动态幅值范围为[xmin, xmax], (xmin, xmax∈[-1, 1])。输入信号x通过与缩放因子S相乘得到:
$ x_{{\mathrm{s}}}=S {\boldsymbol{x}}=\frac{x_{\max }-x_{\min }}{u_{{\mathrm{h}}}-u_{1}} {\boldsymbol{x}}, (0<S<1) $
(24) 因缩放因子的取值范围为0 < S < 1, 信号的幅度值变小。在缩放后的信号上添加直流偏置进行限幅操作,其信号的幅度值都被控制在LED的有限范围内。系统PAPR可表示为:
$ R_{{\mathrm{PAPR}}}=\left[\frac{\left(S {\boldsymbol{x}}+x_{\rm {bias }}\right)_{\max }}{E\left(S {\boldsymbol{x}}+x_{{\mathrm{bias}}}\right)}\right]^{2} $
(25) 自适应缩放因子与离散时域信号峰峰值有着密切关系,预编码首先降低了信号自相关性而影响缩放因子的大小,当自适应缩放因子趋近于0时,系统PAPR会达到最低为1。故而改进的自适应缩放因子可以有效地抑制系统PAPR。
-
通过MATLAB仿真软件对文中所提算法进行验证。仿真参数为10000帧,128个子载波,调制方式为16QAM的DCO-OFDM系统,LED的电压范围为[1, 4],水下可见光通信信道参数如表 1所示。
Table 1. Simulation parameters of underwater visible light communication channel
ηt ηr αt αr θ c λ ρl ρp 0.91 0.91 0.003m 0.003m 0.33mrad 0.151 530nm 5mg/m3 1mg/L 图 5中比较了不同条件下的PAPR性能。当FCCDF=10-3时,与不使用任何抑制算法的DCO-OFDM水下可见光通信系统相比较,只使用VLM预编码的系统PAPR降低了1dB;只使用改进的自适应缩放的系统PAPR降低了2.2dB;而使用所提VLM预编码和改进的自适应缩放联合PAPR抑制算法效果最好,降低了3.2dB。图中, S表示改进的自适应缩放。
在DCO-OFDM水下可见光通信系统中,DCT预编码是一种无失真降低系统PAPR的有效方法。图 6中对此进行了仿真对比,在CCDF为10-3时,VLM预编码比DCT预编码的系统PAPR抑制性能更好,下降了0.5dB,所提VLM预编码和改进的自适应缩放联合算法比DCT预编码和改进的自适应缩放联合算法的系统PAPR下降了0.4dB。
图 7和图 8中仿真了以上两类情况下的误比特率性能。从图 7中可以看出,仅使用VLM预编码的系统与原始系统的误比特率相同,说明VLM预编码不会对系统误比特率有影响。在误比特率为10-4时,单独采用改进的自适应缩放的系统相比于原始的系统高10dB。这是由于对自适应缩放后的信号进行了偏置直流的添加和限幅操作,极大地改变了信号幅度范围,从而导致误比特率增大。因此采用VLM预编码与改进的自适应缩放联合算法,在同一信噪比下,误比特率和单独采用改进的自适应缩放相同。由此可知,VLM预编码与改进的自适应缩放联合算法是牺牲了误比特率的性能,改进了PAPR,为了达到同样的误比特率性能,需要增加信噪比。
从图 8中可以看出,采用VLM预编码和DCT预编码的DCO-OFDM水下可见光通信系统的误比特率与原始系统误比特率相同,说明预编码方式是一种不改变系统性能的PAPR抑制方式。而两种预编码方式与改进的自适应缩放联合算法,在误比特率为10-4时,信噪比比原始系统的高10dB。为了改善误比特率过高,可以增加信噪比。
DCO-OFDM水下可见光通信系统的峰均比抑制算法
Peak-to-average ratio suppression algorithm for DCO-OFDM underwater visible light communication system
-
摘要: 为了解决直流偏置光正交频分复用(DCO-OFDM)水下可见光通信系统中存在的高峰均功率比(PAPR)问题,采用了一种基于范德蒙类矩阵(VLM)预编码与改进的自适应缩放联合PAPR抑制的算法。首先,该算法利用VLM对频域信号进行预编码处理,以此降低其自相关性,再通过比值计算方式改进的自适应缩放因子对时域信号进行缩放,使其更好地减小发光二极管带来的非线性失真,最后实现了抑制系统PAPR的效果。结果表明,在互补累计分布函数为10-3时,相比于原始系统,此联合算法的PAPR下降了3.2dB。由此可知,该研究对抑制DCO-OFDM水下可见光通信系统的PAPR是有帮助的。
-
关键词:
- 光通信 /
- 直流偏置光正交频分复用 /
- 峰均功率比 /
- 范德蒙类矩阵 /
- 改进的自适缩放
Abstract: In order to solve the peak-to-average power ratio (PAPR) problem in the direct current biased optical orthogonal frequency division multiplexing (DCO-OFDM) underwater visible light communication system, a joint PAPR suppression algorithm based on Vandermonde-like matrix (VLM) precoding and improved adaptive scaling was adopted. Firstly, the frequency domain signal was preceded by VLM to reduce the autocorrelation of the signal, and then the time domain signal was adaptively scaled to reduce the nonlinear distortion caused by light-emitting diode (LED). Finally, the PAPR suppression effect of the system was realized. The results show that when the complementary commulative distribution function was 10-3, the PAPR of the joint algorithm was reduced by 3.2dB compared with the original system. Therefore, this research is helpful to suppress the PAPR of DCO-OFDM underwater visible light communication system. -
Table 1. Simulation parameters of underwater visible light communication channel
ηt ηr αt αr θ c λ ρl ρp 0.91 0.91 0.003m 0.003m 0.33mrad 0.151 530nm 5mg/m3 1mg/L -
[1] QIN H R, XIE Zh B, WENG Zh H. Research on routing algorithm in underwater wireless optical communication sensor networks[J]. Laser & Optoelectronics Progress, 2021, 58(11): 1106009(in Chinese). [2] ZENG Z, FU S, ZHANG H, et al. A survey of underwater optical wireless communications[J]. Communications Surveys & Tutorials, 2017, 19(1): 204-238. [3] SPAGNOLO G S, COⅡELLA L, LECCESE F. Underwater optical wireless communication: Overview[J]. Sensor, 2020, 20(8): 2261. doi: 10.3390/s20082261 [4] AKHOUNDI F, SALEHI J A, TASHAKORI A. Cellular underwater wireless optical CDMA network: Performance analysis and implementation concepts[J]. IEEE Transactions on Communications, 2015, 63(3): 882-891. doi: 10.1109/TCOMM.2015.2400441 [5] TANG S, DONG Y, ZHANG X. Receiver design for underwater wireless optical communication link based on APD[C]//7th International Conference on Communications and Networking in China. Kunming, China: IEEE, 2012: 301-305. [6] MI X, DONG Y. Polarized digital pulse interval modulation for underwater wireless optical communications[C]//OCEANS 2016-Shanghai. Shanghai, China: IEEE, 2016: 1-4. [7] WANG Y, XIE S, XIE Z. FISTA-based PAPR reduction method for tone reservation's OFDM system[J]. IEEE Wireless Communications Letters, 2018, 7(3): 300-303. doi: 10.1109/LWC.2017.2773478 [8] NAKAMURA K, MIZUKOSHI I, HANAWA M. Optical wireless transmission of 405nm, 1.45Gbit/s optical IM/DD-OFDM signals through a 4.8m underwater channel[J]. Optics Express, 2015, 23(2): 1558-1566. doi: 10.1364/OE.23.001558 [9] ABDULKAFI A A, ALIAS M Y, HUSSEIN Y S. Performance analysis of DCO-OFDM in VLC system[C]//IEEE Malaysia International Conference on Communications. Kuching, Malaysia: IEEE, 2016: 163-168. [10] ZHANG Y T, ZHAO L, ZHANG F. Performance optimization of visible light OFDM communication system based on wavelet transform[J]. Laser Technology, 2020, 44(2): 261-265(in Chinese). [11] ZHANG J, ZHANG J. Research on partial linear companding transform for reducing the clipping distortion of visible light OFDM system[J]. Journal of Signal Processing, 2016, 32(1): 113-118(in Chinese). [12] MIRIYALA G, MANI V V. A new PAPR reduction technique in DCO-OFDM for visible light communication systems[J]. Optics Communications, 2020, 474: 126064. doi: 10.1016/j.optcom.2020.126064 [13] AZARNIA G, SHARIFI A A. Clipping-based PAPR reduction of optical OFDM signals using compressive sensing: Bayesian signal reconstruction approach[J]. Optical Fiber Technology, 2021, 64: 102527. doi: 10.1016/j.yofte.2021.102527 [14] ZHANG T, ZOU Y, SUN J, et al. Improved companding transform for PAPR reduction in ACO-OFDM-based VLC systems[J]. IEEE Communications Letters, 2018, 22(6): 1180-1193. doi: 10.1109/LCOMM.2018.2827940 [15] IDRIS A, SAPARI N, IDRIS M S, et al. Reduction of PAPR using block coding method and APSK modulation techniques for F-OFDM in 5G system[C]//TENCON 2018-2018 IEEE Region 10 Confe-rence. Jeju, Korea (South): IEEE, 2018: 2456-2460. [16] ZHANG J, ZHANG J, WANG D M. PAPR reduction by DCT-GCS method in visible light communication DCO-OFDM system[J]. Journal of Signal Processing, 2015, 31: 39-44(in Chinese). [17] XIAO Y, CHEN L, LI F, et al. PAPR reduction based on chaos combined with SLM technique in optical OFDM IM/DD system[J]. Optical Fiber Technology, 2015, 21(1): 81-86. [18] SHARIFI A A, EMAMI H. PAPR reduction of asymmetrically clipped optical OFDM signals: Optimizing PTS technique using improved flower pollination algorithm[J]. Optics Communications, 2020, 474: 126057. doi: 10.1016/j.optcom.2020.126057 [19] BAI J R. Research on OFDM system PAPR reduction technology and its applications[D]. Xi'an: Northwestern Polytechnical University, 2018: 82-83(in Chinese). [20] BAI J R, CAO C, YANG Y, et al. Peak-to-average power ratio reduction for DCO-OFDM underwater optical wireless communication system based on an interleaving technique[J]. Optical Engineering, 2018, 57(8): 086110. [21] BAI J R, YANG S. In UOWC systems: A combined PAPR reduction method by PTS approach based on improved particle swarm optimization[J]. Optik, 2021, 232: 166581. doi: 10.1016/j.ijleo.2021.166581 [22] XU Z, YOU Z. A fast algorithm for inversion of confluent Vandermonde-like matrices involving polynomials that satisfy a three-term recurrence relation[J]. Siam Journal on Matrix Analysis & Applications, 1998, 19(3): 797-806. [23] HASAN M M. VLM precoded SLM technique for PAPR reduction in OFDM systems[J]. Wireless Personal Communications, 2013, 73(3): 791-801. doi: 10.1007/s11277-013-1217-6 [24] ABDULKAFI A A, ALIAS M Y, HUSSEIN Y S. A novel approach for PAPR reduction in OFDM-based visible light communications[C]//2017 International Conference on Platform Technology and Service (PlatCon). Busan, Korea (South): IEEE, 2017: 1-4.