INS Error Propagation (IEEE Format)
INS Error Propagation (IEEE Format)
_
_ _ _
+ + + +
, , ,
,
$
B
$
S
$
w
$
M
v
B
meter per second s(uare 0m/s1 or g2. The second is the
bias instabi!ity (for gyro or accelerometer), which is the
random variation in the !ias as computed over specified
finite sample time and averaging time intervals. This non-
stationary (evolutionary) process is characteri#ed !y a ,/f
power spectral density. It is typically expressed in degree
per hour (./hr) or meter per second s(uare 0m/s1 or g2,
respectively.
(.&.( S'ale 2a'tor
Scale factor is the ratio of a change in the input
intended to !e measured. Scale factor is generally evaluated
as the slope of the straight line that can !e fit !y the method
of least s(uares to input-output data. The scale factor error
is deterministic in nature and can !e determined !y
cali!ration. Scale factor asymmetry (for gyro or
accelerometer) is the difference !etween the scale factor
measured with positive input and that measured with
negative input, specified as a fraction of the scale factor
measured over the input range. Scale factor asymmetry
implies that the slope of the input-output function is
discontinuous at #ero input. It must !e distinguished from
other nonlinearities.
Sca!e factor stabi!ity, which is the capa!ility of the
inertial sensor to accurately sense angular velocity (or
acceleration) at different angular rates (or at different
accelerations), can also !e used to descri!e scale factor.
Scale factor sta!ility is presumed to mean the variation of
scale factor with temperature and its repeatabi!ity, which is
expressed as part per million (;;*). 4eviations from the
theoretical scale are due to system imperfections.
or example a gyro with ->> ppm scale factor
error (?>.>-<) the angular error is in the area of >., degree
after one revolution. =ith a laser gyro or high performance
fi!er optical gyro system with ,>> ppm scale factor error
the angular error is >.>>+ degrees if the rotation angle is +>
degrees.
(.&.% A''eleration dependent dri3ts (3or gyros)
5yros (mechanical or *:*S) produce a drift
(angular rate offset) dependent on linear acceleration and
environmental vi!ration influence, a so-called g-
dependent drift expressed in ./hr/g.
$nisoinertia errors arise in spinning mass
gyroscopes and introduce !iases owing to ine(ualities in
gyroscope moments of inertia a!out different axes.
'nisoinertia is fre(uency sensitive if the rotor is driven !y
a hysterisis motor. This is a conse(uence of the elastic
coupling !etween the magnetic ring on the rotor and the
rotating magnetic field. The resulting !iases are
proportional to the product of angular rates applied a!out
pairs of orthogonal axes. The anisoinertia coefficients may
!e expressed in units of ./hr/(rad/sec)1. To eliminate such
an error source re(uires complicated compensation, so the
system designer must understand the interplay !etween
vehicle environment, sensor design, and the resulting
navigation system performance when selecting a
particular type of gyroscope for a strap-down system.
(.&.4 A5is (Sensor -ounting) -isalignments
'xes misalignment is the error resulting from the
imperfection of mounting the sensors. It usually results in
a non-orthogonality of the axes defining the INS !ody
frame. 's a result, each axis is affected !y the
measurements of the other two axes in the !ody frame
(sometimes called cross-coup!ing errors). 'xes
misalignment can, in general, !e compensated or modeled
in INS error e(uation.
' gyro misa!ignment, allows !oth roll and pitch
gyros to sense a fraction of the yaw rate, causing an
attitude error a!out level axes. The attitude error, in turn,
causes an acceleration error along the level axes. The
sensor misalignment factor depends not on the (uality of
the sensors !ut on how well a factory cali!ration is
performed with inertial sensors rigidly mounted together.
or example, a misalignment !etween the gyro
axes (or accelerometer axes) causes a cross- coupling
!etween the measurement axes. ' misalignment of >.,
mrad inside of the system (e.g. residual cali!ration
mismatch) leads to a roll error of >.>-@ degree during one
revolution around the yaw axis (if t%e system is unaided).
The smaller the re(uired misalignment the higher the
re(uirements to sensor performance and cali!ration
e(uipment such as three-axes turn-tab!es.
(.&.6 7andom 8al9
This value, given in deg/s(rt(hr) for gyro or
g/s(rt()#) for accelerometer, shows the noise of the used
sensors. 8andom wal" arises due to the integration of white
noise present in the inertial sensors ('ngle random wal"
('8=) for gyros or Aelocity random wal" for
accelerometers). Noise is an additional signal resulting from
the sensor itself or other electronic e(uipment that interfere
with the output signals trying to measure. Noise is in
general non-systematic and therefore cannot !e removed
from the data using deterministic models. It can only !e
modeled !y stoc%astic process.
or gyros, the higher the noise the more noise is
measured on the angular rates and on the angles. Some
manufacturers also specify it as the noise density in
deg/hr/s(rt()#). 3oth values are e(uivalent - if the second
value is divided !y @>, you get it in deg/s(rt(hr). 'n angular
random wal" of >.>>- deg/s(rt(hr) indicates, that the
angular error (uncertainty) due to random wal" is e.g. >.>>-
deg after , hour (unaided) or >.>>+ deg after -> minutes (all
values one sigma). The angular random wal" is very
important for the accuracy of north see"ing (gyro
compassing), !ecause if the random wal" decreases times +
then the needed duration for north see"ing decreases !y
times four (if the resolution of the gyro is high enough) as
per the e(uation (+.-),
(+.-)
where
is the alignment time,
is the earth rotation rate (? ,B.>C deg/hr),
is the latitude at which the alignment ta"es
place,
A is the a#imuth error during alignment.
(.&.# :i+ration dependent errors
There are some inertial sensor errors arises !ecause
of the effects of vi!ratory and oscillatory motion and
includes
'nisoelastic !iases
Scale factor Non-linearity
Ai!ro-pendulous errors
6oning errors
Sculling errors
Si#e effect errors
;seudo motion errors
(.&.#.& Anisoelasti' +iases (3or gyros)
$nisoe!astic biases, so-called g1-dependent drift
expressed in ./hr/g1 are proportional to the product of
acceleration along orthogonal pairs of axes. Such !iases
arise in spinning mass gyroscopes !ecause the gyroscope
rotor suspension structure, particularly the !earings, has
finite compliances which are une(ual in different directions.
(.&.#.( :i+ro-pendulous errors (3or a''elerometers)
4ynamic cross-coupling in pendulous
accelerometers arises owing to angular displacement of the
pendulum which gives rise to a rectified output when
su!Dected to vi!ratory motion. This type of error can arise in
any pendulous accelerometer depending on the phasing
!etween the vi!ration and the pendulum displacement. The
magnitude of the resulting error is maximi#ed when the
vi!ration acts in a plane normal to the pivot axis at CB. to
the sensitive axis and when the pendulum displacement is in
phase with the vi!ration. This error may !e expressed in
units of g/g1.
(.&.#.% !oning errors
6oning is the conical (or near conical) motion in
inertial space of one of the gyroscope input axes, as
+
cos
a
e
$R&
'
$
_
,
a
'
e
+
, cos
s
/ a/
s
t
x B
,
( )
g,
B
( )
>
v
( )
a/
B
(-.B)
Similarly, a heading error gives rise to a
north position error of$
(-.@)
'n effective angular rate !ias acting a!out the
vertical axis gives rise to a north position
error of$
(-.G)
7ver very short periods of navigation, that is,
navigation over a small fraction of a Schuler period, further
simplifications may !e made to these expressions. The
position error contri!utions in the long and short term are
summari#ed in Ta!le -.,. These e(uations can !e used to
assess system performance or to specify an inertial
navigation system to fulfill a particular application. The
propagation of each error type with time in an unaided 0/S
is illustrated in igures -.C % -.H.
In igure -.C % -.H, plots are given illustrating the
propagation of position errors with time over a --hour
period. The system is assumed to !e stationary and located
on the surface of the :arth at a latitude o3 &>?. The
following errors sources are included.
Initial alignment errors - roll/pitch
? >., mrad
Initial heading error ? ,.> mrad
Initial velocity error ? , m/sec
5yro !ias ? >.>, deg/hr
'ccelerometer !ias ? +>> micro-g
*a+le %.&1 ,rowth o3 INS position errors in the long term
and short term
igure -.C shows the north position error resulting
from the hori#ontal component of gyro !ias. It is seen that
given a residual gyro !ias of >.>, deg/hr, the position error
will grow at the rate of >.G nm/hr. It should !e noted that
this represents the error !uild up along the North and :ast
axes, so that the radial error will !e multiplied !y s=rt(() or
s=rt(north@( A east@() and hence will !e e(ual to & nmBhr.
Thus in order to achieve * nm1%r position accuracy, one
should go for >.>, deg/hr class gyros.
2igure %.41 North position error verses time 'aused +y hori;ontal
gyro +ias
>
sin
s
/ g,
s
t
x B R t
_
,
( )
>
> >
sin
cos
s
/
s
t
x R t
_
,
&
( )
g.
B
+
> +
, cos
cos
+
s
/ g.
s
t t
x B R
_
' ;
,
&
( )
> >
,
( )
>
( )
> /
v
( )
g
B
( )
a
B
igure -.B shows a plot on the north position error
resulting from an initial velocity error. This error is !ounded
!y the effect of the Schuler loop and a figure for the
maximum position error caused !y a , m/sec initial velocity
error.
2igure %.61 North position error verses time 'aused +y initial
velo'ity error
igure -.@ shows a plot of the north position error
resulting from the hori#ontal component of accelerometer
!ias. The propagation of this type of error is also !ounded
!y the effect of the Schuler loop. ' residual accelerometer
!ias of +>> micro-g (corresponding to an initial tilt error of
>.+ mrad) will result in pea" position error of ,.-F nautical
miles.
2igure %.#1 North position error verses time 'aused +y
a''elerometer +ias
igure -.G shows a plot of the north position error
resulting from an initial error in a#imuth alignment. The
resulting position error comprises a Schuler oscillation as a
ramp function, the magnitude of the error varying with
system latitude and speed. Some examples of the resulting
error magnitudes are given for an a#imuth misalignment of
, mrad.
2igure %.>1 North position error verses time 'aused +y initial
a;imuth misalignment
igure -.F shows a plot of the north position error
resulting from an initial misalignment error with respect to
vertical (initial roll or pitch error). The propagation of this
type of error is also !ounded !y the effect of the Schuler
loop. 'n initial roll or pitch error of >., mrad (resulting
from ,>> micro-g accelerometer !ias) will result in pea"
position error of around >.G nautical miles.
2igure %.C1 North position error verses time 'aused +y initial
misalignment w.r.t verti'al
igure -.H shows a plot of north position error
resulting from a vertical gyroscope !ias component. The
position error comprises a (uadratic term with a
superimposed Schuler oscillation& the relative magnitude of
the oscillatory component is small and hence is not readily
apparent in the figure with the scaling adopted here. The
propagation of this error also varies with latitude and
example figures are given in the figure for a >.>,./h
gyroscope !ias.
2igure %.D1 North position error verses time 'aused +y verti'al
gyro +ias 'omponent
's with the north channel, the growth of many of
the errors is !ounded !y the effects of the Schuler tuning.
)owever, this is not the case in the vertical channel where
the errors increase rapidly with time. or example, a
vertical accelerometer !ias, B2az, will give rise to a
position error of B2az3t114# It is for this reason that aircraft
navigation systems commonly operate in conDunction with a
barometric or radar altimeter in order to restrict the growth
of vertical channel errors.
2igure %.&$1 Air'ra3t unaided INS single 'hannel error
propagation
The errors are summari#ed in igure -.,> for an
aircraft navigation system over a - hour period. The figure
illustrates the form of the position errors resulting from the
various error sources whilst the upper curve (dotted !lue
line) represents the com!ined effect of the individual errors.
The upper curve has !een o!tained !y summing the
individual error components (uadratically to give the total
navigation error. It will !e seen that the gyroscope !ias and
a#imuthal misalignment contri!utions grow with time
whilst the other terms are !ounded as a result of Schuler
tuning. It should !e noted that these are simplified results
for a single channel. In a full model, the Foucau!t effect is
noticea!le even at the first Schuler period.
The oucault oscillation, given !y e = sin .
This maintains itself as a modulation of the Schuler
oscillation, the modulation in the two hori#ontal channels
!eing H>. apart in phase. The oucault oscillation has a
period of 4 1 sin . The period of this oscillation is
a!out -> hours for moderate latitudes.
%.% Error propagation during air'ra3t maneuvering
The analysis of inertial navigation system
performance for air!orne applications rapidly !ecomes
complex when account is ta"en of realistic vehicle
traDectories and maneuvers, in which case the analyst will
usually turn to simulation to aid the design process.
In addition to the usual effects of alignment errors
and sensor !iases illustrated in the a!ove error plots, a
num!er of error contri!utions arise as a result of the
acceleration experienced !y an aircraft during ta"e-off.
3oth alignment errors and mounting misalignments of the
accelerometers will give rise to cross-trac5 ve!ocity errors
as the aircraft accelerates during ta"e-off. or example, an
a#imuthal misalignment of will give rise to a velocity
error of , where 6 is the cruise speed of the aircraft.
The velocity error which has !uilt up during ta"e-off then
propagates during the su!se(uent cruise phase of flight in
the same way as an initial velocity error. Sensor scale-factor
inaccuracy and acceleration dependent !iases will give rise
to additional navigation errors in the event of an aircraft
maneuver. 3ecause of applied accelerations, angular rates,
induced vi!rations during maneuvering degrades the
performance of a S4INS.
or example, an aircraft at vertically down
position, the accelerometer will measure , g at that time
(3ut in real case, it was not always at ,g). 't that time with
accelerometer scale factor error of +>> ppm will
corresponds to +>> micro-g. 3ecause of this scale factor
error, the position error also propagates and is bounded
within Schuler oscillation. =hen using an optical gyros, a
coning motion of >.,I at a fre(uency of B> )#, for instance,
can result in a drift in the computed attitude of almost ,>>I/
hr.
3ecause of Sculling error, a steady acceleration
component of $.6A.'os occurs in the x-direction. It is
stressed that this error term can arise even when using
perfect accelerometers, !eing purely a function of the
inaccuracy in the resolution process. If, for example, ' ?
,>g, ? >.,. and the phase difference () is #ero, the
resulting acceleration !ias is around H milli-g.
Si#e effect also introduces errors during the
navigation. or example, in the presence of cyclic angular
motion of amplitude >.,., and of fre(uency B> )#, for
example, a !ias of J ,.B milli-g arises for a displacement of
,> cm, i.e. the distance !etween the center of gravity of the
vehicle and the location of the sensor.
4 7e3eren'es
0,2 4avid ). Titterton K Lohn 9. =eston, Strapdown
Inertial Navigation Technology, +
nd
:dition, 0,, Radar,
Sonar, and /avigation Series, +>>C.
0+2 4r. Ing. :. )inue!er, MIf you investigate in an inertial
measurement system... which technical data you should
analyse and compare !efore ma"ing your decisionN, iM$R
'ec%nica! /otes.
0-2 '. 4. Oing, MInertial Navigation % orty Pears of
:volutionN, G,7 Review, Aol.,-, No. -, ,HHF.
0C2 )aiying )ou, M*odeling Inertial Sensor :rrors
Qsing 'llan AarianceN, 87G, Report R +>>,, 4epartment
of 5eomatics :ngineering, Qniversity of 6algary, +>>C.
0B2 6hing-ang 9in, *odern Navigation, 5uidance, and
6ontrol ;rocessing, -rentice 9a!! Series, ,HH,.
0@2 Sergei '. Lere!ets, M5yro :valuation for the *ission
to LupiterN, 0,,, $erospace 7onference, +>>G.
0G2 8. ;. 5. 6ollinson, Introduction to 'vionics, irst
:dition, 7%apman : 9a!! pub!ications, ,HH@.
0F2 4r. Poung 6. 9ee K 4aniel 5. 7E9aughlin, M'
;erformance 'nalysis of a Tightly 6oupled 5;S/Inertial
System for Two Integrity *onitoring *ethodsN, 7enter for
$dvanced $viation System .eve!opment (6''S4)
Technical ;aper, *IT8: 6orporation, ,HHH.
0H2 MInertial Navigation SystemN, &i!;am F!ig%t
'raining 7enter Technical Notes.
>
6
>