350 likes | 460 Views
ASEN 5070: Statistical Orbit Determination I Fall 2014 Professor Brandon A. Jones Lecture 5: State Deviations and Fundamentals of Linear Algebra. Announcements. Homework 2– Due September 12 Make-up Lecture Today @ 3pm, here. Today’s Lecture. Effects of State Deviations
E N D
ASEN 5070: Statistical Orbit Determination I Fall 2014 Professor Brandon A. Jones Lecture 5: State Deviations and Fundamentals of Linear Algebra
Announcements • Homework 2– Due September 12 • Make-up Lecture • Today @ 3pm, here
Today’s Lecture • Effects of State Deviations • Linear Algebra (Appendix B)
Effects of Small Variations • Quantification of such effects is fundamental to the OD methods discussed in this course! Time
Effects of Small Variations • Let’s think about the effects of small variations in coordinates, and how these impact future states. Example: Propagating a state in the presence of NO forces Final State: (xf, yf, zf, vxf, vyf, vzf) Initial State: (x0, y0, z0, vx0, vy0, vz0)
Effects of Small Variations • What happens if we perturb the value of x0? Force model: 0 Initial State: (x0+Δx, y0, z0, vx0, vy0, vz0) Final State: (xf, yf, zf, vxf, vyf, vzf) Initial State: (x0, y0, z0, vx0, vy0, vz0)
Effects of Small Variations • What happens if we perturb the value of x0? Force model: 0 Final State: (xf+Δx, yf, zf, vxf, vyf, vzf) Initial State: (x0+Δx, y0, z0, vx0, vy0, vz0) Final State: (xf, yf, zf, vxf, vyf, vzf) Initial State: (x0, y0, z0, vx0, vy0, vz0)
Effects of Small Variations • What happens if we perturb the position? Force model: 0 Final State: (xf+Δx, yf+Δy, zf+Δz, vxf, vyf, vzf) Initial State: (x0+Δx, y0+Δy, z0+Δz, vx0, vy0, vz0) Initial State: (x0, y0, z0, vx0, vy0, vz0)
Effects of Small Variations • What happens if we perturb the value of vx0? Force model: 0 Initial State: (x0, y0, z0, vx0-Δvx, vy0, vz0) Final State: (xf, yf, zf, vxf, vyf, vzf) Initial State: (x0, y0, z0, vx0, vy0, vz0)
Effects of Small Variations • What happens if we perturb the value of vx0? Force model: 0 Final State: (xf+tΔvx, yf, zf, vxf+Δvx, vyf, vzf) Initial State: (x0, y0, z0, vx0+Δvx, vy0, vz0) Final State: (xf, yf, zf, vxf, vyf, vzf) Initial State: (x0, y0, z0, vx0, vy0, vz0)
Effects of Small Variations • What happens if we perturb the position and velocity? Force model: 0
Effects of Small Variations • We could have arrived at this easily enough from the equations of motion. Force model: 0
Effects of Small Variations • This becomes more challenging with nonlinear dynamics Force model: two-body
Effects of Small Variations • This becomes more challenging with nonlinear dynamics Final State: (xf, yf, zf, vxf, vyf, vzf) Force model: two-body Initial State: (x0, y0, z0, vx0, vy0, vz0) The partial of one Cartesian parameter wrt the partial of another Cartesian parameter is ugly.
Effects of Small Variations • This becomes more challenging with nonlinear dynamics Final State: (xf, yf, zf, vxf, vyf, vzf) Force model: two-body
Matrix Basics • Matrix A is comprised of elements ai,j • The matrix transpose swaps the indices
Matrix Basics • Matrix inverse A-1 is the matrix such that • For the inverse to exist, A must be square • We will treat vectors as n×1 matrices
2x2 Matrix Inverse Trick • If we have a 2x2, nonsingular matrix: Asking you to invert a full 2x2 matrix on an exam is fair game!
Matrix Determinant • The square matrix determinant, |A|, describes if a solution to a linear system exists: • It also describes the change in area/volume/etc. due to a linear operation:
Linear Independence • A set of vectors are linearly independent if none of them can be expressed as a linear combination of other vectors in the set • In other words, no scalars αi exist such that for some vector vj in the set {vi}, i=1,…,n,
Matrix Rank • The matrix column rank is the number of linearly independent columns of a matrix • The matrix row rank is the number of linearly independent rows of a matrix • rank(A) = min( col. rank of A, row rank of A)
Vector Differentiation • When differentiating a scalar function w.r.t. a vector:
Vector Differentiation • When differentiating a function with vector output w.r.t. a vector:
Matrix Derivative Identities • If A and B are n×1 vectors that are functions of X:
Positive Definite Matrices • The n×n matrix A is positive definite if and only if: • The n×n matrix A is positive semi-definite if and only if:
Minimum of a function • The point x is a minimum if and is positive definite.
Eigenvalues/vectors • Given the n×n matrix A, there are n eigenvalues λ and vectors X≠0 where
Book Appendix B • Other identities/definitions in Appendix B of the book • Matrix Trace • Maximum/Minimum Properties • Matrix Inversion Theorems • Review the appendix and make sure you understand the material