Title: AA 4362 Astrodynamics
1AA 4362Astrodynamics
Solving Kepler's (Barkers) Equation for
parabolic and hyperbolic trajectories
Week 4 Vallado, Chapter 2, sections 2.2, 2.3, 2.4
2Homework 8
Extend Your Kepler Propagation and Solver
Algorithm To Be smart enough to solve all
trajectories. Elliptical, Parabolic, Hyperbolic
Allow For arbitrary M, e, and Starting value
for E Evaluate Algorithm performance for
e0. 5, e1.0, e1.5, e5 Use Methods, 1,
2, 3 as startup algorithms Count number of
iterations for convergence Draw Convergence
Plots for e1, e5 cases
3Generalized Formsfor Keplers Equation
Ellipse Parabola ??????????
4Ellipse
5Ellipse (concluded)
6Parabola
7Parabola (continued)
True Mean Anomaly
Iteration Index
Current estimate
Refined estimate
Current estimate of mean anomaly
Steepest Descent
8Parabola Convergence plot(near perigee)
9Parabola Convergence plot(far from perigee)
10Parabolic Solver, Startup Condition
Approximate Solution as n --gt p
n
11Parabolic Solver, Startup Condition
Example M 10.0
12Parabolic Solver, Startup Condition
Example M 0.1
13Hyperbola
14Hyperbola (concluded)
True Mean Anomaly
Iteration Index
Current estimate
Refined estimate
Current estimate of mean anomaly
Steepest Descent
15Hyperbola, Numerical example
Example
n, Approx 115?
16Hyperbola, convergence plot
true anomaly near perigee
n, Approx 115?
17Hyperbola, convergence plot
18Hyperbola, convergence plot
true anomaly far from perigee
n, Approx 131.8?
19Hyperbola, convergence plot
20Soooo we need a reliablestartup algorithm
So where do we start? with a really Simple
ad-hoc solution .. When M 0, we simply
kick it off zero by adding or subtracting e
(Vallado Algorithm)
F0
F0
21Ad-Hoc Start Up Algorithm
Me 16.148
22Homework 9
Develop a Higher order Taylor Series
method Similar to that Developed for the
Elliptical Kepler Startup Hint
0
0
23Homework 9 (contd)
Find single real root of resolvent cubic
24Homework 9 (contd)
Real Root
Compare convergence properties of this startup
method Versus Ad-Hoc startup algorithm for a
variety of e, M
F0
F0