Sensitivity of the propagated two-body state to the initial state.

Definition

  • State transition matrix (STM): the Jacobian matrix (matrix of partial derivatives) of the orbital state at the final time with respect to the state at the initial time
  • Notation: initial position and velocity subscript 1, final subscript 2
    • Φ(𝑡2,𝑡1)=[𝐴𝐶𝐵𝐷]=[𝜕𝐫2𝜕𝐫1𝜕̇𝐫2𝜕𝐫1𝜕𝐫2𝜕̇𝐫1𝜕̇𝐫2𝜕̇𝐫1]
    • 𝐴, 𝐵, 𝐶, 𝐷 are 3×3 blocks
  • Propagation itself uses the Lagrange coefficients; the STM needs their derivatives too

Battin’s Method

  • References:
    • Battin’s method
    • Arora, Russell, Strange, Ottesen, “Partial derivatives of the solution to the Lambert boundary value problem”, JGCD, 2015, doi: 10.2514/1.G001030
  • Universal formulation: Battin’s technique formulates the STM blocks 𝐴, 𝐵, 𝐶, 𝐷 in universal variables, so one set of formulas covers every conic

Universal Formulation

  • Inverse semimajor axis from the initial state:
    • 𝛼=2𝑟1−𝑣21𝜇
  • Radial-velocity-like terms at both ends:
    • 𝜎1=𝐫𝖳1𝐯1√𝜇,𝜎2=𝐫𝖳2𝐯2√𝜇
  • Universal anomaly from the time of flight (TOF):
    • 𝜒=𝛼√𝜇TOF+𝜎2−𝜎1
  • Auxiliary quantities 𝑈1,…,𝑈5 (with the Lagrange coefficients 𝑓, 𝑔, ̇𝑓 over the transfer):
    • 𝑈1=−𝑟2𝑟1̇𝑓√𝜇,𝑈2=𝑟1(1−𝑓),𝑈3=√𝜇(TOF−𝑔)
    • 𝑈4=𝑈1𝑈3−12(𝑈22−𝛼𝑈23)
    • 𝑈5=16𝜒3−𝑈3𝛼
  • And
    • 𝐶=1√𝜇(3𝑈5−𝜒𝑈4−√𝜇TOF𝑈2)

STM Blocks

  • 𝐼 is the 3×3 identity matrix; bold 𝐫, ̇𝐫 are vectors, 𝑟1, 𝑟2 magnitudes

Position-to-Position Submatrix

  • 𝐴=𝑟2𝜇(̇𝐫2−̇𝐫1)(̇𝐫2−̇𝐫1)𝖳+1𝑟31[𝑟1(1−𝑓)𝐫2𝐫𝖳1+𝐶̇𝐫2𝐫𝖳1]+𝑓𝐼

Velocity-to-Position Submatrix

  • 𝐵=𝑟1𝜇(1−𝑓)[(𝐫2−𝐫1)̇𝐫𝖳1−(̇𝐫2−̇𝐫1)𝐫𝖳1]+𝐶𝜇̇𝐫2̇𝐫𝖳1+𝑔𝐼

Position-to-Velocity Submatrix

  • 𝐶=−1𝑟21(̇𝐫2−̇𝐫1)𝐫𝖳1−1𝑟22𝐫2(̇𝐫2−̇𝐫1)𝖳+̇𝑓[𝐼−1𝑟22𝐫2𝐫𝖳2+1𝜇𝑟2(𝐫2̇𝐫𝖳2−̇𝐫2𝐫𝖳2)𝐫2(̇𝐫2−̇𝐫1)]−𝜇𝐶𝑟31𝑟32𝐫2𝐫𝖳1

Unclear in notes

  • In 𝐶 the term 𝐫2(̇𝐫2−̇𝐫1) inside the brackets is printed without a transpose, which does not give a 3×3 matrix; presumably it is the outer product 𝐫2(̇𝐫2−̇𝐫1)𝖳. Check against Battin / Arora et al.

Velocity-to-Velocity Submatrix

  • 𝐷=𝑟1𝜇(̇𝐫2−̇𝐫1)(̇𝐫2−̇𝐫1)𝖳+1𝑟32[𝑟1(1−𝑓)𝐫2𝐫𝖳1−𝐶𝐫2̇𝐫𝖳1]+̇𝑔𝐼
  • The identity terms (𝑓𝐼, 𝑔𝐼, ̇𝑓𝐼, ̇𝑔𝐼) are the propagation matrix; the rest accounts for 𝑓, 𝑔 changing with the initial state
  • Other STMs in the course: the HCW STM for relative motion about a circular orbit, and the Yamanaka-Ankersen STM for elliptical reference orbits