CN106301593B - Adaptive blind polarization demultiplexing treating method and apparatus - Google Patents
Adaptive blind polarization demultiplexing treating method and apparatus Download PDFInfo
- Publication number
- CN106301593B CN106301593B CN201610650339.5A CN201610650339A CN106301593B CN 106301593 B CN106301593 B CN 106301593B CN 201610650339 A CN201610650339 A CN 201610650339A CN 106301593 B CN106301593 B CN 106301593B
- Authority
- CN
- China
- Prior art keywords
- signal
- moment
- polarisation state
- sign
- solid part
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/60—Receivers
- H04B10/61—Coherent receivers
- H04B10/616—Details of the electronic signal processing in coherent optical receivers
- H04B10/6165—Estimation of the phase of the received optical signal, phase error estimation or phase error correction
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J14/00—Optical multiplex systems
- H04J14/06—Polarisation multiplex systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Optical Communication System (AREA)
Abstract
The present invention relates to a kind of adaptive blind polarization demultiplexing treating method and apparatus, this method includes：Determine auto-adaptive fir filter in the tap coefficient at current time；According to the tap coefficient at the current time, the auto-adaptive fir filter carries out polarization demultiplexing processing to the input signal at current time, and treated that signal exports by polarization demultiplexing；Wherein, according to the tap coefficient of the last moment and two polarization state complex signals, the tap coefficient at the current time is calculated.Become constant modulus signals on the same circle since the present invention is by decomposing, coordinate transform can be such that the multimode signal of different modulating format transforms to, therefore polarization demultiplexing can be carried out to it under the premise of not knowing format modulation signal, therefore method provided by the invention is unrelated with modulation format, modulation format is further flexible suitable for network, flow dynamic, multiplicity and uncertain application scenarios further.
Description
Technical field
The present invention relates to technical field of optical fiber communication, more particularly, to a kind of adaptive blind polarization demultiplexing processing method and
Device.
Background technique
In recent years, sharply increasing with network flow and bandwidth demand, high speed coherent light communication technology has become realization
The key technology of long range large capacity information transmission.By technologies such as palarization multiplexing and wavelength-division multiplex, light can be effectively improved
The bandwidth availability ratio of fiber communication system, lifting system capacity.For high speed coherent light communication, due to fibre-optical dispersion, polarization with
The influence of the factors such as machine crosstalk, it is therefore desirable to which residual dispersion, polarizing coating are carried out using the adaptive butterfly FIR filter of multi-tap
The equilibrium and polarization demultiplexing of dispersion.
Currently, polarization demultiplexing method includes data auxiliary and blind two kinds of processing.Wherein, data householder method passes through transmission
Training sequence carries out channel estimation, occupies partial frequency spectrum resource, reduces system bandwidth utilization rate.However, blind processing side
Method is as a key technology divided in multiplexed optical communications demodulating system partially, can be only without periodically sending training sequence
Restore originally transmitted signal by the statistical property for receiving signal, there is no need to occupy too many bandwidth resources.It can be seen that two kinds of sides
Method is compared, the especially aobvious advantage of blind processing method.Due to the continuous development of data center and cloud computing, network flow will become to move further
State, multiplicity and unpredictable support the elastic optical transceiver of more modulation format to become effectively carrying burst in elastic optical network
Dynamic change network flow, the key for improving the availability of frequency spectrum, optimizing network resource usage.It is therefore desirable to provide a kind of and adjust
The unrelated blind processing method of format processed.
Summary of the invention
For disadvantages described above, the present invention provides a kind of adaptive blind polarization demultiplexing treating method and apparatus, treatment process
It is unrelated with modulation format, it is suitable for network flow dynamic, multiplicity and uncertain application scenarios further.
Adaptive blind polarization demultiplexing processing method provided by the invention includes：
Determine auto-adaptive fir filter in the tap coefficient at current time；
According to the tap coefficient at the current time, the auto-adaptive fir filter to the input signal at current time into
The processing of row polarization demultiplexing, and polarization demultiplexing treated signal is exported；
Wherein, tap coefficient of the determining auto-adaptive fir filter at current time, including：
Each polarization state complex signal of the auto-adaptive fir filter in the output signal of last moment is divided
Solution, obtains the solid part signal and imaginary signals of the polarization state complex signal；
The solid part signal and imaginary signals of the polarization state complex signal are coordinately transformed respectively；
According to the real number signal obtained after the reference modulus value of the solid part signal of the polarization state complex signal and coordinate transform, meter
The error of the solid part signal is calculated, and according to the reality after the reference modulus value and coordinate transform of the imaginary signals of the polarization state complex signal
Number signal calculates the error of the imaginary signals；
The error of the error of the solid part signal and the imaginary signals is synthesized, answering for the polarization state complex signal is obtained
Number error；
According to the tap coefficient of the last moment and each polarization state complex signal, the tap at the current time is calculated
Coefficient.
Optionally, the adaptive butterfly FIR filter includes that x-polarisation state plural number is believed in the output signal at each moment
Number and y-polarisation state complex signal；It is corresponding,
It is coordinately transformed using solid part signal of the following formula to x-polarisation state complex signal：
x′i(k)=xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
Wherein, xiIt (k) is adaptive butterfly FIR filter x-polarisation state complex signal in the output signal at k moment
Solid part signal, x 'i(k) for solid part signal xi(k) real number signal obtained after coordinate transform；And/or
It is coordinately transformed using imaginary signals of the following formula to x-polarisation state complex signal：
x′q(k)=xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
Wherein, xqIt (k) is adaptive butterfly FIR filter x-polarisation state complex signal in the output signal at k moment
Imaginary signals, x 'q(k) for imaginary signals xq(k) real number signal obtained after being coordinately transformed；And/or
It is coordinately transformed using solid part signal of the following formula to y-polarisation state complex signal：
y′i(k)=yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
Wherein, yiIt (k) is adaptive butterfly FIR filter y-polarisation state complex signal in the output signal at k moment
Solid part signal, y 'i(k) for solid part signal yi(k) real number signal obtained after being coordinately transformed；And/or
It is coordinately transformed using imaginary signals of the following formula to y-polarisation state complex signal：
y′q(k)=yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
Wherein, yqIt (k) is adaptive butterfly FIR filter y-polarisation state complex signal in the output signal at k moment
Imaginary signals, y 'q(k) for imaginary signals yq(k) real number signal obtained after being coordinately transformed.
Optionally, the error of the solid part signal of x-polarisation state complex signal is calculated using following formula：
εxi(k)=x 'i(k)(|Rx′i|2-|x′i(k)|2)
Wherein, Rx′iFor solid part signal xi(k) reference modulus value, εxiIt (k) is solid part signal xi(k) error；And/or
The error of the imaginary signals of x-polarisation state complex signal is calculated using following formula：
εxq(k)=x 'q(k)(|Rx′q|2-|x′q(k)|2)
Wherein, Rx′qFor imaginary signals xq(k) reference modulus value, εxqIt (k) is imaginary signals xq(k) error；And/or
The error of the solid part signal of y-polarisation state complex signal is calculated using following formula：
εyi(k)=y 'i(k)(|Ry′i|2-|y′i(k)|2)
Wherein, Ry′iFor solid part signal yi(k) reference modulus value, εyiIt (k) is solid part signal yi(k) error；And/or
The error of y-polarisation state complex signal imaginary signals is calculated using following formula：
εyq(k)=y 'q(k)(|Ry′q|2-|y′q(k)|2)
Wherein, Ry′qFor imaginary signals yq(k) reference modulus value, εyqIt (k) is imaginary signals yq(k) error.
Optionally, the plural error of x-polarisation state is calculated using following formula：
εx(k)=εxi(k)+jεxq(k)
Wherein, εxIt (k) is the plural error of x-polarisation state；And/or
The plural error of y-polarisation state is calculated using following formula：
εy(k)=εyi(k)+jεyq(k)
Wherein, εyIt (k) is the plural error of y-polarisation state.
Optionally, the tap coefficient at+1 moment of kth is calculated using following formula：
Hk+1,xx(m)=Hk,xx(m)+μ·εx(k)·xin(k-m)*
Hk+1,xy(m)=Hk,xy(m)+μ·εx(k)·yin(k-m)*
Hk+1,yx(m)=Hk,yx(m)+μ·εy(k)·xin(k-m)*
Hk+1,yy(m)=Hk,yy(m)+μ·εy(k)·yin(k-m)*
Wherein, μ is the step-length that tap updates, Hk,xx(m)、Hk,xy(m)、Hk,yx(m)、Hk,yy(m) at m-th of the kth moment
The coefficient of tap, Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yy(m) in the coefficient of m-th of tap of+1 moment of kth, xin
(k-m)*For in the conjugation of the input signal of kth moment x-polarisation state, yin(k-m)*For in the input signal of kth moment y-polarisation state
Conjugation.
Optionally, polarization demultiplexing processing is carried out using input signal of the following formula to+1 moment of kth：
Wherein, xin(k+1-m) in the input signal of+1 moment of kth x-polarisation state, yin(k+1-m) in+1 moment of kth y
The input signal of polarization state, x (k+1) are obtained after carrying out polarization demultiplexing processing to+1 moment of kth x-polarisation state input signal
Signal, y (k+1) are the signal for obtain after polarization demultiplexing processing to+1 moment of kth y-polarisation state input signal, and 2N+1 is
Tap number.
Optionally, solid part signal x is calculated using following formulai(k) reference modulus value Rx′i：
Wherein, sxiIt (k) is the real part of complex signal under kth moment x-polarisation state both ideal constellation；And/or
Imaginary signals x is calculated using following formulaq(k) reference modulus value Rx′q：
Wherein, sxqIt (k) is the imaginary part of complex signal under kth moment x-polarisation state both ideal constellation；And/or
Solid part signal y is calculated using following formulai(k) reference modulus value Ry′i：
Wherein, syiIt (k) is the real part of complex signal under kth moment y-polarisation state both ideal constellation；And/or
Imaginary signals y is calculated using following formulaq(k) reference modulus value Ry′q：
Wherein, syqIt (k) is the imaginary part of complex signal under kth moment y-polarisation state both ideal constellation.
Adaptive blind polarization demultiplexing processing unit provided by the invention includes：
Tap coefficient determining module, for determining auto-adaptive fir filter in the tap coefficient at current time；
The auto-adaptive fir filter believes the input at current time for the tap coefficient according to the current time
Number polarization demultiplexing processing is carried out, and polarization demultiplexing treated signal is exported；
Wherein, the tap coefficient determining module is specifically used for：
To each polarization state complex signal of the adaptive butterfly FIR filter in the output signal of last moment into
Row decomposes, and obtains the solid part signal and imaginary signals of the polarization state complex signal；
The solid part signal and imaginary signals of the polarization state complex signal are coordinately transformed respectively；
According to the real number signal obtained after the reference modulus value of the solid part signal of the polarization state complex signal and coordinate transform, meter
The error of the solid part signal is calculated, and according to the reality after the reference modulus value and coordinate transform of the imaginary signals of the polarization state complex signal
Number signal calculates the error of the imaginary signals；
The error of the error of the solid part signal and the imaginary signals is synthesized, answering for the polarization state complex signal is obtained
Number error；
According to the tap coefficient of the last moment and each polarization state complex signal, the tap at the current time is calculated
Coefficient.
Adaptive blind polarization demultiplexing processing method provided by the invention is needed when calculating the tap coefficient of subsequent time
Each polarization state complex signal in the output signal of the last moment of the subsequent time is carried out being decomposed into solid part signal and imaginary part
Signal is coordinately transformed the signal after decomposition, seeks error, and error merges to obtain plural error, further according to plural error meter
The tap coefficient of the subsequent time is calculated, can be utilized when auto-adaptive fir filter receives the input signal of the subsequent time
The tap coefficient of the subsequent time carries out the processing work of polarization demultiplexing.Due to the present invention is by decomposing, coordinate transform can be with
Therefore signal tune can not known as constant modulus signals by transforming to the multimode signal of different modulating format on the same circle
Polarization demultiplexing is carried out under the premise of format processed to it, therefore method provided by the invention is unrelated with modulation format, is suitable for net
Modulation format is further flexible in network, flow dynamic, multiplicity and uncertain application scenarios further.
Detailed description of the invention
In order to illustrate more clearly of the embodiment of the present disclosure or technical solution in the prior art, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Disclosed some embodiments for those of ordinary skill in the art without creative efforts, can be with
Other attached drawings are obtained according to these figures.
Fig. 1 shows the flow diagram of adaptive blind polarization demultiplexing processing method in one embodiment of the invention.
Specific embodiment
Below in conjunction with the attached drawing in the embodiment of the present disclosure, the technical solution in the embodiment of the present disclosure is carried out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the disclosure, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment belongs to the range of disclosure protection.
The present invention provides a kind of adaptive blind polarization demultiplexing processing method, as shown in Figure 1, this method includes：
Determine auto-adaptive fir filter in the tap coefficient at current time；
According to the tap coefficient at the current time, the auto-adaptive fir filter to the input signal at current time into
The processing of row polarization demultiplexing, and polarization demultiplexing treated signal is exported；
Wherein, tap coefficient of the adaptive butterfly FIR filter of the determination at current time, including：
To each polarization state complex signal of the adaptive butterfly FIR filter in the output signal of last moment into
Row decomposes, and obtains the solid part signal and imaginary signals of the polarization state complex signal；
The solid part signal and imaginary signals of the polarization state complex signal are coordinately transformed respectively；
According to the real number signal obtained after the reference modulus value of the solid part signal of the polarization state complex signal and coordinate transform, meter
The error of the solid part signal is calculated, and according to the reality after the reference modulus value and coordinate transform of the imaginary signals of the polarization state complex signal
Number signal calculates the error of the imaginary signals；
The error of the error of the solid part signal and the imaginary signals is synthesized, answering for the polarization state complex signal is obtained
Number error；
According to the tap coefficient of the last moment and each polarization state complex signal, the tap at the current time is calculated
Coefficient.
It will be appreciated that auto-adaptive fir filter has limit for length's unit impulse response filter for adaptive.
It will be appreciated that auto-adaptive fir filter used when carrying out polarization demultiplexing to the input signal at current time
Current tap coefficient is to carry out certain processing to the output signal of last moment auto-adaptive fir filter to obtain, for example, right
The tap coefficient that the input signal of 2nd moment auto-adaptive fir filter use when polarization demultiplexing is according to the 1st moment
What the output signal of auto-adaptive fir filter was handled, to the input signal of the 3rd moment auto-adaptive fir filter
The tap coefficient use when polarization demultiplexing is handled according to the output signal of the 2nd moment auto-adaptive fir filter
It obtains, and so on ....Wherein, it is adopted when the input signal of the 1st moment auto-adaptive fir filter carries out polarization demultiplexing
Initial tap coefficient values are self-settinies.As it can be seen that the tap coefficient at each moment is the output according to previous moment below
What signal obtained.Tap coefficient is constantly updated in this way, it is made to meet the requirement of polarization demultiplexing.
Adaptive blind polarization demultiplexing processing method provided by the invention is needed when calculating the tap coefficient of subsequent time
Each polarization state complex signal in the output signal of the last moment of the subsequent time is carried out being decomposed into solid part signal and imaginary part
Signal is coordinately transformed the signal after decomposition, seeks error, and error merges to obtain plural error, further according to plural error meter
The tap coefficient of the subsequent time is calculated, can be utilized when auto-adaptive fir filter receives the input signal of the subsequent time
The tap coefficient of the subsequent time carries out the processing work of polarization demultiplexing.Due to the present invention by complex signal be decomposed into real part and
Then imaginary part is coordinately transformed, therefore the signal of different modulating format can be made (for example, PM-BPSK, PM-QPSK, PM-
4PAM, PM-8PAM, PM-16QAM, PM-32QAM, PM-64QAM) it transforms on the same circle, therefore, it can not know letter
Polarization demultiplexing is carried out to it under the premise of number modulation format, therefore method provided by the invention is unrelated with modulation format, be applicable in
In network flow dynamic, multiplicity and uncertain application scenarios further.Moreover, method provided by the invention is relative to traditional
RDA has the advantages that simple, anti-noise ability is strong, strong robustness, relative to traditional DD-LMS, have it is easy to accomplish, without doing
Wave carries the advantages of phase recovery.Therefore method provided by the invention is applicable not only to coherent optical communication system, to entire optic communication
System or other specific application scenes are also suitable.
It will be appreciated that auto-adaptive fir filter may include x-polarisation state plural number letter in the output signal at each moment
Number and y-polarisation state complex signal.
In the specific implementation, the solid part signal of x-polarisation state complex signal can be coordinately transformed using following formula：
x′i(k)=xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]} (1)
Wherein, xiIt (k) is the reality of auto-adaptive fir filter x-polarisation state complex signal in the output signal at k moment
Portion's signal, x 'i(k) for solid part signal xi(k) real number signal obtained after coordinate transform.
Likewise, can be coordinately transformed using following formula to the imaginary signals of x-polarisation state complex signal：
x′q(k)=xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]} (2)
Wherein, xqIt (k) is the void of auto-adaptive fir filter x-polarisation state complex signal in the output signal at k moment
Portion's signal, x 'q(k) for imaginary signals xq(k) real number signal obtained after being coordinately transformed.
Likewise, can be coordinately transformed using following formula to the solid part signal of y-polarisation state complex signal：
y′i(k)=yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]} (3)
Wherein, yiIt (k) is the reality of auto-adaptive fir filter y-polarisation state complex signal in the output signal at k moment
Portion's signal, y 'i(k) for solid part signal yi(k) real number signal obtained after being coordinately transformed.
Likewise, can be coordinately transformed using following formula to the imaginary signals of y-polarisation state complex signal：
y′q(k)=yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]} (4)
Wherein, yqIt (k) is the void of auto-adaptive fir filter y-polarisation state complex signal in the output signal at k moment
Portion's signal, y 'q(k) for imaginary signals yq(k) real number signal obtained after being coordinately transformed.
It is understood that sign () is to take sign function.
In the specific implementation, the error of the solid part signal of x-polarisation state complex signal can be calculated using following formula：
εxi(k)=x 'i(k)(|Rx′i|2-|x′i(k)|2) (5)
Wherein, Rx′iFor solid part signal xi(k) reference modulus value, εxiIt (k) is solid part signal xi(k) error.
Likewise, the error of the imaginary signals of x-polarisation state complex signal can be calculated using following formula：
εxq(k)=x 'q(k)(|Rx′q|2-|x′q(k)|2) (6)
Wherein, Rx′qFor imaginary signals xq(k) reference modulus value, εxqIt (k) is imaginary signals xq(k) error.
Likewise, the error of the solid part signal of y-polarisation state complex signal can be calculated using following formula：
εyi(k)=y 'i(k)(|Ry′i|2-|y′i(k)|2) (7)
Wherein, Ry′iFor solid part signal yi(k) reference modulus value, εyiIt (k) is solid part signal yi(k) error.
Likewise, the error of y-polarisation state complex signal imaginary signals can be calculated using following formula：
εyq(k)=y 'q(k)(|Ry′q|2-|y′q(k)|2) (8)
Wherein, Ry′qFor imaginary signals yq(k) reference modulus value, εyqIt (k) is imaginary signals yq(k) error.
R in above formulax′i、Rx′q、Ry′i、Ry′qCalculation method it is as follows：
Wherein it is possible to calculate solid part signal x using following formulai(k) reference modulus value Rx′i：
In formula, sxiIt (k) is the real part of complex signal under kth moment x-polarisation state both ideal constellation；And/or
Wherein, imaginary signals x is calculated using following formulaq(k) reference modulus value Rx′q：
In formula, sxqIt (k) is the imaginary part of complex signal under kth moment x-polarisation state both ideal constellation；And/or
Wherein, solid part signal y is calculated using following formulai(k) reference modulus value Ry′i：
In formula, syiIt (k) is the real part of complex signal under kth moment y-polarisation state both ideal constellation；And/or
Wherein, imaginary signals y is calculated using following formulaq(k) reference modulus value Ry′q：
In formula, syqIt (k) is the imaginary part of complex signal under kth moment y-polarisation state both ideal constellation.
In the specific implementation, the plural error of x-polarisation state can be calculated using following formula：
εx(k)=εxi(k)+jεxq(k) (13)
Wherein, εxIt (k) is the plural error of x-polarisation state.
Likewise, the plural error of y-polarisation state can be calculated using following formula：
εy(k)=εyi(k)+jεyq(k) (14)
Wherein, εyIt (k) is the plural error of y-polarisation state.
In the specific implementation, the tap coefficient at+1 moment of kth can be calculated using following formula：
Hk+1,xx(m)=Hk,xx(m)+μ·εx(k)·xin(k-m)* (15)
Hk+1,xy(m)=Hk,xy(m)+μ·εx(k)·yin(k-m)* (16)
Hk+1,yx(m)=Hk,yx(m)+μ·εy(k)·xin(k-m)* (17)
Hk+1,yy(m)=Hk,yy(m)+μ·εy(k)·yin(k-m)* (18)
Wherein, μ is the step-length that tap updates, Hk,xx(m)、Hk,xy(m)、Hk,yx(m)、Hk,yy(m) at m-th of the kth moment
The coefficient of tap, Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yy(m) in the coefficient of m-th of tap of+1 moment of kth, xin
(k-m)*For in the conjugation of the input signal of kth moment x-polarisation state, yin(k-m)*For in the input signal of kth moment y-polarisation state
Conjugation.
In the specific implementation, polarization demultiplexing processing can be carried out using input signal of the following formula to+1 moment of kth：
Wherein, xin(k+1-m) in the input signal of+1 moment of kth x-polarisation state, yin(k+1-m) in+1 moment of kth y
The input signal of polarization state, x (k+1) are obtained after carrying out polarization demultiplexing processing to+1 moment of kth x-polarisation state input signal
Signal, y (k+1) are the signal for obtain after polarization demultiplexing processing to+1 moment of kth y-polarisation state input signal, and 2N+1 is
Tap number.
It will be appreciated that the value range of m is [- N, N].
In the specific implementation, can also before the signal is input to filter to signal carry out quadrature imbalance compensation,
The pretreatments such as fixed dispersion compensation is synchronous with clock, are then input to filter for pretreated signal and handle, and is defeated
Out, offset estimation can also be carried out and the processing such as compensation, skew estimate and compensation, symbol reflection are penetrated certainly for output signal,
Finally obtain originator original bit sequence.
It is experimentally confirmed, solution vibration processing partially is carried out to signal using method provided by the invention, relative to traditional
The methods of CMA, RDE have the characteristics that superior performance.
Based on identical inventive concept, the present invention also provides a kind of adaptive blind polarization demultiplexing processing unit, the devices
Including：
Tap coefficient determining module, for determining auto-adaptive fir filter in the tap coefficient at current time；
The auto-adaptive fir filter believes the input at current time for the tap coefficient according to the current time
Number polarization demultiplexing processing is carried out, and polarization demultiplexing treated signal is exported；
Wherein, the tap coefficient determining module is specifically used for：
Each polarization state complex signal of the auto-adaptive fir filter in the output signal of last moment is divided
Solution, obtains the solid part signal and imaginary signals of the polarization state complex signal；
The solid part signal and imaginary signals of the polarization state complex signal are coordinately transformed respectively；
According to the real number signal obtained after the reference modulus value of the solid part signal of the polarization state complex signal and coordinate transform, meter
The error of the solid part signal is calculated, and according to the reality after the reference modulus value and coordinate transform of the imaginary signals of the polarization state complex signal
Number signal calculates the error of the imaginary signals；
The error of the error of the solid part signal and the imaginary signals is synthesized, answering for the polarization state complex signal is obtained
Number error；
According to the tap coefficient of the last moment and each polarization state complex signal, the tap at the current time is calculated
Coefficient.
Those of ordinary skill in the art will appreciate that：Realize that all or part of step of above method embodiment can lead to
The relevant hardware of program instruction is crossed to complete, program above-mentioned can store in computer-readable storage medium, the journey
Sequence when being executed, executes step including the steps of the foregoing method embodiments.
In specification of the invention, numerous specific details are set forth.It is to be appreciated, however, that the embodiment of the present invention can be with
It practices without these specific details.In some instances, well known method, structure and skill is not been shown in detail
Art, so as not to obscure the understanding of this specification.
The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although with reference to the foregoing embodiments
Invention is explained in detail, those skilled in the art should understand that；It still can be to aforementioned each implementation
Technical solution documented by example is modified or equivalent replacement of some of the technical features；And these modification or
Replacement, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution.
Claims (7)
1. a kind of adaptive blind polarization demultiplexing processing method, which is characterized in that including：
Determine adaptive butterfly FIR filter in the tap coefficient at current time；
According to the tap coefficient at the current time, the adaptive butterfly FIR filter to the input signal at current time into
The processing of row polarization demultiplexing, and polarization demultiplexing treated signal is exported；
Wherein, tap coefficient of the adaptive butterfly FIR filter of the determination at current time, including：
Each polarization state complex signal of the adaptive butterfly FIR filter in the output signal of last moment is divided
Solution, obtains the solid part signal and imaginary signals of the polarization state complex signal；
The solid part signal and imaginary signals of the polarization state complex signal are coordinately transformed respectively；
According to the real number signal obtained after the reference modulus value of the solid part signal of the polarization state complex signal and coordinate transform, calculating should
The error of solid part signal, and believed according to the real number after the reference modulus value and coordinate transform of the imaginary signals of the polarization state complex signal
Number calculate the error of the imaginary signals；
The error of the error of the solid part signal and the imaginary signals is synthesized, obtains the plural number of the polarization state complex signal accidentally
Difference；
According to the tap coefficient of the last moment and each polarization state complex signal, the tap system at the current time is calculated
Number；
The auto-adaptive fir filter includes x-polarisation state complex signal and y-polarisation state plural number in the output signal at each moment
Signal；It is corresponding,
It is coordinately transformed using solid part signal of the following formula to x-polarisation state complex signal：
x′i(k)=xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
Wherein, xiIt (k) is the reality of adaptive butterfly FIR filter x-polarisation state complex signal in the output signal at k moment
Portion's signal, x 'i(k) for solid part signal xi(k) real number signal obtained after coordinate transform；And/or
It is coordinately transformed using imaginary signals of the following formula to x-polarisation state complex signal：
x′q(k)=xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
Wherein, xq(k) believe for the imaginary part of auto-adaptive fir filter x-polarisation state complex signal in the output signal at k moment
Number, x 'q(k) for imaginary signals xq(k) real number signal obtained after being coordinately transformed；And/or
It is coordinately transformed using solid part signal of the following formula to y-polarisation state complex signal：
y′i(k)=yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
Wherein, yi(k) believe for the real part of auto-adaptive fir filter y-polarisation state complex signal in the output signal at k moment
Number, y 'i(k) for solid part signal yi(k) real number signal obtained after being coordinately transformed；And/or
It is coordinately transformed using imaginary signals of the following formula to y-polarisation state complex signal：
y′q(k)=yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
Wherein, yq(k) believe for the imaginary part of auto-adaptive fir filter y-polarisation state complex signal in the output signal at k moment
Number, y 'q(k) for imaginary signals yq(k) real number signal obtained after being coordinately transformed.
2. the method according to claim 1, wherein
The error of the solid part signal of x-polarisation state complex signal is calculated using following formula：
εxi(k)=x 'i(k)(|Rx′i|2-|x′i(k)|2)
Wherein, Rx′iFor solid part signal xi(k) reference modulus value, εxiIt (k) is solid part signal xi(k) error；And/or
The error of the imaginary signals of x-polarisation state complex signal is calculated using following formula：
εxq(k)=x 'q(k)(|Rx′q|2-|x′q(k)|2)
Wherein, Rx′qFor imaginary signals xq(k) reference modulus value, εxqIt (k) is imaginary signals xq(k) error；And/or
The error of the solid part signal of y-polarisation state complex signal is calculated using following formula：
εyi(k)=y 'i(k)(|Ry′i|2-|y′i(k)|2)
Wherein, Ry′iFor solid part signal yi(k) reference modulus value, εyiIt (k) is solid part signal yi(k) error；And/or
The error of y-polarisation state complex signal imaginary signals is calculated using following formula：
εyq(k)=y 'q(k)(|Ry′q|2-|y′q(k)|2)
Wherein, Ry′qFor imaginary signals yq(k) reference modulus value, εyqIt (k) is imaginary signals yq(k) error.
3. according to the method described in claim 2, it is characterized in that,
The plural error of x-polarisation state is calculated using following formula：
εx(k)=εxi(k)+jεxq(k)
Wherein, εxIt (k) is the plural error of x-polarisation state；And/or
The plural error of y-polarisation state is calculated using following formula：
εy(k)=εyi(k)+jεyq(k)
Wherein, εyIt (k) is the plural error of y-polarisation state.
4. according to the method described in claim 3, it is characterized in that, calculating the tap coefficient at+1 moment of kth using following formula：
Hk+1,xx(m)=Hk,xx(m)+μ·εx(k)·xin(k-m)*
Hk+1,xy(m)=Hk,xy(m)+μ·εx(k)·yin(k-m)*
Hk+1,yx(m)=Hk,yx(m)+μ·εy(k)·xin(k-m)*
Hk+1,yy(m)=Hk,yy(m)+μ·εy(k)·yin(k-m)*
Wherein, μ is the step-length that tap updates, Hk,xx(m)、Hk,xy(m)、Hk,yx(m)、Hk,yy(m) in m-th of tap of kth moment
Coefficient, Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yy(m) in the coefficient of m-th of tap of+1 moment of kth, xin(k-
m)*For in the conjugation of the input signal of kth moment x-polarisation state, yin(k-m)*For in the input signal of kth moment y-polarisation state
Conjugation.
5. according to the method described in claim 4, it is characterized in that, being carried out using input signal of the following formula to+1 moment of kth inclined
Vibration demultiplexing process：
Wherein, xin(k+1-m) in the input signal of+1 moment of kth x-polarisation state, yin(k+1-m) in+1 moment of kth y-polarisation
The input signal of state, x (k+1) are inclined to the x obtained after the progress polarization demultiplexing processing of two polarization state input signal of+1 moment of kth
Polarization state signal, y (k+1) are the y-polarisation state for obtain after polarization demultiplexing processing to two polarization state input signal of+1 moment of kth
Signal, 2N+1 are tap number.
6. according to the method described in claim 2, it is characterized in that,
Solid part signal x is calculated using following formulai(k) reference modulus value Rx′i：
Wherein, sxiIt (k) is the real part of complex signal under kth moment x-polarisation state both ideal constellation；And/or
Imaginary signals x is calculated using following formulaq(k) reference modulus value Rx′q：
Wherein, sxqIt (k) is the imaginary part of complex signal under kth moment x-polarisation state both ideal constellation；And/or
Solid part signal y is calculated using following formulai(k) reference modulus value Ry′i：
Wherein, syiIt (k) is the real part of complex signal under kth moment y-polarisation state both ideal constellation；And/or
Imaginary signals y is calculated using following formulaq(k) reference modulus value Ry′q：
Wherein, syqIt (k) is the imaginary part of complex signal under kth moment y-polarisation state both ideal constellation.
7. a kind of adaptive blind polarization demultiplexing processing unit, which is characterized in that including：
Tap coefficient determining module, for determining auto-adaptive fir filter in the tap coefficient at current time；
The auto-adaptive fir filter, for the tap coefficient according to the current time, to the input signal at current time into
The processing of row polarization demultiplexing, and polarization demultiplexing treated signal is exported；
Wherein, the tap coefficient determining module is specifically used for：
Each polarization state complex signal of the auto-adaptive fir filter in the output signal of last moment is decomposed,
Obtain the solid part signal and imaginary signals of the polarization state complex signal；
The solid part signal and imaginary signals of the polarization state complex signal are coordinately transformed respectively；
According to the real number signal obtained after the reference modulus value of the solid part signal of the polarization state complex signal and coordinate transform, calculating should
The error of solid part signal, and believed according to the real number after the reference modulus value and coordinate transform of the imaginary signals of the polarization state complex signal
Number calculate the error of the imaginary signals；
The error of the error of the solid part signal and the imaginary signals is synthesized, obtains the plural number of the polarization state complex signal accidentally
Difference；
According to the tap coefficient of the last moment and each polarization state complex signal, the tap system at the current time is calculated
Number；
The auto-adaptive fir filter includes x-polarisation state complex signal and y-polarisation state plural number in the output signal at each moment
Signal；It is corresponding,
It is coordinately transformed using solid part signal of the following formula to x-polarisation state complex signal：
x′i(k)=xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
Wherein, xiIt (k) is the reality of adaptive butterfly FIR filter x-polarisation state complex signal in the output signal at k moment
Portion's signal, x 'i(k) for solid part signal xi(k) real number signal obtained after coordinate transform；And/or
It is coordinately transformed using imaginary signals of the following formula to x-polarisation state complex signal：
x′q(k)=xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
Wherein, xq(k) believe for the imaginary part of auto-adaptive fir filter x-polarisation state complex signal in the output signal at k moment
Number, x 'q(k) for imaginary signals xq(k) real number signal obtained after being coordinately transformed；And/or
It is coordinately transformed using solid part signal of the following formula to y-polarisation state complex signal：
y′i(k)=yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
Wherein, yi(k) believe for the real part of auto-adaptive fir filter y-polarisation state complex signal in the output signal at k moment
Number, y 'i(k) for solid part signal yi(k) real number signal obtained after being coordinately transformed；And/or
It is coordinately transformed using imaginary signals of the following formula to y-polarisation state complex signal：
y′q(k)=yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
Wherein, yq(k) believe for the imaginary part of auto-adaptive fir filter y-polarisation state complex signal in the output signal at k moment
Number, y 'q(k) for imaginary signals yq(k) real number signal obtained after being coordinately transformed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610650339.5A CN106301593B (en) | 2016-08-09 | 2016-08-09 | Adaptive blind polarization demultiplexing treating method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610650339.5A CN106301593B (en) | 2016-08-09 | 2016-08-09 | Adaptive blind polarization demultiplexing treating method and apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106301593A CN106301593A (en) | 2017-01-04 |
| CN106301593B true CN106301593B (en) | 2018-11-16 |
Family
ID=57667631
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610650339.5A Active CN106301593B (en) | 2016-08-09 | 2016-08-09 | Adaptive blind polarization demultiplexing treating method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106301593B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110011733B (en) * | 2019-03-25 | 2020-10-16 | 华中科技大学 | Method and system for depolarization multiplexing based on momentum factor |
| CN111884960B (en) * | 2020-09-28 | 2020-12-25 | 烽火通信科技股份有限公司 | Polarization demultiplexing signal processing blind equalization method and polarization demultiplexing device |
| CN113708845B (en) * | 2021-09-08 | 2022-05-06 | 北京科技大学 | Equalization demultiplexing method, apparatus, electronic device and computer readable storage medium |
| CN114337843B (en) * | 2021-12-23 | 2023-11-28 | 武汉邮电科学研究院有限公司 | Self-adaptive channel equalization method and system for coherent optical communication |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103283201A (en) * | 2011-03-31 | 2013-09-04 | 富士通株式会社 | Method and device for setting coefficients of butterfly filter, receiver and receiving method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101895499B (en) * | 2010-06-30 | 2014-11-19 | 北京邮电大学 | Method for demultiplexing polarization by using constant rotation sign training sequence |
| US9154231B2 (en) * | 2013-01-17 | 2015-10-06 | Alcatel Lucent | Generation of an optical local-oscillator signal for a coherent-detection scheme |
| CN105703838B (en) * | 2016-01-26 | 2018-07-03 | 哈尔滨工业大学深圳研究生院 | Coherent optical heterodyne communicatio method for dynamically balancing based on butterfly linear kalman filter |
-
2016
- 2016-08-09 CN CN201610650339.5A patent/CN106301593B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103283201A (en) * | 2011-03-31 | 2013-09-04 | 富士通株式会社 | Method and device for setting coefficients of butterfly filter, receiver and receiving method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106301593A (en) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106301593B (en) | Adaptive blind polarization demultiplexing treating method and apparatus | |
| US8260156B2 (en) | Adaptive crossing frequency domain equalization (FDE) in digital PolMux coherent systems | |
| US9467246B2 (en) | Clock recovery method and clock recovery arrangement for coherent polarization multiplex receivers | |
| CN101552640B (en) | Device and method for altering filter coefficients | |
| US8538278B2 (en) | Superimposed training and digital filtering coherent optical receivers | |
| CN103338171B (en) | A kind of receiving terminal equalization methods based on channel estimation in frequency domain and system | |
| US9252988B2 (en) | System and methods for adaptive equalization for optical modulation formats | |
| JP5120507B2 (en) | Optical receiver, polarization separation device, and polarization separation method | |
| EP2436127A1 (en) | Method for carrier frequency recovery and optical intradyne coherent receiver | |
| US10171177B2 (en) | Digital signal processor, digital optical receiver using the same, and digital signal processing method | |
| CN104168069B (en) | The modulating system of phase controlled rotation in the coherent light communications field and modulator approach | |
| CN106330322B (en) | Frequency deviation skew processing method and processing device | |
| CN105530050A (en) | Joint processing method and device for equilibrium and polarization demultiplexing and phase offset estimation and compensation | |
| CN112929124B (en) | Method and system for realizing depolarization and equalization of polarization multiplexing signal | |
| CN102724151B (en) | Communication system, communication method, multichannel adaptive equalizer and OFDM (orthogonal frequency division multiplexing) signal demultiplexing method | |
| CN107925485A (en) | Relevant optical pickup apparatus | |
| CN106789793A (en) | For the LMS methods without frequency difference difference feedback that adaptive equalizer coefficients update | |
| CN104767702B (en) | The phase hit removing method and device of single carrier coherent optical communication system | |
| CN108199777A (en) | Coherent optical communication system blind balance method based on probability density function fitting and fuzzy logic | |
| JP6378654B2 (en) | Spatial multiplexed optical transmission system and spatial multiplexed optical transmission method | |
| CN105553529B (en) | A kind of less fundamental mode optical fibre Transmission system and its digital signal recovery method | |
| CN110224761A (en) | A kind of method and system of rapid solving polarization spin matrix and equalizer coefficients | |
| CN105282072B (en) | Frequency deviation estimating method and device in optical transport network | |
| JP6387832B2 (en) | Polarization demultiplexing optical communication receiver, polarization demultiplexing optical communication system, and polarization demultiplexing optical communication method | |
| CN108173794A (en) | A kind of novel LMS methods and the real number adaptive equalizer using this method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |