无线光副载波调制系统相位噪声估计补偿算法

摘要：无线光副载波调制与数字调制相比，是一种有效抑制大气湍流影响的调制技术。对于相移键控调制格式，接收机电解调器引入的相位噪声会引起符号相干检测错误判决，导致通信系统误码性能劣化。文章在Malaga大气湍流信道下，分别分析了受大气湍流和相位噪声影响的星座图特征。基于卡尔曼滤波和扩展卡尔曼滤波对无线光副载波QPSK系统相位噪声进行了估计，并对受相位噪声影响的QPSK星座相位点进行了补偿。仿真结果表明，经过两种算法对相位噪声估计补偿后的系统误码率性能均有所改善，扩展卡尔曼滤波算法相对于卡尔曼滤波对相位噪声状态估计偏差小，补偿后误码率低，算法收敛速度快，可以有效改善系统性能。

关键词：副载波调制;相位噪声;Malaga大气湍流;相位噪声估计;误码率

中图分类号：TN929.1        文献标识码：A

文章编号：1009-3044（2022）21-0032-04

开放科学（资源服务）标识码（OSID）：

1 引言

无线光通信具有速率高、带宽大、可自由传输等优点，连续波调制与数字调制相比，具有抑制大气湍流光强起伏影响良好的性能[1-2]。对于无线光副载波多进制相移键控（M-ary Phase Shift Keying，MPSK）调制系统，接收端解调模块内本地振荡器所产生的相位噪声会导致载波恢复时所需要的相位不准确，星座图内的相位调制点偏离其所处的判决区域，引起解调时符号误判，系统性能劣化，因此快速准确预测或估计相位噪声显得尤为重要[3]。为了保证相位信息的正确，在解调前需要先对载波中相位噪声的大小进行估计，然后对相位噪声影响下偏移的相位点进行补偿，使其回到原本所处的判决区域内。

2 国内外研究进展

对于星座图内采样点相位均匀分布的相移键控等相位调制格式，维特比（Viterbi）算法可以直接估计相移键控的载波相位[4]，但算法实现硬件复杂度高。Tao Z等人提出了一种基于预判决的载波相位估计算法，硬件复杂度低，易于实现，且性能与维特比算法相当[5]，但对于高阶M-QAM（[M≥16]）调制，该算法不能直接使用。2014年，Bilal S M等人[6]基于前向反馈载波估计算法，对激光器相位噪声引起的64-QAM星座图相位点旋转后的相位信息进行了恢复。曹国亮等人[7]在偏振复用16-QAM光通信系统中，采用扩展卡尔曼滤波对偏振态和载波相位进行了估计，并对激光器相位噪声引起的星座图中采样点混叠进行了补偿。2014年Pakala L等人[8]运用扩展卡尔曼滤波算法分别实现了对相位噪声和幅度噪声的估计与补偿。2016年Jain A等人[9]提出一种基于直接判决的载波相位估计算法，实现了相干光16-QAM通信系统中由光纤克尔效应产生的非线性相位噪声的估计和补偿。2016年，Xie X M[10]提出了一种改进的多基线相位展开算法，该算法将无迹卡尔曼滤波与增强的联合相位梯度估计器相结合，能够同时进行相位展开和相位噪声抑制。2017年，文献[11]采用无迹卡尔曼滤波，在时域中对高阶M-QAM和大线宽相干光OFDM系统的激光器相位噪声进行估计，补偿了相位噪声引起的星座图采样点偏移。2017年，Shu T等人[12]扩展卡尔曼滤波器与Viterbi-Viterbi相位估计器相结合用于自适应非线性相位噪声抑制的性能。2018年，代亮亮等人[13]以并行卡尔曼滤波对偏振复用16-QAM传输系统中包含的相位噪声进行了估计与补偿。对于光通信系统内相位噪声的研究，文献[13]基于高阶正交幅度调制，以服从Wiener分布的激光器相位噪声为主要研究目标，通过卡尔曼滤波等滤波算法估计激光器相位噪声，并对由相位噪声造成的星座图相位点旋转进行补偿。

目前在无线光通信系统中相位噪声的研究大多是针对光通信系统内的光器件相位噪声进行的。本文基于无线光副载波MPSK调制系统，对接收机电域解调模块锁相环路所产生的相位噪声展开研究。在Malaga大气湍流信道下分析了相位噪声对无线光副载波QPSK星座图的影响，最后采用卡尔曼滤波和扩展卡尔曼滤波对QPSK系统中的相位噪声进行了估计和补偿技术研究。

3 Malaga大气湍流信道

Malaga大气湍流模型是一种通用的信道模型，可以描述目前常见的大气湍流信道。基于电磁波在湍流大气中传播的物理模型，Malaga分布的湍流信道光强I的概率密度函数可表示为[14-15]：

[fI（I）=Ak=1βakIα+k2-1Kα-k2αβIγβ+Ω']          （1）

其中：

[A=2αα2γ1+α2Γ（α）γβγβ+Ω′β+α2]

[ak=β-1k-11（k-1）！2Ω'γk-1αβk2γβ+Ω'1-k2]

[Ω′=Ω+2ρb0+22b0Ωρcos（φA-φB）]

上式中，[α]为正实参数，与散射过程中大尺度涡旋的有效数目有关;[β]为自然数，表示衰落参数的数量;[Kv·]表示第二类[n]阶修正贝塞尔函数;参数[Ω′]为相干平均光功率;[Ω]表示视距衰落分量的平均功率;[φA]和分[φB]别为视距传播分量[UL]与耦合到视距分量[UCS]的相位。[Γ·]表示Gamma函数。

4 相位噪声估计及补偿

4.1 相位噪声模型

无线光副载波通信系统接收端进行载波恢复时，电解调模块相位噪声主要由锁相环（Phase locked loop，PLL）内的本地振荡器输出的随机相位抖动与接收端接收信号相位间的相位差产生，该相位噪声[φ]服从Tikhonov分布，其概率密度函数为[16]：

[p（φ）=exp（αcos（φ））2πI0（α），  φ≤π]                         （2）

式中，[α]为环路信噪比[17]，且[α=A2/N0BL];[A2]表示接收信号功率;[BL]表示一阶锁相环的环路带宽，反映了环路对噪声的抑制作用;[N0]为加性高斯白噪声的单边功率谱密度。通常，[σφ2=N0BL/A2=1/α]，因此，式（2）可以写为：

[p（φ）=exp（cos（φ）/σ2φ）2πI0（1/σ2φ），  φ≤π]                        （3）

式中，[σ2φ]表示相位噪声方差，[I0（·）]表示零阶的第一类修正贝塞尔函数。

4.2相位噪声的卡尔曼滤波估计与补偿

假设经过光电探测器后，进入解调模块前的副载波MPSK电信号可通过式（4）表示：

[yk=akexp（jφk）+nk]                             （4）

式中[k]表示采样到的第[k]个MPSK符号，[ak]表示第[k]个传输符号调制信号，[φk]表示相位噪声，[nk]表示零均值加性高斯白噪声。[φk]服从Tikhonov分布。

卡尔曼滤波的状态方程和观测方程表示为：

[φk=Φφk-1+Γωk]                               （5）

[rk=Hφk+mk]                                        （6）

其中，[Φ]为状态转移矩阵，[Γ]为噪声驱动矩阵，[H]为观测矩阵，[ωk]和[mk]分别表示卡尔曼滤波中的过程噪声和观测噪声，均为零均值加性高斯白噪声。[Q]和[R]分别为[ωk]和[mk]的协方差矩阵。图1所示为卡尔曼滤波估计算法实现框图。

为求[φk]的估计值，卡尔曼滤波算法分为四个步骤：

1） 初始化

[φ00=0]，[P0=Varφ]                          （7）

其中[P0]为误差协方差阵，[Var[φ]]表示相位噪声[φ]的方差。

2）对第k个采样点的状态估计值和协方差进行预测

[φkk-1=Φφk-1k-1]，

[Pkk-1=ΓPk-1k-1Γ'+Q]                   （8）

其中[Q]为过程噪声的协方差阵。

3）求卡尔曼滤波增益

[Kk=Pkk-1H'HPkk-1H'+R-1]                       （9）

其中[R]为观测噪声的协方差矩阵。

4）对第k个采样点的状态估计值和协方差进行更新

[φkk=φkk-1+Kkrk-Hφkk-1]，

[Pkk=1-HKkPkk-1]                                    （10）

通过卡尔曼滤波求得对应状态的最佳估计值的一个计算周期，然后重复该步骤，就可以在误差最小的情况下递归地获取状态值的估计值。补偿后的MPSK信号表示为[14]：

[yk=yke-jφk]                                 （11）

式中，[φk]表示对应状态值的最佳估计值，[yk]表示进入解调模块前未经过相位噪声补偿算法补偿的副载波MPSK信号。

4.3 相位噪声的扩展卡尔曼滤波估计与补偿

扩展卡尔曼滤波的状态方程和观测方程表示为：

[φk=φk-1+ωk]                               （12）

[rk=akexpjφk+mk]                            （13）

其中，[ωk]和[mk]分别表示过程噪声和观测噪声，均为零均值加性高斯白噪声，协方差阵分别为[Q]、[R]。为求得相位噪声[φk]的估计值，所用扩展卡尔曼滤波算法可分为以下六个步骤：