CN114244439B - 一种光纤通信系统的信号补偿方法及装置 - Google Patents

一种光纤通信系统的信号补偿方法及装置 Download PDF

Info

Publication number
CN114244439B
CN114244439B CN202110806092.2A CN202110806092A CN114244439B CN 114244439 B CN114244439 B CN 114244439B CN 202110806092 A CN202110806092 A CN 202110806092A CN 114244439 B CN114244439 B CN 114244439B
Authority
CN
China
Prior art keywords
signal
model
current
compensation
step length
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110806092.2A
Other languages
English (en)
Other versions
CN114244439A (zh
Inventor
赵启涵
张琦
忻向军
姚海鹏
田凤
田清华
张尼
蒋锦坤
王曦朔
王拥军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Posts and Telecommunications
6th Research Institute of China Electronics Corp
Original Assignee
Beijing University of Posts and Telecommunications
6th Research Institute of China Electronics Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Posts and Telecommunications, 6th Research Institute of China Electronics Corp filed Critical Beijing University of Posts and Telecommunications
Priority to CN202110806092.2A priority Critical patent/CN114244439B/zh
Publication of CN114244439A publication Critical patent/CN114244439A/zh
Application granted granted Critical
Publication of CN114244439B publication Critical patent/CN114244439B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/25Arrangements specific to fibre transmission
    • H04B10/2507Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion
    • H04B10/2513Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion due to chromatic dispersion
    • H04B10/2525Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion due to chromatic dispersion using dispersion-compensating fibres

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Optical Communication System (AREA)

Abstract

本发明提供了一种光纤通信系统的信号补偿方法及装置,方法包括:获取光纤通信系统的输出信号为待补偿信号;将待补偿信号输入补偿模型中,得到补偿后的信号;补偿模型包括多个MIMO滤波器模型;MIMO滤波器模型是根据光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的。本发明通过补偿模型能够同时对光纤通信系统的输出信号进行色散效应补偿、偏振模色散效应补偿和非线性效应补偿,提高信号补偿的精度。

Description

一种光纤通信系统的信号补偿方法及装置
技术领域
本发明涉及光纤通信技术领域,特别是涉及一种光纤通信系统的信号补偿方法及装置。
背景技术
随着信息化时代的到来,光纤通信系统朝着长距离、大容量的方向不断发展。偏振复用技术的引入使得光纤通信系统的传输容量得到成倍提升,但信道中产生的色散效应、偏振模色散效应和非线性效应会对传输的信号造成损伤,使系统的性能受到很大的限制。因此,对接收到的信号进行色散、偏振模色散、非线性效应的补偿是当前研究的热点。目前的补偿方案实现简单,不能同时对光纤通信系统色散、偏振模色散和非线性效应进行补偿,补偿效果差,难以精确恢复偏振复用光纤通信系统的原始信号,使得偏振复用光纤通信系统的使用受到限制。
发明内容
本发明的目的是提供一种光纤通信系统的信号补偿方法及装置,能够同时对光纤通信系统的输出信号进行色散效应补偿、偏振模色散效应补偿和非线性效应补偿,提高信号补偿的精度。
为实现上述目的,本发明提供了如下方案:
一种光纤通信系统的信号补偿方法,包括:
获取光纤通信系统的输出信号为待补偿信号;
将所述待补偿信号输入补偿模型中,得到补偿后的信号;所述补偿模型包括多个MIMO滤波器模型;所述MIMO滤波器模型是根据所述光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的。
可选的,在所述将所述待补偿信号输入补偿模型中,得到补偿后的信号之后,还包括:
利用相偏恢复算法,对所述补偿后的信号进行相偏补偿。
可选的,在所述获取光纤通信系统的输出信号为待补偿信号之前,还包括:
将所述光纤通信系统的传输距离划分为多个随机步长;
确定任一随机步长为当前随机步长;
根据所述当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型;
获取当前随机步长的光纤历史传输数据;所述当前随机步长的光纤历史传输数据包括原始信号,以及原始信号经当前随机步长的光纤传输后的输出信号;
以所述原始信号经当前随机步长的光纤传输后的输出信号为输入,以所述原始信号为输出,对所述当前MIMO滤波器模型进行训练,得到所述当前MIMO滤波器模型中多个抽头系数的抽头系数响应;
对多个所述抽头系数响应进行量化处理,得到多个抽头系数组;所述抽头系数组中的抽头系数的数量与所述当前MIMO滤波器模型中抽头系数的数量相等。
可选的,所述将所述光纤通信系统的传输距离划分为多个随机步长,具体包括:
将所述光纤通信系统的传输距离划分为多个跨段;
确定任一跨段为当前跨段;
利用公式
Figure GDA0004054617380000021
将所述当前跨段划分为多个对数步长;
根据所述当前跨段的多个对数步长,利用公式
Figure GDA0004054617380000022
得到所述当前跨段的多个随机步长;
其中,hn表示第n个对数步长,Γ表示光纤损耗,σ表示参量,σ=[1-exp(-2ΓL)]/n,exp(·)表示自然对数e的(·)次幂,L表示跨段长度,hr_n表示第n个随机步长,h1和hn-1分别表示第1个和第n-1个对数步长,r1、rn-1和rn分别表示第1个、第n-1个和第n个随机步长的随机系数,
Figure GDA0004054617380000023
rand(·)表示在(·)范围内进行随机取值,R表示最大扰动范围。
可选的,所述根据所述当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型,具体包括:
根据所述当前随机步长的色散效应,建立色散效应信号传播模型;
根据所述当前随机步长的偏振模色散效应,建立偏振模色散效应信号传播模型;
根据所述当前随机步长的非线性效应,建立非线性效应信号传播模型;
根据所述色散效应信号传播模型,确定第一线性算子;
根据所述偏振模色散效应信号传播模型,确定第二线性算子;
根据所述非线性效应信号传播模型,确定非线性算子;
根据所述第一线性算子、所述第二线性算子和所述非线性算子,确定当前随机步长的偏振复用信号传播模型;
根据所述偏振复用信号传播模型建立所述当前随机步长的反向传播模型;所述反向传播模型与所述偏振复用信号传播模型的结构相同;所述反向传播模型与所述偏振复用信号传播模型的信号传输方向相反;
根据所述反向传播模型,建立当前MIMO滤波器模型;所述当前MIMO滤波器模型的卷积用于描述所述反向传播模型的信号传输过程;所述当前MIMO滤波器模型的卷积是通过将色散补偿滤波器与差分群时延滤波器进行卷积处理后,乘以主偏振态旋转矩阵得到的;所述色散补偿滤波器、所述差分群时延滤波器和所述主偏振态旋转矩阵均是根据所述反向传播模型确定的。
可选的,
所述色散效应信号传播模型为:
Figure GDA0004054617380000031
所述偏振模色散效应信号传播模型为:
Figure GDA0004054617380000032
所述非线性效应信号传播模型为:
Figure GDA0004054617380000033
所述偏振复用信号传播模型为:
Figure GDA0004054617380000034
其中,
Figure GDA0004054617380000035
Figure GDA0004054617380000036
分别表示色散效应信号传播模型、偏振模色散效应信号传播模型和非线性效应信号传播模型,U表示位置为z、频率为ω的信号,j为虚数单位,β2表示色散效应的群速度时延,h表示步长,βk_PMD表示偏振模色散效应的群速度时延,
Figure GDA0004054617380000037
Rk表示主偏振态旋转矩阵,
Figure GDA0004054617380000038
τk表示差分群时延,σ1表示泡利矩阵,
Figure GDA0004054617380000039
u表示位置为z、时间为t的信号,γ表示克尔效应参数,
Figure GDA00040546173800000310
H表示u的厄密转置,u(tn,z+h)表示原始信号经当前随机步长的光纤传输后的输出信号,
Figure GDA00040546173800000311
表示第一线性算子,
Figure GDA00040546173800000312
表示第二线性算子,
Figure GDA00040546173800000313
表示非线性算子,u(tn,z)表示原始信号。
可选的,所述对多个所述抽头系数响应进行量化处理,得到多个抽头系数组,具体包括:
分别对每个所述抽头系数响应进行时隙划分,得到每个抽头系数的不同时隙内的抽头系数响应;
利用公式
Figure GDA0004054617380000041
对不同时隙内的抽头系数响应的实时值进行量化,得到多个抽头系数量化值;
将同一时刻的不同抽头系数响应的实时值对应的抽头系数量化值确定为同一抽头系数组,得到多个所述抽头系数组;
其中,
Figure GDA0004054617380000042
表示抽头系数量化值,cj1,m(k)表示第m个时隙内抽头系数的实时响应,
Figure GDA0004054617380000043
表示第m个时隙内抽头系数的响应最大值,j1表示抽头系数的序数,Δj表示量化因子,
Figure GDA0004054617380000044
表示时隙的个数,<·>表示四舍五入取整。
可选的,所述将所述待补偿信号输入补偿模型中,得到补偿后的信号,具体包括:
令迭代次数为i;i=1;
将所述待补偿信号确定为第i次输入信号;
将所述第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个所述抽头系数组设置为所述第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号;
将第i次补偿后的信号作为第i+1次迭代的输入信号,令i的数值增加1并返回步骤“将所述第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个所述抽头系数组设置为所述第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号”直至遍历所有MIMO滤波器模型,得到所述补偿后的信号。
一种光纤通信系统的信号补偿装置,包括:
待补偿信号获取模块,用于获取光纤通信系统的输出信号为待补偿信号;
信号补偿模块,用于将所述待补偿信号输入补偿模型中,得到补偿后的信号;所述补偿模型包括多个MIMO滤波器模型;所述MIMO滤波器模型是根据所述光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的。
可选的,所述装置,还包括:
相偏补偿模块,用于利用相偏恢复算法,对所述补偿后的信号进行相偏补偿。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种光纤通信系统的信号补偿方法及装置,方法包括:获取光纤通信系统的输出信号为待补偿信号;将待补偿信号输入补偿模型中,得到补偿后的信号;补偿模型包括多个MIMO(Multiple-InputMultiple-Output,多入多出技术)滤波器模型;MIMO滤波器模型是根据光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的。本发明通过补偿模型能够同时对光纤通信系统的输出信号进行色散效应补偿、偏振模色散效应补偿和非线性效应补偿,提高信号补偿的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中光纤通信系统的信号补偿方法流程图;
图2为本发明实施例中的信号补偿方法流程图;
图3为本发明实施例中MIMO滤波器模型的卷积示意图;
图4为本发明实施例中光纤通信系统的跨段划分示意图;
图5为本发明实施例中光纤通信系统的信号补偿装置结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种光纤通信系统的信号补偿方法及装置,能够同时对光纤通信系统的输出信号进行色散效应补偿、偏振模色散效应补偿和非线性效应补偿,提高信号补偿的精度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例中光纤通信系统的信号补偿方法流程图,图2为本发明实施例中的信号补偿方法流程图,如图1-2所示,本发明提供了一种光纤通信系统的信号补偿方法,包括:
步骤101:获取光纤通信系统的输出信号为待补偿信号;
步骤102:将待补偿信号输入补偿模型中,得到补偿后的信号;补偿模型包括多个MIMO滤波器模型;MIMO滤波器模型是根据光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的。
在步骤102之后,还包括:
利用相偏恢复算法,对补偿后的信号进行相偏补偿。
在步骤101之前,还包括:
将光纤通信系统的传输距离划分为多个随机步长;
确定任一随机步长为当前随机步长;
根据当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型;
获取当前随机步长的光纤历史传输数据;当前随机步长的光纤历史传输数据包括原始信号,以及原始信号经当前随机步长的光纤传输后的输出信号;
以原始信号经当前随机步长的光纤传输后的输出信号为输入,以原始信号为输出,对当前MIMO滤波器模型进行训练,得到当前MIMO滤波器模型中多个抽头系数的抽头系数响应;
对多个抽头系数响应进行量化处理,得到多个抽头系数组;抽头系数组中的抽头系数的数量与当前MIMO滤波器模型中抽头系数的数量相等。
其中,将光纤通信系统的传输距离划分为多个随机步长,具体包括:
将光纤通信系统的传输距离划分为多个跨段;
确定任一跨段为当前跨段;
利用公式
Figure GDA0004054617380000061
将当前跨段划分为多个对数步长;
根据当前跨段的多个对数步长,利用公式
Figure GDA0004054617380000071
得到当前跨段的多个随机步长;
其中,hn表示第n个对数步长,Γ表示光纤损耗,σ表示参量,σ=[1-exp(-2ΓL)]/n,exp(·)表示自然对数e的(·)次幂,L表示跨段长度,hr_n表示第n个随机步长,h1和hn-1分别表示第1个和第n-1个对数步长,r1、rn-1和rn分别表示第1个、第n-1个和第n个随机步长的随机系数,
Figure GDA0004054617380000072
rand(·)表示在(·)范围内进行随机取值,R表示最大扰动范围。本发明的跨段划分示意图如图4所示,图中z、z-L和z-kL表示不同跨段,虚曲线表示多个跨段的光功率随位置变化趋势。h1、h2、h3和h5表示第1个、第2个、第3个和第5个对数步长,hr-1、hr-2、hr-3和hr-5表示第1个、第2个、第3个和第5个对数步长。
根据当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型,具体包括:
根据当前随机步长的色散效应,建立色散效应信号传播模型;
根据当前随机步长的偏振模色散效应,建立偏振模色散效应信号传播模型;
根据当前随机步长的非线性效应,建立非线性效应信号传播模型;
根据色散效应信号传播模型,确定第一线性算子;
根据偏振模色散效应信号传播模型,确定第二线性算子;
根据非线性效应信号传播模型,确定非线性算子;
根据第一线性算子、第二线性算子和非线性算子,确定当前随机步长的偏振复用信号传播模型;
根据偏振复用信号传播模型建立当前随机步长的反向传播模型;反向传播模型与偏振复用信号传播模型的结构相同;反向传播模型与偏振复用信号传播模型的信号传输方向相反;
根据反向传播模型,建立当前MIMO滤波器模型;当前MIMO滤波器模型的卷积用于描述反向传播模型的信号传输过程;当前MIMO滤波器模型的卷积是通过将色散补偿滤波器与差分群时延滤波器进行卷积处理后,乘以主偏振态旋转矩阵得到的;色散补偿滤波器、差分群时延滤波器和主偏振态旋转矩阵均是根据反向传播模型确定的,MIMO滤波器模型的卷积如图3所示。
色散补偿滤波器为两个奇长复值对称滤波器,在两个偏振态中采用如下的相同的滤波器抽头系数:
Figure GDA0004054617380000081
其中,FCD为该滤波器的长度;cn表示色散补偿滤波器的抽头系数,
Figure GDA0004054617380000082
Figure GDA0004054617380000083
分别表示长度为
Figure GDA0004054617380000084
和0时色散补偿滤波器的抽头系数。
差分群时延滤波器为两个实值非对称滤波器,在两个偏振态中具有翻转的抽头系数
Figure GDA0004054617380000085
其中,FDGD为该滤波器长度。dn表示色散补偿滤波器的抽头系数,
Figure GDA0004054617380000086
Figure GDA0004054617380000087
分别表示长度为
Figure GDA0004054617380000088
和0时差分群时延滤波器的抽头系数。
在进行模型训练时,滤波器可进行两种参数化:一种为对全部抽头系数进行自由优化,具有与抽头个数相同的自由度;另一种为将滤波器约束为抽头系数如下的拉格朗日分数延迟滤波器:
Figure GDA0004054617380000089
其中,fs为中心频率,τk为差分群时延,F为非对称滤波器长度,根据该滤波器可以对τk直接优化,在每步只有一个自由度,与滤波器长度无关;
主偏振态旋转矩阵,为2x2复值琼斯矩阵,可进行两种参数化:一种为对主偏振态旋转矩阵中的全部项自由优化,具有8个自由度;另一种为采用如下二次酉矩阵对主偏振态旋转矩阵进行约束,具有4个自由度:
Figure GDA00040546173800000810
具体的,MIMO滤波器训练过程如下:将训练序列分成多个批次输入MIMO滤波器,对建立的基于机器学习的滤波器模型以有效信噪比的倒数为收敛目标使用小批次随机梯度下降法进行训练,输出当前信道匹配的抽头系数;损失函数为如下的有效信噪比的倒数:
Figure GDA00040546173800000811
其中,s表示发送符号序列,
Figure GDA00040546173800000812
表示接收符号序列。在此时的模型训练中,模型以有效信噪比最大化的目标进行收敛。最后通过引入动量和自适应速率的小批次随机梯度下降法对循行进行优化训练,得到抽头系数相应。
具体的,
色散效应信号传播模型为:
Figure GDA0004054617380000091
偏振模色散效应信号传播模型为:
Figure GDA0004054617380000092
非线性效应信号传播模型为:
Figure GDA0004054617380000093
偏振复用信号传播模型为
Figure GDA0004054617380000094
其中,
Figure GDA0004054617380000095
Figure GDA0004054617380000096
分别表示色散效应信号传播模型、偏振模色散效应信号传播模型和非线性效应信号传播模型,U表示位置为z、频率为ω的信号,j为虚数单位,β2表示色散效应的群速度时延,h表示步长,βk_PMD表示偏振模色散效应的群速度时延,
Figure GDA0004054617380000097
Rk表示主偏振态旋转矩阵,
Figure GDA0004054617380000098
τk表示差分群时延,σ1表示泡利矩阵,
Figure GDA0004054617380000099
u表示位置为z、时间为t的信号,γ表示克尔效应参数,
Figure GDA00040546173800000910
H表示u的厄密转置,u(tn,z+h)表示原始信号经当前随机步长的光纤传输后的输出信号,
Figure GDA00040546173800000911
表示第一线性算子,
Figure GDA00040546173800000912
表示第二线性算子,
Figure GDA00040546173800000913
表示非线性算子,u(tn,z)表示原始信号。
具体的,偏振复用光纤通信系统的信号传播模型,包括:
色散信号传播模型:设定位置为z,频率为ω的信号为U=U(z,ω),步长为h,(色散效应的)群速度时延为β2,根据以下表达式建立色散对光信号的频域影响模型:
Figure GDA00040546173800000914
非线性效应信号传播模型:设定位置为z,时间为t的信号为u=u(t,z),克尔效应参数为γ,根据以下表达式建立色非线性效应对光信号的时域影响模型:
Figure GDA00040546173800000915
其中,
Figure GDA00040546173800000916
上标H表示该值的厄密转置。
偏振模色散信号传播模型:根据以下表达式建立偏振模色散对光信号的影响模型:
Figure GDA0004054617380000101
其中,由偏振模色散引起的
Figure GDA0004054617380000102
可分解为:主偏振态旋转矩阵和差分群时延矩阵。
主偏振态旋转矩阵
Figure GDA0004054617380000103
为2阶特殊酉矩阵,表示信号的偏振态旋转;τk表示在两偏振态之间引入的差分群时延,
Figure GDA0004054617380000104
为泡利矩阵,定义了差分群时延对信号的影响。
综合上述三种模型,即得到如下的偏振复用信号传播模型:
Figure GDA0004054617380000105
其中,
Figure GDA0004054617380000106
为第一线性算子,与色散效应相关,根据色散效应信号传播模型对其后信号进行长度为h2的受色散效应影响下信号变化的运算得到;
Figure GDA0004054617380000107
为第二线性算子,与偏振模色散效应相关,根据偏振模色散效应信号传播模型对其后信号进行长度为h2的受偏振模色散效应影响下信号变化的运算得到;
Figure GDA0004054617380000108
为与非线性效应相关的非线性算子,根据非线性效应信号传播模型对其后信号进行长度为h的受非线性效应影响下信号变化的运算得到的。该模型表述了偏振复用信号在从z位置开始,经过长度为h光纤传输后的变化。
此外,对多个抽头系数响应进行量化处理,得到多个抽头系数组,具体包括:
分别对每个抽头系数响应进行时隙划分,得到每个抽头系数的不同时隙内的抽头系数响应;
利用公式
Figure GDA0004054617380000109
对不同时隙内的抽头系数响应的实时值进行量化,得到多个抽头系数量化值;
将同一时刻的不同抽头系数响应的实时值对应的抽头系数量化值确定为同一抽头系数组,得到多个所述抽头系数组;
其中,
Figure GDA00040546173800001010
表示抽头系数量化值,cj1,m(k)表示第m个时隙内抽头系数的实时响应,
Figure GDA00040546173800001011
表示第m个时隙内抽头系数的响应最大值,j1表示抽头系数的序数,Δj表示量化因子,
Figure GDA0004054617380000111
表示时隙的个数,<·>表示四舍五入取整。
具体的,对多个抽头系数响应进行量化过程如下:
时隙划分:根据幅值将滤波器抽头系数的脉冲响应划分为多个幅值差异较大的时隙,时隙之间的最大幅值通常相差5倍以上;
量化滤波器抽头系数的脉冲响应:将各时隙的滤波器抽头脉冲响应按照不同精度量化为离散值,完成滤波器的抽头系数配置。根据不同的响应幅度,可通过调整量化因子来改变量化精度,更大幅值的时隙需要更大量化精度来对脉冲响应进行准确地量化。
步骤102,具体包括:
令迭代次数为i;i=1;
将待补偿信号确定为第i次输入信号;
将第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个抽头系数组设置为第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号;
将第i次补偿后的信号作为第i+1次迭代的输入信号,令i的数值增加1并返回步骤“将第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个抽头系数组设置为第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号”直至遍历所有MIMO滤波器模型,得到补偿后的信号。
具体的,本发明包括以下步骤:
S1.对接收到的包含训练数据(光纤历史传输数据)和实际数据的信号进行低通滤波及下采样处理,得到携带信道损伤的待处理信号。其中,训练数据为接收端已知的固定序列,利用该序列在信道中产生的变化可以对当前的信道环境进行评估,并训练得到信号处理算法(即反向传播模型)的参数。
S2.通过信号中的训练数据对基于机器学习模型的MIMO滤波器进行训练,得到滤波器抽头系数脉冲响应。
S3.对滤波器的抽头系数脉冲响应进行分块量化,得到在不同时隙中精度不同的滤波器抽头系数。
S4.用已完成抽头系数配置的MIMO滤波器对接收到的实际数据进行滤波处理,即采用滤波器抽头系数对接收到的实际数据进行与MIMO滤波器结构相同的矩阵运算,补偿信号的色散效应、偏振模色散效应和非线性效应;
S5.对S4中得到的信号采用相偏恢复算法(例如盲相位恢复算法)进行相偏恢复,补偿激光器线宽引起的相位噪声,得到初始信号。
本发明提供的光纤通信系统的信号补偿方法及系统,基于现有补偿方案建模精度低,复杂度过高,量化误差相干累积,补偿性能低、不适用于偏振复用系统的问题。通过对光纤通信系统的进行跨段划分,每一个跨段划分为多个对数步长,在每个对数步长上添加上随机扰动,得到不同的随机步长,并针对每一个随机步长都建立一个MIMO滤波器模型,经总传输距离划分为多个小范围传输数据进行补偿,提高了补偿精度。
图5为本发明实施例中光纤通信系统的信号补偿装置结构示意图,如图5所示,本发明提供了一种光纤通信系统的信号补偿装置,包括:
待补偿信号获取模块501,用于获取光纤通信系统的输出信号为待补偿信号;
信号补偿模块502,用于将待补偿信号输入补偿模型中,得到补偿后的信号;补偿模型包括多个MIMO滤波器模型;MIMO滤波器模型是根据光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的。
此外,本发明提供的光纤通信系统的信号补偿装置,还包括:
相偏补偿模块,用于利用相偏恢复算法,对补偿后的信号进行相偏补偿。
随机步长划分模块,用于将光纤通信系统的传输距离划分为多个随机步长;
当前随机步长确定模块,用于确定任一随机步长为当前随机步长;
当前MIMO滤波器模型建立模块,用于根据当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型;
光纤历史传输数据获取模块,用于获取当前随机步长的光纤历史传输数据;当前随机步长的光纤历史传输数据包括原始信号,以及原始信号经当前随机步长的光纤传输后的输出信号;
当前MIMO滤波器模型训练模块,用于以原始信号经当前随机步长的光纤传输后的输出信号为输入,以原始信号为输出,对当前MIMO滤波器模型进行训练,得到当前MIMO滤波器模型中多个抽头系数的抽头系数响应;
抽头系数组确定模块,用于对多个抽头系数响应进行量化处理,得到多个抽头系数组;抽头系数组中的抽头系数的数量与当前MIMO滤波器模型中抽头系数的数量相等。
随机步长划分模块,具体包括:
跨段划分单元,用于将光纤通信系统的传输距离划分为多个跨段;
当前跨段单元,用于确定任一跨段为当前跨段;
对数步长确定单元,用于利用公式
Figure GDA0004054617380000131
将当前跨段划分为多个对数步长;
随机步长确定单元,用于根据当前跨段的多个对数步长,利用公式
Figure GDA0004054617380000132
得到当前跨段的多个随机步长;
其中,hn表示第n个对数步长,Γ表示光纤损耗,σ表示参量,σ=[1-exp(-2ΓL)]/n,exp(·)表示自然对数e的(·)次幂,L表示跨段长度,hr_n表示第n个随机步长,h1和hn-1分别表示第1个和第n-1个对数步长,r1、rn-1和rn分别表示第1个、第n-1个和第n个随机步长的随机系数,
Figure GDA0004054617380000133
rand(·)表示在(·)范围内进行随机取值,R表示最大扰动范围。
当前MIMO滤波器模型建立模块,具体包括:
色散效应信号传播模型建立单元,用于根据当前随机步长的色散效应,建立色散效应信号传播模型;
偏振模色散效应信号传播模型建立单元,用于根据当前随机步长的偏振模色散效应,建立偏振模色散效应信号传播模型;
非线性效应信号传播模型建立单元,用于根据当前随机步长的非线性效应,建立非线性效应信号传播模型;
第一线性算子确定单元,用于根据色散效应信号传播模型,确定第一线性算子;
第二线性算子确定单元,用于根据偏振模色散效应信号传播模型,确定第二线性算子;
非线性算子确定单元,用于根据非线性效应信号传播模型,确定非线性算子;
偏振复用信号传播模型确定单元,用于根据第一线性算子、第二线性算子和非线性算子,确定当前随机步长的偏振复用信号传播模型;
反向传播模型确定单元,用于根据偏振复用信号传播模型建立当前随机步长的反向传播模型;反向传播模型与偏振复用信号传播模型的结构相同;反向传播模型与偏振复用信号传播模型的信号传输方向相反;
当前MIMO滤波器模型确定单元,用于根据反向传播模型,建立当前MIMO滤波器模型;当前MIMO滤波器模型的卷积用于描述反向传播模型的信号传输过程;当前MIMO滤波器模型的卷积是通过将色散补偿滤波器与差分群时延滤波器进行卷积处理后,乘以主偏振态旋转矩阵得到的;色散补偿滤波器、差分群时延滤波器和主偏振态旋转矩阵均是根据反向传播模型确定的。
具体的,
色散效应信号传播模型为:
Figure GDA0004054617380000141
偏振模色散效应信号传播模型为:
Figure GDA0004054617380000142
非线性效应信号传播模型为:
Figure GDA0004054617380000143
偏振复用信号传播模型为:
Figure GDA0004054617380000144
其中,
Figure GDA0004054617380000145
Figure GDA0004054617380000146
分别表示色散效应信号传播模型、偏振模色散效应信号传播模型和非线性效应信号传播模型,U表示位置为z、频率为ω的信号,j为虚数单位,β2表示色散效应的群速度时延,h表示步长,βk_PMD表示偏振模色散效应的群速度时延,
Figure GDA0004054617380000147
Rk表示主偏振态旋转矩阵,
Figure GDA0004054617380000148
τk表示差分群时延,σ1表示泡利矩阵,
Figure GDA0004054617380000149
u表示位置为z、时间为t的信号,γ表示克尔效应参数,
Figure GDA00040546173800001410
H表示u的厄密转置,u(tn,z+h)表示原始信号经当前随机步长的光纤传输后的输出信号,
Figure GDA00040546173800001411
表示第一线性算子,
Figure GDA00040546173800001412
表示第二线性算子,
Figure GDA00040546173800001413
表示非线性算子,u(tn,z)表示原始信号。
抽头系数组确定模块,具体包括:
时隙划分单元,用于分别对每个抽头系数响应进行时隙划分,得到每个抽头系数的不同时隙内的抽头系数响应;
抽头系数量化值确定单元,用于利用公式
Figure GDA0004054617380000151
对不同时隙内的抽头系数响应的实时值进行量化,得到多个抽头系数量化值;
抽头系数组确定单元,用于将同一时刻的不同抽头系数响应的实时值对应的抽头系数量化值确定为同一抽头系数组,得到多个抽头系数组;
其中,
Figure GDA0004054617380000152
表示抽头系数量化值,cj1,m(k)表示第m个时隙内抽头系数的实时响应,
Figure GDA0004054617380000153
表示第m个时隙内抽头系数的响应最大值,j1表示抽头系数的序数,Δj表示量化因子,Mj1表示时隙的个数,<·>表示四舍五入取整。
信号补偿模块,具体包括:
迭代次数确定单元,用于令迭代次数为i;i=1;
第i次输入信号确定单元,用于将待补偿信号确定为第i次输入信号;
第i次补偿后的信号确定单元,用于将第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个抽头系数组设置为第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号;
输入信号更新单元,用于将第i次补偿后的信号作为第i+1次迭代的输入信号,令i的数值增加1并返回步骤“将第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个抽头系数组设置为第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号”直至遍历所有MIMO滤波器模型,得到补偿后的信号。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种光纤通信系统的信号补偿方法,其特征在于,所述方法,包括:
获取光纤通信系统的输出信号为待补偿信号;
将所述待补偿信号输入补偿模型中,得到补偿后的信号;所述补偿模型包括多个MIMO滤波器模型;所述MIMO滤波器模型是根据所述光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的;
在所述获取光纤通信系统的输出信号为待补偿信号之前,还包括:
将所述光纤通信系统的传输距离划分为多个随机步长;
确定任一随机步长为当前随机步长;
根据所述当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型;
获取当前随机步长的光纤历史传输数据;所述当前随机步长的光纤历史传输数据包括原始信号,以及原始信号经当前随机步长的光纤传输后的输出信号;
以所述原始信号经当前随机步长的光纤传输后的输出信号为输入,以所述原始信号为输出,对所述当前MIMO滤波器模型进行训练,得到所述当前MIMO滤波器模型中多个抽头系数的抽头系数响应;
对多个所述抽头系数响应进行量化处理,得到多个抽头系数组;所述抽头系数组中的抽头系数的数量与所述当前MIMO滤波器模型中抽头系数的数量相等。
2.根据权利要求1所述的光纤通信系统的信号补偿方法,其特征在于,在所述将所述待补偿信号输入补偿模型中,得到补偿后的信号之后,还包括:
利用相偏恢复算法,对所述补偿后的信号进行相偏补偿。
3.根据权利要求1所述的光纤通信系统的信号补偿方法,其特征在于,所述将所述光纤通信系统的传输距离划分为多个随机步长,具体包括:
将所述光纤通信系统的传输距离划分为多个跨段;
确定任一所述跨段为当前跨段;
利用公式
Figure FDA0004075779550000011
将所述当前跨段划分为多个对数步长;
根据所述当前跨段的多个对数步长,利用公式
Figure FDA0004075779550000012
得到所述当前跨段的多个随机步长;
其中,hn表示第n个对数步长,Γ表示光纤损耗,σ表示参量,σ=[1-exp(-2ΓL)]/n,exp(·)表示自然对数e的(·)次幂,L表示跨段长度,hr_n表示第n个随机步长,h1和hn-1分别表示第1个和第n-1个对数步长,r1、rn-1和rn分别表示第1个、第n-1个和第n个随机步长的随机系数,
Figure FDA0004075779550000021
rand(·)表示在(·)范围内进行随机取值,R表示最大扰动范围。
4.根据权利要求1所述的光纤通信系统的信号补偿方法,其特征在于,所述根据所述当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型,具体包括:
根据所述当前随机步长的色散效应,建立色散效应信号传播模型;
根据所述当前随机步长的偏振模色散效应,建立偏振模色散效应信号传播模型;
根据所述当前随机步长的非线性效应,建立非线性效应信号传播模型;
根据所述色散效应信号传播模型,确定第一线性算子;
根据所述偏振模色散效应信号传播模型,确定第二线性算子;
根据所述非线性效应信号传播模型,确定非线性算子;
根据所述第一线性算子、所述第二线性算子和所述非线性算子,确定当前随机步长的偏振复用信号传播模型;
根据所述偏振复用信号传播模型建立所述当前随机步长的反向传播模型;所述反向传播模型与所述偏振复用信号传播模型的结构相同;所述反向传播模型与所述偏振复用信号传播模型的信号传输方向相反;
根据所述反向传播模型,建立当前MIMO滤波器模型;所述当前MIMO滤波器模型的卷积用于描述所述反向传播模型的信号传输过程;所述当前MIMO滤波器模型的卷积是通过将色散补偿滤波器与差分群时延滤波器进行卷积处理后,乘以主偏振态旋转矩阵得到的;所述色散补偿滤波器、所述差分群时延滤波器和所述主偏振态旋转矩阵均是根据所述反向传播模型确定的。
5.根据权利要求4所述的光纤通信系统的信号补偿方法,其特征在于,
所述色散效应信号传播模型为:
Figure FDA0004075779550000022
所述偏振模色散效应信号传播模型为:
Figure FDA0004075779550000023
所述非线性效应信号传播模型为:
Figure FDA0004075779550000031
所述偏振复用信号传播模型为:
Figure FDA0004075779550000032
其中,
Figure FDA0004075779550000033
Figure FDA0004075779550000034
分别表示色散效应信号传播模型、偏振模色散效应信号传播模型和非线性效应信号传播模型,U表示位置为z、频率为ω的信号,j为虚数单位,β2表示色散效应的群速度时延,h表示步长,βk_PMD表示偏振模色散效应的群速度时延,
Figure FDA0004075779550000035
Rk表示主偏振态旋转矩阵,
Figure FDA0004075779550000036
τk表示差分群时延,σ1表示泡利矩阵,
Figure FDA0004075779550000037
u表示位置为z、时间为t的信号,γ表示克尔效应参数,
Figure FDA0004075779550000038
H表示u的厄密转置,u(tn,z+h)表示原始信号经当前随机步长的光纤传输后的输出信号,
Figure FDA0004075779550000039
表示第一线性算子,
Figure FDA00040757795500000310
表示第二线性算子,
Figure FDA00040757795500000311
表示非线性算子,u(tn,z)表示原始信号。
6.根据权利要求4所述的光纤通信系统的信号补偿方法,其特征在于,所述对多个所述抽头系数响应进行量化处理,得到多个抽头系数组,具体包括:
分别对每个所述抽头系数响应进行时隙划分,得到每个抽头系数的不同时隙内的抽头系数响应;
利用公式
Figure FDA00040757795500000312
对不同时隙内的抽头系数响应的实时值进行量化,得到多个抽头系数量化值;
将同一时刻的不同抽头系数响应的实时值对应的抽头系数量化值确定为同一抽头系数组,得到多个所述抽头系数组;
其中,
Figure FDA00040757795500000313
表示抽头系数量化值,
Figure FDA00040757795500000314
表示第m个时隙内抽头系数的实时响应,
Figure FDA00040757795500000315
表示第m个时隙内抽头系数的响应最大值,j1表示抽头系数的序数,Δj表示量化因子,
Figure FDA00040757795500000316
表示时隙的个数,<·>表示四舍五入取整。
7.根据权利要求6所述的光纤通信系统的信号补偿方法,其特征在于,所述将所述待补偿信号输入补偿模型中,得到补偿后的信号,具体包括:
令迭代次数为i;i=1;
将所述待补偿信号确定为第i次输入信号;
将所述第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个所述抽头系数组设置为所述第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号;
将第i次补偿后的信号作为第i+1次迭代的输入信号,令i的数值增加1并返回步骤“将所述第i次输入信号输入第i个随机步长对应MIMO滤波器模型中,按照时间顺序将第i个随机步长对应MIMO滤波器模型的多个所述抽头系数组设置为所述第i个随机步长对应MIMO滤波器模型的抽头系数,得到第i次补偿后的信号”直至遍历所有MIMO滤波器模型,得到所述补偿后的信号。
8.一种光纤通信系统的信号补偿装置,其特征在于,所述装置,包括:
待补偿信号获取模块,用于获取光纤通信系统的输出信号为待补偿信号;
信号补偿模块,用于将所述待补偿信号输入补偿模型中,得到补偿后的信号;所述补偿模型包括多个MIMO滤波器模型;所述MIMO滤波器模型是根据所述光纤通信系统的色散效应、偏振模色散效应和非线性效应建立的;
随机步长划分模块,用于将光纤通信系统的传输距离划分为多个随机步长;
当前随机步长确定模块,用于确定任一随机步长为当前随机步长;
当前MIMO滤波器模型建立模块,用于根据当前随机步长的色散效应、偏振模色散效应和非线性效应,建立当前MIMO滤波器模型;
光纤历史传输数据获取模块,用于获取当前随机步长的光纤历史传输数据;当前随机步长的光纤历史传输数据包括原始信号,以及原始信号经当前随机步长的光纤传输后的输出信号;
当前MIMO滤波器模型训练模块,用于以原始信号经当前随机步长的光纤传输后的输出信号为输入,以原始信号为输出,对当前MIMO滤波器模型进行训练,得到当前MIMO滤波器模型中多个抽头系数的抽头系数响应;
抽头系数组确定模块,用于对多个抽头系数响应进行量化处理,得到多个抽头系数组;抽头系数组中的抽头系数的数量与当前MIMO滤波器模型中抽头系数的数量相等。
9.根据权利要求8所述的光纤通信系统的信号补偿装置,其特征在于,所述装置,还包括:
相偏补偿模块,用于利用相偏恢复算法,对所述补偿后的信号进行相偏补偿。
CN202110806092.2A 2021-07-16 2021-07-16 一种光纤通信系统的信号补偿方法及装置 Active CN114244439B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110806092.2A CN114244439B (zh) 2021-07-16 2021-07-16 一种光纤通信系统的信号补偿方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110806092.2A CN114244439B (zh) 2021-07-16 2021-07-16 一种光纤通信系统的信号补偿方法及装置

Publications (2)

Publication Number Publication Date
CN114244439A CN114244439A (zh) 2022-03-25
CN114244439B true CN114244439B (zh) 2023-03-21

Family

ID=80742831

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110806092.2A Active CN114244439B (zh) 2021-07-16 2021-07-16 一种光纤通信系统的信号补偿方法及装置

Country Status (1)

Country Link
CN (1) CN114244439B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106452593A (zh) * 2016-05-30 2017-02-22 北京邮电大学 滤波器的构建方法和装置、及非线性噪声抑制方法和系统
CN106788740A (zh) * 2017-02-21 2017-05-31 武汉邮电科学研究院 光纤通信系统中的自适应非线性均衡器及方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102983910B (zh) * 2012-11-19 2017-06-23 中兴通讯股份有限公司 相干光通信系统中色散和非线性补偿方法及系统
CN107809282B (zh) * 2017-10-13 2019-08-16 北京邮电大学 相干光纤偏分复用系统中极端场景下的均衡方法及设备
CN111181655A (zh) * 2020-01-06 2020-05-19 同济大学 一种基于bp算法和非线性dfe算法的光接收机

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106452593A (zh) * 2016-05-30 2017-02-22 北京邮电大学 滤波器的构建方法和装置、及非线性噪声抑制方法和系统
CN106788740A (zh) * 2017-02-21 2017-05-31 武汉邮电科学研究院 光纤通信系统中的自适应非线性均衡器及方法

Also Published As

Publication number Publication date
CN114244439A (zh) 2022-03-25

Similar Documents

Publication Publication Date Title
CN101965698B (zh) 色散均衡装置和方法以及数字相干光接收机
Diniz Adaptive filtering
EP1430675B1 (en) Fast computation of mimo decision feedback equalizer coefficients
CN104579243B (zh) 一种dft调制滤波器组的快速设计方法
US20120114341A1 (en) System and method for frequency-domain chromatic dispersion and polarization mode dispersion compensation with time-domain channel estimation
CN109450552B (zh) 一种基于cma的信号均衡方法及系统
CN102160336B (zh) 仿真装置和仿真方法
CN107534484B (zh) 光接收装置及光接收方法
CN102835040B (zh) 微波mimo系统中接收信号的均衡方法和均衡器
CN103384901B (zh) 改进变换域中的滤波
CN114244439B (zh) 一种光纤通信系统的信号补偿方法及装置
GB1578938A (en) Frequency domain equalizer
Fareed et al. Channel impulse response equalization based on genetic algorithm in mode division multiplexing
Guiomar et al. Simplified Volterra series nonlinear equalizer by intra-channel cross-phase modulation oriented pruning
EP3143708B1 (en) Non-linear propagation impairment equalisation in optical transmission
Qiu et al. Fourier Neural Operator Based Fibre Channel Modelling for Optical Transmission
Glentis et al. Electronic dispersion compensation of fiber links using sparsity induced volterra equalizers
CN107005307B (zh) 一种设置均衡装置的方法及均衡装置
Bai et al. Preconditioning multichannel adaptive filtering algorithms using EVD-and SVD-based signal prewhitening and system decoupling
CN1312867C (zh) 偏振模式色散补偿器及其自适应调节方法
Zhang et al. Machine learning for long-haul optical systems
Shulkind et al. Estimating the Volterra series transfer function over coherent optical OFDM for efficient monitoring of the fiber channel nonlinearity
Zheng Implementation of an adaptive equalizer based on machine learning for real-time coherent optical receivers
CN102347921A (zh) 接收器以及对接收信号进行均衡处理的方法
Foster et al. A novel algorithm for calculating the QR decomposition of a polynomial matrix

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant