激光与光电子学进展, 2020, 57 (21): 210601, 网络出版: 2020-10-24  

基于互补斯托克斯向量实现偏振态旋转的追踪和预测 下载: 985次

Tracking and Prediction of Polarization State Rotation Using Complementary Stokes Vectors
作者单位
暨南大学光子技术研究院, 广东 广州 511486
摘要
针对偏振复用相干光通信系统链路中存在的高速偏振态旋转造成的信道损伤、相干接收机的动态均衡算法模块在实际系统中存在的反馈延迟导致的偏振态追踪性能不足的问题,提出一种基于互补斯托克斯向量的偏振态旋转的追踪和预测算法。所提算法将测得的琼斯空间中的偏振旋转矩阵和其逆矩阵变换到斯托克斯空间,并计算得到一对互补向量,并对斯托克斯空间中的互补向量的运动轨迹进行追踪和预测,以实现链路偏振态旋转的补偿。仿真结果表明,在均衡反馈延迟为5帧时,通过对互补向量进行循迹,预测的琼斯矩阵的平方误差可比原始方案降低10%。通过采用所提算法,系统的偏振旋转速度容限从1.3 MHz提升至1.5 MHz。
Abstract
The high-speed polarization state rotation associated with the polarization multiplexing coherent optical communication system link causes channel damage. Further, the feedback delay associated with the dynamic equalization algorithm module of the coherent receiver in the actual system results in insufficient polarization state tracking performance. In this study, a tracking and prediction algorithm for polarization state rotation is proposed based on complementary Stokes vectors to resolve the aforementioned problems. The proposed algorithm initially transforms the measured polarization rotation matrix of the Jones space and its inverse matrix into the Stokes space, obtains a pair of complementary vectors, and tracks and predicts the motion trajectories of the complementary vectors in Stokes space. Finally, the compensation of the polarization state rotation of the link is realized. The simulation results show that the squared error of the predicted Jones matrix can be reduced by 10% compared with that in the original scheme by tracking the complementary vector when the equalization feedback delay is 5 frames. Furthermore, the polarization rotation speed tolerance of the system increases from 1.3 MHz to 1.5 MHz when using the proposed algorithm.

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为厄米特转置,则UV可写为

U=a'+jb'c'+jd'-c'+jd'a'-jb',V=a'-jb'-c'-jd'c'-jd'a'+jb',(1)

式中: a'2+b'2+c'2+d'2=1,因此,UV的空间为四维空间中的单位球面。偏振态旋转在斯托克斯空间中表示为一个向量。令uv分别表示UV的斯托克斯空间向量,uv分别为uv的模长,则Uu的转换公式[15]

U=Icosu-jsinu=e-j,(2)

式中:I为二维单位矩阵;σ为泡利旋转矩阵。

σ=σ1σ2σ3,(3)

其中,

σ1=100-1,σ2=0110,σ3=0-jj0,(4)

式中:uv的空间R为半径为π的三维球内部。通过(2)、(3)式不难推出,uv方向相反,且u+v=π,因此称uv为互补斯托克斯向量。

3 偏振旋转矩阵在斯托克斯空间的运动轨迹

RSOP在斯托克斯空间中可以表示为一个点的旋转,目前其他工作对RSOP转速的仿真都是恒定转速的。仿真模型中使用的是多段合成的随机旋转速度,在其中的每一段中,所提RSOP模型是随机方向的静态旋转模型,每段的转速服从高斯分布,方差为σ2。所以,为了分析在斯托克斯空间中RSOP的运动方式,需要先推导出RSOP的等效转速。球表面积S的表达式为

S=20π22πrcosθrdθ=4πr2,(5)

式中:r为球的半径;θ为旋转角度。球面上任一点绕任意轴旋转一周的长度均值为

L=20π2(2πrcosθ)24πr2rdθ=2πr·π4,(6)

即针对赤道上一点的旋转长度,系数为 π4。先求球面上任一点绕任意轴旋转一周的周长的平方,再对其求均值可以得到

L2=20π2(2πrcosθ)34πr2rdθ=4π2r2·23,(7)

即针对赤道上一点的旋转长度的平方均值,系数为 23。在时间上将信号分段,每段的符号数为N,每段在该点的速度贡献之和呈Rayleigh分布,则每个维度方差为

σ0=12Nσ2·23=13Nσ2(8)

所以速度均值为

V=Nσ2/3·π/2=πNσ2/6(9)

对于单位转速的等效平均速度,表达式为

V=πNσ26·4π=8πNσ23(10)

从(6)~(10)式得出了在等效的斯托克斯空间中的运动速度和实际的RSOP转速之间的关系。结合(1)、(2)式,即可仿真出RSOP,图1为RSOP矩阵在斯托克斯空间的轨迹图,S1、S2、S3为斯托克斯空间的三个维度的坐标,RSOP在斯托克斯空间的运动相对平滑,可以对其进行预测。

图1看出,RSOP轨迹可以在局部拟合成一个圆,所以可以通过圆拟合方案来对RSOP进行追踪和预测,即通过已知的三个点的位置来确定一个圆形,然后通过圆的轨迹来预测接下来RSOP点的位置。圆拟合方案如图2(a)所示,其中d为延迟时间。

图 1. 斯托克斯空间中RSOP矩阵运动轨迹图

Fig. 1. Motion trajectory of RSOP matrix in Stokes space

下载图片 查看所有图片

图 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的运动轨迹,如图2(b)所示,根据向量的前进位移和速度就可以计算出向量的前进速度,从而预测出RSOP点的位置。

为了测试圆拟合方案和向量追踪方案的效果,估计使用和不使用互补向量算法的差别,仿真了1000个连续RSOP矩阵,测试了延迟为5帧时的效果,使用琼斯矩阵参数的平方误差来衡量算法的效果。图3为圆拟合方案和向量追踪方案的琼斯矩阵平方误差对比图,可以看到,在大多数时间内,圆拟合方案的平方误差大于向量追踪方案,然而,向量追踪方案在某些时刻的平方误差会特别大,这种误差可以用互补向量追踪算法来修正。

图 3. 圆拟合方案和向量追踪方案的琼斯矩阵平方误差曲线

Fig. 3. Square error curves of Jones matrix obtained by circle fitting scheme and vector tracking scheme

下载图片 查看所有图片

uv两个向量都是RSOP向量,分别使用向量追踪算法对两个互补向量进行了追踪。发现两个向量得到的平方误差曲线,一条误差较大时,另一条的误差则较小。因为uv两个向量为互补向量,将uv构成的空间R进一步分为半径为π/2的球内部和半径在π/2至π之间的空心球,分别称它们为内球和外球,则对于任意旋转矩阵,uv必然一个在内球一个在外球。空间R为非线性空间,R的表面均对应U=-I,而其球心位置附近可近似为线性空间,因此,对于时变的U和对应的V,选择位于内球的一个矩阵进行RSOP的追踪和预测。由此提出了一种互补向量RSOP算法。

使用和不使用自动选择互补向量算法的平方误差曲线如图4所示,其中,连续RSOP旋转矩阵的仿真点数为1000个,延迟为5帧。可以看到,互补向量追踪算法在大部分时间都可以得到具有较小平方误差的展开点,对比单独使用一个展开点的算法,此算法的平方误差始终是一个相对较小的值,而单独使用一个展开点的算法则存在突然出现巨大误差的情况。实验结果表明,互补向量算法是明显优于不使用互补向量算法的。

图 4. 互补向量追踪算法和向量追踪算法的平方误差曲线

Fig. 4. Square error curves of vector tracking algorithm and complementary vector tracking algorithm

下载图片 查看所有图片

4 算法原理与流程

对LMS算法和互补向量RSOP矩阵预测算法进行对比。LMS算法的流程如图5所示,该算法结合了载波相位恢复(CPE)、快速傅里叶变换(FFT)、快速傅里叶逆变换(IFFT)、最小二乘(LS)等算法。其中,ri为输入的信号,StSf经过傅里叶逆变换后的时域信号,So为判决处理后的输出信号。

图 5. LMS算法流程图

Fig. 5. Flow chart of LMS algorithm

下载图片 查看所有图片

图5可以看到,算法流程为信号经过频域LMS进行均衡,判决后对得到的信号再进行LS算法计算,得到RSOP矩阵和载波相位回复需要的旋转角后再进行载波相位恢复。

LS算法表达式为

M=SHofrf*rHfrf*,(11)

式中:Sf为经LMS算法后的频域信号;rf为经LMS算法之前的频域信号;M为二阶矩阵。

首先任意一个二阶矩阵可以表示为

M=ejφM1M2,(12)

式中:φ为相位偏转角;M2为形如(1)式的偏振旋转矩阵;厄米特矩阵M1可以表示为

M1=ac+jdc-jdb(13)

显然,M1与其厄米特共轭矩阵相等,所以将M与其厄米特共轭矩阵相乘,可得

MMH=M12,(14)

再将(14)式代入(12)式,即可得

ejφM2=(MMH)-12M(15)

即得到了用于RSOP预测的偏振旋转矩阵M2,用琼斯矩阵来进行RSOP的预测,另外因为琼斯矩阵的模为1,所以可以用模值计算出相位旋转角,用来进行载波相位恢复。

RSOP预测算法流程如图6所示,相比图5,在LS算法后加入了选择是否使用互补向量的算法,在LMS算法前使用预测得到的矩阵进行RSOP均衡。RSOP矩阵可以用斯托克斯空间的三维向量V表示,且它的分布在斯托克斯空间中构成了一个三维球体,则追踪和预测RSOP的问题就转换为了追踪三维空间中点的运动轨迹的问题。因此提出了通过追踪RSOP在斯托克斯空间上的速度和位移来对其进行预测的方法。

图 6. RSOP预测算法流程图

Fig. 6. Flow chart of RSOP prediction algorithm

下载图片 查看所有图片

RSOP在斯托克斯空间中的位移可以分为两个方向来计算,分别是沿着当前速度方向的位移X0和旋转方向的位移X1,总位移为两个方向的位移之和。先计算单位延时后的RSOP在斯托克斯空间的位移,假设已知当前时间和上一时间的RSOP点的位置,则可以计算出当前时间k的速度 vk,并可以得到当前速度的模值 vk和方向向量 vkvk,对当前时间的位移向量和上一时间的位移向量进行叉乘,再计算出旋转轴l,则通过速度位移公式可以计

算下一时间的两个方向的位移,取d个时间的位移进行累加,即可计算得到d时间之后点的轨迹。为了简化编程,使用泰勒展开对计算过程进行化简,化简后的计算公式为

Tk+d=dvkcos(n2ta)+n1avcos(n4ta)·vkvk+dvksinn1tad+n1avsin(n3ta)u+vk,(16)

式中:T为RSOP在斯托克斯空间的轨迹点;vk为速度;av为加速度;ta为两个速度向量之间的旋转角度;n1,n2,n3n4是为了简化积分计算,通过泰勒展开得到的参数,计算方法为

n1=d(d+1)2,n2=1dd2d,n3=1dd2n1,n4=1dd3n1(17)

计算d个延迟后轨迹点的位置后,再对更新系数进行更新,并计算下一时间的速度、加速度、旋转方向、移动方向等。通过旋转方向和速度的预测,可以计算出偏振态旋转的琼斯矩阵,即可对偏振态旋转进行均衡。

5 仿真结果

绘制如图7所示的斯托克斯空间中RSOP轨迹预测的效果图,以直观地观察互补向量算法对RSOP的追踪效果,可以看出,在斯托克斯空间中,预测点的轨迹始终跟随实际的RSOP轨迹。结果表明,所提算法可以较好地对斯托克斯空间中的RSOP运动轨迹进行追踪,使用所提算法预测RSOP是可行的。

图 7. 斯托克斯空间RSOP轨迹预测

Fig. 7. RSOP trajectory prediction in Stokes space

下载图片 查看所有图片

在不同的RSOP旋转速度下,对使用互补向量的RSOP矩阵预测算法和不使用RSOP只使用LMS的算法的BER和EVM进行对比。仿真调制格式为16QAM,Block长度为256,结果如图8图9所示。

图 8. EVM对比

Fig. 8. EVM comparison

下载图片 查看所有图片

图 9. BER对比

Fig. 9. BER comparison

下载图片 查看所有图片

图8可知,LMS算法的EVM从RSOP转速大于1.3 MHz后开始有明显的上升,而RSOP预测算法则在RSOP转速大于1.5 MHz时才开始有明显的上升。从图9可知,在RSOP转速超过1.3 MHz前,LMS算法的曲线始终小于16QAM的BER容限(3.8×10-3),但是,大部分都在10-3量级,而RSOP预测算法的曲线在1.5 MHz前都是10-5量级,即通过采用互补向量的循迹追踪算法,系统的偏振态旋转速度容限从1.3 MHz提升至1.5 MHz。预测得EVM可以降低约10%,BER可以从10-3量级降低到10-5量级。

6 结论

首先分析RSOP在琼斯空间和斯托克斯空间之间的转换,提出一种将RSOP矩阵转换到斯托克斯空间中的算法,然后,在斯托克斯空间中使用多段合成的随机旋转速度对RSOP进行建模。依据此模型,提出了两种圆拟合和计算旋转点的向量速度加速度的预测方案,并对RSOP矩阵进行预测,最终经过对比,选择了计算向量运动的方式,并且为了提高算法精度,提出了互补向量算法。经过对比,使用互补向量算法的平方误差明显小于不使用互补向量的算法。从RSOP轨迹的追踪图中可以看出,所提算法可以较好地预测RSOP矩阵的运动轨迹,使解决实际通信系统中的延迟问题成为可能。最后,对比了互补向量预测算法和LMS算法的EVM和BER曲线。经过计算可得,互补向量的寻迹追踪算法的EVM可以降低10%,且BER可从10-3量级降低到10-5量级,偏振态旋转的速度容限从1.3 MHz提升至1.5 MHz。由此可知,所提算法可以较好地预测RSOP矩阵,通信效果有明显提升,具有良好的应用前景。

参考文献

[1] Ghazisaeidi A, Sab O A, et al. 25.4-Tb/s transmission over transpacific distances using truncated probabilistically shaped PDM-64QAM[J]. Journal of Lightwave Technology, 2018, 36(6): 1354-1361.

[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.

[10] Yu Z M, Yi X W, Zhang J, et al. Modified constant modulus algorithm with polarization demultiplexing in Stokes space in optical coherent receiver[J]. Journal of Lightwave Technology, 2013, 31(19): 3203-3209.

[11] Noe R, Hinz S, Sandel D, et al. Crosstalk detection schemes for polarization division multiplex transmission[J]. Journal of Lightwave Technology, 2001, 19(10): 1469-1475.

[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.

    Cao G L, Yang Y F, Wang F, et al. Extended Kalman based polarization and carrier phase quickly tracking for PDM-16QAM[J]. Acta Optica Sinica, 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.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!