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

Double closed-loop control of integrated optical resonance gyroscope with mean-square exponential stability

Open Access Open Access

Abstract

A new double closed-loop control system with mean-square exponential stability is firstly proposed to optimize the detection accuracy and dynamic response characteristic of the integrated optical resonance gyroscope (IORG). The influence mechanism of optical nonlinear effects on system detection sensitivity is investigated to optimize the demodulation gain, the maximum sensitivity and the linear work region of a gyro system. Especially, we analyze the effect of optical parameter fluctuation on the parameter uncertainty of system, and investigate the influence principle of laser locking-frequency noise on the closed-loop detection accuracy of angular velocity. The stochastic disturbance model of double closed-loop IORG is established that takes the unfavorable factors such as optical effect nonlinearity, disturbed disturbance, optical parameter fluctuation and unavoidable system noise into consideration. A robust control algorithm is also designed to guarantee the mean-square exponential stability of system with a prescribed H performance in order to improve the detection accuracy and dynamic performance of IORG. The conducted experiment results demonstrate that the IORG has a dynamic response time less than 76us, a long-term bias stability 7.04°/h with an integration time of 10s over one-hour test, and the corresponding bias stability 1.841°/h based on Allan deviation, which validate the effectiveness and usefulness of the proposed detection scheme.

© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Optical gyros based on Sagnac effect are the key components of inertial measurement units which are widely used in the fields of inertial navigation and attitude determination [1,2]. Especially, the integrated optical resonance gyroscopes (IORGs) have many advantages such as high robustness, theoretical sensitivity and superior reliability due to its inherent characteristics of miniaturized structure, all-solid-state and the combination between integrated optic fabrication and MEMS fabrication [3,4]. Thus, it has been considered as the next generation of micro-sized optical gyro and a promising candidate in the field of inertial navigation [5,6]. However, the numerous noise resulted from complex mechanism of IORG causes the drift of gyro bias and influences the detection accuracy of gyro [7–9]. Therefore, the closed-loop detection technology is important to suppress noise and disturbance and further to optimize the performance of IORG.

Recently, a great deal of attention has been paid to closed-loop detection methods with the aim of achieving precise and reliable IORGs [10–16]. In order to lock the output frequency of laser, Hotate et al. proposed a feedback method that utilized digital serrodyne modulation to track the resonant frequency by LiNbO3 integrated phase modulator [10]. However, the frequent resets of serrodyne waves bring high-frequency harmonics and deteriorate waveform quality of resonator output. The triangular phase modulation method was also designed to control the laser frequency and effectively suppress the backreflection error in the IORGs [12,13]. Feng et al. described a low-delay, high-bandwidth frequency-locking loop implemented on a single filed-programmable gate array with triangular phase modulation, and the signal processing delay-time was reduced to less than 1μs and the bandwidth of 10 kHz was obtained under the optimized loop parameters [14]. In order to reduce backscatter-induced errors, the integer period sampling method [15] and the three lasers scheme [16] using two phase-locked lasers were proposed to improve the scale factor of resonant optic gyro, respectively. These closed-loop detection schemes are achieved by locking the resonant frequency of one propagating lightwave at light source frequency, but the closed-loop tracking of angular velocity is not implemented which consequently leads to IORGs with the inherent limited linearity and dynamic range as well as the sensitivities of existing disturbance. Recently, a double closed-loop IORG is experimentally demonstrated by using a controlled hybrid digital phase modulation scheme [17]. However, we found that the influence mechanism of locking-frequency accuracy from laser frequency locking loop (LFLL) on main closed-loop performance is neglected, which is a challenging problem that confines the engineering process of IORG. Furthermore, other influence factors on angular rate tracking precision are unclear and neglected in the related researches on IORG. Therefore, one of our motivations behind this work is to propose a double closed-loop detection method for IORG to accurately extract the angular velocity while various noise sources are substantially suppressed.

In this work, a novel double closed-loop control scheme with mean-square exponential stability is proposed to optimize detection precision and dynamic performance of the integrated optical resonance gyroscope. Firstly, we investigate the influence mechanism of the nonlinearity of optical effect on detection sensitivity, and the maximum sensitivity of a gyro system is obtained through the optimization of modulation parameters. Secondly, we present the influence theory of laser locking-frequency noise on the performance of angular velocity tracking loop (AVTL), and obtain the stochastic disturbance model of the double closed-loop IORG. Thirdly, a control algorithm based on the double closed-loop model is proposed to ensure that the closed-loop IORG has mean-square exponential stability with a prescribed H∞ performance. Finally, the relevant experiment results are given to demonstrate the effectiveness of proposed closed-loop control method.

2. Problem description

The schematic diagram of the double closed-loop system of IORG with triangular phase modulation is shown in Fig. 1. When the IORG is in a rotary state, the CW and CCW lightwaves propagate different cycles in the resonant cavity, which brings in different resonant frequencies [12].

 figure: Fig. 1

Fig. 1 The principle scheme of the double closed-loop IORG based on Sagnac effect. The triangular wave for phase modulation is applied on the upper arm and bottom arm of the IOPM. The feedback sawtooth wave is applied on the bottom arm of the IOPM. Meanwhile, the phase modulation triangular wave and digital feedback sawtooth wave are all differentially applied on the arms of the IOPM. The IOPM achieves the phase modulation and also enables the closed-loop control of AVTL.

Download Full Size | PDF

The rotation angular velocity is obtained by detecting the resonant frequency difference ΔfΩ. Proper laser frequency is chosen and locked at resonant frequency of CCW path through laser frequency locking loop (LFLL). Meanwhile, the frequency difference ΔfΩ can be demodulated through the difference between the demodulated signal of two loops, and compensated by optical phase modulator in the AVTL to lock the resonant frequency of CW path. The data output of AVTL controller is proportional to the slope coefficient of stair-like digital sawtooth wave imposed on IOPM. And the control update time of controller is constant, which is same as the period of phase modulation triangular wave. Therefore, the data output of AVTL controller and slope of sawtooth wave are equal to the frequency difference ΔfΩ.

The frequency difference ΔfΩ is proportional to the rotation angular velocity Ω according to the principle of Sagnac effect and is taken as the output of an IORG system, which can be presented as [18]

ΔfΩ=DneλΩ
where D is the diameter of ring resonator, ne is the refractive index of resonator host material, λ is the central wavelength of light, c is the speed of light in free space, and Ω is the rotation angular velocity.

However, the interference light intensity representing the difference frequency information is a weak signal containing numerous noise in practical applications [7], which is sensitive to external disturbance and unavoidable optical parameter fluctuation. Meanwhile, the detection accuracy of the IORG system often suffers from the frequency-locking noise of LFLL, but its influence mechanism has rarely been analyzed and reported yet. In this work, the double closed-loop detection model of IORG is analyzed that includes the laser frequency locking and angular velocity tracking. Besides, a robust control algorithm with mean-square exponential stability is also proposed to optimize the measurement accuracy and dynamic performance of IORGs and to promote the process of their engineering applications.

3. Analysis of closed-loop detection system

3.1 Optimization of sensitivity level

Because the resonant frequency difference containing the rotational angular velocity information is a weak signal, the high-frequency carrier technology is adopted in IORG to suppress noise based on weak signal detection theory. Since the interference light intensity is a nonlinear function of resonant frequency, we first analyze the optimal modulation frequency to achieve maximum signal to noise (SNR) and linear work region of IORG.

In the IORG system depicted in Fig. 1, the interference intensity detected by the photoelectric detector can be written as [19]

IT=I0ηβR1+q22qcos2πfτ
where I0 stands for the light intensity of laser source, f is the central frequency of lightwave outputting from the resonator, η is the total loss coefficient of optical path, β and R are the conversion efficiency and trans-impedance of photoelectric detector, τ is the resonator’s round-trip time given by τ = πneD/c, and q is the loss coefficient of ring resonator defined by F=π/cos1(2q/(1+q2)) [19] with F the resonator finesse in the experimental test.

As shown in Fig. 2(a), when the IORG is in the state of rotation, the output interference intensity of resonator IT, is a function of light frequency. The resonant frequency containing information of rotational angular velocity is delivered into the high frequency carrier by triangular phase modulation shown as Fig. 2(b). The relationship between the resonator output and frequency of triangular wave fT can be analyzed by the equivalent frequency bias method, where the equivalent square modulation frequency bias are ± fbias. Then, the time-domain output of photoelectric detector is described in Fig. 2(c), and the closed error of tracking angular velocity is synchronizing demodulated by a digital square wave with same frequency and phase as those of modulation wave. Thus, the demodulated signal can be expressed as

Id(f)=I0ηβRKqd1+q22qcos2π(f+fbias)τI0ηβRKqd1+q22qcos2π(ffbias)τ=2msin2πfbiasτ1+m22cos4πfbiasτ+m22cos4πfτ2mcos2πfτcos2πfbiasτI0ηβRKqd1+q2sin(2πfτ)
where Kqd is the circuit gain of the forward channel in the closed-loop IORG. Thus, we can simplify the Eq. (3) as follows
Id(f)=(k1+Δk1)sin(2πfτ)
where f=f0+Δfwith f0 the center frequency of resonator curve and is an integral multiple of 1/τ and Δf the frequency difference between real optical frequency and resonant frequency, and m = 2q/(1 + q2). Considering the closed-loop system operates at zero point, we can obtain cos2πΔfτ1. Then, k1 is the gain of demodulation gain while Δf=0, and also the forward chain of double closed-loop system. Further, the optimal square frequency bias fbias can be obtained by the Eq. (5).

 figure: Fig. 2

Fig. 2 The flow diagram of detected signals during the modulation and demodulation processes including (a) the interference intensity IT as a function of light frequency (b) the phase modulated carrier ϕ and equal frequency bias fbias; (c) the time-domain output Im of photoelectric detector after modulation; (d) the demodulated signal Id .

Download Full Size | PDF

fbias|k1fbias=0=arccos8m2+112m2πτ

Therefore, based on the optimal square frequency bias fbias, the optimal frequency of phase modulation triangular wave can be obtained. In our designed IORG system, a buried-type Ge-doped silica-based waveguide-type ring resonator with a diameter of 60mm and a finesse of 82 is employed by which we obtain q = 0.9624 and τ = 0.911ns. The integrated optical phase modulator (IOPM) made by LiNbO3 waveguide, has two independent modulation arms, a bandwidth 200MHz and a half-wave voltage Vπ 5.6V. The phase modulation triangular wave has peak voltage VPP 22.4V.

Based on the above analysis, the optimized frequency bias 3.8657MHz is obtained, and then the optimal frequency of phase modulation triangular wavefT=fbiasVπ/VPPis 0.966MHz. Therefore, the largest demodulation gain and the least nonlinearity can be realized as shown in Fig. 3(a). Meanwhile, as illustrated in Fig. 3(b), the multiple-beam interference effect of waveguide ring resonator brings in system nonlinearity, however the proposed modulation frequency guarantees the IORG system operates at largest work region and optimal detection sensitivity.

 figure: Fig. 3

Fig. 3 (a) the relationship between demodulation gain k1 and modulation frequency; (b)the demodulated signals Id shown as functions of frequency deviation for different modulation frequencies.

Download Full Size | PDF

In this part, the influencing mechanism of optical effect nonlinearity on detection sensitivity of the IORG system is investigated, and the optimized modulation frequency fbias is derived to obtain the maximum demodulation gain k1. Consequently, the closed-loop error sensitivity before controller can be improved with a precise ΔfΩ. However, the IORG is still a nonlinear system, and it is necessary to design a control algorithm that takes nonlinearity into consideration for ensuring the fast dynamic response of the IORG system.

3.2 Analysis of detection accuracy

The high dynamic model of double closed-loop IORG is established and analyzed in the following. In order to obtain ΔfΩ, the input of AVTL controller is the differential mode signal of PD1 and PD2. In ideal conditions, the frequency difference ΔfΩ can be accurately demodulated through the digital closed-loop processor in AVTL. However, the unavoidable frequency-locking accuracy Δfn existing in LFLL causes considerable difficulty in precisely extracting frequency difference ΔfΩ. We analyze the working principle of frequency-locking accuracy on AVTL. According to the frequency relationship, Δf=ΔfΩ+Δfn, we can obtain the demodulated frequency difference in AVTL as

Id2=(k1+Δk1)sin(2π(ΔfΩ+Δfn)τ)
which can be also rewritten as

Id2=(k1+Δk1)cos2πΔfnτsin2πΔfΩτ+(k1+Δk1)cos2πΔfΩτsin2πΔfnτ

Noting that sin2πΔfΩτ2πΔfΩτ, we find that the first term of demodulation result, k1cos2πΔfnτsin2πΔfΩτ, is the closed-loop error of AVTL. Due to the real-time locking of laser frequency, the closed-loop error in LFLL Δfn is within a small range, so that we have cos2πΔfnτ1, which obviously leads to the parameter uncertainty in AVTL. As for the frequency-locking noise signal g(x1(k))=sin2πΔfnτ existing in the second term of demodulation result k1cos2πΔfΩτsin2πΔfnτ, we know that its mean E(g(x1(k))) is 0 and variance of noise σ satisfying E(g2(x1(k)))=σ2 is finite. Therefore, the frequency-locking noise signal can be considered as the noise term with zero mean and finite variance. Meanwhile, the closed-loop mean error ΔfΩ in AVTL is 0, so that v2(k)=(k1+Δk1)cos2πΔfΩτ can be supposed to be the finite-energy disturbance signal which belonging to L2[0, + ∞).

If the error of laser frequency locking Δfn is zero, we can see that the disturbing signal g(x1(k))v2(k) will disappear in AVTL. However, the error of laser frequency locking does exist in a practical system, and influences the AVTL of IORGs. Thus, the closed-loop model of AVTL is a stochastic dynamical system introduced by the frequency locking noise.

The system gain parameters can be defined as k1 and k2, which satisfy k1 = k2. We define the gain variation of the forward channels of two loops as Δk1 and Δk2 which satisfyk2+Δk2=(k1+Δk1)cos2πΔfnτ. Meanwhile, the output light power is considerably affected by environmental temperature fluctuation, which also leads to the gain variation of forward channels of two loops. Based on the above analysis, we suppose that Δk1 and Δk2 are less than ± 10%. Thus, the synchronizing demodulated closed-loop signal of AVTL can be obtained respectively as

Id2(Δf)=(k2+Δk2)sin2πΔfΩτ+g(x1(k))v2(k)

4. Design of closed-loop detection

In this section, a robust control algorithm based on the closed-loop signal detection method discussed above is proposed to establish a fast tracking system with high detection accuracy.

Because the frequency locking accuracy can only be demodulated from PD1, the tracking result of the laser frequency is free from the influence of angular closed-loop error. However, the parameter uncertainty Δk1 of LFLL exists based on Eq. (4). According to controllability canonical form of automatic control theory, the dynamic equation of LFLL of IORG system can be deduced as

x1(k+1)=A1x1(k)+B1(k1+Δk1)sin(kf1KC1x1(k))+D1w1(k)
Meanwhile, the mathematic model for the AVTL of IORG system is given by
x2(k+1)=A2x2(k)+B2(k2+Δk2)sin(kf2KC2x2(k))+B2g(x1(k))v2(k)+D2w2(k)
whereA2=A1=[1100011000110001],B2=B1=[0001].w1(k)andw2(k)are the one-dimensional finite-energy disturbance inputs caused by optical error and noise, respectively.D1 and D2Rn are the disturbance intensity constant vectors. x1(k) and x2(k)Rn are the state variables, kf1 and kf2 are the gain of feedback chain of LFLL and AVTL, KC1 and KC2R1 × n are the feedback gain matrices of controller of LFLL and AVTL, respectively.

In order to take the system nonlinearity, the inevitable parameter drift, the external disturbance and the frequency-locking noise into consideration, Eq. (9) and (10) are combined to establish the stochastic disturbance model of double closed loop system given by

x(k+1)=Ax(k)+B¯f(kf1KC1x1(k),kf2KC2x2(k))+Ig(x1(k))v2(k)+Dw(k)
whereA=(A100A2),I¯=(0B2),D=(D1D2)and B¯=(B1(k1+Δk1)00B2(k2+Δk2)). The matrix ΔB=(B1Δk100B2Δk2) can be transformed into ΔB=HF(k)E due to the bounded parameter Δk1 and Δk2, where H and E are constant matrices of appropriate dimensions that describe the variation intensity of the gain, F(k) is an uncertain matrix of appropriate dimension that satisfies the condition F(k)TF(k)I.

A desired control algorithm should be designed as robust as possible in practical engineering to overcome the previously-mentioned disadvantages. Therefore, it is of great significance to introduce the performance index for enhancing the disturbance rejection attenuation level of control algorithm. In addition to the detection accuracy, the exponential stability with desired H performance is necessary to ensure IORG system with a good tracking and a fast dynamic response. Thus, we present the method for designing the gain matrices, KC1 and KC2 of controller in the following.

In order to derive the main results, we first give some definitions and lemmas.

Definition 1 [20]:The system (11) with w(k) = 0 is said to be exponentially stable in a mean-square sense, if there exist some scalars κ>0 and 0<α<1, such that the solution x(k) of system (11) satifies E{x(k)}καkk0x(k0) for all any x(k0)Rn and kk0.

Definition 2 [21]: For a given scalar γ>0, the system (11) is said to be exponentially stable in a mean-square sense with an expinential H performance γ, if it is exponentially stable and there holdss=0+E{eT(s)e(s)}s=0+γ2wT(s)w(s) for all non-zero w(k)l2[0,) under zero initial condition.

Lemma 1 (Schur Complement) [22]: For a given symmetric matrix S with the form S = [Sij],S11Rr×r,S12Rr×(nr), S22R(n-r) × (n-r), then, S<0 if and only if S11<0,S22S21S111S12<0 or S22<0,S11S12S221S21<0.

Lemma 2 [23]: Q, H and E are known constant matrixes with appropriate dimensions, where Q is a symmetrical matrix stisfying Q = QT, the following inequality Q + HFE + ETFTHT <0 holds, for all satisfying FTFI, if and only if there exists a scalar ε>0, such that Q+ε1HHT+εETE<0.

Then, we first consider the exponential stability of system (11) with w(k) = 0 by a Lyapunov-based approach. The following theorem provides a sufficient condition to guarantee that IORG is mean-square exponentially stable for obtaining the desirable dynamic performance.

Theorem 1: For a given scalar 0<α<1, the system (11) is locally exponential stable in the mean square sense with w(k) = 0, if there exists positive definite matrix PRn×n, feedback gain matrices Kc1,Kc2R1×n, and positive scalars ε1, ε2, ε3 such that

(Φ+ε3ETEϕ1TPϕ2TP0*P0PH**P0***ε3I)<0
where Φ=(αPε1kf12H1KC1ε2kf22H2KC2*ε1I0**ε2I), ϕ1=(AB), ϕ2=(σG0H10).

Proof: See Appendix A.

Next, based on the results obtained in Theorem 1, we now consider the H performance for system(11) with noise w(k) to improve the detection accuracy of IORG. The following theorem guarantees that IORG is mean-square exponentially stabile with the desired H performance.

Theorem 2: For a given scalar 0<α<1 and the H performance index γ, the system (11) is mean-square exponential stable with a prescribed H performance index γ, if there exist symmetric positive definite matrix P,Q,RRn×n, feedback gain matrices Kc1,Kc2R1×n, and positive scalars ε1, ε2, ε3, such that

(Φ˜+ε3ETEϕ^1TPϕ˜2P0P0PHP0ε3I)<0
where

Φ˜=(αP+Iε1kf12H1KC1ε2kf22H2KC20ε1I00ε2I0γ2I),ϕ^1=(ABD),ϕ˜2=(σG0H100).

Proof: See Appendix B.

Theorem 2 provides the design rule of closed-loop feedback gain matrices of IORGs considering the nonlinearity, optical parameter uncertainty, external noise and frequency-locking disturbance of the LFLL. It guarantees that the system has mean-square exponential stability and optimal detection sensitivity with a prescribed H performance index γ. It is noted that, the perturbation parameters and frequency-locking disturbance of the LFLL in closed-loop control system of IORG vary dynamically due to the uncertain external environment, which cause the angular velocity tracking loop of IORG to be a parameter-uncertainty existing and noise-perturbed stochastic dynamical system.

We herewith obtain the model of double closed-loop detection system as shown in Fig. 4. When the control parameters of two loops are designed, the closed-loop error of AVTL does not affect the LFLL accuracy. Therefore, the control parameter KC1 of the LFLL, whose design method of control parameter in LFLL is simple, can be designed first, and its closed-loop control model is similar with the model in [24]. Based on the obtained KC1 referred to Theorems 1 in [24], we can calculate the control parameter KC2 based on Theorems 2 with the designed KC1.

 figure: Fig. 4

Fig. 4 The block diagram of the double closed-loop detection scheme of IORG.

Download Full Size | PDF

The entire forward-path circuit consists of interference, photoelectric detector, analog amplifier, A/D converter and demodulation. The analog amplifier, drive circuit and linear amplifier are analog. The modulation, demodulation, controller and sawtooth generating are digital. During the demodulation, after the sampling and accumulating in the positive and negative half-period of the modulation triangular wave, the closed-loop error of each loop can be demodulated by subtraction of the positive and negative accumulation. And, Section 3.1 and 3.2 deduce that the forward-path circuit can be supposed as a proportional link multiplied by a sine function.

The circuit gain of the forward channel Kqd=KGaKadNdem where KGa is the gain of pre-amplifier, Kad is the gain of A/D converter and Ndem is the sampling number of half demodulation period.

In accordance with the composition of forward channel of closed-loop system as shown in Fig. 4, the demodulation gain is given by

k1=21+NadmI0ηRβKGaNdem(1+q2)Vref
where Kad=2Nad/Vrefad and Nad and Vrefad are the resolution and reference voltage of A/D converter respectively. The values of the key parameters are shown in Table 1. Given that n=2,A1=(1101),B1=(01), D1=(01), P=I, k1 = 0.0265 and kf1=kGdVrefdakla/2Nda=3.0469×104 where kGd is the gain of drive circuit, and kla is the gain of linear amplifier. Then, the matrix parameter of KC1 is calculated asKC1=(0.31250.00488) by linear matrix inequality (LMI) [25]. The feedback gain of AVTLkf2=kGdfclkVrefVpp2Nda/4VπΔHsaw where ΔHsaw is the digital amplitude of sawtooth wave, and fclk is the working clock of FPGA. In our designed IORG system, k2 = k1 = 0.0265 and kf2 = 0.0571, Then given that α = 0.05, Q = I, the feedback gain matrix of AVTL can be calculated as KC2=(0.20.00078) based on Theorem 2 with the designed KC1.

Tables Icon

Table 1. The value of parameters of block diagram

The realization of the control algorithm illustrated in Fig. 4 provides an efficient and robust method for optimizing the detection precision and fast tracking performance of IORG. Theorem 2 provides the method for designing the AVTL controller, which is robust for suppressing the influence of locking frequency accuracy and the main closed-loop disturbance.

5. Experiments and results

Some experiments are conducted in this section to verify the validity and effectiveness of the designed signal processing scheme of IORG. The optical configuration of the IORG experimental setup as previously presented in Fig. 1 is composed of the laser, the optical isolator (ISO), the IOPM and waveguide-type ring resonator. A piezoelectric low phase-noise fiber laser operating around 1550nm is utilized with linewidth less than 1kHz and outputting power 36.2mW. The type and parameter of IOPM and waveguide-type ring resonator is shown in Section 3.1.

The gyro output is at millisecond level and the closed-loop step of tracking angular velocity is at dozens of nano-second level, which brings difficulty in verifying the influence mechanism analysis of laser locking noise on main closed-loop performance and the validity of the double closed-loop controller. In this work, a high speed sampling design platform based on Fig. 4 is developed. The digital signal processing tasks, performed by the reconfigurable FPGA of Kintex-7 XC7K410T with working clock 100MHz, includes the modulation signal generation, the demodulation and control algorithm process. The FPGA is connected to host computer by PXI-e bus. Therefore, the designed platform not only monitors the nano-second state variables and analyze the frequency spectrum of the key signal, but also assist us in verifying the influence mechanism analysis of laser locking noise on the main closed-loop performance.

The step response experiment is first conducted to demonstrate that the dynamic performance of IORG has fast exponential stability in Fig. 5. A ramp signal with digital stair height hstair 93.75 is generated by FPGA with working clock fclk 100MHz, then converted into digital amplitude 65536 with 16 bits D/A, and applied to phase modulator along with the feedback signal. The resonant frequency difference of ramp signal ΔfΩ=hstairfclk/2nD/A is obtained as 143.05kHz and the corresponding slope ratio is equivalent to 300°/s of the system input. The time domain output curves of CW, CCW loop and AVTL controller are observed by an oscilloscope in Fig. 5(a). It can be observed that the resulting output of digital controller is a standard sawtooth signal, which is imposed on the IOPM to feedback the resonant frequency difference ΔfΩ. The output light intensity of resonator is further converted into voltage signal by photoelectric detector, which is depicted in Fig. 5(b).

 figure: Fig. 5

Fig. 5 Implementation of the closed-loop IORG system. (a) The oscilloscope display of output signals under 300°/s rotational speed measurement. Curve A is the CW output which detected by PD2 and displayed on oscilloscope, Curve B is the CCW output which detected by PD1, and curve C is AVTL controller output; (b) the PD1 curve of LFLL; (c) the frequency locking noise of LFLL after the frequency of laser is locked to the center frequency of curve resonator (d) the step response of IORG system.

Download Full Size | PDF

The control update time of frequency-locking loop is 0.2ms, and the experimental data depicted in Fig. 5(c) is the demodulated signal of LFLL. Since the demodulated signal of LFLL is mainly the laser's locking noise, so that we may estimate the maximum variance of the laser's locking noise by the experimental data shown in Fig. 5(c). The estimated variance of locking-frequency noise σ is 44.56Hz illustrated in Fig. 5(c), which is used to optimize the closed-loop control algorithm of AVTL.

In consideration of the influence of laser locking-frequency noise, the simulation and experiment results are presented in Fig. 5 about the dynamic performance of IORG by using the proposed control algorithm. And, the simulation result of the time-domain model shows a fast rise time of 73.265us. The experimental result demonstrates that the response time of IORG is 76us, which can be sampled by FPGA and illustrated in Fig. 5(d). The exponential stability coefficient α stays consistency with the theoretical value in Section 4, and the experimental result 76us is roughly consistent with the theoretical result 73.265us. It clearly shows that the improved dynamic performance of IORG with the influence of laser locking-frequency noise can meet the exponential stability using the proposed control algorithm.

Further, the measuring accuracy experiment is conducted to validate the effectiveness of our proposed closed-loop algorithm with a prescribed H performance. As shown in Fig. 6, a bias stability of 7.04°/h is demonstrated with an integration time of 10s over one-hour test, and the Allan deviation analysis of typical one-hour test shows that the bias stability has reached 1.841°/h. In contrast, without adopting our closed-loop detection method, the bias stability of the IORG system with an integration time of 10s is 46.8°/h adopting the same optical system [12]. Meanwhile, the IORG has the best long-term bias stability that ever been reported for an IORG based on a silica waveguide ring resonator. The experiment results verify that the proposed algorithm can effectively suppress the locking frequency noise, the system nonlinearity and the external disturbance, and therefore improve the detection accuracy and dynamic performance of the IORG system.

 figure: Fig. 6

Fig. 6 The results of measuring accuracy experiment of the IORG system. a) 1h test of bias stability; c) Allan variance of the 1h test data.

Download Full Size | PDF

Finally, we conduct the measurement experiment of the scale factor to verify the linearity of IORG by using the proposed closed-loop controller, which is illustrated in Fig. 7. The obtained scale factor (SF) of IORG prototype is 7.6 (°/s)−1 and the nonlinearity of scale factor is 344.71ppm. These results demonstrate that the proposed signal processing method of IORG successfully corrects the system nonlinearity.

 figure: Fig. 7

Fig. 7 Results of scale factor experiment of the IORG system. a) the test of scale factor; b) The out curve and corresponding relative fitting error of the IORG system.

Download Full Size | PDF

6. Conclusions

In a practical IORG system, the performance of main angular velocity tracking loop is greatly affected by the locking-frequency accuracy of the laser frequency-locking loop. Meanwhile, the detection accuracy of IORG is degenerated by several adverse factors such as the nonlinear optical effect, the optical parameter fluctuation, and the unavoidable external noise existing in practical environment. In this paper, a novel closed-loop signal detection scheme is designed to extract the angular velocity of IORG system with high sensitivity and high precision of measurement. Firstly, we derive the optimal modulation frequency to obtain the maximum linear work region and demodulation gain for suppressing the nonlinearity of optical effect. Secondly, we explored the influence mechanism of the above four factors on system sensitivity and measurement accuracy. The mathematical model of the main angular rate tracking loop is also established. Thirdly, a robust control algorithm of the IORG for suppressing the frequency-locking disturbance and the existing optical noise is then proposed to ensure mean-square exponential stability with a prescribed H performance. The experimental results demonstrate the effectiveness and usefulness of our proposed method, by which the detection accuracy and dynamic performance of IORG have been greatly improved. The result of our research is of significance in promoting the wide application of IORGs in inertial navigation.

Appendix A

Proof of Theorem 1:

We select the following Lyapunov function for system (11) as V(k)=xT(k)Px(k), where P is defined as that in Theorem 1. It is noted that sin(x)is a monotone increasing and differentiable function in the rangeπxπ, and sin() locally satisfies the Lipschitz condition0<sin(x(k))/x(k)1 and g(0)0. It is easy to see that there exists positive scalar 0<G01 such thatsin(x1(k))G0x1(k). And, we can obtain the following condition

fT(KCixi(k))(f(KCixi(k))KCixi(k))0,wherei=1,2.
Thus, we have
E{V(k+1)αV(k)}E{xT(k+1)Px(k+1)αxT(k)Px(k)ε1sinT(KC1x1(k))(sin(KC1x1(k))KC1x1(k))ε2sinT(KC2x2(k))(sin(KC2x2(k))KC2x2(k))}
In consideration ofsin(x1(k))G0x1(k) and E(vT(k)v(k))=σ2, the inequality holds:
E{(I¯g(x1(k))v(k))TP(I¯g(x1(k))v(k))}E{xT(k)(σG0H1)TP(σG0H1)x(k)}.
Noting that E(v(k))=0 and the dynamics of System (11), we have
E{V(k+1)αV(k)}=E{ς1T(k)(ϕ¯1TPϕ¯1)ς1(k)+ς1T(k)Φς1(k)+xT(k)(σG0H1)TP(σG0H1)x(k)}
where ς1T(k)=(xT(k)sinT(KC1x1(k))sinT(KC2x2(k))),x1(k)=H1x(k)=(I0)x(k),x2(k)=H2x(k)=(0I)x(k), ϕ¯1=(AB¯) and Φ=(αPε12H1KC1ε22H2KC2*ε1I0**ε2I).

Thus, we obtain

E{V(k+1)αV(k)}E{ς1T(k)(ϕ¯1TPϕ¯1+ϕ2TPϕ2+Φ)ς1(k)}

We can see that the sufficient condition for E{V(k+1)αV(k)}<0is if and only if ϕ¯1TPϕ¯1++ϕ2TPϕ2+Φ<0.

Then by using Lemma 1, ϕ¯1TPϕ¯1++ϕ2TPϕ2+Φ<0 holds if and only if

(Φϕ¯1TPϕ2TPPϕ¯1P0Pϕ20P)<0

Because the matrixΔB is bounded, from Lemma 2, one has

(Φϕ¯1TPϕ¯2TPPϕ¯1P0Pϕ¯20P)=(Φϕ1TPϕ2TPPϕ1TP0Pϕ20P)+MTF(k)N+NTFT(k)M<0
where ϕ1=(AB),B=(B¯1k100B¯2k2),M=(0HTP0),N=(E00).

By Lemma 1, Eq. (18) is equivalent to

(Φ+ε3ETEϕ1TPϕ2TP0*P0PH**P0***ε3I)<0

Then, we can see thatE{v(k)}(α)kk0V(k0) such that E{x(k)}β2β1(α12)kk0x(k0), where β1=minλ(P) and β2=maxλ(P) are the minimum and maximum eigenvalues of matrix P, respectivly. According to Definition 1, we conclude that the system(12) is exponentially stable in the mean-square sense with w(k)=0.

Appendix B

Proof of Theorem 2

We consider the performance Γ(k)=xT(k)x(k)γ2wT(k)w(k) for any nonzero w(k)L2[0,). Following the similar proof of Theorem 1 with the same Lyapunov-Krasovskii function, we have

E{V(k+1)αV(k)+Γ(k)}E{V(k+1)αV(k)+Γ(k)ε1fT(KC1x1(k))(sin(KC1x1(k))KC1x1(k))ε2fT(KC2x2(k))(sin(KC2x2(k))KC2x2(k))}=E{ξ2T(k)(ϕ˜1TPϕ˜1++ϕ˜2TPϕ˜2+Φ˜)ξ2(k)}
where

ξ2(k)=(x(k)Tx(k1)Tsin(Δφ(k))Tsin(Δφ(k1))Tw(k)T)T

According to the above analysis, we can see that a sufficient condition for E{V(k+1)αV(k)}<0 is equivalent to ϕ˜1TPϕ˜1++ϕ˜2TPϕ˜2+Φ˜<0,where ϕ˜1=(AB¯D),ϕ˜2=(σG0H100),Φ˜=(αP+Iε12H1KC1ε22H2KC20ε1I00ε2I0γ2I)

It follows from Lemma 1, ϕ˜1TPϕ˜1++ϕ˜2TPϕ˜2+Φ˜<0 holds if and only if

(Φ˜ϕ˜1TPϕ˜2PP0P)<0

By using Lemma 2, Eq. (20) can be rewritten as

(Φ˜ϕ^1TPϕ˜2PP0P)+M2TF2(k)N2+N2TF2(k)M2(Φ˜+ε3ETEϕ^1TPϕ˜2P0P0PHP0ε3I)<0

where M2=(0HTP0) and N2=(E00).

Finally, we obtain

E{V(k+1)}<E{αV(k)Γ(k)}<E{αkk0V(k0)s=k0kαksΓ(s)}

On the other hand, considering the definition of the Lyapunov function V(k+1)>0and the zero initial condition x(k0)=0, we can see thats=k0kαksE(xT(s)x(s))<s=k0kαksE(γ2wT(s)w(s)). By suming both sides of the equation from k=k0 to k=, one gets s=k0+E{xT(s)x(s)}<γ2s=k0+wT(s)w(s). Then, from Theorem 2, it follows that the closed-loop system (11) can satisfy H performance. This completes the proof.

Funding

National Natural Science Foundation of China (No. 61405005).

References and links

1. C. Ciminelli, F. Dell’Olio, C. E. Campanella, and M. N. Armenise, “Photonic technologies for angular velocity sensing,” Adv. Opt. Photonics 2(3), 370–404 (2010). [CrossRef]  

2. J. Haavisto and G. A. Pajer, “Resonance effects in low-loss ring waveguides,” Opt. Lett. 5(12), 510–512 (1980). [CrossRef]   [PubMed]  

3. M. Á. Guillén-Torres, E. Cretu, N. A. F. Jaeger, and L. Chrostowski, “Ring resonator optical gyroscopes-Parameter optimization and robustness analysis,” J. Lightwave Technol. 30(12), 1802–1817 (2012). [CrossRef]  

4. N. Barbour, “Inertial Navigation Sensors [R],” Charles Stark Draper Lab Inc Cambridge ME, 2010.

5. N. Barbour and G. Schmidt, “Inertial sensor technology trends,” IEEE Sensors 1(4), 332–339 (2001). [CrossRef]  

6. W. Liang, V. S. Ilchenko, A. A. Savchenkov, E. Dale, D. Eliyahu, A. B. Matsko, and L. Maleki, “Resonant micro photonic gyroscope,” Optica 4(1), 114–117 (2017). [CrossRef]  

7. J. Wang, L. Feng, Y. Tang, and Y. Zhi, “Resonator integrated optic gyro employing trapezoidal phase modulation technique,” Opt. Lett. 40(2), 155–158 (2015). [CrossRef]   [PubMed]  

8. X. Wang, Z. He, and K. Hotate, “Reduction of polarization-fluctuation induced drift in resonator fiber optic gyro by a resonator with twin 90 ° polarization-axis rotated splices,” Opt. Express 18(2), 1677–1683 (2010). [CrossRef]   [PubMed]  

9. J. I. Thorpe, K. Numata, and J. Livas, “Laser frequency stabilization and control through offset sideband locking to optical cavities,” Opt. Express 16(20), 15980–15990 (2008). [CrossRef]   [PubMed]  

10. K. Hotate and M. Harumoto, “Resonator fiber optic gyro using digital serrodyne modulation,” J. Lightwave Technol. 15(3), 466–473 (1997). [CrossRef]  

11. M. Lei, L. Feng, Y. Zhi, and H. Liu, “Effect of intensity variation of laser in resonator integrated optic gyro,” Appl. Opt. 52(19), 4576–4581 (2013). [CrossRef]   [PubMed]  

12. L. Hong, C. Zhang, L. Feng, M. Lei, and H. Yu, “Effect of phase modulation nonlinearity in resonator micro-optic gyro,” Opt. Eng. 50(9), 094404 (2011). [CrossRef]  

13. J. Wang, L. Feng, Q. Wang, H. Jiao, and X. Wang, “Suppression of backreflection error in resonator integrated optic gyro by the phase difference traversal method,” Opt. Lett. 41(7), 1586–1589 (2016). [CrossRef]   [PubMed]  

14. Y. Zhi, L. Feng, M. Lei, and K. Wang, “Low-delay, high-bandwidth frequency-locking loop of resonator integrated optic gyro with triangular phase modulation,” Appl. Opt. 52(33), 8024–8031 (2013). [CrossRef]   [PubMed]  

15. J. Wang, L. Feng, Y. Zhi, H. Liu, W. Wang, and M. Lei, “Reduction of backreflection noise in resonator micro-optic gyro by integer period sampling,” Appl. Opt. 52(32), 7712–7717 (2013). [CrossRef]   [PubMed]  

16. J. Wu, M. Smiciklas, L K. Strandjord, T. Qiu, W. Ho and G. A. Sanders, “Resonator fiber optic gyro with high backscatter-error suppression using two independent phase-locked lasers”, International Society for Optics and Photonics 2015: 96341O–96341O–4.

17. H. Ma, J. Zhang, L. Wang, and Z. Jin, “Double closed-loop resonant micro optic gyro using hybrid digital phase modulation,” Opt. Express 23(12), 15088–15097 (2015). [CrossRef]   [PubMed]  

18. H. C. Lefevre, “The fiber-optic gyroscope,”. In French, Artech house, 1993, pp. 20–23.

19. L. Feng, J. Wang, Y. Zhi, Y. Tang, Q. Wang, H. Li, and W. Wang, “Transmissive resonator optic gyro based on silica waveguide ring resonator,” Opt. Express 22(22), 27565–27575 (2014). [CrossRef]   [PubMed]  

20. D. Zhang, L. Yu, Q. G. Wang, C. J. Ong, and Z. G. Wu, “Exponential H filtering for discrete-time singular switched system with time-varying delays,” J. Franklin Inst. 349(7), 2323 (2012). [CrossRef]  

21. T. Liu, J. Zhao, and J. H. David, “Exponential synchronization of complex delayed dynamical networks with switching topology,” IEEE Trans. Circuits Syst. I, Reg. Papers 57(11), 2967–2980 (2010).

22. S. Boyd, L. E. Ghaoui, E. Feron, and V. Balakrishinan, “Linear Matrix Inequalities in Systems and Control Theory,” Philadelphia, PA, USA: SIAM, (1994).

23. B. Shen, Z. Wang, and H. Shu, “Distributed H1-consensus Filtering in Sensor Networks with Multiple Missing Measurements: The Finite-horizon Case,” Automatica 64(46), 1682–1688 (2010). [CrossRef]  

24. W. Deng, H. Li, C. Zhang, and P. Wang, “Optimization of Detection Accuracy of Closed-Loop Optical Voltage Sensors Based on Pockels Effect,” Sensors (Basel) 17(8), 1723 (2017). [CrossRef]   [PubMed]  

25. S. Boyd, L. E. Ghaoui, E. Feron, and V. Balakrishinan, Linear Matrix Inequalities in Systems and Control Theory. Philadelphia, PA, USA: SIAM, 1994.

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 The principle scheme of the double closed-loop IORG based on Sagnac effect. The triangular wave for phase modulation is applied on the upper arm and bottom arm of the IOPM. The feedback sawtooth wave is applied on the bottom arm of the IOPM. Meanwhile, the phase modulation triangular wave and digital feedback sawtooth wave are all differentially applied on the arms of the IOPM. The IOPM achieves the phase modulation and also enables the closed-loop control of AVTL.
Fig. 2
Fig. 2 The flow diagram of detected signals during the modulation and demodulation processes including (a) the interference intensity IT as a function of light frequency (b) the phase modulated carrier ϕ and equal frequency bias fbias; (c) the time-domain output Im of photoelectric detector after modulation; (d) the demodulated signal Id .
Fig. 3
Fig. 3 (a) the relationship between demodulation gain k1 and modulation frequency; (b)the demodulated signals Id shown as functions of frequency deviation for different modulation frequencies.
Fig. 4
Fig. 4 The block diagram of the double closed-loop detection scheme of IORG.
Fig. 5
Fig. 5 Implementation of the closed-loop IORG system. (a) The oscilloscope display of output signals under 300°/s rotational speed measurement. Curve A is the CW output which detected by PD2 and displayed on oscilloscope, Curve B is the CCW output which detected by PD1, and curve C is AVTL controller output; (b) the PD1 curve of LFLL; (c) the frequency locking noise of LFLL after the frequency of laser is locked to the center frequency of curve resonator (d) the step response of IORG system.
Fig. 6
Fig. 6 The results of measuring accuracy experiment of the IORG system. a) 1h test of bias stability; c) Allan variance of the 1h test data.
Fig. 7
Fig. 7 Results of scale factor experiment of the IORG system. a) the test of scale factor; b) The out curve and corresponding relative fitting error of the IORG system.

Tables (1)

Tables Icon

Table 1 The value of parameters of block diagram

Equations (28)

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

Δ f Ω = D n e λ Ω
I T = I 0 ηβR 1+ q 2 2qcos2πfτ
I d ( f )= I 0 ηβR K qd 1+ q 2 2qcos2π( f+ f bias )τ I 0 ηβR K qd 1+ q 2 2qcos2π(f f bias )τ = 2msin2π f bias τ 1+ m 2 2 cos4π f bias τ+ m 2 2 cos4πfτ2mcos2πfτcos2π f bias τ I 0 ηβR K qd 1+ q 2 sin( 2πfτ )
I d ( f )=( k 1 +Δ k 1 )sin( 2πfτ )
f bias | k 1 f bias =0 = arccos 8 m 2 +1 1 2m 2πτ
I d 2 =( k 1 +Δ k 1 )sin( 2π( Δ f Ω +Δ f n )τ )
I d 2 =( k 1 +Δ k 1 )cos2πΔ f n τsin2πΔ f Ω τ+( k 1 +Δ k 1 )cos2πΔ f Ω τsin2πΔ f n τ
I d 2 ( Δf )=( k 2 +Δ k 2 )sin2πΔ f Ω τ+g( x 1 ( k ) ) v 2 ( k )
x 1 ( k+1 )= A 1 x 1 ( k )+ B 1 ( k 1 +Δ k 1 )sin( k f 1 K C 1 x 1 ( k ) )+ D 1 w 1 ( k )
x 2 ( k+1 )= A 2 x 2 ( k )+ B 2 ( k 2 +Δ k 2 )sin( k f 2 K C 2 x 2 ( k ) )+ B 2 g( x 1 ( k ) ) v 2 ( k )+ D 2 w 2 ( k )
x(k+1)=Ax(k)+ B ¯ f( k f 1 K C 1 x 1 (k), k f 2 K C 2 x 2 (k))+Ig( x 1 (k)) v 2 (k)+Dw(k)
( Φ+ ε 3 E T E ϕ 1 T P ϕ 2 T P 0 * P 0 PH * * P 0 * * * ε 3 I )<0
( Φ ˜ + ε 3 E T E ϕ ^ 1 T P ϕ ˜ 2 P 0 P 0 PH P 0 ε 3 I )<0
Φ ˜ =( αP+I ε 1 k f 1 2 H 1 K C 1 ε 2 k f 2 2 H 2 K C 2 0 ε 1 I 0 0 ε 2 I 0 γ 2 I ), ϕ ^ 1 =( A B D ), ϕ ˜ 2 =( σ G 0 H 1 0 0 ).
k 1 = 2 1+ N ad m I 0 ηRβ K Ga N dem ( 1+ q 2 ) V ref
f T ( K C i x i ( k ) )( f( K C i x i ( k ) ) K C i x i ( k ) )0, where i=1, 2.
E{ V( k+1 )αV( k ) }E{ x T ( k+1 )Px( k+1 )α x T ( k )Px( k ) ε 1 sin T ( K C 1 x 1 ( k ) )( sin( K C 1 x 1 ( k ) ) K C 1 x 1 ( k ) ) ε 2 sin T ( K C 2 x 2 ( k ) )( sin( K C 2 x 2 ( k ) ) K C 2 x 2 ( k ) ) }
E{ ( I ¯ g( x 1 ( k ) )v(k) ) T P( I ¯ g( x 1 ( k ) )v(k) ) }E{ x T ( k ) ( σ G 0 H 1 ) T P( σ G 0 H 1 )x( k ) }
E{ V( k+1 )αV( k ) }=E{ ς 1 T ( k ) ( ϕ ¯ 1 T P ϕ ¯ 1 ) ς 1 ( k )+ ς 1 T ( k )Φ ς 1 ( k ) + x T ( k ) ( σ G 0 H 1 ) T P( σ G 0 H 1 )x( k ) }
E{ V( k+1 )αV( k ) }E{ ς 1 T ( k )( ϕ ¯ 1 T P ϕ ¯ 1 + ϕ 2 T P ϕ 2 +Φ ) ς 1 ( k ) }
( Φ ϕ ¯ 1 T P ϕ 2 T P P ϕ ¯ 1 P 0 P ϕ 2 0 P )<0
( Φ ϕ ¯ 1 T P ϕ ¯ 2 T P P ϕ ¯ 1 P 0 P ϕ ¯ 2 0 P )=( Φ ϕ 1 T P ϕ 2 T P P ϕ 1 T P 0 P ϕ 2 0 P )+ M T F( k )N+ N T F T ( k )M<0
( Φ+ ε 3 E T E ϕ 1 T P ϕ 2 T P 0 * P 0 PH * * P 0 * * * ε 3 I )<0
E{ V( k+1 )αV( k )+Γ( k ) } E{ V( k+1 )αV( k )+Γ( k ) ε 1 f T ( K C 1 x 1 ( k ) )( sin( K C 1 x 1 ( k ) ) K C 1 x 1 ( k ) ) ε 2 f T ( K C 2 x 2 ( k ) )( sin( K C 2 x 2 ( k ) ) K C 2 x 2 ( k ) ) } =E{ ξ 2 T ( k )( ϕ ˜ 1 T P ϕ ˜ 1 ++ ϕ ˜ 2 T P ϕ ˜ 2 + Φ ˜ ) ξ 2 ( k ) }
ξ 2 (k)= ( x (k) T x (k1) T sin (Δφ(k)) T sin (Δφ(k1)) T w (k) T ) T
( Φ ˜ ϕ ˜ 1 T P ϕ ˜ 2 P P 0 P )<0
( Φ ˜ ϕ ^ 1 T P ϕ ˜ 2 P P 0 P )+ M 2 T F 2 ( k ) N 2 + N 2 T F 2 ( k ) M 2 ( Φ ˜ + ε 3 E T E ϕ ^ 1 T P ϕ ˜ 2 P 0 P 0 PH P 0 ε 3 I )<0
E{V(k+1)}<E{αV(k)Γ(k)}<E{ α k k 0 V( k 0 ) s= k 0 k α ks Γ(s) }
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.