CN106301593A

CN106301593A - Adaptive blind polarization demultiplexing treating method and apparatus

Info

Publication number: CN106301593A
Application number: CN201610650339.5A
Authority: CN
Inventors: 陈雪; 杨桃; 王立芊
Original assignee: Beijing University of Posts and Telecommunications
Current assignee: Beijing University of Posts and Telecommunications
Priority date: 2016-08-09
Filing date: 2016-08-09
Publication date: 2017-01-04
Anticipated expiration: 2036-08-09
Also published as: CN106301593B

Abstract

The present invention relates to an adaptive blind polarization demultiplexing processing method and device, the method comprising: determining the tap coefficient of the adaptive FIR filter at the current moment; according to the tap coefficient at the current moment, the adaptive FIR filter Perform polarization demultiplexing processing on the input signal at the current moment, and output the signal after polarization demultiplexing processing; wherein, according to the tap coefficient at the previous moment and the two polarization state complex signals, calculate the current moment tap coefficient. Since the present invention can convert multi-mode signals of different modulation formats to the same circle to become constant-mode signals through decomposition and coordinate transformation, it can perform polarization demultiplexing on the premise of not knowing the signal modulation format, so the present invention The method provided has nothing to do with the modulation format, and is suitable for application scenarios where the modulation format in the network becomes more flexible and the traffic becomes more dynamic, diverse and unpredictable.

Description

Adaptive blind polarization demultiplexing processing method and device

技术领域technical field

本发明涉及光纤通信技术领域，尤其是涉及一种自适应盲偏振解复用处理方法和装置。The invention relates to the technical field of optical fiber communication, in particular to an adaptive blind polarization demultiplexing processing method and device.

背景技术Background technique

近年来，随着网络流量和带宽需求的急剧增加，高速相干光通信技术已成为实现长距离大容量信息传输的关键技术。通过偏振复用和波分复用等技术，可以有效的提高光纤通信系统的带宽利用率，提升系统容量。对高速相干光通信而言，由于光纤色散、偏振随机串扰等因素的影响，因此需要采用多抽头的自适应蝶形FIR滤波器进行残余色散、偏振膜色散的均衡和偏振解复用。In recent years, with the sharp increase in network traffic and bandwidth requirements, high-speed coherent optical communication technology has become a key technology for long-distance and large-capacity information transmission. Through technologies such as polarization multiplexing and wavelength division multiplexing, the bandwidth utilization rate of the optical fiber communication system can be effectively improved, and the system capacity can be increased. For high-speed coherent optical communication, due to the influence of factors such as fiber dispersion and polarization random crosstalk, it is necessary to use a multi-tap adaptive butterfly FIR filter for residual dispersion, equalization of polarization film dispersion, and polarization demultiplexing.

目前，偏振解复用方法包括数据辅助和盲处理两种。其中，数据辅助方法通过发送训练序列来进行信道估计，占用了部分频谱资源，降低了系统带宽利用率。然而，盲处理方法作为偏分复用光通信解调系统中的一项关键技术，无需周期性地发送训练序列，可仅仅依靠接收信号的统计特性来恢复原始发送信号，因此无需占用太多带宽资源。可见两种方法相比，盲处理方法尤显优势。由于数据中心和云计算的不断发展，网络流量将变得愈加动态、多样和不可预测，支持多种调制格式的弹性光收发机成为弹性光网络中有效承载突发动态变化网络流量、提高频谱利用率、优化网络资源利用的关键。因此有必要提供一种与调制格式无关的盲处理方法。Currently, polarization demultiplexing methods include data-assisted and blind processing. Among them, the data-assisted method performs channel estimation by sending training sequences, which occupies part of spectrum resources and reduces system bandwidth utilization. However, as a key technology in the polarization division multiplexing optical communication demodulation system, the blind processing method does not need to periodically send the training sequence, and can only rely on the statistical characteristics of the received signal to restore the original sent signal, so it does not need to occupy too much bandwidth resource. It can be seen that compared with the two methods, the blind processing method is particularly superior. Due to the continuous development of data centers and cloud computing, network traffic will become more dynamic, diverse and unpredictable. Elastic optical transceivers that support multiple modulation formats have become an effective way to carry sudden dynamic changes in network traffic and improve spectrum utilization in elastic optical networks. The key to efficiency and optimal utilization of network resources. Therefore, it is necessary to provide a blind processing method independent of the modulation format.

发明内容Contents of the invention

针对以上缺陷，本发明提供一种自适应盲偏振解复用处理方法和装置，处理过程与调制格式无关，适用于网络流量愈加动态、多样和不可预测的应用场景。In view of the above defects, the present invention provides an adaptive blind polarization demultiplexing processing method and device, the processing process has nothing to do with the modulation format, and is suitable for application scenarios where network traffic becomes more dynamic, diverse and unpredictable.

本发明提供的自适应盲偏振解复用处理方法包括：The adaptive blind polarization demultiplexing processing method provided by the present invention includes:

确定自适应FIR滤波器在当前时刻的抽头系数；Determine the tap coefficient of the adaptive FIR filter at the current moment;

根据所述当前时刻的抽头系数，所述自适应FIR滤波器对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；According to the tap coefficient at the current moment, the adaptive FIR filter performs polarization demultiplexing processing on the input signal at the current moment, and outputs the signal after polarization demultiplexing processing;

其中，所述确定自适应FIR滤波器在当前时刻的抽头系数，包括：Wherein, said determining the tap coefficient of the adaptive FIR filter at the current moment includes:

对所述自适应FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；Decomposing each polarization state complex signal in the output signal of the adaptive FIR filter at a previous moment, to obtain a real part signal and an imaginary part signal of the polarization state complex signal;

分别对该偏振态复数信号的实部信号和虚部信号进行坐标变换；Coordinate transformation is performed on the real part signal and the imaginary part signal of the polarization state complex signal respectively;

根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；According to the reference modulus of the real part signal of the polarization state complex signal and the real signal obtained after coordinate transformation, calculate the error of the real part signal, and according to the reference modulus and coordinate transformation of the imaginary part signal of the polarization state complex signal Calculate the error of the imaginary part signal of the real number signal;

将该实部信号的误差和该虚部信号的误差进行合成，得到该偏振态复数信号的复数误差；Combining the error of the real part signal with the error of the imaginary part signal to obtain the complex error of the polarization state complex signal;

根据所述上一时刻的抽头系数和各个偏振态复数信号，计算所述当前时刻的抽头系数。The tap coefficient at the current time is calculated according to the tap coefficient at the last time and each polarization state complex signal.

可选的，所述自适应蝶形FIR滤波器在每一时刻的输出信号中包括x偏振态复数信号和y偏振态复数信号；对应的，Optionally, the output signal of the adaptive butterfly FIR filter at each moment includes an x polarization state complex signal and a y polarization state complex signal; correspondingly,

采用下式对x偏振态复数信号的实部信号进行坐标变换：Use the following formula to transform the coordinates of the real part of the x-polarized complex signal:

x′_i(k)＝x_i(k)-4·sign[x_i(k)]-2·sign{x_i(k)-4·sign[x_i(k)]}x′ _i (k)= _xi (k)-4·sign[ _xi (k)]-2·sign{ _xi (k)-4·sign[ _xi (k)]}

其中，x_i(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中x偏振态复数信号的实部信号，x′_i(k)为对实部信号x_i(k)坐标变换后得到的实数信号；和/或Wherein, x _i (k) is the real part signal of the x polarization state complex signal in the output signal of the adaptive butterfly FIR filter at time k, and x' _i (k) is the real part signal of the real part signal x _i (k) the real signal obtained after coordinate transformation; and/or

采用下式对x偏振态复数信号的虚部信号进行坐标变换：Use the following formula to transform the imaginary part signal of the x-polarization state complex signal:

x′_q(k)＝x_q(k)-4·sign[x_q(k)]-2·sign{x_q(k)-4·sign[x_q(k)]}x′ _q (k)＝x _q (k)-4·sign[x _q (k)]-2·sign{x _q (k)-4·sign[x _q (k)]}

其中，x_q(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中x偏振态复数信号的虚部信号，x′_q(k)为对虚部信号x_q(k)进行坐标变换后得到的实数信号；和/或Wherein, x _q (k) is the imaginary part signal of the x polarization state complex signal in the output signal of the adaptive butterfly FIR filter at time k, and x′ _q (k) is the imaginary part signal x _q (k) the real signal obtained after coordinate transformation; and/or

采用下式对y偏振态复数信号的实部信号进行坐标变换：Use the following formula to transform the real part of the y-polarized complex signal:

y′_i(k)＝y_i(k)-4·sign[y_i(k)]-2·sign{y_i(k)-4·sign[y_i(k)]}y′ _i (k)=y _i (k)-4·sign[y _i (k)]-2·sign{y _i (k)-4·sign[y _i (k)]}

其中，y_i(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中y偏振态复数信号的实部信号，y′_i(k)为对实部信号y_i(k)进行坐标变换后得到的实数信号；和/或Wherein, y _i (k) is the real part signal of the y polarization state complex signal in the output signal of the adaptive butterfly FIR filter at k moment, and y' _i (k) is the real part signal y _i (k) of the real part signal the real signal obtained after coordinate transformation; and/or

采用下式对y偏振态复数信号的虚部信号进行坐标变换：Use the following formula to transform the imaginary part signal of the y-polarized complex signal:

y′_q(k)＝y_q(k)-4·sign[y_q(k)]-2·sign{y_q(k)-4·sign[y_q(k)]}y′ _q (k)=y _q (k)-4·sign[y _q (k)]-2·sign{y _q (k)-4·sign[y _q (k)]}

其中，y_q(k)为所述自适应蝶形FIR滤波器在k时刻的输出信号中y偏振态复数信号的虚部信号，y′_q(k)为对虚部信号y_q(k)进行坐标变换后得到的实数信号。Wherein, y _q (k) is the imaginary part signal of the y polarization state complex signal in the output signal of the adaptive butterfly FIR filter at k moment, and y′ _q (k) is the imaginary part signal y _q (k) Real signal obtained after coordinate transformation.

可选的，采用下式计算x偏振态复数信号的实部信号的误差：Optionally, the following formula is used to calculate the error of the real part signal of the x polarization state complex signal:

ε_xi(k)＝x′_i(k)(|R_x′i|²-|x′_i(k)|²)ε _xi (k)=x′ _i (k)(|R _x′i | ² -|x′ _i (k)| ² )

其中，R_x′i为实部信号x_i(k)的参考模值，ε_xi(k)为实部信号x_i(k)的误差；和/或Wherein, R _x'i is the reference modulus value of the real part signal x _i (k), ε _xi (k) is the error of the real part signal x _i (k); and/or

采用下式计算x偏振态复数信号的虚部信号的误差：Use the following formula to calculate the error of the imaginary part signal of the x polarization state complex signal:

ε_xq(k)＝x′_q(k)(|R_x′q|²-|x′_q(k)|²)ε _xq (k)=x′ _q (k)(|R _x′q | ² -|x′ _q (k)| ² )

其中，R_x′q为虚部信号x_q(k)的参考模值，ε_xq(k)为虚部信号x_q(k)的误差；和/或Wherein, R _x'q is the reference modulus value of the imaginary part signal x _q (k), ε _xq (k) is the error of the imaginary part signal x _q (k); and/or

采用下式计算y偏振态复数信号的实部信号的误差：The error of the real part signal of the y-polarized complex signal is calculated using the following formula:

ε_yi(k)＝y′_i(k)(|R_y′i|²-|y′_i(k)|²)ε _yi (k)=y′ _i (k)(|R _y′i | ² -|y′ _i (k)| ² )

其中，R_y′i为实部信号y_i(k)的参考模值，ε_yi(k)为实部信号y_i(k)的误差；和/或Wherein, R _y'i is the reference modulus value of the real part signal y _i (k), ε _yi (k) is the error of the real part signal y _i (k); and/or

采用下式计算y偏振态复数信号虚部信号的误差：Use the following formula to calculate the error of the imaginary part signal of the y-polarized complex signal:

ε_yq(k)＝y′_q(k)(|R_y′q|²-|y′_q(k)|²)ε _yq (k)=y′ _q (k)(|R _y′q | ² -|y′ _q (k)| ² )

其中，R_y′q为虚部信号y_q(k)的参考模值，ε_yq(k)为虚部信号y_q(k)的误差。Among them, R _y′q is the reference modulus value of the imaginary part signal y _q (k), and ε _yq (k) is the error of the imaginary part signal y _q (k).

可选的，采用下式计算x偏振态的复数误差：Optionally, the complex error of the x polarization state is calculated using the following formula:

ε_x(k)＝ε_xi(k)+jε_xq(k)ε _x (k) = ε _xi (k) + jε _xq (k)

其中，ε_x(k)为x偏振态的复数误差；和/或where ε _x (k) is the complex error of the x polarization state; and/or

采用下式计算y偏振态的复数误差：The complex error for the y-polarization state is calculated using the following equation:

ε_y(k)＝ε_yi(k)+jε_yq(k)ε _y (k) = ε _yi (k) + jε _yq (k)

其中，ε_y(k)为y偏振态的复数误差。where ε _y (k) is the complex error of the y polarization state.

可选的，采用下式计算第k+1时刻的抽头系数：Optionally, the tap coefficient at the k+1th moment is calculated using the following formula:

H_k+1,xx(m)＝H_k,xx(m)+μ·ε_x(k)·x_in(k-m)^* H _k+1,xx (m)＝H _k,xx (m)+μ·ε _x (k)·x _in (km) ^*

H_k+1,xy(m)＝H_k,xy(m)+μ·ε_x(k)·y_in(k-m)^* H _k+1,xy (m)＝H _k,xy (m)+μ·ε _x (k)·y _in (km) ^*

H_k+1,yx(m)＝H_k,yx(m)+μ·ε_y(k)·x_in(k-m)^* H _k+1,yx (m)＝H _k,yx (m)+μ·ε _y (k)·x _in (km) ^*

H_k+1,yy(m)＝H_k,yy(m)+μ·ε_y(k)·y_in(k-m)^* H _k+1,yy (m)＝H _k,yy (m)+μ·ε _y (k)·y _in (km) ^*

其中，μ为抽头更新的步长，H_k,xx(m)、H_k,xy(m)、H_k,yx(m)、H_k,yy(m)为在第k时刻第m个抽头的系数，H_k+1,xx(m)、H_k+1,xy(m)、H_k+1,yx(m)、H_k+1,yy(m)为在第k+1时刻第m个抽头的系数，x_in(k-m)^*为在第k时刻x偏振态的输入信号的共轭，y_in(k-m)^*为在第k时刻y偏振态的输入信号的共轭。Among them, μ is the step size of the tap update, H _k,xx (m), H _k,xy (m), H _k,yx (m), H _k,yy (m) is the mth tap at the kth moment The coefficients of H _k+1,xx (m), H _k+1,xy (m), H _k+1,yx (m), H _k+1,yy (m) are The coefficients of m taps, x _in (km) ^* is the conjugate of the input signal of the x polarization state at the kth moment, and y _in (km) ^* is the conjugate of the input signal of the y polarization state at the kth moment.

可选的，采用下式对第k+1时刻的输入信号进行偏振解复用处理：Optionally, the following formula is used to perform polarization demultiplexing processing on the input signal at the k+1th moment:

$[\begin{matrix} x x ((k k + + 11)) \\ y the y ((k k + + 11)) \end{matrix}] = = [\begin{matrix} {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, x x x x} ((m m)) \cdot \cdot {x x}_{i i n no} ((k k + + 11 - - m m)) + + {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, x x y the y} ((m m)) \cdot \cdot {y the y}_{i i n no} ((k k + + 11 - - m m)) \\ {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, y the y x x} ((m m)) \cdot \cdot {x x}_{i i n no} ((k k + + 11 - - m m)) + + {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, y the y y the y} ((m m)) \cdot \cdot {y the y}_{i i n no} ((k k + + 11 - - m m)) \end{matrix}]$

其中，x_in(k+1-m)为在第k+1时刻x偏振态的输入信号，y_in(k+1-m)为在第k+1时刻y偏振态的输入信号，x(k+1)为对第k+1时刻x偏振态输入信号进行偏振解复用处理后得到的信号，y(k+1)为对第k+1时刻y偏振态输入信号进行偏振解复用处理后得到的信号，2N+1为抽头个数。Among them, x _in (k+1-m) is the input signal of the x polarization state at the k+1 time, y _in (k+1-m) is the input signal of the y polarization state at the k+1 time, x( k+1) is the signal obtained by performing polarization demultiplexing on the input signal of the x polarization state at the k+1 time, and y(k+1) is the polarization demultiplexing of the input signal of the y polarization state at the k+1 time For the signal obtained after processing, 2N+1 is the number of taps.

可选的，采用下式计算实部信号x_i(k)的参考模值R_x′i：Optionally, the following formula is used to calculate the reference modulus R _x'i of the real part signal x _i (k):

${R R}_{{x x}^{' '} i i} = = \sqrt{\frac{E E. {{{| | {s the s}_{x x i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x i i} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{x x i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x i i} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{x x i i} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{x x i i} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{x x i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x i i} ((k k))]]}} | |}^{22}}}}}$

其中，s_xi(k)为第k时刻x偏振态理想星座下复数信号的实部；和/或Wherein, s _xi (k) is the real part of the complex signal under the ideal constellation of polarization state x at the kth moment; and/or

采用下式计算虚部信号x_q(k)的参考模值R_x′q：Use the following formula to calculate the reference modulus R _x′q of the imaginary part signal x _q (k):

${R R}_{{x x}^{' '} q q} = = \sqrt{\frac{E E. {{{| | {s the s}_{x x q q} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x q q} ((k k))]] - - 22 \cdot \cdot s the s i i g g n no {{{x x}_{x x q q} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x q q} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{x x q q} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x q q} ((k k))]] - - 22 \cdot \cdot s the s i i g g n no {{{x x}_{x x q q} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x q q} ((k k))]]}} | |}^{22}}}}}$

其中，s_xq(k)为第k时刻x偏振态理想星座下复数信号的虚部；和/或Wherein, s _xq (k) is the imaginary part of the complex signal under the ideal constellation of polarization state x at the kth moment; and/or

采用下式计算实部信号y_i(k)的参考模值R_y′i：Use the following formula to calculate the reference modulus R _y′i of the real part signal y _i (k):

${R R}_{{y the y}^{' '} i i} = = \sqrt{\frac{E E. {{{| | {s the s}_{y the y i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]] - - 22 \cdot \cdot s the s i i g g n no {{{s the s}_{y the y i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{y the y i i} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{y the y i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]]}} | |}^{22}}}}}$

其中，s_yi(k)为第k时刻y偏振态理想星座下复数信号的实部；和/或Wherein, s _yi (k) is the real part of the complex signal under the ideal constellation of y polarization state at the kth moment; and/or

采用下式计算虚部信号y_q(k)的参考模值R_y′q：Use the following formula to calculate the reference modulus R _y′q of the imaginary part signal y _q (k):

${R R}_{{y the y}^{' '} q q} = = \sqrt{\frac{E E. {{{| | {s the s}_{y the y q q} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]] - - 22 \cdot \cdot s the s i i g g n no {{{s the s}_{y the y q q} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{y the y q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{y the y q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]]}} | |}^{22}}}}}$

其中，s_yq(k)为第k时刻y偏振态理想星座下复数信号的虚部。Wherein, s _yq (k) is the imaginary part of the complex signal under the ideal constellation of the y-polarized state at the kth moment.

本发明提供的自适应盲偏振解复用处理装置包括：The adaptive blind polarization demultiplexing processing device provided by the present invention includes:

抽头系数确定模块，用于确定自适应FIR滤波器在当前时刻的抽头系数；The tap coefficient determination module is used to determine the tap coefficient of the adaptive FIR filter at the current moment;

所述自适应FIR滤波器，用于根据所述当前时刻的抽头系数，对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；The adaptive FIR filter is used to perform polarization demultiplexing processing on the input signal at the current time according to the tap coefficient at the current time, and output the signal after the polarization demultiplexing process;

其中，所述抽头系数确定模块具体用于：Wherein, the tap coefficient determination module is specifically used for:

对所述自适应蝶形FIR滤波器在上一时刻的输出信号中的每一偏振态复数信号进行分解，得到该偏振态复数信号的实部信号和虚部信号；Decomposing each polarization state complex signal in the output signal of the adaptive butterfly FIR filter at the previous moment, to obtain the real part signal and the imaginary part signal of the polarization state complex signal;

根据该偏振态复数信号的实部信号的参考模值和坐标变换后得到的实数信号，计算该实部信号的误差，及根据该偏振态复数信号的虚部信号的参考模值和坐标变换后的实数信号计算该虚部信号的误差；Calculate the error of the real part signal according to the reference modulus of the real part signal of the polarization state complex signal and the real signal obtained after coordinate transformation, and calculate the error of the real part signal according to the reference modulus and coordinate transformation of the imaginary part signal of the polarization state complex signal Calculate the error of the imaginary part signal of the real number signal;

本发明提供的自适应盲偏振解复用处理方法，在计算下一时刻的抽头系数时需要对该下一时刻的上一时刻的输出信号中每一偏振态复数信号进行分解为实部信号和虚部信号，对分解后的信号进行坐标变换、求误差，误差合并得到复数误差，再根据复数误差计算该下一时刻的抽头系数，自适应FIR滤波器接收到该下一时刻的输入信号时便可以利用该下一时刻的抽头系数进行偏振解复用的处理工作。由于本发明通过分解、坐标变换可以使不同调制格式的多模信号变换到同一个圆上成为恒模信号，因此，可以在不知道信号调制格式的前提下对其进行偏振解复用，因此本发明提供的方法与调制格式无关，适用于网络中调制格式愈加灵活，流量愈加动态、多样和不可预测的应用场景。In the adaptive blind polarization demultiplexing processing method provided by the present invention, when calculating the tap coefficient at the next moment, it is necessary to decompose each polarization state complex signal in the output signal at the previous moment at the next moment into real part signals and For the imaginary part signal, coordinate transformation is performed on the decomposed signal, the error is calculated, the error is combined to obtain a complex error, and then the tap coefficient at the next moment is calculated according to the complex error, and when the adaptive FIR filter receives the input signal at the next moment Then, the tap coefficient at the next moment can be used to perform polarization demultiplexing processing. Since the present invention can convert multi-mode signals of different modulation formats to the same circle to become constant-mode signals through decomposition and coordinate transformation, it can perform polarization demultiplexing on the premise of not knowing the signal modulation format, so the present invention The method provided by the invention has nothing to do with the modulation format, and is suitable for application scenarios where the modulation format in the network becomes more flexible and the traffic becomes more dynamic, diverse and unpredictable.

附图说明Description of drawings

为了更清楚地说明本公开实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present disclosure or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only These are some embodiments of the present disclosure. For those skilled in the art, other drawings can also be obtained according to these drawings on the premise of not paying creative efforts.

图1示出了本发明一实施例中自适应盲偏振解复用处理方法的流程示意图。Fig. 1 shows a schematic flowchart of an adaptive blind polarization demultiplexing processing method in an embodiment of the present invention.

具体实施方式detailed description

下面将结合本公开实施例中的附图，对本公开实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present disclosure with reference to the accompanying drawings in the embodiments of the present disclosure. Obviously, the described embodiments are only some of the embodiments of the present invention, not all of them. Based on the embodiments in the present disclosure, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present disclosure.

本发明提供一种自适应盲偏振解复用处理方法，如图1所示，该方法包括：The present invention provides a kind of adaptive blind polarization demultiplexing processing method, as shown in Figure 1, this method comprises:

其中，所述确定自适应蝶形FIR滤波器在当前时刻的抽头系数，包括：Wherein, said determining the tap coefficients of the adaptive butterfly FIR filter at the current moment includes:

可理解的是，自适应FIR滤波器为自适应的有限长单位冲激响应滤波器。It can be understood that the adaptive FIR filter is an adaptive finite-length unit impulse response filter.

可理解的是，自适应FIR滤波器对当前时刻的输入信号进行偏振解复用时采用的当前抽头系数是对上一时刻自适应FIR滤波器的输出信号进行一定的处理得到的，例如，对第2个时刻自适应FIR滤波器的输入信号进行偏振解复用时采用的抽头系数是根据第1时刻自适应FIR滤波器的输出信号进行处理得到的，对第3个时刻自适应FIR滤波器的输入信号进行偏振解复用时采用的抽头系数是根据第2时刻自适应FIR滤波器的输出信号进行处理得到的，依次类推……。其中，在第1时刻自适应FIR滤波器的输入信号进行偏振解复用时采用的初始抽头系数是自行设置的。可见，后面每一时刻的抽头系数是根据前一时刻的输出信号得到的。通过这种方式对抽头系数不断的更新，使其满足偏振解复用的要求。It can be understood that the current tap coefficient used when the adaptive FIR filter performs polarization demultiplexing on the input signal at the current moment is obtained by performing certain processing on the output signal of the adaptive FIR filter at the previous moment, for example, for The tap coefficients used in polarization demultiplexing of the input signal of the adaptive FIR filter at the second moment are obtained by processing the output signal of the adaptive FIR filter at the first moment, and the adaptive FIR filter at the third moment The tap coefficients used when performing polarization demultiplexing of the input signal are obtained by processing the output signal of the adaptive FIR filter at the second moment, and so on. . . . Wherein, the initial tap coefficient used when the input signal of the adaptive FIR filter is subjected to polarization demultiplexing at the first moment is set by itself. It can be seen that the tap coefficient at each subsequent moment is obtained according to the output signal at the previous moment. In this way, the tap coefficients are continuously updated to meet the requirements of polarization demultiplexing.

本发明提供的自适应盲偏振解复用处理方法，在计算下一时刻的抽头系数时需要对该下一时刻的上一时刻的输出信号中每一偏振态复数信号进行分解为实部信号和虚部信号，对分解后的信号进行坐标变换、求误差，误差合并得到复数误差，再根据复数误差计算该下一时刻的抽头系数，自适应FIR滤波器接收到该下一时刻的输入信号时便可以利用该下一时刻的抽头系数进行偏振解复用的处理工作。由于本发明将复数信号分解为实部和虚部，然后进行坐标变换，因此可以使不同调制格式的信号(例如，PM-BPSK、PM-QPSK、PM-4PAM、PM-8PAM、PM-16QAM、PM-32QAM、PM-64QAM)变换到同一个圆上，因此，可以在不知道信号调制格式的前提下对其进行偏振解复用，因此本发明提供的方法与调制格式无关，适用于网络流量愈加动态、多样和不可预测的应用场景。而且，本发明提供的方法相对于传统的RDA，具有简单、抗噪能力强、鲁棒性强的优点，相对于传统的DD-LMS，具有容易实现、无需做波载相位恢复的优点。因此本发明提供的方法不仅适用于相干光通信系统，对整个光通信系统或其他特定应用场景也适用。In the adaptive blind polarization demultiplexing processing method provided by the present invention, when calculating the tap coefficient at the next moment, it is necessary to decompose each polarization state complex signal in the output signal at the previous moment at the next moment into a real part signal and For the imaginary part signal, coordinate transformation is performed on the decomposed signal, the error is calculated, the error is combined to obtain a complex error, and then the tap coefficient at the next moment is calculated according to the complex error, and when the adaptive FIR filter receives the input signal at the next moment Then, the tap coefficient at the next moment can be used to perform polarization demultiplexing processing. Since the present invention decomposes the complex signal into a real part and an imaginary part, and then performs coordinate transformation, signals of different modulation formats (for example, PM-BPSK, PM-QPSK, PM-4PAM, PM-8PAM, PM-16QAM, PM-32QAM, PM-64QAM) to the same circle, therefore, it can be demultiplexed under the premise of not knowing the signal modulation format, so the method provided by the present invention has nothing to do with the modulation format, and is suitable for network traffic More dynamic, diverse and unpredictable application scenarios. Moreover, compared with traditional RDA, the method provided by the present invention has the advantages of simplicity, strong anti-noise ability and strong robustness. Compared with traditional DD-LMS, it has the advantages of easy implementation and no need for carrier phase recovery. Therefore, the method provided by the present invention is not only applicable to coherent optical communication systems, but also applicable to the entire optical communication system or other specific application scenarios.

可理解的是，自适应FIR滤波器在每一时刻的输出信号中可以包括x偏振态复数信号和y偏振态复数信号。It can be understood that the output signal of the adaptive FIR filter at each moment may include an x polarization state complex signal and a y polarization state complex signal.

在具体实施时，可以采用下式对x偏振态复数信号的实部信号进行坐标变换：During specific implementation, the real part signal of the x polarization state complex signal can be coordinate transformed by using the following formula:

x′_i(k)＝x_i(k)-4·sign[x_i(k)]-2·sign{x_i(k)-4·sign[x_i(k)]} (1)x' _i (k)= _xi (k)-4·sign[ _xi (k)]-2·sign{ _xi (k)-4·sign[ _xi (k)]} (1)

其中，x_i(k)为所述自适应FIR滤波器在k时刻的输出信号中x偏振态复数信号的实部信号，x′_i(k)为对实部信号x_i(k)坐标变换后得到的实数信号。Wherein, x _i (k) is the real part signal of the x polarization state complex signal in the output signal of the adaptive FIR filter at time k, and x' _i (k) is the coordinate transformation of the real part signal x _i (k) The resulting real signal.

同样的，可以采用下式对x偏振态复数信号的虚部信号进行坐标变换：Similarly, the following formula can be used to transform the imaginary part signal of the x polarization state complex signal:

x′_q(k)＝x_q(k)-4·sign[x_q(k)]-2·sign{x_q(k)-4·sign[x_q(k)]} (2)x′ _q (k)＝x _q (k)-4·sign[x _q (k)]-2·sign{x _q (k)-4·sign[x _q (k)]} (2)

其中，x_q(k)为所述自适应FIR滤波器在k时刻的输出信号中x偏振态复数信号的虚部信号，x′_q(k)为对虚部信号x_q(k)进行坐标变换后得到的实数信号。Wherein, x _q (k) is the imaginary part signal of the x polarization state complex signal in the output signal of the adaptive FIR filter at time k, and x′ _q (k) is the coordinate of the imaginary part signal x _q (k) The transformed real signal.

同样的，可以采用下式对y偏振态复数信号的实部信号进行坐标变换：Similarly, the following formula can be used to transform the real part of the y-polarized complex signal:

y′_i(k)＝y_i(k)-4·sign[y_i(k)]-2·sign{y_i(k)-4·sign[y_i(k)]} (3)y′ _i (k)=y _i (k)-4·sign[y _i (k)]-2·sign{y _i (k)-4·sign[y _i (k)]} (3)

其中，y_i(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的实部信号，y′_i(k)为对实部信号y_i(k)进行坐标变换后得到的实数信号。Wherein, y _i (k) is the real part signal of the y-polarized complex signal in the output signal of the adaptive FIR filter at time k, and y' _i (k) is the coordinate of the real part signal y _i (k) The transformed real signal.

同样的，可以采用下式对y偏振态复数信号的虚部信号进行坐标变换：Similarly, the coordinate transformation of the imaginary part signal of the y-polarized complex signal can be carried out by using the following formula:

y′_q(k)＝y_q(k)-4·sign[y_q(k)]-2·sign{y_q(k)-4·sign[y_q(k)]} (4)y′ _q (k)=y _q (k)-4·sign[y _q (k)]-2·sign{y _q (k)-4·sign[y _q (k)]} (4)

其中，y_q(k)为所述自适应FIR滤波器在k时刻的输出信号中y偏振态复数信号的虚部信号，y′_q(k)为对虚部信号y_q(k)进行坐标变换后得到的实数信号。Wherein, y _q (k) is the imaginary part signal of the y polarization state complex signal in the output signal of the adaptive FIR filter at k moment, and y′ _q (k) is the coordinate of the imaginary part signal y _q (k) The transformed real signal.

可以理解的是，sign()为取符号函数。It can be understood that sign() is a sign function.

在具体实施时，可以采用下式计算x偏振态复数信号的实部信号的误差：In specific implementation, the following formula can be used to calculate the error of the real part signal of the x polarization state complex signal:

ε_xi(k)＝x′_i(k)(|R_x′i|²-|x′_i(k)|²) (5)ε _xi (k)=x′ _i (k)(|R _x′i | ² -|x′ _i (k)| ² ) (5)

其中，R_x′i为实部信号x_i(k)的参考模值，ε_xi(k)为实部信号x_i(k)的误差。Wherein, R _x'i is the reference modulus value of the real part signal _xi (k), and ε _xi (k) is the error of the real part signal _xi (k).

同样的，可以采用下式计算x偏振态复数信号的虚部信号的误差：Similarly, the following formula can be used to calculate the error of the imaginary part signal of the x polarization state complex signal:

ε_xq(k)＝x′_q(k)(|R_x′q|²-|x′_q(k)|²) (6)ε _xq (k)=x′ _q (k)(|R _x′q | ² -|x′ _q (k)| ² ) (6)

其中，R_x′q为虚部信号x_q(k)的参考模值，ε_xq(k)为虚部信号x_q(k)的误差。Among them, R _x'q is the reference modulus value of the imaginary part signal x _q (k), and ε _xq (k) is the error of the imaginary part signal x _q (k).

同样的，可以采用下式计算y偏振态复数信号的实部信号的误差：Similarly, the error of the real part signal of the y-polarized complex signal can be calculated by the following formula:

ε_yi(k)＝y′_i(k)(|R_y′i|²-|y′_i(k)|²) (7)ε _yi (k)=y′ _i (k)(|R _y′i | ² -|y′ _i (k)| ² ) (7)

其中，R_y′i为实部信号y_i(k)的参考模值，ε_yi(k)为实部信号y_i(k)的误差。Among them, R _y′i is the reference modulus value of the real part signal y _i (k), and ε _yi (k) is the error of the real part signal y _i (k).

同样的，可以采用下式计算y偏振态复数信号虚部信号的误差：Similarly, the following formula can be used to calculate the error of the imaginary part signal of the y-polarized complex signal:

ε_yq(k)＝y′_q(k)(|R_y′q|²-|y′_q(k)|²) (8)ε _yq (k)=y′ _q (k)(|R _y′q | ² -|y′ _q (k)| ² ) (8)

上式中的R_x′i、R_x′q、R_y′i、R_y′q的计算方法如下：The calculation methods of R x′i , R _x′q , R _y′i , and R _y′q in the _above formula are as follows:

其中，可以采用下式计算实部信号x_i(k)的参考模值R_x′i：Among them, the reference modulus R _x′i of the real part signal x _i (k) can be calculated by the following formula:

${R R}_{{x x}^{' '} i i} = = \sqrt{\frac{E E. {{{| | {s the s}_{x x i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x i i} ((k k))]] - - 22 \cdot \cdot s the s i i g g n no {{{s the s}_{x x i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{x x i i} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{x x i i} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{x x i i} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{x x i i} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{x x i i} ((k k))]]}} | |}^{22}}}}} - - - - - - ((99))$

式中，s_xi(k)为第k时刻x偏振态理想星座下复数信号的实部；和/或In the formula, s _xi (k) is the real part of the complex signal under the ideal constellation of polarization state x at the kth moment; and/or

其中，采用下式计算虚部信号x_q(k)的参考模值R_x′q：Among them, the reference modulus R _x′q of the imaginary part signal x _q (k) is calculated by the following formula:

${R R}_{{x x}^{' '} q q} = = \sqrt{\frac{E E. {{{| | {s the s}_{x x q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{x x q q} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{x x q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{x x q q} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{x x q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{x x q q} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{x x q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{x x q q} ((k k))]]}} | |}^{22}}}}} - - - - - - ((1010))$

式中，s_xq(k)为第k时刻x偏振态理想星座下复数信号的虚部；和/或In the formula, s _xq (k) is the imaginary part of the complex signal under the ideal constellation of x polarization state at the kth moment; and/or

其中，采用下式计算实部信号y_i(k)的参考模值R_y′i：Among them, the reference modulus R _y′i of the real part signal y _i (k) is calculated by the following formula:

${R R}_{{y the y}^{' '} i i} = = \sqrt{\frac{E E. {{{| | {s the s}_{y the y i i} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]] - - 22 \cdot \cdot s the s i i g g n no {{{s the s}_{y the y i i} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{y the y i i} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]] - - 22 \cdot \cdot s the s i i g g n no {{{s the s}_{y the y i i} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{y the y i i} ((k k))]]}} | |}^{22}}}}} - - - - - - ((1111))$

式中，s_yi(k)为第k时刻y偏振态理想星座下复数信号的实部；和/或In the formula, s _yi (k) is the real part of the complex signal under the ideal constellation of y polarization state at the kth moment; and/or

其中，采用下式计算虚部信号y_q(k)的参考模值R_y′q：Among them, the reference modulus R _y′q of the imaginary part signal y _q (k) is calculated by the following formula:

${R R}_{{y the y}^{' '} q q} = = \sqrt{\frac{E E. {{{| | {s the s}_{y the y q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{y the y q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]]}} | |}^{44}}}}{E E. {{{| | {s the s}_{y the y q q} ((k k)) - - 44 \cdot \cdot s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]] - - 22 \cdot &Center Dot; s the s i i g g n no {{{s the s}_{y the y q q} ((k k)) - - 44 \cdot &Center Dot; s the s i i g g n no [[{s the s}_{y the y q q} ((k k))]]}} | |}^{22}}}}} - - - - - - ((1212))$

式中，s_yq(k)为第k时刻y偏振态理想星座下复数信号的虚部。In the formula, s _yq (k) is the imaginary part of the complex signal under the ideal constellation of y-polarization state at the kth moment.

在具体实施时，可以采用下式计算x偏振态的复数误差：In specific implementation, the complex error of the x polarization state can be calculated by the following formula:

ε_x(k)＝ε_xi(k)+jε_xq(k) (13)ε _x (k) = ε _xi (k) + jε _xq (k) (13)

其中，ε_x(k)为x偏振态的复数误差。where ε _x (k) is the complex error of the x polarization state.

同样的，可以采用下式计算y偏振态的复数误差：Similarly, the complex error of the y polarization state can be calculated using the following formula:

ε_y(k)＝ε_yi(k)+jε_yq(k) (14)ε _y (k) = ε _yi (k) + jε _yq (k) (14)

在具体实施时，可以采用下式计算第k+1时刻的抽头系数：In specific implementation, the following formula can be used to calculate the tap coefficient at the k+1th moment:

H_k+1,xx(m)＝H_k,xx(m)+μ·ε_x(k)·x_in(k-m)^* (15)H _k+1,xx (m)＝H _k,xx (m)+μ·ε _x (k)·x _in (km) ^* (15)

H_k+1,xy(m)＝H_k,xy(m)+μ·ε_x(k)·y_in(k-m)^* (16)H _k+1,xy (m)＝H _k,xy (m)+μ·ε _x (k)·y _in (km) ^* (16)

H_k+1,yx(m)＝H_k,yx(m)+μ·ε_y(k)·x_in(k-m)^* (17)H _k+1,yx (m)＝H _k,yx (m)+μ·ε _y (k)·x _in (km) ^* (17)

H_k+1,yy(m)＝H_k,yy(m)+μ·ε_y(k)·y_in(k-m)^* (18)H _k+1,yy (m)＝H _k,yy (m)+μ·ε _y (k)·y _in (km) ^* (18)

在具体实施时，可以采用下式对第k+1时刻的输入信号进行偏振解复用处理：In specific implementation, the following formula can be used to perform polarization demultiplexing processing on the input signal at the k+1th moment:

$[\begin{matrix} x x ((k k + + 11)) \\ y the y ((k k + + 11)) \end{matrix}] = = [\begin{matrix} {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, x x x x} ((m m)) \cdot \cdot {x x}_{i i n no} ((k k + + 11 - - m m)) + + {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, x x y the y} ((m m)) \cdot \cdot {y the y}_{i i n no} ((k k + + 11 - - m m)) \\ {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, y the y x x} ((m m)) \cdot \cdot {x x}_{i i n no} ((k k + + 11 - - m m)) + + {Σ Σ}_{m m = = - - N N}^{N N} {H h}_{k k + + 11,, y the y y the y} ((m m)) \cdot &Center Dot; {y the y}_{i i n no} ((k k + + 11 - - m m)) \end{matrix}] - - - - - - ((1919))$

可理解的是，m的取值范围为[-N，N]。Understandably, the value range of m is [-N, N].

在具体实施时，还可以在将信号输入至滤波器之前对信号进行正交不平衡补偿、固定色散补偿和时钟同步等预处理，然后将预处理后的信号输入至滤波器进行处理，并输出，当然对于输出信号还可以进行频偏估计与补偿、相偏估计与补偿、符号反映射等处理，最终得到发端原始比特序列。In specific implementation, preprocessing such as quadrature imbalance compensation, fixed dispersion compensation, and clock synchronization can also be performed on the signal before the signal is input to the filter, and then the preprocessed signal is input to the filter for processing and output , of course, frequency offset estimation and compensation, phase offset estimation and compensation, symbol inverse mapping and other processing can also be performed on the output signal, and finally the original bit sequence of the transmitting terminal can be obtained.

通过实验证明，采用本发明提供的方法对信号进行偏解振处理，相对于传统的CMA、RDE等方法，都具有性能优越的特点。It is proved by experiments that using the method provided by the invention to perform polarization and devibration processing on signals has the characteristics of superior performance compared with traditional methods such as CMA and RDE.

基于相同的发明构思，本发明还提供一种自适应盲偏振解复用处理装置，该装置包括：Based on the same inventive concept, the present invention also provides an adaptive blind polarization demultiplexing processing device, which includes:

所述自适应FIR滤波器，用于根据所述当前时刻的抽头系数，对当前时刻的输入信号进行偏振解复用处理，并将偏振解复用处理后的信号输出；The adaptive FIR filter is used to perform polarization demultiplexing processing on the input signal at the current time according to the tap coefficient at the current time, and output the signal after the polarization demultiplexing processing;

本领域普通技术人员可以理解：实现上述方法实施例的全部或者部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储在计算机可读取的存储介质中，该程序在执行时，执行包括上述方法实施例的步骤。Those of ordinary skill in the art can understand that all or part of the steps for realizing the above-mentioned method embodiments can be completed by hardware related to program instructions, and the aforementioned program can be stored in a computer-readable storage medium. When the program is executed, the The steps of the above method embodiments are included.

本发明的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。In the description of the invention, numerous specific details are set forth. However, it is understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure the understanding of this description.

以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解；其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand; Modifications are made to the recorded technical solutions, or equivalent replacements are made to some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. An adaptive blind polarization demultiplexing method, comprising:

determining a tap coefficient of the self-adaptive butterfly FIR filter at the current moment;

according to the tap coefficient at the current moment, the self-adaptive butterfly FIR filter carries out polarization demultiplexing processing on the input signal at the current moment and outputs the signal after the polarization demultiplexing processing;

wherein the determining the tap coefficient of the adaptive butterfly FIR filter at the current time comprises:

decomposing each polarization state complex signal in the output signal of the self-adaptive butterfly FIR filter at the last moment to obtain a real part signal and an imaginary part signal of the polarization state complex signal;

respectively carrying out coordinate transformation on a real part signal and an imaginary part signal of the polarization state complex signal;

calculating the error of the real part signal according to the reference module value of the real part signal of the polarization state complex signal and the real signal obtained after coordinate transformation, and calculating the error of the imaginary part signal according to the reference module value of the imaginary part signal of the polarization state complex signal and the real signal obtained after coordinate transformation;

synthesizing the error of the real part signal and the error of the imaginary part signal to obtain a complex error of the polarization state complex signal;

and calculating the tap coefficient of the current moment according to the tap coefficient of the previous moment and each polarization state complex signal.

2. The method of claim 1, wherein the output signal of the adaptive FIR filter at each time instant comprises an x-polarization state complex signal and a y-polarization state complex signal; in a corresponding manner, the first and second optical fibers are,

and (3) carrying out coordinate transformation on the real part signal of the x polarization state complex signal by adopting the following formula:

x′_i(k)＝x_i(k)-4·sign[x_i(k)]-2·sign{x_i(k)-4·sign[x_i(k)]}

wherein x is_i(k) Is the real part signal x 'of the x polarization state complex signal in the output signal of the adaptive butterfly FIR filter at the time k'_i(k) For the real part signal x_i(k) Obtaining a real number signal after coordinate transformation; and/or

And (3) carrying out coordinate transformation on the imaginary part signal of the x polarization state complex signal by adopting the following formula:

x′_q(k)＝x_q(k)-4·sign[x_q(k)]-2·sign{x_q(k)-4·sign[x_q(k)]}

wherein x is_q(k) For the adaptive FIR filter at time kImaginary signal, x 'of complex signal of x polarization state in output signal'_q(k) For the imaginary signal x_q(k) Carrying out coordinate transformation to obtain a real number signal; and/or

And (3) carrying out coordinate transformation on the real part signal of the y polarization state complex signal by adopting the following formula:

y′_i(k)＝y_i(k)-4·sign[y_i(k)]-2·sign{y_i(k)-4·sign[y_i(k)]}

wherein, y_i(k) Is the real part signal, y 'of the y polarization state complex signal in the output signal of the adaptive FIR filter at the time point of k'_i(k) For the real part signal y_i(k) Carrying out coordinate transformation to obtain a real number signal; and/or

And (3) carrying out coordinate transformation on the imaginary part signal of the y polarization state complex signal by adopting the following formula:

y′_q(k)＝y_q(k)-4·sign[y_q(k)]-2·sign{y_q(k)-4·sign[y_q(k)]}

wherein, y_q(k) Is the imaginary part signal y 'of the y polarization state complex signal in the output signal of the adaptive FIR filter at the time point of k'_q(k) For the imaginary signal y_q(k) And carrying out coordinate transformation to obtain a real number signal.

3. The method of claim 2,

calculating the error of the real part signal of the x-polarization state complex signal by adopting the following formula:

_xi(k)＝x′_i(k)(|R_x′i|²-|x′_i(k)|²)

wherein R is_x′iIs a real part signal x_i(k) Is determined by the reference modulus value of (a),_xi(k) is a real part signal x_i(k) An error of (2); and/or

The error of the imaginary signal of the complex signal of x polarization state is calculated using the following formula:

_xq(k)＝x′_q(k)(|R_x′q|²-|x′_q(k)|²)

wherein R is_x′qIs an imaginary partSignal x_q(k) Is determined by the reference modulus value of (a),_xq(k) as an imaginary signal x_q(k) An error of (2); and/or

And calculating the error of the real part signal of the y polarization state complex signal by adopting the following formula:

_yi(k)＝y′_i(k)(|R_y′i|²-|y′_i(k)|²)

wherein R is_y′iAs a real part signal y_i(k) Is determined by the reference modulus value of (a),_yi(k) as a real part signal y_i(k) An error of (2); and/or

And calculating the error of the imaginary part signal of the y polarization state complex signal by adopting the following formula:

_yq(k)＝y′_q(k)(|R_y′q|²-|y′_q(k)|²)

wherein R is_y′qIs an imaginary signal y_q(k) Is determined by the reference modulus value of (a),_yq(k) is an imaginary signal y_q(k) The error of (2).

4. The method of claim 3,

the complex error for the x-polarization state is calculated using the following equation:

_x(k)＝_xi(k)+j_xq(k)

wherein,_x(k) complex error for x polarization; and/or

The complex error of the y polarization state is calculated using the following equation:

_y(k)＝_yi(k)+j_yq(k)

wherein,_y(k) is the complex error of the y polarization state.

5. The method of claim 4, wherein the tap coefficient at time k +1 is calculated using the following equation:

H_k+1,xx(m)＝H_k,xx(m)+μ·_x(k)·x_in(k-m)^*

H_k+1,xy(m)＝H_k,xy(m)+μ·_x(k)·y_in(k-m)^*

H_k+1,yx(m)＝H_k,yx(m)+μ·_y(k)·x_in(k-m)^*

H_k+1,yy(m)＝H_k,yy(m)+μ·_y(k)·y_in(k-m)^*

where μ is the step size of the tap update, H_k,xx(m)、H_k,xy(m)、H_k,yx(m)、H_k,yy(m) is the coefficient of the m tap at time k, H_k+1,xx(m)、H_k+1,xy(m)、H_k+1,yx(m)、H_k+1,yy(m) is the coefficient of the m tap at time k +1, x_in(k-m)^*Is the conjugate of the input signal in the x-polarization state at time k, y_in(k-m)^*Is the conjugate of the input signal at the y-polarization state at time instant k.

6. The method of claim 5, wherein the polarization demultiplexing is performed on the input signal at time k +1 using the following equation:

[\begin{matrix} x (k + 1) \\ y (k + 1) \end{matrix}] = [\begin{matrix} Σ_{m = - N}^{N} H_{k + 1, x x} (m) \cdot x_{i n} (k + 1 - m) + Σ_{m = - N}^{N} H_{k + 1, x y} (m) \cdot y_{i n} (k + 1 - m) \\ Σ_{m = - N}^{N} H_{k + 1, y x} (m) \cdot x_{i n} (k + 1 - m) + Σ_{m = - N}^{N} H_{k + 1, y y} (m) \cdot y_{i n} (k + 1 - m) \end{matrix}]

wherein x is_in(k +1-m) is the input signal for the x polarization state at time k +1, y_inAnd (k +1-m) is an input signal of y polarization state at the k +1 moment, x (k +1) is an x polarization state signal obtained by performing polarization demultiplexing on the input signals of the two polarization states at the k +1 moment, y (k +1) is a y polarization state signal obtained by performing polarization demultiplexing on the input signals of the two polarization states at the k +1 moment, and 2N +1 is the number of taps.

7. The method of claim 3,

calculating the real part signal x using the following equation_i(k) Reference modulus value R_x′i：

R_{x^{'} i} = \sqrt{\frac{E {| s_{x i} (k) - 4 \cdot s i g n [s_{x i} (k)] - 2 \cdot s i g n {s_{x i} (k) - 4 \cdot s i g n [s_{x i} (k)]} |^{4}}}{E {| s_{x i} (k) - 4 \cdot s i g n [s_{x i} (k)] - 2 \cdot s i g n {s_{x i} (k) - 4 \cdot s i g n [s_{x i} (k)]} |^{2}}}}

Wherein s is_xi(k) The real part of the complex signal under the x polarization state ideal constellation at the k time; and/or

The imaginary signal x is calculated using the equation_q(k) Reference modulus value R_x′q：

R_{x^{'} q} = \sqrt{\frac{E {| s_{x q} (k) - 4 \cdot s i g n [s_{x q} (k)] - 2 \cdot s i g n {s_{x q} (k) - 4 \cdot s i g n [s_{x q} (k)]} |^{4}}}{E {| s_{x q} (k) - 4 \cdot s i g n [s_{x q} (k)] - 2 \cdot s i g n {s_{x q} (k) - 4 \cdot s i g n [s_{x q} (k)]} |^{2}}}}

Wherein s is_xq(k) The imaginary part of the complex signal under the ideal constellation of the x polarization state at the k time; and/or

Calculating the real part signal y using the following equation_i(k) Reference modulus value R_y′i：

R_{y^{'} i} = \sqrt{\frac{E {| s_{y i} (k) - 4 \cdot s i g n [s_{y i} (k)] - 2 \cdot s i g n {s_{y i} (k) - 4 \cdot s i g n [s_{y i} (k)]} |^{4}}}{E {| s_{y i} (k) - 4 \cdot s i g n [s_{y i} (k)] - 2 \cdot s i g n {s_{y i} (k) - 4 \cdot s i g n [s_{y i} (k)]} |^{2}}}}

Wherein s is_yi(k) The real part of the complex signal under the ideal constellation of the y polarization state at the kth moment; and/or

The imaginary signal y is calculated using the formula_q(k) Reference modulus value R_y′q：

R_{y^{'} q} = \sqrt{\frac{E {| s_{y q} (k) - 4 \cdot s i g n [s_{y q} (k)] - 2 \cdot s i g n {s_{y q} (k) - 4 \cdot s i g n [s_{y q} (k)]} |^{4}}}{E {| s_{y q} (k) - 4 \cdot s i g n [s_{y q} (k)] - 2 \cdot s i g n {s_{y q} (k) - 4 \cdot s i g n [s_{y q} (k)]} |^{2}}}}

Wherein s is_yq(k) The imaginary part of the complex signal in the ideal constellation of the y polarization state at the k time.

8. An adaptive blind polarization demultiplexing device, comprising:

the tap coefficient determining module is used for determining the tap coefficient of the self-adaptive FIR filter at the current moment;

the self-adaptive FIR filter is used for carrying out polarization demultiplexing processing on the input signal at the current moment according to the tap coefficient at the current moment and outputting the signal after the polarization demultiplexing processing;

wherein the tap coefficient determining module is specifically configured to:

decomposing each polarization state complex signal in the output signal of the self-adaptive FIR filter at the last moment to obtain a real part signal and an imaginary part signal of the polarization state complex signal;