CN106301593B - 自适应盲偏振解复用处理方法和装置 - Google Patents

自适应盲偏振解复用处理方法和装置 Download PDF

Info

Publication number
CN106301593B
CN106301593B CN201610650339.5A CN201610650339A CN106301593B CN 106301593 B CN106301593 B CN 106301593B CN 201610650339 A CN201610650339 A CN 201610650339A CN 106301593 B CN106301593 B CN 106301593B
Authority
CN
China
Prior art keywords
signal
moment
polarisation state
sign
solid part
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
CN201610650339.5A
Other languages
English (en)
Other versions
CN106301593A (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
Original Assignee
Beijing University of Posts and Telecommunications
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 filed Critical Beijing University of Posts and Telecommunications
Priority to CN201610650339.5A priority Critical patent/CN106301593B/zh
Publication of CN106301593A publication Critical patent/CN106301593A/zh
Application granted granted Critical
Publication of CN106301593B publication Critical patent/CN106301593B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60—Receivers
    • H04B10/61—Coherent receivers
    • H04B10/616—Details of the electronic signal processing in coherent optical receivers
    • H04B10/6165—Estimation of the phase of the received optical signal, phase error estimation or phase error correction
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04J—MULTIPLEX COMMUNICATION
    • H04J14/00—Optical multiplex systems
    • H04J14/06—Polarisation multiplex systems

Landscapes

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

Abstract

本发明涉及一种自适应盲偏振解复用处理方法和装置，该方法包括：确定自适应FIR滤波器在当前时刻的抽头系数；根据所述当前时刻的抽头系数，所述自适应FIR滤波器对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；其中，根据所述上一时刻的抽头系数和两个偏振态复数信号，计算所述当前时刻的抽头系数。由于本发明通过分解、坐标变换可以使不同调制格式的多模信号变换到同一个圆上成为恒模信号，因此可以在不知道信号调制格式的前提下对其进行偏振解复用，因此本发明提供的方法与调制格式无关，适用于网络中调制格式愈加灵活，流量愈加动态、多样和不可预测的应用场景。

Description

自适应盲偏振解复用处理方法和装置
技术领域
本发明涉及光纤通信技术领域，尤其是涉及一种自适应盲偏振解复用处理方法和装置。
背景技术
近年来，随着网络流量和带宽需求的急剧增加，高速相干光通信技术已成为实现长距离大容量信息传输的关键技术。通过偏振复用和波分复用等技术，可以有效的提高光纤通信***的带宽利用率，提升***容量。对高速相干光通信而言，由于光纤色散、偏振随机串扰等因素的影响，因此需要采用多抽头的自适应蝶形FIR滤波器进行残余色散、偏振膜色散的均衡和偏振解复用。
目前，偏振解复用方法包括数据辅助和盲处理两种。其中，数据辅助方法通过发送训练序列来进行信道估计，占用了部分频谱资源，降低了***带宽利用率。然而，盲处理方法作为偏分复用光通信解调***中的一项关键技术，无需周期性地发送训练序列，可仅仅依靠接收信号的统计特性来恢复原始发送信号，因此无需占用太多带宽资源。可见两种方法相比，盲处理方法尤显优势。由于数据中心和云计算的不断发展，网络流量将变得愈加动态、多样和不可预测，支持多种调制格式的弹性光收发机成为弹性光网络中有效承载突发动态变化网络流量、提高频谱利用率、优化网络资源利用的关键。因此有必要提供一种与调制格式无关的盲处理方法。
发明内容
针对以上缺陷，本发明提供一种自适应盲偏振解复用处理方法和装置，处理过程与调制格式无关，适用于网络流量愈加动态、多样和不可预测的应用场景。
本发明提供的自适应盲偏振解复用处理方法包括：
确定自适应FIR滤波器在当前时刻的抽头系数；
根据所述当前时刻的抽头系数，所述自适应FIR滤波器对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；
其中，所述确定自适应FIR滤波器在当前时刻的抽头系数，包括：
对所述自适应FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；
分别对该偏振态复数信号的实部信号和虚部信号进行坐标变换；
根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；
将该实部信号的误差和该虚部信号的误差进行合成，得到该偏振态复数信号的复数误差；
根据所述上一时刻的抽头系数和各个偏振态复数信号，计算所述当前时刻的抽头系数。
可选的，所述自适应蝶形FIR滤波器在每一时刻的输出信号中包括x偏振态复数信号和y偏振态复数信号；对应的，
采用下式对x偏振态复数信号的实部信号进行坐标变换：
x′i(k)＝xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
其中，xi(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中x偏振态复数信号的实部信号，x′i(k)为对实部信号xi(k)坐标变换后得到的实数信号；和/或
采用下式对x偏振态复数信号的虚部信号进行坐标变换：
x′q(k)＝xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
其中，xq(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中x偏振态复数信号的虚部信号，x′q(k)为对虚部信号xq(k)进行坐标变换后得到的实数信号；和/或
采用下式对y偏振态复数信号的实部信号进行坐标变换：
y′i(k)＝yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
其中，yi(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中y偏振态复数信号的实部信号，y′i(k)为对实部信号yi(k)进行坐标变换后得到的实数信号；和/或
采用下式对y偏振态复数信号的虚部信号进行坐标变换：
y′q(k)＝yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
其中，yq(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中y偏振态复数信号的虚部信号，y′q(k)为对虚部信号yq(k)进行坐标变换后得到的实数信号。
可选的，采用下式计算x偏振态复数信号的实部信号的误差：
εxi(k)＝x′i(k)(|Rx′i|2-|x′i(k)|2)
其中，Rx′i为实部信号xi(k)的参考模值，εxi(k)为实部信号xi(k)的误差；和/或
采用下式计算x偏振态复数信号的虚部信号的误差：
εxq(k)＝x′q(k)(|Rx′q|2-|x′q(k)|2)
其中，Rx′q为虚部信号xq(k)的参考模值，εxq(k)为虚部信号xq(k)的误差；和/或
采用下式计算y偏振态复数信号的实部信号的误差：
εyi(k)＝y′i(k)(|Ry′i|2-|y′i(k)|2)
其中，Ry′i为实部信号yi(k)的参考模值，εyi(k)为实部信号yi(k)的误差；和/或
采用下式计算y偏振态复数信号虚部信号的误差：
εyq(k)＝y′q(k)(|Ry′q|2-|y′q(k)|2)
其中，Ry′q为虚部信号yq(k)的参考模值，εyq(k)为虚部信号yq(k)的误差。
可选的，采用下式计算x偏振态的复数误差：
εx(k)＝εxi(k)+jεxq(k)
其中，εx(k)为x偏振态的复数误差；和/或
采用下式计算y偏振态的复数误差：
εy(k)＝εyi(k)+jεyq(k)
其中，εy(k)为y偏振态的复数误差。
可选的，采用下式计算第k+1时刻的抽头系数：
Hk+1,xx(m)＝Hk,xx(m)+μ·εx(k)·xin(k-m)*
Hk+1,xy(m)＝Hk,xy(m)+μ·εx(k)·yin(k-m)*
Hk+1,yx(m)＝Hk,yx(m)+μ·εy(k)·xin(k-m)*
Hk+1,yy(m)＝Hk,yy(m)+μ·εy(k)·yin(k-m)*
其中，μ为抽头更新的步长，Hk,xx(m)、Hk,xy(m)、Hk,yx(m)、Hk,yy(m)为在第k时刻第m个抽头的系数，Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yy(m)为在第k+1时刻第m个抽头的系数，xin(k-m)*为在第k时刻x偏振态的输入信号的共轭，yin(k-m)*为在第k时刻y偏振态的输入信号的共轭。
可选的，采用下式对第k+1时刻的输入信号进行偏振解复用处理：
其中，xin(k+1-m)为在第k+1时刻x偏振态的输入信号，yin(k+1-m)为在第k+1时刻y偏振态的输入信号，x(k+1)为对第k+1时刻x偏振态输入信号进行偏振解复用处理后得到的信号，y(k+1)为对第k+1时刻y偏振态输入信号进行偏振解复用处理后得到的信号，2N+1为抽头个数。
可选的，采用下式计算实部信号xi(k)的参考模值Rx′i：
其中，sxi(k)为第k时刻x偏振态理想星座下复数信号的实部；和/或
采用下式计算虚部信号xq(k)的参考模值Rx′q：
其中，sxq(k)为第k时刻x偏振态理想星座下复数信号的虚部；和/或
采用下式计算实部信号yi(k)的参考模值Ry′i：
其中，syi(k)为第k时刻y偏振态理想星座下复数信号的实部；和/或
采用下式计算虚部信号yq(k)的参考模值Ry′q：
其中，syq(k)为第k时刻y偏振态理想星座下复数信号的虚部。
本发明提供的自适应盲偏振解复用处理装置包括：
抽头系数确定模块，用于确定自适应FIR滤波器在当前时刻的抽头系数；
所述自适应FIR滤波器，用于根据所述当前时刻的抽头系数，对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；
其中，所述抽头系数确定模块具体用于：
对所述自适应蝶形FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；
分别对该偏振态复数信号的实部信号和虚部信号进行坐标变换；
根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；
将该实部信号的误差和该虚部信号的误差进行合成，得到该偏振态复数信号的复数误差；
根据所述上一时刻的抽头系数和各个偏振态复数信号，计算所述当前时刻的抽头系数。
本发明提供的自适应盲偏振解复用处理方法，在计算下一时刻的抽头系数时需要对该下一时刻的上一时刻的输出信号中每一偏振态复数信号进行分解为实部信号和虚部信号，对分解后的信号进行坐标变换、求误差，误差合并得到复数误差，再根据复数误差计算该下一时刻的抽头系数，自适应FIR滤波器接收到该下一时刻的输入信号时便可以利用该下一时刻的抽头系数进行偏振解复用的处理工作。由于本发明通过分解、坐标变换可以使不同调制格式的多模信号变换到同一个圆上成为恒模信号，因此，可以在不知道信号调制格式的前提下对其进行偏振解复用，因此本发明提供的方法与调制格式无关，适用于网络中调制格式愈加灵活，流量愈加动态、多样和不可预测的应用场景。
附图说明
为了更清楚地说明本公开实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些图获得其他的附图。
图1示出了本发明一实施例中自适应盲偏振解复用处理方法的流程示意图。
具体实施方式
下面将结合本公开实施例中的附图，对本公开实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。
本发明提供一种自适应盲偏振解复用处理方法，如图1所示，该方法包括：
确定自适应FIR滤波器在当前时刻的抽头系数；
根据所述当前时刻的抽头系数，所述自适应FIR滤波器对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；
其中，所述确定自适应蝶形FIR滤波器在当前时刻的抽头系数，包括：
对所述自适应蝶形FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；
分别对该偏振态复数信号的实部信号和虚部信号进行坐标变换；
根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；
将该实部信号的误差和该虚部信号的误差进行合成，得到该偏振态复数信号的复数误差；
根据所述上一时刻的抽头系数和各个偏振态复数信号，计算所述当前时刻的抽头系数。
可理解的是，自适应FIR滤波器为自适应的有限长单位冲激响应滤波器。
可理解的是，自适应FIR滤波器对当前时刻的输入信号进行偏振解复用时采用的当前抽头系数是对上一时刻自适应FIR滤波器的输出信号进行一定的处理得到的，例如，对第2个时刻自适应FIR滤波器的输入信号进行偏振解复用时采用的抽头系数是根据第1时刻自适应FIR滤波器的输出信号进行处理得到的，对第3个时刻自适应FIR滤波器的输入信号进行偏振解复用时采用的抽头系数是根据第2时刻自适应FIR滤波器的输出信号进行处理得到的，依次类推……。其中，在第1时刻自适应FIR滤波器的输入信号进行偏振解复用时采用的初始抽头系数是自行设置的。可见，后面每一时刻的抽头系数是根据前一时刻的输出信号得到的。通过这种方式对抽头系数不断的更新，使其满足偏振解复用的要求。
本发明提供的自适应盲偏振解复用处理方法，在计算下一时刻的抽头系数时需要对该下一时刻的上一时刻的输出信号中每一偏振态复数信号进行分解为实部信号和虚部信号，对分解后的信号进行坐标变换、求误差，误差合并得到复数误差，再根据复数误差计算该下一时刻的抽头系数，自适应FIR滤波器接收到该下一时刻的输入信号时便可以利用该下一时刻的抽头系数进行偏振解复用的处理工作。由于本发明将复数信号分解为实部和虚部，然后进行坐标变换，因此可以使不同调制格式的信号(例如，PM-BPSK、PM-QPSK、PM-4PAM、PM-8PAM、PM-16QAM、PM-32QAM、PM-64QAM)变换到同一个圆上，因此，可以在不知道信号调制格式的前提下对其进行偏振解复用，因此本发明提供的方法与调制格式无关，适用于网络流量愈加动态、多样和不可预测的应用场景。而且，本发明提供的方法相对于传统的RDA，具有简单、抗噪能力强、鲁棒性强的优点，相对于传统的DD-LMS，具有容易实现、无需做波载相位恢复的优点。因此本发明提供的方法不仅适用于相干光通信***，对整个光通信***或其他特定应用场景也适用。
可理解的是，自适应FIR滤波器在每一时刻的输出信号中可以包括x偏振态复数信号和y偏振态复数信号。
在具体实施时，可以采用下式对x偏振态复数信号的实部信号进行坐标变换：
x′i(k)＝xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]} (1)
其中，xi(k)为所述自适应FIR滤波器在k时刻的输出信号中x偏振态复数信号的实部信号，x′i(k)为对实部信号xi(k)坐标变换后得到的实数信号。
同样的，可以采用下式对x偏振态复数信号的虚部信号进行坐标变换：
x′q(k)＝xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]} (2)
其中，xq(k)为所述自适应FIR滤波器在k时刻的输出信号中x偏振态复数信号的虚部信号，x′q(k)为对虚部信号xq(k)进行坐标变换后得到的实数信号。
同样的，可以采用下式对y偏振态复数信号的实部信号进行坐标变换：
y′i(k)＝yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]} (3)
其中，yi(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的实部信号，y′i(k)为对实部信号yi(k)进行坐标变换后得到的实数信号。
同样的，可以采用下式对y偏振态复数信号的虚部信号进行坐标变换：
y′q(k)＝yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]} (4)
其中，yq(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的虚部信号，y′q(k)为对虚部信号yq(k)进行坐标变换后得到的实数信号。
可以理解的是，sign()为取符号函数。
在具体实施时，可以采用下式计算x偏振态复数信号的实部信号的误差：
εxi(k)＝x′i(k)(|Rx′i|2-|x′i(k)|2) (5)
其中，Rx′i为实部信号xi(k)的参考模值，εxi(k)为实部信号xi(k)的误差。
同样的，可以采用下式计算x偏振态复数信号的虚部信号的误差：
εxq(k)＝x′q(k)(|Rx′q|2-|x′q(k)|2) (6)
其中，Rx′q为虚部信号xq(k)的参考模值，εxq(k)为虚部信号xq(k)的误差。
同样的，可以采用下式计算y偏振态复数信号的实部信号的误差：
εyi(k)＝y′i(k)(|Ry′i|2-|y′i(k)|2) (7)
其中，Ry′i为实部信号yi(k)的参考模值，εyi(k)为实部信号yi(k)的误差。
同样的，可以采用下式计算y偏振态复数信号虚部信号的误差：
εyq(k)＝y′q(k)(|Ry′q|2-|y′q(k)|2) (8)
其中，Ry′q为虚部信号yq(k)的参考模值，εyq(k)为虚部信号yq(k)的误差。
上式中的Rx′i、Rx′q、Ry′i、Ry′q的计算方法如下：
其中，可以采用下式计算实部信号xi(k)的参考模值Rx′i：
式中，sxi(k)为第k时刻x偏振态理想星座下复数信号的实部；和/或
其中，采用下式计算虚部信号xq(k)的参考模值Rx′q：
式中，sxq(k)为第k时刻x偏振态理想星座下复数信号的虚部；和/或
其中，采用下式计算实部信号yi(k)的参考模值Ry′i：
式中，syi(k)为第k时刻y偏振态理想星座下复数信号的实部；和/或
其中，采用下式计算虚部信号yq(k)的参考模值Ry′q：
式中，syq(k)为第k时刻y偏振态理想星座下复数信号的虚部。
在具体实施时，可以采用下式计算x偏振态的复数误差：
εx(k)＝εxi(k)+jεxq(k) (13)
其中，εx(k)为x偏振态的复数误差。
同样的，可以采用下式计算y偏振态的复数误差：
εy(k)＝εyi(k)+jεyq(k) (14)
其中，εy(k)为y偏振态的复数误差。
在具体实施时，可以采用下式计算第k+1时刻的抽头系数：
Hk+1,xx(m)＝Hk,xx(m)+μ·εx(k)·xin(k-m)* (15)
Hk+1,xy(m)＝Hk,xy(m)+μ·εx(k)·yin(k-m)* (16)
Hk+1,yx(m)＝Hk,yx(m)+μ·εy(k)·xin(k-m)* (17)
Hk+1,yy(m)＝Hk,yy(m)+μ·εy(k)·yin(k-m)* (18)
其中，μ为抽头更新的步长，Hk,xx(m)、Hk,xy(m)、Hk,yx(m)、Hk,yy(m)为在第k时刻第m个抽头的系数，Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yy(m)为在第k+1时刻第m个抽头的系数，xin(k-m)*为在第k时刻x偏振态的输入信号的共轭，yin(k-m)*为在第k时刻y偏振态的输入信号的共轭。
在具体实施时，可以采用下式对第k+1时刻的输入信号进行偏振解复用处理：
其中，xin(k+1-m)为在第k+1时刻x偏振态的输入信号，yin(k+1-m)为在第k+1时刻y偏振态的输入信号，x(k+1)为对第k+1时刻x偏振态输入信号进行偏振解复用处理后得到的信号，y(k+1)为对第k+1时刻y偏振态输入信号进行偏振解复用处理后得到的信号，2N+1为抽头个数。
可理解的是，m的取值范围为[-N，N]。
在具体实施时，还可以在将信号输入至滤波器之前对信号进行正交不平衡补偿、固定色散补偿和时钟同步等预处理，然后将预处理后的信号输入至滤波器进行处理，并输出，当然对于输出信号还可以进行频偏估计与补偿、相偏估计与补偿、符号反映射等处理，最终得到发端原始比特序列。
通过实验证明，采用本发明提供的方法对信号进行偏解振处理，相对于传统的CMA、RDE等方法，都具有性能优越的特点。
基于相同的发明构思，本发明还提供一种自适应盲偏振解复用处理装置，该装置包括：
抽头系数确定模块，用于确定自适应FIR滤波器在当前时刻的抽头系数；
所述自适应FIR滤波器，用于根据所述当前时刻的抽头系数，对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；
其中，所述抽头系数确定模块具体用于：
对所述自适应FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；
分别对该偏振态复数信号的实部信号和虚部信号进行坐标变换；
根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；
将该实部信号的误差和该虚部信号的误差进行合成，得到该偏振态复数信号的复数误差；
根据所述上一时刻的抽头系数和各个偏振态复数信号，计算所述当前时刻的抽头系数。
本领域普通技术人员可以理解：实现上述方法实施例的全部或者部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储在计算机可读取的存储介质中，该程序在执行时，执行包括上述方法实施例的步骤。
本发明的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解；其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种自适应盲偏振解复用处理方法，其特征在于，包括：
确定自适应蝶形FIR滤波器在当前时刻的抽头系数；
根据所述当前时刻的抽头系数，所述自适应蝶形FIR滤波器对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；
其中，所述确定自适应蝶形FIR滤波器在当前时刻的抽头系数，包括：
对所述自适应蝶形FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；
分别对该偏振态复数信号的实部信号和虚部信号进行坐标变换；
根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；
将该实部信号的误差和该虚部信号的误差进行合成，得到该偏振态复数信号的复数误差；
根据所述上一时刻的抽头系数和各个偏振态复数信号，计算所述当前时刻的抽头系数；
所述自适应FIR滤波器在每一时刻的输出信号中包括x偏振态复数信号和y偏振态复数信号；对应的，
采用下式对x偏振态复数信号的实部信号进行坐标变换：
x′i(k)＝xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
其中，xi(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中x偏振态复数信号的实部信号，x′i(k)为对实部信号xi(k)坐标变换后得到的实数信号；和/或
采用下式对x偏振态复数信号的虚部信号进行坐标变换：
x′q(k)＝xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
其中，xq(k)为所述自适应FIR滤波器在k时刻的输出信号中x偏振态复数信号的虚部信号，x′q(k)为对虚部信号xq(k)进行坐标变换后得到的实数信号；和/或
采用下式对y偏振态复数信号的实部信号进行坐标变换：
y′i(k)＝yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
其中，yi(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的实部信号，y′i(k)为对实部信号yi(k)进行坐标变换后得到的实数信号；和/或
采用下式对y偏振态复数信号的虚部信号进行坐标变换：
y′q(k)＝yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
其中，yq(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的虚部信号，y′q(k)为对虚部信号yq(k)进行坐标变换后得到的实数信号。
2.根据权利要求1所述的方法，其特征在于，
采用下式计算x偏振态复数信号的实部信号的误差：
εxi(k)＝x′i(k)(|Rx′i|2-|x′i(k)|2)
其中，Rx′i为实部信号xi(k)的参考模值，εxi(k)为实部信号xi(k)的误差；和/或
采用下式计算x偏振态复数信号的虚部信号的误差：
εxq(k)＝x′q(k)(|Rx′q|2-|x′q(k)|2)
其中，Rx′q为虚部信号xq(k)的参考模值，εxq(k)为虚部信号xq(k)的误差；和/或
采用下式计算y偏振态复数信号的实部信号的误差：
εyi(k)＝y′i(k)(|Ry′i|2-|y′i(k)|2)
其中，Ry′i为实部信号yi(k)的参考模值，εyi(k)为实部信号yi(k)的误差；和/或
采用下式计算y偏振态复数信号虚部信号的误差：
εyq(k)＝y′q(k)(|Ry′q|2-|y′q(k)|2)
其中，Ry′q为虚部信号yq(k)的参考模值，εyq(k)为虚部信号yq(k)的误差。
3.根据权利要求2所述的方法，其特征在于，
采用下式计算x偏振态的复数误差：
εx(k)＝εxi(k)+jεxq(k)
其中，εx(k)为x偏振态的复数误差；和/或
采用下式计算y偏振态的复数误差：
εy(k)＝εyi(k)+jεyq(k)
其中，εy(k)为y偏振态的复数误差。
4.根据权利要求3所述的方法，其特征在于，采用下式计算第k+1时刻的抽头系数：
Hk+1,xx(m)＝Hk,xx(m)+μ·εx(k)·xin(k-m)*
Hk+1,xy(m)＝Hk,xy(m)+μ·εx(k)·yin(k-m)*
Hk+1,yx(m)＝Hk,yx(m)+μ·εy(k)·xin(k-m)*
Hk+1,yy(m)＝Hk,yy(m)+μ·εy(k)·yin(k-m)*
其中，μ为抽头更新的步长，Hk,xx(m)、Hk,xy(m)、Hk,yx(m)、Hk,yy(m)为在第k时刻第m个抽头的系数，Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yy(m)为在第k+1时刻第m个抽头的系数，xin(k-m)*为在第k时刻x偏振态的输入信号的共轭，yin(k-m)*为在第k时刻y偏振态的输入信号的共轭。
5.根据权利要求4所述的方法，其特征在于，采用下式对第k+1时刻的输入信号进行偏振解复用处理：
其中，xin(k+1-m)为在第k+1时刻x偏振态的输入信号，yin(k+1-m)为在第k+1时刻y偏振态的输入信号，x(k+1)为对第k+1时刻两偏振态输入信号进行偏振解复用处理后得到的x偏振态信号，y(k+1)为对第k+1时刻两偏振态输入信号进行偏振解复用处理后得到的y偏振态信号，2N+1为抽头个数。
6.根据权利要求2所述的方法，其特征在于，
采用下式计算实部信号xi(k)的参考模值Rx′i：
其中，sxi(k)为第k时刻x偏振态理想星座下复数信号的实部；和/或
采用下式计算虚部信号xq(k)的参考模值Rx′q：
其中，sxq(k)为第k时刻x偏振态理想星座下复数信号的虚部；和/或
采用下式计算实部信号yi(k)的参考模值Ry′i：
其中，syi(k)为第k时刻y偏振态理想星座下复数信号的实部；和/或
采用下式计算虚部信号yq(k)的参考模值Ry′q：
其中，syq(k)为第k时刻y偏振态理想星座下复数信号的虚部。
7.一种自适应盲偏振解复用处理装置，其特征在于，包括：
抽头系数确定模块，用于确定自适应FIR滤波器在当前时刻的抽头系数；
所述自适应FIR滤波器，用于根据所述当前时刻的抽头系数，对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；
其中，所述抽头系数确定模块具体用于：
对所述自适应FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；
分别对该偏振态复数信号的实部信号和虚部信号进行坐标变换；
根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；
将该实部信号的误差和该虚部信号的误差进行合成，得到该偏振态复数信号的复数误差；
根据所述上一时刻的抽头系数和各个偏振态复数信号，计算所述当前时刻的抽头系数；
所述自适应FIR滤波器在每一时刻的输出信号中包括x偏振态复数信号和y偏振态复数信号；对应的，
采用下式对x偏振态复数信号的实部信号进行坐标变换：
x′i(k)＝xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
其中，xi(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中x偏振态复数信号的实部信号，x′i(k)为对实部信号xi(k)坐标变换后得到的实数信号；和/或
采用下式对x偏振态复数信号的虚部信号进行坐标变换：
x′q(k)＝xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
其中，xq(k)为所述自适应FIR滤波器在k时刻的输出信号中x偏振态复数信号的虚部信号，x′q(k)为对虚部信号xq(k)进行坐标变换后得到的实数信号；和/或
采用下式对y偏振态复数信号的实部信号进行坐标变换：
y′i(k)＝yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
其中，yi(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的实部信号，y′i(k)为对实部信号yi(k)进行坐标变换后得到的实数信号；和/或
采用下式对y偏振态复数信号的虚部信号进行坐标变换：
y′q(k)＝yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
其中，yq(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的虚部信号，y′q(k)为对虚部信号yq(k)进行坐标变换后得到的实数信号。
CN201610650339.5A 2016-08-09 2016-08-09 自适应盲偏振解复用处理方法和装置 Active CN106301593B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610650339.5A CN106301593B (zh) 2016-08-09 2016-08-09 自适应盲偏振解复用处理方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610650339.5A CN106301593B (zh) 2016-08-09 2016-08-09 自适应盲偏振解复用处理方法和装置

Publications (2)

Publication Number Publication Date
CN106301593A CN106301593A (zh) 2017-01-04
CN106301593B true CN106301593B (zh) 2018-11-16

Family

ID=57667631

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610650339.5A Active CN106301593B (zh) 2016-08-09 2016-08-09 自适应盲偏振解复用处理方法和装置

Country Status (1)

Country Link
CN (1) CN106301593B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110011733B (zh) * 2019-03-25 2020-10-16 华中科技大学 一种基于动量因子的解偏振复用方法及***
CN111884960B (zh) * 2020-09-28 2020-12-25 烽火通信科技股份有限公司 一种偏振解复用信号处理盲均衡方法及偏振解复用装置
CN113708845B (zh) * 2021-09-08 2022-05-06 北京科技大学 均衡解复用方法、装置、电子设备及计算机可读存储介质
CN114337843B (zh) * 2021-12-23 2023-11-28 武汉邮电科学研究院有限公司 一种用于相干光通信的自适应信道均衡方法及***

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103283201A (zh) * 2011-03-31 2013-09-04 富士通株式会社 蝶形滤波器系数设定方法和装置、接收机及接收方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101895499B (zh) * 2010-06-30 2014-11-19 北京邮电大学 用恒定转动符号训练序列解偏振复用的方法
US9154231B2 (en) * 2013-01-17 2015-10-06 Alcatel Lucent Generation of an optical local-oscillator signal for a coherent-detection scheme
CN105703838B (zh) * 2016-01-26 2018-07-03 哈尔滨工业大学深圳研究生院 基于蝶形线性卡尔曼滤波器的相干光接收机动态均衡方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103283201A (zh) * 2011-03-31 2013-09-04 富士通株式会社 蝶形滤波器系数设定方法和装置、接收机及接收方法

Also Published As

Publication number Publication date
CN106301593A (zh) 2017-01-04

Similar Documents

Publication Publication Date Title
CN106301593B (zh) 自适应盲偏振解复用处理方法和装置
US8260156B2 (en) Adaptive crossing frequency domain equalization (FDE) in digital PolMux coherent systems
US9467246B2 (en) Clock recovery method and clock recovery arrangement for coherent polarization multiplex receivers
US9515742B2 (en) System and methods for adaptive equalization for optical modulation formats
CN103338171B (zh) 一种基于频域信道估计的接收端均衡方法和***
US20110194855A1 (en) Superimposed training and digital filtering coherent optical receivers
JP5120507B2 (ja) 光受信器、偏光分離装置および偏光分離方法
EP2436127A1 (en) Method for carrier frequency recovery and optical intradyne coherent receiver
CN104579621B (zh) 时钟恢复方法和装置
US10171177B2 (en) Digital signal processor, digital optical receiver using the same, and digital signal processing method
CN104168069B (zh) 相干光通信领域中相位受控旋转的调制***及调制方法
CN105530050A (zh) 均衡与偏振解复用和相偏估计与补偿的联合处理方法及装置
CN112929124B (zh) 一种实现偏振复用信号的解偏与均衡的方法及***
CN107925485A (zh) 相干光接收装置
CN102724151B (zh) 通信***、通信方法以及多通道自适应均衡器和ofdm信号解复用的方法
CN106789793A (zh) 用于自适应均衡器系数更新的无需频差相差反馈的lms方法
CN108199777A (zh) 基于概率密度函数拟合和模糊逻辑的相干光通信***盲均衡方法
JP6378654B2 (ja) 空間多重光伝送システム及び空間多重光伝送方法
CN105553529B (zh) 一种少模光纤传输***及其数字信号恢复方法
CN110224761A (zh) 一种快速求解偏振旋转矩阵和均衡器系数的方法及***
CN105282072B (zh) 光传输网络中的频偏估计方法与装置
JP6387832B2 (ja) 偏光多重分離光通信受信機、偏光多重分離光通信システム、および偏光多重分離光通信方法
CN108173794A (zh) 一种新型lms方法及使用该方法的实数自适应均衡器
JPH0621762A (ja) デイジタル形トランスバーサル自動等化器
EP3076568B1 (en) Method for eliminating nonlinear effect, transmitter and receiver

Legal Events

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