Visual walkthrough — Modified Rodrigues parameters — singularity-free, compact

3.5.11 · D2 · Physics › Guidance, Navigation & Control (GNC) › Modified Rodrigues parameters — singularity-free, compact

---

Step 1 — Ek rotation sirf ek axis aur ek angle hai

KYA HAI. Koi bhi rigid object uthao — ek phone, ek spacecraft. Tum use ek pose se doosre pose mein jaise bhi ghuma lo, hamesha ek single seedhi line (the axis) hoti hai jiske around tum use ek baar mein ghuma sakte ho, aur ghoomne ki ek single matra (the angle) hoti hai. Yahi ek rotation ka poora matlab hai.

Hum in dono cheezoon ko naam dete hain:

• $\hat{e}$ — the axis, ek aisi arrow ke roop mein jis ki length exactly $1$ hai (ek "unit" arrow, isliye choti hat $\hat{}$). Yeh sirf batata hai ki spin-line kis direction mein point karti hai, kitni door nahi.
• $\Phi$ — the angle, degrees mein maapa gaya, jo batata hai ki hum $\hat{e}$ ke around kitna ghoomte hain.

KYUN. Har doosra attitude description (matrices, Euler angles, quaternions, MRPs) secretly isi $(\hat{e},\Phi )$ pair ki repackaging hai. Agar hum yahan se shuru karein, toh hum sachchi bottom se shuru karte hain.

Sign ke baare mein ek note — turn ka sense. Axis arrow yeh bhi fix karta hai ki kaun sa direction positive turn hai: apne right hand ki ungliyon ko $\hat{e}$ ke around curl karo jab thumb arrow ke saath point kar raha ho, aur tumhari ungliyan positive $\Phi$ ki direction dikhayengi. Toh $+\Phi$ ka $\hat{e}$ ke around rotation wahi motion hai jaise $-\Phi$ ka $-\hat{e}$ ke around rotation — axis ko flip karna aur angle ka sign flip karna cancel ho jaate hain. Yeh twin fact yaad rakho; yeh aage har sign decide karta hai.

PICTURE. Neeche, amber arrow $\hat{e}$ hai; object $\Phi$ angle se uske around swing karta hai (right-hand rule dikhaya gaya hai).

---

Step 2 — Sine aur cosine: ek right triangle padhna

KYA HAI. Isse pehle ki hum $(\hat{e},\Phi )$ ko numbers mein pack karein, humein ek angle ko length mein convert karne ka tarika chahiye. Yahi sine aur cosine ka kaam hai. Ek right triangle banao (ek corner ek square $90^{\circ }$ hai). Ek tilted corner chuno; uska angle $\theta$ hai.

• $\cos \theta =\frac{\text{adjacent}}{\text{hypotenuse}}$ — woh side jo $\theta$ ko touch karti hai, long slanted side se divided.
• $\sin \theta =\frac{\text{opposite}}{\text{hypotenuse}}$ — woh side jo $\theta$ ke across hai, long side par.

YEH TOOLS KYUN? Humein aise functions chahiye jo angle ke saath smoothly badhein aur hamesha jaane jaate hon. Sine aur cosine exactly "kitna aage" aur "kitna upar" gaye ho ek circle ke around — angle aur length ke beech natural dictionary hain.

PICTURE. Radius $1$ ke ek circle mein, angle $\theta$ tumhe point $(\cos \theta ,\sin \theta )$ par pohnchaata hai: cosine horizontal reach hai, sine vertical reach.

---

Step 3 — Half-angle packing: the quaternion

KYA HAI. Ab hum apne rotation ke liye pehla compact code banate hain. Angle ka aadha, $\Phi /2$, lo aur use cosine part aur sine part mein split karo:

$q_0=\cos \frac{\Phi }{2},\mathbf{q}=\hat{e}\sin \frac{\Phi }{2}$

• $q_0$ — ek single number (scalar part). Yeh aadhe angle ka cosine hai.
• $\mathbf{q}$ — ek chota arrow (vector part): axis $\hat{e}$ ko $\sin (\Phi /2)$ se stretch kiya gaya. Yeh axis ki tarah same direction mein point karta hai, lekin uski length ab angle kitna bada hai woh record karti hai.

Saath mein yeh chaar numbers $(q_0,\mathbf{q})$ quaternion banate hain — dekho Quaternions (Euler symmetric parameters).

HALF-ANGLE KYUN — aur yeh rotations ko sahi combine kyun karta hai. Yeh woh picture hai jis par ruk ke dhyan dene laayak hai. Ek rotation karo, phir doosra identical rotation karo — total usi axis ke around $2\Phi$ ka turn hai. Dekho vector part is doubling ke under kya karta hai. Ek single quaternion $\sin (\Phi /2)$ carry karta hai; lekin trig doubling rule kehta hai $\sin \Phi =2\sin (\Phi /2)\cos (\Phi /2)$ — poora angle $\Phi$ woh hai jo tum do half-angle pieces ko multiply karke paate ho. Doosre shabdon mein, half-angle exactly woh exponent hai jo turns compose karne par add hota hai: do rotations stack karo aur unke half-angles add hote hain ($\Phi /2+\Phi /2=\Phi$), jo sine/cosine machinery sahi combined rotation mein convert karti hai. Agar hum full angle ki jagah half store karte, toh do rotations compose karna sirf codes multiply karne ke barabar nahi hota — half woh unique choice hai jo "quaternions multiply karna" ka matlab "ek turn phir doosra karna" banata hai. Aur kyunki half-angle ke $\cos$ aur $\sin$ dono kabhi ek saath zero nahi hote, code kabhi blow up nahi hota: kahi koi singularity nahi.

PICTURE. Jab $\Phi$ $0\to 360^{\circ }$ sweep karta hai, half-angle sirf $0\to 180^{\circ }$ sweep karta hai. Dekho $q_0=\cos (\Phi /2)$ $+1$ se start hota hai, $\Phi =180^{\circ }$ par $0$ cross karta hai, aur $\Phi =360^{\circ }$ par $-1$ tak pohnchhta hai; doubling rule $\sin \Phi =2\sin (\Phi /2)\cos (\Phi /2)$ ko amber curve ke roop mein draw kiya gaya hai jo do cyan ones se bana hai.

---

Step 4 — Pehla attempt: $q_0$ se divide karo (Gibbs vector, aur yeh kyun fail hota hai)

KYA HAI. Chaar numbers ek zyada hain. Teen tak pohnchhne ke liye, hum arrow $\mathbf{q}$ ko scalar $q_0$ se divide karte hain:

$\mathbf{g}=\frac{\mathbf{q}}{q_0}=\frac{\hat{e}\sin (\Phi /2)}{\cos (\Phi /2)}=\hat{e}\tan \frac{\Phi }{2}$

Yahan ek naya tool aata hai: tangent.

• $\tan \theta =\frac{\sin \theta }{\cos \theta }$ — sine over cosine, "kitna steep." Jab cosine (neeche wala) $0$ tak shrink hota hai, tangent infinity ki taraf shoot karta hai.

Yeh $\mathbf{g}=\hat{e}\tan (\Phi /2)$ classical Rodrigues parameters (Gibbs vector) hain — dekho Classical Rodrigues parameters (Gibbs vector).

DIVIDE KYUN? Divide karne se $\hat{e}$-arrow ki length ka hypotenuse par dependence khatam ho jaata hai aur axis ke saath ek clean 3-number vector milta hai, jis ki length angle ka pure function hai. Teen numbers achieve ho gaaye.

YEH FAIL KYUN HOTA HAI. Denominator $q_0=\cos (\Phi /2)$ dekho. $\Phi =180^{\circ }$ par, half-angle $90^{\circ }$ hai, aur $\cos 90^{\circ }=0$. Zero se divide — Gibbs vector sirf $180^{\circ }$ par explode karta hai. Ek spacecraft ke liye jo half turn se aage tumble kar sakta hai, yeh bahut jaldi hai.

PICTURE. Tangent curve infinity ki taraf rocket karti hai jab uska angle $90^{\circ }$ ke paas aata hai; woh vertical wall hi $\Phi =180^{\circ }$ par Gibbs singularity hai.

---

Step 5 — Fix: $q_0$ ki jagah $1+q_0$ se divide karo

KYA HAI. Teen-number idea rakho, lekin divisor ko $q_0$ se $1+q_0$ mein badlo:

$\mathbf{\sigma }=\frac{\mathbf{q}}{1+q_0}$

• $\mathbf{\sigma }$ — Modified Rodrigues Parameters, hamare target teen numbers.
• $1+q_0$ — naya denominator. Yeh $1+\cos (\Phi /2)$ ke barabar hai.

YEH DIVISOR KYUN? Poochho: $1+q_0$ kab zero hit karta hai? Sirf tab jab $q_0=-1$, yaani $\cos (\Phi /2)=-1$, yaani $\Phi /2=180^{\circ }$, yaani $\Phi =360^{\circ }$ — ek full turn. Gibbs se compare karo (jo $180^{\circ }$ par mara) aur Euler angles se (gimbal lock $90^{\circ }$ par). $1$ add karna blow-up ko $360^{\circ }$ tak shift kar deta hai, jo utna door ho sakta hai jitna possible hai.

PICTURE. Denominators saath mein plot karo: $q_0$ $\Phi =180^{\circ }$ par zero cross karta hai (Gibbs wahan marta hai), jabki $1+q_0$ sirf $\Phi =360^{\circ }$ par zero touch karta hai (MRP poori normal range survive karta hai).

---

Step 6 — $\mathbf{q}/(1+q_0)$ exactly $\hat{e}\tan (\Phi /4)$ kyun hai

KYA HAI. Quaternion pieces substitute karo aur simplify karo:

$\mathbf{\sigma }=\frac{\hat{e}\sin (\Phi /2)}{1+\cos (\Phi /2)}=\hat{e}\tan \frac{\Phi }{4}$

Middle-to-right jump mein ek trigonometric identity use hoti hai:

$\frac{\sin \theta }{1+\cos \theta }=\tan \frac{\theta }{2}$

Identity earn karo — ek do-line algebraic sketch. $\theta =2\alpha$ likho, toh $\alpha =\theta /2$ woh $\theta$ ka half hai (jo $\Phi /4$ ban jaata hai jab $\theta =\Phi /2$). Ab do schoolbook doubling rules use karo:

$\sin \theta =\sin 2\alpha =2\sin \alpha \cos \alpha ,1+\cos \theta =1+\cos 2\alpha =2{\cos }^2\alpha .$

Pehla kehta hai "ek full sine do half-pieces multiply hain"; doosra standard $\cos 2\alpha =2{\cos }^2\alpha -1$ ko rearrange karta hai toh $-1$, $+1$ cancel ho jaata hai. Ek ko doosre se divide karo — ek factor $2\cos \alpha$ upar aur neeche cancel ho jaata hai:

$\frac{\sin \theta }{1+\cos \theta }=\frac{2\sin \alpha \cos \alpha }{2{\cos }^2\alpha }=\frac{\sin \alpha }{\cos \alpha }=\tan \alpha =\tan \frac{\theta }{2}.$

$\theta =\Phi /2$ set karo. Tab $\alpha =\theta /2=\Phi /4$ — quarter angle. Yahi se famous $\tan (\Phi /4)$ aata hai.

• $\hat{e}$ — abhi bhi sirf axis direction (scalar se divide karna arrow ko rotate nahi kar sakta).
• $\tan (\Phi /4)$ — $\mathbf{\sigma }$ ki length, ab quarter angle ka tangent.

YEH MATTER KYUN KARTA HAI. Quarter-angle hi woh reason hai ki singularity $360^{\circ }$ par hai: $\tan (\Phi /4)$ tab blow up karta hai jab $\Phi /4=90^{\circ }$, yaani $\Phi =360^{\circ }$. "Half of the half angle" = doosri halving = singularity Gibbs se do baar zyada door push ho gayi.

PICTURE. Half-angle right triangle ek circle ke andar baitha hai; identity $\sin \theta /(1+\cos \theta )$ literally point $(-1,0)$ se $(\cos \theta ,\sin \theta )$ tak chord ki slope hai — ek chord jo $\theta$ ka half subtend karti hai, $\tan (\theta /2)$ deti hai. Cancelled factor $2\cos \alpha$ bhi draw kiya gaya hai.

---

Step 6b — $\tan (\Phi /4)$ ka sign kya matlab rakhta hai

KYA HAI. Working range $0^{\circ }\le \Phi <360^{\circ }$ mein quarter angle $0^{\circ }\le \Phi /4<90^{\circ }$ mein rehta hai — first quadrant — jahan $\tan (\Phi /4)$ $\ge 0$ hota hai. Toh forward turn ke liye scalar length $\tan (\Phi /4)$ kabhi negative nahi hoti, aur poori direction (including which way you spin) axis $\hat{e}$ carry karta hai, exactly jaisa Step 2 mein sine ne sign carry kiya tha.

YEH CLEAN DESIGN KYUN HAI. Kyunki quarter angle normal use mein kabhi quadrant I nahi chodta, hum tangent ke sign-flipping se kabhi nahi ladte (tangent quadrant II mein negative ho jaata hai, III mein phir positive — wahi repeat-every-$180^{\circ }$ trouble jo arctan ko kaatta hai). MRPs use sidestep karte hain: $\tan (\Phi /4)\ge 0$ rakho aur $\hat{e}$ ko signing karne do.

Negative angles / reverse spin. $-\Phi$ ka turn (doosri taraf spin karo) Step 1 ke twin fact se $-\hat{e}$ ke around $+\Phi$ ke barabar hai. Woh feed karo: $\mathbf{\sigma }=(-\hat{e})\tan (\Phi /4)=-\hat{e}\tan (\Phi /4)$ — MRP simply sign flip kar deta hai. Toh rotation ka sense reverse karna $\mathbf{\sigma }$ ko origin ke through flip karta hai, jo exactly wahi hai jo tumhara intuition chahta hai: opposite turn ⇒ opposite chota arrow.

PICTURE. Do arrows: $\hat{e}$ ke around $+\Phi$ turn $\mathbf{\sigma }$ deta hai; reverse turn $-\mathbf{\sigma }$ deta hai, origin ke through mirror image.

---

Step 7 — $\Phi$ badhne par $\Vert \mathbf{\sigma }\Vert$ ka har case

KYA HAI. Kyunki $\Vert \mathbf{\sigma }\Vert =\tan (\Phi /4)$ (range mein non-negative), $\Phi$ ka poora span dekho aur dekhlo ki length kya karti hai:

$\Phi$	$\Phi /4$	$\Vert \mathbf{\sigma }\Vert =\tan (\Phi /4)$	matlab
$0^{\circ }$	$0^{\circ }$	$0$	identity, $\mathbf{\sigma }=0$
$120^{\circ }$	$30^{\circ }$	$0.577$	chota, safe
$180^{\circ }$	$45^{\circ }$	$1$	switch threshold
$300^{\circ }$	$75^{\circ }$	$3.73$	dangerously large
$\to 360^{\circ }$	$\to 90^{\circ }$	$\to \infty$	singularity

THRESHOLD $1$ KYUN HAI. $\Phi =180^{\circ }$ par quarter angle $45^{\circ }$ hai, aur $\tan 45^{\circ }=1$ exactly. Toh $\Vert \mathbf{\sigma }\Vert =1$ natural halfway warning line hai: iske neeche safe ho, iske upar blow-up ki taraf ja rahe ho.

Degenerate case $\Phi =0$. Tab $\Phi /4=0$, $\tan 0=0$, toh $\mathbf{\sigma }=0$ axis $\hat{e}$ chahey jo bhi ho — kyunki kuch nahi karna ka koi meaningful axis nahi hota. Origin identity rotation hai. Achha: wahan koi division problem nahi kyunki $1+q_0=2$.

PICTURE. $\Vert \mathbf{\sigma }\Vert$ ki number line jismein $=1$ gate $180^{\circ }$ par hai aur wall $360^{\circ }$ par.

---

Step 8 — Shadow set: last wall se bachna

KYA HAI. Step 3 se yaad karo ki $q$ aur $-q$ same rotation hain. Toh negated quaternion $-q=(-q_0,-\mathbf{q})$ ko usi same recipe $\mathbf{\sigma }=\mathbf{q}/(1+q_0)$ mein feed karo. Denominator dhyan se dekho: scalar part ab $-q_0$ hai, toh "$1+q_0$" "$1+(-q_0)=1-q_0$" ban jaata hai, aur vector part $-\mathbf{q}$ hai:

${\mathbf{\sigma }}^S=\frac{-\mathbf{q}}{1+(-q_0)}=\frac{-\mathbf{q}}{1-q_0}.$

Isliye denominator $1-q_0$ mein flip hota hai — yeh same formula hai jisme sign-twin quaternion feed kiya gaya. Ab ise sirf $\mathbf{\sigma }$ mein formula banao. ${\mathbf{\sigma }}^S$ ke upar aur neeche $(1+q_0)$ se multiply karo aur $(1-q_0)(1+q_0)=1-q_0^2=\Vert \mathbf{q}{\Vert }^2$ use karo (kyunki $q_0^2+\Vert \mathbf{q}{\Vert }^2=1$), saath mein $\mathbf{\sigma }=\mathbf{q}/(1+q_0)$ aur $\Vert \mathbf{\sigma }{\Vert }^2=\Vert \mathbf{q}{\Vert }^2/(1+q_0{)}^2$:

${\mathbf{\sigma }}^S=\frac{-\mathbf{q}(1+q_0)}{\Vert \mathbf{q}{\Vert }^2}=-\frac{\mathbf{\sigma }}{\Vert \mathbf{\sigma }{\Vert }^2}=-\frac{\mathbf{\sigma }}{{\mathbf{\sigma }}^{\top }\mathbf{\sigma }}.$

• ${\mathbf{\sigma }}^{\top }\mathbf{\sigma }=\Vert \mathbf{\sigma }{\Vert }^2$ — $\mathbf{\sigma }$ ki length squared.
• Minus sign aur divide-by-length-squared vector ko unit sphere ke through invert karta hai: lamba $\mathbf{\sigma }$ chota ${\mathbf{\sigma }}^S$ ban jaata hai.

KYUN. Agar $\Vert \mathbf{\sigma }\Vert >1$ (past $180^{\circ }$, $360^{\circ }$ wall ki taraf chadh raha hai), ${\mathbf{\sigma }}^S$ par switch karo. Kyunki $\Vert {\mathbf{\sigma }}^S\Vert =1/\Vert \mathbf{\sigma }\Vert <1$, tum instantly safe zone mein wapas aa jaate ho — same pose describe karte hue. Shadow $\Phi \mapsto 360^{\circ }-\Phi$ encode karta hai, toh $300^{\circ }$ turn ek comfortable $60^{\circ }$ wala ban jaata hai. Isliye MRPs "practically singularity-free" hain.

PICTURE. Unit circle ke bahar ek lamba red $\mathbf{\sigma }$ ek chote cyan ${\mathbf{\sigma }}^S$ mein map hota hai iske andar — origin se wohi line, flipped side.

---

Ek-picture summary

Upar sab kuch, compressed: axis–angle (with sign) → half-angle quaternion (composition ke under add hota hai) → $1+q_0$ se divide karo → quarter-angle tangent → unit circle ke baad shadow-set flip (origin par singular).

Recall Feynman retelling — seedhe words mein bolo

Ek rotation kuch nahi hai bas ek spin-line aur tum uspar kitna ghoomte ho, aur us line ke around right-hand rule batata hai kaun sa direction positive hai — toh doosri taraf ghoomna waisa hi hai jaise line flip karna. Use numbers mein store karne ke liye, main pehle angle ko aadha karta hoon aur likhhta hoon "kitna flat" aur "kitna tall" woh half-angle hai — yahi quaternion hai, chaar safe numbers. Half KYUN? Kyunki jab main ek turn do baar karta hoon, half-angles simply add ho jaate hain, aur sine-doubling rule $\sin \Phi =2\sin (\Phi /2)\cos (\Phi /2)$ exactly wahi hai jo "codes multiply karna" ka matlab "ek turn phir doosra karna" banata hai — half woh magic exponent hai jo add hota hai. Chaar ek zyada hai, toh main teen tak pohnchhne ke liye tall part ko flat part se divide karne ki koshish karta hoon; lekin flat part half turn par zero hit karta hai aur poori cheez explode ho jaati hai — yahi purana Gibbs vector hai, useless. Clever fix hai "one plus flat part" se divide karna; woh denominator sirf full turn par zero hota hai. Do doubling rules — $\sin \theta =2\sin \alpha \cos \alpha$ aur $1+\cos \theta =2{\cos }^2\alpha$ — se algebra clean karne par, ek factor cancel ho jaata hai aur main axis times quarter angle ka tangent pe aa jaata hoon. Length woh quarter-tangent hai: rest par zero, half turn par one, full turn par infinity, aur uska sign axis mein rehta hai, tangent mein nahi. Finally, kyunki quaternion ka sign flip karna same rotation mean karta hai, sign-twin ko usi recipe mein feed karne se "$1+q_0$" "$1-q_0$" ban jaata hai aur ek backup "shadow" version milta hai; jab bhi mera vector length one se bada ho jaata hai toh main shadow par swap kar leta hoon, jo chota aur safe land karta hai — lekin main swap kabhi origin ke paas nahi karta, kyunki "no turn" ka shadow khud infinite hai. Isliye teen numbers enough hain aur main practically koi wall kabhi nahi milta.

One-line summary ::: $\mathbf{\sigma }=\hat{e}\tan (\Phi /4)$ axis–angle se, $1+q_0$ se divide kiye gaye quaternion ke zariye; shadow-set flip $\Vert \mathbf{\sigma }\Vert =1$ se aage lone $360^{\circ }$ singularity se bachata hai. Quarter-angle half-angle kyun nahi ::: $q_0$ ki jagah $1+q_0$ se divide karne se (not $q_0$) $\tan (\Phi /2)$ $\tan (\Phi /4)$ ban jaata hai, singularity $180^{\circ }$ se $360^{\circ }$ shift ho jaati hai. Shadow denominator $1-q_0$ kyun hai ::: yeh same formula $\mathbf{q}/(1+q_0)$ hai jisme sign-twin $-q$ feed kiya gaya, toh $q_0\to -q_0$; yeh identity par singular hai, toh sirf tab switch karo jab $\Vert \sigma \Vert >1$.

Dekho bhi: Direction Cosine Matrix (DCM), Attitude kinematics and $\boldsymbol\omega$, Spacecraft attitude control laws.