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.
- 2. The method of claim 1 further comprising determining an estimate of density {circumflex over (ρ)}(H) in real-time for use in the computation of the acceleration of the satellite due to atmospheric drag during satellite trajectory propagation, wherein the satellite has a current height “H”, and wherein the determination of the estimate of density {circumflex over (ρ)}(H) comprises: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 D^(HP)=Δ ρ^(HP)ρ_(HP);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 D^(HP)=Δ ρ^(HP)ρ_(HP)comprises applying a Gauss-Markov transition function Φ(tk+1, tk) to propagate {circumflex over (D)}(HP) from time tk to time tk+1 where:{circumflex over (D)}k+1(HP)=Φ(tk+1, tk){circumflex over (D)}k(HP); ((tk+1−tk)→0)(Φ(tk+1, tk)→1) and ((tk+1−tk)→∞)(Φ(tk+1, tk)→0).
- 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.
CROSS REFERENCE TO RELATED APPLICATIONS
This application claims priority under 35 U.S.C. § 119(e) from provisional application No. 60/337,809, filed Nov. 13, 2001. The 60/337,809 provisional application is incorporated by reference herein, in its entirety, for all purposes.
US Referenced Citations (3)
Number |
Name |
Date |
Kind |
6515617 |
Demers et al. |
Feb 2003 |
B1 |
6567712 |
Rog et al. |
May 2003 |
B1 |
6570530 |
Gaal et al. |
May 2003 |
B2 |
Provisional Applications (1)
|
Number |
Date |
Country |
|
60/337809 |
Nov 2001 |
US |