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

Polarization self-compensation in a laser-driven interferometric fiber optic gyroscope with high long-term stability

Open Access Open Access

Abstract

We present a laser-driven interferometric fiber optic gyroscope (IFOG) with polarization self-compensation to achieve high scale-factor stability, sensitivity, and long-term stability. Coherent light with 200kHz linewidth is employed to keep the scale factor stable. The optical scheme ensures polarization reciprocity as well as the optimal working point for good sensitivity. Furthermore, a hybrid machine learning loop (MLL) method, combining the advantages of PID fast response and artificial neural network (ANN) dynamic search, can control a liquid crystal rotator (LCR) to dynamically compensate for slow drift induced by polarization coupling. In open environment, when the sensitivity is 0.005$^{\circ }/\sqrt{\mathrm{h}}$, the bias instability (BI) is significantly optimized from 0.6723°/h at 60s (PID) to 0.3869°/h at 200s (MLL), which is close to the Sagnac interferometric limit (SIL). Such IFOG can meet the real-time and robust requirements for inertial navigation systems in long-term measurement.

© 2024 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Optical rotational sensing has applications in navigation as an optical gyroscope and is usually done with a Sagnac interferometer through phase measurement. At present, interferometric fiber optic gyroscopes (IFOGs) are widely utilized owing to their high sensitivity and portability [15]. In general, the performance of IFOGs can be evaluated in terms of short-term sensitivity and long-term drift [6,7]. Short-term noise directly limits the detection sensitivity. Long-term drift leads to instability. Therefore, high-sensitivity IFOGs with low bias instability (BI) and angular random walk (ARW) have always drawn lots of interest and attention around the world [819], and becoming a goal of many researchers.

Commonly, IFOGs based on the minimal scheme consist of a broadband source, two couplers, a polarizer, and a polarization-maintaining (PM) fiber. It can effectively mitigate the main error sources, namely Kerr-induced drift, Rayleigh backscattering, and polarization error [2023], which significantly promotes short-term sensitivity. However, the broadband source and PM fiber bring some negative effects in long-term drift [24,25]. For one thing, the mean wavelength of the broadband source is difficult to stabilize which yields a worse stable scale factor for the IFOGs [26]. For another thing, the PM fiber can hardly maintain a single polarization in practice due to external environmental factors, such as temperature, magnetic field, and vibration [23,2732], which induces random coupling between two orthogonally polarized light (we refer to it as "polarization coupling"). As a result, there exists an obvious slow drift in the output of IFOGs. Therefore, compared with short-term sensitivity, long-term stability is the core bottleneck of the current engineering application.

In order to reduce slow drift and improve the environmental adaptability of IFOGs, software compensation methods have been widely used. The existing software compensation methods can be classified into two types: Traditional compensation [3335] and artificial intelligence (AI) algorithms [30,3642]. PID control is the most commonly used traditional fitting method in gyroscopes due to its quick response. However, it is difficult for a PID controller [43], whose structure and parameters are fixed in the whole control process, to achieve good performance in a time-varying system. By comparison, AI methods, such as artificial neural networks (ANN) [4345], can fully learn the nonlinear characteristics of the data via using different construction ways, complex network structure, and appropriate training methods, to adaptively search for optimal parameters but usually have to take a lot of time to obtain training samples and construct the model, limiting its application in continuous measurement situations [46,47]. Hence, combining PID control with ANN can dynamically compensate for polarization error and achieve great long-term stability in IFOG.

In this paper, we devote to developing a laser-driven IFOG with polarization self-compensation to achieve high scale-factor stability, sensitivity, and long-term stability. Coherent light with 200$\mathrm {kHz}$ linewidth is employed to keep the scale-factor stability. The combination of a polarization beam splitter (PBS), Faraday rotator (FR), and half-wave plate (HWP) can ensure polarization reciprocity. Benefiting from the combination, the bias phase of an electro-optic modulator (EOM) can be fixed at optimal working point for good sensitivity. Furthermore, a liquid crystal rotator (LCR) is controlled by a hybrid machine learning loop (MLL) method based on an artificial neural network (ANN) and a PID to dynamically compensate for the slow drift induced by external environment fluctuations. The experimental results show that, in the open environment, when the sensitivity is achieved to 0.005$\mathrm {^{\circ }/\sqrt {h}}$, the BI is improved from 0.6723$\mathrm {^{\circ }/h}$ at 60$\mathrm {s}$ using traditional PID method to 0.3869$\mathrm {^{\circ }/h}$ at 200$\mathrm {s}$ via MLL. Moreover, under the same light intensity, the ARW and BI of the current IFOG is nearly equal to the Sagnac interferometric limit (SIL). Such IFOG can meet the real-time and robust requirements for inertial navigation systems in long-term measurement.

2. Principle and theory

2.1 Theoretical model with noise and drift analysis

The schematic diagram is shown in Fig. 1(a). The light is emitted from a distributed feedback laser with 200 $\rm {kHz}$ linewidth and passes a Glan polarizer to generate a vertically polarized light ($y$). Different from these IFOGs based on the conventional minimal scheme, a polarization beam splitter (PBS1), Faraday rotators (FR1 and FR2), and half-wave plates (HWP1 and HWP2) are the key components to keep polarization reciprocity in the current IFOG.

 figure: Fig. 1.

Fig. 1. (a) Experiment setup of the polarization self-compensation IFOG. Laser: distributed feedback laser with 200kHz linewidth, and the central wavelength is 795nm; GL: Glan prism with extinction ratio $\epsilon \thickapprox −50 \mathrm {dB}$; HWP: half-wave plate; PBS: polarization beam splitter with extinction ratio $\epsilon \leq -30 \mathrm {dB}$; FR: Faraday rotator; EOM: electro-optic modulator for shifting phase with half-wave voltage 230V; LCR: liquid crystal rotator for adjusting polarization; PMF: 300m polarization-maintaining fiber, radius 10.1 cm, intrinsic frequency 333kHz; M: mirror; PD1: a photoelectric detector for signal measurement; DAQ: digital acquisition card with $50 \rm {kps}$ sample rate; PID: proportion integration differentiation controller; AWG: arbitrary waveform generator; VC: voltage controller, used to generate high voltage. SUM: addition operator. The red arrow shows the polarization direction of the incident light, while the blue arrow shows the polarization direction of the returned light. (b) Machine learning loop (MLL) schematic diagram. It’s a loop between the fiber optic gyroscope hardware system and the hybrid machine learning algorithm. PD1: Photoelectric detection. The ANN and PID parts are socked together and run separately depending on the cost value. ANN: contains an input layer, an output layer, and five hidden layers with 64 neurons in a single layer. Each hidden layer contains a linear layer and a dropout layer. The dropout layer is a regularization method [48] that prevents the ANN from falling into a local minimum or overfitting during the prediction process. The dropout rate is set to $\eta =0.2$.

Download Full Size | PDF

With the help of a half-wave plate (HWP0), the vertical ($y$) light and the horizontal ($x$) light go through PBS1 with the same powers to ignore the Kerr-induced drift. Then, the vertical ($y$) light and the horizontal ($x$) light propagate in CW and CCW direction respectively. After passing through an electro-optic modulator (EOM), there is a phase difference $\phi _{0}$ between the vertical ($y$) light and the horizontal ($x$) light. Here it can be fixed at an optimal working point to obtain great measurement sensitivity and a linear scale factor for IFOG. In CW direction, after passing through FR1 and HWP1, the direction of the polarization remains unchanged. In CCW direction, with the help of FR2 and HWP2, the direction of polarization rotates counterclockwise by 90$^{\circ }$. In this way, the polarization of CW and CCW directions are both vertical($y$) when coupling into the fiber coil. Benefiting from the optical scheme, the light coupled into orthogonal (unwanted) polarization cannot enter the photodetector (PD1), which partially reduces the noise of the proposed IFOG.

Subsequently, the lights in CW and CCW direction pass through a 300$\rm {m}$ PM fiber coil. However, in practice, the temperature variation makes the PM fiber can hardly maintain the single polarization. As a result, there is a rotation angle $\theta _{n}$ caused by polarization changes. In this case, a liquid crystal rotation (LCR) modulator with the rotation angle $\theta _{c}$ is employed to compensate for $\theta _{n}$. Hence, the lights returning to PBS1 are

$$E_{cw}=A_{0}e^{i\omega _{0}t}e^{i\phi _{cw}}\left( \begin{array}{c} \sin (\theta_{n} +\theta_{c} ) \\ \cos (\theta_{n} +\theta_{c} ) \end{array} \right),$$
$$E_{ccw}=A_{0}e^{i\omega _{0}t}e^{i(\phi _{ccw}+\phi_{0} )}\left( \begin{array}{c} -\cos (\theta_{n} +\theta_{c} ) \\ -\sin (\theta_{n} +\theta_{c} ) \end{array} \right).$$

Here $A_{0}$ and $\omega _{0}$ are the amplitude and the central frequency of the light source. $\phi _{cw}$ and $\phi _{ccw}$ are the induced phase by the Sagnac effect in CW and CCW, respectively. Finally, the horizontal ($x$) light goes through PBS2 and is detected by PD1, which is

$$\begin{aligned} I_{1}(t) &=\left\vert E_{x}^{{\ast} }E_{x}\right\vert \\ &=I_{0}\cos ^{2}(\theta_{c} +\theta_{n} )[1+\cos (\phi_{0} -\phi )]. \end{aligned}$$

Here $\phi =\phi _{cw}-\phi _{ccw}$ is the phase difference related to the rotational angular velocity $\Omega =\frac {c\lambda }{4\pi Lr}\phi$. $c$ is the speed of light in vacuum. $\lambda$ is the laser wavelength. $L$ is the length of the Sagnac loop. $r$ is the radius of the optical fiber loop. Under static conditions $\phi \approx 0$, the interference signal $I_{1}(t)$ is stabilized at $I_{0}=|A_{0}|^2$ , which is defined as the target value. When the phase difference induced by EOM is $\phi _{0}=-\frac {\pi }{2}$ nearby, the bias phase of the IFOG can be fixed at the optimal working point to obtain great measurement sensitivity and a linear scale factor. According to Eq. (3), when $\theta _{c} =\theta _{n}$, the polarization coupling angle $\theta _{n}$ caused by environmental fluctuations can be compensated via the rotation angle $\theta _{c}$ of LCR. The relationship between polarization rotation angle $\theta _{c}$ and LCR modulation amplitude is given in Fig. 2(f). So in the case of small angles, the gyroscope output and the LCR modulation voltage approximately satisfy the cosine function squared relation. Below we briefly introduce how to control LCR via a hybrid machine learning loop (MLL) method to realize polarization self-compensation.

 figure: Fig. 2.

Fig. 2. MLL training and predicted results. (a) Cost versus run numbers. The purple dots are generated by the differential evolution algorithm during iterations, and the red dots are iterated and generated by the neural network. (b) Parameters versus run numbers. The modulation frequency of the LCR [10kHz- 50kHz], which is normalized. The voltage of the LCR [3.75V- 4V], which is normalized, too. It is positively correlated with the polarization rotation angle $\theta _{n}$, but not completely proportional to it. (c) The relationship between the cost and voltage predicted by MLL when the modulation frequency is kept at the optimum (26kHz). (d) The predicted relationship between the cost and modulation frequency when the voltage of LCR is 3.915V. The dashed lines indicate the maximum and minimum boundaries of the ANN prediction. (e) ANN predicted parameters and PID total output vs time. The sampling average interval is 1s. (f) LCR modulation amplitude vs polarization rotation angle $\theta _{c}$.

Download Full Size | PDF

2.2 Hybrid machine learning loop method for polarization self-compensation

Based on the machine-learning package [49], the working principle of the hybrid machine learning loop (MLL) method is shown in Fig. 1(b). The output $I_{1}(t)$ detected by PD1 is sent to the computer via the digital acquisition (DAQ) card. The MLL algorithm program in a computer can actively search the optimal voltage amplitude and modulation frequency for LCR to compensate for the polarization coupling angle $\theta _{n}$ caused by environmental fluctuations. In order to avoid the accumulation of charge drift in the LCR, square-wave modulation with variable parameters is adopted. The whole MLL is divided into 4 processes as follows.

•Initialization. The differential evolution algorithm [50] built-in ANN stochastically generates two sets of parameters, the voltage amplitude and modulation frequency, within the range we set to initialize the IFOG. These parameters are passed through LCR to the IFOG, then the corresponding outputs are fed back to the MLL to optimize the parameters which eventually approach the target value. After the beginning twenty loops, the ANN obtains the initial training data set and constructs the error model between the output values and the target value, also called the $cost$ function, given as

$$cost(t)=|I_{1}(t)-I_{0}|,$$
where $I_{0}$ is the target value discussed in Eq. (3).

•ANN search. We set $k$ as a threshold for $cost(t)$. When $cost(t)$ is larger than the threshold $k$, the ANN starts working. The voltage amplitude and modulation frequency of LCR are continuously iterated by the built-in gradient descent and backpropagation algorithms [51] in ANN. In this process, the ANN can update the voltage amplitude and modulation frequency parameters for the next loop and predict whether the cost reaches the optimal minimum. The cost function updates in the direction of minimizing until it is smaller than the threshold $k$. Meanwhile, in order to find the minimum of the cost function, the mean-square error loss function [52] is adapted in ANN to determine the "distance" between the predicted value $cost_{pre}(t_{i})$ and the measured value $cost(t_{i})$ of the $i$ loop.

$$loss=\frac{1}{N}\sum_{i=1}^{N}[cost_{pre}(t_{i})-cost(t_{i})]^{2},$$
where $N$ is the number of training samples. Besides, to prevent overfitting, the differential evolution iteration is added to the loop every three times shown by the purple dots in Fig. 2(a).

•Fast feedback. As shown in Fig. 1(b), the ANN pauses until it has found the global minimum $cost(t) \rightarrow 0$. The corresponding parameters are the final control values for LCR. To ensure the compensation effect of LCR, the optimal $cost$ is fed to the PID controller (model SIM960) as its locking setpoint via the communication line for a small range of fast feedback. Finally, the outputs of ANN together with PID are fed back to LCR for dynamic adjustment.

•Real-time monitoring. When the ANN completes the prediction in the previous loop, it doesn’t stop working. On the contrary, it continuously monitors the output of our IFOG. When a sudden external event causes the PID to unlock and the condition $cost>k$ is satisfied, a fast search is carried out by ANN to pull back to the optimal parameters. But if PID locks ($cost<k$ is always satisfied at any time ), ANN will identify and receive this PID total output into the training dataset. So the prediction optimal parameters are adaptively updated over time.

Figure 2 shows the training process and prediction results of the MLL. The target value $I_{0}$ is $40.0 {\mathrm{\mu} \mathrm{W}}$. Here, the maximum running number is 80. The threshold $k$ we set is 0.02$\rm {\mu {W}}$, which is usually set as 3 times intensity fluctuations. In Fig. 2(a)(b), the cost function and the parameters, including voltage and frequency, gradually converge to stable values after about 30 iterations. This indicates that ANN has found the best cost and the corresponding optimal parameters. In Fig. 2(c) and (d), the optimal frequency and voltage predicted by ANN are about 26kHz and 3.915V. The corresponding optimal cost is 0.016$\rm {\mu {W}}$. It can be fed to PID as the locking setpoint. The small feedback output of PID is added to the output of the ANN through the adder shown in Fig. 1(a). This is the PID total output in Fig. 2(e). Finally, the voltage of LCR is locked at about 3.915V nearby and the optimal voltage varies with time as shown in Fig. 2(e), which is equivalent to the optimal polarization rotation angle $\theta _{c}=0.40^{\circ }$ nearby. Hence, the optimal voltage and frequency for LCR are adaptively locked by MLL to improve long-term stability.

3. Experimental setup and results

Based on the experimental setup shown in Fig. 1(a), the short-term sensitivity and long-term drift are analyzed. The performance of the current IFOG is analyzed under static conditions and without intrinsic frequency modulation. At a room temperature of $21.0\rm {^{\circ }{C}}$, the average intensity before entering PBS1 is locked at 120$\rm {\mu {W}}$. Intensity stability root meam square $\mathrm {{(RMS)\approx 10^{-4}}}$. The final light intensity at the PD1 is 40.0$\mathrm {\mu {W}}$. With the help of HWP0, the lights go through PBS1 with the same powers. The combination of FR and HWP remains the polarization in CW and CCW direction the same in the PM fiber coil while perpendicularity at PBS1. The output is detected by PD1 and transmitted to the hybrid machine learning algorithm program by DAQ. The updated voltage amplitude and modulation frequency obtained by the ANN are transmitted to AWG. The resultant output light intensity of the proposed IFOG goes back to the software algorithm to predict whether the result is within the given boundary. The optimal minimum cost is passed to PID as a setpoint. The output of PID together with the optimal modulated signal predicted by ANN are added and then fed back to LCR. Finally, LCR compensates for the polarization drift without destroying the reciprocity of the path. The response time of the LCR for $\mathrm {{90^{\circ }}}$ rotation is about 0.15 ms, which can also be replaced by a Lithium niobate electro-optical polarization modulator. In addition, when the proposed IFOG without fiber coil, it is a typical Sagnac interferometer and provides a good benchmark for demonstrating the effect of polarization compensation, called the Sagnac interferometric limit (SIL). Below we analyze the performance of the polarization self-compensation IFOG in detail and compare it with SIL.

3.1 Scale factor

In IFOG, the relation between the angular velocity $\Omega$ and the interference signal $I_1(t)$ can be calibrated via scale factor. Therefore, $\Omega$ can be obtained from $I_1(t)$ via scale factor.

In this experiment, the interference signal $I_{1}(t)$ detected by PD1 is shown in Fig. 3(a). When the EOM is working at $\phi _{0} =-\frac {\pi }{2}$, the interference signal slope is the largest. Under this situation, the IFOG can be fixed at the optimal working point to obtain great measurement sensitivity. When a tiny phase bias is applied to this IFOG, the intensity of the interference signal is shown in Fig. 3(b). Hereby, the relation between the phase shift $\phi$ and the change of output light intensity $\Delta I=I-I_{0}$ satisfies $\phi =K_{1}\Delta I$ with $K_{1}={0.025204\pm 0.000085}\mathrm { \ (rad/\mu W)}$. As $\Omega =\frac {c\lambda }{4\pi Lr}\phi$, the scaling factor between $\Omega$ and $\Delta I$ is $K=K_{1}K_{2}=3.245\times 10^{4} \mathrm {^{\circ }/(h\cdot \mu W)}$. Here $K_{2}=\frac {c\lambda }{4\pi Lr} =1.292\times 10^{6}\mathrm {{^{\circ }{/(h\cdot rad)}}}$ is the conversion factor between $\Omega$ and $\phi$.

 figure: Fig. 3.

Fig. 3. Scale factor calibration of the IFOG. (a) The interference signal of the gyroscope when the triangular wave scanning signal with an amplitude of $460\mathrm {{V_{pp}}}$ and frequency of $0.4\mathrm {Hz}$ drives the EOM. The visibility reaches 99.3%. (b) The calibration factor between the gyroscope output light intensity and phase, where the black dots are the experimental data and the red line is the fitting line.

Download Full Size | PDF

3.2 Short-term measurement sensitivity

To obtain the sensitivity of the proposed IFOG, we employ the relation between signal-to-noise ratio (SNR) and sensitivity as follow [5355]:

$$\delta \Omega =K_{2}\frac{\phi}{\sqrt{10^{\frac{SNR}{10}}B}},$$
where $\delta \Omega$ is the short-term angular velocity sensitivity. $\phi$ is the small phase shift we introduce. The unit of SNR is dB. $B$ is the detection bandwidth. $K_{2}$ is the conversion factor between $\Omega$ and $\phi$ above.

In Fig. 4, the tiny phase shift $\rm {\phi =4.8\times 10^{-4}rad}$ is applied by a sinusoidal voltage ($\mathrm {900Hz, 0.05V_{pp}}$) on EOM. It shows under the same output light intensity of 40.0 $\mathrm { \mu W}$, the noise power spectrum of SIL is 8.9dB lower than the IFOG with a 300m fiber coil, which indicates that the fiber coil indeed introduces extra noise. Besides, the SNRs of IFOG with software compensation are respectively 61.97dB (PID) and 63.04dB (MLL), which are better than the IFOG without feedback (56.49dB). Meanwhile, the MLL algorithm (red line) is about 1dB better than the single PID control (blue line). Based on Eq. (6), the equivalent angular velocity sensitivity is achieved at 0.00582$\mathrm {{^{\circ }/\sqrt {h}}}$ (PID), 0.00513$\mathrm {{^{\circ }/\sqrt {h}}}$ (MLL), and 0.01093$\mathrm {{^{\circ }/\sqrt {h}}}$ (without feedback). Hence, although the utilization of the fiber coil induces extra noise compared with SIL, the sensitivity of our polarization self-compensation IFOG is still maintained at 0.005$\rm {^{\circ }\sqrt {h}}$ level.

 figure: Fig. 4.

Fig. 4. The power spectra obtained by different methods. A tiny phase shift $\mathrm {{\phi =4.8\times 10^{-4}rad}}$ $\mathrm {(900Hz,0.05V_{pp})}$ is introduced through EOM. The four lines are the signal-to-noise ratio(SNR) corresponding to the case PID, MLL, Without feedback, and SIL, respectively, with a resolution bandwidth of 2Hz and output light intensity $40 {\mathrm{\mu} \mathrm{W}}$.

Download Full Size | PDF

3.3 Long-term drift and robustness

In addition to the short-term sensitivity, long-term drift plays a crucial role in inertial navigation systems. At the optimal bias point $\phi _{0} =- \frac {\pi }{2}$, the static output for 2 hours is shown in Fig. 5. The angular velocity drift without software feedback can reach $58.633\mathrm {^{\circ }/h}$. By comparison, the one-hour drifts with PID and MLL can be lower than $5.062\rm {^{\circ }/h}$ and $2.980\rm {^{\circ }/h}$, respectively. Under the same conversion factor $K_{2}=1.292\times 10^{6}\mathrm {{^{\circ }{/(h\cdot rad)}}}$, the equivalent one-hour drift of the angular velocity in Sagnac interferometer is $7.197\rm {^{\circ }/h}$. Here we define the suppression ratio:

$$\gamma =1-Z_{1}/Z_{0},$$
where $Z_{1}$ is the drift with a software algorithm optimization (PID and MLL). $Z_{0}$ is the initial drift without software feedback. The suppression ratio of MLL and PID algorithms reach 91.3% and 94.9%, respectively. The results show that the software algorithms, including PID and MLL, can effectively compensate for the polarization change.

 figure: Fig. 5.

Fig. 5. The long-time output of gyroscope corresponding to different methods. The output interval is 0.25s, and the sample rate is $50\rm {kps}.$

Download Full Size | PDF

Figure 6 shows the Allan variance calculated from the data sampled in Fig. 5. The angular random walk (ARW) and bias instability (BI) obtained by Allan variance are shown in Table 1. Compared with the IFOG without feedback, the utilization of the MLL method optimizes the BI from $11.4256\rm {^{\circ }/h}$ at 120s to $0.3869\rm {^{\circ }/h}$ at 200s and optimizes ARW from $0.2098\rm {^{\circ }/\sqrt {h}}$ to $0.0502 \rm {^{\circ }/\sqrt {h}}$. The corresponding BI suppression rate reaches 96.61%. Besides, both the ARW and BI of the proposed IFOG exceed the performance of a single PID (ARW $0.0540\rm {^{\circ }/\sqrt {h}}$, BI $0.6723\rm {^{\circ }/h}$) and approaches Sagnac interferometric limit (SIL, ARW $0.0454\rm {^{\circ }/\sqrt {h}}$, BI $0.3883\rm {^{\circ }/h}$), which is measured by the proposed IFOG without fiber coil but with the same conversion factor $K_{2}=1.292\times 10^{6}\mathrm {{^{\circ }{/(h\cdot rad)}}}$. Hence, the MLL method effectively decreases the long-term slow drift. It is important to point out that the optimization results are based on the current experimental setup. If the $K_{2}$ factor is increased or the laser linewidth is larger in the future, the optimization limit of the MLL method will be further raised [23,24].

 figure: Fig. 6.

Fig. 6. Allan variance. The logarithmic coordinate is adopted. The output time interval is 0.25s.

Download Full Size | PDF

Tables Icon

Table 1. Allan variance results when the IFOG is in static sensing.

Furthermore, the robustness is also a very important indicator for IFOG. In this experiment, we artificially apply a $0.13 \rm {rad}$ phase bias to EOM at about $130\rm {s}$ to simulate a sudden vibration. The experimental results are shown in Fig. 7. The recovery time of MLL is only about 22s. In contrast, the oscillation duration of PID is much longer than 180s and is difficult to reset. Hence, the MLL method, combining the advantages of PID fast response and ANN dynamic search, can control the LCR to dynamically compensate for slow drift induced by polarization coupling. In actual testing, the recovery time could be further shortened by reducing the number of iterations or narrowing the range of parameters.

 figure: Fig. 7.

Fig. 7. Robustness comparison under the sudden external fluctuation. The red line represents the MLL method, and the blue line corresponds to the single PID method. Sample rate 25 points per second. The red part of the oscillation is caused by the parameter search and iteration of the ANN. The blue part of the oscillation is PID oscillation caused by a sudden external interference.

Download Full Size | PDF

Therefore, although the PM fiber coil introduces extra drift, the proposed polarization Self-compensation IFOG can effectively compensate for slow drift while ensuing short-term sensitivity at the same time.

4. Conclusion

In conclusion, we propose a polarization self-compensation IFOG to achieve high scale-factor stability, sensitivity, and long-term stability. Coherent light with 200kHz linewidth is employed to keep the scale factor stable. The optical scheme together with the software algorithm decreases the intensity fluctuation induced by polarization coupling. Here the combination of polarization beam splitter (PBS), Faraday rotator (FR), and half-wave plate (HWP) allow the polarizations to maintain consistency in the PM fiber coil, which ensures the polarization reciprocity. Furthermore, an LCR modulator is controlled via the hybrid machine learning loop (MLL) method based on the ANN and PID to compensate for the polarization coupling drift dynamically. Experiment results show that the proposed IFOG can effectively improve the long-term stability with good short-term sensitivity ($\mathrm {0.005^{\circ }/\sqrt {h}}$). Compared with the single PID method, the BI of our IFOG is reduced from $\rm {0.6723^{\circ }/h}$ at 60s to $\rm {0.3869^{\circ }/h}$ at 200s. More importantly, the ARW and BI of the current IFOG are close to the Sagnac interferometric limit (SIL). The good long-term stability and robustness of this IFOG are extremely useful in complex and volatile environments, as well as suitable for inertial navigation systems and fiber-optic systems.

Funding

National Natural Science Foundation of China (12104161, U23A2075, 12274132, 12174409); E-Institutes of Shanghai Municipal Education Commission (202101070008E00099); Science and Technology Commission of Shanghai Municipality (19YF1414300); China Postdoctoral Science Foundation (2023M741187); Fundamental Research Funds for the Central Universities.

Disclosures

The authors declare no conflicts of interest.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

References

1. C. K. Kao, Optical fiber systems: technology, design, and applications (McGraw-Hill New York, 1982).

2. E. J. Post, “Sagnac effect,” Rev. Mod. Phys. 39(2), 475–493 (1967). [CrossRef]  

3. K.-X. Sun, M. M. Fejer, E. Gustafson, et al., “Sagnac Interferometer for Gravitational-Wave Detection,” Phys. Rev. Lett. 76(17), 3053–3056 (1996). [CrossRef]  

4. H. C. Lefèvre, “The fiber-optic gyroscope, a century after sagnac’s experiment: The ultimate rotation-sensing technology?” C. R. Phys. 15(10), 851–858 (2014). [CrossRef]  

5. R. Ulrich, “Fiber-optic rotation sensing with low drift,” Opt. Lett. 5(5), 173–175 (1980). [CrossRef]  

6. B. Culshaw, “The optical fibre sagnac interferometer: an overview of its principles and applications,” Meas. Sci. Technol. 17(1), R1–R16 (2006). [CrossRef]  

7. Y. Cao, L. Zhu, and F. Shi, “Dual-polarization interferometric fiber optic gyroscope with shupe effect compensation,” Appl. Phys. Lett. 123(1), 1 (2023). [CrossRef]  

8. M. Mehmet, T. Eberle, and S. Steinlechner, “Demonstration of a quantum-enhanced fiber sagnac interferometer,” Opt. Lett. 35(10), 1665–1667 (2010). [CrossRef]  

9. M. R. Grace, C. N. Gagatsos, Q. Zhuang, et al., “Quantum-enhanced fiber-optic gyroscopes using quadrature squeezing and continuous-variable entanglement,” Phys. Rev. Appl. 14(3), 034065 (2020). [CrossRef]  

10. T. Kim, M. Fiorentino, and F. N. Wong, “Phase-stable source of polarization-entangled photons using a polarization sagnac interferometer,” Phys. Rev. A 73(1), 012316 (2006). [CrossRef]  

11. M. Fink, F. Steinlechner, and J. Handsteiner, “Entanglement-enhanced optical gyroscope,” New J. Phys. 21(5), 053010 (2019). [CrossRef]  

12. W. Burns, C.-L. Chen, and R. Moeller, “Fiber-optic gyroscopes with broad-band sources,” J. Lightwave Technol. 1(1), 98–105 (1983). [CrossRef]  

13. Y. Wu, J. Guo, and X. Feng, “Atom-light hybrid quantum gyroscope,” Phys. Rev. Appl. 14(6), 064023 (2020). [CrossRef]  

14. D. T. Mead and S. Mosor, “Progress with interferometric fiber optic gyro at honeywell,” in Optical Waveguide and Laser Sensors, vol. 11405 (SPIE, 2020), pp. 7–22.

15. G. T. Schmidt and R. E. Phillips, “Ins/gps integration architecture performance comparisons,” Low-Cost Navig. Sens. Integr. Technol 1, 1–21 (2011).

16. X. Wang, C. He, and Z. Wang, “Method for suppressing the bias drift of interferometric all-fiber optic gyroscopes,” Opt. Lett. 36(7), 1191–1193 (2011). [CrossRef]  

17. J. Dong and R. Lu, “Dual-loop sagnac interferometer with a geometric phase shifter for quadrature phase bias locking,” Opt. Lett. 44(22), 5422–5425 (2019). [CrossRef]  

18. A. Rickman, “The commercialization of silicon photonics,” Nat. Photonics 8(8), 579–582 (2014). [CrossRef]  

19. A. A. Berni, T. Gehring, and B. M. Nielsen, “Ab initio quantum-enhanced optical phase estimation using real-time feedback control,” Nat. Photonics 9(9), 577–581 (2015). [CrossRef]  

20. S. W. Lloyd, V. Dangui, and M. J. Digonnet, “Measurement of reduced backscattering noise in laser-driven fiber optic gyroscopes,” Opt. Lett. 35(2), 121–123 (2010). [CrossRef]  

21. G. Pavlath and H. J. Shaw, “Birefringence and polarization effects in fiber gyroscopes,” Appl. Opt. 21(10), 1752–1757 (1982). [CrossRef]  

22. S. W. Lloyd, M. J. Digonnet, and S. Fan, “Modeling coherent backscattering errors in fiber optic gyroscopes for sources of arbitrary line width,” J. Lightwave Technol. 31(13), 2070–2078 (2013). [CrossRef]  

23. J. N. Chamoun and M. J. Digonnet, “Noise and bias error due to polarization coupling in a fiber optic gyroscope,” J. Lightwave Technol. 33(13), 2839–2847 (2015). [CrossRef]  

24. S. W. Lloyd, S. Fan, and M. J. Digonnet, “Experimental observation of low noise and low drift in a laser-driven fiber optic gyroscope,” J. Lightwave Technol. 31(13), 2079–2085 (2013). [CrossRef]  

25. S. W. Lloyd, M. J. Digonnet, and S. Fan, “Near shot-noise limited performance of an open loop laser-driven interferometric fiber optic gyroscope,” in 21st International Conference on Optical Fiber Sensors, vol. 7753 (SPIE, 2011), pp. 1462–1465.

26. J. Chamoun and M. J. Digonnet, “Aircraft-navigation-grade laser-driven fog with gaussian-noise phase modulation,” Opt. Lett. 42(8), 1600–1603 (2017). [CrossRef]  

27. S. Carrara, B. Y. Kim, and H. Shaw, “Bias drift reduction in polarization-maintaining fiber gyroscope,” Opt. Lett. 12(3), 214–216 (1987). [CrossRef]  

28. E. C. Kintner, “Polarization control in optical-fiber gyroscopes,” Opt. Lett. 6(3), 154–156 (1981). [CrossRef]  

29. I. A. Andronova and G. B. Malykin, “Physical problems of fiber gyroscopy based on the sagnac effect,” Phys.-Usp. 45(8), 793–817 (2002). [CrossRef]  

30. C. Zhang, S. Du, J. Jin, et al., “Thermal analysis of the effects of thermally induced nonreciprocity in fiber optic gyroscope sensing coils,” Optik 122(1), 20–23 (2011). [CrossRef]  

31. K. Zhou, S. Pan, S. Liu, et al., “Fiber gyroscope with a double sensitivity employing a polarization splitter,” Opt. Lett. 38(8), 1337–1339 (2013). [CrossRef]  

32. Z. Wang, Y. Yang, and P. Lu, “Dual-polarization interferometric fiber-optic gyroscope with an ultra-simple configuration,” Opt. Lett. 39(8), 2463–2466 (2014). [CrossRef]  

33. R. Peesapati, S. L. Sabat, and K. Karthik, “Efficient hybrid kalman filter for denoising fiber optic gyroscope signal,” Optik 124(20), 4549–4556 (2013). [CrossRef]  

34. R. M. Rajulapati, J. Nayak, and S. Naseema, “Modeling and simulation of signal processing for a closed loop fiber optic gyro’s using FPGA,” Int. J. Eng. Sci & Techn 4, 947–959 (2012).

35. M. A. Johnson and M. H. Moradi, PID control (Springer, 2005).

36. X. Chen, R. Song, C. Shen, et al., “Application of a genetic algorithm elman network in temperature drift modeling for a fiber-optic gyroscope,” Appl. Opt. 53(26), 6043–6050 (2014). [CrossRef]  

37. S. Li, X. Zuo, Z. Li, et al., “Applying deep learning to continuous bridge deflection detected by fiber optic gyroscope for damage detection,” Sensors 20(3), 911 (2020). [CrossRef]  

38. D. Silver, A. Huang, and C. J. Maddison, “Mastering the game of go with deep neural networks and tree search,” Nature 529(7587), 484–489 (2016). [CrossRef]  

39. Y. Cao, W. Xu, and B. Lin, “Long short-term memory network of machine learning for compensating temperature error of a fiber optic gyroscope independent of the temperature sensor,” Appl. Opt. 61(28), 8212–8222 (2022). [CrossRef]  

40. X. Chen and C. Shen, “Study on temperature error processing technique for fiber optic gyroscope,” Optik 124(9), 784–792 (2013). [CrossRef]  

41. W. Zhang, D. Zhang, P. Zhang, et al., “A new fusion fault diagnosis method for fiber optic gyroscopes,” Sensors 22(8), 2877 (2022). [CrossRef]  

42. C. Shen, H. Cao, and J. Li, “Hybrid de-noising approach for fiber optic gyroscopes combining improved empirical mode decomposition and forward linear prediction algorithms,” Rev. Sci. Instrum. 87(3), 1 (2016). [CrossRef]  

43. D. Pogorelaya, M. Smolovik, and S. Volkovskiy, “Adjustment of pid controller in fiber-optic gyro feedback loop,” Gyroscopy Navig. 8(3), 235–239 (2017). [CrossRef]  

44. V. K. Ojha, A. Abraham, and V. Snášel, “Metaheuristic design of feedforward neural networks: A review of two decades of research,” Eng. Appl. Artif. Intell. 60, 97–116 (2017). [CrossRef]  

45. A. Venketeswaran, N. Lalam, and J. Wuenschell, “Recent advances in machine learning for fiber optic sensor applications,” Adv. Intell. Syst. 4(1), 2100067 (2022). [CrossRef]  

46. G. Carleo, I. Cirac, and K. Cranmer, “Machine learning and the physical sciences,” Rev. Mod. Phys. 91(4), 045002 (2019). [CrossRef]  

47. J. Fang, R. Yin, and X. Lei, “An adaptive decoupling control for three-axis gyro stabilized platform based on neural networks,” Mechatronics 27, 38–46 (2015). [CrossRef]  

48. N. Srivastava, G. Hinton, A. Krizhevsky, et al., “Dropout: a simple way to prevent neural networks from overfitting,” The journal of machine learning research 15, 1929–1958 (2014).

49. P. B. Wigley, P. J. Everitt, and A. van den Hengel, “Fast machine-learning online optimization of ultra-cold-atom experiments,” Sci. Rep. 6(1), 25890 (2016). [CrossRef]  

50. J. Chen, Z. Wu, and G. Bao, “Design of coaxial coils using hybrid machine learning,” Rev. Sci. Instrum. 92(4), 1 (2021). [CrossRef]  

51. R. H. Byrd, P. Lu, J. Nocedal, et al., “A limited memory algorithm for bound constrained optimization,” SIAM J. Sci. Comput. 16(5), 1190–1208 (1995). [CrossRef]  

52. K. P. Körding and D. M. Wolpert, “The loss function of sensorimotor learning,” Proc. Natl. Acad. Sci. 101(26), 9839–9842 (2004). [CrossRef]  

53. A. J. Stevenson, M. B. Gray, H.-A. Bachor, et al., “Quantum-noise-limited interferometric phase measurements,” Appl. Opt. 32(19), 3481–3493 (1993). [CrossRef]  

54. S. Liu, Y. Chen, and J. Jiang, “Quantum enhanced electro-optic sensor for e-field measurement,” Opt. Express 29(21), 32865–32874 (2021). [CrossRef]  

55. W. Huang, X. Liang, and B. Zhu, “Protection of noise squeezing in a quantum interferometer with optimal resource allocation,” Phys. Rev. Lett. 130(7), 073601 (2023). [CrossRef]  

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

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. (a) Experiment setup of the polarization self-compensation IFOG. Laser: distributed feedback laser with 200kHz linewidth, and the central wavelength is 795nm; GL: Glan prism with extinction ratio $\epsilon \thickapprox −50 \mathrm {dB}$; HWP: half-wave plate; PBS: polarization beam splitter with extinction ratio $\epsilon \leq -30 \mathrm {dB}$; FR: Faraday rotator; EOM: electro-optic modulator for shifting phase with half-wave voltage 230V; LCR: liquid crystal rotator for adjusting polarization; PMF: 300m polarization-maintaining fiber, radius 10.1 cm, intrinsic frequency 333kHz; M: mirror; PD1: a photoelectric detector for signal measurement; DAQ: digital acquisition card with $50 \rm {kps}$ sample rate; PID: proportion integration differentiation controller; AWG: arbitrary waveform generator; VC: voltage controller, used to generate high voltage. SUM: addition operator. The red arrow shows the polarization direction of the incident light, while the blue arrow shows the polarization direction of the returned light. (b) Machine learning loop (MLL) schematic diagram. It’s a loop between the fiber optic gyroscope hardware system and the hybrid machine learning algorithm. PD1: Photoelectric detection. The ANN and PID parts are socked together and run separately depending on the cost value. ANN: contains an input layer, an output layer, and five hidden layers with 64 neurons in a single layer. Each hidden layer contains a linear layer and a dropout layer. The dropout layer is a regularization method [48] that prevents the ANN from falling into a local minimum or overfitting during the prediction process. The dropout rate is set to $\eta =0.2$.
Fig. 2.
Fig. 2. MLL training and predicted results. (a) Cost versus run numbers. The purple dots are generated by the differential evolution algorithm during iterations, and the red dots are iterated and generated by the neural network. (b) Parameters versus run numbers. The modulation frequency of the LCR [10kHz- 50kHz], which is normalized. The voltage of the LCR [3.75V- 4V], which is normalized, too. It is positively correlated with the polarization rotation angle $\theta _{n}$, but not completely proportional to it. (c) The relationship between the cost and voltage predicted by MLL when the modulation frequency is kept at the optimum (26kHz). (d) The predicted relationship between the cost and modulation frequency when the voltage of LCR is 3.915V. The dashed lines indicate the maximum and minimum boundaries of the ANN prediction. (e) ANN predicted parameters and PID total output vs time. The sampling average interval is 1s. (f) LCR modulation amplitude vs polarization rotation angle $\theta _{c}$.
Fig. 3.
Fig. 3. Scale factor calibration of the IFOG. (a) The interference signal of the gyroscope when the triangular wave scanning signal with an amplitude of $460\mathrm {{V_{pp}}}$ and frequency of $0.4\mathrm {Hz}$ drives the EOM. The visibility reaches 99.3%. (b) The calibration factor between the gyroscope output light intensity and phase, where the black dots are the experimental data and the red line is the fitting line.
Fig. 4.
Fig. 4. The power spectra obtained by different methods. A tiny phase shift $\mathrm {{\phi =4.8\times 10^{-4}rad}}$ $\mathrm {(900Hz,0.05V_{pp})}$ is introduced through EOM. The four lines are the signal-to-noise ratio(SNR) corresponding to the case PID, MLL, Without feedback, and SIL, respectively, with a resolution bandwidth of 2Hz and output light intensity $40 {\mathrm{\mu} \mathrm{W}}$.
Fig. 5.
Fig. 5. The long-time output of gyroscope corresponding to different methods. The output interval is 0.25s, and the sample rate is $50\rm {kps}.$
Fig. 6.
Fig. 6. Allan variance. The logarithmic coordinate is adopted. The output time interval is 0.25s.
Fig. 7.
Fig. 7. Robustness comparison under the sudden external fluctuation. The red line represents the MLL method, and the blue line corresponds to the single PID method. Sample rate 25 points per second. The red part of the oscillation is caused by the parameter search and iteration of the ANN. The blue part of the oscillation is PID oscillation caused by a sudden external interference.

Tables (1)

Tables Icon

Table 1. Allan variance results when the IFOG is in static sensing.

Equations (7)

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

E c w = A 0 e i ω 0 t e i ϕ c w ( sin ( θ n + θ c ) cos ( θ n + θ c ) ) ,
E c c w = A 0 e i ω 0 t e i ( ϕ c c w + ϕ 0 ) ( cos ( θ n + θ c ) sin ( θ n + θ c ) ) .
I 1 ( t ) = | E x E x | = I 0 cos 2 ( θ c + θ n ) [ 1 + cos ( ϕ 0 ϕ ) ] .
c o s t ( t ) = | I 1 ( t ) I 0 | ,
l o s s = 1 N i = 1 N [ c o s t p r e ( t i ) c o s t ( t i ) ] 2 ,
δ Ω = K 2 ϕ 10 S N R 10 B ,
γ = 1 Z 1 / Z 0 ,
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.