CN101718557B - United signal processing method of double-polarization interference type optical fiber gyroscope - Google Patents

United signal processing method of double-polarization interference type optical fiber gyroscope Download PDF

Info

Publication number
CN101718557B
CN101718557B CN2009102374215A CN200910237421A CN101718557B CN 101718557 B CN101718557 B CN 101718557B CN 2009102374215 A CN2009102374215 A CN 2009102374215A CN 200910237421 A CN200910237421 A CN 200910237421A CN 101718557 B CN101718557 B CN 101718557B
Authority
CN
China
Prior art keywords
polarization
optical fiber
light
double
polarized light
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN2009102374215A
Other languages
Chinese (zh)
Other versions
CN101718557A (en
Inventor
王玉杰
王子南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University
Original Assignee
Peking University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Peking University filed Critical Peking University
Priority to CN2009102374215A priority Critical patent/CN101718557B/en
Publication of CN101718557A publication Critical patent/CN101718557A/en
Application granted granted Critical
Publication of CN101718557B publication Critical patent/CN101718557B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a united signal processing method of a double-polarization interference type optical fiber gyroscope, which belongs to the technical field of optical fiber communication. The method comprises the following steps of: (1) modulating the double-polarization interference type optical fiber gyroscope and outputting signal light at intervals by two linearly-polarized light inputlight paths; (2) measuring light power output on the two polarized light input light paths at different temperatures by a phase modulator; (3) according to a recorded temperature value and corresponding light power, determining the error ratio values rho (T) of two paths of equivalent rotation rates caused by temperature transition; (4) modulating the double-polarization interference type optical fiber gyroscope and outputting the signal light simultaneously by the two linearly-polarized light input light paths; (5) measuring the light power output on the two linearly-polarized light input light paths; (6) calculating rotation rates omega 1 and omega 2 respectively according to the light power measured in the step (5); and (7) calculating a modified rotation speed omega according to a formula omega=(omega1-rho(T)omega2)/(1-rho(T)). The invention carries out united analysis on noise in the double-polarization optical fiber gyroscope to obtain an accurate rotation angle speed.

Description

A kind of united signal processing method of double-polarization interference type optical fiber gyroscope
Technical field
The present invention relates to a kind of signal processing method of double polarizing light fiber gyroscope, relate in particular to a kind of united signal processing method of double-polarization interference type optical fiber gyroscope, belong to technical field of optical fiber communication.
Background technology
Optical fibre gyro is a kind of device that can the detection system space rotates rotation, is a kind of simple in structure, cost is lower, precision is high novel all solid state inertia device.The optical fibre gyro application is extensive, comprising: tactical weapon guidance, orientation, tracking and navigation in the Aero-Space, optics compass, automobile navigation instrument, dual-use fields such as robot control.
The ultimate principle of optical fibre gyro is based on Sagnac (sagnac) effect and the ring interferometer made.In the closed optical loop of a random geometry, send in opposite direction the two-beam of transmission from any point, if closed light path relative inertness space rotate along a certain direction, when two-beam returns starting point, relative variation can take place in phase place, and this phenomenon is called the Sagnac effect.Therefore, can be by detecting the interference light of two-way, the counter rotational angular velocity of releasing closed light path system.
The relative phase variation of the two-way light that rotation causes is called the Sagnac phase shift,
φ s = 2 π LD λ ‾ c Ω - - - ( 1 )
Wherein, L is a fiber optic coils length, and D is the fiber optic coils diameter, and λ is the centre wavelength of flashlight, and c is the light velocity in the vacuum, and Ω is system's rotational angular velocity.By formula (1) as can be seen, Sagnac effect and phase shift φ thereof sSize is irrelevant with the refractive index of medium.
The structure of Sagnac fibre optic interferometer needs satisfying reciprocity, makes two beam reversals transmission light wave have identical transport property.Reciprocity can play the effect of " common mode inhibition ", eliminates various ghost effects, thereby can measure delicately because the nonreciprocal Sagnac phase shift that rotation causes.
Particularly, reciprocity comprises the reciprocity of fiber coupler, single mode reciprocity and polarization reciprocity.
In traditional guarantor's bias tyre optical fibre gyro scheme,, only taked the single shaft mode of operation of fast axle or slow axis as the minimum reciprocal structure among Fig. 1.Just because of φ sWith the irrelevant character of the refractive index of dielectric material, adopt the dual-polarization optical fiber Gyroscope Design scheme that is operated in two axles, by the two-way output signal is carried out check and analysis, can further improve the precision and the stability of system.
Summary of the invention
The object of the present invention is to provide a kind of united signal processing method of double-polarization interference type optical fiber gyroscope, it has obtained angular velocity of rotation accurately by the noise in the dual-polarization optical fiber gyro is carried out Conjoint Analysis.
Technical scheme of the present invention is:
A kind of united signal processing method of double-polarization interference type optical fiber gyroscope, wherein said double-polarization interference type optical fiber gyroscope is: two linearly polarized light input light paths are protected shift/beam splitter through one and are connected with a coupling mechanism, described coupling mechanism opposite side is connected with a fiber optic loop, be connected a phase-modulator between described fiber optic loop and the described coupling mechanism, the steps include:
1) modulation double-polarization interference type optical fiber gyroscope makes the two described linearly polarized light input light paths output signal light of property at interval;
2) measure described phase-modulator when different temperatures, the luminous power of exporting on the two described polarized light input light paths;
3), determine because the two-way equivalence speed of rotation error ratio that temperature transient causes according to the temperature value and the corresponding luminous power of record ρ ( T ) = ( n 1 dn 1 dT ) / ( n 2 dn 2 dT ) ; N wherein 1Be fast axle refractive index, n 2Be the slow axis refractive index, T is a temperature;
4) modulation double-polarization interference type optical fiber gyroscope makes two described linearly polarized light input light paths output signal light simultaneously;
5) measure the luminous power of exporting on the two described linearly polarized light input light paths;
6) calculate speed of rotation Ω respectively according to step 5) institute photometry power 1And Ω 2
7) according to formula Ω = Ω 1 - ρ ( T ) Ω 2 1 - ρ ( T ) Calculate revised rotational speed Ω.
Further, described coupling mechanism is a polarization-maintaining coupler, and described fiber optic loop is a polarization-maintaining fiber coil.
Further, described linearly polarized light input light path comprises a signal modulator, and the linear polarization light source is connected with described guarantor's shift/beam splitter through described signal modulator.
Further, described linearly polarized light input light path comprises a circulator, and described signal modulator is connected with described guarantor's shift/beam splitter through described circulator, and described circulator is connected with a detector.
Further, described signal modulator is a phase-modulator.
Further, described signal modulator is an amplitude modulator.
Further, described ratio ρ ( T ) = ( n 1 dn 1 dT ) / ( n 2 dn 2 dT ) Computing method be: at first the expression formula of interference signal luminous power is launched according to Bessel's function; Then according to step 2) luminous power surveyed, utilize that wave filter takes out wherein one, two, the amplitude of four-time harmonic; At last according to extracted one, two, φ when the amplitude of four-time harmonic solves different temperatures Im, obtain the n under the different temperatures 1/ n 2Value, thus extrapolate dn iThe curve of/dT is determined described ratio ρ (T); Wherein i=1 or 2 corresponds respectively to fast axle or slow axis.
Further, adopt kalman filter method to handle the signal of two described linearly polarized light input light paths respectively, eliminate temperature and float, obtain the speed of rotation Ω on two branch roads 1' and Ω 2'; Then according to formula Ω '=(η Ω 2'-Ω 1')/(η-1) calculate revised speed of rotation Ω ', η=ξ n wherein 2/ n 1, ξ = | E 1 cw | 2 - | E 1 ccw | 2 | E 2 cw | 2 - | E 2 ccw | 2 , | E i Cw| 2Be clockwise light intensity, | E i Ccw| 2Be counterclockwise light intensity, i=1 or 2 corresponds respectively to fast axle or slow axis.
Further, utilize formula Ω s=(Ω 2Δ n 11Δ n 2)/(Δ n 1-Δ n 2) calculate revised speed of rotation Ω s, Δ n wherein 1For causing that by stress or vibration fast axle local indexes of refraction changes Δ n 2For causing that by stress or vibration the slow axis local indexes of refraction changes.
Fig. 2 is a dual-polarization optical fiber Gyroscope Design scheme synoptic diagram.Wide spectrum light source produces flashlight, becomes the light (rise and be biased on fast axle or the slow axis) of linear polarization after the polarizer rises partially, is divided into through the 3db beam splitter then that two-way power equates, the identical light in polarization direction.Two-way light by circulator 1 and circulator 2, passes through bundling device respectively then through identical or different phase places or amplitude premodulation mode.After bundling device can be realized one road light signal polarization rotated 90 ° with another road combiner.The light that comes out from bundling device becomes the fast axial light and the slow axial light of quadrature.Polarization-maintaining coupler and polarization maintaining optical fibre can be so that fast axial light and slow axial light transmit as the flashlight in the optical fibre gyro of minimum reciprocal structure respectively.Whole optical path after the circulator is reversible.Since the identical noise factor of experience in transmission course, and the Sagnac phase place is identical.If fast axle and slow axis refractive index are respectively n 1And n 2, then phase place can be expressed as form:
Δφ 1=Δφ s+Δφ n(n 1)
Δφ 2=Δφ s+Δφ n(n 2) (2)
By different modulation formats, we can carry out input at the port 3 of two circulators and the nonreciprocal port of coupling mechanism.According to experiment Device Design, detect the interference light intensity of fast axle and slow axis respectively at circulator 1 and circulator 2, carry out nonreciprocal signal light power detection at the nonreciprocal port of polarization-maintaining coupler.
Amplitude modulation(PAM) sign indicating number type as shown in Figure 3.Only on fast axle (or slow axis), there is the light signal input when amplitude modulation(PAM) makes, promptly gets back to single reciprocity structure.But the depolarized characteristic of device is undesirable to make the light on the fast and slow axis have certain conversion.The signal section of slow axis is mainly changed over by following mode from fast axle: (1) forward direction is (2) afterwards to by the 3db coupling mechanism time (coupling mechanism accomplish two-way polarization property all very height ratio is difficult) during by the 3db polarization-maintaining coupler; (3) back is to by beam splitter the time.Still adopt two-way to detect, just can obtain being transformed into the splitting ratio of slow axis (or fast axle) from fast axle (or slow axis).Obtain the overall characteristic of these several light-splitting device groups, method is provided for the correction of the angular velocity of rotation of surveying.
Make fast and slow axis have the process that clear signal is handled the orthogonal signal light simultaneously with modulation format below.
Can be expressed as from the interference signal of circulator output:
I i=I 0i{1+cos[φ S+Δφ i(t)]} (3)
I=1 wherein, 2, expression fast axle pattern and slow axis pattern respectively.The φ of two patterns sIdentical, Δ φ i(t) phase differential of introducing for the PZT phase modulation (PM), because fast and slow axis refractive index difference, the phase differential that the PZT phase-modulator is introduced has minute differences.When utilizing sine wave modulation,
Δφ i(t)=2φ 0isin(ω mτ i/2)cos[ω m(t-τ i/2)] (4)
Wherein, ω mExpression PZT phase-modulator modulating frequency, m is the meaning of modulation, the modulation of expression PZT phase-modulator.(4) are brought into equation (3), can obtain according to the interference signal that Bessel's function launches be:
I i=I 0i{1+[J 0im)+2∑(-1) nJ 2nim)cos2nω m(t-τ i/2)]cos(φ s)+ (5)
2∑(-1) n+1J 2n-1m)sin[(2n-1)ω m(t-τ i/2)]sin(φ s)}
Wherein, φ Im=2 φ 0iSin (ω mτ i/ 2).
Utilize that wave filter takes out wherein one, two, the amplitude of four-time harmonic is respectively:
B i1=2I 0iJ 1im)sinφ s
B i2=2I 0iJ 2im)cosφ s
B i4=2I 0iJ 4im)cosφ s
By B I2/ B I4=J 2Im)/J 4Im), can try to achieve φ ImValue; Again by B I1/ B I2=J 1Im)/J 2Im) tan (φ s), can further try to achieve Sagnac phase shift φ Si, obtain angular velocity of rotation Ω according to the geometric properties of ring i
The φ that obtains for synchronization S1And φ S2, it comprises because the phase shift that pure rotation causes, and the phase shift that causes of various noise.
Can obtain the data of two output ports of different polarization states a moment simultaneously, analyze and obtain angular velocity of rotation.During this state, suppose to adopt above analytical approach to obtain containing the angular velocity of rotation Ω of different noises respectively from two polarization directions 1And Ω 2As shown in Figure 4, at a time, Ω 1And Ω 2Two-dimensional coordinate on, actual angular velocity of rotation value should be on 45 ° of lines.With a bit projecting on 45 ° of lines on the coordinate, pass through Ω 1And Ω 2Unite judgement, compare single Ω 1Value can obtain angular velocity of rotation more accurately.
Utilize the dual-polarization optical fiber gyro structure of Fig. 2, can further explore and analyze the common noise characteristic, to reach the purpose of eliminating the noise item.
In the fibre optic gyroscope, thermonoise is an important noise source, and the Approximate Equivalent speed of rotation error that the thermal transient process causes is:
Ω E ( t ) = 1 DL n dn dT ∫ 0 L dT dt ( z ) ( L - 2 z ) dz - - - ( 6 )
Wherein, D is the fiber optic loop diameter, and L is a fiber optic loop length, and T is a temperature, and n is an effective refractive index.
When using the different polarized light of two-way, can try to achieve φ respectively 1mAnd φ 2m,
φ im=2φ 0isin(ω mτ i/2)
=2n iΔlsin(ω mn iL/2c)
When the PZT modulator is operated in the pi/2 working point, sin (ω mn iL/2c) get maximal value, be in n iChange insensitive zone.So φ 1m/ φ 2m=n 1/ n 2
Use temperature control circuit to regulate the PZT temperature, gather the interference light signal data, { the φ when finding the solution different temperatures Im, obtain the n under the different temperatures 1/ n 2Value, thus extrapolate dn iThe curve of/dT is determined because the two-way equivalence speed of rotation error ratio that temperature transient causes is
Figure G2009102374215D00052
Be designated as ρ (T).Detailed process is as described in the embodiment of back.
Adopt modulation format 2, promptly on two branch roads flashlight is arranged all, establish the speed of rotation that obtains according to two paths of signals and be respectively Ω 1And Ω 2, real rotational speed is Ω, so
1-Ω)/(Ω 2-Ω)=ρ(T) (7)
Can instead solve with this: Ω = Ω 1 - ρ ( T ) Ω 2 1 - ρ ( T ) - - - ( 8 )
Like this,, can carry out Conjoint Analysis, reach more accurate angular velocity of rotation value to reduce influence for the noise that the temperature drift of two-way polarized light causes.This mode can be rejected the influence that temperature drift brings effectively, obtains angular velocity of rotation value accurately, and has guaranteed stability.
Other effects that cause noise have different forms with the relation of refractive index.If have difference power between the light wave of backpropagation in the fiber optic loop, will there be non-linear Kerr effect, a nonreciprocity phase differential that causes.Nonreciprocity refringence between the reverse light wave of this two-way is
Δn k = δn NL 1 - δn NL 2 = χ e ( 3 ) ( | E cw | 2 - | E ccw | 2 ) 2 n - - - ( 9 )
Wherein, χ e (3)Being three rank electric susceptibilities, is the main cause that causes third harmonic, four-wave mixing and nonlinear refractive index, | E Cw| 2With | E Ccw| 2It is respectively clockwise and counterclockwise light intensity.
For different polarization states, (n can be similar to and think the n difference n = ϵ r = 1 + χ e )。The variation that temporary transient consideration is caused by Sagnac effect and non-linear Kerr effect, then
φ 1 = φ s + Δn 1 · φ = φ s + χ e ( 3 ) ( | E 1 cw | 2 - | E 1 ccw | 2 ) 2 n 1 · φ
φ 2 = φ s + Δn 2 · φ = φ s + χ e ( 3 ) ( | E 2 cw | 2 - | E 2 ccw | 2 ) 2 n 2 · φ - - - ( 10 )
Wherein φ = 2 π λ L , (L refers to integration lengths here, and for wideband light source, this integration lengths is positioned at very short one section of coil mid point).
For the nonreciprocity effect that non-linear Kerr effect causes, phase interference item and the polarization state institute corresponding refractive index relation of being inversely proportional to be not so its signal processing method of taking just can be the same.Than temperature drift, the noise ratio that non-linear Kerr effect causes is less, so need (float as controlling temperature under at experiment condition, or adopt Kalman filtering to handle two tributary signals respectively and eliminate temperature and float, what obtain be respectively Ω on two branch roads eliminating temperature noise influence back i'), consider to eliminate non-linear Kerr effect influence again.
Adopt amplitude modulation(PAM) sign indicating number type 2, promptly two branch roads all have flashlight, in the process that the experiment light path is built, record two-way light and close the extinction ratio δ that restraints to after the orthogonal optical on fast axle and the slow axis through bundling device, can extrapolate beam intensity ratio ξ=(1+ δ)/(the 1-δ) on fast axle and the slow axis.In equation (10),
ξ = | E 1 cw | 2 - | E 1 ccw | 2 | E 2 cw | 2 - | E 2 ccw | 2
Under the condition of modulation format 1, measured the n under the different temperatures in front 1/ n 2Value.Note η=ξ n 2/ n 1
Release the formula of similar (7) by equation (10):
1’-Ω’)/(Ω 2’-Ω’)=η (11)
Wherein, Ω 1' and Ω 2' be respectively two branch roads to be carried out angular velocity of rotation value after temperature is floated processing.
Thereby obtain Ω '=(η Ω 2'-Ω 1')/(η-1) (12)
For the noise that other factors cause, analyze the difference of its noisiness on fast and slow axis, finding algorithm pointedly just is to reduce The noise.
For thermal stress and vibration, the variation that this can cause the local birefringence rate causes the nonreciprocity biased error.Its noise error that causes is directly proportional with the variations in refractive index rate, can consider that following manner handles.
φ 1=φ s+Δn 1·φ (13)
φ 2=φ s+Δn 2·φ
Ω is easy to get s=(Ω 2Δ n 11Δ n 2)/(Δ n 1-Δ n 2) (14)
But, need on statistical significance, carry out possible correction, or the factor that causes thermal stress and vibration of taking measures as far as possible to reduce to the result because thermal stress and vibration be difficult to determine.
Because the influence of cross polarization effect, high to the characteristic requirement of polarization-maintaining coupler in the light path, to forward and reverse, each polarization axle all will reach the extinction ratio of 30db.Therefore can consider depolarized scheme, both reduce, also can reduce the cost of polarization maintaining optical fibre the device property requirement.
The light channel structure of depolarized scheme as shown in Figure 5.
The power of two-way light can be different, carrying out different premodulation by two light paths behind the polarization-maintaining beam splitter, can be phase adjusted or amplitude adjusted promptly.
Good effect of the present invention is:
For dual-polarization optical fiber gyro scheme, carry out different amplitude sign indicating number type modulation, according to noisiness relevant with refractive index on the fast and slow axis, The noise is eliminated in Conjoint Analysis, to obtain angular velocity of rotation accurately.
Description of drawings
Fig. 1 is traditional single shaft polarization maintaining optical fibre gyro minimum reciprocal structure;
Fig. 2 is the double-polarization interference type optical fiber gyroscope structural drawing;
Fig. 3 is to dual-polarization light path amplitude modulation(PAM) sign indicating number type;
Fig. 4 is simple signal Processing schematic diagram;
Fig. 5 is the double-polarization interference type optical fiber gyroscope that adopts the depolarized structure of single mode.
Embodiment
Below in conjunction with the drawings and specific embodiments the present invention is described in further detail:
The present invention adopts ASE wide spectrum light source output signal light, becomes linearly polarized light behind the polarizer, and the polarization-maintaining coupler that adopts 50: 50 is as beam splitter.Amplitude premodulation on the two-way that is divided into adopts optical attenuator, and circuit is controlled the two-way attenuator state simultaneously, to reach required modulation format.Two kinds of amplitude modulation(PAM) sign indicating number types commonly used in experiment are: type 1---have only on the branch road has light signal to pass through to sign indicating number, and no light signal on another; Sign indicating number type 2---all there is light signal to pass through on two branch roads.Protect inclined to one side circulator and realize that forward is by (port one->port 2), oppositely by (port 2->port 3).Protect shift/beam splitter and the light polarization direction on the branch road can be rotated 90 °, and with two branch roads photosynthetic be the flashlight of quadrature on fast and the slow axis.50: 50 polarization-maintaining couplers with High Extinction Ratio are input to the polarization-maintaining fiber coil two ends with the flashlight of quadrature.
Port 3 at two circulators of reverse optical path is converted to electric signal to the interference light signal of different polarization states with photoelectric commutator, and gathers by the NI data collecting card.Utilize LabView that the data of gathering are carried out spectrum analysis and data storage.
Before experiment, want earlier to make having only that (as above branch road) has flashlight on the branch road by modulation format 1, by the temperature that temperature-control circuit is regulated PZT, store the interference signal luminous power of this branch road circulator port 3 after the PD conversion with data collecting card, utilize second harmonic B 2With four-time harmonic B 4Calculate corresponding { φ 1m; Under this temperature, making to have only has flashlight (following branch road) on another branch road, obtain { φ according to same step 2m.So just can obtain the n under the different temperatures 1/ n 2Value, and extrapolate dn iThe curve of/dT.For the polarization property of entire system, need be at this modulation format 1 time, the signal data of circulator port 3 on another branch road of collection analysis.
When measuring the Sagnac effect, use modulation format 2, obtain two interference signal light signals on the branch road respectively, demodulation gets speed of rotation Ω 1And Ω 2After, utilize equation (8) to try to achieve Ω value more accurately.
The depolarizer that adopts in the depolarized structure of single mode can be the Lyot depolarizer of 45 ° of formation of polarization maintaining optical fibre deflection of 1: 2: 4 with three segment length ratio.Fiber optic loop adopts the general single mode fiber ring to get final product.Behind polarization-maintaining beam splitter, be divided into different polarization states in the reverse optical path, and carry out check and analysis.Because the various noises that two polarization states that are divided into experience in polarization-maintaining fiber coil more become identical with refractive index, its signal processing method has the need for further discussion.

Claims (9)

1. the united signal processing method of a double-polarization interference type optical fiber gyroscope, wherein said double-polarization interference type optical fiber gyroscope is:
Two linearly polarized light input light paths are protected shift/beam splitter through one and are connected with a coupling mechanism, and described coupling mechanism opposite side is connected with a fiber optic loop, are connected a phase-modulator between described fiber optic loop and the described coupling mechanism, the steps include:
1) modulation double-polarization interference type optical fiber gyroscope makes the two described linearly polarized light input light paths output signal light of property at interval;
2) measure described phase-modulator when different temperatures, the luminous power of exporting on the two described polarized light input light paths;
3), determine because the two-way equivalence speed of rotation error ratio that temperature transient causes according to the temperature value and the corresponding luminous power of record
Figure FSB00000496643800011
N wherein 1Be fast axle refractive index, n 2Be the slow axis refractive index, T is a temperature;
4) modulation double-polarization interference type optical fiber gyroscope makes two described linearly polarized light input light paths output signal light simultaneously;
5) measure the luminous power of exporting on the two described linearly polarized light input light paths;
6) calculate speed of rotation Ω respectively according to step 5) institute photometry power 1And Ω 2
7) according to formula
Figure FSB00000496643800012
Calculate revised rotational speed Ω.
2. the method for claim 1 is characterized in that described coupling mechanism is a polarization-maintaining coupler, and described fiber optic loop is a polarization-maintaining fiber coil.
3. method as claimed in claim 1 or 2 is characterized in that described linearly polarized light input light path comprises a signal modulator, and the linear polarization light source is connected with described guarantor's shift/beam splitter through described signal modulator.
4. method as claimed in claim 3 is characterized in that described linearly polarized light input light path comprises a circulator, and described signal modulator is connected with described guarantor's shift/beam splitter through described circulator, and described circulator is connected with a detector.
5. method as claimed in claim 4 is characterized in that described signal modulator is a phase-modulator.
6. method as claimed in claim 4 is characterized in that described signal modulator is an amplitude modulator.
7. method as claimed in claim 1 or 2 is characterized in that described ratio
Figure FSB00000496643800013
Computing method be: at first the expression formula of interference signal luminous power is launched according to Bessel's function; Then according to step 2) luminous power surveyed, utilize that wave filter takes out wherein one, two, the amplitude of four-time harmonic; Then according to extracted one, two, φ when the amplitude of four-time harmonic solves different temperatures Im, φ wherein Im=2 φ 0iSin (ω mτ i/ 2); Obtain the n under the different temperatures at last 1/ n 2Value, thus extrapolate dn iThe curve of/dT is determined described ratio ρ (T); I=1 or 2 wherein, during i=1 corresponding to fast axle, during i=2 corresponding to slow axis.
8. method as claimed in claim 1 or 2 is characterized in that adopting kalman filter method to handle the signal of two described linearly polarized light input light paths respectively, eliminates temperature and floats, and obtains the speed of rotation Ω on two branch roads 1' and Ω 2'; Then according to formula Ω '=(η Ω 2'-Ω 1')/(η-1) calculate revised speed of rotation Ω ', η=ξ n wherein 2/ n 1,
Figure FSB00000496643800021
| E i Cw| 2Be clockwise light intensity, | E i Ccw| 2Be counterclockwise light intensity, i=1 or 2, during i=1 corresponding to fast axle, during i=2 corresponding to slow axis.
9. method as claimed in claim 1 or 2 is characterized in that utilizing formula Ω s=(Ω 2Δ n 11Δ n 2)/(Δ n 1-Δ n 2) calculate revised speed of rotation Ω s, Δ n wherein 1For causing that by stress or vibration fast axle local indexes of refraction changes Δ n 2For causing that by stress or vibration the slow axis local indexes of refraction changes.
CN2009102374215A 2009-11-06 2009-11-06 United signal processing method of double-polarization interference type optical fiber gyroscope Expired - Fee Related CN101718557B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009102374215A CN101718557B (en) 2009-11-06 2009-11-06 United signal processing method of double-polarization interference type optical fiber gyroscope

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009102374215A CN101718557B (en) 2009-11-06 2009-11-06 United signal processing method of double-polarization interference type optical fiber gyroscope

Publications (2)

Publication Number Publication Date
CN101718557A CN101718557A (en) 2010-06-02
CN101718557B true CN101718557B (en) 2011-07-27

Family

ID=42433164

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009102374215A Expired - Fee Related CN101718557B (en) 2009-11-06 2009-11-06 United signal processing method of double-polarization interference type optical fiber gyroscope

Country Status (1)

Country Link
CN (1) CN101718557B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607547A (en) * 2012-02-17 2012-07-25 北京大学 Digital domain balanced detecting method and device for optical fiber gyroscope

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101886925B (en) * 2010-05-27 2012-06-27 北京大学 Multi-wavelength interference type optical fiber gyro based on carrier modulation
CN102042849B (en) * 2010-12-23 2012-03-21 北京大学 Truth value estimation method based on multi-dimensional vector information and device thereof
CN102288388B (en) * 2011-05-09 2013-04-10 哈尔滨工程大学 Device and method for improving polarization-maintaining optical fiber polarization coupling measurement precision and symmetry
CN102393198B (en) * 2011-10-20 2015-04-08 武汉虹拓新技术有限责任公司 Optical gyroscope
CN102494681B (en) * 2011-12-06 2014-03-12 北京航空航天大学 Difference double-interference type optical fiber gyroscope based on birefringence modulation
CN103399191B (en) * 2013-07-25 2016-06-22 哈尔滨理工大学 FBG-GMM electric current sensing method based on sideband demodulation
CN103674070B (en) * 2013-12-26 2017-05-31 北京航天时代光电科技有限公司 A kind of fiber-optic gyroscope light path mismachining tolerance tests update the system and method
CN104075705A (en) * 2014-06-26 2014-10-01 中航捷锐(北京)光电技术有限公司 Optical fiber gyroscope for improving velocity sensitivity
CN105674975B (en) * 2016-01-07 2019-02-01 西安邮电大学 A kind of single axis fiber gyro instrument
CN106338283A (en) * 2016-09-27 2017-01-18 哈尔滨工程大学 High-precision interference type optical fiber gyroscope temperature drift error model optimizing method
CN107576341B (en) * 2017-08-09 2021-06-04 武汉昊衡科技有限公司 Device and method for eliminating polarization fading in OFDR (offset frequency domain digital radiography)
CN109556593B (en) * 2018-10-09 2020-08-04 清华大学深圳研究生院 Angular velocity measuring device, method and carrier thereof
CN109724585A (en) * 2019-01-28 2019-05-07 北京航空航天大学 Device is inhibited using the light source relative intensity noise of double polarization-maintaining couplers
CN110006418B (en) * 2019-04-15 2020-12-22 清华大学深圳研究生院 Angular velocity measuring method, device and application system
CN111457911B (en) * 2020-05-29 2022-02-01 北京航空航天大学 Bionic polarization compass calibration method based on polarization two-dimensional residual error information
CN112051411A (en) * 2020-08-07 2020-12-08 北京航天控制仪器研究所 Angular motion measuring device based on nonlinear optical wave beam splitting and beam combining
CN112051606B (en) * 2020-09-10 2024-04-02 北京大学 Six-component seismograph
CN112083477B (en) * 2020-09-10 2024-03-19 北京大学 Three-component rotary seismograph
CN112797971B (en) * 2020-12-11 2022-05-27 浙江大学 Differential fiber optic gyroscope based on characteristic of inhibiting temperature drift of double-core fiber

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101482412A (en) * 2009-02-24 2009-07-15 北京航天时代光电科技有限公司 Closed-loop optic fiber gyroscope light path structure with low polarization error

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101482412A (en) * 2009-02-24 2009-07-15 北京航天时代光电科技有限公司 Closed-loop optic fiber gyroscope light path structure with low polarization error

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP特开2002-54933A 2002.02.20

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607547A (en) * 2012-02-17 2012-07-25 北京大学 Digital domain balanced detecting method and device for optical fiber gyroscope
CN102607547B (en) * 2012-02-17 2014-07-09 北京大学 Digital domain balanced detecting method and device for optical fiber gyroscope

Also Published As

Publication number Publication date
CN101718557A (en) 2010-06-02

Similar Documents

Publication Publication Date Title
CN101718557B (en) United signal processing method of double-polarization interference type optical fiber gyroscope
CN101629825B (en) Dual-polarization interferometric fiber-optic gyro
US11268811B2 (en) Non-interferometric optical gyroscope based on polarization sensing
CN101660910B (en) Low-noise dual-polarization interference optic fiber gyroscope
CN108534798B (en) Polarization nonreciprocal error elimination method in dual-polarization fiber-optic gyroscope and dual-polarization fiber-optic gyroscope
CN101825465B (en) Reentrant interference optical fiber gyroscope
CN101886925B (en) Multi-wavelength interference type optical fiber gyro based on carrier modulation
CN102128621A (en) Interference optical fiber gyroscope for measuring multidimensional vector
CN102706340B (en) A kind of interferometric fiber optic gyroscope
CN109883412A (en) A kind of double path optical fiber gyroscopes
KR101978444B1 (en) Optical fiber Sagnac interferometer using a polarizing beam splitter
CN103438882B (en) A kind of optical fibre gyro of low scale factor error
US20030123064A1 (en) Symmetrical depolarized fiber optic gyroscope
Li et al. Test and analysis of the optical Kerr-effect in resonant micro-optic gyros
CN102607547B (en) Digital domain balanced detecting method and device for optical fiber gyroscope
CN103292976A (en) Interferometry-based optical fiber Verdet constant measuring system
CN104075705A (en) Optical fiber gyroscope for improving velocity sensitivity
CN102607589A (en) Method and device for measuring angular speed of optical fiber gyroscope based on double-frequency modulating signals
CN103411600B (en) A kind of two interferometric fiber optic gyroscopes ± pi/2 measurement method of parameters based on light path difference
Wang et al. Optimization of the sinusoidal phase modulation technique in resonant fiber optic gyro
CN107976300B (en) Method for measuring beat length of polarization maintaining optical fiber
CN103033178B (en) Laser gyroscope offset frequency method based on external cavity feedback
Qi et al. An ultra-short coil fiber optic gyroscope
Kim et al. Fiber-optic gyroscopes: In harsh, confining environments this advanced gyroscope, a close cousin to the ring laser gyro, offers great advantages
CN107806981B (en) Measuring device for beat length of polarization maintaining optical fiber

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110727

Termination date: 20131106