基于互补斯托克斯向量实现偏振态旋转的追踪和预测 下载: 985次
1 引言
近年来,偏振复用相干光通信系统得到了越来越多的研究,并已产品化[1-2]。然而在一些极端情况下,链路中的高速偏振态旋转(RSOP)可能由于振动或大风等因素分别在地缆或天缆中产生。过高的偏振态旋转速度可能令相干接收机中的动态均衡模块追踪失效。追踪偏振态旋转的另一个问题在于动态均衡模块的反馈延迟。目前,高速相干系统大多采用频域恒模算法(CMA)或最小方均(LMS)算法等实现动态均衡[3-5],相比于时域均衡技术,在动态均衡系统中,串流信号以帧的形式进行并行处理,并利用快速傅里叶变换和逆变换实现卷积,以降低均衡的计算复杂度[6-7]。在实际系统中,帧长一般为256至512符号,而由于并行结构,动态均衡算法的反馈延迟大约需要几十帧。这些延迟给偏振态旋转的追踪带来了双重的困难,一方面,更大的反馈延迟导致追踪变得更困难,另一方面,延迟的存在使得系统实际需要的不是偏振态旋转的追踪,而是对应几十帧延迟的偏振预测。关于延迟带来的偏振态旋转追踪问题及偏振态旋转追踪的预测方案,目前鲜有报道。
目前,偏振态旋转追踪算法一般在琼斯空间对琼斯矩阵的参量变化进行追踪。算法包括恒模算法、为了弥补恒模算法中高阶调制信号噪声模不为零而改进的半径指导的均衡算法和多模算法、为了改进恒模算法中的奇异点问题而提出的改进算法等[8-10]。另外卡尔曼算法也是追踪琼斯矩阵变量的常用算法,包括扩展卡尔曼算法、线性卡尔曼算法、为了解决非线性问题的非线性卡尔曼算法等[11-14]。
本文提出了一种偏振态旋转预测算法,该算法将频域均衡输出得到的琼斯空间偏振旋转矩阵转换为斯托克斯空间向量,再进行循迹追踪和预测。相比关于琼斯空间的追踪方案,所提方案中,在斯托克斯空间的偏振态旋转轨迹更接近于线性空间,有助于通过简单算法实现高效的偏振态旋转追踪和预测。对偏振态旋转进行了理论分析,并在存在延迟的情况下,使用所提算法对偏振态旋转轨迹进行了预测,分析了所提算法的仿真效果。通过采用对互补向量进行循迹追踪的方法,预测的偏振态旋转结果误差向量幅度(EVM)可比原始方案低10%,误码率(BER)可从10-3量级降低到10-5量级。通过采用所提算法,系统的偏振态旋转速度容限从1.3 MHz提升至1.5 MHz。
2 琼斯空间旋转矩阵和斯托克斯空间向量的转换
为了对偏振态旋转进行分析,需要分别分析偏振态旋转在琼斯空间和斯托克斯空间的表达和相互转换的方法。偏振态旋转在琼斯空间中可以表示为一个琼斯矩阵,令U为琼斯空间中的偏振旋转矩阵,则U为酉矩阵,有U-1=UH,令V=UH,其中上角标H为厄米特转置,则U和V可写为
式中:
式中:I为二维单位矩阵;σ为泡利旋转矩阵。
其中,
式中:u和v的空间R为半径为π的三维球内部。通过(2)、(3)式不难推出,u和v方向相反,且u+v=π,因此称u和v为互补斯托克斯向量。
3 偏振旋转矩阵在斯托克斯空间的运动轨迹
RSOP在斯托克斯空间中可以表示为一个点的旋转,目前其他工作对RSOP转速的仿真都是恒定转速的。仿真模型中使用的是多段合成的随机旋转速度,在其中的每一段中,所提RSOP模型是随机方向的静态旋转模型,每段的转速服从高斯分布,方差为σ2。所以,为了分析在斯托克斯空间中RSOP的运动方式,需要先推导出RSOP的等效转速。球表面积S的表达式为
式中:r为球的半径;θ为旋转角度。球面上任一点绕任意轴旋转一周的长度均值为
即针对赤道上一点的旋转长度,系数为
即针对赤道上一点的旋转长度的平方均值,系数为
所以速度均值为
对于单位转速的等效平均速度,表达式为
从(6)~(10)式得出了在等效的斯托克斯空间中的运动速度和实际的RSOP转速之间的关系。结合(1)、(2)式,即可仿真出RSOP,
从
图 2. RSOP预测方案示意图。(a)圆拟合;(b)向量追踪
Fig. 2. Schematic of RSOP prediction scheme. (a) Curve fitting; (b) vector tracking
已知前三个点的位置S1,S2,S3,由于三个点可以确定一个圆,可以计算出这个圆的圆心O和半径,然后假设在较短的时间内旋转速度不变且RSOP点都是在同一个圆上运动,根据每次前进的角度θ1和θ2,就可以计算出延迟时间后点的前进角度和位置。然而圆拟合方案在每次预测时都要计算圆的位置,包括计算圆心、半径、前进角度,运算量比较大且只考虑了在一定时间内以固定的角速度旋转的情况,因此提出了另一种更加合理的仿真方法,即通过向量追踪RSOP的运动轨迹,如
为了测试圆拟合方案和向量追踪方案的效果,估计使用和不使用互补向量算法的差别,仿真了1000个连续RSOP矩阵,测试了延迟为5帧时的效果,使用琼斯矩阵参数的平方误差来衡量算法的效果。
图 3. 圆拟合方案和向量追踪方案的琼斯矩阵平方误差曲线
Fig. 3. Square error curves of Jones matrix obtained by circle fitting scheme and vector tracking scheme
u和v两个向量都是RSOP向量,分别使用向量追踪算法对两个互补向量进行了追踪。发现两个向量得到的平方误差曲线,一条误差较大时,另一条的误差则较小。因为u和v两个向量为互补向量,将u和v构成的空间R进一步分为半径为π/2的球内部和半径在π/2至π之间的空心球,分别称它们为内球和外球,则对于任意旋转矩阵,u和v必然一个在内球一个在外球。空间R为非线性空间,R的表面均对应U=-I,而其球心位置附近可近似为线性空间,因此,对于时变的U和对应的V,选择位于内球的一个矩阵进行RSOP的追踪和预测。由此提出了一种互补向量RSOP算法。
使用和不使用自动选择互补向量算法的平方误差曲线如
图 4. 互补向量追踪算法和向量追踪算法的平方误差曲线
Fig. 4. Square error curves of vector tracking algorithm and complementary vector tracking algorithm
4 算法原理与流程
对LMS算法和互补向量RSOP矩阵预测算法进行对比。LMS算法的流程如
从
LS算法表达式为
式中:Sf为经LMS算法后的频域信号;rf为经LMS算法之前的频域信号;M为二阶矩阵。
首先任意一个二阶矩阵可以表示为
式中:φ为相位偏转角;M2为形如(1)式的偏振旋转矩阵;厄米特矩阵M1可以表示为
显然,M1与其厄米特共轭矩阵相等,所以将M与其厄米特共轭矩阵相乘,可得
再将(14)式代入(12)式,即可得
即得到了用于RSOP预测的偏振旋转矩阵M2,用琼斯矩阵来进行RSOP的预测,另外因为琼斯矩阵的模为1,所以可以用模值计算出相位旋转角,用来进行载波相位恢复。
RSOP预测算法流程如
RSOP在斯托克斯空间中的位移可以分为两个方向来计算,分别是沿着当前速度方向的位移X0和旋转方向的位移X1,总位移为两个方向的位移之和。先计算单位延时后的RSOP在斯托克斯空间的位移,假设已知当前时间和上一时间的RSOP点的位置,则可以计算出当前时间k的速度
算下一时间的两个方向的位移,取d个时间的位移进行累加,即可计算得到d时间之后点的轨迹。为了简化编程,使用泰勒展开对计算过程进行化简,化简后的计算公式为
式中:T为RSOP在斯托克斯空间的轨迹点;vk为速度;av为加速度;ta为两个速度向量之间的旋转角度;n1,n2,n3及n4是为了简化积分计算,通过泰勒展开得到的参数,计算方法为
计算d个延迟后轨迹点的位置后,再对更新系数进行更新,并计算下一时间的速度、加速度、旋转方向、移动方向等。通过旋转方向和速度的预测,可以计算出偏振态旋转的琼斯矩阵,即可对偏振态旋转进行均衡。
5 仿真结果
绘制如
在不同的RSOP旋转速度下,对使用互补向量的RSOP矩阵预测算法和不使用RSOP只使用LMS的算法的BER和EVM进行对比。仿真调制格式为16QAM,Block长度为256,结果如
从
6 结论
首先分析RSOP在琼斯空间和斯托克斯空间之间的转换,提出一种将RSOP矩阵转换到斯托克斯空间中的算法,然后,在斯托克斯空间中使用多段合成的随机旋转速度对RSOP进行建模。依据此模型,提出了两种圆拟合和计算旋转点的向量速度加速度的预测方案,并对RSOP矩阵进行预测,最终经过对比,选择了计算向量运动的方式,并且为了提高算法精度,提出了互补向量算法。经过对比,使用互补向量算法的平方误差明显小于不使用互补向量的算法。从RSOP轨迹的追踪图中可以看出,所提算法可以较好地预测RSOP矩阵的运动轨迹,使解决实际通信系统中的延迟问题成为可能。最后,对比了互补向量预测算法和LMS算法的EVM和BER曲线。经过计算可得,互补向量的寻迹追踪算法的EVM可以降低10%,且BER可从10-3量级降低到10-5量级,偏振态旋转的速度容限从1.3 MHz提升至1.5 MHz。由此可知,所提算法可以较好地预测RSOP矩阵,通信效果有明显提升,具有良好的应用前景。
[2] Shen S Y, Yan J H, Peng P C, et al. Polarization-tracking-free PDM supporting hybrid digital-analog transport for fixed-mobile systems[J]. IEEE Photonics Technology Letters, 2019, 31(1): 54-57.
[3] Chen C, Zhuge Q B, Plant D V. Zero-guard-interval coherent optical OFDM with overlapped frequency-domain CD and PMD equalization[J]. Optics Express, 2011, 19(8): 7451-7467.
[4] ChenC, ZhugeQ, Plant DV. Reduced-guard-interval CO-OFDM with overlapped frequency-domain CD and PMD equalization[C]∥Optical Fiber Communication Conference/National Fiber Optic Engineers Conference 2011, March 6-10, 2011, Los Angeles, California, USA. Washington, DC: OSA, 2011: OWE7.
[5] Zheng Z B, Cui N, Xu H Y, et al. Window-split structured frequency domain Kalman equalization scheme for large PMD and ultra-fast RSOP in an optical coherent PDM-QPSK system[J]. Optics Express, 2018, 26(6): 7211-7226.
[6] Kudo R, Kobayashi T, Ishihara K, et al. PMD compensation in optical coherent single carrier transmission using frequency-domain equalisation[J]. Electronics Letters, 2009, 45(2): 124-125.
[7] Ishihara K, Kobayashi T, Kudo R, et al. Frequency-domain equalization for coherent optical single-carrier transmission systems[J]. IEICE Transactions on Communications, 2009, 92(12): 3736-3743.
[8] Godard D. Self-recovering equalization and carrier tracking in two-dimensional data communication systems[J]. IEEE Transactions on Communications, 1980, 28(11): 1867-1875.
[9] Ready MJ, Gooch RP. Blind equalization based on radius directed adaptation[C]∥International Conference on Acoustics, April 3-6, 1990, Albuquerque, NM, USA. New York: IEEE, 1990: 1699- 1702.
[12] Yang Y F, Cao G L, Zhong K P, et al. Fast polarization-state tracking scheme based on radius-directed linear Kalman filter[J]. Optics Express, 2015, 23(15): 19673-19680.
[13] Pakala L, Schmauss B. Joint tracking of polarization state and phase noise using adaptive cascaded Kalman filtering[J]. IEEE Photonics Technology Letters, 2017, 29(16): 1297-1300.
[14] 曹国亮, 杨彦甫, 王非, 等. 基于扩展卡尔曼的PDM-16QAM偏振态和载波相位快速跟踪[J]. 光学学报, 2014, 34(12): 1206005.
[15] Gordon J P, Kogelnik H. PMD fundamentals: polarization mode dispersion in optical fibers[J]. Proceedings of the National Academy of Sciences of the United States of America, 2000, 97(9): 4541-4550.
常雪峰, 黄斌, 张黎杰, 李霞, 隋琪. 基于互补斯托克斯向量实现偏振态旋转的追踪和预测[J]. 激光与光电子学进展, 2020, 57(21): 210601. Chang Xuefeng, Huang Bin, Zhang Lijie, Li Xia, Sui Qi. Tracking and Prediction of Polarization State Rotation Using Complementary Stokes Vectors[J]. Laser & Optoelectronics Progress, 2020, 57(21): 210601.