300 likes | 390 Views
ASEN 5070: Statistical Orbit Determination I Fall 2013 Professor Brandon A. Jones Professor George H. Born Lecture 7: Generating the STM. Announcements. Lecture Quiz – Due Wednesday Homework 2 – Due September 13 Bring pen and paper on Friday!. Today’s Lecture. Previous Lecture
E N D
ASEN 5070: Statistical Orbit Determination I Fall 2013 Professor Brandon A. Jones Professor George H. Born Lecture 7: Generating the STM
Announcements • Lecture Quiz – Due Wednesday • Homework 2 – Due September 13 • Bring pen and paper on Friday!
Today’s Lecture • Previous Lecture • Linearization – State Transition Matrix (STM) • Alternative Derivation • Solution Methods
Least Squares: Derivation of Normal Form For now, assume G is a linear function • Differentiate with respect to x. What is the answer?
State Transition Matrix • Since x is linear (note lower case!) then there exists a solution to the linear, first order system of differential equations: • The solution is of the form: • Φ(t,ti) is the state transition matrix (STM) that maps x(ti) to the state x(t) at time t.
STM Differential Equation Constant! • What is the differential equation?
Methods to Generate the STM • There are four methods to generate the STM: • Solve from the direct equation (next lecture) • If A is constant, use the Laplace Transform or eigenvector/value analysis • Analytically integrate the differential equation directly • Numerically integrate the equations (ode45)
STM – Alternative Derivation • Expand X(t) in a Taylor series about X*(t):
Laplace Transforms • Laplace Transforms are useful for analysis of linear time-invariant systems: • electrical circuits, • harmonic oscillators, • optical devices, • mechanical systems, • even orbit problems. • Transformation from the time domain into the frequency domain. • Inverse Laplace Transform converts the system back.
Example • Solve the ODE • We can solve this using “traditional” calculus:
Example • Solve the ODE • Or, we can solve this using Laplace Transforms:
Applied to Stat OD • Solve the ODE
‘Direct Approach’ to Solving STM • Leverage the differential equation and combine it with analytic diff. eq. methods • Compatible with simple equations, but not with larger estimated state vectors or complicated dynamics
Numeric Integration • For more complicated dynamics, must integrate X*(t) and Φ(t,t0) simultaneously in propagator • Up to n+n2 propagated states • Derivative function must include the evaluation of the [A(t)]* matrix in addition to F(X,t)
Example Implementation • Use the MATLAB reshape() command to turn matrix into a vector • v = reshape( V, nrows*ncols, 1 ); • MATLAB Demo…