Method and apparatus for orbit determination
First Claim
1. A method of estimating in real-time an orbit of a satellite using a recursive filter and wherein the recursive filter first receives an a priori state estimate and an a priori state error covariance matrix for the satellite, and wherein the method comprises:
- inputting tracking measurements for the satellite received in real-time from tracking stations to the recursive filter at filter time update “
T”
;
inputting filter tracking station platform data to the recursive at filter time update “
T”
;
computing recursively the current state estimate, the current error covariance and the current filter measurement update of the satellite at filter time update “
T”
;
moving the state estimate and error covariance from filter time update “
T”
to filter time update “
T+1”
;
incorporating into the filter measurement update the information provided by the tracking data measurement at filter time update T; and
determining a current estimate of the orbit of the satellite.
8 Assignments
0 Petitions
Accused Products
Abstract
A method and apparatus for estimating the orbits of spacecraft or other objects relative to primary celestial bodies, given applicable measurements. The present invention uses the variation of parameters in universal variables for spacecraft trajectory propagation during operation of its optimal sequential orbit estimation. This method provides a universal capability for all three conics: elliptic, parabolic, and hyperbolic, forward-time and backward-time numerical integration for filtering and smoothing, respectively. A method for determining local and global estimates of atmospheric density is also provided. A method for realistic orbit simulation provides a means for validating an orbit estimate.
13 Citations
11 Claims
-
1. A method of estimating in real-time an orbit of a satellite using a recursive filter and wherein the recursive filter first receives an a priori state estimate and an a priori state error covariance matrix for the satellite, and wherein the method comprises:
-
inputting tracking measurements for the satellite received in real-time from tracking stations to the recursive filter at filter time update “
T”
;
inputting filter tracking station platform data to the recursive at filter time update “
T”
;
computing recursively the current state estimate, the current error covariance and the current filter measurement update of the satellite at filter time update “
T”
;
moving the state estimate and error covariance from filter time update “
T”
to filter time update “
T+1”
;
incorporating into the filter measurement update the information provided by the tracking data measurement at filter time update T; and
determining a current estimate of the orbit of the satellite. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11)
computing from an atmospheric density model (ADM) atmospheric density {overscore (ρ
)} at height H ({overscore (ρ
)}(H)), wherein the ADM derives {overscore (ρ
)}(H) in part from current observatory measured solar and geomagnetic indices;
propagating through time an estimate of the relative density correction at the perigee height {circumflex over (D)}(HP) of the satellite, wherein HP is the satellite height at perigee and deweighting the variance estimate of {circumflex over (D)}(HP); computing an estimate of the correction to atmospheric density at a perigee height, Δ
{circumflex over (ρ
)}(HP);
computing from Δ
{circumflex over (ρ
)}(HP) and the current height of the satellite “
H”
the value of a mapping function, SΔ
ρ
;
computing an estimate of the density correction at the current height, Δ
{circumflex over (ρ
)}(H), by multiplying SΔ
ρ
by {circumflex over (D)}(HP); and
adding the density correction estimate Δ
{circumflex over (ρ
)}(H) to {overscore (ρ
)}(H) to determine {circumflex over (ρ
)}(H).
-
-
3. The method of claim 2, wherein propagating through time an estimate of the relative density correction at the perigee height {circumflex over (D)}(HP) of the satellite, wherein HP is the satellite height at perigee and
-
( H P ) = Δ ρ ^ ( H P ) ρ _ ( H P ) comprises applying a Gauss-Markov transition function Φ
(tk+1, tk) to propagate {circumflex over (D)}(HP) from time tk to time tk+1 where;
-
-
4. The method of claim 2, wherein deweighting the variance estimate of {circumflex over (D)}(HP) comprises applying a Gauss-Markov deweighting function for variance on the filter time update where:
-
{overscore (σ
)}k+12({circumflex over (D)}k(HP))=└
1−
Φ
2(tk+1, tk)┘
{overscore (σ
)}Δ
ρ
2(HP)where {overscore (σ
)}Δ
ρ
2(HP) is a positive time constant.
-
-
5. The method of claim 2, wherein the current observatory measured solar and geomagnetic indices comprises the ten centimeter solar flux measurement F10, its mean across several months {overscore (F)}10, and a geomagnetic measurement.
-
6. The method of claim 5 wherein geomagnetic measurement is selected from the group consisting of the three-hourly geomagnetic measurement KP and aP.
-
7. The method of claim 6 wherein the recursive filter gain is determined, in part, by the value of KP.
-
8. The method of claim 6 wherein the recursive filter gain is determined, in part, by the value of aP.
-
9. The method of claim 1 further comprising validating an orbit determination using simulated data wherein the simulated data is generated by;
-
obtaining a priori state of the satellite;
computing an initial simulated orbit estimate and a simulated truth orbit;
determining the initial orbit error by computing the difference between simulated truth orbit and the simulated orbit estimate; and
adjusting the initial orbit error using error estimates reflected in the state error covariance matrix.
-
-
10. The method of claim 9 wherein the error estimates reflected in the state error covariance matrix are selected from the group consisting of serially correlated transponder biases for two-way measurements, serially correlated measurement biases for two-way measurements, serially correlated clock biases for one-way measurements;
- and serially correlated force model error biases for all orbits.
-
11. The method of claim 1 wherein computing recursively at filter time update “
- T”
the current state estimate, the current error covariance and the current filter measurement update of the satellite, moving the state estimate and error covariance from filter time update “
T”
to filter time update “
T+1”
, incorporating into the filter measurement update the information provided by the tracking data measurement at filter time update T, and arriving at a current estimate of the orbit of the satellite are performed in an object-oriented computer program.
- T”
Specification