Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kalman Filtering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71

TRMM Re-Entry Planning: Attitude Determination and Control During Thruster Modes

The Tropical Rainfall Measuring Mission (TRMM) spacecraft has been undergoing design for a controlled re-entry to Earth. During simulation of the re-entry plan, there was evidence of errors in the attitude determination algorithms during thruster modes. These errors affected the bum efficiency, and thus planning, during re-entry. During thruster modes, the spacecraft attitude is controlled off of integrated Gyro Error Angles that were designed to closely follow the nominal spacecraft pointing frame (Tip Frame). These angles, however, were not exactly mapped to the Tip Frame from the Body Frame. Additionally, in the initial formulation of the thruster mode attitude determination algorithms, several assumptions and approximations were made to conserve processor speed. These errors became noticeable and significant when simulating bums of much longer duration (-10 times) than had been produced in flight. A solution is proposed that uses attitude determination information from a propagated extended Kalman filter that already exists in the TRMM thruster modes. This attitude information is then used to rotate the Gyro Error Angles into the Tip Frame. An error analysis is presented that compares the two formulations. The new algorithm is tested using the TRMM High-Fidelity Simulator and verified with the TRMM Software Testing and Training Facility. Simulation results for both configurations are also presented.

DeWeese, Keith↗

Combinations of Earth Orientation Measurements: SPACE2004, COMB2004, and POLE2004

Independent Earth orientation measurements taken by the space-geodetic techniques of lunar and satellite laser ranging, very long baseline interferometry, and the global positioning system have been combined using a Kalman filter. The resulting combined Earth orientation series, SPACE2004, consists of values and uncertainties for Universal Time, polar motion, and their rates that span from September 28, 1976, to January 22, 2005, at daily intervals and is available in versions whose epochs are given at either midnight or noon. The space-geodetic measurements used to generate SPACE2004 have then been combined with optical astrometric measurements to form two additional combined Earth orientation series: (1) COMB2004, consisting of values and uncertainties for Universal Time, polar motion, and their rates that span from January 20, 1962, to January 22, 2005, at daily intervals and which is also available in versions whose epochs are given at either midnight or noon, and (2) POLE2004, consisting of values and uncertainties for polar motion and its rate that span from January 20, 1900, to January 20, 2005, at 30.4375-day intervals.

Gross, Richard R.↗

Iterative Magnetometer Calibration

This paper presents an iterative method for three-axis magnetometer (TAM) calibration that makes use of three existing utilities recently incorporated into the attitude ground support system used at NASA's Goddard Space Flight Center. The method combines attitude-independent and attitude-dependent calibration algorithms with a new spinning spacecraft Kalman filter to solve for biases, scale factors, nonorthogonal corrections to the alignment, and the orthogonal sensor alignment. The method is particularly well-suited to spin-stabilized spacecraft, but may also be useful for three-axis stabilized missions given sufficient data to provide observability.

Sedlak, Joseph↗

Addition of random run FM noise to the KPW time scale algorithm

The KPW (Kalman plus weights) time scale algorithm uses a Kalman filter to provide frequency and drift information to a basic time scale equation. This paper extends the algorithm to three-state clocks nd gives results for a simulated eight-clock ensemble.

time scale clock ensemble clock noise↗

Combinations of Earth Orientation Observations: SPACE94, COMB94, and POLE94

A Kalman filter has been used to combine all publicly available, independently determined measurements of the Earth's orientation taken by the modern, space-geodetic techniques of very long baseline interferometry, satellite laser ranging, lunar laser ranging, and the global positioning system. Prior to combining the data, tidal terms were removed from the UT1 measurements, outlying data points were deleted, series-specific corrections were applied for bias and rate, and the stated uncertainties of the measurements were adjusted by multiplying them by series-specific scale factors. Values for these bias- rate corrections and uncertainty scale factors were determined by an iterative, round-robin procedure wherein each data set is compared, in turn, to a combination of all other data sets. When applied to the measurements, the bias-rate corrections thus determined make the data sets agree with each other in bias and rate, and the uncertainty scale factors thus determined make the residual of each series (when differenced with a combination of all others) have a reduced chi-square of one. The corrected and adjusted series are then placed within an IERS reference frame by aligning them with the IERS Earth orientation series EOP (IERS)90C04. The result of combining these corrected, adjusted and aligned series is designated SPCE94 and spans October 6.0, 1976 to January 27.0, 1995 at daily intervals.

SPACE94 COMB94↗

Stochastic Modeling of Length of Day

We propose a general stochastic model for the UT1-LOD system and derive the corresponding Kalman filter model.

stochastic model autoregressive moving average tra↗

The Short-Term Prediction of Universal Time and Length-Of-Day Using Atmospheric Angular Momemtum

Methods to combine the diverse set of geodetic measurements of Earth orientation that are currently available, and to interpolate and extrapolate these data to generate an optimal estimate of Earth orientation, have been under development at JPL for a number of years. The strategy currently in use is a Kalman filtering scheme based on Earth orentation parameters and their excitation functions that incorporates stochastic models of the important physical processes and takes into account the variable form, quality, and temporal density of the data provided by diffent measurement services.

Universal Time↗

Interplanetary navigation using a continental baseline large antenna arrays

Navigation is a key component of interplanetary missions and must continue to be precise with the changing landscape of antenna design. Improvements for the Deep Space Network (DSN) may include the use of antenna arrays to simulate the power of a larger single antenna at much lower operating and construction costs. Therefore, it is necessary to test the performance of arrayed antennas from a navigational point-of-view. This initial investigation focuses on the performance of arrayed antennas from a navigational point-of-view. This initial investigation focuses on the performance of delta one-way range measurements using a shorter baseline with more data collection then current systems use. With all other parameter equal, the longer the baseline, the better the accuracy for navigation making the number of data packets very important. This trade study compares baseline distances ranging from 1 to 1000km with an in use baseline, looking at a due east baseline, a due north baseline at 45 degrees East of North. The precision of the baseline systems can be found through a simulated created for this purpose using the Jet Propulsion Lab based Monte navigation and mission design tool. The simulation combines the delta one-way range measurements with two-range and two-way Doppler measurements and puts the measurements through a Kalman filter to determine an orbit solution. Noise is added along with initial errors to give the simulation realism. This study is an important step towards the assessment of the utility of arrays for navigational purposes. The preliminary results have showed a decrease in reliability as the baseline is shortened but the larger continental baselines show comparable results t that of the current Goldstone to Canberra.

very long baseline interferometry (VLBI)↗

Automated daily processing of more than 1000 ground-based GPS receivers for studying intense ionospheric storms

To take advantage of the vast amount of GPS data, researchers use a number of techniques to estimate satellite and receiver interfrequency biases and the total electron content (TEC) of the ionosphere. Most techniques estimate vertical ionospheric structure and, simultaneously, hardware-related biases treated as nuisance parameters. These methods often are limited to 200 GPS receivers and use a sequential least squares or Kalman filter approach. The biases are later removed from the measurements to obtain unbiased TEC. In our approach to calibrating GPS receiver and transmitter interfrequency biases we take advantage of all available GPS receivers using a new processing algorithm based on the Global Ionospheric Mapping (GIM) software developed at the Jet Propulsion Laboratory. This new capability is designed to estimate receiver biases for all stations. We solve for the instrumental biases by modeling the ionospheric delay and removing it from the observation equation using precomputed GIM maps. The precomputed GIM maps rely on 200 globally distributed GPS receivers to establish the ''background'' used to model the ionosphere at the remaining 800 GPS sites.

total electron content (TEC)↗

Relative Attitude Determination of Earth Orbiting Formations Using GPS Receivers

Satellite formation missions require the precise determination of both the position and attitude of multiple vehicles to achieve the desired objectives. In order to support the mission requirements for these applications, it is necessary to develop techniques for representing and controlling the attitude of formations of vehicles. A generalized method for representing the attitude of a formation of vehicles has been developed. The representation may be applied to both absolute and relative formation attitude control problems. The technique is able to accommodate formations of arbitrarily large number of vehicles. To demonstrate the formation attitude problem, the method is applied to the attitude determination of a simple leader-follower along-track orbit formation. A multiplicative extended Kalman filter is employed to estimate vehicle attitude. In a simulation study using GPS receivers as the attitude sensors, the relative attitude between vehicles in the formation is determined 3 times more accurately than the absolute attitude.

Lightsey, E. Glenn↗

Convective Weather Forecast Quality Metrics for Air Traffic Management Decision-Making

Since numerical weather prediction models are unable to accurately forecast the severity and the location of the storm cells several hours into the future when compared with observation data, there has been a growing interest in probabilistic description of convective weather. The classical approach for generating uncertainty bounds consists of integrating the state equations and covariance propagation equations forward in time. This step is readily recognized as the process update step of the Kalman Filter algorithm. The second well known method, known as the Monte Carlo method, consists of generating output samples by driving the forecast algorithm with input samples selected from distributions. The statistical properties of the distributions of the output samples are then used for defining the uncertainty bounds of the output variables. This method is computationally expensive for a complex model compared to the covariance propagation method. The main advantage of the Monte Carlo method is that a complex non-linear model can be easily handled. Recently, a few different methods for probabilistic forecasting have appeared in the literature. A method for computing probability of convection in a region using forecast data is described in Ref. 5. Probability at a grid location is computed as the fraction of grid points, within a box of specified dimensions around the grid location, with forecast convection precipitation exceeding a specified threshold. The main limitation of this method is that the results are dependent on the chosen dimensions of the box. The examples presented Ref. 5 show that this process is equivalent to low-pass filtering of the forecast data with a finite support spatial filter. References 6 and 7 describe the technique for computing percentage coverage within a 92 x 92 square-kilometer box and assigning the value to the center 4 x 4 square-kilometer box. This technique is same as that described in Ref. 5. Characterizing the forecast, following the process described in Refs. 5 through 7, in terms of percentage coverage or confidence level is notionally sound compared to characterizing in terms of probabilities because the probability of the forecast being correct can only be determined using actual observations. References 5 through 7 only use the forecast data and not the observations. The method for computing the probability of detection, false alarm ratio and several forecast quality metrics (Skill Scores) using both the forecast and observation data are given in Ref. 2. This paper extends the statistical verification method in Ref. 2 to determine co-occurrence probabilities. The method consists of computing the probability that a severe weather cell (grid location) is detected in the observation data in the neighborhood of the severe weather cell in the forecast data. Probabilities of occurrence at the grid location and in its neighborhood with higher severity, and with lower severity in the observation data compared to that in the forecast data are examined. The method proposed in Refs. 5 through 7 is used for computing the probability that a certain number of cells in the neighborhood of severe weather cells in the forecast data are seen as severe weather cells in the observation data. Finally, the probability of existence of gaps in the observation data in the neighborhood of severe weather cells in forecast data is computed. Gaps are defined as openings between severe weather cells through which an aircraft can safely fly to its intended destination. The rest of the paper is organized as follows. Section II summarizes the statistical verification method described in Ref. 2. The extension of this method for computing the co-occurrence probabilities in discussed in Section HI. Numerical examples using NCWF forecast data and NCWD observation data are presented in Section III to elucidate the characteristics of the co-occurrence probabilities. This section also discusses the procedure for computing throbabilities that the severity of convection in the observation data will be higher or lower in the neighborhood of grid locations compared to that indicated at the grid locations in the forecast data. The probability of coverage of neighborhood grid cells is also described via examples in this section. Section IV discusses the gap detection algorithm and presents a numerical example to illustrate the method. The locations of the detected gaps in the observation data are used along with the locations of convective weather cells in the forecast data to determine the probability of existence of gaps in the neighborhood of these cells. Finally, the paper is concluded in Section V.

Chatterji, Gano B.↗

Attitude Issues on the Huygens Probe: Balloon Dropped Mock up Role in Determining Reconstruction Strategies During Descent in Lower Atmosphere

As part of the collaboration with Italian Space Agency on HASI instrument for Huygens mission, University of Padova has been conducting since 2001 scientific activity on Stratospheric Balloon Launches from the Trapani base in Sicily. The most recent boomerang flight in July 2003 has successfully flown a mock up of the Huygens probe hosting spares of flight scientific units and extra housekeeping and scientific sensors on a parachuted descent from 33 kilometre altitude. This work presents the studies conducted on attitude reconstruction of the probe, as well as the utilisation of iterative extended Kalman filtering in investigating vanes induced spin rate and in providing a baseline for the performance evaluation of Huygens accelerometers operations. Finally some possible contributions on the reconstruction of the lower part of Titan descent for Huygens probe are suggested based on the confrontation of sensor data for 2003 flight.

Bettanini, C.↗

Real-time Sub-cm Differential Orbit Determination of two Low-Earth Orbiters with GPS Bias Fixing

An effective technique for real-time differential orbit determination with GPS bias fixing is formulated. With this technique, only real-time GPS orbits and clocks are needed (available from the NASA Global Differential GPS System with 10-20 cm accuracy). The onboard, realtime orbital states of user satellites (few meters in accuracy) are used for orbit initialization and integration. An extended Kalman filter is constructed for the estimation of the differential orbit between the two satellites as well as a reference orbit, together with their associating dynamics parameters. Due to close proximity of the two satellites and of similar body shapes, the differential dynamics are highly common and can be tightly constrained which, in turn, strengthens the orbit estimation. Without explicit differencing of GPS data, double-differenced phase biases are formed by a transformation matrix. Integer-valued fixing of these biases are then performed which greatly strengthens the orbit estimation. A 9-day demonstration between GRACE orbits with baselines of approx.200 km indicates that approx.80% of the double-differenced phase biases can successfully be fixed and the differential orbit can be determined to approx.7 mm as compared to the results of onboard K-band ranging.

global positioning system (GPS)↗

Toward a Real-Time Measurement-Based System for Estimation of Helicopter Engine Degradation Due to Compressor Erosion

This paper presents a preliminary demonstration of an automated health assessment tool, capable of real-time on-board operation using existing engine control hardware. The tool allows operators to discern how rapidly individual turboshaft engines are degrading. As the compressor erodes, performance is lost, and with it the ability to generate power. Thus, such a tool would provide an instant assessment of the engine s fitness to perform a mission, and would help to pinpoint any abnormal wear or performance anomalies before they became serious, thereby decreasing uncertainty and enabling improved maintenance scheduling. The research described in the paper utilized test stand data from a T700-GE-401 turboshaft engine that underwent sand-ingestion testing to scale a model-based compressor efficiency degradation estimation algorithm. This algorithm was then applied to real-time Health Usage and Monitoring System (HUMS) data from a T700-GE-701C to track compressor efficiency on-line. The approach uses an optimal estimator called a Kalman filter. The filter is designed to estimate the compressor efficiency using only data from the engine s sensors as input.

Litt, Jonathan S.↗

Reduced-Order Modeling and Wavelet Analysis of Turbofan Engine Structural Response Due to Foreign Object Damage "FOD" Events

The development of a wavelet-based feature extraction technique specifically targeting FOD-event induced vibration signal changes in gas turbine engines is described. The technique performs wavelet analysis of accelerometer signals from specified locations on the engine and is shown to be robust in the presence of significant process and sensor noise. It is envisioned that the technique will be combined with Kalman filter thermal/ health parameter estimation for FOD-event detection via information fusion from these (and perhaps other) sources. Due to the lack of high-frequency FOD-event test data in the open literature, a reduced-order turbofan structural model (ROM) was synthesized from a finite-element model modal analysis to support the investigation. In addition to providing test data for algorithm development, the ROM is used to determine the optimal sensor location for FOD-event detection. In the presence of significant noise, precise location of the FOD event in time was obtained using the developed wavelet-based feature.

Turso, James A.↗

An Ensemble-Based Smoother with Retrospectively Updated Weights for Highly Nonlinear Systems

Monte Carlo computational methods have been introduced into data assimilation for nonlinear systems in order to alleviate the computational burden of updating and propagating the full probability distribution. By propagating an ensemble of representative states, algorithms like the ensemble Kalman filter (EnKF) and the resampled particle filter (RPF) rely on the existing modeling infrastructure to approximate the distribution based on the evolution of this ensemble. This work presents an ensemble-based smoother that is applicable to the Monte Carlo filtering schemes like EnKF and RPF. At the minor cost of retrospectively updating a set of weights for ensemble members, this smoother has demonstrated superior capabilities in state tracking for two highly nonlinear problems: the double-well potential and trivariate Lorenz systems. The algorithm does not require retrospective adaptation of the ensemble members themselves, and it is thus suited to a streaming operational mode. The accuracy of the proposed backward-update scheme in estimating non-Gaussian distributions is evaluated by comparison to the more accurate estimates provided by a Markov chain Monte Carlo algorithm.

Monte Carlo↗

Hierarchically Structured Non-Intrusive Sign Language Recognition

This work presents a hierarchically structured approach at the nonintrusive recognition of sign language from a monocular frontal view. Robustness is achieved through sophisticated localization and tracking methods, including a combined EM/CAMSHIFT overlap resolution procedure and the parallel pursuit of multiple hypotheses about hands position and movement. This allows handling of ambiguities and automatically corrects tracking errors. A biomechanical skeleton model and dynamic motion prediction using Kalman filters represents high level knowledge. Classification is performed by Hidden Markov Models. 152 signs from German sign language were recognized with an accuracy of 97.6%.

Zieren, Jorg↗

SPARTAN: A High-Fidelity Simulation for Automated Rendezvous and Docking Applications

bd Systems (a subsidiary of SAIC) has developed the Simulation Package for Autonomous Rendezvous Test and ANalysis (SPARTAN), a high-fidelity on-orbit simulation featuring multiple six-degree-of-freedom (6DOF) vehicles. SPARTAN has been developed in a modular fashion in Matlab/Simulink to test next-generation automated rendezvous and docking guidance, navigation,and control algorithms for NASA's new Vision for Space Exploration. SPARTAN includes autonomous state-based mission manager algorithms responsible for sequencing the vehicle through various flight phases based on on-board sensor inputs and closed-loop guidance algorithms, including Lambert transfers, Clohessy-Wiltshire maneuvers, and glideslope approaches The guidance commands are implemented using an integrated translation and attitude control system to provide 6DOF control of each vehicle in the simulation. SPARTAN also includes high-fidelity representations of a variety of absolute and relative navigation sensors that maybe used for NASA missions, including radio frequency, lidar, and video-based rendezvous sensors. Proprietary navigation sensor fusion algorithms have been developed that allow the integration of these sensor measurements through an extended Kalman filter framework to create a single optimal estimate of the relative state of the vehicles. SPARTAN provides capability for Monte Carlo dispersion analysis, allowing for rigorous evaluation of the performance of the complete proposed AR&D system, including software, sensors, and mechanisms. SPARTAN also supports hardware-in-the-loop testing through conversion of the algorithms to C code using Real-Time Workshop in order to be hosted in a mission computer engineering development unit running an embedded real-time operating system. SPARTAN also contains both runtime TCP/IP socket interface and post-processing compatibility with bdStudio, a visualization tool developed by bd Systems, allowing for intuitive evaluation of simulation results. A description of the SPARTAN architecture and capabilities is provided, along with details on the models and algorithms utilized and results from representative missions.

Turbe, Michael A.↗