3.5.4 · D5Guidance, Navigation & Control (GNC)
Question bank — DCM kinematics — Ċ = −[ω×]C
Recall Quick symbol refresher before you start
- = the Direction Cosine Matrix, maps a vector's components from inertial frame to body frame : .
- = the skew-symmetric matrix (cross-product operator) built from vector , so that . See Skew-symmetric matrices & cross-product operator.
- = angular velocity of the body, written in body coordinates.
- "Orthonormal" = rows are unit-length and mutually perpendicular, i.e. .
True or false — justify
Each answer must give the reason, not a bare verdict.
says the rate of change of is proportional to itself.
True — the equation is linear in : the derivative is left-multiplied by a fixed (at that instant) skew matrix. That is why the solution is a matrix exponential, not a polynomial.
The matrix is itself a rotation matrix.
False — lives in the tangent space to SO(3), not on the group. It is (skew) (rotation), which is generally neither orthogonal nor of determinant 1.
If then is constant in time.
True — zero spin gives , so freezes. Physically the body is not turning, so the map between frames cannot change.
The product is skew-symmetric at every instant.
True — it follows directly from differentiating the constraint , which holds for all . That single fact is what forces the whole skew structure.
is an assumption we impose at and hope survives.
False for the exact ODE — the equation guarantees , so orthonormality is conserved automatically once true initially. It is only numerical integration that must re-impose it.
Swapping to track removes the minus sign.
True — for the inverse map you get . The sign is a bookkeeping consequence of which direction the matrix maps, not a physical fact.
The three independent numbers in the skew matrix are exactly the three components of .
True — a skew matrix has zeros on the diagonal and one sign-flipped copy of each off-diagonal, leaving 3 free entries. These map one-to-one to via the vee-map.
Because is linear in , you can integrate it like a scalar: always.
False in general — the clean exponential holds only when is constant. If varies in time and the skew matrices at different times don't commute, you need a time-ordered exponential, not .
Spot the error
Find the flaw in each stated claim.
"I'll plug the inertial angular velocity into ."
Wrong frame — with the skew must use body-frame components , so the output lands consistently in the body frame. Using mismatches the frames on the two sides.
" has a plus, so ."
The plus formula describes a vector seen from the inertial frame; maps inertial into body, so from the body's viewpoint inertial vectors appear to spin backward — hence the minus. The two formulas describe opposite viewpoints.
"Euler-integrating keeps orthogonal since the ODE preserves it."
The continuous ODE preserves exactly; a discrete Euler step takes a straight line off the curved surface SO(3), so drifts. You must renormalize (Gram–Schmidt / SVD) or propagate a quaternion.
"The vee-map gives ."
Sign error — since , the correct extraction carries a minus: . Dropping it flips the recovered spin direction.
" works for any invertible ."
Only for proper rotations (). A reflection () flips handedness and introduces an extra minus sign, breaking the identity used in the sign derivation.
" has an inverse, so I can always solve for ."
A skew matrix is singular (odd-dimensional skew matrices have determinant 0), so it has no inverse. Recovery uses the transpose trick , exploiting , not any inverse of the skew.
"Since maps N→B, its rows are the inertial basis vectors in body coordinates."
Backwards — the rows of are the body-frame basis vectors expressed in inertial coordinates. The entry makes this explicit.
Why questions
Why must the derivative have exactly a skew-symmetric-times- form and not some arbitrary structure?
Because is locked onto the curved surface ; the only velocities that keep it on that surface point along the surface's tangent directions, which are precisely {skew matrix} . Any other form would drag off the rotation group.
Why does angular velocity — a 3-component vector — control the change of a 9-component matrix without over- or under-determining it?
A rotation matrix, despite 9 entries, has only 3 true degrees of freedom (orthonormality removes 6). So 3 numbers in exactly match the 3 ways can change — a perfect count from the Lie-algebra dimension, see Rotation group SO(3) and Lie algebra so(3).
Why is the minus sign a convention and not a physics choice you could get wrong experimentally?
The physical spin is fixed; the sign only records whether you wrote the map as (minus) or (plus). Both predict identical physical motion — you just must declare which convention you use.
Why is called the "parent" of quaternion and Euler-angle kinematics?
Quaternions and Euler angles are just different coordinates for the same rotation; their rate equations (Quaternion kinematics — q̇ = ½ Ω(ω) q, Euler angle kinematics & gimbal lock) are re-parametrizations of this single DCM equation. Same physics, different bookkeeping.
Why does this equation appear under the name "Poisson's equation" in some texts?
Because Poisson's equation for rotating frames is the general statement that a body-fixed quantity's inertial rate equals ; applied to each column/basis vector of the DCM it collapses to exactly .
Why can't you extract from a single snapshot of alone?
is a rate, encoded in how is changing, not in its instantaneous value. You need (or two nearby samples of ) to form and read the spin off.
Edge cases
What is at the instant , and does have to equal then?
: the map is momentarily frozen. But itself can be any rotation — a non-spinning body can still be pointing in any fixed orientation.
If is a pure spin about a body axis, what special property does the resulting have?
It is a planar rotation about that fixed axis, given by the matrix exponential . The spin axis is an eigenvector with eigenvalue 1, so that axis's direction is preserved.
What happens to the solution if points along a time-varying axis whose direction keeps shifting?
The simple formula fails because skew matrices at different times don't commute; you must integrate the ODE step by step (or use a time-ordered product). Physically this is why real tumbling attitude has no closed form.
At the "gimbal lock" configuration of Euler angles, does the DCM equation also break down?
No — the DCM equation is perfectly well-behaved everywhere; only the Euler-angle re-parametrization becomes singular (rates blow up), see Euler angle kinematics & gimbal lock. This robustness is a key reason to propagate the DCM or a quaternion instead.
For an angular velocity of enormous magnitude, does anything in cap or saturate?
Nothing in the exact equation saturates — simply rotates faster, staying on SO(3). The practical limit is numerical: large times a finite step size causes big drift, forcing smaller steps or renormalization.
If you feed in a that is accidentally not orthonormal, will the equation fix it over time?
No — the equation only conserves ; it does not restore it. If , that error is carried forward unchanged (its rate of change is zero, but so is its correction). Start clean.
Can ever flip to under this equation?
No — the flow stays continuously on the connected component with . Since is continuous and starts at , and the ODE never leaves , it cannot jump to a reflection.