CN106525390A - Dispersion compensation method for optical fiber polarization maintaining device with ultra high distributed birefringence dispersion - Google Patents
Dispersion compensation method for optical fiber polarization maintaining device with ultra high distributed birefringence dispersion Download PDFInfo
- Publication number
- CN106525390A CN106525390A CN201610985383.1A CN201610985383A CN106525390A CN 106525390 A CN106525390 A CN 106525390A CN 201610985383 A CN201610985383 A CN 201610985383A CN 106525390 A CN106525390 A CN 106525390A
- Authority
- CN
- China
- Prior art keywords
- dispersion
- data
- compensation
- optical path
- head end
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/30—Testing of optical devices, constituted by fibre optics or optical waveguides
Landscapes
- Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Optical Communication System (AREA)
Abstract
The invention belongs to the technical field of optical coherence polarization measurement and particularly relates to a dispersion compensation method for an optical fiber polarization maintaining device with ultra high distributed birefringence dispersion, and the performance of a testing system is improved through inhibiting the influence of the ultra high distributed birefringence dispersion on the output signal of an optical coherence domain polarization measurement device. The method comprises a step of carrying out zero fill on the two ends of a distribution polarization crosstalk original data obtained by the distribution polarization crosstalk test of a polarization maintaining device to be tested by using the optical coherence polarization measurement device to obtain an initial data to be compensated, and a step of extracting the data segment of the data to be compensated in an optical path difference range as head end data. According to the method, the ultra high distributed birefringence dispersion compensation of the OCDP output signal can be realized without external hardware configuration, the flexibility is high, two-order and three-order birefringence dispersion can be compensated, and higher order birefringence dispersion can be compensated according to needs.
Description
Technical field
The invention belongs to optical coherence domain polarimetry technology (Optical Coherence Domain
Polarization, referred to as OCDP) field, and in particular to it is a kind of by suppressing the distributed birefringence dispersion of superelevation to optics phase
Dry domain polarimeter output signal impact come improve test system performance for the distributed birefringence color of superelevation
The dispersion compensation method of scattered optical fibre polarization-maintaining device.
Background technology
Optical coherence domain polarimetry technology (OCDP) is based on White Light Interferometer (White Light
Interferometry, referred to as WLI) a kind of technical scheme for measuring the distribution polarization interference for protecting inclined device.OCDP
Wide spectrum light source, such as super-radiance light emitting diode (SLD) are adopted typically, and in order to avoid interference peaks have intrinsic secondary lobe, typically
From the wide spectrum light source of Gaussian spectrum；Wide spectrum optical injects device under test through the polarizer, according to the different choice of device under test
The countershaft angle of the different polarizers, such as the detection of fiber optic loop, we generally use 0 ° of polarizer and wide spectrum optical are noted by slow axis
Enter, and for the detection of Y waveguide, generally use 45 ° of polarizers by wide spectrum optical by fast axle and slow axis while injection；Device under test
The light of fast and slow axis output is coupled to its single-mode output tail optical fiber by 45 ° of analyzers, is then entered by 2 × 2 bonders
One Mach increases Dare interferometer (Mach-Zehnder Interferometer, abbreviation MZI) or Michelson's interferometer
(Michelson Interferometer, abbreviation MI)；The signal of last interferometer output is entered by a difference detecting circuit
Row detection.The defect point protected in inclined device under test can cause the light in its fast and slow axis to couple to another axle, and treat unequally due to protecting
There is mode birefringence (Modal Birefringence) between the fast and slow axis for surveying device, thus it is defeated from device under test fast and slow axis
Can there is certain optical path difference between the light for going out, and the optical path difference is typically greater than the coherence length of laser.By interferometer
In displacement platform can carry out light path matching, the output light of device under test fast and slow axis is interfered by realization, and final output one
Serial interference peaks, these interference peaks can reflect that the parameter information such as the position of defect point, extinction ratio is (referring to first in device under test
Technology [1], Chinese patent CN103900680A).
Due to the birefringence dispersion (Birefringence Dispersion) existed between the fast and slow axis of device under test, meeting
Cause the deterioration of OCDP output signals, it is double that (referring to first technology [2], Chinese patent CN 102332956B) specifically includes second order
Birefringence dispersion can cause the reduction of the broadening and amplitude of interference peaks envelope, that is, reduce spatial resolution and make extinction ratio
Measured value distortion.Even in some cases, can there is division due to the presence of birefringence dispersion in two adjacent interference peaks
Situation, cause to obtain the test result (referring to first technology [3], Chinese patent CN 104006948B) of mistake.Therefore, it is double
Birefringence dispersion compensation becomes a key technology for realizing high-resolution, the OCDP systems of High Extinction Ratio certainty of measurement.
Birefringence dispersion compensation method is broadly divided into two class of hardware compensating and software compensation.Hardware compensating method typically exists
Dispersion compensation device is inserted in reference arm, matches the dispersion of interferometer two-arm, and during for birefringence linear change, can
Second order, third-order dispersion compensation are realized (referring to first technology with by a pair of parallel balzed grating, is embedded in displacement platform
[4], Chinese patent CN 100398057C).However, hardware compensating method needs, using complicated hardware unit, to improve displacement
The debugging difficulty of platform, and the birefringence dispersion amount of device under test known to needing.Software compensation method is typically without to existing system
The hardware of system is modified, but the signal to gathering carries out post processing to realize that birefringence dispersion is compensated, although can not be as hard
Part compensation method carries out real-Time Compensation like that, but software compensation method is more flexible, and typically without the need for known device under test
Birefringence dispersion amount.For the interference data exported using space optical path interferometer, can be by using rectangular window function point
The data of main peak data and a certain interference peaks for needing to carry out birefringence dispersion compensation are not intercepted out, then by calculating both
The ratio of envelope width carries out birefringence dispersion to carry out the estimation of birefringence dispersion coefficient, and then construction phase compensating factor
Compensation is (referring to first technology [2], Chinese patent CN 102332956B).But the technology is only effective to space optical path interferometer,
And need the interference peaks for carrying out birefringence dispersion compensation carry out envelope width measurement, for those are flooded by noise
Interference peaks cannot then realize birefringence dispersion compensate.
The content of the invention
It is an object of the invention to provide a kind of compensation is caused due to the presence of the distributed birefringence dispersion of superelevation
The deterioration of OCDP test system performances, including optical coherence domain polarization measurement device is improved for being distributed the position of polarization interference test
Put precision, spatial resolution, and the amplitude precision of extinction ratio test for the light with the distributed birefringence dispersion of superelevation
The fine dispersion compensation method for protecting inclined device.
The object of the present invention is achieved like this：
The present invention comprises the following steps：
(1) optical coherence domain polarization measurement device is protected inclined device and carries out being distributed that polarization interference test obtains point to be measured
The two ends of cloth polarization interference initial data carry out zero padding and obtain initial data I to be compensated1D (), wherein d >=0 are optical path differences, and
I1D () is in optical path difference scope d ∈ [0, Δ d1] and d ∈ [dmax-Δd1,dmax] in be added data 0, Δ d1It is taken as
Δd1=0.5mm, dmaxFor the maximum of optical path difference d；
(2) data I to be compensated are extracted1D () is in optical path difference scope d ∈ [Δ d1,d1+Δd1] in data segment as head end
Data Is1(d), d1It is taken as d1=4mm；
(3) to head end data Is1D () measures dispersion parameters, i.e. 2nd order chromatic dispersion G using optimization chromatic dispersion measurement technology1With three ranks
Dispersion T1, and determine head end data Is1D () is using dispersion parameters G1With T1Carry out the peak value of the result after FFT dispersion compensation modules
Optical path difference l of position1, then calculate distributed second order dispersion DGj=Gj/lj；
(4) using dispersion parameters G obtained in previous step1With T1, treat offset data I1D () is using FFT dispersion compensatings
Block carries out dispersion compensation, obtains local offset data Ilocal(d)；
(5) to local offset data IlocalD () carries out data segmentation, extract wherein optical path difference scope d ∈ [Δ d1,d1+Δ
d1] in data segment as the head end data I after compensationc1(d), and extract wherein optical path difference scope d ∈ [d1,dmax-Δd1] in
Data segment as endian data I after compensatione1(d), IlocalD ∈ [0, Δ d in (d)1] and d ∈ [dmax-Δd1,dmax]
Part is dropped, and head end data Ic1(d) and endian data Ie1Then there is length for optical path difference Δ d in (d)1Overlapping part；
(6) endian data I after inspection compensatione1Whether the data length of (d) meets δ d=dmax-d1≥d1It is if meeting, first
Range error Δ I caused by first remaining according to the dispersion for allowingdB=0.5dB calculates a parameter relevant with dispersionUpdate d1And Δ d1RespectivelyWith
WhereinFor wide spectrum light source coherence length, λ0, the centre wavelength and half Gao Quan of Δ λ respectively light sources
Width, c are the light velocity in vacuum, and SNR=90 is the dynamic range of system, by endian data I after compensatione1D the end of () is carried out
Zero padding is used as new data I to be compensated2D (), zero-filled data length are Δ d2, update optical path difference d to make its minima is 0；
(7) endian data I of the repeat step (2) to (6) after when compensationejD the data length of () is unsatisfactory for δ d=
dmax-d1≥d1Relation；
(8) by the head end data I after all compensationcjD () in turn carries out head and the tail splicing and is fully compensated for data Iout
(d)。
Described optimization chromatic dispersion measurement includes：
One object function S=Obj (I of constructionsj(d),Gj,Tj), which is with head end data Is1(d) and 2nd order chromatic dispersion G1With
Third-order dispersion T1As input, and FFT dispersion compensation modules are used to head end data I using dispersion parameterss1D () carries out dispersion benefit
Trial is repaid, and acutance S of compensation result is finally calculated as output；Object function S=Obj (I are searched for using optimization algorithmsj
(d),Gj,Tj) obtain maximum when dispersion parameters Gj,Tj。
Described FFT dispersion compensation modules include：
Treat offset data IjD () carries out fast Fourier transform, obtain its amplitude spectrum A (ω) and phase spectrumI.e.Wherein ω is light angular frequency；Dispersion compensation phase place is calculated using dispersion parameters simultaneouslyWherein ω0For light source center angular frequency；Then calculate the phase after dispersion compensation
Compose positionLocal offset data I is obtained finally by inverse fast Fourier transform IFFTlocal(d)
=IFFT [A (ω) exp [i φ (ω)]].
Compared with prior art, it is an advantage of the current invention that：
1st, the present invention is capable of achieving to carry out OCDP output signals the distributed birefringence of superelevation without the need for extra hardware configuration
Dispersion compensation, and flexibility ratio is high, can compensate second order, three rank birefringence dispersions, and can compensate higher order according to demand
Birefringence dispersion.
2nd, birefringence dispersion this priori not only without the need for device under test, and for those as birefringence is accumulated
The excessive and caused interference peaks flooded by noise of amount can also realize that birefringence dispersion is compensated.
3rd, compared to data to be directly divided into multiple segment individual processing, the present invention has faster processing speed.
Description of the drawings
Fig. 1 is OCDP test principle figures；
Fig. 2 is the theory diagram of the distributed birefringence dispersion compensation of superelevation；
Fig. 3 is the emulation signal of OCDP test systems and its output result after birefringence dispersion compensation.
Fig. 4 is the result that there is to some polarization-maintaining fiber coil measured datas for demarcating solder joint to carry out dispersion compensation.
Specific embodiment
Below in conjunction with the accompanying drawings the present invention is described in more detail：
A kind of dispersion compensation method for the optical fibre polarization-maintaining device with the distributed birefringence dispersion of superelevation.The method is led to
The dispersion parameters for intercepting that head end data carries out obtaining herein based on optimized chromatic dispersion measurement technology are crossed, and is compensated with this and is needed
Offset data, after intercepting and preserving the head end data of perfect compensation, remaining endian data is carried out above-mentioned steps again,
Until all data are by perfect compensation, the distributed birefringence dispersion of the superelevation that is eliminated after finally splice all of which affects
Polarization interference data；Amount of redundancy of the method when operation and setting data segmentation of data padding is come after ensureing to compensate
Data do not receive termination effects such that it is able to directly carry out head and the tail splicing.The method using chromatic dispersion measurement technology make dispersion
Parameter determination does not need high s/n ratio interference peaks, and the suitability is higher；And gradually dispersion compensation the characteristics of those can be made to be made an uproar
The interference peaks that sound floods gradually appear with iterative process, finally by perfect compensation.
The present invention includes：
1st, optical coherence domain polarization measurement device is carried out being distributed what polarization interference test was obtained to the inclined device of guarantor to be measured
The two ends of distribution polarization interference initial data carry out zero padding and obtain initial data I to be compensated1D (), wherein d >=0 are optical path differences,
And I1D () is in optical path difference scope d ∈ [0, Δ d1] and d ∈ [dmax-Δd1,dmax] in be added data 0, Δ here
d1It is taken as Δ d1=0.5mm, dmaxFor the maximum of optical path difference d；
2nd, data I to be compensated are extracted1D () is in optical path difference scope d ∈ [Δ d1,d1+Δd1] in data segment as head end
Data Is1(d), here d1It is taken as d1=4mm；
3rd, to head end data Is1D () measures dispersion parameters, i.e. 2nd order chromatic dispersion G using optimization chromatic dispersion measurement technology1With three
Rank dispersion T1, and determine head end data Is1D () is using dispersion parameters G1With T1Carry out the peak of the result after FFT dispersion compensation modules
Optical path difference l of value position1, then calculate distributed second order dispersion DGj=Gj/lj；
4th, using dispersion parameters G obtained in previous step1With T1, treat offset data I1D () is using FFT dispersion compensations
Module carries out dispersion compensation, obtains local offset data Ilocal(d)；
5th, to local offset data IlocalD () carries out data segmentation, extract wherein optical path difference scope d ∈ [Δ d1,d1+
Δd1] in data segment as the head end data I after compensationc1(d), and extract wherein optical path difference scope d ∈ [d1,dmax-Δd1]
Interior data segment is used as endian data I after compensatione1(d), here IlocalD ∈ [0, Δ d in (d)1] and d ∈ [dmax-Δd1,
dmax] part be dropped, and head end data Ic1(d) and endian data Ie1Then there is length for optical path difference Δ d in (d)1Overlap
Point；
6th, endian data I after inspection compensatione1Whether full -0 foot δ d=d of the data length of (d)max-d1≥d1If, it is full
Foot, then it is assumed that data length is sufficiently large, the then dispersion for being first depending on allowing remain caused range error Δ IdB=0.5dB is calculated
One parameter relevant with dispersionThen update d1And Δ d1Respectively
WithWhereinFor wide spectrum light source coherence length, λ0, Δ λ point
Not Wei light source centre wavelength and full width at half maximum, c be vacuum in the light velocity, and SNR=90 for system dynamic range, then will
Endian data I after compensatione1D the end of () carries out zero padding as new data I to be compensated2D () (zero-filled data length is Δ
d2), and update optical path difference d to make its minima to be 0；
7th, repeat step 2 to 6 (the lower footnote of wherein all variables adds 1), endian data I after when compensationej(d)
Data length be unsatisfactory for δ d=dmax-d1≥d1Relation；
8th, by the head end data I after all compensationcjD () in turn carries out head and the tail splicing and is fully compensated for data
Iout(d)。
Described optimization chromatic dispersion measurement is：
An object function S=Obj (I is constructed firstsj(d),Gj,Tj), which is with head end data Is1(d) and 2nd order chromatic dispersion
G1With third-order dispersion T1As input, and FFT dispersion compensation modules are used to head end data I using dispersion parameterss1D () carries out color
Dissipate compensation to attempt, acutance S of compensation result is finally calculated as output；Secondly, object function S=is searched for using optimization algorithm
Obj(Isj(d),Gj,Tj) obtain maximum when dispersion parameters Gj,Tj。
Described FFT dispersion compensation modules are：
Treat offset data IjD () carries out fast Fourier transform (FFT), obtain its amplitude spectrum A (ω) and phase spectrumI.e.Wherein ω is light angular frequency；Dispersion is calculated using dispersion parameters simultaneously to mend
Repay phase placeWherein ω0For light source center angular frequency；Then calculate dispersion compensation
Phase spectrum afterwards isLocal offset data is obtained finally by inverse fast Fourier transform IFFT
Ilocal(d)=IFFT [A (ω) exp [i φ (ω)]].
It is OCDP test principle figures shown in Fig. 1, it is 1550nm that super continuous luminous diode 101 sends centre wavelength,
Wide spectrum optical of half spectral width for the Gaussian spectrum shape of 50nm, through one 2:98 1 × 2 bonder 102, wherein 2% light is entered
Photodiode 103 is used for real-time monitoring light source works state, and remaining 98% light enters isolator 104 to be used to intercept back instead
Light, then via the polarizer 105 and the injection device under test 107 of polarizer tail optical fiber 106, determines according to the type of device under test 107
Determine the countershaft angle of the polarizer 105, if device under test 105 is polarization-maintaining fiber coil, the countershaft angle of the polarizer 105 is 0 °, if
Device under test 105 is Y waveguide, then the countershaft angle of the polarizer 105 is 45 °, and the light for injecting device under test 105 can be due to wherein depositing
Defect point and there is the coupling between fast and slow axis, the light exported from the fast and slow axis of device under test 105 is via analyzer tail optical fiber
108 and analyzer 109,2 × 2 bonders 112 are then injected into, the operation wavelength of 2 × 2 bonder 112 is 1310nm/
1550nm, and output ratio is 50:The light of 50,2 × 2 bonder, 112 two-arm output respectively enters reference arm 113 and scanning
Arm one end 114, at the same time, the centre wavelength of the output of distributed feedback laser 110 is the laser of 1310nm, through isolator
2 × 2 bonders 112 are also passed through after 111 and enters the reference arm 113 and scan arm 114 that Mach increases Dare interferometer, wherein entering
The light of scan arm one end 114 can initially enter a circulator 115, and through 116 freedom of entry space of GRIN Lens,
Return into GRIN Lens 116 after displacement platform reflecting mirror 117 reflects, again pass by circulator 115 another into scan arm
End 118, the laser for turning finally to the wide spectrum optical and 1310nm wavelength of 1550nm centre wavelengths in arm 113 and scan arm 118 are equal
Interfere at 2 × 2 bonders 119, interference signal is input into wavelength division multiplexer 120 and 121 respectively by 2 × 2 bonders, wherein
Wide spectrum optical is detected via photodiode 122 and photodiode 123, then carries out calculus of differences by difference engine 125, and
The laser of 1310nm wavelength is then detected by the photodiode 124 of 121 one end of wavelength division multiplexer, last all of detection letter
Number gather via data acquisition unit 126, the data acquisition unit 126 includes a data collecting card and a computer.
Fig. 2 is the theory diagram of the distributed birefringence dispersion compensation of superelevation, and data acquisition unit 126 is by sync bit moving stage
The motion of reflecting mirror 117 and the signal acquisition process of its own are distributed polarization interference output signal as initial data to obtain
The head and the tail two ends of initial data 201 are carried out data padding 202 and obtain data to be compensated 203 by 201 (wide range optical interference signals),
Head end data 204 is subsequently extracted to which, and the second order of head end data 204, third-order dispersion is obtained using optimization chromatic dispersion measurement 205
206, the position of peak value in the result of dispersion compensation is further carried out according to second order, third-order dispersion 206 to head end data 204, is counted
Distributed second order dispersion 207 is calculated, then using second order, third-order dispersion parameter 206, offset data 203 is treated and is used FFT dispersions
Compensating module 208 carries out dispersion compensation, obtains local offset data 209, and next step carries out data segmentation 210, obtains head to which
End data 211 and endian data 212, during data are split, have been thrown aside local offset data head and the tail two ends and have been mended with data
Institute's zero padding number identical data volume in 0 202 operations, and the head end data 211 for obtaining there is also this with endian data 212
The overlapping part of individual data volume, then carries out the whether sufficiently large judgement 213 of data length, if data to endian data 212
Length is more than 1 times of head end data, that is, think data long enough, then endian data 212 is carried out data padding 202 again
Aforementioned 202 to 213 operation is carried out, until judging that 213 obtain result for "No", finally splices all head end datas successively from beginning to end
214, that is, obtain required full remuneration data 215.
The specific implementation step of this method is described in detail with reference to Fig. 3.It is by computer journey shown in the first width of Fig. 3
The OCDP test system output signals of sequence emulation, it can be seen that wherein with 4 interference peaks 301,302,303, and 304, he
Receive a certain degree of dispersive influence respectively, and if in fact, there is no the presence of birefringence dispersion, the width of this 4 interference peaks
Degree should be consistent；Shown in the second width of Fig. 3 is to carry out the result after an iteration using this method, specially by Fig. 3
Data head and the tail zero padding shown in one width, then intercepts part that head end optical path difference is 0.5～4.5mm (namely comprising interference peaks 301
Part) chromatic dispersion measurement is optimized, it is used for compensating all data in the first width of Fig. 3 after obtaining dispersion parameters, compensation result goes
Fall head and the tail it is a bit of after as shown in the second width of Fig. 3, it can be seen that 4 after compensation corresponding with 4 interference peaks in the first width of Fig. 3
The amplitude and width of individual interference peaks makes moderate progress, and especially the birefringence dispersion of interference peaks 305 is compensated completely,
So far first time iterative compensation is completed；The part that optical path difference in the second width of Fig. 3 is 3.5～20mm is taken out, second is carried out repeatedly
For dispersion compensation, as a result as shown in the 3rd width of Fig. 3；Third time iteration dispersion compensation is carried out to the data shown in the 3rd width of Fig. 3
As a result as shown in the 4th width of Fig. 3；Data shown in the 4th width of Fig. 3 are carried out with result such as Fig. 3 of the 4th iteration dispersion compensation
It is shown in five width, convenient herein for contrast, when the result of each wheel birefringence dispersion compensation is shown, will be processed good
Signal splices before current signal to be compensated.The method through us be can be seen that from the result shown in the 5th width of last Fig. 3
After carrying out birefringence dispersion compensation, the signal to noise ratio of 4 interference peaks is realized and is obviously improved, under effect of noise, 4 interference peaks
Amplitude it is basically identical.
Finally, Fig. 4 is illustrated and is processed polarization-maintaining fiber coil measured datas with some demarcation solder joints using this method
Carry out the result of dispersion compensation, the first width initial data, after the second width dispersion compensation.
Claims (3)
1. a kind of dispersion compensation method for the optical fibre polarization-maintaining device with the distributed birefringence dispersion of superelevation, its feature exist
In comprising the following steps：
(1) optical coherence domain polarization measurement device is protected inclined device to carry out being distributed the distribution that polarization interference test obtains inclined to be measured
The shake two ends of crosstalk initial data carry out zero padding and obtain initial data I to be compensated1D (), wherein d >=0 are optical path differences, and I1(d)
In optical path difference scope d ∈ [0, Δ d1] and d ∈ [dmax-Δd1,dmax] in be added data 0, Δ d1It is taken as Δ d1=
0.5mm, dmaxFor the maximum of optical path difference d；
(2) data I to be compensated are extracted1D () is in optical path difference scope d ∈ [Δ d1,d1+Δd1] in data segment as head end data
Is1(d), d1It is taken as d1=4mm；
(3) to head end data Is1D () measures dispersion parameters, i.e. 2nd order chromatic dispersion G using optimization chromatic dispersion measurement technology1And third-order dispersion
T1, and determine head end data Is1D () is using dispersion parameters G1With T1Carry out the peak of the result after FFT dispersion compensation modules
Optical path difference l1, then calculate distributed second order dispersion DGj=Gj/lj；
(4) using dispersion parameters G obtained in previous step1With T1, treat offset data I1D () is entered using FFT dispersion compensation modules
Row dispersion compensation, obtains local offset data Ilocal(d)；
(5) to local offset data IlocalD () carries out data segmentation, extract wherein optical path difference scope d ∈ [Δ d1,d1+Δd1] in
Data segment as the head end data I after compensationc1(d), and extract wherein optical path difference scope d ∈ [d1,dmax-Δd1] in number
According to section as endian data I after compensatione1(d), IlocalD ∈ [0, Δ d in (d)1] and d ∈ [dmax-Δd1,dmax] part
It is dropped, and head end data Ic1(d) and endian data Ie1Then there is length for optical path difference Δ d in (d)1Overlapping part；
(6) endian data I after inspection compensatione1Whether the data length of (d) meets δ d=dmax-d1≥d1If meeting, being first depending on
Range error Δ I caused by the dispersion residual of permissiondB=0.5dB calculates a parameter relevant with dispersion
Update d1And Δ d1RespectivelyWithWhereinFor wide spectrum light source coherence length, λ0, the centre wavelength and full width at half maximum of Δ λ respectively light sources, c
For the light velocity in vacuum, dynamic ranges of the SNR=90 for system, by endian data I after compensatione1D the end of () carries out zero padding
As new data I to be compensated2D (), zero-filled data length are Δ d2, update optical path difference d to make its minima is 0；
(7) endian data I of the repeat step (2) to (6) after when compensationejD the data length of () is unsatisfactory for δ d=dmax-d1
≥d1Relation；
(8) by the head end data I after all compensationcjD () in turn carries out head and the tail splicing and is fully compensated for data Iout(d)。
2. a kind of dispersion for the optical fibre polarization-maintaining device with the distributed birefringence dispersion of superelevation according to claim 1
Compensation method, it is characterised in that described optimization chromatic dispersion measurement includes：
One object function S=Obj (I of constructionsj(d),Gj,Tj), which is with head end data Is1(d) and 2nd order chromatic dispersion G1With three ranks
Dispersion T1As input, and FFT dispersion compensation modules are used to head end data I using dispersion parameterss1D () carries out dispersion compensation and tastes
Examination, finally calculates acutance S of compensation result as output；Object function S=Obj (I are searched for using optimization algorithmsj(d),Gj,
Tj) obtain maximum when dispersion parameters Gj,Tj。
3. a kind of dispersion for the optical fibre polarization-maintaining device with the distributed birefringence dispersion of superelevation according to claim 1
Compensation method, it is characterised in that described FFT dispersion compensation modules include：
Treat offset data IjD () carries out fast Fourier transform, obtain its amplitude spectrum A (ω) and phase spectrumI.e.Wherein ω is light angular frequency；Dispersion compensation phase place is calculated using dispersion parameters simultaneouslyWherein ω0For light source center angular frequency；Then calculate the phase after dispersion compensation
Compose positionLocal offset data I is obtained finally by inverse fast Fourier transform IFFTlocal(d)
=IFFT [A (ω) exp [i φ (ω)]].
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610985383.1A CN106525390B (en) | 2016-11-09 | 2016-11-09 | A kind of dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610985383.1A CN106525390B (en) | 2016-11-09 | 2016-11-09 | A kind of dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106525390A true CN106525390A (en) | 2017-03-22 |
| CN106525390B CN106525390B (en) | 2018-10-26 |
Family
ID=58350363
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610985383.1A Active CN106525390B (en) | 2016-11-09 | 2016-11-09 | A kind of dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106525390B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109697270A (en) * | 2019-01-16 | 2019-04-30 | 中国工程物理研究院激光聚变研究中心 | A kind of light beam dispersion characteristics inversion algorithm based on spatial spectral interference |
| CN113804406A (en) * | 2021-08-16 | 2021-12-17 | 广东工业大学 | Birefringence dispersion iterative compensation method for optical frequency domain polarization crosstalk measurement |
| CN113804405A (en) * | 2021-08-16 | 2021-12-17 | 广东工业大学 | Micro optical fiber dispersion measuring device based on double-coupler ring optical path structure |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1249813A (en) * | 1997-11-28 | 2000-04-05 | 富士通株式会社 | Polarization mode dispersion measuring method, and dispersion compensation control device and dispersion compensation control method |
| JP4008470B2 (en) * | 2002-07-19 | 2007-11-14 | 株式会社フジクラ | Measuring method and apparatus for measuring polarization mode dispersion of optical fiber |
| CN101819086A (en) * | 2010-05-19 | 2010-09-01 | 中国人民解放军国防科学技术大学 | Optical fiber dispersion measurement system and use method thereof |
| JP2011007806A (en) * | 2004-12-14 | 2011-01-13 | Luna Innovations Inc | Compensation method for time-varying phase in interferometric measurement |
| CN102332956A (en) * | 2011-08-23 | 2012-01-25 | 天津大学 | Dispersion compensation method for broadband light source |
| CN102914421A (en) * | 2012-10-19 | 2013-02-06 | 苏州光环科技有限公司 | Method and device for measuring polarization crosstalk in optical double-refraction medium |
-
2016
- 2016-11-09 CN CN201610985383.1A patent/CN106525390B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1249813A (en) * | 1997-11-28 | 2000-04-05 | 富士通株式会社 | Polarization mode dispersion measuring method, and dispersion compensation control device and dispersion compensation control method |
| JP4008470B2 (en) * | 2002-07-19 | 2007-11-14 | 株式会社フジクラ | Measuring method and apparatus for measuring polarization mode dispersion of optical fiber |
| JP2011007806A (en) * | 2004-12-14 | 2011-01-13 | Luna Innovations Inc | Compensation method for time-varying phase in interferometric measurement |
| CN101819086A (en) * | 2010-05-19 | 2010-09-01 | 中国人民解放军国防科学技术大学 | Optical fiber dispersion measurement system and use method thereof |
| CN102332956A (en) * | 2011-08-23 | 2012-01-25 | 天津大学 | Dispersion compensation method for broadband light source |
| CN102914421A (en) * | 2012-10-19 | 2013-02-06 | 苏州光环科技有限公司 | Method and device for measuring polarization crosstalk in optical double-refraction medium |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109697270A (en) * | 2019-01-16 | 2019-04-30 | 中国工程物理研究院激光聚变研究中心 | A kind of light beam dispersion characteristics inversion algorithm based on spatial spectral interference |
| CN109697270B (en) * | 2019-01-16 | 2022-04-01 | 中国工程物理研究院激光聚变研究中心 | Light beam dispersion characteristic inversion algorithm based on spatial spectrum interference |
| CN113804406A (en) * | 2021-08-16 | 2021-12-17 | 广东工业大学 | Birefringence dispersion iterative compensation method for optical frequency domain polarization crosstalk measurement |
| CN113804405A (en) * | 2021-08-16 | 2021-12-17 | 广东工业大学 | Micro optical fiber dispersion measuring device based on double-coupler ring optical path structure |
| CN113804406B (en) * | 2021-08-16 | 2023-12-26 | 广东工业大学 | Birefringent dispersion iterative compensation method for optical frequency domain polarization crosstalk measurement |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106525390B (en) | 2018-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1705471B1 (en) | Apparatus for measuring differential mode delay of multimode optical fiber | |
| US5712704A (en) | Appliance for measuring polarization mode dispersion and corresponding measuring process | |
| WO2019029163A1 (en) | Apparatus and method for eliminating polarization fading in ofdr | |
| US6204924B1 (en) | Method and apparatus for measuring polarization mode dispersion of optical devices | |
| US9041935B2 (en) | Measuring polarization crosstalk in optical birefringent materials and devices based on reduction of line broadening caused by birefringent dispersion | |
| ITPD970025A1 (en) | REFLECTOMETRIC INSTRUMENT FOR THE MEASUREMENT OF BIRIFRANGENCE DISTRIBUTED IN SINGLE-MODE FIBER OPTICS | |
| CN102928199A (en) | Device and method for improving polarization crosstalk measurement performance of optical device | |
| CN106092078B (en) | A kind of birefringence dispersion real-time compensation polarization beam splitting interferometer | |
| CA2709779A1 (en) | System and method to determine chromatic dispersion in short lengths of waveguides using a 3-wave interference pattern and a single-arm interferometer | |
| Calvani et al. | Polarization measurements on single-mode fibers | |
| CN106525390A (en) | Dispersion compensation method for optical fiber polarization maintaining device with ultra high distributed birefringence dispersion | |
| CN104792503A (en) | Device for enhancing sensitivity of distributed crosstalk measurement of optical polarizing device | |
| JP4241252B2 (en) | Optical fiber characteristic measuring apparatus and optical fiber characteristic measuring method | |
| CN103900680A (en) | Device and detecting method for restraining polarization crosstalk measuring noise by the adoption of light source | |
| CN106768877A (en) | A kind of Larger Dynamic scope scaling method for optical coherence domain polarimeter | |
| CN112082735B (en) | Optical fiber sensing ring bidirectional synchronous measurement device and method based on Sagnac structure | |
| CN107271147B (en) | A kind of optical fibre gyro polarization-maintaining fiber coil distribution polarization interference data processing method | |
| CN105066899B (en) | A kind of optical-fiber laser static strain sensing demodulating system of mutual reference | |
| CN104280217A (en) | Dual-channel optical performance measuring device for Y waveguide | |
| US6671038B2 (en) | Measurement of polarization dependent characteristic of optical components | |
| US10363101B2 (en) | Optical shape sensing system, medical apparatus and method for optical shape sensing | |
| CN208805343U (en) | The spectral measurement device of optical cable jumper head internal optical fiber breakpoint | |
| CN107764517B (en) | Method for eliminating second-order pseudo coupling point of interference signal of white light interferometer | |
| Fang et al. | Tuning Rate Estimation for Improving Positioning Accuracy in OFDR System | |
| US6600148B2 (en) | Polarization mode dispersion measuring method and polarization mode dispersion measuring system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |