Helio Koiti Kuga

Single frequency GPS measurements in real-time artificial satellite orbit determination ☆

Abstract A simplified and compact algorithm with low computational cost providing an accuracy around tens of meters for artificial satellite orbit determination in real-time and on-board is developed in this work. The state estimation method is the extended Kalman filter. The Cowells method is used to propagate the state vector, through a simple Runge–Kutta numerical integrator of fourth order with fixed step size. The modeled forces are due to the geopotential up to 50th order and degree of JGM-2 model. To time-update the state error covariance matrix, it is considered a simplified force model. In other words, in computing the state transition matrix, the effect of J2 (Earth flattening) is analytically considered, which unloads dramatically the processing time. In the measurement model, the single frequency GPS pseudorange is used, considering the effects of the ionospheric delay, clock offsets of the GPS and user satellites, and relativistic effects. To validate this model, real live data are used from Topex/Poseidon satellite and the results are compared with the Topex/Poseidon Precision Orbit Ephemeris (POE) generated by NASA/JPL, for several test cases. It is concluded that this compact algorithm enables accuracies of tens of meters with such simplified force model, analytical approach for computing the transition matrix, and a cheap GPS receiver providing single frequency pseudorange measurements.

Unscented Kalman Filter Applied to the Spacecraft Attitude Estimation with Euler Angles

The aim of this work is to test an algorithm to estimate, in real time, the attitude of an artificial satellite using real data supplied by attitude sensors that are on board of the CBERS-2 satellite (China Brazil Earth Resources Satellite). The real-time estimator used in this work for attitude determination is the Unscented Kalman Filter. This filter is a new alternative to the extended Kalman filter usually applied to the estimation and control problems of attitude and orbit. This algorithm is capable of carrying out estimation of the states of nonlinear systems, without the necessity of linearization of the nonlinear functions present in the model. This estimation is possible due to a transformation that generates a set of vectors that, suffering a nonlinear transformation, preserves the same mean and covariance of the random variables before the transformation. The performance will be evaluated and analyzed through the comparison between the Unscented Kalman filter and the extended Kalman filter results, by using real onboard data.

A Discussion Related to Orbit Determination Using Nonlinear Sigma Point Kalman Filter

Herein, the purpose is to present a Kalman filter based on the sigma point unscented transformation development, aiming at real-time satellite orbit determination using GPS measurements. First, a brief review of the extended Kalman filter will be done. After, the sigma point Kalman filter will be introduced as well as the basic idea of the unscented transformation, in which this filter is based. Following, the unscented Kalman filter applied to orbit determination will be explained. Such explanation encloses formulations about the orbit determination through GPS; the dynamic model; the observation model; the unmodeled acceleration estimation; also an application of this new filter approaches on orbit determination using GPS measurements discussion.

Spin-Stabilized Spacecrafts: Analytical Attitude Propagation Using Magnetic Torques

An analytical approach for spin-stabilized satellites attitude propagation is presented, considering the influence of the residual magnetic torque and eddy currents torque. It is assumed two approaches to examine the influence of external torques acting during the motion of the satellite, with the Earths magnetic field described by the quadripole model. In the first approach is included only the residual magnetic torque in the motion equations, with the satellites in circular or elliptical orbit. In the second approach only the eddy currents torque is analyzed, with the satellite in circular orbit. The inclusion of these torques on the dynamic equations of spin stabilized satellites yields the conditions to derive an analytical solution. The solutions show that residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spacecrafts spin axis and the eddy currents torque causes an exponential decay of the angular velocity magnitude. Numerical simulations performed with data of the Brazilian Satellites (SCD1 and SCD2) show the period that analytical solution can be used to the attitude propagation, within the dispersion range of the attitude determination system performance of Satellite Control Center of Brazil National Research Institute.

Onboard and real-time artificial satellite orbit determination using GPS

An algorithm for real-time and onboard orbit determination applying the Extended Kalman Filter (EKF) method is developed. Aiming at a very simple and still fairly accurate orbit determination, an analysis is performed to ascertain an adequacy of modeling complexity versus accuracy. The minimum set of to-be-estimated states to reach the level of accuracy of tens of meters is found to have at least the position, velocity, and user clock offset components. The dynamical model is assessed through several tests, covering force model, numerical integration scheme and step size, and simplified variational equations. The measurement model includes only relevant effects to the order of meters. The EKF method is chosen to be the simplest real-time estimation algorithm with adequate tuning of its parameters. In the developed procedure, the obtained position and velocity errors along a day vary from 15 to 20 m and from 0.014 to 0.018 m/s, respectively, with standard deviation from 6 to 10 m and from 0.006 to 0.008 m/s, respectively, with the SA either on or off. The results, as well as analysis of the final adopted models used, are presented in this work.

Comparison between Two Methods to Calculate the Transition Matrix of Orbit Motion

Two methods to evaluate the state transition matrix are implemented and analyzed to verify the computational cost and the accuracy of both methods. This evaluation represents one of the highest computational costs on the artificial satellite orbit determination task. The first method is an approximation of the Keplerian motion, providing an analytical solution which is then calculated numerically by solving Keplers equation. The second one is a local numerical approximation that includes the effect of 𝐽2. The analysis is performed comparing these two methods with a reference generated by a numerical integrator. For small intervals of time (1 to 10 s) and when one needs more accuracy, it is recommended to use the second method, since the CPU time does not excessively overload the computer during the orbit determination procedure. For larger intervals of time and when one expects more stability on the calculation, it is recommended to use the first method.

Analysis of Attitude Determination Methods Using GPS Carrier Phase Measurements

If three or more GPS antennas are mounted properly on a platform and differences of GPS signals measurements are collected simultaneously, the baselines vectors between antennas can be determined and the platform orientation defined by these vectors can be calculated. Thus, the prerequisite for attitude determination technique based on GPS is to calculate baselines between antennas to millimeter level of accuracy. For accurate attitude solutions to be attained, carrier phase double differences are used as main type of measurements. The use of carrier phase measurements leads to the problem of precise determination of the ambiguous integer number of cycles in the initial carrier phase (integer ambiguity). In this work two algorithms (LSAST and LAMBDA) were implemented and tested for ambiguity resolution allowing accurate real-time attitude determination using measurements given by GPS receivers in coupled form. Platform orientation was obtained using quaternions formulation, and the results showed that LSAST method performance is similar to LAMBDA as far as the number of epochs which are necessary to resolve ambiguities is concerned, but with processing time significantly higher. The final result accuracy was similar for both methods, better than 0.1° to 0.2°, when baselines are considered in decoupled form.

Orbital propagation for Brazilian satellites using NORAD models

Abstract The performance of the NORAD models for near Earth satellites (SGP, SGP4, SGP8) using two Brazilian flying satellites, SCD-1 and SCD-2, and the Chinese-Brazilian satellite CBERS-1 is compared. The achievable accuracy of such models is compared against the predicted 2-lines elements set for the satellites. Every week an updated fresh set of 2-lines elements for these satellites is made available through the Internet. About ten years of observations of the SCD-1 satellite are available and therefore solar activity influences on the orbit can be analyzed. Data are selected considering also orbital (for CBERS-1) and attitude (for SCD-2) maneuvers.

Simulations under Ideal and Nonideal Conditions for Characterization of a Passive Doppler Geographical Location System Using Extension of Data Reception Network

This work presents a (Data Reception Network) DRN software investigation considering simulated conditions inserting purposely errors into the Doppler measurements, satellites ephemeris, and time stamp, to characterize the geographical location software (GEOLOC) developed by Sousa (2000) and Sousa et al. (2003). The extension of reception stations in Brazilian territory can result in more precise locations if the network is considered in the GEOLOC. The results and analyses were first obtained considering the ground stations separately, to characterize their effects in the geographical location (GL) result. Six conditions were investigated: ideal simulated conditions, random and bias errors in the Doppler measurements, errors in the satellite ephemeris, and errors in the time stamp in order to investigate the DRN importance to get more accurate locations; an analysis was performed considering the random errors of 1 Hz in the Doppler measurements. The results are quite satisfactory and also show good compatibility between the simulator and the GEOLOC using the DRN.

Nonrecursive algorithm for remote geolocation using ranging measurements

The objective of this work is to present a new algebraic solution for the problem of remote determination of geographic coordinates of a target, using a new remote geopositioning system being developed in Brazil. It can be useful for double-check measurements obtained with other methods, for certain critical applications, being capable to perform independently from them. This system requires three-reference bases on the surface of the earth with synchronized clocks and a repeater in space. Calculations are derived from measurements of propagation time of clock signal transmitted by one base to all bases and target, via a transponder in space. The algorithm also provides the “instantaneous” determination of the repeater position in space and brings other applications in navigation and remote clock synchronization. The algorithm has been successfully tested through computational software.