CN110034827A - A kind of depolarization multiplexing method and system based on reverse observation error - Google Patents
A kind of depolarization multiplexing method and system based on reverse observation error Download PDFInfo
- Publication number
- CN110034827A CN110034827A CN201910227443.7A CN201910227443A CN110034827A CN 110034827 A CN110034827 A CN 110034827A CN 201910227443 A CN201910227443 A CN 201910227443A CN 110034827 A CN110034827 A CN 110034827A
- Authority
- CN
- China
- Prior art keywords
- error
- matrix
- module
- demultiplexing
- depolarization
- 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
-
- 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/6162—Compensation of polarization related effects, e.g., PMD, PDL
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Optical Communication System (AREA)
Abstract
The invention discloses a kind of depolarization multiplexing methods based on reverse observation error, belong to coherent light communication technical field.The method of the present invention initializes demultiplexing matrix A (n), the matrix that the demultiplexing matrix is one 2 × 2 first；Two-way polarization state signal X'(n), Y'(n) is inputted into demultiplexing matrix again and obtains demultiplexed signal；The constellation point X nearest apart from demultiplexed signal is found out later#(n)、Y#(n), demultiplexed signal and constellation point X are calculated#(n)、Y#(n) the preliminary error ε between1(n)；Again by constellation point X#(n)、Y#(n) the inverse matrix A of input demultiplexing matrix A (n)‑1(n) pseudo- observation signal X " (n), Y " (n) are obtained, constellation point X is calculated#(n)、Y#(n) the reverse observation error ε between puppet observation signal X " (n), Y " (n)2(n)；Gradient descent method is reused, preliminary error ε is substituted into1(n) and reverse observation error ε2(n) multiplex matrices A (n) is updated；Above step is finally repeated, solution obtains demultiplexed signal and updates demultiplexing matrix A (n).The invention also achieves a kind of depolarization multiplex systems based on reverse observation error.The present invention effectively raises the convergence rate of depolarization multiplexing algorithm.
Description
Technical field
The invention belongs to coherent light communication technical fields, inclined more particularly, to a kind of solution based on reverse observation error
Vibration multiplexing method and system.
Background technique
With the rapidly growth of communication service, the promotion of power system capacity is imperative.Polarization multiplexing
The polarization characteristic that (Polarization Multiplexing, PM) utilizes light to transmit in single mode optical fiber is equal by transmission wavelength
Two independent and mutually orthogonal polarization states as independent channel distinguish transmitting two paths of signals, significantly improved power system capacity,
Increase the availability of frequency spectrum.
Since to will receive polarization mode dispersion, Polarization Dependent Loss and input signal in transmission process non-for two polarization states
The influence of the factors such as orthogonal, so that polarization state changes to form mutual crosstalk.Under atrocious weather environment, thunder and lightning
Faraday effect caused by the high-intensity magnetic field that the heavy current of generation induces will cause the quick polarization state rotation of light, this is to reception
The demultiplexing convergence speed of the algorithm at end is a stern challenge.
In conventional method the most commonly used is palarization multiplexing algorithm be permanent mould (Constant Modulus Algorithm,
CMA) algorithm, constant modulus algorithm utilize the structure of butterfly, constantly structural parameters are adjusted by error function, to reach depolarization
The purpose of multiplexing, algorithm complexity is low, realizes simply, but its error function is only concerned the modulus value of complex signal, for phase noise
It is extremely sensitive, and it is only applicable to 4QAM modulation format.Although there is the improvement constant modulus algorithm based on this algorithm to solve phase noise
Problem, and the multi-modulus algorithm based on judgement are suitable for the modulation format of higher order, but its algorithm the convergence speed excessively slowly can not be
Stable demultiplexing is realized under lightning environment.Therefore we need a kind of faster algorithm of convergence rate to cause to realize to lightning
Quick polarization state rotation tracking.
Summary of the invention
Aiming at the above defects or improvement requirements of the prior art, the present invention provides a kind of solutions based on reverse observation error
Palarization multiplexing method and system, its object is to effectively raise solution by introducing reverse this new dimension of observation error
Palarization multiplexing convergence speed of the algorithm.Thus the technical issues of polarization state is tracked under lightning environment is solved.
To achieve the above object, described the present invention provides a kind of depolarization multiplexing method based on reverse observation error
Method specifically includes the following steps:
(1) initialization demultiplexing matrix A (n), the matrix that demultiplexing matrix A (n) is one 2 × 2；
(2) two-way polarization state signal X'(n), Y'(n) input A (n) is obtained into demultiplexed signal X#'(n)、Y#'(n)；
(3) distance X is found out#'(n)、Y#'(n) nearest constellation point X#(n)、Y#(n), X is calculated#'(n)、Y#'(n) and X#(n)、
Y#(n) the preliminary error ε between1(n)；
(4) by constellation point X#(n)、Y#(n) the inverse matrix A of input demultiplexing matrix A (n)-1(n) pseudo- observation signal X " is obtained
(n), Y " (n) calculates X#(n)、Y#(n) the reverse observation error ε between X " (n), Y " (n)2(n)；
(5) gradient descent method is used, error ε is substituted into1(n) and error ε2(n) multiplex matrices A (n) is updated；
(6) n=n+1 repeats step (2)~(5), exports demultiplexed signal X#'(n)、Y#'(n) and demultiplexing matrix is updated
A(n)。
Further, the error ε1(n) calculation formula is as follows:
Wherein,WithRespectively indicate X#'(n) real and imaginary parts；WithIt respectively indicates
X#(n) real and imaginary parts.
Further, the error ε2(n) calculation formula is as follows:
ε2(n)=(X "r(n)2-X′r(n)2)2+(X"i(n)2-X′i(n)2)2
Wherein, X "r(n) and X "i(n) " the real and imaginary parts of (n) that respectively indicate X；X′r(n) and X 'i(n) X' is respectively indicated
(n) real and imaginary parts.
Further, described (5) specifically:
Solve parameter
Solve parameter
It updates matrix A (n):
Wherein, α, β are weight factor, μ1、μ2For Learning Step.
It is another aspect of this invention to provide that the present invention provides a kind of, the depolarization based on reverse observation error is multiplexed system
System, the system comprises successively run with lower module:
Initialization module, for initializing demultiplexing matrix A (n), the square that demultiplexing matrix A (n) is one 2 × 2
Battle array；
Demultiplexing module, for two-way polarization state signal X'(n), Y'(n) input A (n) to be obtained demultiplexed signal X#'
(n)、Y#'(n)；
Preliminary error solves module, for finding out distance X#'(n)、Y#'(n) nearest constellation point X#(n)、Y#(n), it calculates
X#'(n)、Y#'(n) and X#(n)、Y#(n) the preliminary error ε between1(n)；
Reverse observation error solves module, is used for constellation point X#(n)、Y#(n) inverse matrix of input demultiplexing matrix A (n)
A-1(n) pseudo- observation signal X " (n), Y " (n) are obtained, X is calculated#(n)、Y#(n) the reverse observation error between X " (n), Y " (n)
ε2(n)；
Matrix update module substitutes into error ε for using gradient descent method1(n) and error ε2(n) multiplex matrices A is updated
(n)；
Iterative solution palarization multiplexing module repeats successively to execute demultiplexing module, preliminary error solution mould for setting n=n+1
Block, reverse observation error solve module, matrix update module, export demultiplexed signal X#'(n)、Y#'(n) and demultiplexing square is updated
Battle array A (n).
Further, the preliminary error solves error ε in module1(n) calculation formula is as follows:
Wherein,WithRespectively indicate X#'(n) real and imaginary parts；WithIt respectively indicates
X#(n) real and imaginary parts.
Further, the reverse observation error solves error ε in module2(n) calculation formula is as follows:
ε2(n)=(X "r(n)2-X′r(n)2)2+(X"i(n)2-X′i(n)2)2
Wherein, X "r(n) and X "i(n) " the real and imaginary parts of (n) that respectively indicate X；X′r(n) and X 'i(n) X' is respectively indicated
(n) real and imaginary parts.
Further, the iterative solution palarization multiplexing module is specifically used for:
Solve parameter
Solve parameter
It updates matrix A (n):
Wherein, α, β are weight factor, μ1、μ2For Learning Step.
In general, through the invention it is contemplated above technical scheme is compared with the prior art, have below beneficial to effect
Fruit:
(1) the application dimension that this is new by reverse observation error introduces depolarization multiplexing method, effectively raises solution
Thus palarization multiplexing convergence speed of the algorithm solves the technical issues of polarization state is tracked under lightning environment；
(2) to error ε in the application1(n) and ε2(n) method for using real and imaginary parts separate computations, can effectively correct
The damage of phase caused by channel；
(3) gradient descent method used herein updates multiplex matrices A (n), and this method belongs to single order optimisation technique, structure
Simply, facilitate calculating, spend the time less.
Detailed description of the invention
Fig. 1 is the signal processing flow schematic diagram of depolarization multiplexing algorithm in the embodiment of the present application；
Fig. 2 is the signal constellation (in digital modulation) schematic diagram of depolarization multiplexing algorithm in the embodiment of the present application.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that described herein, specific examples are only used to explain the present invention, not
For limiting the present invention.As long as in addition, technical characteristic involved in the various embodiments of the present invention described below that
Not constituting conflict between this can be combined with each other.
As shown in Figure 1, being that the signal processing flow of the depolarization multiplexing algorithm in embodiment based on reverse observation error shows
It is intended to；The signal wherein respectively carried for two mutually orthogonal polarization states of transmitting terminal X, Y, by the transmission of fiber channel, letter
Road available channel matrix T is simulated, the Jones matrix that matrix T is one 2 × 2, obtains drawing since polarization state rotates in receiving end
Signal X', Y' after two-way crosstalk are played, depolarization multiplex matrices A is input to, the matrix that matrix A is one 2 × 2 can use one
A 2 × 2 matrixTo describe；
Demultiplexed signal X is obtained by demultiplexing matrix#'、Y#'；
By being adjudicated apart from nearest principle, find from X#'、Y#'Nearest constellation point remembers that the constellation point is X#、Y#, by X#、Y#
Input demultiplexing inverse of a matrix matrix A-1, also referred to as analog channel matrix obtains pseudo- observation signal X ", Y ".
As shown in Fig. 2, being the signal constellation (in digital modulation) schematic diagram of the depolarization multiplexing algorithm in embodiment based on reverse observation error；
Wherein demultiplexed signal X#'With judgement constellation point X#Between error be defined as ε1, between pseudo- observation signal X " and observation signal X'
Error be defined as ε2；As seen from the figure, the calculation method of error we carried out using the mode that real and imaginary part is opened, it is specific to calculate
Formula is as follows:
ε2(n)=(X "r(n)2-X′r(n)2)2+(X"i(n)2-X′i(n)2)2
After obtaining two errors, we seek local derviation to error using gradient descent method, are finally demultiplexed using following formula to update
With matrix:
Wherein α and β is weight factor, is responsible for weight when adjustment updates between two errors；μ1、μ2For Learning Step, adjust
Whole convergence speed of the algorithm；Specific value needs to be adjusted according to the actual situation, is Learning Step, and Learning Step difference is right
Whole system also has an impact, and lesser Learning Step is unfavorable for the scene that high speed changes, and biggish Learning Step is then difficult to surely
It is fixed, it is easy to miss optimal value, specific different parameters needs are manually adjusted according to actual scene；Finally we are by introducing reverse see
Error is surveyed, demultiplexing algorithm is made to accelerate convergence rate, realizes the tracking of the polarization state rotation of higher speed.
The above content as it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention,
It is not intended to limit the invention, any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention,
It should all be included in the protection scope of the present invention.
Claims (8)
1. a kind of depolarization multiplexing method based on reverse observation error, which is characterized in that the described method comprises the following steps:
(1) initialization demultiplexing matrix A (n), the matrix that demultiplexing matrix A (n) is one 2 × 2；
(2) two-way polarization state signal X'(n), Y'(n) input A (n) is obtained into demultiplexed signal X#'(n)、Y#'(n)；
(3) distance X is found out#'(n)、Y#'(n) nearest constellation point X#(n)、Y#(n), X is calculated#'(n)、Y#'(n) and X#(n)、Y#
(n) the preliminary error ε between1(n)；
(4) by constellation point X#(n)、Y#(n) the inverse matrix A of input demultiplexing matrix A (n)-1(n) obtain pseudo- observation signal X " (n),
Y " (n) calculates X#(n)、Y#(n) the reverse observation error ε between X " (n), Y " (n)2(n)；
(5) gradient descent method is used, error ε is substituted into1(n) and error ε2(n) multiplex matrices A (n) is updated；
(6) n=n+1 repeats step (2)~(5), exports demultiplexed signal X#'(n)、Y#' (n) and update demultiplexing matrix A
(n)。
2. a kind of depolarization multiplexing method based on reverse observation error according to claim 1, which is characterized in that described
Error ε1(n) calculation formula is as follows:
Wherein,WithRespectively indicate X#'(n) real and imaginary parts；WithRespectively indicate X#(n)
Real and imaginary parts.
3. a kind of depolarization multiplexing method based on reverse observation error according to claim 1, which is characterized in that described
Error ε2(n) calculation formula is as follows:
ε2(n)=(X "r(n)2-X'r(n)2)2+(X”i(n)2-X'i(n)2)2
Wherein, X "r(n) and X "i(n) real and imaginary parts of X " (n) are respectively indicated；X'r(n) and X'i(n) X'(n is respectively indicated)
Real and imaginary parts.
4. a kind of depolarization multiplexing method based on reverse observation error according to claim 1, which is characterized in that described
(5) specifically:
Solve parameter
Solve parameter
It updates matrix A (n):
Wherein, α, β are weight factor, μ1、μ2For Learning Step.
5. a kind of depolarization multiplex system based on reverse observation error, which is characterized in that the system comprises what is successively run
With lower module:
Initialization module, for initializing demultiplexing matrix A (n), the matrix that demultiplexing matrix A (n) is one 2 × 2；
Demultiplexing module, for two-way polarization state signal X'(n), Y'(n) input A (n) to be obtained demultiplexed signal X#'(n)、Y#'
(n)；
Preliminary error solves module, for finding out distance X#'(n)、Y#'(n) nearest constellation point X#(n)、Y#(n), X is calculated#'
(n)、Y#'(n) and X#(n)、Y#(n) the preliminary error ε between1(n)；
Reverse observation error solves module, is used for constellation point X#(n)、Y#(n) the inverse matrix A of input demultiplexing matrix A (n)-1
(n) pseudo- observation signal X " (n), Y " (n) are obtained, X is calculated#(n)、Y#(n) the reverse observation error ε between X " (n), Y " (n)2
(n)；
Matrix update module substitutes into error ε for using gradient descent method1(n) and error ε2(n) multiplex matrices A (n) is updated；
Iterative solution palarization multiplexing module, for setting n=n+1, repetition successively executes demultiplexing module, preliminary error solves module,
Reverse observation error solves module, matrix update module, exports demultiplexed signal X#'(n)、Y#'(n) and demultiplexing matrix A is updated
(n)。
6. a kind of depolarization multiplex system based on reverse observation error according to claim 5, which is characterized in that described
Preliminary error solves error ε in module1(n) calculation formula is as follows:
Wherein,WithRespectively indicate X#'(n) real and imaginary parts；WithRespectively indicate X#(n)
Real and imaginary parts.
7. a kind of depolarization multiplex system based on reverse observation error according to claim 5, which is characterized in that described
Reverse observation error solves error ε in module2(n) calculation formula is as follows:
Wherein, X "r(n) and X "i(n) real and imaginary parts of X " (n) are respectively indicated；X'r(n) and X'i(n) X'(n is respectively indicated)
Real and imaginary parts.
8. a kind of depolarization multiplex system based on reverse observation error according to claim 5, which is characterized in that described
Iterative solution palarization multiplexing module is specifically used for:
Solve parameter
Solve parameter
It updates matrix A (n):
Wherein, α, β are weight factor, μ1、μ2For Learning Step.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910227443.7A CN110034827B (en) | 2019-03-25 | 2019-03-25 | Depolarization multiplexing method and system based on reverse observation error |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910227443.7A CN110034827B (en) | 2019-03-25 | 2019-03-25 | Depolarization multiplexing method and system based on reverse observation error |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110034827A true CN110034827A (en) | 2019-07-19 |
| CN110034827B CN110034827B (en) | 2020-11-17 |
Family
ID=67236552
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910227443.7A Active CN110034827B (en) | 2019-03-25 | 2019-03-25 | Depolarization multiplexing method and system based on reverse observation error |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110034827B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114650097A (en) * | 2020-12-21 | 2022-06-21 | 富士通株式会社 | Polarization state tracking recovery method, device and optical receiver |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160315637A1 (en) * | 2015-04-24 | 2016-10-27 | Zte Corporation | Turbo product coded modulation |
| CN106875002A (en) * | 2017-02-20 | 2017-06-20 | 中国石油大学(华东) | Complex value neural network training method based on gradient descent method Yu generalized inverse |
| CN107809282A (en) * | 2017-10-13 | 2018-03-16 | 北京邮电大学 | Equalization methods and equipment in the inclined division multiplexing system of coherent optical-fiber under extreme scenes |
| US10103799B2 (en) * | 2015-03-16 | 2018-10-16 | Lockheed Martin Corporation | Apparatus and method for increased data rates in underwater communications using orbital angular momentum |
| CN108933626A (en) * | 2017-05-24 | 2018-12-04 | 中兴通讯股份有限公司 | A kind of signal processing method and device |
-
2019
- 2019-03-25 CN CN201910227443.7A patent/CN110034827B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10103799B2 (en) * | 2015-03-16 | 2018-10-16 | Lockheed Martin Corporation | Apparatus and method for increased data rates in underwater communications using orbital angular momentum |
| US20160315637A1 (en) * | 2015-04-24 | 2016-10-27 | Zte Corporation | Turbo product coded modulation |
| CN106875002A (en) * | 2017-02-20 | 2017-06-20 | 中国石油大学(华东) | Complex value neural network training method based on gradient descent method Yu generalized inverse |
| CN108933626A (en) * | 2017-05-24 | 2018-12-04 | 中兴通讯股份有限公司 | A kind of signal processing method and device |
| CN107809282A (en) * | 2017-10-13 | 2018-03-16 | 北京邮电大学 | Equalization methods and equipment in the inclined division multiplexing system of coherent optical-fiber under extreme scenes |
Non-Patent Citations (2)
| Title |
|---|
| KAZURO KIKUCHI: "Performance analyses of polarization demultiplexing based on constant-modulus algorithm in digital coherent optical receivers", 《OPTICS EXPRESS》 * |
| 童程: "一种改进的光偏分复用16-QAM***盲解复用算法的研究", 《光电子·激光》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114650097A (en) * | 2020-12-21 | 2022-06-21 | 富士通株式会社 | Polarization state tracking recovery method, device and optical receiver |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110034827B (en) | 2020-11-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Eriksson et al. | 56 Gbaud probabilistically shaped PAM8 for data center interconnects | |
| CN107809282B (en) | Equalization methods and equipment in the inclined division multiplexing system of coherent optical-fiber under extreme scenes | |
| US9531480B2 (en) | Optical receiver | |
| US9843388B1 (en) | Laser communications | |
| CN103973614A (en) | Adaptive equalizer tap coefficient correction method and optical receiver | |
| US10148465B2 (en) | Training assisted joint equalization | |
| CN110034827A (en) | A kind of depolarization multiplexing method and system based on reverse observation error | |
| CN110011733A (en) | A kind of depolarization multiplexing method and system based on factor of momentum | |
| Bajaj et al. | 54.5 Tb/s WDM transmission over field deployed fiber enabled by neural network-based digital pre-distortion | |
| CN108933626A (en) | A kind of signal processing method and device | |
| Arık et al. | High-dimensional modulation for mode-division multiplexing | |
| CN114124289B (en) | Fast polarization state rotation compensation method based on cubature Kalman filter | |
| CN110224761A (en) | A kind of method and system of rapid solving polarization spin matrix and equalizer coefficients | |
| CN108521385A (en) | A kind of coherent light communication equalization methods using cross-mixing state | |
| Che et al. | Blind polarization de-multiplexing for stokes vector direct detection | |
| CN112468223A (en) | Operation simulation method and system for broadband communication constellation | |
| Omidi et al. | Optimizing geometric constellations for phase noise channels using deep learning | |
| Schädler et al. | Machine learning in digital signal processing for optical transmission systems | |
| Xian et al. | PDM PAM-4 with IM-DD using a simple MIMO DSP-based receiver for short reach communications | |
| CN115242308B (en) | Signal phase reconstruction method and system based on digital subcarrier multiplexing | |
| Buelow et al. | Experimental analysis of transmission and soft decoding of optimized 4D constellations | |
| Liu et al. | In-service blind transceiver calibration for long-haul non-dispersion managed systems | |
| Li et al. | A CGAN-aided Autoencoder Supporting Joint Geometric Probabilistic Shaping for Optical Fiber Communication System | |
| CN114866145B (en) | Method, apparatus, device and computer readable medium for optical communication | |
| Ip et al. | Temporal variations in the channel matrix in few-mode fiber recirculating loop transmission |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |