3.5.4 · D3 · Physics › Guidance, Navigation & Control (GNC) › [[3.5.04 DCM kinematics — Ċ = −[ω×]C|DCM kinematics — Ċ = −[ω×]C]]
Yeh page ek drill hai. [[3.5.04 DCM kinematics — Ċ = −[ω×]C (index 3.5.4)|Parent note]] ne yeh equation prove ki thi:
C ˙ = − [ ω B × ] C .
Yahan hum isme har tarah ka input daalenge — har spin axis, ek general tilted axis, zero spin, sign-changing spin, ek limiting long-time case, time-varying spin, ek real spacecraft word problem, aur ek exam twist. Agar tum yeh sab kar sako, toh exam mein kuch bhi surprise nahi karega.
Shuru karne se pehle, ek reminder plain words mein. Symbol [ w × ] ("skew matrix of w ") ek 3 × 3 box of numbers hai jo vector w = ( w 1 , w 2 , w 3 ) se banta hai, aur ise kisi bhi vector a se multiply karne par cross product w × a milta hai. Dekho Skew-symmetric matrices & cross-product operator . Likha hua:
[ w × ] = 0 w 3 − w 2 − w 3 0 w 1 w 2 − w 1 0 .
Dhyaan do ki diagonal hamesha zero hai aur ise diagonal ke paas flip karne par har sign flip ho jaata hai — yahi "skew-symmetric" ka matlab hai, aur yehi wajah hai ki C ek rotation rehta hai.
Definition Is poori page par use hone waala Notation
Hum bracket ke andar vector par arrow rakhte hain: [ ω B × ] ka matlab hamesha "vector ω B se bana skew matrix" hai. Jab space tight ho toh hum [ ω × ] likh sakte hain, lekin bracket ke andar object hamesha vector hota hai, koi scalar nahi.
Row/column indices 1 , 2 , 3 se chalte hain (one-based), jo subscripts b ^ 1 , b ^ 2 , b ^ 3 se match karte hain. "Index 1" = pehli row/column, yaani x -row.
Vee-map ∨ (formally Ex 7 mein introduce kiya gaya) skew bracket ka inverse hai: diya hua skew matrix S , S ∨ woh vector w padhta hai jis se [ w × ] = S banta hai.
Is topic ke har problem ka jawab in mein se kisi ek cell mein aata hai. Neeche ke worked examples mein cell ka tag laga hai.
#
Case class
Kya khaas hai
Example
A
Kisi ek body axis (x , y , ya z ) ke baare mein Spin
ek 2-D rotation ban jaata hai
Ex 1, Ex 2
B
Zero spin ω = 0
degenerate: C ˙ = 0 , kuch nahi hilta
Ex 3
C
Sign flip / reverse spin (− ω )
rotation ki direction ulti ho jaati hai
Ex 4
D
Limiting / long-time behaviour (t → ∞ , periodicity)
solution periodic hai, kabhi blow up nahi hoti
Ex 5
E
General (tilted) constant axis (Rodrigues / matrix exponential)
full 3-component ω , ek integration
Ex 6
F
Extract ω from data (vee-map)
inverse problem
Ex 7
G
Time-varying ω ( t )
simple exponential nahi; dhyan se integrate karo
Ex 8
H
Real-world word problem (spacecraft rate)
physical scenario se [ ω × ] banao
Ex 9
I
Exam twist : wrong-sign trap / frame trap
plus/minus aur N-vs-B error pakdo
Ex 10
J
Degenerate check : kya integration det C = + 1 preserve karti hai?
orthonormality + handedness
Ex 11
Worked example Ex 1 — Body
z ke baare mein Constant spin, C ( t ) nikalo
Ek rigid body inertial frame ke saath aligned shuru hota hai (C ( 0 ) = I ) aur apne body z -axis ke baare mein constant rate ω = 2 rad/s se ghoomta hai, toh ω B = ( 0 , 0 , 2 ) . C ( t ) nikalo aur ise t = 4 π s par evaluate karo.
Forecast: pehle guess karo — t = 4 π ≈ 0.79 s par body ω t = 2 π (ek quarter turn) ghoom chuki hogi. Kya top-left 2 × 2 block + ya − rotation jaisa dikhega? Equation mein minus sign yaad rakho.
Skew matrix banao. ω B = ( 0 , 0 , 2 ) ke saath,
[ ω × ] = 0 2 0 − 2 0 0 0 0 0 .
Yeh step kyun? Har DCM-kinematics problem ω ko uske skew box mein convert karke shuru hoti hai; woh box hi woh operator hai jo C ko "twist" karta hai.
ODE type pehchano. Kyunki ω constant hai, C ˙ = − [ ω × ] C ek linear ODE with constant coefficients hai. Iska solution matrix exponential hai:
C ( t ) = e − [ ω × ] t C ( 0 ) .
Yeh step kyun? Jaise numbers ke liye x ˙ = a x ⇒ x = e a t x 0 , matrix version bhi kaam karta hai kyunki [ ω × ] kabhi nahi badlta.
Exponential evaluate karo. Ek z -axis skew matrix ka exponential ek rotation hota hai. Minus sign ke saath yeh z ke baare mein − ω t angle ki rotation hai:
C ( t ) = cos ω t − sin ω t 0 sin ω t cos ω t 0 0 0 1 .
Yeh step kyun? Minus + ω t ko clockwise turn mein badal deta hai — Lie-algebra picture dekho jahan skew matrix I ke tangent space mein rehta hai.
t = π /4 , ω = 2 plug karo: ω t = π /2 , toh cos = 0 , sin = 1 :
C ( 4 π ) = 0 − 1 0 1 0 0 0 0 1 .
Verify: Check karo C C T = I (hota hai, columns orthonormal hain), aur det C = + 1 . t = 0 par, ω t = 0 deta hai C = I ✔.
Figure 1 (neeche): white arrows fixed inertial axes n ^ 1 , n ^ 2 hain; blue aur yellow arrows quarter turn ke baad body axes b ^ 1 , b ^ 2 hain. Red arc dikhata hai rotation clockwise (− ω t ) hai, aur yehi wajah hai ki − sin ω t C ke lower-left mein hai — equation mein minus sign apne aap clockwise sweep ki tarah dikh raha hai.
Intuition Figure 1 caption
Body z -spin at ω t = π /2 : inertial axes (white) fixed hain, body axes (blue/yellow) clockwise − ω t se sweep hue hain (red arc). Clockwise sweep C ˙ = − [ ω × ] C mein minus sign ka geometric chehra hai.
Worked example Ex 2 — Body
x ke baare mein Spin: kaun si entries hilti hain?
Same setup (C ( 0 ) = I ) lekin ω B = ( 3 , 0 , 0 ) rad/s. C ( t ) likho aur t = 6 π s par C do.
Forecast: kaun sa 2 × 2 block frozen rehega, aur kaun sa rotate karega?
Skew matrix: [ ω × ] = 0 0 0 0 0 3 0 − 3 0 .
Yeh step kyun? Sirf w 1 = 3 nonzero hai, toh sirf w 1 entries survive karti hain — x -axis untouched rehta hai.
x -axis skew ka exponential (minus ke saath) x ke baare mein − ω t ki rotation hai:
C ( t ) = 1 0 0 0 cos ω t − sin ω t 0 sin ω t cos ω t .
Yeh step kyun? Pehli row aur pehla column (index 1, yaani hamare one-based convention mein x -row/column) identity ke roop mein rehte hain; y –z block turning karta hai.
ω t = 3 ⋅ π /6 = π /2 plug karo: cos = 0 , sin = 1 :
C = 1 0 0 0 0 − 1 0 1 0 .
Verify: det C = + 1 , C C T = I , aur pehli row/column ( 1 , 0 , 0 ) unhili hai — x -spin ke liye sahi ✔.
Worked example Ex 3 — Bilkul bhi rotation nahi
ω B = ( 0 , 0 , 0 ) kisi starting C ( 0 ) = C 0 ke saath. C ( t ) kya hai?
Forecast: trivial — lekin paka hona chahiye ki kyun trivial hai.
Skew matrix zero matrix hai: [ ω × ] = 0 .
Yeh step kyun? Box ki har entry ω ke ek component se banti hai, aur sab zero hain.
ODE ban jaati hai C ˙ = 0 , toh C ( t ) = C 0 har t ke liye.
Yeh step kyun? Zero derivative ka matlab hai attitude kabhi nahi badlta — physically body spin nahi kar raha.
Exponential ke saath consistency: e − 0 ⋅ t = I , toh C ( t ) = I ⋅ C 0 = C 0 ✔.
Verify: C ( t ) = C 0 constant hai, iska rate zero hai, orthonormality trivially rakhi gayi hai. Yeh poore topic ka degenerate corner hai.
Worked example Ex 4 — Spin direction palto
Ex 1 lo lekin ω B = ( 0 , 0 , − 2 ) ke saath (dusri taraf ghoomna). C ( t ) nikalo aur Ex 1 se t = π /4 par compare karo.
Forecast: kya Ex 1 ke compare mein off-diagonal signs swap honge?
Skew matrix sign flip karta hai: [ ω × ] = 0 − 2 0 2 0 0 0 0 0 .
Yeh step kyun? [( − ω ) × ] = − [ ω × ] — skew map linear hai.
Exponential z ke baare mein − ω t = + 2 t ki rotation deta hai:
C ( t ) = cos ω t sin ω t 0 − sin ω t cos ω t 0 0 0 1 ,
cos ( − x ) = cos x , sin ( − x ) = − sin x use karke.
Yeh step kyun? ω reverse karna turn ki handedness reverse karta hai; sirf sin terms ka sign badlta hai.
t = π /4 , ω t = π /2 par:
C = 0 1 0 − 1 0 0 0 0 1 .
Verify: Yeh bilkul Ex 1 ke answer ka transpose hai — spin reverse karna inverse rotation deta hai ✔. C Ex4 = C Ex1 T , aur C Ex1 C Ex4 = I .
Worked example Ex 5 — Kya
C kabhi blow up karta hai? Periodicity.
Constant ω B = ( 0 , 0 , ω ) , ω = 2 rad/s ke saath, t → ∞ hone par C ( t ) describe karo. Sabse chhota T > 0 nikalo jahan C ( T ) = C ( 0 ) .
Forecast: unstable growth, ya hamesha bounded?
Entries hain cos ω t , sin ω t — hamesha − 1 aur + 1 ke beech bounded.
Yeh step kyun? Real a ke liye e a t ke unlike, skew matrix ka exponential purely imaginary eigenvalues rakhta hai, toh yeh barhne ki jagah oscillate karta hai. Yeh S O ( 3 ) ki pahchaan hai: rotations kabhi blow up nahi hote.
Period: C tab repeat karta hai jab ω t mein 2 π ka izaafa hota hai, yaani T = ω 2 π = 2 2 π = π s.
Yeh step kyun? Ek poora revolution frame ko apni shuru wali jagah wapas le aata hai.
Limit: t → ∞ ke saath koi ek limit nahi hai; C ( t ) π s period se cycle karta hai.
Verify: ω T = 2 ⋅ π = 2 π , toh cos = 1 , sin = 0 , deta hai C ( T ) = I = C ( 0 ) ✔. Bounded aur periodic — "stable corner."
Worked example Ex 6 — Ek tilted body axis ke baare mein Spin
Ab axis ek coordinate axis nahi hai. Maano constant body rate ω B = ω e ^ hai jahan unit axis e ^ = 2 1 ( 1 , 1 , 0 ) aur ω = 2 π rad/s hai. C ( 0 ) = I ke saath, C ( t ) closed form mein nikalo aur t = 1 s par evaluate karo.
Forecast: tilted axis ke saath tum sirf ek 2 × 2 block nahi bhar sakte — lekin phir bhi ek clean formula hai. Guess karo kaun se trig functions aate hain.
Angular velocity ko magnitude aur axis mein split karo: ω B = ω e ^ , ∣ e ^ ∣ = 1 . Iska skew hai [ ω × ] = ω [ e ^ × ] .
Yeh step kyun? Scalar rate ko aage nikaalne se hum unit-axis skew [ e ^ × ] use kar sakte hain, jinki powers simple hain: [ e ^ × ] 3 = − [ e ^ × ] .
Rodrigues' formula se exponential likho. Kyunki C ( t ) = e − [ ω × ] t = e − ϕ [ e ^ × ] hai jahan turn angle ϕ = ω t hai, matrix-exponential identity deta hai:
C ( t ) = I − sin ϕ [ e ^ × ] + ( 1 − cos ϕ ) [ e ^ × ] 2 , ϕ = ω t .
Yeh step kyun? C ˙ = − [ ω × ] C mein minus sign − sin ϕ term ban jaata hai — e ^ ke baare mein clockwise turn. Kisi bhi constant axis ke liye yahi general answer hai; single-axis Ex 1/Ex 2 special cases hain.
[ e ^ × ] aur [ e ^ × ] 2 banao. e ^ = 2 1 ( 1 , 1 , 0 ) ke saath:
[ e ^ × ] = 2 1 0 0 − 1 0 0 1 1 − 1 0 , [ e ^ × ] 2 = 2 1 − 1 1 0 1 − 1 0 0 0 − 2 .
Yeh step kyun? Sirf yeh do matrices Rodrigues ko chahiye hain — koi infinite series nahi.
ϕ = ω t = 2 π ⋅ 1 = 2 π plug karo: sin ϕ = 1 , 1 − cos ϕ = 1 :
C ( 1 ) = I − [ e ^ × ] + [ e ^ × ] 2 = 2 1 2 1 2 1 2 1 2 1 − 2 1 − 2 1 2 1 0 .
Verify: C C T = I aur det C = + 1 (VERIFY mein check kiya gaya). Saath hi axis e ^ unhila hai: C e ^ = e ^ , "e ^ ke baare mein rotation" ki fingerprint ✔.
∨
Skew bracket [ w × ] ek vector ko matrix mein badalta hai. Vee-map ∨ ulta karta hai — yeh ek skew matrix ko wapas uske vector mein badalta hai. Agar
S = 0 s 3 − s 2 − s 3 0 s 1 s 2 − s 1 0 , toh S ∨ = s 1 s 2 s 3 .
Words mein: s 1 ko ( 3 , 2 ) entry se padho, s 2 ko ( 1 , 3 ) entry se, s 3 ko ( 2 , 1 ) entry se. Yeh exactly skew bracket ka inverse hai: ( [ w × ] ) ∨ = w .
Worked example Ex 7 — Measured
C ˙ C T se ω padho
Ek sensor tumhe product deta hai:
C ˙ C T = 0 − 5 4 5 0 − 3 − 4 3 0 .
ω B nikalo.
Forecast: teen independent off-diagonal numbers hi ω hain signs tak — compute karne se pehle guess karo.
Inverse relation use karo. Parent se, [ ω B × ] = − C ˙ C T :
[ ω × ] = 0 5 − 4 − 5 0 3 4 − 3 0 .
Yeh step kyun? C ˙ = − [ ω × ] C ko right mein C T se multiply karke aur C C T = I use karke skew matrix isolate ho jaata hai.
Vee-map apply karo ∨ (abhi upar define kiya gaya) skew template se components padhne ke liye:
ω 1 = ( ⋅ ) 32 = 3 , ω 2 = ( ⋅ ) 13 = 4 , ω 3 = ( ⋅ ) 21 = 5.
Yeh step kyun? Standard [ w × ] template se compare karo: ( 3 , 2 ) entry w 1 hai, ( 1 , 3 ) entry w 2 hai, ( 2 , 1 ) entry w 3 hai.
Result: ω B = ( 3 , 4 , 5 ) rad/s.
Verify: ( 3 , 4 , 5 ) se [ ω × ] dobara banao — yeh step 1 se match karta hai. Magnitude ∣ ω ∣ = 9 + 16 + 25 = 50 ≈ 7.071 rad/s ✔. Yeh attitude determination ko mirror karta hai jahan tum data se rates recover karte ho.
Worked example Ex 8 — Fixed body axis ke baare mein Spin-up,
ω ( t ) ramp karta hai
Body apne z -axis ke baare mein ghoomti hai lekin rate linearly ramp karta hai : ω B ( t ) = ( 0 , 0 , α t ) jahan α = 2 rad/s 2 , aur C ( 0 ) = I . C ( t ) nikalo aur t = 1 s par evaluate karo.
Forecast: shayad e − [ ω × ] t likhne ka mann kare — lekin ω ab constant nahi hai. Kya naive exponential yahan bhi kaam karta hai?
Check karo ki [ ω ( t ) × ] alag-alag times par khud se commute karta hai ya nahi. Yahan saare skew matrices [ ω ( t ) × ] ek hi same axis z ke baare mein hain, bas α t se scale kiye gaye hain. Matrices jo ek fixed matrix ke scalar multiples hain woh hamesha commute karte hain.
Yeh step kyun? C ˙ = − [ ω ( t ) × ] C ka general solution ek time-ordered exponential ki zaroorat hai; yeh ordinary exponential of the integral mein tab collapse hota hai jab alag-alag times par skew matrices commute karein. Yahan karte hain (single fixed axis), toh hum angle integrate kar sakte hain.
Rate ko turn angle mein integrate karo. Define karo ϕ ( t ) = ∫ 0 t ω ( s ) d s = ∫ 0 t α s d s = 2 1 α t 2 .
Yeh step kyun? Total angle ghuma hua rate curve ke area under se aata hai, rate× time se nahi, kyunki rate badhta rehta hai.
C ( t ) ko − ϕ ( t ) se z -rotation ke roop mein likho:
C ( t ) = cos ϕ ( t ) − sin ϕ ( t ) 0 sin ϕ ( t ) cos ϕ ( t ) 0 0 0 1 , ϕ ( t ) = 2 1 α t 2 .
Yeh step kyun? Ex 1 jaisi z -block structure, lekin angle ab integrated ϕ ( t ) hai, ω t nahi.
t = 1 , α = 2 plug karo: ϕ = 2 1 ⋅ 2 ⋅ 1 2 = 1 rad, toh cos 1 ≈ 0.5403 , sin 1 ≈ 0.8415 :
C ( 1 ) = cos 1 − sin 1 0 sin 1 cos 1 0 0 0 1 ≈ 0.5403 − 0.8415 0 0.8415 0.5403 0 0 0 1 .
Yeh maanga hua numerical matrix C ( 1 ) hai.
Verify: C C T = I aur det C = + 1 har t ke liye (VERIFY mein check kiya gaya). Warning: agar ω ( t ) ka axis khud tumble karta, toh skew matrices commute nahi karte aur yeh shortcut fail ho jaata — tumhe numerically ODE integrate karni padti (aur yehi wajah hai ki practice mein quaternions use kiye jaate hain).
Worked example Ex 9 — Spacecraft nadir spin
Ek imaging satellite apne body y -axis (pitch) ke baare mein ω = 0.10 rad/s se rotate karta hai taaki uska camera ground ki taraf point karta rahe. t = 0 par uska DCM C ( 0 ) = I hai. Ek body vector b ^ = ( 0 , 0 , 1 ) (camera boresight, shuru mein inertial z ke saath) t = 5 s ke baad inertial frame mein express karke kahan point karta hai?
Forecast: boresight x –z plane mein sweep karta hai; predict karo kis taraf.
y -spin ke liye Skew: ω B = ( 0 , 0.10 , 0 ) deta hai:
[ ω × ] = 0 0 − 0.10 0 0 0 0.10 0 0 .
Yeh step kyun? Sirf w 2 = 0.10 nonzero hai; y -axis spin axis hai.
C ( t ) = y ke baare mein − ω t ki rotation:
C ( t ) = cos ω t 0 sin ω t 0 1 0 − sin ω t 0 cos ω t .
Yeh step kyun? Cell A ka standard body-axis result, y -block ke liye adapted.
Inertial coords mein boresight hai b N = C T b ^ (kyunki C N→B map karta hai, iska transpose B→N map karta hai). ω t = 0.10 ⋅ 5 = 0.5 rad ke saath:
b N = C T ( 0 , 0 , 1 ) T = ( − sin 0.5 , 0 , cos 0.5 ) = ( − 0.4794 , 0 , 0.8776 ) .
Yeh step kyun? Hum boresight ground/inertial frame se dekhna chahte hain, toh body se inertial mein transform karte hain.
Verify: ∣ b N ∣ = 0.479 4 2 + 0.877 6 2 = 1.000 (unit length preserved) ✔. Camera 0.5 rad (≈ 28. 6 ∘ ) − x ki taraf tilt hua hai — ek sensible pitch sweep.
Worked example Ex 10 — Do galtiyan pakdo
Ek student likhta hai: "C : N → B aur ω N = ( 0 , 0 , ω ) ke liye, kinematics hai C ˙ = + [ ω N × ] C ." Dono galtiyan nikalo aur theek karo, phir sahi equation do.
Forecast: padhne se pehle galtiyan count karo — bilkul do hain.
Galti 1 — sign. N→B map karne wale matrix ke liye, sahi sign minus hai: C ˙ = − [ ⋅ ] C . Plus C ˙ T = [ ω × ] C T (B→N map) se belong karta hai. Kyun: forward body rotation ki wajah se inertial vectors body frame mein backward rotate hote hue dikhte hain.
Galti 2 — frame. Skew ko body components ω B use karne chahiye, ω N nahi, kyunki output side (C ˙ ka left index) body frame mein rehta hai.
Corrected equation:
C ˙ = − [ ω B × ] C , ω B = C ω N .
Yeh step kyun? Sign convention aur ω ka frame dono ko C ki mapping direction se match karna chahiye.
Verify: Agar hum (galat tarike se) plus rakhte, toh d t d ( C C T ) nonzero hoti, orthonormality violate hoti — ek contradiction jo sign error flag karta hai. Sahi minus ke saath, d t d ( C C T ) = 0 ✔. Quaternion kinematics — q̇ = ½ Ω(ω) q aur Euler angle kinematics & gimbal lock mein analogous convention care se compare karo.
Worked example Ex 11 — Kya flow
det C = + 1 rakhta hai?
Ex 1 ke constant z -spin solution use karke, confirm karo ki det C ( t ) = + 1 har t ke liye hai (naki − 1 , jo ek reflection hoti).
Forecast: rotations ko handedness preserve karni chahiye — lekin prove karo, assume mat karo.
Determinant compute karo Ex 1 se C ( t ) ka. Last row ke baad expand karte hue:
det C = 1 ⋅ ( cos 2 ω t + sin 2 ω t ) = cos 2 ω t + sin 2 ω t .
Yeh step kyun? z -row ( 0 , 0 , 1 ) hai, toh determinant 2 × 2 block mein reduce ho jaata hai.
Pythagorean identity: cos 2 + sin 2 = 1 , toh det C = + 1 har t ke liye.
Yeh step kyun? Yeh "proper rotation, no reflection" ki algebraic fingerprint hai.
General reason: d t d det C = det C ⋅ tr ( − [ ω × ]) = 0 kyunki har skew matrix ka trace zero hota hai. Toh det C apni starting value + 1 se kabhi nahi badlta.
Yeh step kyun? Jacobi's formula aur tr [ ω × ] = 0 (zero diagonal) determinant ko lock kar deta hai.
Verify: t = π /4 par (ω t = π /2 ): cos 2 + sin 2 = 0 + 1 = 1 ✔. Flow S O ( 3 ) mein rehta hai — koi reflection nahi, kabhi bhi. Yeh Poisson's equation structure guarantee karta hai.
Recall Kya har cell cover hua?
Single-axis (A) ::: Ex 1 (z ), Ex 2 (x ), Ex 9 (y ).
Zero spin (B) ::: Ex 3 — C ˙ = 0 , C frozen.
Sign flip (C) ::: Ex 4 — reverse spin transpose deta hai.
Limiting/periodic (D) ::: Ex 5 — bounded, period 2 π / ω .
General tilted axis (E) ::: Ex 6 — Rodrigues' formula, axis 2 1 ( 1 , 1 , 0 ) .
Extract ω , vee-map (F) ::: Ex 7 — ω = ( 3 , 4 , 5 ) .
Time-varying ω ( t ) (G) ::: Ex 8 — rate ko angle ϕ = 2 1 α t 2 mein integrate karo.
Word problem (H) ::: Ex 9 — 5 s ke baad boresight.
Exam twist (I) ::: Ex 10 — sign + frame traps.
Handedness (J) ::: Ex 11 — det C = + 1 hamesha.
ω B = ( 0 , 0 , ω ) aur C ( 0 ) = I ke liye, C ( t ) kya hai?z ke baare mein − ω t ki rotation: top-left block [[ cos ω t , sin ω t ] , [ − sin ω t , cos ω t ]] , C 33 = 1 .
Jab ω = 0 toh C ˙ ka kya hota hai? [ ω × ] = 0 , toh
C ˙ = 0 aur
C ( t ) = C ( 0 ) hamesha ke liye (degenerate case).
Spin reverse karna ω → − ω C ( t ) ke saath kya karta hai? Inverse/transpose rotation deta hai; sirf sin terms ka sign flip hota hai.
Constant-spin C ( t ) kabhi blow up kyun nahi karta? Skew matrices ke imaginary eigenvalues hote hain, toh
e − [ ω × ] t oscillate karta hai (period
2 π / ω ) — bounded,
S O ( 3 ) mein.
General constant axis e ^ ke liye, kaun sa closed form C ( t ) deta hai? Rodrigues: C = I − sin ϕ [ e ^ × ] + ( 1 − cos ϕ ) [ e ^ × ] 2 jahan ϕ = ω t .
Jab ω ( t ) vary karta hai lekin fixed axis rakhta hai, toh turn angle kaise nikalte hain? Rate integrate karo: ϕ ( t ) = ∫ 0 t ω ( s ) d s (jaise linear ramp ke liye 2 1 α t 2 ).
Time-varying ω ( t ) ke liye e − [ ω × ] t hamesha kyun use nahi kar sakte? Integral ka exponential valid hai sirf tab jab alag-alag times ke skew matrices commute karein (fixed axis); warna time-ordered exponential ya numerical integration chahiye.
− C ˙ C T se ω B kaise padhte hain?Vee-map: ω 1 = ( ⋅ ) 32 , ω 2 = ( ⋅ ) 13 , ω 3 = ( ⋅ ) 21 .
DCM flow ke under det C hamesha + 1 kyun rehta hai? d t d det C = det C ⋅ tr ( − [ ω × ]) = 0 kyunki har skew matrix ka trace zero hota hai.