CN106301593A - Adaptive blind polarization demultiplexing treating method and apparatus - Google Patents
Adaptive blind polarization demultiplexing treating method and apparatus Download PDFInfo
- Publication number
- CN106301593A CN106301593A CN201610650339.5A CN201610650339A CN106301593A CN 106301593 A CN106301593 A CN 106301593A CN 201610650339 A CN201610650339 A CN 201610650339A CN 106301593 A CN106301593 A CN 106301593A
- Authority
- CN
- China
- Prior art keywords
- signal
- centerdot
- moment
- error
- 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.)
- Granted
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, 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
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)=xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
Wherein, xiK () 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 xiThe 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)=xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
Wherein, xqK () 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 xqK 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)=yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
Wherein, yiK () 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 yiK 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)=yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
Wherein, yqK () 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 yqK 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)(|Rx′i|2-|x′i(k)|2)
Wherein, Rx′iFor solid part signal xiThe reference modulus value of (k), εxiK () is solid part signal xiThe 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)(|Rx′q|2-|x′q(k)|2)
Wherein, Rx′qFor imaginary signals xqThe reference modulus value of (k), εxqK () is imaginary signals xqThe 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)(|Ry′i|2-|y′i(k)|2)
Wherein, Ry′iFor solid part signal yiThe reference modulus value of (k), εyiK () is solid part signal yiThe error of (k);And/or
The error of employing following formula calculating y-polarisation state complex signal imaginary signals:
εyq(k)=y 'q(k)(|Ry′q|2-|y′q(k)|2)
Wherein, Ry′qFor imaginary signals yqThe reference modulus value of (k), εyqK () is imaginary signals yqThe 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:
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,yyM () is in kth moment m-th
The coefficient of tap, Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yyM () is the coefficient in kth+1 moment m-th tap, xin
(k-m)*For the conjugation of the input signal in kth moment x-polarisation state, yin(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:
Wherein, xin(k+1-m) it is the input signal in kth+1 moment x-polarisation state, yin(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 xiReference modulus value R of (k)x′i:
Wherein, sxiK () 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 xqReference modulus value R of (k)x′q:
Wherein, sxqK () 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 yiReference modulus value R of (k)y′i:
Wherein, syiK () 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 yqReference modulus value R of (k)y′q:
Wherein, syqK () 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)=xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]} (1)
Wherein, xiK () 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 xiThe 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)=xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]} (2)
Wherein, xqK () 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 xqK 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)=yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]} (3)
Wherein, yiK () 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 yiK 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)=yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]} (4)
Wherein, yqK () 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 yqK 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)(|Rx′i|2-|x′i(k)|2) (5)
Wherein, Rx′iFor solid part signal xiThe reference modulus value of (k), εxiK () is solid part signal xiThe 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)(|Rx′q|2-|x′q(k)|2) (6)
Wherein, Rx′qFor imaginary signals xqThe reference modulus value of (k), εxqK () is imaginary signals xqThe 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)(|Ry′i|2-|y′i(k)|2) (7)
Wherein, Ry′iFor solid part signal yiThe reference modulus value of (k), εyiK () is solid part signal yiThe 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)(|Ry′q|2-|y′q(k)|2) (8)
Wherein, Ry′qFor imaginary signals yqThe reference modulus value of (k), εyqK () is imaginary signals yqThe error of (k).
R in above formulax′i、Rx′q、Ry′i、Ry′qComputational methods as follows:
Wherein it is possible to use following formula to calculate solid part signal xiReference modulus value R of (k)x′i:
In formula, sxiK () 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 xqReference modulus value R of (k)x′q:
In formula, sxqK () 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 yiReference modulus value R of (k)y′i:
In formula, syiK () 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 yqReference modulus value R of (k)y′q:
In formula, syqK () 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:
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,yyM () is in kth moment m-th
The coefficient of tap, Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yyM () is the coefficient in kth+1 moment m-th tap, xin
(k-m)*For the conjugation of the input signal in kth moment x-polarisation state, yin(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:
Wherein, xin(k+1-m) it is the input signal in kth+1 moment x-polarisation state, yin(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)=xi(k)-4·sign[xi(k)]-2·sign{xi(k)-4·sign[xi(k)]}
Wherein, xiK () 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 xiThe 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)=xq(k)-4·sign[xq(k)]-2·sign{xq(k)-4·sign[xq(k)]}
Wherein, xqK () 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 xqK 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)=yi(k)-4·sign[yi(k)]-2·sign{yi(k)-4·sign[yi(k)]}
Wherein, yiK () 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 yiK 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)=yq(k)-4·sign[yq(k)]-2·sign{yq(k)-4·sign[yq(k)]}
Wherein, yqK () 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 yqK 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)(|Rx′i|2-|x′i(k)|2)
Wherein, Rx′iFor solid part signal xiThe reference modulus value of (k), εxiK () is solid part signal xiThe 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)(|Rx′q|2-|x′q(k)|2)
Wherein, Rx′qFor imaginary signals xqThe reference modulus value of (k), εxqK () is imaginary signals xqThe 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)(|Ry′i|2-|y′i(k)|2)
Wherein, Ry′iFor solid part signal yiThe reference modulus value of (k), εyiK () is solid part signal yiThe error of (k);And/or
The error of employing following formula calculating y-polarisation state complex signal imaginary signals:
εyq(k)=y 'q(k)(|Ry′q|2-|y′q(k)|2)
Wherein, Ry′qFor imaginary signals yqThe reference modulus value of (k), εyqK () is imaginary signals yqThe 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:
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,yyM () is in the m-th tap of kth moment
Coefficient, Hk+1,xx(m)、Hk+1,xy(m)、Hk+1,yx(m)、Hk+1,yyM () is the coefficient in kth+1 moment m-th tap, xin(k-
m)*For the conjugation of the input signal in kth moment x-polarisation state, yin(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:
Wherein, xin(k+1-m) it is the input signal in kth+1 moment x-polarisation state, yin(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 xiReference modulus value R of (k)x′i:
Wherein, sxiK () 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 xqReference modulus value R of (k)x′q:
Wherein, sxqK () 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 yiReference modulus value R of (k)y′i:
Wherein, syiK () 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 yqReference modulus value R of (k)y′q:
Wherein, syqK () 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.
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 true CN106301593A (en) | 2017-01-04 |
CN106301593B 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) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110011733A (en) * | 2019-03-25 | 2019-07-12 | 华中科技大学 | A kind of depolarization multiplexing method and system based on factor of momentum |
CN111884960A (en) * | 2020-09-28 | 2020-11-03 | 烽火通信科技股份有限公司 | Polarization demultiplexing signal processing blind equalization method and polarization demultiplexing device |
CN113708845A (en) * | 2021-09-08 | 2021-11-26 | 北京科技大学 | Equalization demultiplexing method, apparatus, electronic device and computer readable storage medium |
CN114337843A (en) * | 2021-12-23 | 2022-04-12 | 武汉邮电科学研究院有限公司 | Self-adaptive channel equalization method and system for coherent optical communication |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101895499A (en) * | 2010-06-30 | 2010-11-24 | 北京邮电大学 | Method for demultiplexing polarization by using constant rotation sign training sequence |
CN103283201A (en) * | 2011-03-31 | 2013-09-04 | 富士通株式会社 | Method and device for setting coefficients of butterfly filter, receiver and receiving method |
US20140199063A1 (en) * | 2013-01-17 | 2014-07-17 | Alcatel-Lucent Usa Inc. | Generation of an optical local-oscillator signal for a coherent-detection scheme |
CN105703838A (en) * | 2016-01-26 | 2016-06-22 | 哈尔滨工业大学深圳研究生院 | A coherent light receiver dynamic balancing method based on a butterfly linear Kalman filter |
-
2016
- 2016-08-09 CN CN201610650339.5A patent/CN106301593B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101895499A (en) * | 2010-06-30 | 2010-11-24 | 北京邮电大学 | Method for demultiplexing polarization by using constant rotation sign training sequence |
CN103283201A (en) * | 2011-03-31 | 2013-09-04 | 富士通株式会社 | Method and device for setting coefficients of butterfly filter, receiver and receiving method |
US20140199063A1 (en) * | 2013-01-17 | 2014-07-17 | Alcatel-Lucent Usa Inc. | Generation of an optical local-oscillator signal for a coherent-detection scheme |
CN105703838A (en) * | 2016-01-26 | 2016-06-22 | 哈尔滨工业大学深圳研究生院 | A coherent light receiver dynamic balancing method based on a butterfly linear Kalman filter |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110011733A (en) * | 2019-03-25 | 2019-07-12 | 华中科技大学 | A kind of depolarization multiplexing method and system based on factor of momentum |
CN111884960A (en) * | 2020-09-28 | 2020-11-03 | 烽火通信科技股份有限公司 | Polarization demultiplexing signal processing blind equalization method and polarization demultiplexing device |
CN111884960B (en) * | 2020-09-28 | 2020-12-25 | 烽火通信科技股份有限公司 | Polarization demultiplexing signal processing blind equalization method and polarization demultiplexing device |
CN113708845A (en) * | 2021-09-08 | 2021-11-26 | 北京科技大学 | Equalization demultiplexing method, apparatus, electronic device and computer readable storage medium |
CN113708845B (en) * | 2021-09-08 | 2022-05-06 | 北京科技大学 | Equalization demultiplexing method, apparatus, electronic device and computer readable storage medium |
CN114337843A (en) * | 2021-12-23 | 2022-04-12 | 武汉邮电科学研究院有限公司 | Self-adaptive channel equalization method and system for coherent optical communication |
CN114337843B (en) * | 2021-12-23 | 2023-11-28 | 武汉邮电科学研究院有限公司 | Self-adaptive channel equalization method and system for coherent optical communication |
Also Published As
Publication number | Publication date |
---|---|
CN106301593B (en) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8260156B2 (en) | Adaptive crossing frequency domain equalization (FDE) in digital PolMux coherent systems | |
CN106301593A (en) | Adaptive blind polarization demultiplexing treating method and apparatus | |
CN101599801B (en) | Filter coefficient regulator and method thereof | |
CN103609040A (en) | Optical receiver for multimode communications | |
CN103338171B (en) | A kind of receiving terminal equalization methods based on channel estimation in frequency domain and system | |
US20110194855A1 (en) | Superimposed training and digital filtering coherent optical receivers | |
CN103812806B (en) | A kind of channel equalization method estimated based on time domain channel and system | |
CN103329497A (en) | Digital demodulator architecture | |
EP2436127A1 (en) | Method for carrier frequency recovery and optical intradyne coherent receiver | |
CN107204945B (en) | Frequency difference estimation device, channel interval estimation device, method and system | |
CN102420658B (en) | Average length detection device and method as well as light coherent receiver | |
CN104168069B (en) | The modulating system of phase controlled rotation in the coherent light communications field and modulator approach | |
CN105432029A (en) | Method, device and communication system for reducing optical transmission impairments | |
CN106992835A (en) | The construction method of mode division multiplexing optical fiber telecommunications system and the optical fiber telecommunications system of structure | |
CN110291728B (en) | Method and apparatus for dual polarized optical communications | |
CN103684601A (en) | Coefficient determining device, balancer, receiver and transmitter | |
CN105530050A (en) | Joint processing method and device for equilibrium and polarization demultiplexing and phase offset estimation and compensation | |
CN102724151B (en) | Communication system, communication method, multichannel adaptive equalizer and OFDM (orthogonal frequency division multiplexing) signal demultiplexing method | |
CN106789793A (en) | For the LMS methods without frequency difference difference feedback that adaptive equalizer coefficients update | |
CN106998229A (en) | It is a kind of based on variable step without constraint FD LMS mode division multiplexing system Deplexing method | |
CN102160336A (en) | Simulation device and simulation method | |
CN105553529B (en) | A kind of less fundamental mode optical fibre Transmission system and its digital signal recovery method | |
CN101764649A (en) | MGDM multimode optical fiber communication system and method based on ICA signal separation algorithm | |
Qin et al. | An improved CMA for dispersion compensation in 100 Gbps DP-QPSK optical signal transmission system | |
CN103179060A (en) | Balancer coefficient updating device and 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 |