CN106525390B - A kind of dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion - Google Patents
A kind of dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion Download PDFInfo
- Publication number
- CN106525390B CN106525390B CN201610985383.1A CN201610985383A CN106525390B CN 106525390 B CN106525390 B CN 106525390B CN 201610985383 A CN201610985383 A CN 201610985383A CN 106525390 B CN106525390 B CN 106525390B
- Authority
- CN
- China
- Prior art keywords
- dispersion
- data
- compensation
- polarization
- distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- 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
Abstract
The invention belongs to optical coherence domain polarimetry technical fields, and in particular to a kind of to improve the dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion of test system performance by inhibiting superelevation distribution birefringence dispersion to the influence of optical coherence domain polarization measurement device output signal.The present invention includes:The both ends progress zero padding that optical coherence domain polarization measurement device carries out polarization-maintaining device to be measured to the distribution polarization interference initial data that distribution polarization interference is tested obtains initial data to be compensated;Data segment of the data to be compensated within the scope of optical path difference is extracted as head end data etc..The present invention can be realized without additional hardware configuration and carry out the compensation of superelevation distribution birefringence dispersion to OCDP output signals, and flexibility ratio is high, can compensate second order, three rank birefringence dispersions, and the birefringence dispersion of higher order can be compensated according to demand.
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 inhibiting superelevation distribution birefringence dispersion to optics phase
The influence of dry domain polarimeter output signal come improve test system performance for having superelevation distribution birefringence color
The dispersion compensation method of scattered optical fibre polarization-maintaining device.
Background technology
Optical coherence domain polarimetry technology (OCDP) is to be based on White Light Interferometer (White Light
Interferometry, referred to as WLI) for measure polarization-maintaining device distribution polarization interference a kind of technical solution.OCDP
Wide spectrum light source is generally used, such as super-radiance light emitting diode (SLD), and in order to avoid interference peaks have intrinsic secondary lobe, generally
Select the wide spectrum light source of Gaussian spectrum;Wide spectrum optical injects device under test by the polarizer, the different selections according to device under test
The different polarizers is to shaft angle degree, and such as the detection of fiber optic loop, we are generally noted wide spectrum optical by slow axis using 0 ° of polarizer
Enter, and for the detection of Y waveguide, generally wide spectrum optical is injected simultaneously by fast axle and slow axis using 45 ° of polarizers;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 couplers
One Mach increases Dare interferometer (Mach-Zehnder Interferometer, abbreviation MZI) or Michelson's interferometer
(Michelson Interferometer, abbreviation MI);The signal of last interferometer output by a difference detecting circuit into
Row detection.The defects of polarization-maintaining device under test point can cause the light in its fast and slow axis to be coupled to another axis, and since polarization-maintaining waits for
It surveys between the fast and slow axis of device there are mode birefringence (Modal Birefringence), therefore defeated from device under test fast and slow axis
There can be certain optical path difference between the light gone out, and the optical path difference is typically greater than the coherence length of laser.Pass through interferometer
In displacement platform can carry out light path matching, the output light of device under test fast and slow axis is interfered in realization, and final output one
Serial interference peaks, these interference peaks can reflect that the parameter informations such as the position of defect point, extinction ratio are (referring to first in device under test
Ji Shu [1], Chinese patent CN103900680A).
Due to existing birefringence dispersion (Birefringence Dispersion), meeting between the fast and slow axis of device under test
Lead to the deterioration of OCDP output signals, (referring to first Ji Shu [2], Chinese patent CN 102332956B) and specifically include second order pair
Birefringence dispersion can lead to the broadening of interference peaks envelope and the reduction of amplitude, that is, reduce spatial resolution and make extinction ratio
Measured value is distorted.Even in some cases, two adjacent interference peaks can be divided due to the presence of birefringence dispersion
The case where, the test result for causing to obtain mistake is (referring to first Ji Shu [3], Chinese patent CN 104006948B).Therefore, double
Birefringence dispersion compensation, which becomes, realizes high-resolution, a key technology of the OCDP systems of High Extinction Ratio measurement accuracy.
Birefringence dispersion compensation method is broadly divided into two class of hardware compensating and software compensation.Hardware compensating method typically exists
It is inserted into dispersion compensation device in reference arm, the dispersion of interferometer two-arm is made to match, and the case where for birefringence linear change, it can
To realize second order, third-order dispersion compensation (referring to first technology by a pair of parallel balzed grating, of the insertion in displacement platform
[4], Chinese patent CN 100398057C).However, hardware compensating method needs, using complicated hardware device, to improve displacement
The debugging difficulty of platform, and need the birefringence dispersion amount of known device under test.Software compensation method is typically without to existing system
The hardware of system is modified, but is post-processed to the signal of acquisition to realize that birefringence dispersion compensates, although cannot be as hard
Part compensation method carries out real-time compensation like that, but software compensation method wants more flexible, and is generally not necessarily to known device under test
Birefringence dispersion amount.It, can be by using rectangular window function point for the interference data of use space light path interferometer output
Main peak data are not intercepted out and need to carry out the data of a certain interference peaks of birefringence dispersion compensation, then by calculating the two
The ratio of envelope width carries out the estimation of birefringence dispersion coefficient, and then constructs phase compensating factor to carry out birefringence dispersion
Compensation is (referring to first Ji Shu [2], Chinese patent CN 102332956B).However the technology is only effective to space optical path interferometer,
And it needs the interference peaks for carrying out birefringence dispersion compensation that can carry out envelope width measurement, those is flooded by noise
Interference peaks then cannot achieve birefringence dispersion compensation.
Invention content
The purpose of the present invention is to provide a kind of compensation caused by the presence of superelevation distribution birefringence dispersion
OCDP tests the deterioration of system performance, including improves the position that optical coherence domain polarization measurement device is used to be distributed polarization interference test
Set the light for having superelevation distribution birefringence dispersion of the amplitude precision of precision, spatial resolution and extinction ratio test
The dispersion compensation method of fine polarization-maintaining device.
The object of the present invention is achieved like this:
The present invention includes the following steps:
(1) optical coherence domain polarization measurement device is subjected to point that distribution polarization interference is tested to polarization-maintaining device to be measured
The both ends of cloth polarization interference initial data carry out zero padding and obtain initial data I to be compensated1(d), wherein d >=0 is optical path difference, and
I1(d) in optical path difference range d ∈ [0,Δd1]And d ∈ [dmax-Δd1,dmax]It is inside added data 0, Δ d1It is taken as
Δd1=0.5mm, dmaxFor the maximum value of optical path difference d;
(2) data I to be compensated is extracted1(d) in optical path difference range d ∈ [Δd1,d1+Δd1]Interior data segment is as head end
Data Is1(d), d1It is taken as d1=4mm;
(3) to head end data Is1(d) dispersion parameters, i.e. second-order dispersion G are measured using optimization chromatic dispersion measurement technology1With three ranks
Dispersion T1, and determine head end data Is1(d) dispersion parameters G is used1With T1Carry out the peak value of the result after FFT dispersion compensation modules
The optical path difference l of position1, then calculate distributed second order dispersion DGj=Gj/lj;
(4) the dispersion parameters G obtained in previous step is utilized1With T1, treat offset data I1(d) FFT dispersion compensatings are used
Block carries out dispersion compensation, obtains local offset data Ilocal(d);
(5) to local offset data Ilocal(d) data segmentation is carried out, wherein optical path difference range d ∈ [ are extracted;Δd1,d1+Δ
d1]Interior data segment is as the head end data I after compensationc1(d), and wherein optical path difference range d ∈ [ are extracted;d1,dmax-Δd1]It is interior
Data segment as compensation after endian data Ie1(d), Ilocal(d) d ∈ [ in;0,Δd1]And d ∈ [dmax-Δd1,dmax]'s
Part is dropped, and head end data Ic1(d) with endian data Ie1(d) then there is length is optical path difference Δ d1Overlapping part;
(6) the endian data I after compensation is examinede1(d) whether data length meets δ d=dmax-d1≥d1, if satisfied,
Then it is first depending on range error Δ I caused by the dispersion residual of permissiondB=0.5dB calculates a parameter related with dispersionUpdate d1And Δ d1RespectivelyWith
WhereinFor wide spectrum light source coherence length, λ0, Δ λ is respectively the centre wavelength and half Gao Quan of light source
Width, c are the light velocity in vacuum, and SNR=90 is the dynamic range of system, by the endian data I after compensatione1(d) end carries out
Zero padding is as new data I to be compensated2(d), zero-filled data length is Δ d2, updating optical path difference d makes its minimum value be 0;
(7) endian data I of the step (2) to (6) after when compensation is repeatedej(d) data length is unsatisfactory for δ d=
dmax-d1≥d1Relationship;
(8) by the head end data I after all compensationcj(d) it in turn carries out head and the tail splicing and is fully compensated for data Iout
(d)。
The optimization chromatic dispersion measurement includes:
Construct an object function S=Obj (Isj(d),Gj,Tj), with head end data Is1(d) and second-order dispersion G1With
Third-order dispersion T1Use FFT dispersion compensation modules to head end data I as input, and using dispersion parameterss1(d) dispersion benefit is carried out
Trial is repaid, finally calculates the acutance S of compensation result as output;Use optimization algorithm search object function S=Obj (Isj
(d),Gj,Tj) obtain maximum value when dispersion parameters Gj,Tj。
The FFT dispersion compensation modules include:
Treat offset data Ij(d) Fast Fourier Transform (FFT) is carried out, its amplitude spectrum A (ω) and phase spectrum are obtainedI.e.Wherein ω is optic angle frequency;Simultaneously dispersion compensation phase is calculated using dispersion parametersWherein ω0For light source center angular frequency;Then the phase after dispersion compensation is calculated
It composes positionLocal offset data I is obtained finally by inverse fast Fourier transform IFFTlocal(d)
=IFFT[A(ω)exp[iφ(ω)]].
Compared with the prior art, the advantages of the present invention are as follows:
1, the present invention can be realized without additional hardware configuration and carry out superelevation distribution birefringence to OCDP output signals
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.
2, not only without this priori of the birefringence dispersion of device under test, but also for those since birefringence is accumulated
Measure it is excessive caused by by the interference peaks that noise floods can also realize birefringence dispersion compensate.
3, it individually handles compared to data to be directly divided into multiple segments, the present invention has faster processing speed.
Description of the drawings
Fig. 1 is OCDP test principle figures;
Fig. 2 is the functional block diagram of superelevation distribution birefringence dispersion compensation;
The emulation signal and its output result after birefringence dispersion compensates that Fig. 3 is OCDP test systems.
Fig. 4 is that have the result that the polarization-maintaining fiber coil measured data of several calibration solder joints carries out dispersion compensation to one.
Specific implementation mode
The present invention is described in more detail below in conjunction with the accompanying drawings:
A kind of dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion.This method is logical
It crosses interception head end data to carry out the chromatic dispersion measurement technology based on optimization and obtain dispersion parameters herein, and is compensated and needed with this
Remaining endian data is carried out above-mentioned steps by offset data again after intercepting and preserving the head end data of perfect compensation,
Until all data are compensated by perfection, the superelevation that is eliminated after finally splicing all of which distribution birefringence dispersion influences
Polarization interference data;Amount of redundancy when this method is divided by the operation of data padding and setting data is come after ensureing compensation
Data are not by end effects, so as to directly carry out head and the tail splicing.The chromatic dispersion measurement technology that this method uses makes dispersion
Parameter determination does not need high s/n ratio interference peaks, and applicability is stronger;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, are finally compensated by perfection.
The present invention includes:
1st, optical coherence domain polarization measurement device is carried out what distribution polarization interference was tested to polarization-maintaining device to be measured
The both ends of distribution polarization interference initial data carry out zero padding and obtain initial data I to be compensated1(d), wherein d >=0 is optical path difference,
And I1(d) in optical path difference range d ∈ [0,Δd1]And d ∈ [dmax-Δd1,dmax]It is inside added data 0, Δ here
d1It is taken as Δ d1=0.5mm, dmaxFor the maximum value of optical path difference d;
2nd, data I to be compensated is extracted1(d) in optical path difference range d ∈ [Δd1,d1+Δd1]Interior data segment is as head end
Data Is1(d), d here1It is taken as d1=4mm;
3rd, to head end data Is1(d) dispersion parameters, i.e. second-order dispersion G are measured using optimization chromatic dispersion measurement technology1With three
Rank dispersion T1, and determine head end data Is1(d) dispersion parameters G is used1With T1Carry out the peak of the result after FFT dispersion compensation modules
It is worth the optical path difference l of position1, then calculate distributed second order dispersion DGj=Gj/lj;
4th, the dispersion parameters G obtained in previous step is utilized1With T1, treat offset data I1(d) FFT dispersion compensations are used
Module carries out dispersion compensation, obtains local offset data Ilocal(d);
5th, to local offset data Ilocal(d) data segmentation is carried out, wherein optical path difference range d ∈ [ are extracted;Δd1,d1+
Δd1]Interior data segment is as the head end data I after compensationc1(d), and wherein optical path difference range d ∈ [ are extracted;d1,dmax-Δd1]
Interior data segment is as the endian data I after compensatione1(d), I herelocal(d) d ∈ [ in;0,Δd1]And d ∈ [dmax-Δd1,
dmax]Part be dropped, and head end data Ic1(d) with endian data Ie1(d) then there is length is optical path difference Δ d1Overlap
Point;
6th, the endian data I after compensation is examinede1(d) whether full -0 foot δ d=d of data lengthmax-d1≥d1If full
Foot, then it is assumed that data length is sufficiently large, then is first depending on range error Δ I caused by the dispersion residual of permissiondB=0.5dB is calculated
One parameter related with dispersionThen d is updated1And Δ d1Respectively
WithWhereinFor wide spectrum light source coherence length, λ0, λ points of Δ
Not Wei light source centre wavelength and full width at half maximum, c be vacuum in the light velocity, and SNR=90 be system dynamic range, then will
Endian data I after compensatione1(d) end carries out zero padding as new data I to be compensated2(d) (zero-filled data length is Δ
d2), and updating optical path difference d makes its minimum value be 0;
7th, step 2 is repeated to 6 (the lower footnote of wherein all variables adds 1), until endian data I after compensatingej(d)
Data length be unsatisfactory for δ d=dmax-d1≥d1Relationship;
8th, by the head end data I after all compensationcj(d) it in turn carries out head and the tail splicing and is fully compensated for data
Iout(d)。
The optimization chromatic dispersion measurement is:
An object function S=Obj (I is constructed firstsj(d),Gj,Tj), with head end data Is1(d) and second-order dispersion
G1With third-order dispersion T1Use FFT dispersion compensation modules to head end data I as input, and using dispersion parameterss1(d) color is carried out
It dissipates compensation to attempt, finally calculates the acutance S of compensation result as output;Secondly, using optimization algorithm search object function S=
Obj(Isj(d),Gj,Tj) obtain maximum value when dispersion parameters Gj,Tj。
The FFT dispersion compensation modules are:
Treat offset data Ij(d) Fast Fourier Transform (FFT) (FFT) is carried out, its amplitude spectrum A (ω) and phase spectrum are obtainedI.e.Wherein ω is optic angle frequency;Dispersion is calculated using dispersion parameters simultaneously to mend
Repay phaseWherein ω0For light source center angular frequency;Then dispersion compensation is calculated
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, which sends out centre wavelength,
Half spectral width is the wide spectrum optical of the Gaussian spectrum shape of 50nm, by one 2:98 1 × 2 coupler 102, wherein 2% light enters
Photodiode 103 is used to monitor light source working state in real time, and the light of residue 98% enters isolator 104 for obstructing back instead
Then light injects device under test 107 via the polarizer 105 and polarizer tail optical fiber 106, is determined according to the type of device under test 107
Determine the polarizer 105 to shaft angle degree, if device under test 105 is polarization-maintaining fiber coil, the polarizer 105 is 0 ° to shaft angle degree, if
Device under test 105 is Y waveguide, then the polarizer 105 is 45 ° to shaft angle degree, and the light of injection device under test 105 can be due to wherein depositing
Defect point and the coupling between fast and slow axis occurs, from the fast and slow axis of device under test 105 export light via analyzer tail optical fiber
108 and analyzer 109,2 × 2 couplers 112 are then injected into, the operation wavelength of 2 × 2 coupler 112 is 1310nm/
1550nm, and output power ratio is 50:The light of 50,2 × 2 coupler, 112 two-arm output respectively enters reference arm 113 and scanning
Arm one end 114, at the same time, the centre wavelength that distributed feedback laser 110 exports are the laser of 1310nm, by isolator
2 × 2 couplers 112 are also passed through after 111 and enter 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 pass through a 116 freedom of entry space of GRIN Lens,
It is returned into GRIN Lens 116 after being reflected by displacement platform speculum 117, it is another into scan arm to again pass by circulator 115
End 118, the laser of the wide spectrum optical and 1310nm wavelength that turn finally to the 1550nm centre wavelengths in arm 113 and scan arm 118 is equal
It being interfered at 2 × 2 couplers 119, interference signal inputs wavelength division multiplexer 120 and 121 respectively by 2 × 2 couplers, 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, finally all detections letter
Number via data acquisition unit 126 acquire, the data acquisition unit 126 include a data collecting card and a computer.
Fig. 2 is the functional block diagram of superelevation distribution birefringence dispersion compensation, and data acquisition unit 126 passes through sync bit moving stage
The movement of speculum 117 and the signal acquisition process of its own are distributed polarization interference output signal as initial data to obtain
201 (wide range optical interference signals) carry out data padding 202 to the head and the tail both ends of initial data 201 and obtain data 203 to be compensated,
Head end data 204 is then extracted to it, and obtains the second order of head end data 204, third-order dispersion using optimization chromatic dispersion measurement 205
206, the position of peak value in the result of dispersion compensation, meter are further carried out to head end data 204 according to second order, third-order dispersion 206
Distributed second order dispersion 207 is calculated, second order, third-order dispersion parameter 206 are then utilized, offset data 203 is treated and uses FFT dispersions
Compensating module 208 carries out dispersion compensation, obtains local offset data 209, carries out data segmentation 210 to it in next step, obtains head
End data 211 and endian data 212 have been thrown aside local offset data head and the tail both ends and have been mended with data during data are divided
The identical data volume of institute's zero padding number in 0 202 operations, and there is also this with endian data 212 for obtained head end data 211
Then the overlapping part of a data volume carries out the whether sufficiently large judgement 213 of data length, if data to endian data 212
Length is 1 times of head end data or more, that is, thinks data long enough, then endian data 212 is carried out data padding 202 and again
Aforementioned 202 to 213 operation is carried out, until judging that 213 obtain result as "No", finally splices all head end datas from beginning to end successively
214 to get to required full remuneration data 215.
To describe the specific implementation step of this method in detail with reference to Fig. 3.It is by computer journey shown in the first width of Fig. 3
The OCDP of sequence emulation tests system output signal, it can be seen that wherein there are 4 interference peaks 301,302,303 and 304, he
Respectively by a degree of dispersive influence, and in fact, if the not presence of birefringence dispersion, the width of this 4 interference peaks
Degree should be consistent;Be shown in the second width of Fig. 3 using this method carry out after an iteration as a result, being specially by Fig. 3 the
Then data head and the tail zero padding shown in one width intercepts the part that head end optical path difference is 0.5~4.5mm and (namely includes interference peaks 301
Part) chromatic dispersion measurement is optimized, it is used to compensate all data in the first width of Fig. 3 after obtaining dispersion parameters, compensation result is gone
It is a bit of rear as shown in the second width of Fig. 3 to fall head and the tail, 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 a interference peaks make 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, carries out second repeatedly
For dispersion compensation, as a result as shown in Fig. 3 third width;Third time iteration dispersion compensation is carried out to data shown in Fig. 3 third width
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,, will be processed when showing the result of each round birefringence dispersion compensation here in order to compare conveniently shown in five width
Signal splices before current signal to be compensated.The result shown in the 5th width of last Fig. 3 can be seen that the method by us
The signal-to-noise ratio of 4 interference peaks, which realizes, after progress birefringence dispersion compensation is obviously improved, under the influence of noise, 4 interference peaks
Amplitude it is almost the same.
Finally, Fig. 4 is illustrated handles a polarization-maintaining fiber coil measured data with several calibration solder joints using this method
Carry out dispersion compensation as a result, 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 superelevation distribution birefringence dispersion, feature exists
In including the following steps:
(1) that optical coherence domain polarization measurement device is carried out the distribution tested of distribution polarization interference to polarization-maintaining device to be measured is inclined
The both ends for crosstalk initial data of shaking carry out zero padding and obtain initial data I to be compensatedj(d), wherein d >=0 is optical path difference, j=1,
And Ij(d) in optical path difference range d ∈ [0,Δd1]And d ∈ [dmax-Δd1,dmax]It is inside added data 0, Δ d1It takes
For Δ d1=0.5mm, dmaxFor the maximum value of optical path difference d;
(2) data I to be compensated is extractedj(d) in optical path difference range d ∈ [Δd1,d1+Δd1]Interior data segment is as head end data
Isj(d), d1It is taken as d1=4mm;
(3) to head end data Isj(d) dispersion parameters, i.e. second-order dispersion G are measured using optimization chromatic dispersion measurement technology1And third-order dispersion
T1, and determine head end data Isj(d) dispersion parameters G is used1With T1Carry out the peak position of the result after FFT dispersion compensation modules
Optical path difference l1, then calculate distributed second order dispersion DGj=Gj/lj;
(4) the dispersion parameters G obtained in previous step is utilized1With T1, treat offset data Ij(d) use FFT dispersion compensation modules into
Row dispersion compensation obtains local offset data Ilocal(d);
(5) to local offset data Ilocal(d) data segmentation is carried out, wherein optical path difference range d ∈ [ are extracted;Δd1,d1+Δd1]It is interior
Data segment as compensation after head end data Icj(d), and wherein optical path difference range d ∈ [ are extracted;d1,dmax-Δd1]Interior number
According to section as the endian data I after compensationej(d), Ilocal(d) d ∈ [ in;0,Δd1]And d ∈ [dmax-Δd1,dmax]Part
It is dropped, and head end data Icj(d) with endian data Iej(d) then there is length is optical path difference Δ d1Overlapping part;
(6) the endian data I after compensation is examinedej(d) whether data length meets δ d=dmax-d1≥d1, if satisfied, then first
Range error Δ I caused by being remained according to the dispersion alloweddB=0.5dB calculates a parameter related with dispersionUpdate d1And Δ d1RespectivelyWith
WhereinFor wide spectrum light source coherence length, λ0, Δ λ is respectively the centre wavelength and half Gao Quan of light source
Width, c are the light velocity in vacuum, and SNR=90 is the dynamic range of system, by the endian data I after compensationej(d) end carries out
Zero padding is as new data I to be compensatedj+1(d), zero-filled data length is Δ d2, updating optical path difference d makes its minimum value be 0;
(7) accordingly increase j=j+1, repeat endian data I of the step (2) to (6) after when compensationej(d) data length is not
Meet δ d=dmax-d1≥d1Relationship;
(8) by the head end data I after all compensationcj(d) it 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 superelevation distribution birefringence dispersion according to claim 1
Compensation method, which is characterized in that the optimization chromatic dispersion measurement includes:
Construct an object function S=Obj (Isj(d),Gj,Tj), with head end data Isj(d) and second-order dispersion G1With three ranks
Dispersion T1Use FFT dispersion compensation modules to head end data I as input, and using dispersion parameterssj(d) dispersion compensation is carried out to taste
Examination finally calculates the acutance of compensation result as output;Use optimization algorithm search object function S=Obj (Isj(d),Gj,
Tj) obtain maximum value when dispersion parameters Gj,Tj。
3. a kind of dispersion for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion according to claim 1
Compensation method, which is characterized in that the FFT dispersion compensation modules include:
Treat offset data Ij(d) Fast Fourier Transform (FFT) is carried out, its amplitude spectrum A (ω) and phase spectrum are obtainedI.e.Wherein ω is optic angle frequency;Simultaneously dispersion compensation phase is calculated using dispersion parametersWherein ω0For light source center angular frequency;Then the phase after dispersion compensation is calculated
It composes 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 CN106525390A (en) | 2017-03-22 |
CN106525390B true 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) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109697270B (en) * | 2019-01-16 | 2022-04-01 | 中国工程物理研究院激光聚变研究中心 | Light beam dispersion characteristic inversion algorithm based on spatial spectrum interference |
CN113804406B (en) * | 2021-08-16 | 2023-12-26 | 广东工业大学 | Birefringent 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 (5)
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 |
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 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8004686B2 (en) * | 2004-12-14 | 2011-08-23 | Luna Innovations Inc. | Compensating for time varying phase changes in interferometric measurements |
-
2016
- 2016-11-09 CN CN201610985383.1A patent/CN106525390B/en active Active
Patent Citations (5)
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 |
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 |
Also Published As
Publication number | Publication date |
---|---|
CN106525390A (en) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1705471B1 (en) | Apparatus for measuring differential mode delay of multimode optical fiber | |
US7728985B2 (en) | Polarization-sensitive common path optical coherence reflectometry/tomography device | |
WO2019029163A1 (en) | Apparatus and method for eliminating polarization fading in ofdr | |
CN102279095B (en) | Device for reducing influence of birefringent chromatic dispersion on polarization coupling measurement of polarization maintaining optical fiber | |
US9041935B2 (en) | Measuring polarization crosstalk in optical birefringent materials and devices based on reduction of line broadening caused by birefringent dispersion | |
CN106525390B (en) | A kind of dispersion compensation method for the optical fibre polarization-maintaining device with superelevation distribution birefringence dispersion | |
CN102928199A (en) | Device and method for improving polarization crosstalk measurement performance of optical device | |
US8300230B2 (en) | System and method to determine chromatic dispersion in short lengths of waveguides using a 3-wave interference pattern and a single-arm interferometer | |
CN103900680B (en) | A kind of device utilizing light source to suppress polarization crosstalk to measure noise and detection method | |
US20070109554A1 (en) | Polarization sensitive optical coherence device for obtaining birefringence information | |
CN107515017A (en) | A kind of optical frequency domain reflectometer of light wave frequency shift modulation | |
CN104792503A (en) | Device for enhancing sensitivity of distributed crosstalk measurement of optical polarizing device | |
Calvani et al. | Polarization measurements on single-mode fibers | |
CN105784336A (en) | Fiber device transmission and reflection performance test device and method | |
Qu et al. | Improvement of strain measurement range via image processing methods in OFDR system | |
CN106989904A (en) | A kind of measuring method of polarization maintaining optical fibre extinction ratio | |
US5654793A (en) | Method and apparatus for high resolution measurement of very low levels of polarization mode dispersion (PMD) in single mode optical fibers and for calibration of PMD measuring instruments | |
CN113804299A (en) | Optical fiber device distributed bidirectional polarization measuring device based on optical frequency domain interference | |
CN111964873B (en) | High-precision distributed extinction ratio measuring method for polarization maintaining optical fiber | |
CN113804412A (en) | Optical fiber device micro-chromaticity dispersion measuring device based on ring light path structure | |
WO2004005973A2 (en) | Polarization diversity detection without a polarizing beam splitter | |
CN113804298B (en) | Distributed bidirectional polarization measurement device based on matching correction light frequency domain interference | |
CN108398098A (en) | A kind of optical surface spacing non-contact measurement apparatus and method | |
CN113804405A (en) | Micro optical fiber dispersion measuring device based on double-coupler ring optical path structure | |
Zhang et al. | Polarization characteristic parameters measurement of Y-waveguide multifunction integrated optic circuit based on noise reduction processing |
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 |