Gauss's method for Lambert's problem: Question bank

3.2.29 · D5 · Physics › Orbital Mechanics & Astrodynamics › Gauss's method for Lambert's problem

Shuru karne se pehle, ek symbol pantry taaki neeche koi bhi symbol unexplained na ho:

Formula Is page par refer ki gayi constants aur functions (reference ke liye)

Geometry / time constants (sirf given data se bante hain):

m = \frac{\mu\,\Delta t^2}{\big(2\sqrt{r_1 r_2}\,\cos(\theta/2)\big)^3}.$$ $\ell$ **shape-constant** hai (pure geometry); $m$ **time-constant** hai (ismein $\Delta t$ hai, $\Delta t^2$ note karo). **Specific angular momentum** $h=\sqrt{\mu p}$ = swept-area ka *rate* per unit mass; Kepler's second law padhta hai $\Delta t = 2\,(\text{sector area})/h$, yaani area steady rate $h/2$ par accumulate hoti hai. **Semi-latus rectum** $p$ = orbit ka "width" scale (focus par ellipse ke across chord ka aadha); ye size/energy fix karta hai aur iteration ke *baad* hi pata chalta hai: $$p = \frac{r_1 r_2 (1-\cos\theta)}{r_1+r_2-2\sqrt{r_1 r_2}\cos(\theta/2)\,(1-2x)}.$$ **Eccentric anomaly** $E$ = ellipse ko circumscribe karne wale circle par ek auxiliary angle jo Kepler's equation ko linear banata hai; $E_1,E_2$ do points par iske values hain, aur $x=\sin^2\!\big(\tfrac{E_2-E_1}{4}\big)$ measure karta hai ki $E$ kitna badla. **The series** (elliptic sector area ko expand karne se): $$X(x)=\frac{4}{3}\left(1+\frac{6}{5}x+\frac{6\cdot 8}{5\cdot 7}x^2+\frac{6\cdot 8\cdot 10}{5\cdot 7\cdot 9}x^3+\cdots\right).$$ **[[Lagrange f and g functions]]** = closed-form coefficients jo ek Kepler orbit propagate karte hain, $\mathbf r_2=f\,\mathbf r_1+g\,\mathbf v_1$ satisfy karte hue: $$f = 1-\frac{r_2}{p}(1-\cos\theta),\quad g=\frac{r_1 r_2 \sin\theta}{\sqrt{\mu p}},\quad \dot g = 1-\frac{r_1}{p}(1-\cos\theta).$$ Inki key property: kyunki wo relation $\mathbf v_1$ mein *linear* hai, ise invert karne se Kepler dobara solve kiye bina velocity mil jaati hai.

Do pictures jinpar poora method tika hai — swept sector versus flat triangle, aur guess-and-correct loop:

Figure — Gauss's method for Lambert's problem

True or false — justify karo

r_{1}

aur

r_{2}

ko given time mein connect karne wala orbit unique hota hai.

False as stated — unique tabhi hoga jab "short way vs long way" branch aur full revolutions ki count bhi fix karo; geometry +

Δ t

akele ek directional ambiguity chhodh dete hain.

Agar transfer almost flat hai (chhota

θ

), to ratio

y

1

ke kareeb hoga.

True — ek patla sector aur uska triangle almost same area enclose karte hain, isliye sector/triangle

\approx 1

; yahi reason hai ki iteration

y = 1

se start hoti hai.

r_{1}

aur

r_{2}

ke beech chord

c

is baat par depend karta hai ki unhe kaun sa orbit connect karta hai.

False —

c^{2} = r_{1}^{2} + r_{2}^{2} - 2 r_{1} r_{2} cos θ

do points ki pure geometry hai; koi orbit choose karne se pehle ye fixed hota hai.

Δ t

constant

m

mein linearly enter karta hai.

False — ye

Δ t^{2}

ke roop mein enter karta hai, jo Kepler's third law (time-squared scaling) ki echo hai; ise linearly treat karne se galat

m

aur bekaari

y

milti hai.

Gauss ki do equations redundant hain — koi ek bhi kaam kar leti.

False — ek geometry ko auxiliary variable se link karti hai, doosri time ko Kepler's equation ke zariye link karti hai; sahi orbit wahan milti hai jahan dono ek saath hold karein, isliye dono zaroori hain.

Auxiliary variable

x = sin^{2} (\frac{E _{2} - E _{1}}{4})

kabhi negative nahi ho sakta.

False — hyperbolic ya fast transfers ke liye "eccentric anomaly" analogue

x < 0

banata hai (jaise worked example mein

x \approx - 0.01

); series

X (x)

wahan bhi valid rehti hai.

Gauss's method

17 9^{\circ}

transfer ke liye utna hi acha kaam karta hai jitna

4 0^{\circ}

wale ke liye.

False — jaise

θ \to 18 0^{\circ}

cos (θ /2) \to 0

, jo

ℓ

aur

m

ke denominators mein hai, isliye ye blow up karte hain aur iteration fail hoti hai; near-antipodal transfers ke liye Battin/Izzo use karo.

Jab

y

converge ho jaata hai, tab bhi velocities ke liye alag se Kepler's equation solve karni padti hai.

False — Lagrange f and g functions

p

v_{1}, v_{2}

closed form mein dete hain;

f, g

use karne ka yahi poora point hai.

Galti dhundho

θ = arccos \frac{r _{1} \cdot r _{2}}{r _{1} r _{2}}

transfer angle deta hai, ho gaya."

Galti ye hai:

arccos

sirf

0

–

18 0^{\circ}

return karta hai, isliye ye short way aur long way (

θ

36 0^{\circ} - θ

) mein fark nahi kar sakta; branch choose karne ke liye orbit-normal / prograde flag use karna padega.

"Kyunki

m

time carry karta hai aur

Δ t

bada hai, main

Δ t

double kar ke

m

double kar lunga."

Galat:

m \propto Δ t^{2}

, isliye

Δ t

double karne se

m

double nahi, balki chauguna ho jaata hai.

y

triangle area ka sector area se ratio hai."

Ulta hai — ye sector ÷ triangle hai; curved sector bada region hota hai, isliye elliptic short transfers ke liye

y \geq 1

hota hai.

"Main iteration

y = 0

se start karunga safe rehne ke liye."

Ye kharab seed hai:

y = 0

x = m / y^{2} - ℓ

ko diverge karta hai (division by zero) aur koi physical zero-area sector nahi hota; sahi, contraction-friendly seed

y = 1

hai.

θ > 18 0^{\circ}

ke liye main

ℓ

aur

m

ke signs same rakhta hun."

Galti:

cos (θ /2)

18 0^{\circ}

ke baad negative ho jaata hai,

ℓ

aur

m

ke andar ke signs flip kar deta hai; us sign ko track karna zaroori hai, ye assume mat karo ki wo positive rehta hai.

"Lagrange coefficient

g = \frac{r _{1} r _{2} sin θ}{μ p}

hai, isliye agar

θ = 0

hai to bhi mujhe

v_{1} = (r_{2} - f r_{1}) / g

se finite

v_{1}

milega."

Galti:

θ = 0

par,

g \to 0

aur

r_{2} ∥ r_{1}

, isliye division singular hai — do points directionally coincide karte hain aur problem degenerate ho jaati hai.

p

(semi-latus rectum) ko iterate kiye bina seedha geometry se padha ja sakta hai."

Galat:

p

ke formula mein

x

hai, jo sirf

y

converge hone ke baad pata chalta hai; geometry akele

p

under-determine karta hai kyunki wo shape fix karta hai lekin energy/time nahi.

Why questions

Sector ko triangle se kyun compare karte hain na ki poori ellipse se?

Triangle sirf given data se computable hai (

\frac{1}{2} r_{1} r_{2} sin θ

), ek fixed known reference deta hai; sector unknown swept area hai, isliye unka ratio exactly wahi cheez isolate karta hai jo time-of-flight constrain karta hai.

Kepler's second law Gauss ki derivation mein aati hi kyun hai?

Kepler's second law kehta hai equal areas equal times mein sweep hoti hain (

Δ t = 2 (sector area) / h

, jahan

h = μ p

steady area-rate hai), jo ek time input ko ek area condition mein convert karta hai — clock aur geometry ke beech ka bridge.

Typical transfers ke liye iteration itni fast contraction kyun hoti hai?

Ordinary short-way arcs ke liye

y

1

ke bahut kareeb hota hai aur update map ka slope wahan magnitude mein

1

se kafi neeche hota hai, isliye har step error ko bade factor se shrink karta hai — 3–5 steps mein

1 0^{- 8}

aa jaata hai.

m

\propto Δ t^{2}

divided by a length cubed ke roop mein kyun express karte hain?

Ye Kepler's-third-law scaling hai (period

^{2} \propto

size

^{3}

); time ko is tarah package karne se

m

dimensionless ho jaata hai aur poori time dependence ek number mein aa jaati hai.

Velocities get karne ke liye Lagrange f and g functions use karte hain instead of orbit dobara solve karne ke?

Kyunki

r_{2} = f r_{1} + g v_{1}

ek linear relation hai, isliye ise invert karne se

v_{1}

algebraically recover ho jaata hai — Kepler's equation ka second pass nahi chahiye.

Series

X (x)

mein wo specific fractions

\frac{6}{5}, \frac{6 \cdot 8}{5 \cdot 7}, \dots

kyun hain?

Ye elliptic sector-area integral ke

x

ki powers mein term-by-term expansion se aate hain; har fraction us convergent series ka agla coefficient hai, valid kyunki

∣ x ∣

chhota rehta hai.

Gauss's method ek Porkchop plot ka engine kyun hai?

Ek porkchop plot hazaron departure/arrival date pairs sweep karta hai, aur har pair ek Lambert solve hai; Gauss jaisi fast, few-iteration method us poore grid ko scan karna sasta banati hai.

Edge cases

Jaise

θ \to 18 0^{\circ}

ℓ

aur

m

ka kya hota hai?

Dono diverge karte hain kyunki

cos (θ /2) \to 0

unke denominators mein hai; method accuracy khota hai aur converge fail karta hai, isliye Universal variable formulation solver use karo.

y = 1

exactly physically kya represent karta hai?

Sector area triangle area ke barabar hai — ek degenerate flat (straight-line) transfer jahan swept region triangle par collapse ho jaata hai; real ellipses mein

y

thoda

1

se upar hota hai.

Agar

r_{1}

aur

r_{2}

bilkul same direction mein point karein (

θ = 0

) to kya hoga?

Chord aur

sin θ

vanish ho jaate hain,

g \to 0

, aur "orbit" undetermined hai — endpoints radially aligned hain, ek degenerate case jo Gauss ke formulas resolve nahi kar sakte.

Agar dono radii equal hों,

r_{1} = r_{2}

Bilkul theek hai — geometry constant

ℓ

simplify ho jaata hai (

r_{1} + r_{2} = 2 r_{1}

) aur method normally chalta hai; equal radii ka matlab sirf symmetric transfer hai, koi singularity nahi.

Agar required

Δ t

bahut bada ho (multiple revolutions)?

Basic Gauss single arc assume karta hai; multi-revolution transfers extra solution branches introduce karte hain, aur single-branch iteration unhe nahi dhundh payegi — har revolution count ko alag handle karna padega, jaise Izzo Lambert solver karta hai.

Hyperbolic (bahut fast) transfer ke liye method kaisa behave karta hai?

Eccentric-anomaly analogue aur

x

negative ho jaate hain, lekin series

X (x)

chhote

∣ x ∣

ke liye abhi bhi converge karti hai, isliye Gauss mildly hyperbolic cases handle kar sakta hai; strongly hyperbolic ones ke liye Universal variable formulation better hai.

Agar pehle iteration ka

x

negative aaye, jaise worked example mein?

Ye normal hai aur koi error nahi hai — iska matlab sirf ye hai ki

E_{2} - E_{1}

chhota hai ya arc "doosri taraf" bend karta hai; iterate karte raho aur

y

phir bhi settle ho jaata hai (wahan ye

1.027

ke kareeb converge hua).

Recall Quick self-check

Connecting orbit unique hone ke liye kaun si do cheezein fix karni padti hain? ::: Transfer direction (short/long way) aur revolution count.
Gauss's method kahan break karta hai, aur kya replace karta hai ise? ::: $θ = 18 0^{\circ}$ ke kareeb; Battin's method ya Izzo Lambert solver use karo.
Kaun sa ek physical law time input ko area condition mein convert karta hai? ::: Kepler's second law.