3.5.16 · D5Guidance, Navigation & Control (GNC)
Question bank — Mechanization equations — integrating IMU to get position, velocity, attitude
This bank hunts the five classic misconceptions: (1) accelerometers "feel" gravity, (2) rotations add like numbers, (3) Coriolis is a real force, (4) the flat-Earth core is always safe, (5) drift is linear in time.
True or false — justify
A stationary, level accelerometer reads zero on its Down axis.
False. It reads (in NED, Down , the reading is opposite the support). The table pushes it up, and specific force is the contact force, so a static reading is magnitude opposite the support — zero happens only in free fall.
An IMU in free fall (dropped) measures its own falling acceleration.
False. It reads . Gravity accelerates every part of the proof mass equally, so nothing squeezes the sensor — there is no specific force to detect during free fall.
Gyro bias causes velocity error that grows linearly with time.
False. Step it out: a gyro bias tilts attitude linearly, . A tilt of misrotates gravity so a horizontal acceleration leaks into the level channel. Integrating that once gives velocity and again gives position — so position error from a gyro bias grows like , worse than the from an accelerometer bias.
You can compute velocity before attitude if you integrate the accelerometer first.
False. The accelerometer output lives in body axes. You need (attitude) to rotate it into nav axes () before you can add gravity, so attitude must come first.
The gravity vector added in the velocity equation is pure Newtonian gravitation.
False. It is plumb-bob (apparent) gravity = gravitational attraction minus centrifugal acceleration from Earth's spin. That is why it points along the local vertical (Down), not exactly at Earth's centre.
Over a short indoor run you may safely drop the Earth-rate and transport-rate terms.
True. For short baselines the and corrections are negligible; the flat-Earth trio (in NED) is the 80/20 core. The NED convention fixes every sign in that trio.
Rotation rates from the gyro can simply be added over time to get the total rotation angle.
False. Finite 3-D rotations do not commute; orientation lives on the curved manifold . That is exactly why we propagate a matrix/quaternion (, using the skew matrix defined above) instead of summing angles.
carries a minus sign because of a computational convention, not physics.
True (it is a convention). In NED the "Down" axis points toward Earth, so positive downward velocity decreases height; the sign is purely from the choice of axis direction. (In ENU it would be .)
Spot the error
"The accelerometer measures acceleration, so I double-integrate it directly to get position."
The device measures specific force , not acceleration. Integrating raw while static injects a fake ; you must rotate to nav axes and use the full so gravity is added back first.
"Since gyros read , I plug that straight into ."
You need the rate relative to the nav frame, , and then form its skew matrix . Using raw makes attitude rotate with the Earth and the moving frame.
"The Coriolis term is a real physical force acting on the vehicle."
It is a fictitious (frame) correction, not a real force. It appears only because we express Newton's law in a rotating/curving nav frame; there is no physical push behind it.
"."
Missing the : correct is , where is the East velocity component and the transverse radius. Longitude circles shrink toward the poles, so the same eastward speed sweeps more longitude at high latitude.
"Both position rates use the same Earth radius."
They use different radii of curvature: latitude rate uses the meridian radius , longitude rate uses the transverse radius . The ellipsoid curves differently north–south vs east–west (more sharply north–south), so — see Radii of Curvature of the Earth Ellipsoid.
"Velocity error and position error are the same size after a while."
No — velocity is one integration of acceleration error, position is two. Position error accumulates faster ( from accel bias) than velocity error ().
"I subtract gravity from the accelerometer reading to remove it."
You add : the equation is , so the term is where gravity re-enters. The accelerometer already missed true gravitational acceleration (it can't feel it), so gravity is put back in software rather than removed.
Why questions
Why must be rotated into the nav frame before gravity is applied?
Gravity is naturally expressed along the local vertical (a nav-frame direction). To add two vectors they must share a frame, so the body-frame force is rotated by first ().
Why does the velocity equation contain but the position and attitude equations do not?
Coriolis arises from differentiating a velocity observed in a rotating frame — it is inherently a first-derivative-of-velocity effect. Attitude and position kinematics don't involve that cross term the same way (see Coriolis and Centrifugal Effects — Coriolis is , largest when velocity is perpendicular to the rotation axis).
Why is a tilt error the dominant INS error path rather than a raw accelerometer bias?
A small attitude tilt leaks the full into the horizontal channel as . Since , even a tilt injects a larger horizontal acceleration than most accelerometer biases.
Why do we prefer a quaternion or DCM over roll-pitch-yaw angles for propagation?
Euler angles suffer gimbal lock (singularities) and their rates blow up near pitch. The DCM/quaternion stays well-behaved everywhere — but you must give it a correct initial condition (an initial alignment or unit quaternion) and, because numerical integration drifts off the constraint, renormalize each step (orthonormalize the DCM, or scale the quaternion back to unit norm). See Direction Cosine Matrix and Quaternions.
Why does the transport rate exist at all?
Because the local "down" direction physically re-tilts as you travel over a curved Earth. The nav frame rotates to keep pointing down, and that rotation must be removed from the gyro measurement.
Why is a strapdown INS more demanding on the mechanization software than a gimbaled one?
In strapdown the sensors are bolted to the body, so all the frame rotations (attitude layer) must be done numerically at high rate; a gimbaled platform mechanically holds the frame, doing the rotation in hardware (see Strapdown vs Gimbaled INS).
Why does dividing arc-rate by radius give an angular rate?
Arc length , so and . Northward speed is the meridian arc rate, so .
Why can the Kalman filter estimate the sensor biases that cause drift?
External aids (like GNSS) reveal the growing INS error, and because bias-driven error has a known time signature (, , ), the filter can attribute drift to the underlying bias and correct it (see Kalman Filter for INS-GNSS Integration — the filter blends a drifting-but-smooth INS with a noisy-but-bounded GNSS fix).
Edge cases
At the geographic pole, what breaks in the position rates, and what stays fine?
The longitude rate diverges because — longitude is ill-defined at the pole. The latitude rate stays perfectly finite (the meridian radius is a normal positive length there), so only the east-west channel is singular; wander-azimuth or a different frame fixes it.
For a perfectly stationary vehicle, what should be, and why isn't it just ?
It should be . The upward support force ( Down in ) rotates into nav axes and exactly cancels the added Down gravity , so the two terms in sum to zero.
What happens to attitude propagation when the gyro reads exactly zero?
, so orientation is held constant — correct for a non-rotating body, and it also means any zero-rate bias still silently tilts the frame.
Is the flat-Earth model safe for a slow-moving vehicle that runs for hours?
No. Even at low speed, Earth rate () and transport effects accumulate over long time, so long-duration navigation needs the full Earth terms regardless of speed.
What is the limiting behaviour of position error from a constant accelerometer bias over time?
It grows as without bound — the unaided INS drift signature. This is why bounded external aiding (GNSS) is required for long missions.
Does an accelerometer at rest on the equator read exactly gravitational attraction?
No — it reads apparent gravity, which is gravitational attraction minus the equatorial centrifugal term ( smaller). The mechanization uses this apparent value (see Gravity Model and Geoid).
Recall One-line summary of every trap here
The accelerometer feels specific force not gravity; rotations don't add (propagate with a correct initial alignment and renormalization); Coriolis is a frame term not a force; the flat-Earth core (NED) fails over long time not just distance; and drift grows like (accel bias) or (gyro tilt) — never linearly.