Satellite Orbit and Ephemeris Determination using Inter Satellite Links
Satellite Orbit and Ephemeris Determination using Inter Satellite Links Satellite Orbit and Ephemeris Determination using Inter Satellite Links
Software DescriptionInter Satellite Linksthan one error has occurred and if at a given point of time m (m > 4 for FD and m > 5 for FDI)range type measurements are available, then linearization yields the linear modely = Gx + εEq. 5.6-1But what, if no measurement is faulty but the ephemeris or clock state? To allow the isolationof the satellites own faulty clock or ephemeris, these parameters are introduced as pseudoobservations and the above equation replaced bywithy = G * ⋅x+ ε ∗* Eq. 5.6-2and⎛ 0 ⎞⎜ ⎟⎜ 0 ⎟⎜ ⎟y * =⎜ ⎟⎜0 ,⎟⎜ ⎟⎜ 0 ⎟⎜ ⎟⎝yM ⎠T Te , R ATTCT⎛T⎜e R⎜⎜e⎜G*=T⎜e⎜T⎜ 0⎜⎝ GTATCT0⎞⎟⎟0⎟⎟0⎟⎟1⎟⎟g⎠and ε∗ =⎛ ε⎜⎜ε⎜⎜⎜ε⎜⎜ ε⎜⎝ εe , e unit vector in radial, along track and cross track directionRATCTCM⎞⎟⎟⎟⎟⎟⎟⎟⎟⎠Eq. 5.6-3The residuals are given as zero, i.e. "no ephemeris fault" and "no clock fault". The RAIMalgorithm is now capable of removing the bad assumption of "no radial error" for example, ifthe removal of this row in the system of observation equations minimises the sample varianceof residuals.A major draw back of that kind of snapshot algorithm based on the sample variance is theneed for a sophisticated pre-processing of the raw data. The sample variance taken from theraw measurements is still too noisy, thus leading to lots of false alarms.Another possible way of monitoring the integrity of a satellites position and clock is byseparating satellite dynamics / errors / and observation noise by their dynamic behaviour.This can be achieved using the Kalman filter with the following state vector⎡ ∆x⎤⎢ ⎥⎢∆y⎥⎢ ∆z⎥⎢ ⎥⎢∆Tx = ⎥⎢ ∆x⎥⎢ ⎥⎢ ∆y⎥⎢ ∆z⎥⎢ ⎥⎢⎣∆T⎥⎦Eq. 5.6-4which is kept very adaptive by adding high process noise. The Kalman filter is then inert withrespect to noise, but reacts immediately on real errors. The dynamic model which has beenPage 90R. Wolf
Inter Satellite LinksSoftware Descriptionimplemented and used successfully assumes the errors in x-y-z direction (ECEF referenceframe) as well as the satellite clock to be composed of a step and a ramp, expressed by to thefollowing transition matrix:⎡1⎢⎢0⎢⎢0⎢⎢0Φ = ⎢⎢0⎢⎢0⎢⎢⎢0⎢⎣00100000000100000000100001000100001000100001000100⎤⎥0⎥⎥0⎥⎥1⎥⎥0⎥⎥0⎥⎥⎥0⎥⎥1⎦The observation matrix is given by⎡1⎢⎢0H = ⎢⎢0⎢⎢⎣00100001000010000000000000⎤⎥0⎥⎥0⎥⎥0⎥⎦Note that this dynamic system is very close to the one used to estimate the orbit corrections,but much less smoothing character. Moreover, the same filter tuning cannot be used forintegrity monitoring and orbit estimation. This approach is more suited for onboardprocessing and therefore elaborated in more detail in chapter 7.R. Wolf Page 91
- Page 56 and 57: Orbit ComputationInter Satellite Li
- Page 58 and 59: Orbit ComputationInter Satellite Li
- Page 60 and 61: Orbit ComputationInter Satellite Li
- Page 62 and 63: Orbit ComputationInter Satellite Li
- Page 64 and 65: Orbit ComputationInter Satellite Li
- Page 66 and 67: Orbit ComputationInter Satellite Li
- Page 68 and 69: Orbit ComputationInter Satellite Li
- Page 70 and 71: Orbit ComputationInter Satellite Li
- Page 72 and 73: Orbit ComputationInter Satellite Li
- Page 74 and 75: Orbit ComputationInter Satellite Li
- Page 76 and 77: Orbit ComputationInter Satellite Li
- Page 78 and 79: Orbit ComputationInter Satellite Li
- Page 80 and 81: Orbit ComputationInter Satellite Li
- Page 82 and 83: Orbit ComputationInter Satellite Li
- Page 84 and 85: Orbit ComputationInter Satellite Li
- Page 86 and 87: Software DescriptionInter Satellite
- Page 88 and 89: Software DescriptionInter Satellite
- Page 90 and 91: Software DescriptionInter Satellite
- Page 92 and 93: Software DescriptionInter Satellite
- Page 94 and 95: Software DescriptionInter Satellite
- Page 96 and 97: Software DescriptionInter Satellite
- Page 98 and 99: Software DescriptionInter Satellite
- Page 100 and 101: Software DescriptionInter Satellite
- Page 102 and 103: Software DescriptionInter Satellite
- Page 104 and 105: Software DescriptionInter Satellite
- Page 108 and 109: Simulations and ResultsInter Satell
- Page 110 and 111: Simulations and ResultsInter Satell
- Page 112 and 113: Simulations and ResultsInter Satell
- Page 114 and 115: Simulations and ResultsInter Satell
- Page 116 and 117: Simulations and ResultsInter Satell
- Page 118 and 119: Simulations and ResultsInter Satell
- Page 120 and 121: Simulations and ResultsInter Satell
- Page 122 and 123: Simulations and ResultsInter Satell
- Page 124 and 125: Simulations and ResultsInter Satell
- Page 126 and 127: Simulations and ResultsInter Satell
- Page 128 and 129: Simulations and ResultsInter Satell
- Page 130 and 131: Simulations and ResultsInter Satell
- Page 132 and 133: Simulations and ResultsInter Satell
- Page 134 and 135: Simulations and ResultsInter Satell
- Page 136 and 137: Simulations and ResultsInter Satell
- Page 138 and 139: Simulations and ResultsInter Satell
- Page 140 and 141: Simulations and ResultsInter Satell
- Page 142 and 143: Simulations and ResultsInter Satell
- Page 144 and 145: Simulations and ResultsInter Satell
- Page 146 and 147: Simulations and ResultsInter Satell
- Page 148 and 149: Simulations and ResultsInter Satell
- Page 150 and 151: Simulations and ResultsInter Satell
- Page 152 and 153: Simulations and ResultsInter Satell
- Page 154 and 155: Simulations and ResultsInter Satell
Software Description<strong>Inter</strong> <strong>Satellite</strong> <strong>Links</strong>than one error has occurred <strong>and</strong> if at a given point of time m (m > 4 for FD <strong>and</strong> m > 5 for FDI)range type measurements are available, then linearization yields the linear modely = Gx + εEq. 5.6-1But what, if no measurement is faulty but the ephemeris or clock state? To allow the isolationof the satellites own faulty clock or ephemeris, these parameters are introduced as pseudoobservations <strong>and</strong> the above equation replaced bywithy = G * ⋅x+ ε ∗* Eq. 5.6-2<strong>and</strong>⎛ 0 ⎞⎜ ⎟⎜ 0 ⎟⎜ ⎟y * =⎜ ⎟⎜0 ,⎟⎜ ⎟⎜ 0 ⎟⎜ ⎟⎝yM ⎠T Te , R ATTCT⎛T⎜e R⎜⎜e⎜G*=T⎜e⎜T⎜ 0⎜⎝ GTATCT0⎞⎟⎟0⎟⎟0⎟⎟1⎟⎟g⎠<strong>and</strong> ε∗ =⎛ ε⎜⎜ε⎜⎜⎜ε⎜⎜ ε⎜⎝ εe , e unit vector in radial, along track <strong>and</strong> cross track directionRATCTCM⎞⎟⎟⎟⎟⎟⎟⎟⎟⎠Eq. 5.6-3The residuals are given as zero, i.e. "no ephemeris fault" <strong>and</strong> "no clock fault". The RAIMalgorithm is now capable of removing the bad assumption of "no radial error" for example, ifthe removal of this row in the system of observation equations minimises the sample varianceof residuals.A major draw back of that kind of snapshot algorithm based on the sample variance is theneed for a sophisticated pre-processing of the raw data. The sample variance taken from theraw measurements is still too noisy, thus leading to lots of false alarms.Another possible way of monitoring the integrity of a satellites position <strong>and</strong> clock is byseparating satellite dynamics / errors / <strong>and</strong> observation noise by their dynamic behaviour.This can be achieved <strong>using</strong> the Kalman filter with the following state vector⎡ ∆x⎤⎢ ⎥⎢∆y⎥⎢ ∆z⎥⎢ ⎥⎢∆Tx = ⎥⎢ ∆x⎥⎢ ⎥⎢ ∆y⎥⎢ ∆z⎥⎢ ⎥⎢⎣∆T⎥⎦Eq. 5.6-4which is kept very adaptive by adding high process noise. The Kalman filter is then inert withrespect to noise, but reacts immediately on real errors. The dynamic model which has beenPage 90R. Wolf