Apparatus and Method for Integrating Continuous and Discontinuous Inertial Instrument

Abstract

The invention is related to a method and an inertial navigation system for combining continuous signal output from a first inertial sensor (14) with discontinuous signal output from a second inertial sensor (12). The first inertial sensor (14) acquires continuous data with respect to a navigation frame of reference for a parameter used in inertial navigation and the continuous data is processed to produce estimated values of the parameter. The second inertial sensor (12) acquires discontinuous data with respect to a case frame of reference indicative of the parameter with respect to a case (25) containing the second inertial sensor (12). The discontinuous data is processed to produce measurements of the parameter at selected times, —and the estimated values of the parameter and the measurements of the parameter are processed at selected times with a Kalman filter to provide corrections to the estimated values of the parameter at the selected times.

Claims

1. A method for combining continuous signals output from a first inertial sensor (14) with discontinuous signals output from a second inertial sensor (12), comprising the steps of: using the first inertial sensor (14) to acquire continuous data with respect to a navigation frame of reference fixed with respect to the first inertial sensor (14) for a parameter used in inertial navigation; processing the continuous data to produce estimated values of the parameter; arranging the second inertial sensor (12) to acquire discontinuous data with respect to a case frame of reference indicative of the parameter with respect to a case (25) that contains the second inertial sensor (12); processing the discontinuous data to produce measurements of the parameter at selected times; and processing the estimated values of the parameter and the measurements of the parameter at selected times with a Kalman filter to provide corrections to the estimated values of the parameter the selected times.

2. The method of claim 1, further comprising the step of arranging the second inertial sensor (12) to have a greater accuracy than the first inertial sensor (14).

3. The method of claim 1, further comprising the steps of: forming the second inertial sensor (12) to comprise an atom optic instrument (18) that is arranged to establish a positional reference in inertial space using a cloud of ultra-cold atoms that move within a (25) for a finite time period; obtaining positional data in the navigation frame for the cloud of atoms as a function of time; determining a position vector for the atom cloud in the navigation reference frame at a time when the atom cloud is released into the (25); forming a scalar product of the positional data in the navigation frame or the cloud of atoms as a function of time.

4. An inertial navigation system (10), comprising: a central processing unit (16); a first inertial sensor (14) arranged to provide a continuous signal indicative of a parameter to the central processing unit (16), the central processing unit (16) being arranged to process the continuous signal to obtain estimated values of the parameter; and a second inertial sensor (12) arranged to provide discontinuous signal indicative of the parameter to the central processing unit (16), the central processing unit (16) being arranged to process the discontinuous signal to obtain measurements of the parameter at selected times, the central processing unit (16) being further arranged to process the estimated values of the parameter and the measurements of the parameter at selected times with a Kalman filter to provide corrections to the estimated values of the parameter the selected times.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0009] FIG. 1 is a block diagram showing a conventional navigation system, a set of atom optic instruments and a central processing unit that may be included in the invention;

[0010] FIG. 2 shows the basic structure of an atom optic interferometer that may be included in the set of atom optic instruments of FIG. 1;

[0011] FIG. 3 shows coordinate systems that may be used for inertial navigation; and

[0012] FIG. 4 is a functional block diagram showing integration of measurements output from an atom optic instrument with a conventional inertial sensing instruments.

DETAILED DESCRIPTION OF THE INVENTION

[0013] Referring to FIG. 1, the present invention is directed to an inertial navigation system 10 in which highly accurate inertial sensors, such as an atom optic instruments 12 are used to provide periodic corrections to a conventional navigation system 14 to provide a continuously operating navigation system having improved accuracy.

[0014] The atom optic instruments 12 provides periodic, high accuracy inertial measurements to a central processing unit (CPU) 16. The CPU 16 also receives navigation information from the conventional navigation system 14. The CPU processes signals received from the atom optic instruments 12 and from the conventional navigation system 14 in a Kalman filter to provide a correction to data received from the conventional navigation system 14.

[0015] Atom optic instruments are disclosed in U.S. Pat. No. 5,274,232, issued Dec. 28, 1993 to Chu et al. and U.S. Pat. No. 4,992,656 issued Feb. 12, 1991 to Clauser, the disclosures of which are hereby incorporated by reference into this disclosure.

[0016] As disclosed in U.S. Pat. No. 5,274,232 an atom optic instrument typically includes an evacuated chamber having a plurality of windows. Cooling laser beams are admitted into the chamber through a first group of the windows. The cooling laser beams converge on a target region within the cell and are included in a laser cooling system. As is known in the art a laser cooling system reduces the velocity, and hence the temperature of an atom, by reducing its momentum as the atom absorbs photons in Raman transitions. The cooling laser beams cooperate with a magnetic field to trap atoms, typically Cesium atoms, within the target region for a time of sufficient duration to permit measurements to be made.

[0017] FIG. 2 illustrates an atom optic interferometer 18 that may be included in the set of atom optic instruments 12. The invention should not be understood to be limited to any particular high accuracy navigation system and should not be understood to be limited to the specific atom optic interferometer 18 structure shown in FIG. 2.

[0018] Still referring to FIG. 2, a Cesium oven 20 produces a beam of Cesium atoms that is collimated by two 1 mm diameter apertures (not shown) and then cooled using laser cooling techniques. A laser 22 provides a Raman population transfer beam to the Cesium atoms, which undergo stimulated Raman transitions between the Cs 6S.sub.1/2, F=3, m.sub.F=0 and the 6S.sub.1/2, F=4, m.sub.F=0 magnetic field insensitive ground state hyperfine levels. The resulting atomic beam has a most probable longitudinal velocity of 290 m/s and a transverse velocity of 10 cm/s. State preparation of the Cesium atoms is accomplished by first optically pumping the atoms into the F=4, m.sub.F=0 level. The atoms are then transferred from the F=4, m.sub.F=0 level to the F=3, m.sub.F=0 level with a Doppler insensitive it Raman pulse. A laser 24 provides a blaster beam that removes atoms that remain in the F=4 state from the Cesium beam.

[0019] After state preparation is completed, atoms enter a magnetically shielded case 25 that encloses an interrogation region 26. Magnetic shielding may be provided by passing electric current through a pair of conducting coils (not shown). In the interrogation region 26 three pairs of counterpropagating laser beams are used to drive additional Raman transitions. The laser beams are oriented to be perpendicular to the atomic beam in a horizontal plane and are provided by lasers 28-33. The Cesium atoms enter the interrogation region 26 in the |F=3, m.sub.F=0; p.sub.0 state, where p.sub.0 is the initial momentum per atom. A first set of π/2 Raman pulses from the lasers 28 and 31 puts most of the Cesium atoms in a superposition of the |F=3, m.sub.F=0; p.sub.0 state and the |F=4, m.sub.F=0; p.sub.0+2k state where 2k is the recoil momentum of each atom that has absorbed two photons and experienced a Raman transition. This recoil momentum corresponds to a 7 mm/s velocity transverse to the initial velocity of the Cesium atoms, which causes these atoms to separate from atoms remaining in their initial states upon entering the interrogation region as the atoms continue their trajectories. In a free-flight distance of about 1 m, the wave packets corresponding to the two atomic states separate by about 23 μm.

[0020] Next the atoms pass through a pair of π Raman beams provided by the lasers 29 and 32, which interact with the Cesium atom states and exchange the two ground states and momenta, which reverses the direction of the atoms that were deflected by the π/2 Raman beam. The atom beams then combine after traveling another meter where a second pair of π/2 Raman pulses overlaps the two wave packets. Atoms that made the transition from the F=3 to F=4 state are detected by a resonant probe laser beam from the laser 33 that is tuned to the 6S.sub.1/2, F=4.fwdarw.6S.sub.1/2, F=5 cycling transition. A lens system 36 focuses the fluorescence resulting from this transition onto a photomultiplier tube, which produces a corresponding electrical signal.

[0021] Dividing the input beam of Cesium atoms into two beams that are subsequently combined as described above produces a Mach-Zehnder interferometer, which is well-known in both bulk optics and fiber optics. Rotation of the atom optic interferometer about an axis perpendicular to the plane defined by the two paths followed by the wave packets produces phase shift between the wave packets in accordance with the Sagnac effect. The prior art is replete with examples of ring laser and fiber optic sensors that use the Sagnac effect to measure rotations.

[0022] FIG. 3 shows a navigation coordinate system having orthogonal axes X.sub.N, Y.sub.N, Z.sub.N that are transported with a vehicle relative to the earth where X.sub.N, Y.sub.N are local-level at the wander angle α relative to the east and north geodetic axes and Z.sub.N is along the local vertical fixed in an object such as an aircraft or a missile. FIG. 3 also shows an earth-fixed coordinate system having orthogonal axes X.sub.E, Y.sub.E, Z.sub.E that are fixed relative to the earth. An inertial coordinate system has orthogonal axes X.sub.1, Y.sub.1, Z.sub.1 that are fixed in inertial space. In a typical navigation situation, the navigation coordinate system moves relative to the earth-fixed coordinate system while at the same time, the earth-fixed coordinate system moves relative to the inertial coordinate system. Inertial navigation typically involves transforming measurements made in one coordinate system such as the vehicle body coordinate system into one or more of the three coordinate systems shown.

[0023] FIG. 4 is a functional block diagram of the invention, which is based on assuming that the periodic measurement from the atom optic (AO) instrument is a phase difference measurement of the form:

Δφ=[φ.sub.1−φ.sub.2+φ.sub.3],  (1)

where the phases φ.sub.i (i=1, 2, 3) are estimated as the scalar product of the k vector, which is the propagation direction for the laser field with the vectors representing the difference in position x.sub.P(t) of the atom cloud falling under the influence of gravity and its initial velocity and the position of a reference point fixed to the atom optic instrument case x.sub.R (t), such reference point being determined by the optical configuration of the instrument, at the times t.sub.i=t.sub.0, t.sub.0+T, t.sub.0+2T, where i=1, 2, 3 and where 2T is a fraction of a second. The phase may be written explicitly as:

φ.sub.i=k(t.sub.i).Math.[x.sub.P(t.sub.i)−x.sub.R(t.sub.i)]  (2)

[0024] The atom optic instrument case 25 has an orthogonal reference set of axes denoted [x.sub.C,y.sub.C,z.sub.C]. The origin of this reference frame is the case-fixed reference point referred to above. The point of release of the atom cloud during the measurement periods is assumed to be fixed in the atom optic instrument case reference axes and is denoted by the fixed “lever arm.”

L=[L.sub.xc,L.sub.yc,L.sub.zc].  (3)

[0025] The vector k is fixed in the case 25 of the atom optic instrument with coordinates [k.sub.xc,k.sub.yc,k.sub.zc] for all time indices i.

[0026] The navigation coordinate system maintained by the conventional inertial navigation system (INS) can be defined to be coincident with the local East, North and Vertical [E, N, V] axes at the present position (latitude and longitude of the INS with respect to the earth. The vertical can be defined to be normal to the ellipsoidal model of the earth. The origin of this navigation coordinate system is a point with respect to which the measurements of force by the conventional accelerometers and angular rotation by the conventional gyros are referred. In general this origin will be separated from the case-fixed reference point of the atom optic instrument by a fixed lever arm that can be expressed in the case coordinate axes of the atom optic instrument or in an orthogonal reference set of axes defined by the conventional inertial instruments. These two orthogonal reference coordinate systems will be related to each other by a fixed transformation.

[0027] To simplify the following description, it is assumed that the transformation between the coordinate axes defined by the INS sensors and the atom optic instrument case 25 is the identity matrix. It is further assumed that the origins of these two coordinate systems are identical such that the position of the reference point of the atom optic instrument will also be the reference point for navigation using the conventional INS measurements.

[0028] The following equation is a simplified representation of the trajectory of the atom cloud with respect to an inertial frame during the measurement period:

[00001]$\begin{array}{cc}{x}_P\left(t\right)=x_0+\mathrm{Vt}+\left[\frac{1}{2}\right]g{t}^2& \left(4\right)\end{array}$

where: [0029] x.sub.o is the initial position of the atom cloud; [0030] V is the initial velocity of the atom cloud with respect to the atom optic case 25 where it is assumed for simplification that the atom optic case 25 has zero velocity with respect to inertial space; and [0031] g is the acceleration of gravity for a simplified flat earth, “no rotation” approximation. These simplifying assumptions are now relaxed for the practical situation below.

[0032] The frame of reference for the trajectory of the atom cloud can be the reference navigation coordinates [X.sub.N, Y.sub.N, Z.sub.N] defined in FIG. 3 above in which the trajectory of the reference point of the atom optic instrument is expressed or the atom optic case coordinate axes [x.sub.c,y.sub.c,z.sub.c]. Without loss of generality the wander angle α can be set to zero such that the X.sub.N axis points east, E and the Y.sub.N axis points north, N. The trajectories of the atom cloud and the optical facet can be expressed in either reference frame using a transformation .sub.NT.sub.C between them. The transformation .sub.NT.sub.C is continuously computed by the INS using the difference of the conventional gyro measured angular rates ω.sub.g with respect to inertial space measured along case coordinates, and the INS computed angular rate of the navigation coordinates with respect to inertial space, ω.sub.Nc expressed in case coordinates.

[0033] Since the relative position of the atom cloud upon release at t=0, is known in the case coordinates as is the k-vector, then the phase φ.sub.i is a fixed number known a′ priori, namely,

φ.sub.1=k.Math.L=[k.sub.xc,k.sub.yc,k.sub.zc].Math.[L.sub.xc,L.sub.yc,L.sub.zc].  (5)

[0034] FIG. 4 shows how to develop an expression of the atom cloud trajectory in navigation coordinate axes. The gravity vector g.sub.N is fixed in the navigation coordinate axes along the local vertical axis, Z.sub.N except for gravity disturbances with respect to the earth. The gravity disturbance vector is ignored at this point but can be readily introduced into the mechanization later if it is known. The variable g.sub.N is defined to be the value of gravity along the local vertical at a current position of the vehicle transporting the INS. The gravity vector g.sub.N is the combination of the acceleration of gravity due to mass attraction and centripetal acceleration due to earth rotation Q with respect to inertial space.

[0035] The initial velocity V.sub.c of the atom cloud relative to the atom optic case 25 upon release is defined in case coordinate axes as:

V.sub.C=[V.sub.X,V.sub.Y,V.sub.Z]  (6)

[0036] The velocity V.sub.c/E of the origin of the atom optic instrument case axes relative to the earth expressed in the navigation coordinate axes as a function of time is defined as:

V.sub.C/E=[V.sub.XN(t),V.sub.YN(t),V.sub.ZN(t)]  (7)

The relative velocity of the point of release of the atom cloud with respect to the origin (reference point) of the atom optic instrument case 25 with respect to the earth is expressed as a vector cross product of the lever arm between the two points L, and the relative angular rate of the case 25 with respect to the earth (ω.sub.g−ω.sub.Nc). This term is written as L×(ω.sub.g−ω.sub.Nc) as shown in FIG. 4.

[0037] Consequently the initial velocity V.sub.Ac/E(0).sub.N of the atom cloud with respect to the earth expressed in the navigation coordinate axes is

V.sub.AC/E(0)=.sub.NT.sub.C(0)+V.sub.C/E(0).sub.N+L×(ω.sub.NC(0)−ω.sub.g(0)),  (8)

where θ.sub.N, is the angular rate of the navigation coordinate axes with respect to the earth but expressed in case coordinates and ω.sub.g is the angular rate of the case 25 with respect to inertial space measured in case coordinates.

[0038] The subsequent position of the atom cloud with respect to the earth expressed in the navigation coordinate axes is then estimated using the conventional inertial navigation solution information over the 2T measurement period as:

P.sub.AC/N(t).sub.N=∫{∫[g.sub.N+C.sub.N]dμ+V.sub.AC/E(0).sub.N}dν+P.sub.AC/E(0).sub.N  (9)

[0039] The initial position of the atom cloud with respect to the earth is obtained from the initial position of the reference point of the atom optic instrument P.sub.C/E(0).sub.N with respect to the earth in the navigation coordinate axes plus the lever arm L between the two positions expressed in the navigation coordinate axes.

P.sub.AC/E(0).sub.N=P.sub.C/E(0).sub.N+.sub.NT.sub.C(0)[L].sub.C.  (10)

[0040] Since g.sub.N represents acceleration with respect to inertial space expressed in the navigation coordinate axes, a Coriolis correction C.sub.N, is required to obtain the time derivative of velocity with respect to the earth with respect to the navigation coordinate axes in which the integration takes place. This Coriolis term in the navigation system coordinate frame is expressed as:

C.sub.N=−[ω.sub.NC+Ω].sub.N×V.sub.AC/E(t)  (11)

where:

[0041] Ω is the Earth's Rate of Rotation with respect to inertial Space, and ω.sub.Nc is the rate of rotation of the navigation coordinate axes with respect to inertial space.

[0042] Therefore, the velocity of the atom cloud relative to the earth is estimated as:

V.sub.AC/E(t).sub.N=∫[g.sub.N+C.sub.N]dμ+V.sub.AC/E(0).sub.N  (12)

[0043] From the above, the INS provides the navigation solution for the reference point of the atom optic instrument P.sub.C/E(t).sub.N from its initial position P.sub.C/E(0).sub.N at the start of the measurement period, 2T.

[0044] FIG. 4 also shows an expression for the k-vector in navigation coordinate axes. Since the transformation between the navigation coordinate axes and the atom optic case coordinate axes is known, the k-vectors at the times T and 2T in the navigation coordinate axes are:

k.sub.2(T).sub.N=.sub.NT.sub.C(T)└k.sub.xc,k.sub.yc,k.sub.zc┘  (13)

and

k.sub.3(2T).sub.N=.sub.NT.sub.C(2T)└k.sub.xc,k.sub.yc,k.sub.zc┘  (14)

[0045] The phase angles may be estimated as

φ.sub.2=k.sub.2(T).sub.N.Math.[P.sub.AC/E(T)−P.sub.OF/E(T)].sub.N  (15)

and

φ.sub.3=k.sub.2(2T).sub.N.Math.[P.sub.AC/E(T)−P.sub.OF/E(T)].sub.N  (16)

[0046] A block diagram depicting the information flow for this mechanization is shown in FIG. 4. The Kalman filter processes a Kalman observation that may be expressed as

Kalman observation=Δφ.sub.AO instrument−Δφ.sub.Estimated  (17)

which is the difference of the atom optic instrument measurement Δϕ.sub.AO instrument, and the estimate of this measurement Δϕ.sub.Estimated, that is computed using the information from the conventional inertial navigation solution.

[0047] The means for correction of the high bandwidth conventional navigation solution obtained with conventional inertial instruments is also shown in FIG. 4. The periodic discontinuous measurement of the phase difference provided by the atom optic instrument Δφ.sub.AO Instrument is compared with an estimate of the phase difference computed using the conventional inertial navigation solution measurements Δφ.sub.Estimated. The difference of these two measurements is the discrete observation that is processed by a Kalman filter to correct all error states modeled in the Kalman filter. These error states pertain to all significant errors in the conventional inertial navigation solution and the conventional inertial instruments and also all significant error states in the atom optic instrument measurements. For example, the Kalman filter may provide corrections to the conventional inertial gyro and accelerometer biases, to attitude errors in .sub.NT.sub.C, to velocity errors, position errors, and other modeled error sources.

[0048] The embodiment of the invention described is meant as an example. Other similar embodiments can also be arrived at by those skilled in the art. For example, the method described could be used with sensors having a different form of output as that given in equations (1) and (2). In this case the conventional navigation solution would be used to compute an estimate of the appropriate mathematical form of the atom optic instrument measurement. Also, navigation could be performed in a wide variety of coordinate systems.