CN106301593A - Adaptive blind polarization demultiplexing treating method and apparatus - Google Patents

Abstract

The present invention relates to a kind of adaptive blind polarization demultiplexing treating method and apparatus, the method comprises determining that the auto-adaptive fir filter tap coefficient at current time；According to the tap coefficient of described current time, described auto-adaptive fir filter carries out polarization demultiplexing process to the input signal of current time, and by the signal output after polarization demultiplexing process；Wherein, according to tap coefficient and two polarization state complex signals in a described upper moment, the tap coefficient of described current time was calculated.Due to the present invention by decomposing, coordinate transform can make the multimode signal of different modulating form transform on same circle becomes constant modulus signals, therefore on the premise of not knowing format modulation signal, it can be carried out polarization demultiplexing, therefore the method that the present invention provides is unrelated with modulation format, be applicable to network, modulation format is further flexible, dynamic, the various and uncertain application scenarios of flow.

Description

Adaptive blind polarization demultiplexing treating method and apparatus

Technical field

The present invention relates to technical field of optical fiber communication, especially relate to a kind of adaptive blind polarization demultiplexing processing method and Device.

Background technology

In recent years, along with sharply increasing of network traffics and bandwidth demand, high speed coherent light communication technology has become realization The key technology of distance Large Copacity information transmission.By the technology such as palarization multiplexing and wavelength-division multiplex, can effectively improve light The bandwidth availability ratio of fiber communication system, promotes power system capacity.For high speed coherent light communication, due to fibre-optical dispersion, polarization with The impact of the factors such as machine crosstalk, it is therefore desirable to use multitap self adaptation butterfly FIR filter to carry out residual dispersion, polarizing coating The equilibrium of dispersion and polarization demultiplexing.

At present, polarization demultiplexing method includes data auxiliary and blind process two kinds.Wherein, data householder method is by sending Training sequence carries out channel estimation, occupies partial frequency spectrum resource, reduces system bandwidth utilization rate.But, blind process side Method is as a key technology in point multiplexed optical communications demodulating system partially, it is not necessary to periodically send training sequence, can be only The statistical property receiving signal is relied on to recover originally transmitted signal, therefore without taking too many bandwidth resources.Visible two kinds of sides Method is compared, and blind processing method shows advantage especially.Due to data center and the development of cloud computing, network traffics will become further to move State, various and unpredictable, support that the elastic optical transceiver of more modulation form becomes and elastic optical network effectively carries burst Dynamically change network traffics, improve the availability of frequency spectrum, the key of optimization network resource usage.It is therefore desirable to provide a kind of and adjust The blind processing method that form processed is unrelated.

Summary of the invention

For disadvantages described above, the present invention provides a kind of adaptive blind polarization demultiplexing treating method and apparatus, processing procedure Unrelated with modulation format, it is adaptable to dynamic, the various and uncertain application scenarios of network traffics.

The adaptive blind polarization demultiplexing processing method that the present invention provides includes:

Determine the auto-adaptive fir filter tap coefficient at current time；

According to the tap coefficient of described current time, the input signal of current time is entered by described auto-adaptive fir filter Row polarization demultiplexing processes, and by the signal output after polarization demultiplexing process；

Wherein, described determine the auto-adaptive fir filter tap coefficient at current time, including:

The described auto-adaptive fir filter each polarization state complex signal in the output signal in a upper moment was carried out point Solve, obtain solid part signal and the imaginary signals of this polarization state complex signal；

Solid part signal and imaginary signals to this polarization state complex signal carry out coordinate transform respectively；

The real number signal obtained after the reference modulus value of the solid part signal according to this polarization state complex signal and coordinate transform, meter Calculate the error of this solid part signal, and according to the reality after the reference modulus value of the imaginary signals of this polarization state complex signal and coordinate transform Number signal calculates the error of this imaginary signals；

The error of the error of this solid part signal He this imaginary signals is synthesized, obtains answering of this polarization state complex signal Number error；

Tap coefficient according to a described upper moment and each polarization state complex signal, calculate the tap of described current time Coefficient.

Optionally, the described self adaptation butterfly FIR filter output signal in each moment includes that x-polarisation state plural number is believed Number and y-polarisation state complex signal；Corresponding,

Use following formula that the solid part signal of x-polarisation state complex signal carries out coordinate transform:

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

Wherein, x_iK () is described self adaptation butterfly FIR filter x-polarisation state complex signal in the output signal in k moment Solid part signal, x '_iK () is to solid part signal x_iThe real number signal obtained after (k) coordinate transform；And/or

Use following formula that the imaginary signals of x-polarisation state complex signal carries out coordinate transform:

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

Wherein, x_qK () is described self adaptation butterfly FIR filter x-polarisation state complex signal in the output signal in k moment Imaginary signals, x '_qK () is to imaginary signals x_qK real number signal that () obtains after carrying out coordinate transform；And/or

Use following formula that the solid part signal of y-polarisation state complex signal carries out coordinate transform:

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

Wherein, y_iK () is described self adaptation butterfly FIR filter y-polarisation state complex signal in the output signal in k moment Solid part signal, y '_iK () is to solid part signal y_iK real number signal that () obtains after carrying out coordinate transform；And/or

Use following formula that the imaginary signals of y-polarisation state complex signal carries out coordinate transform:

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

Wherein, y_qK () is described self adaptation butterfly FIR filter y-polarisation state complex signal in the output signal in k moment Imaginary signals, y '_qK () is to imaginary signals y_qK real number signal that () obtains after carrying out coordinate transform.

Optionally, the error of the solid part signal of employing following formula calculating x-polarisation state complex signal:

ε_xi(k)=x '_i(k)(|R_x′i|²-|x′_i(k)|²)

Wherein, R_x′iFor solid part signal x_iThe reference modulus value of (k), ε_xiK () is solid part signal x_iThe error of (k)；And/or

The error of the imaginary signals of employing following formula calculating x-polarisation state complex signal:

ε_xq(k)=x '_q(k)(|R_x′q|²-|x′_q(k)|²)

Wherein, R_x′qFor imaginary signals x_qThe reference modulus value of (k), ε_xqK () is imaginary signals x_qThe error of (k)；And/or

The error of the solid part signal of employing following formula calculating y-polarisation state complex signal:

ε_yi(k)=y '_i(k)(|R_y′i|²-|y′_i(k)|²)

Wherein, R_y′iFor solid part signal y_iThe reference modulus value of (k), ε_yiK () is solid part signal y_iThe error of (k)；And/or

The error of employing following formula calculating y-polarisation state complex signal imaginary signals:

ε_yq(k)=y '_q(k)(|R_y′q|²-|y′_q(k)|²)

Wherein, R_y′qFor imaginary signals y_qThe reference modulus value of (k), ε_yqK () is imaginary signals y_qThe error of (k).

Optionally, the plural error of employing following formula calculating x-polarisation state:

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

Wherein, ε_xK () is the plural error of x-polarisation state；And/or

The plural error of employing following formula calculating y-polarisation state:

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

Wherein, ε_yK () is the plural error of y-polarisation state.

Optionally, the tap coefficient in employing following formula calculating kth+1 moment:

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)^*

Wherein, μ is the step-length that tap updates, H_k,xx(m)、H_k,xy(m)、H_k,yx(m)、H_k,yyM () is in kth moment m-th The coefficient of tap, H_k+1,xx(m)、H_k+1,xy(m)、H_k+1,yx(m)、H_k+1,yyM () is the coefficient in kth+1 moment m-th tap, x_in (k-m)^*For the conjugation of the input signal in kth moment x-polarisation state, y_in(k-m)^*For the input signal in kth moment y-polarisation state Conjugation.

Optionally, use following formula that the input signal in kth+1 moment carries out polarization demultiplexing process:

$\left[\begin{array}{c}x(k+1)\\ y(k+1)\end{array}\right]=\left[\begin{array}{c}\underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,xx}\left(m\right)\·x_{in}(k+1-m)+\underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,xy}\left(m\right)\·y_{in}(k+1-m)\\ \underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,yx}\left(m\right)\·x_{in}(k+1-m)+\underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,yy}\left(m\right)\·y_{in}(k+1-m)\end{array}\right]$

Wherein, x_in(k+1-m) it is the input signal in kth+1 moment x-polarisation state, y_in(k+1-m) it is at kth+1 moment y The input signal of polarization state, x (k+1) is for obtaining after kth+1 moment x-polarisation state input signal is carried out polarization demultiplexing process Signal, y (k+1) be kth+1 moment y-polarisation state input signal is carried out polarization demultiplexing process after the signal that obtains, 2N+1 is Tap number.

Optionally, following formula is used to calculate solid part signal x_iReference modulus value R of (k)_x′i:

$R_{x^{\′}i}=\sqrt{\frac{E\left\{{|s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]-2\·sign\{s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]-2\·sign\{s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]\left\}\right|}^2\right\}}}$

Wherein, s_xiK () is the real part of complex signal under kth moment x-polarisation state both ideal constellation；And/or

Following formula is used to calculate imaginary signals x_qReference modulus value R of (k)_x′q:

$R_{x^{\′}q}=\sqrt{\frac{E\left\{{|s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]-2\·sign\{x_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]-2\·sign\{x_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]\left\}\right|}^2\right\}}}$

Wherein, s_xqK () is the imaginary part of complex signal under kth moment x-polarisation state both ideal constellation；And/or

Following formula is used to calculate solid part signal y_iReference modulus value R of (k)_y′i:

$R_{y^{\′}i}=\sqrt{\frac{E\left\{{|s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]-2\·sign\{s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]-2\·sign\{s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]\left\}\right|}^2\right\}}}$

Wherein, s_yiK () is the real part of complex signal under kth moment y-polarisation state both ideal constellation；And/or

Following formula is used to calculate imaginary signals y_qReference modulus value R of (k)_y′q:

$R_{y^{\′}q}=\sqrt{\frac{E\left\{{|s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]-2\·sign\{s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]-2\·sign\{s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]\left\}\right|}^2\right\}}}$

Wherein, s_yqK () is the imaginary part of complex signal under kth moment y-polarisation state both ideal constellation.

The adaptive blind polarization demultiplexing processing means that the present invention provides includes:

Tap coefficient determines module, for determining the auto-adaptive fir filter tap coefficient at current time；

Described auto-adaptive fir filter, for the tap coefficient according to described current time, the input to current time is believed Number carry out polarization demultiplexing process, and the signal output after polarization demultiplexing is processed；

Wherein, described tap coefficient determine module specifically for:

The described self adaptation butterfly FIR filter each polarization state complex signal in the output signal in a upper moment is entered Row decomposes, and obtains solid part signal and the imaginary signals of this polarization state complex signal；

Solid part signal and imaginary signals to this polarization state complex signal carry out coordinate transform respectively；

The real number signal obtained after the reference modulus value of the solid part signal according to this polarization state complex signal and coordinate transform, meter Calculate the error of this solid part signal, and according to the reality after the reference modulus value of the imaginary signals of this polarization state complex signal and coordinate transform Number signal calculates the error of this imaginary signals；

The error of the error of this solid part signal He this imaginary signals is synthesized, obtains answering of this polarization state complex signal Number error；

Tap coefficient according to a described upper moment and each polarization state complex signal, calculate the tap of described current time Coefficient.

The adaptive blind polarization demultiplexing processing method that the present invention provides, needs when calculating the tap coefficient of subsequent time Each polarization state complex signal in the output signal in a upper moment of this subsequent time is decomposed into solid part signal and imaginary part Signal, carries out coordinate transform to the signal after decomposing, seeks error, and error merging obtains plural number error, further according to plural number error meter Calculate the tap coefficient of this subsequent time, just can utilize when auto-adaptive fir filter receives the input signal of this subsequent time The tap coefficient of this subsequent time carries out the process work of polarization demultiplexing.Owing to the present invention is permissible by decomposition, coordinate transform The multimode signal making different modulating form transforms to becomes constant modulus signals on same circle, therefore, it can not knowing that signal is adjusted On the premise of form processed, it being carried out polarization demultiplexing, the method that therefore present invention provides is unrelated with modulation format, it is adaptable to net In network, modulation format is further flexible, dynamic, the various and uncertain application scenarios of flow.

Accompanying drawing explanation

In order to be illustrated more clearly that disclosure embodiment or technical scheme of the prior art, below will be to embodiment or existing In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this Discloseder embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to Other accompanying drawing is obtained according to these figures.

Fig. 1 shows the schematic flow sheet of adaptive blind polarization demultiplexing processing method in one embodiment of the invention.

Detailed description of the invention

Below in conjunction with the accompanying drawing in disclosure embodiment, the technical scheme in disclosure embodiment is carried out clear, complete Describe, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments wholely.Based on Embodiment in the disclosure, it is every other that those of ordinary skill in the art are obtained under not making creative work premise Embodiment, broadly falls into the scope of disclosure protection.

The present invention provides a kind of adaptive blind polarization demultiplexing processing method, as it is shown in figure 1, the method includes:

Determine the auto-adaptive fir filter tap coefficient at current time；

According to the tap coefficient of described current time, the input signal of current time is entered by described auto-adaptive fir filter Row polarization demultiplexing processes, and by the signal output after polarization demultiplexing process；

Wherein, described determine the self adaptation butterfly FIR filter tap coefficient at current time, including:

The described self adaptation butterfly FIR filter each polarization state complex signal in the output signal in a upper moment is entered Row decomposes, and obtains solid part signal and the imaginary signals of this polarization state complex signal；

Solid part signal and imaginary signals to this polarization state complex signal carry out coordinate transform respectively；

The real number signal obtained after the reference modulus value of the solid part signal according to this polarization state complex signal and coordinate transform, meter Calculate the error of this solid part signal, and according to the reality after the reference modulus value of the imaginary signals of this polarization state complex signal and coordinate transform Number signal calculates the error of this imaginary signals；

The error of the error of this solid part signal He this imaginary signals is synthesized, obtains answering of this polarization state complex signal Number error；

Tap coefficient according to a described upper moment and each polarization state complex signal, calculate the tap of described current time Coefficient.

It will be appreciated that auto-adaptive fir filter is adaptive to have limit for length's unit impulse response wave filter.

It will be appreciated that auto-adaptive fir filter carries out employing during polarization demultiplexing to the input signal of current time Current tap coefficient is that the output signal to a upper moment auto-adaptive fir filter carries out certain process and obtains, such as, right When the input signal of the 2nd moment auto-adaptive fir filter carries out polarization demultiplexing, the tap coefficient of employing is according to the 1st moment The output signal of auto-adaptive fir filter carries out processing and obtains, the input signal to the 3rd moment auto-adaptive fir filter Carrying out the tap coefficient of employing during polarization demultiplexing is that the output signal according to the 2nd moment auto-adaptive fir filter processes Obtain, the like ....Wherein, the input signal at the 1st moment auto-adaptive fir filter is adopted when carrying out polarization demultiplexing Initial tap coefficient values arrange voluntarily.Visible, after the tap coefficient in each moment be the output according to previous moment Signal obtains.Renewal continuous to tap coefficient in this way so that it is meet the requirement of polarization demultiplexing.

The adaptive blind polarization demultiplexing processing method that the present invention provides, needs when calculating the tap coefficient of subsequent time Each polarization state complex signal in the output signal in a upper moment of this subsequent time is decomposed into solid part signal and imaginary part Signal, carries out coordinate transform to the signal after decomposing, seeks error, and error merging obtains plural number error, further according to plural number error meter Calculate the tap coefficient of this subsequent time, just can utilize when auto-adaptive fir filter receives the input signal of this subsequent time The tap coefficient of this subsequent time carries out the process work of polarization demultiplexing.Due to the present invention complex signal is decomposed into real part and Imaginary part, then carries out coordinate transform, therefore can make signal (such as, PM-BPSK, PM-QPSK, PM-of different modulating form 4PAM, PM-8PAM, PM-16QAM, PM-32QAM, PM-64QAM) transform on same circle, therefore, it can not knowing letter On the premise of number modulation format, it being carried out polarization demultiplexing, the method that therefore present invention provides is unrelated with modulation format, applicable Dynamic in network traffics, various and uncertain application scenarios.And, the method that the present invention provides is relative to traditional RDA, have simple, anti-noise ability strong, the advantage of strong robustness, relative to traditional DD-LMS, has easily realization, without doing The advantage that ripple carries phase recovery.Therefore the method that the present invention provides is applicable not only to coherent optical communication system, to whole optic communication System or other application-specific scenes are also suitable.

It will be appreciated that auto-adaptive fir filter can include that in the output signal in each moment x-polarisation state plural number is believed Number and y-polarisation state complex signal.

In the specific implementation, can use following formula that the solid part signal of x-polarisation state complex signal carries out coordinate transform:

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

Wherein, x_iK () is described auto-adaptive fir filter reality of x-polarisation state complex signal in the output signal in k moment Portion's signal, x '_iK () is to solid part signal x_iThe real number signal obtained after (k) coordinate transform.

Same, can use following formula that the imaginary signals of x-polarisation state complex signal carries out coordinate transform:

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

Wherein, x_qK () is described auto-adaptive fir filter void of x-polarisation state complex signal in the output signal in k moment Portion's signal, x '_qK () is to imaginary signals x_qK real number signal that () obtains after carrying out coordinate transform.

Same, can use following formula that the solid part signal of y-polarisation state complex signal carries out coordinate transform:

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

Wherein, y_iK () is described auto-adaptive fir filter reality of y-polarisation state complex signal in the output signal in k moment Portion's signal, y '_iK () is to solid part signal y_iK real number signal that () obtains after carrying out coordinate transform.

Same, can use following formula that the imaginary signals of y-polarisation state complex signal carries out coordinate transform:

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

Wherein, y_qK () is described auto-adaptive fir filter void of y-polarisation state complex signal in the output signal in k moment Portion's signal, y '_qK () is to imaginary signals y_qK real number signal that () obtains after carrying out coordinate transform.

It is understood that sign () is for taking sign function.

In the specific implementation, following formula can be used to calculate the error of solid part signal of x-polarisation state complex signal:

ε_xi(k)=x '_i(k)(|R_x′i|²-|x′_i(k)|²) (5)

Wherein, R_x′iFor solid part signal x_iThe reference modulus value of (k), ε_xiK () is solid part signal x_iThe error of (k).

Same, following formula can be used to calculate the error of imaginary signals of x-polarisation state complex signal:

ε_xq(k)=x '_q(k)(|R_x′q|²-|x′_q(k)|²) (6)

Wherein, R_x′qFor imaginary signals x_qThe reference modulus value of (k), ε_xqK () is imaginary signals x_qThe error of (k).

Same, following formula can be used to calculate the error of solid part signal of y-polarisation state complex signal:

ε_yi(k)=y '_i(k)(|R_y′i|²-|y′_i(k)|²) (7)

Wherein, R_y′iFor solid part signal y_iThe reference modulus value of (k), ε_yiK () is solid part signal y_iThe error of (k).

Same, following formula can be used to calculate the error of y-polarisation state complex signal imaginary signals:

ε_yq(k)=y '_q(k)(|R_y′q|²-|y′_q(k)|²) (8)

Wherein, R_y′qFor imaginary signals y_qThe reference modulus value of (k), ε_yqK () is imaginary signals y_qThe error of (k).

R in above formula_x′i、R_x′q、R_y′i、R_y′qComputational methods as follows:

Wherein it is possible to use following formula to calculate solid part signal x_iReference modulus value R of (k)_x′i:

$R_{x^{\′}i}=\sqrt{\frac{E\left\{{|s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]-2\·sign\{s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]-2\·sign\{s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]\left\}\right|}^2\right\}}}---\left(9\right)$

In formula, s_xiK () is the real part of complex signal under kth moment x-polarisation state both ideal constellation；And/or

Wherein, following formula is used to calculate imaginary signals x_qReference modulus value R of (k)_x′q:

$R_{x^{\′}q}=\sqrt{\frac{E\left\{{|s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]-2\·sign\{s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]-2\·sign\{s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]\left\}\right|}^2\right\}}}---\left(10\right)$

In formula, s_xqK () is the imaginary part of complex signal under kth moment x-polarisation state both ideal constellation；And/or

Wherein, following formula is used to calculate solid part signal y_iReference modulus value R of (k)_y′i:

$R_{y^{\′}i}=\sqrt{\frac{E\left\{{|s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]-2\·sign\{s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]-2\·sign\{s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]\left\}\right|}^2\right\}}}---\left(11\right)$

In formula, s_yiK () is the real part of complex signal under kth moment y-polarisation state both ideal constellation；And/or

Wherein, following formula is used to calculate imaginary signals y_qReference modulus value R of (k)_y′q:

$R_{y^{\′}q}=\sqrt{\frac{E\left\{{|s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]-2\·sign\{s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]\left\}\right|}^4\right\}}{E\left\{{|s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]-2\·sign\{s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]\left\}\right|}^2\right\}}}---\left(12\right)$

In formula, s_yqK () is the imaginary part of complex signal under kth moment y-polarisation state both ideal constellation.

In the specific implementation, can use following formula calculate x-polarisation state plural error:

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

Wherein, ε_xK () is the plural error of x-polarisation state.

Same, following formula can be used to calculate the plural error of y-polarisation state:

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

Wherein, ε_yK () is the plural error of y-polarisation state.

In the specific implementation, can use following formula calculate kth+1 moment tap coefficient:

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

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

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

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

Wherein, μ is the step-length that tap updates, H_k,xx(m)、H_k,xy(m)、H_k,yx(m)、H_k,yyM () is in kth moment m-th The coefficient of tap, H_k+1,xx(m)、H_k+1,xy(m)、H_k+1,yx(m)、H_k+1,yyM () is the coefficient in kth+1 moment m-th tap, x_in (k-m)^*For the conjugation of the input signal in kth moment x-polarisation state, y_in(k-m)^*For the input signal in kth moment y-polarisation state Conjugation.

In the specific implementation, can use following formula that the input signal in kth+1 moment carries out polarization demultiplexing process:

$\left[\begin{array}{c}x(k+1)\\ y(k+1)\end{array}\right]=\left[\begin{array}{c}\underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,xx}\left(m\right)\·x_{in}(k+1-m)+\underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,xy}\left(m\right)\·y_{in}(k+1-m)\\ \underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,yx}\left(m\right)\·x_{in}(k+1-m)+\underset{m=-N}{\overset{N}{\Σ}}{H}_{k+1,yy}\left(m\right)\·y_{in}(k+1-m)\end{array}\right]---\left(19\right)$

Wherein, x_in(k+1-m) it is the input signal in kth+1 moment x-polarisation state, y_in(k+1-m) it is at kth+1 moment y The input signal of polarization state, x (k+1) is for obtaining after kth+1 moment x-polarisation state input signal is carried out polarization demultiplexing process Signal, y (k+1) be kth+1 moment y-polarisation state input signal is carried out polarization demultiplexing process after the signal that obtains, 2N+1 is Tap number.

It will be appreciated that the span of m is [-N, N].

In the specific implementation, it is also possible to input a signal into signal is carried out before wave filter quadrature imbalance compensation, The fixing pretreatment such as dispersion compensation and clock synchronization, then processes the input of pretreated signal to wave filter, and defeated Go out, frequency deviation can also be carried out certainly for output signal and estimate to process with compensate, skew is estimated and compensates, symbol reflection is penetrated etc., Finally give original bit sequence of making a start.

Being experimentally confirmed, the method using the present invention to provide carries out inclined solution and shakes process, relative to traditional signal The methods such as CMA, RDE, all have the feature of superior performance.

Based on identical inventive concept, the present invention also provides for a kind of adaptive blind polarization demultiplexing processing means, this device Including:

Tap coefficient determines module, for determining the auto-adaptive fir filter tap coefficient at current time；

Described auto-adaptive fir filter, for the tap coefficient according to described current time, the input to current time is believed Number carry out polarization demultiplexing process, and the signal output after polarization demultiplexing is processed；

Wherein, described tap coefficient determine module specifically for:

The described auto-adaptive fir filter each polarization state complex signal in the output signal in a upper moment was carried out point Solve, obtain solid part signal and the imaginary signals of this polarization state complex signal；

Solid part signal and imaginary signals to this polarization state complex signal carry out coordinate transform respectively；

The real number signal obtained after the reference modulus value of the solid part signal according to this polarization state complex signal and coordinate transform, meter Calculate the error of this solid part signal, and according to the reality after the reference modulus value of the imaginary signals of this polarization state complex signal and coordinate transform Number signal calculates the error of this imaginary signals；

The error of the error of this solid part signal He this imaginary signals is synthesized, obtains answering of this polarization state complex signal Number error；

Tap coefficient according to a described upper moment and each polarization state complex signal, calculate the tap of described current time Coefficient.

One of ordinary skill in the art will appreciate that: all or part of step realizing said method embodiment can be led to The hardware crossing programmed instruction relevant completes, and aforesaid program can be stored in the storage medium of embodied on computer readable, this journey Sequence upon execution, performs to include the step of said method embodiment.

In the description of the present invention, illustrate a large amount of detail.It is to be appreciated, however, that embodiments of the invention are permissible Put into practice in the case of there is no these details.In some instances, it is not shown specifically known method, structure and skill Art, in order to do not obscure the understanding of this description.

Above example only in order to technical scheme to be described, is not intended to limit；Although with reference to previous embodiment The present invention is described in detail, it will be understood by those within the art that；It still can be to aforementioned each enforcement Technical scheme described in example is modified, or wherein portion of techniques feature is carried out equivalent；And these amendment or Replace, do not make the essence of appropriate technical solution depart from the spirit and scope of various embodiments of the present invention technical scheme.

Claims (8)

1. an adaptive blind polarization demultiplexing processing method, it is characterised in that including:

Determine the self adaptation butterfly FIR filter tap coefficient at current time；

According to the tap coefficient of described current time, the input signal of current time is entered by described self adaptation butterfly FIR filter Row polarization demultiplexing processes, and by the signal output after polarization demultiplexing process；

Wherein, described determine the self adaptation butterfly FIR filter tap coefficient at current time, including:

The described self adaptation butterfly FIR filter each polarization state complex signal in the output signal in a upper moment was carried out point Solve, obtain solid part signal and the imaginary signals of this polarization state complex signal；

Solid part signal and imaginary signals to this polarization state complex signal carry out coordinate transform respectively；

The real number signal obtained after the reference modulus value of the solid part signal according to this polarization state complex signal and coordinate transform, calculating should The error of solid part signal, and according to the real number letter after the reference modulus value of the imaginary signals of this polarization state complex signal and coordinate transform Number calculate this imaginary signals error；

The error of the error of this solid part signal He this imaginary signals is synthesized, obtains the plural number of this polarization state complex signal by mistake Difference；

Tap coefficient according to a described upper moment and each polarization state complex signal, calculate the tap system of described current time Number.

Method the most according to claim 1, it is characterised in that described auto-adaptive fir filter is in the output in each moment Signal includes x-polarisation state complex signal and y-polarisation state complex signal；Corresponding,

Use following formula that the solid part signal of x-polarisation state complex signal carries out coordinate transform:

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

Wherein, x_iK () is described self adaptation butterfly FIR filter reality of x-polarisation state complex signal in the output signal in k moment Portion's signal, x '_iK () is to solid part signal x_iThe real number signal obtained after (k) coordinate transform；And/or

Use following formula that the imaginary signals of x-polarisation state complex signal carries out coordinate transform:

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

Wherein, x_qK () is the imaginary part letter of described auto-adaptive fir filter x-polarisation state complex signal in the output signal in k moment Number, x '_qK () is to imaginary signals x_qK real number signal that () obtains after carrying out coordinate transform；And/or

Use following formula that the solid part signal of y-polarisation state complex signal carries out coordinate transform:

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

Wherein, y_iK () is the real part letter of described auto-adaptive fir filter y-polarisation state complex signal in the output signal in k moment Number, y '_iK () is to solid part signal y_iK real number signal that () obtains after carrying out coordinate transform；And/or

Use following formula that the imaginary signals of y-polarisation state complex signal carries out coordinate transform:

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

Wherein, y_qK () is the imaginary part letter of described auto-adaptive fir filter y-polarisation state complex signal in the output signal in k moment Number, y '_qK () is to imaginary signals y_qK real number signal that () obtains after carrying out coordinate transform.

Method the most according to claim 2, it is characterised in that

The error of the solid part signal of employing following formula calculating x-polarisation state complex signal:

ε_xi(k)=x '_i(k)(|R_x′i|²-|x′_i(k)|²)

Wherein, R_x′iFor solid part signal x_iThe reference modulus value of (k), ε_xiK () is solid part signal x_iThe error of (k)；And/or

The error of the imaginary signals of employing following formula calculating x-polarisation state complex signal:

ε_xq(k)=x '_q(k)(|R_x′q|²-|x′_q(k)|²)

Wherein, R_x′qFor imaginary signals x_qThe reference modulus value of (k), ε_xqK () is imaginary signals x_qThe error of (k)；And/or

The error of the solid part signal of employing following formula calculating y-polarisation state complex signal:

ε_yi(k)=y '_i(k)(|R_y′i|²-|y′_i(k)|²)

Wherein, R_y′iFor solid part signal y_iThe reference modulus value of (k), ε_yiK () is solid part signal y_iThe error of (k)；And/or

The error of employing following formula calculating y-polarisation state complex signal imaginary signals:

ε_yq(k)=y '_q(k)(|R_y′q|²-|y′_q(k)|²)

Wherein, R_y′qFor imaginary signals y_qThe reference modulus value of (k), ε_yqK () is imaginary signals y_qThe error of (k).

Method the most according to claim 3, it is characterised in that

The plural error of employing following formula calculating x-polarisation state:

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

Wherein, ε_xK () is the plural error of x-polarisation state；And/or

The plural error of employing following formula calculating y-polarisation state:

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

Wherein, ε_yK () is the plural error of y-polarisation state.

Method the most according to claim 4, it is characterised in that the tap coefficient in employing following formula calculating kth+1 moment:

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)^*

Wherein, μ is the step-length that tap updates, H_k,xx(m)、H_k,xy(m)、H_k,yx(m)、H_k,yyM () is in the m-th tap of kth moment Coefficient, H_k+1,xx(m)、H_k+1,xy(m)、H_k+1,yx(m)、H_k+1,yyM () is the coefficient in kth+1 moment m-th tap, x_in(k- m)^*For the conjugation of the input signal in kth moment x-polarisation state, y_in(k-m)^*For the input signal in kth moment y-polarisation state Conjugation.

Method the most according to claim 5, it is characterised in that use following formula that the input signal in kth+1 moment is carried out partially Shake demultiplexing process:

$\left[\begin{array}{c}x\left(k+1\right)\\ y\left(k+1\right)\end{array}\right]=\left[\begin{array}{c}\underset{m=-N}{\overset{N}{\mathrm{\Σ}}}{H}_{k+1,xx}\left(m\right)\·x_{in}\left(k+1-m\right)+\underset{m=-N}{\overset{N}{\mathrm{\Σ}}}{H}_{k+1,xy}\left(m\right)\·y_{in}\left(k+1-m\right)\\ \underset{m=-N}{\overset{N}{\mathrm{\Σ}}}{H}_{k+1,yx}\left(m\right)\·x_{in}\left(k+1-m\right)+\underset{m=-N}{\overset{N}{\mathrm{\Σ}}}{H}_{k+1,yy}\left(m\right)\·y_{in}\left(k+1-m\right)\end{array}\right]$

Wherein, x_in(k+1-m) it is the input signal in kth+1 moment x-polarisation state, y_in(k+1-m) it is in kth+1 moment y-polarisation The input signal of state, x (k+1) is that after kth+1 moment two polarization state input signal is carried out polarization demultiplexing process, the x that obtains is inclined Polarization state signal, y (k+1) be kth+1 moment two polarization state input signal is carried out polarization demultiplexing process after the y-polarisation state that obtains Signal, 2N+1 is tap number.

Method the most according to claim 3, it is characterised in that

Following formula is used to calculate solid part signal x_iReference modulus value R of (k)_x′i:

$R_{x^{\′}i}=\sqrt{\frac{E\left\{\right|s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]-2\·sign\{s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]\}{|}^4\}}{E\left\{\right|s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]-2\·sign\{s_{xi}\left(k\right)-4\·sign\[s_{xi}\left(k\right)\]\}{|}^2\}}}$

Wherein, s_xiK () is the real part of complex signal under kth moment x-polarisation state both ideal constellation；And/or

Following formula is used to calculate imaginary signals x_qReference modulus value R of (k)_x′q:

$R_{x^{\′}q}=\sqrt{\frac{E\left\{\right|s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]-2\·sign\{s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]\}{|}^4\}}{E\left\{\right|s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]-2\·sign\{s_{xq}\left(k\right)-4\·sign\[s_{xq}\left(k\right)\]\}{|}^2\}}}$

Wherein, s_xqK () is the imaginary part of complex signal under kth moment x-polarisation state both ideal constellation；And/or

Following formula is used to calculate solid part signal y_iReference modulus value R of (k)_y′i:

$R_{y^{\′}i}=\sqrt{\frac{E\left\{\right|s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]-2\·sign\{s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]\}{|}^4\}}{E\left\{\right|s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]-2\·sign\{s_{yi}\left(k\right)-4\·sign\[s_{yi}\left(k\right)\]\}{|}^2\}}}$

Wherein, s_yiK () is the real part of complex signal under kth moment y-polarisation state both ideal constellation；And/or

Following formula is used to calculate imaginary signals y_qReference modulus value R of (k)_y′q:

$R_{y^{\′}q}=\sqrt{\frac{E\left\{\right|s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]-2\·sign\{s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]\}{|}^4\}}{E\left\{\right|s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]-2\·sign\{s_{yq}\left(k\right)-4\·sign\[s_{yq}\left(k\right)\]\}{|}^2\}}}$

Wherein, s_yqK () is the imaginary part of complex signal under kth moment y-polarisation state both ideal constellation.

8. an adaptive blind polarization demultiplexing processing means, it is characterised in that including:

Tap coefficient determines module, for determining the auto-adaptive fir filter tap coefficient at current time；

Described auto-adaptive fir filter, for the tap coefficient according to described current time, enters the input signal of current time Row polarization demultiplexing processes, and by the signal output after polarization demultiplexing process；

Wherein, described tap coefficient determine module specifically for:

The described auto-adaptive fir filter each polarization state complex signal in the output signal in a upper moment is decomposed, Obtain solid part signal and the imaginary signals of this polarization state complex signal；

Solid part signal and imaginary signals to this polarization state complex signal carry out coordinate transform respectively；

The real number signal obtained after the reference modulus value of the solid part signal according to this polarization state complex signal and coordinate transform, calculating should The error of solid part signal, and according to the real number letter after the reference modulus value of the imaginary signals of this polarization state complex signal and coordinate transform Number calculate this imaginary signals error；

The error of the error of this solid part signal He this imaginary signals is synthesized, obtains the plural number of this polarization state complex signal by mistake Difference；

Tap coefficient according to a described upper moment and each polarization state complex signal, calculate the tap system of described current time Number.