Converting between orbital elements and state vectors (r, v): Question bank

3.2.17 · D5 · Physics › Orbital Mechanics & Astrodynamics › Converting between orbital elements and state vectors (r, v)

Definition Is page pe kaam aane waale symbols aur tools (taaki aap kahin aur jump na karo)

$μ = GM$ — gravitational parameter: central body ka mass $M$ times Newton's constant $G$ , ek number mein bundle kiya hua (km³/s²). Ye gravity ki pull ki strength set karta hai; neeche har speed, energy, aur period formula isko carry karta hai.
$ε = \frac{v ^{2}}{2} - \frac{μ}{r}$ — specific energy (energy per kg): kinetic minus gravitational potential. Bound orbits ke liye negative, parabolic escape ke liye zero, hyperbolic ke liye positive.
$h = r \times v$ — specific angular momentum: orbit plane ke perpendicular.
$e$ — eccentricity vector: periapsis ki taraf point karta hai, magnitude $= e$ .
$n = \overset{z}{^} \times h$ — node vector: ascending node ki taraf point karta hai.
Rotation matrices (recipe ke step-B tools). $R_{3} (θ)$ ek vector ko $\overset{z}{^}$ axis ke baare mein angle $θ$ se rotate karta hai; $R_{1} (θ)$ $\overset{x}{^}$ axis ke baare mein rotate karta hai. Figure s01 dekho ki "us axis ke baare mein" aur "positive angle" ka matlab kya hai.

Ek positive rotation counter-clockwise hota hai jab aap axis ke neeche origin ki taraf dekhte ho (right-hand rule). Recipe negative-angle versions use karti hai kyunki ye frame build-up ko undo karta hai (neeche "Spot the error" note aur figure s01 dekho).

Figure — Converting between orbital elements and state vectors (r, v)

True or false — justify

Chhe numbers, ek ellipse. Decide karo, phir defend karo.

State vectors

(r, v)

aur elements

(a, e, i, Ω, ω, ν)

bilkul same information carry karte hain.

True in general — dono 6 numbers hain jo same orbit describe karte hain, aur conversion koi physics nahi khoता. Caveat: circular (

e = 0

) ya equatorial (

i = 0

) singularities pe kuch angles (

ω

Ω

) undefined ho jaate hain, isliye classical element set momentarily incomplete ho jaata hai bhale hi

(r, v)

bilkul valid rehta hai.

Eccentricity vector

e

hamesha focus se apoapsis ki taraf point karta hai.

False — ye periapsis ki taraf point karta hai, yaani sabse paas ka point. Periapsis pe

r

e

ke along hota hai, isliye

cos ν = e \cdot r / (er) = 1

aur

ν = 0

Kyunki

h = r \times v

constant hai, two-body problem mein orbit plane kabhi nahi badalta.

True — constant

h

ka matlab hai ek fixed perpendicular direction, isliye

r

aur

v

dono us ek plane mein rehte hain jiska normal

h

hai; yahi wajah hai ki

h

i

aur

Ω

fix karta hai.

Specific angular momentum

h

orbit ki motion ke saath point karta hai.

False —

h

cross product

r \times v

hai, isliye ye motion ke plane ke perpendicular hota hai (ye orbit ke "bahar" point karta hai), path ke along nahi.

Bada semi-major axis

a

hamesha faster satellite ka matlab hota hai.

False — vis-viva se, bada

a

matlab average pe lower total speed (unchi, dheemi orbits); energy

ε = - μ / (2 a)

a

badhne ke saath zero ki taraf jaati hai, matlab diye gaye

r

ke liye kam kinetic energy.

Agar kisi instant pe

r \cdot v = 0

ho, toh body periapsis ya apoapsis pe hai.

True —

r \cdot v

distance ke change ki rate hai; zero matlab

r

momentarily na badh raha hai na ghad raha hai, jo sirf ellipse ke do turning points pe hota hai.

Node vector

n = \overset{z}{^} \times h

periapsis ki taraf point karta hai.

False —

n

ascending node ki taraf point karta hai (jahan orbit upar jaate hue reference plane cross karta hai). Periapsis

e

se fix hota hai, aur

n

aur

e

ke beech ka angle hi argument of periapsis

ω

hai.

Hyperbolic trajectory ke liye semi-major axis

a

negative aata hai.

True — hyperbolas ka positive specific energy

ε > 0

hota hai, aur

a = - μ / (2 ε)

tab

a < 0

force karta hai; formula bina kisi special case ke saare conics handle karta hai.

Do orbits jinke identical

(a, e, i, Ω, ω)

hain lekin alag

ν

hai, alag orbits hain.

False — same orbit, sirf usme alag positions;

ν

ek aisa element hai jo time ke saath badalta hai, isliye ye path nahi, point batata hai.

Spot the error

Har line mein ek plausible-sounding mistake hai. Use name karo aur fix karo. (Rotation-axis convention ke liye figure s01 dekho.)

"Perifocal se inertial jaane ke liye,

R_{3} (+ ω) R_{1} (+ i) R_{3} (+ Ω)

se rotate karo."

Galat signs aur order. Perifocal frame ko inertial se build karna ek swivel

Ω

, tilt

i

, spin

ω

hai; doosri taraf jaana inverse hai, aur rotation product ka inverse order reverse karta hai aur har angle negate karta hai. Toh correct transform hai

R_{3} (- Ω) R_{1} (- i) R_{3} (- ω)

"Maine

ν = cos^{- 1} (\frac{e \cdot r}{er}) = 4 8^{\circ}

use kiya, ho gaya."

Quadrant check missing hai —

cos^{- 1}

sirf

[0, π]

return karta hai. Agar

r \cdot v < 0

ho (periapsis ki taraf ja raha hai) toh sahi value

36 0^{\circ} - 4 8^{\circ}

hai. Zyada robust:

ν = atan2 (r \cdot v / (scaled), e \cdot r)

style two-argument arctangent compute karo, jo ek hi step mein bina kisi sign patch ke poora

[0, 2 π)

return karta hai. Figure s03 dekho.

e = 0.02

hai isliye argument of periapsis by definition

ω = 0

hai."

Galat reasoning — chota

e

ω

ko zero nahi fix karta;

ω

abhi bhi node se periapsis tak ek well-defined angle hai. Sirf jab

e = 0

exactly ho tabhi

ω

undefined hota hai.

cos i = h_{z} / h

, aur kyunki cosine negative ho sakta hai isliye mujhe

i

ka sign alag se check karna hoga."

i

ke liye koi quadrant check nahi chahiye — inclination

[0, π]

mein rehta hai, jo exactly wahi range hai jo

cos^{- 1}

return karta hai, isliye

cos i = h_{z} / h

i

ko uniquely aur bina kisi ambiguity ke deta hai.

v_{pq w}

ke liye maine

μ / p (- sin ν, cos ν, 0)

likha."

e

drop ho gaya — tangential component

e + cos ν

hai, na ki sirf

cos ν

. Ise bhoolne se velocity khaas points ke alaawa har jagah galat ho jaati hai, aur energy conservation violate hoti hai.

a = μ / (2 ε)

orbit ka size energy se deta hai."

Sign slip — sahi relation hai

a = - μ / (2 ε)

. Minus kyun: bound orbit ke liye specific energy negative hoti hai (

ε < 0

, body potential well mein trapped hai), isliye

a = - μ / (2 ε)

use positive size mein flip karta hai. Plus sign ek ordinary ellipse ke liye nonsensical negative

a

deta.

"RAAN

Ω

periapsis se measure kiya jaata hai."

Galat reference. $\overset{x}{^}$ se kyun:

Ω

describe karta hai ki orbit plane kaise swivel hoti hai, jo poori tarah se yahan set hoti hai ki plane equator kahan cross karti hai — ascending node. Woh crossing reference plane mein fixed inertial direction

\overset{x}{^}

(reference/vernal direction) se measure ki jaati hai. Periapsis plane ke andar rehta hai aur alag se

ω

(node se) aur

ν

(periapsis se) se measure hota hai.

Why questions

Parent recipe cheezein karti hai — yahan bataya gaya hai ki inhe skip kyun nahi kar sakte. (Step numbers upar di gayi one-glance recipe se hain.)

Steps 6, 7, 8 ko companion-sign check kyun chahiye lekin step 5 (inclination) ko nahi?

Kyunki

i \in [0, π]

exactly

cos^{- 1}

ke output range se match karta hai, lekin

Ω, ω, ν

poore

[0, 2 π)

mein phele hain, isliye ye jaanne ke liye ki angle circle ke kis half mein hai, ek doosri quantity ka sign chahiye. Figure s03 dekho.

True anomaly

ν

ke liye

r \cdot v

sahi sign-flag kyun hai?

Kyunki

r \cdot v > 0

matlab body focus se door ja rahi hai, yaani periapsis se apoapsis ki taraf chadh rahi hai, toh

ν \in (0, π)

; negative value matlab utar rahi hai, toh

ν \in (π, 2 π)

Elements → state convert karna "PQW mein build karo, phir rotate karo" mein kyun split hota hai?

Perifocal frame mein periapsis ek axis pe hota hai, isliye

r

aur

v

ke dead-simple

cos ν / sin ν

forms hote hain; phir saara messy orientation (

Ω, i, ω

) baad mein ek single rotation matrix se milta hai.

Eccentricity vector ko derive karna memorise karne se zyada zaroori kyun hai?

Iska derivation dikhata hai ki ye genuinely ek constant of motion hai (two-body equation ko

h

se cross karne se ek fixed vector integrate hota hai), yahi wajah hai ki ye reliably

e

aur periapsis direction

ω

dono encode kar sakta hai.

h

do elements (

i

aur

Ω

) fix karta hai jabki

ε

sirf ek (

a

) kyun fix karta hai?

h

ek vector hai — iska direction plane ki orientation pin karta hai (tilt

i

aur swivel

Ω

dono deta hai), jabki

ε

ek single scalar hai aur sirf ek size number,

a

, set kar sakta hai.

Vis-viva relation ek single speed measurement (

r

ke saath) se

a

kyun reveal karne deta hai?

Kyunki

v^{2} = μ (2/ r - 1/ a)

speed, position, aur size ko saath bandh karta hai;

v

aur

r

jaanne ke baad sirf

a

unknown rehta hai — Two-body problem and the vis-viva equation dekho.

Same rotation matrix

\tilde{R}

ko dono

r

aur

v

pe kyun apply karna padta hai?

Ye dono vectors same perifocal frame mein express hain; frames change karna ek geometric re-orientation hai, isliye ye us frame mein rehne waale har vector pe identically act karta hai.

Edge cases

Woh scenarios jahan naive formulas quietly break ho jaate hain. (Figure s04 do singular geometries dikhata hai.)

Circular orbit (

e = 0

): kaunsa element undefined ho jaata hai, aur ise kaise patch kiya jaata hai?

ω

undefined ho jaata hai kyunki measure karne ke liye koi periapsis nahi hai (

e = 0

). Patch: argument of latitude

u = ω + ν

use karo, node se body tak ka angle. Orbital singularities and equinoctial elements dekho.

Equatorial orbit (

i = 0

): kya fail hota hai, aur kyun?

n = \overset{z}{^} \times h = 0

kyunki

h

\overset{z}{^}

ke parallel hai, isliye koi ascending node nahi hai aur

Ω

undefined hai; fix ye hai ki

Ω

aur

ω

alag se track karne ki jagah true longitude track karo.

Dono

e = 0

aur

i = 0

(circular equatorial): kitne elements collapse ho jaate hain?

ω

aur

Ω

dono disappear ho jaate hain, aur sirf

a

aur ek single angle (true longitude

ℓ = Ω + ω + ν

) meaningful rehte hain; equinoctial elements precisely inhi dual singularities se bachne ke liye bane hain.

Parabolic escape (

e = 1

a

ka kya hota hai?

Energy

ε = 0

hai, isliye

a = - μ / (2 ε)

blow up ho jaata hai (infinite) — ek parabola ka koi finite semi-major axis nahi hota, aur iske bajaye semi-latus rectum

p = h^{2} / μ

use kiya jaata hai.

Exact apoapsis pe,

r \cdot v = 0

— recipe

ν = 0

ke bajaye

ν = π

kyun choose karti hai?

Dot product dono turning points pe zero hota hai, isliye sign rule inconclusive hai; aap ise

cos ν = e \cdot r / (er)

se resolve karte ho, jo apoapsis pe

- 1

hai (

ν = π

deta hai) aur periapsis pe

+ 1

Radial (straight-line) trajectory,

r ∥ v

: poora scheme kyun break ho jaata hai?

Tab

h = r \times v = 0

hai, isliye orbit plane,

i

Ω

, aur

p = h^{2} / μ

sab undefined hain — ek degenerate "orbit" jisme koi angular momentum nahi hai aur jise element set represent hi nahi kar sakta.

Retrograde orbit (

i

9 0^{\circ}

aur

18 0^{\circ}

ke beech): kya

cos i = h_{z} / h

abhi bhi kaam karta hai?

Haan —

h_{z} < 0

cosine ko negative banata hai aur

cos^{- 1}

9 0^{\circ}

se upar ka angle return karta hai, bina kisi extra sign fix ke correctly retrograde motion flag karta hai (yahi wajah hai ki

i

ko quadrant check nahi chahiye).

Nearly-zero eccentricity jahan numerical noise

e \cdot r / (er)

1

se thoda exceed kara deta hai: safe habit kya hai?

cos^{- 1}

call karne se pehle argument ko

[- 1, 1]

mein clamp karo; floating-point round-off ise valid range se bahar push kar sakta hai aur ek otherwise correct computation crash kar sakta hai.

Recall Har trap ki one-sentence summary

Baar baar aane wale lessons: $h$ plane ke bahar point karta hai (along nahi), $e$ periapsis ki taraf point karta hai, semi-major axis minus sign use karta hai ( $a = - μ /2 ε$ ), $cos^{- 1}$ ko $i$ ke alaawa companion sign chahiye, PQW→ECI transform negative angles use karta hai, aur jab bhi $e = 0$ ya $i = 0$ ho koi element dissolve ho jaata hai aur use $u$ , longitude, ya equinoctial elements se swap karna padta hai.