Solving Kepler's equation — Newton-Raphson iteration: Foundations

3.2.15 · D1 · Physics › Orbital Mechanics & Astrodynamics › Solving Kepler's equation — Newton-Raphson iteration

Is central equation ko padhne se pehle, us mein har letter ka kuch visual matlab hona chahiye. Yeh page har letter ko ek ek karke, build-up order mein introduce karta hai, aur unhe define karne ke baad hi hum finally equation likhte hain. Upar se neeche padho — har symbol apne pehle waale par depend karta hai.

1. Ellipse aur uske do secret numbers: $a$ aur $b$

Ek ellipse ek aisa circle hai jise ek direction mein doosre se zyada stretch kiya gaya ho. Socho ek gol balloon ko dono haathon ke beech dabaa rahe ho.

Picture: ellipse ke centre par khade ho. Door wale end tak chalo — woh distance $a$ hai. 90° mudo aur paas wale edge tak chalo — woh distance $b$ hai. Kyunki ellipse dabaa hua hai, $b < a$ (circle ke liye, $b = a$ ).

Figure 1 (neeche) exactly yahi dikhata hai: ek violet ellipse jis mein magenta arrow $a$ ko long axis ke saath mark karta hai aur orange arrow $b$ ko short axis ke saath, plus centre aur off-centre focus jo hum aage milenge.

Is topic ko iski zaroorat kyun hai: poori orbit us box ke andar rehti hai jo $a$ aur $b$ set karte hain. Woh clever trick jo Kepler's equation derive karti hai woh ellipse ko uske around drawn ek perfect circle of radius $a$ se compare karti hai — toh $a$ aur $b$ pehle tools hain jo hum use karte hain.

Figure — Solving Kepler's equation — Newton-Raphson iteration

2. Focus aur offset $c$ — jahan Sun baitha hai

Ek ellipse mein ek special point hota hai jise focus kehte hain. Sun (ya Earth, satellite ke liye) ek focus par baitha hai — centre par nahin. Yahi off-centre-ness saari mushkil ki jad hai.

Picture: Figure 1 mein centre woh navy dot hai jo geometric middle mein hai; "mote" end ki taraf $c$ distance (dashed navy segment) slide karo aur tum orange star tak pahuncho — wahi Sun baitha hai. Agar $c = 0$ toh Sun bilkul centre mein hai aur orbit ek perfect circle hai.

Is topic ko iski zaroorat kyun hai: kyunki Sun $c$ se offset hai, body kabhi uske paas hoti hai aur kabhi door. Paas = fast, door = slow. Woh uneven speed exactly wahi reason hai ki hum time se position simple algebra se compute nahin kar sakte.

3. Eccentricity $e$ — "kitna dabaa hua" ka dial

Ab hum squash ko ek single dimensionless number mein pack karte hain.

Rearrange karne par milta hai $c = a e$ — focal offset $e$ use karke likha gaya. Pichle section ke triangle law $c^{2} = a^{2} - b^{2}$ ke saath combine karne par, yeh handy relation bhi milta hai $b = a 1 - e^{2}$ ( $c = a e$ substitute karo aur $b$ ke liye solve karo).

Dial padho:

$e = 0$ → focus centre par → perfect circle (koi squash nahin).
$e$ near $0$ (jaise $0.017$ Earth ke liye) → almost circle.
$e$ near $1$ (jaise $0.97$ kuch comets ke liye) → lamba, patla, cigar-shaped orbit.
$e = 1$ → ellipse ek parabola mein khul jaata hai (ab orbit nahin raha).

4. Radians mein angles — woh unit jo $sin$ ko addable banati hai

Topic ki central equation ek angle ko us angle ki sine ke saath add karegi. Yeh sense banane ke liye, angle radians mein hona chahiye, degrees mein nahin.

Picture: radius $1$ ka ek circle lo. Length $1$ ka ek string ka tukda uske rim par rakho. Us string ko span karne waala centre se angle exactly ek radian hai ( $\approx 57.3°$ ).

Is topic ko iski zaroorat kyun hai: aage mile teen anomalies radian angles hain, aur topic mein har $sin$ / $cos$ radians assume karta hai.

5. Teen anomalies: $M$ , $E$ , $ν$ — "kahan hai" kehne ke teen tarike

"Anomaly" ek purana astronomy word hai jiska matlab sirf angle jo bataye body kahan hai hai. Teen hain kyunki teen natural reference points hain.

Figure 2 (neeche) teeno ek saath dikhata hai: magenta ellipse (real orbit) apne violet dashed auxiliary circle ke andar. Body magenta dot hai; uska "shadow" seedha circle par push kiya gaya woh violet dot hai. Violet angle $E$ navy centre par measure kiya gaya hai, aur orange angle $ν$ orange Sun par measure kiya gaya hai — tum literally dekh sakte ho woh thodi alag directions mein point karte hain.

Is topic ko iski zaroorat kyun hai: mission hai time $\to$ position, yani $M \to E \to ν$ . $M$ time se aasaan hai, $ν$ $E$ se aasaan hai, lekin $M \to E$ woh transcendental knot hai jise whole topic untangle karne ke liye exist karta hai.

6. Mass, gravity, aur parameter $μ$ — pull ki strength

Time ko angle mein badalne se pehle, humein ek number chahiye jo capture kare ki central body kitni strongly pull karti hai.

Picture: $μ$ ko Sun ke liye ek single "gravity dial" ki tarah socho. Orbit work mein tumhe $G$ aur $M_{⋆}$ alag alag kabhi nahin chahiye — woh hamesha is ek product ke roop mein saath aate hain, toh astronomers inhe ek baar mein $μ$ mein bundle kar dete hain.

Is topic ko iski zaroorat kyun hai: mean angle ki steady rate jis par woh grow karta hai (aage define) depend karta hai gravity kitni strong hai, toh $μ$ table par hona chahiye pehle jab hum woh rate likh saken.

7. Mean motion $n$ , period $T$ , aur periapsis time $t_{p}$ — time ko $M$ mein convert karna

Picture: agar ek poora lap ( $2 π$ radians) time $T$ leta hai, toh har second $2 π / T$ radians milte hain. Woh rate $n$ hai. Us se, $M = n (t - t_{p}),$ elapsed radians since body ne last time time $t_{p}$ par closest approach pass kiya. Zyada Mean Motion and Orbital Period mein.

$a$ , $T$ aur $μ$ ke beech link jo $n = μ / a^{3} = 2 π / T$ mein chhupa hai woh exactly Kepler's third law hai.

8. Moving tools: $sin$ , $cos$ , derivative $f^{'}$ , aur root

Final ingredients woh mathematical machines hain jo topic run karta hai.

$sin E$ Kepler mein kyun aata hai? Derivation mein centre-to-focus-to-circle-point ka triangle height $a sin E$ rakhta hai. Woh height hi $E$ ko steady $M$ se door kheenchti hai — toh correction term $- e sin E$ ban jaata hai.

9. Central equation aur uska root function assemble karna

Ab jab $M$ , $E$ aur $e$ sab ke paas pictures hain, hum finally topic ka central relation likh sakte hain:

Hum ise rearrange karke $E$ akela nahin pa sakte, isliye hum woh $E$ dhundte hain jo ise balance kare. Sab kuch ek taraf move karo:

Derivative kyun bilkul? Newton-Raphson is slope ko use karta hai predict karne ke liye ki curve zero par kahan cross karegi: tangent line ko axis tak follow karo aur root ke paas land karo. Figure 3 (neeche) exactly yahi dikhata hai — violet curve $f (E)$ , uska magenta root, ek orange guess, aur dashed orange tangent jiska slope $f^{'}$ hai jo guess ko crossing ki taraf slide karta hai. Ise slower slope-free cousin se compare karo Fixed-Point Iteration for Kepler mein.

Yeh sab topic mein kaise feed karta hai

Equipment checklist

Khud test karo — right side cover karo aur zor se jawab do.

Ellipse par

a

aur

b

kya measure karte hain?

Sabse lambe (

a

) aur sabse chote (

b

) diameters ka aadha — semi-major aur semi-minor axes.

a

b

, aur

c

kaise related hain?

Right-triangle law

c^{2} = a^{2} - b^{2}

se; focal offset, minor axis, aur major axis ek right triangle banate hain jisme

a

hypotenuse hai.

Orbit par Sun kahan baitha hai?

Ek focus par, centre se

c = a e

se offset — kabhi centre par nahin.

Eccentricity

e

tumhe kya batata hai?

Orbit kitna squashed hai:

e = 0

circle hai,

e

near

1

ek lamba patla ellipse hai; yeh har correction term ko scale karta hai.

Har angle radians mein kyun hona chahiye?

Kyunki Kepler's equation ek raw angle ko ek sine ke saath add karti hai, aur sirf radians mein ek angle ek plain addable number ki tarah behave karta hai.

M

E

, aur

ν

mein kya difference hai?

M

= steady fake-body angle (time ka proxy);

E

= auxiliary-circle shadow tak ka angle centre par;

ν

= Sun par real angle.

Plain

arctan

galat

ν

kyun de sakta hai, aur tum ise kaise fix karte ho?

Kyunki

arctan

sirf

(- \frac{π}{2}, \frac{π}{2})

return karta hai aur top half bottom half se tell nahin kar sakta; atan2 use karo ya

sin E

ka sign carry karo taaki

E

aur

ν

same half-plane mein rahen.

t_{p}

kya hai aur yeh kyun hai?

Closest approach par clock reading;

t - t_{p}

measure karne se

M = 0

exactly perihelion par hota hai.

Gravitational parameter

μ

kya hai?

Product

μ = G M_{⋆}

Newton's constant aur central mass ka — ek single "gravity dial" jo set karta hai body ko kitna fast move karna hai.

Time

t

M

kaise nikalte ho?

M = n (t - t_{p})

jisme mean motion

n = 2 π / T = μ / a^{3}

hai.

Unit circle par

sin E

aur

cos E

kya hain?

Angle

E

par point ki height aur rightward reach; dono

[- 1, 1]

mein rehte hain.

f

ka "root" kya hai, aur

f^{'}

kya hai?

Root woh input hai jo

f = 0

banaye;

f^{'} = 1 - e cos E

f

ki local slope hai, Newton dwara root aim karne ke liye use kiya jaata hai.

Taiyaar ho? Main topic par wapas jao aur har symbol ab ek purane dost ki tarah feel hona chahiye.