Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Quasi-distributed birefringence dispersion measurement for polarization maintain device with high accuracy based on white light interferometry

Open Access Open Access

Abstract

A white light temporal interferometric technique for measurement of the quasi-distributed birefringence dispersion (BD) in a polarization maintain (PM) device with high accuracy based on weighted least square (WLS) method is presented. It is verified theoretically and experimentally that the accuracy of WLS method and the conventional ordinary least square (OLS) method both are proportional to the signal-to-noise ratio (SNR) of interferogram, whereas the WLS method holds a higher scaling factor because it is more suitable for heteroscedastic model that has unequal error variance. The experiment results show a repeatability of ~4.6 × 10−5 ps/nm @ 1550 nm with WLS method for 100 sets of data, and ~4.3 × 10−4 ps/nm with OLS method, for an interferogram with SNR of 30 dB. Besides, WLS method without iterative operation is carried out by using power spectrum of interferogram as weight value. The feasibility of this technique is demonstrated by distinguishing the quasi-distributed BD of every part for a packaged Y-waveguide with two 1m-long PM pigtails.

© 2016 Optical Society of America

1. Introduction

As an intrinsic Mach-Zehnder interferometer, polarization maintain (PM), or rather high birefringence, components have various applications in the fields of microwave photonic filters [1], interferometric fiber optic sensors [2]. The presence of birefringence dispersion (BD), which is the differential chromatic dispersion of the two eigenmodes, will decrease the sensitivity and spatial resolution for distributed sensor [3, 4 ] and may lead to reflection spectrum distortion for microwave photonic filters [5]. There are two types of white light interferometry (WLI), spectral [6, 7 ] and temporal [8, 9 ], to measure the BD for PM components. Spectral WLI can carry out broad spectral range measurement of BD using supercontinuum source, but it is not suitable for distributed measurement of BD for complicated devices. Temporal WLI that typically requires a Gaussian source, has the capability for distributed measurement by lateral force method for PM fibers [9] or splice points induced crosstalk for other PM components.

There are several temporal WLI to measure the BD for PM components. Okamoto and Hosaka [10] measured the BD of PM fiber by subtracting the chromatic dispersion of its two eigenmodes. This technique employs commercial instrument to obtain chromatic dispersion, which is quite convenient. The accuracy however is not sufficient for PM fiber nowadays with less BD. The technique for BD measurement proposed by Tang et al. [9] holds the advantages of high sensitivity and accuracy, and demonstrates the distributed measurement ability. This technique however is able to obtain the BD merely at central wavelength, while the BD is wavelength dependent. Flavin et al. [8] presented a method for BD measurement that deals with the phase spectrum of interferogram using polynomial fitting method based on ordinary least square (OLS) algorithm. Unavoidable temporal white Gaussian noise will erode the measurement results, and the signal-to-noise ratio (SNR) of interferogram is related to the coupling intensity of splice points [11] which is a fixed number. When making up a packaged PM device, manufacturers expect less coupling intensity [12] of splice points between operating component and its pigtails, which will thus make it difficult for measurement of BD with high accuracy.

In this paper, we propose a novel white light interferometric technique to measure the quasi-distributed BD for a PM device with high accuracy. The noise in conventional temporal WLI based on polynomial fitting technique is analyzed and results show that the variance of phase spectrum error is inversely proportional to the power spectrum of interferogram. Accordingly, the weighted least square (WLS) method, which specially deals with heteroscedastic model [13] that has unequal error variance, should be utilized to estimate the dispersion coefficient from phase spectrum and then the local BD curve for a interferogram could be reconstructed. Besides, using power spectrum of interferogram as weight value instead of employing iterative operation to carry out WLS method, the technique we proposed is more accurate and timesaving. Combining local BD curves of several interferograms, we obtain the quasi-distributed BD of every part for a packaged Y-waveguide with two 1m-long PM pigtails.

2. Theory

2.1 Quasi-distributed birefringence dispersion measurement

The experiment setup for recording of a series of temporal interferograms from which we characterize the quasi-distributed BD in the multifunctional integrated waveguide modulator (MFIWM) is illustrated in Fig. 1 . Passing through a coupler, an isolator and a polarizer in turn, the broadband light will be coupled in the fast and slow axes of the MFIWM under test, respectively. Subsequently, the output light from these two axes will be coupled in a single mode fiber using an analyzer. Finally, a Mach-Zehnder interferometer, which is considered as an tunable birefringent device for simplicity, is exploited to obtain interferograms by varying the optical path difference between the probe and reference arms. The MFIWM under test consists of input extend fiber (PMF1), input pigtail (PMF2), LiNbO3 Y-waveguide chip, output pigtail (PMF3) and output extend fiber (PMF4). These two extend fibers are actually pigtails of the polarizer and the analyzer, respectively. The splice points A, B, C, D between these components all can result in a pair of interferograms because of 1st order polarization coupling. In addition, the polarization extinction capability of Y-waveguide chip will also lead to a pair of interferograms, and the fiber Mach-Zehnder interferometer will arise a singleinterferogram. All these interferograms are shown in Fig. 2 . We can see the full width at half maximum (FWHM) of interferograms C and D are narrow and the FWHM of other interferograms located left or right of them are gradually increasing. This is because of the positions of interferograms C and D are close to the zero BD area, given the BD of interferometer is negative, and negative or positive BD will both broaden the interferograms.

 figure: Fig. 1

Fig. 1 Experiment setup of the measurement system.

Download Full Size | PDF

 figure: Fig. 2

Fig. 2 The output signal of the measurement system. The four pairs of interferograms A, A’, B, B’, C, C’, D, and D’ are first order coupling crosstalk induced by the four splice points A, B, C, and D in Fig. 1, respectively. Polarization extinction capability of the Y-waveguide chip leads to the interferograms Y and Y’. The interferogram M is primarily resulted by the interference of excited mode in the probe arm and in the reference arm. The interferograms between B and Y, as well as interferograms between B’ and Y’, are second order coupling crosstalk which are beyond this paper.

Download Full Size | PDF

One can employ fast Fourier transform (FFT) technique and a simple program to obtain the amplitude spectrum and unwrapped phase spectrum from an interferogram. Ignoring terms more than three order, the phase spectrum φ(ω) can be represented by the Taylor expansion about the source central frequency ω0 as [14]

φ(ω)=φ(ω0)+L[LB(ωω0)+12SBD(ωω0)2+16TBD(ωω0)3]
where ω is the angular frequency, L is the length of device under test, and the dispersion coefficients are the linear birefringence (LB), the second order birefringence dispersion (SBD), and the third order birefringence dispersion (TBD), respectively. The second derivative of phase spectrum φ(ω) with respect to wavelength is the BD of interferogram D(λ)=(2πc/λ2)(d2φ/dλ2)/L, the BD of interferogram therefore can simplicity be expressed with dispersion coefficient as
D(λ)=[SBD+TBD(2πcλ2πcλ0)]2πcλ2
where λ0 is the central wavelength of source and c is the velocity of light in vacuum. Consequently, the quasi-distributed BD of every component in the MFIWM can be calculated as [11]
{ΔDl1=(DMDA)/l1=(DA'DM)/l1,ΔDl2=(DADB)/l2=(DB'DA')/l2,ΔDl3=(DDDC)/l3=(DC'DD')/l3,ΔDl4=(DMDD)/l4=(DD'DM)/l4,ΔDY=(DBDMDY+DC)/lY=(DY'DB'DC'+DM)/lY
where ΔDl1, ΔDl2, ΔDl3 ΔDl4 and ΔDY are the quasi-distributed BD of PMF1, PMF2, PMF3, PMF4 and Y-waveguide chip, respectively; DA, DA, DB, DB, DC, DC, DD, DD, DY, DY, and DM are the local BD of interferograms A, A’, B, B’, C, C’, D, D’, Y, Y’, and M, respectively.

2.2 Dispersion coefficient estimation with weighted least square method

The discrete phase spectrum obtained by means of FFT technique can be expressed as

φ=[φ(ω1)φ(ω2)...φ(ωk)]T
where []T represents to matrix transposition, φ(ωi) is the phase value at differential angular frequency ωj(we replace Ωj=ωjω0 with ωj for simplicity), and k is the sample channel number of effective phase spectrum. Without regard to dispersion with order higher than three, one can employ the WLS algorithm to fit the discrete phase spectrum to 3rd order polynomial curve that expressed as
φ=WωB+ε
here
ω=[1ω1ω12ω131ω2ω22ω231ωkωk2ωk3],W=[w10w20wk],B=[B0B1B2B3]T,ε=[ε1ε2εk]T
where εj(j=0,1,,k) and wj is the error and weight value of every phase channel, respectively. The optimal value of weight value is wj=1/σεj2 [13], where σεj2 is the variance of εj. Bi(i=0,1,2,3) is the i th coefficient of fitting phase curve and represent initial phase, LB, SBD/2, TBD/6, respectively.

If the interferogram is subject to white Gaussian noise, the effective phase spectrum obtained by FFT technique will be eroded. In addition, the variance of phase spectrum error is [15]

σε2(ω)=1S[a(ω)]2
where a(ω) is the normalized amplitude spectrum of interferogram, S is the SNR of [a(ω)]2, i.e. Spectral SNR. The weight value therefore can be set as wj=[a(ωj)]2, or equivalently normalized power spectrum, for simplicity. Moreover, the weighted least squares estimator of dispersion coefficient matrix B is
B^wls=(ωTWω)1ωTWφ
And the covariance matrix of estimator B^wls is [13]
V=1S(ωTWω)1
where the diagonal elements of covariance matrix V, are the variance of Bi in turn, i.e. Vi,i=σBi2. Accordingly, the variance of measurement error for SBD and TBD is σSBD2=4V3,3 and σTBD2=36V4,4, respectively. On the other hand, if OLS method is utilized to estimate the dispersion coefficient matrix B, in other word let the weight value matrix W be equal to unit matrix I, the estimator will be B^ols=(ωTω)1ωTφ, and the corresponding covariance matrix is
Vols=1S(ωTω)1ωTW1ω(ωTω)1
where the meaning of diagonal elements of Vols are the same as Eq. (9). We can know from [16] that the the matrix VolsV is nonnegative definite, therefore all the diagonal elements of matrix VolsV are 0. Accordingly, OLS estimator is no longer a minimum-variance estimator for the unequal variance model described as Eq. (5) when the source spectrum is not ideal flat-top, given the variance of WLS estimator is less than it at least. Whereas the WLS method is equivalent to OLS method when the normalized amplitude spectrum of interferogram is constant, because the weight value matrix degenerate to unit matrix. However, the temporal interferometry typically employs Gaussian source for its symmetry between time and frequency domain. Consequently, the WLS method will be more suitable than OLS method since it is still a minimum-variance estimator for Gaussian source.

3. Experiments and results

3.1 Simulation for dispersion coefficient measurement

For given broadband source and Spectral SNR value S, one can utilizing Eqs. (9) and (10) to calculate the theoretical variance of SBD and TBD for WLS and OLS method, respectively. The measurement error variance of WLS and OLS method are both inversely proportional to spectral SNR, and this relationship is verified by simulation experiment. The interferogram in simulation experiment is with SBD of 100 fs2 and TBD of 104 fs3, which is added by white Gaussian noise with different noise level. The BD of interferogram with noise is characterized by WLS and OLS method, respectively. Figure 3 shows the experiment results measured via 1000-time simulations, which demonstrates the relationship between noise power and the standard deviation of SBD and TBD. Notice that the noise power is in time domain, whereas the S in Eq. (7) is spectral SNR. The noise power is proportional to spectral SNR [15] andwe don’t care about what the exact factor of proportionality is. As depicted in Fig. 3, the WLS method is roughly an order of magnitude improvement than OLS method in standard deviation of TBD and SBD at spectral SNR of 50 dB to 65 dB. When the noise is >50 dB, the OLS method is seriously affected by outliers and no longer meet the linear relationship.

 figure: Fig. 3

Fig. 3 Relationship between noise power and the measurement error of TBD (black) and SBD (blue) with OLS method (dash line) and WLS method (solid line) for 1000 sets simulation data, respectively. Inset: zoom out to show the results of OLS method at high noise power from −45 dB to −30 dB.

Download Full Size | PDF

3.2 Error variance of birefringence dispersion for interferograms

The measurement error variance of BD can be calculated by utilizing Eq. (2) as

σD2(λ)=(2πcλ2)2[4V3,3+36(2πcλ2πcλ0)V4,4],
given the SBD and TBD are statistical independent. The packaged Y-waveguide was measured for 100 times with the technique proposed in this paper. The experiment setup and the MFIWM under test are demonstrated in Fig. 1. The broadband source we used is a super luminescent diode (SLD) emitting at 1550 nm with FWHM of approximately 50 nm. A distributed feedback (DFB) laser is used to eliminate the influence of mechanical vibration of scanning mirror. The 2/98 coupler C1 and photo diode PD1 followed the SLD are used for monitoring the source. The splitting ratio of coupler C2 and C3 both are 50:50, and actually C2 and C3 are wavelength division multiplexers. Photo diode PD2 and PD3 are followed by a differential circuit. The length of PMF1, PMF2, PMF3, PMF4 and Y-waveguide chip in the MFIWM under test are 15.36 m, 1.41 m, 0.71 m, 20.16 m and 0.02 m, respectively. In addition, all the four PM fibers are made by YOFC. The product type of PMF1 and PMF4 are PMF 1310/125-18/25-Y, PMF2 is PMF 1550/125-18/250-Y, and PMF3 is PMF 1550/80-18/165-Y.

Figure 4 demonstrates the measurement repeatability of interferograms M and Y’, and corresponding theoretical error profile with almost the same spectral SNR for each interferogram. Using FFT technique, we obtain the spectral SNR of M and Y’ interferogram of approximate 60 dB and 30 dB, respectively, and the central wavelength of them is about 1552 nm. For interferogram M, the WLS method is roughly a fact of 5 improvement than OLS method from Fig. 4(a) and (b). The improvement less than the simulation value (nearly an order of magnitude) and the large deviations between theoretical and experimental results away central wavelength in Fig. 4(b) perhaps due to the ripple of light source which is not white Gaussian noise (it can be seen from Fig. 2). The region of ripple varies with the properties of the working medium in SLD and its cavity length, and the amplitude of ripple is approximately 60 dB less than the corresponding interferogram [17]. We will possibly obtain better results if utilizing the sidelobe-free SLD in [17]. For interferograms Y and Y’, which are not affected by the ripple of first order interferograms induced by SLD source, the WLS method holds indeed roughly an order of magnitude improvement than OLS method.

 figure: Fig. 4

Fig. 4 Comparison of measurement error variance between theory and 100 sets experimental data. (a) Measurement error of interferogram M with OLS method (red solid line) and corresponding theoretical error curve with spectral SNR of 60 dB (blue dash line), (b) error of M with WLS method and corresponding 60 dB theoretical error, (c) error of Y’ with OLS method and corresponding 30 dB theoretical error, (d) error of Y’ with WLS method and corresponding 30 dB theoretical error.

Download Full Size | PDF

3.3 Measurement results of birefringence dispersion for the interferometer

It has already been stated that the interferogram M is caused by the interferometer, the BD of interferogram M is therefore the BD of the Mach-Zehnder interferometer in Fig. 1 in which the length difference of single-mode fiber for two arms is ~0.12 m. According to Eq. (3), the measurement accuracy of BD for PMF1, PMF4, and Y waveguide is bound up with that for the interferometer. Although the spectral domain (SD) WLI is not appropriate for distributed measurement of BD for complicated devices, we could utilize it to measure the BD of interferometer, and compare the measurement result with the Flavin’s OLS method and our WLS method that are both based on temporal WLI. In order to use the SD method proposed by Lee and Kim in [7] to measure the BD of interferometer, we merely require splicing the point S and E in Fig. 1 directly, without the polarizer and analyzer of course, and employing an optical spectrum analyzer (YOKOGAWA AQ6370C) in stead of the differential signal acquisition device to obtain the channel spectrum, i.e. the spectral interference pattern. We adjust the SNR to ~60 dB via varying the power of light source and measure the interferometer 100 times as well as the temporal methods for the sake of consistency.

Figure 5(a) demonstrates the measurement results of BD for the interferogram with Flavin’s OLS method, WLS method, and Lee’s SD method, in which the maximum deviation of these three methods is less than 5×105 ps/nm that is approximately equivalent to a single- mode fiber with length of 3 mm. Figure 5(b) shows the corresponding measurement errors calculated with 100 sets data, respectively. Actually, the measurement error curves of OLS method and WLS method are identical with the experimental results of interferogram M in Fig. 4(a) and (b). It is demonstrated that the measurement error of WLS method is better than that of SD method. The measurement error of SD method is primarily induced by the phase noise in the process of data acquisition, and the suppression of phase noise is beyond this paper. Whereas, the phase noise induced at the time of data acquisition in WLS method is calibrated with a DFB laser. Besides, the phase noise induced by the amplitude noise of raw data via FFT technique is suppressed by the WLS method.

 figure: Fig. 5

Fig. 5 Comparison of the measurement results of BD for the interferometer. (a) Results of BD for the interferometer measured with two temporal methods and a spectral domain method. (b) The measurement error of these three methods calculated with 100 sets data. The error curves of OLS method and WLS method are actually the same as the experimental error curves in Fig. 4(a) and (b), respectively.

Download Full Size | PDF

3.4 Measurement results of quasi-distributed birefringence dispersion for packaged device

One can obtain the quasi-distributed BD of every component under test from Eq. (3) after measuring the local BD of every labeled interferogram in Fig. 2. It is demonstrated in Fig. 6 that the measurement results for quasi-distributed BD of the four PM fibers. PMF1, PMF3 and PMF4 have the similar dispersion characteristics, whereas the PMF2 much differs from them. We have mentioned before that the zero BD point is located in the vicinity of interferograms C and D, so the interferograms on the right side of M have larger dispersion than those on the left. Figure 6 also compares the results calculated from left and right interferograms, which confirms our proposed method can measure positive and negative BD that quite close to zero. It is shown in Fig. 6(b) that PMF2 has large 3rd order BD, and if we don’t distinguish the BD of PMF2 from Y-waveguide chip, it will greatly influences the measurement result of BD for Y-waveguide chip. Besides, varying the length of pigtails when using the device in practice will change the dispersion characterization of this whole device. Figure 7 depicts the comparison of measurement results for distributed BD of Y-waveguide with and without the two pigtails.

 figure: Fig. 6

Fig. 6 Measurement results for quasi-distributed BD of the four PM fibers. (a) PMF1, PMF3, and PMF4; (b) PMF2. Where, left means the curve is calculated by the data of interferograms on the left hand of interferogram M as shown in Fig. 2, and right means the right hand of M.

Download Full Size | PDF

 figure: Fig. 7

Fig. 7 Measurement results for quasi-distributed BD of Y-waveguide chip with and without two pigtails of it.

Download Full Size | PDF

4. Conclusion

In conclusion, we propose a novel white light interferometric technique for quasi-distributed BD measurement with high accuracy based on WLS method. The higher accuracy than OLS method is because of WLS method is more suitable for unequal error variance model. Furthermore, employing power spectrum of interferogram as weight value is more accurate and cancels iterative operation which is time-consuming. The quasi-distributed BD of a packaged Y-waveguide with two ~1 m-long PM pigtails is measured with this technique. The measurement results show the BD of pigtails are comparable with 0.02 m-long Y-waveguide chip, thus distinguishing them becomes very important.

Acknowledgment

This word is supported by the National Natural Science Foundation of China (61422505, 61227013, 61307104, 61405044); The Program for New Century Excellent Talents in University (NCET-12-0623); National Key Scientific Instrument and Equipment Development Project (No. 2013YQ040815); and the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20122304110022).

References and links

1. W. Muguang and Y. Jianping, “Multitap microwave photonic filter with negative coefficients based on the inherent birefringence in a LiNbO3 phase modulator,” IEEE Photonics J. 5, 5500709 (2013). [CrossRef]  

2. J. Noda, K. Okamoto, and Y. Sasaki, “Polarization-maintaining fibers and their applications,” J. Lightwave Technol. 4(8), 1071–1089 (1986). [CrossRef]  

3. W. Urbanczyk and W. J. Bock, “Analysis of dispersion effects for white-light interferometric fiber-optic sensors,” Appl. Opt. 33(1), 124–129 (1994). [CrossRef]   [PubMed]  

4. T. Xu, W. Jing, H. Zhang, K. Liu, D. Jia, and Y. Zhang, “Influence of birefringence dispersion on a distributed stress sensor using birefringent optical fiber,” Opt. Fiber Technol. 15(1), 83–89 (2009). [CrossRef]  

5. Y. Xiaoke, L. Chao, Y. Xiufeng, Z. Wen-De, W. Fang, D. Lei, and W. Yixin, “Continuously tunable microwave-photonic filter design using high-birefringence linear chirped grating,” IEEE Photonics Technol. Lett. 15(5), 754–756 (2003). [CrossRef]  

6. D. Deng, D. Sega, T. Cheng, W. Gao, X. Xue, T. Suzuki, and Y. Ohishi, “Dispersion characterization of two orthogonal modes in a birefringence tellurite microstructured optical fiber,” Opt. Express 22(20), 23920–23927 (2014). [CrossRef]   [PubMed]  

7. J. Y. Lee and D. Y. Kim, “Versatile chromatic dispersion measurement of a single mode fiber using spectral white light interferometry,” Opt. Express 14(24), 11608–11615 (2006). [CrossRef]   [PubMed]  

8. D. A. Flavin, R. McBride, and J. D. C. Jones, “Dispersion of birefringence and differential group delay in polarization-maintaining fiber,” Opt. Lett. 27(12), 1010–1012 (2002). [CrossRef]   [PubMed]  

9. F. Tang, X. Z. Wang, Y. Zhang, and W. Jing, “Characterization of birefringence dispersion in polarization-maintaining fibers by use of white-light interferometry,” Appl. Opt. 46(19), 4073–4080 (2007). [CrossRef]   [PubMed]  

10. K. Okamoto and T. Hosaka, “Polarization-dependent chromatic dispersion in birefringent optical fibers,” Opt. Lett. 12(4), 290–292 (1987). [CrossRef]   [PubMed]  

11. Y. Jun, Y. Yonggui, Z. Ai, C. Jun, L. Chuang, Y. Dekai, H. Sheng, P. Feng, W. Bing, Z. Yu, L. Zhihai, and Y. Libo, “Full Evaluation of Polarization Characteristics of Multifunctional Integrated Optic Chip With High Accuracy,” J. Lightwave Technol. 32(22), 4243–4252 (2014). [CrossRef]  

12. Z. Wenzin, “Automated fusion-splicing of polarization maintaining fibers,” J. Lightwave Technol. 15(1), 125–134 (1997). [CrossRef]  

13. D. C. Montgomery, E. A. Peck, and G. G. Vining, Introduction to linear regression analysis (John Wiley & Sons, 2012), Chap. 5.

14. T. Grósz, A. P. Kovács, M. Kiss, and R. Szipőcs, “Measurement of higher order chromatic dispersion in a photonic bandgap fiber: comparative study of spectral interferometric methods,” Appl. Opt. 53(9), 1929–1937 (2014). [CrossRef]   [PubMed]  

15. A. Gosteva, M. Haiml, R. Paschotta, and U. Keller, “Noise-related resolution limit of dispersion measurements with white-light interferometers,” J. Opt. Soc. Am. B 22(9), 1868–1874 (2005). [CrossRef]  

16. S. Puntanen, G. P. Styan, and J. Isotalo, Matrix tricks for linear statistical models: our personal top twenty (Springer Science & Business Media, 2011), Chap. 15.

17. J. E. Fouquet, G. R. Trott, W. V. Sorin, M. J. Ludowise, and D. M. Braun, “High-power semiconductor edge-emitting light-emitting diodes for optical low coherence reflectometry,” IEEE J. Quantum Electron. 31(8), 1494–1503 (1995). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1
Fig. 1 Experiment setup of the measurement system.
Fig. 2
Fig. 2 The output signal of the measurement system. The four pairs of interferograms A, A’, B, B’, C, C’, D, and D’ are first order coupling crosstalk induced by the four splice points A, B, C, and D in Fig. 1, respectively. Polarization extinction capability of the Y-waveguide chip leads to the interferograms Y and Y’. The interferogram M is primarily resulted by the interference of excited mode in the probe arm and in the reference arm. The interferograms between B and Y, as well as interferograms between B’ and Y’, are second order coupling crosstalk which are beyond this paper.
Fig. 3
Fig. 3 Relationship between noise power and the measurement error of TBD (black) and SBD (blue) with OLS method (dash line) and WLS method (solid line) for 1000 sets simulation data, respectively. Inset: zoom out to show the results of OLS method at high noise power from −45 dB to −30 dB.
Fig. 4
Fig. 4 Comparison of measurement error variance between theory and 100 sets experimental data. (a) Measurement error of interferogram M with OLS method (red solid line) and corresponding theoretical error curve with spectral SNR of 60 dB (blue dash line), (b) error of M with WLS method and corresponding 60 dB theoretical error, (c) error of Y’ with OLS method and corresponding 30 dB theoretical error, (d) error of Y’ with WLS method and corresponding 30 dB theoretical error.
Fig. 5
Fig. 5 Comparison of the measurement results of BD for the interferometer. (a) Results of BD for the interferometer measured with two temporal methods and a spectral domain method. (b) The measurement error of these three methods calculated with 100 sets data. The error curves of OLS method and WLS method are actually the same as the experimental error curves in Fig. 4(a) and (b), respectively.
Fig. 6
Fig. 6 Measurement results for quasi-distributed BD of the four PM fibers. (a) PMF1, PMF3, and PMF4; (b) PMF2. Where, left means the curve is calculated by the data of interferograms on the left hand of interferogram M as shown in Fig. 2, and right means the right hand of M.
Fig. 7
Fig. 7 Measurement results for quasi-distributed BD of Y-waveguide chip with and without two pigtails of it.

Equations (11)

Equations on this page are rendered with MathJax. Learn more.

φ ( ω ) = φ ( ω 0 ) + L [ LB ( ω ω 0 ) + 1 2 SBD ( ω ω 0 ) 2 + 1 6 TBD ( ω ω 0 ) 3 ]
D ( λ ) = [ SBD + TBD ( 2 π c λ 2 π c λ 0 ) ] 2 π c λ 2
{ Δ D l 1 = ( D M D A ) / l 1 = ( D A ' D M ) / l 1 , Δ D l 2 = ( D A D B ) / l 2 = ( D B ' D A ' ) / l 2 , Δ D l 3 = ( D D D C ) / l 3 = ( D C ' D D ' ) / l 3 , Δ D l 4 = ( D M D D ) / l 4 = ( D D ' D M ) / l 4 , Δ D Y = ( D B D M D Y + D C ) / l Y = ( D Y ' D B ' D C ' + D M ) / l Y
φ = [ φ ( ω 1 ) φ ( ω 2 ) ... φ ( ω k ) ] T
φ = W ω B + ε
ω = [ 1 ω 1 ω 1 2 ω 1 3 1 ω 2 ω 2 2 ω 2 3 1 ω k ω k 2 ω k 3 ] , W = [ w 1 0 w 2 0 w k ] , B = [ B 0 B 1 B 2 B 3 ] T , ε = [ ε 1 ε 2 ε k ] T
σ ε 2 ( ω ) = 1 S [ a ( ω ) ] 2
B ^ wls = ( ω T W ω ) 1 ω T W φ
V = 1 S ( ω T W ω ) 1
V ols = 1 S ( ω T ω ) 1 ω T W 1 ω ( ω T ω ) 1
σ D 2 ( λ ) = ( 2 π c λ 2 ) 2 [ 4 V 3,3 + 36 ( 2 π c λ 2 π c λ 0 ) V 4 , 4 ] ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.