Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimal estimation”

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 199 records · Page 11

Relationships of earthquakes (and earthquake-associated mass movements) and polar motion as determined by Kalman filtered, Very-Long-Baseline-Interferometry

A Kalman filter was designed to yield optimal estimates of geophysical parameters from Very Long Baseline Interferometry (VLBI) group delay data. The geophysical parameters are the polar motion components, adjustments to nutation in obliquity and longitude, and a change in the length of day parameter. The VLBI clock (and clock rate) parameters and atmospheric zenith delay parameters are estimated simultaneously. Filter background is explained. The IRIS (International Radio Interferometric Surveying) VLBI data are Kalman filtered. The resulting polar motion estimates are examined. There are polar motion signatures at the times of three large earthquakes occurring in 1984 to 1986: Mexico, 19 September, 1985 (Magnitude M sub s = 8.1); Chile, 3 March, 1985 (M sub s = 7.8); and Taiwan, 14 November, 1986 (M sub s = 7.8). Breaks in polar motion occurring about 20 days after the earthquakes appear to correlate well with the onset of increased regional seismic activity and a return to more normal seismicity (respectively). While the contribution of these three earthquakes to polar motion excitations is small, the cumulative excitation due to earthquakes, or seismic phenomena over a Chandler wobble damping period may be significant. Mechanisms for polar motion excitation due to solid earth phenomena are examined. Excitation functions are computed, but the data spans are too short to draw conclusions based on these data.

Preisig, Joseph Richard Mark↗

Stable Kalman filters for processing clock measurement data

Kalman filters have been used for some time to process clock measurement data. Due to instabilities in the standard Kalman filter algorithms, the results have been unreliable and difficult to obtain. During the past several years, stable forms of the Kalman filter have been developed, implemented, and used in many diverse applications. These algorithms, while algebraically equivalent to the standard Kalman filter, exhibit excellent numerical properties. Two of these stable algorithms, the Upper triangular-Diagonal (UD) filter and the Square Root Information Filter (SRIF), have been implemented to replace the standard Kalman filter used to process data from the Deep Space Network (DSN) hydrogen maser clocks. The data are time offsets between the clocks in the DSN, the timescale at the National Institute of Standards and Technology (NIST), and two geographically intermediate clocks. The measurements are made by using the GPS navigation satellites in mutual view between clocks. The filter programs allow the user to easily modify the clock models, the GPS satellite dependent biases, and the random noise levels in order to compare different modeling assumptions. The results of this study show the usefulness of such software for processing clock data. The UD filter is indeed a stable, efficient, and flexible method for obtaining optimal estimates of clock offsets, offset rates, and drift rates. A brief overview of the UD filter is also given.

Clements, P. A.↗

The shape of Eros

Monte Carlo simulations are presently used to optimize estimation, ascertain associated errors, and guide bias-correction procedures, for the Eros polar silhouette convex hull that has been estimated from radar echo spectra. This hull is trapezoidal; this nonaxisymmetric shape may account for odd harmonics in Eros' echo spectral signature as a function of rotation phase. Additional constraints have been obtained for the figure of Eros through the inversion of the optical lightcurve to estimate the asteroid's two-dimensional average of the three-dimensional shape. This 'mean cross-section' and the polar silhouette exhibit similar elongations.

Ostro, S. J.↗

Relative-Error-Covariance Algorithms

Two algorithms compute error covariance of difference between optimal estimates, based on data acquired during overlapping or disjoint intervals, of state of discrete linear system. Provides quantitative measure of mutual consistency or inconsistency of estimates of states. Relative-error-covariance concept applied, to determine degree of correlation between trajectories calculated from two overlapping sets of measurements and construct real-time test of consistency of state estimates based upon recently acquired data.

Bierman, Gerald J.↗

An experimental study of nonlinear dynamic system identification

A technique for robust identification of nonlinear dynamic systems is developed and illustrated using both simulations and analog experiments. The technique is based on the Minimum Model Error optimal estimation approach. A detailed literature review is included in which fundamental differences between the current approach and previous work is described. The most significant feature of the current work is the ability to identify nonlinear dynamic systems without prior assumptions regarding the form of the nonlinearities, in constrast to existing nonlinear identification approaches which usually require detailed assumptions of the nonlinearities. The example illustrations indicate that the method is robust with respect to prior ignorance of the model, and with respect to measurement noise, measurement frequency, and measurement record length.

Stry, Greselda I.↗

An experimental study of nonlinear dynamic system identification

A technique based on the Minimum Model Error optimal estimation approach is employed for robust identification of a nonlinear dynamic system. A simple harmonic oscillator with quadratic position feedback was simulated on an analog computer. With the aid of analog measurements and an assumed linear model, the Minimum Model Error Algorithm accurately identifies the quadratic nonlinearity. The tests demonstrate that the method is robust with respect to prior ignorance of the nonlinear system model and with respect to measurement record length, regardless of initial conditions.

Stry, Greselda I.↗

A drag measurement technique for free piston shock tunnels

A new technique is described for measuring drag with 100-microsecond rise time on a nonlifting model in a free piston shock tunnel. The technique involves interpretation of the stress waves propagating within the model and its support. A finite element representation and spectral methods are used to obtain a mean square optimal estimate of the time history of the aerodynamic loading. Thus, drag is measured instantaneously and the previous restriction caused by the mechanical time constant of balances is overcome. The effectiveness of the balance is demonstrated by measuring the drag on cones with 5 and 15 deg semi-vertex angles in nominally Mach 5.6 flow with stagnation enthalpies from 2.6 to 33 MJ/kg.

Sanderson, S. R.↗

Robust Kalman filter design for predictive wind shear detection

Severe, low-altitude wind shear is a threat to aviation safety. Airborne sensors under development measure the radial component of wind along a line directly in front of an aircraft. In this paper, optimal estimation theory is used to define a detection algorithm to warn of hazardous wind shear from these sensors. To achieve robustness, a wind shear detection algorithm must distinguish threatening wind shear from less hazardous gustiness, despite variations in wind shear structure. This paper presents statistical analysis methods to refine wind shear detection algorithm robustness. Computational methods predict the ability to warn of severe wind shear and avoid false warning. Comparative capability of the detection algorithm as a function of its design parameters is determined, identifying designs that provide robust detection of severe wind shear.

Stratton, Alexander D.↗

Ionospheric gravity wave measurements with the USU dynasonde

A method for the measurement of ionospheric Gravity Wave (GW) using the USU Dynasonde is outlined. This method consists of a series of individual procedures, which includes functions for data acquisition, adaptive scaling, polarization discrimination, interpolation and extrapolation, digital filtering, windowing, spectrum analysis, GW detection, and graphics display. Concepts of system theory are applied to treat the ionosphere as a system. An adaptive ionogram scaling method was developed for automatically extracting ionogram echo traces from noisy raw sounding data. The method uses the well known Least Mean Square (LMS) algorithm to form a stochastic optimal estimate of the echo trace which is then used to control a moving window. The window tracks the echo trace, simultaneously eliminating the noise and interference. Experimental results show that the proposed method functions as designed. Case studies which extract GW from ionosonde measurements were carried out using the techniques described. Geophysically significant events were detected and the resultant processed results are illustrated graphically. This method was also developed for real time implementation in mind.

Berkey, Frank T.↗

Altimeter and gravity data analysis

The studies carried out under this grant fell into two broad areas. The first area was the analysis of surface gravity data with the ultimate goal of providing normal equations that could be used in combination with normal equations from the analysis of satellite orbit perturbations to obtain an optimal estimate of the gravitational potential coefficients of the Earth. The second main research activity was the estimation of gravity anomalies in ocean areas from satellite altimeter data. Such anomalies could enable the improved calibration of potential coefficient models derived solely from the analysis of orbital perturbation information. The studies in these two areas are discussed.

Rapp, Richard H.↗

Global average concentration and trend for hydroxyl radicals deduced from ALE/GAGE trichloroethane (methyl chloroform) data for 1978-1990

An optimal estimation inversion scheme is utilized with atmospheric data and emission estimates to determined the globally averaged CH3CCl3 tropospheric lifetime and OH concentration. The data are taken from atmospheric measurements from surface stations of 1,1,1-trichloroethane and show an annual increase of 4.4 +/- 0.2 percent. Industrial emission estimates and a small oceanic loss rate are included, and the OH concentration for the same period (1978-1990) are incorporated at 1.0 +/- 0.8 percent/yr. The positive OH trend is consistent with theories regarding OH and ozone trends with respect to land use and global warming. Attention is given to the effects of the ENSO on the CH3CCl3 data and the assumption of continuing current industrial anthropogenic emissions. A novel tropical atmospheric tracer-transport mechanism is noted with respect to the CH3CCl3 data.

Prinn, R.↗

Attitude determination using vector observations: A fast optimal matrix algorithm

The attitude matrix minimizing Wahba's loss function is computed directly by a method that is competitive with the fastest known algorithm for finding this optimal estimate. The method also provides an estimate of the attitude error covariance matrix. Analysis of the special case of two vector observations identifies those cases for which the TRIAD or algebraic method minimizes Wahba's loss function.

Markley, F. Landis↗

Attitude determination using vector observations - A fast optimal matrix algorithm

The attitude matrix minimizing Wahba's loss function is computed directly by a method that is competitive with the fastest known algorithm for finding this optimal estimate. The method also provides an estimate of the attitude error covariance matrix. Analysis of the special case of two vector observations identifies those cases for which the TRIAD or algebraic method minimizes Wahba's loss function.

Markley, F. L.↗

Advances in time-scale algorithms

The term clock is usually used to refer to a device that counts a nearly periodic signal. A group of clocks, called an ensemble, is often used for time keeping in mission critical applications that cannot tolerate loss of time due to the failure of a single clock. The time generated by the ensemble of clocks is called a time scale. The question arises how to combine the times of the individual clocks to form the time scale. One might naively be tempted to suggest the expedient of averaging the times of the individual clocks, but a simple thought experiment demonstrates the inadequacy of this approach. Suppose a time scale is composed of two noiseless clocks having equal and opposite frequencies. The mean time scale has zero frequency. However if either clock fails, the time-scale frequency immediately changes to the frequency of the remaining clock. This performance is generally unacceptable and simple mean time scales are not used. First, previous time-scale developments are reviewed and then some new methods that result in enhanced performance are presented. The historical perspective is based upon several time scales: the AT1 and TA time scales of the National Institute of Standards and Technology (NIST), the A.1(MEAN) time scale of the US Naval observatory (USNO), the TAI time scale of the Bureau International des Poids et Measures (BIPM), and the KAS-1 time scale of the Naval Research laboratory (NRL). The new method was incorporated in the KAS-2 time scale recently developed by Timing Solutions Corporation. The goal is to present time-scale concepts in a nonmathematical form with as few equations as possible. Many other papers and texts discuss the details of the optimal estimation techniques that may be used to implement these concepts.

Stein, S. R.↗

Ground-based measurements of water vapor in the middle atmosphere

We present measurements of the middle atmospheric water vapor mixing ratio profile obtained using the ground-based Naval Research Laboratory (NRL) water vapor millimeter-wave spectrometer (WVMS) instrument at the Jet Propulsion Laboratory (JPL) Table Mountain Observatory. The measurements cover a period of 262 days from January 23, 1992, to October 13, 1992. During this campaign it was possible to retrieve useful daily mixing ratio profiles for 186 days. We thus have a nearly continuous record of water vapor mixing ratios for altitudes from approximately equals 35 to 75 km. The retrievals are obtained using the optimal estimation method. Details of the error analysis are presented, and a technique is introduced that reduces baseline effects and helps to estimate the baseline error. The high-altitude (greater than or approx. equal to 65 km) data show a sharp rise prior to the expected maximum near the summer solstice and a gradual decline in the following months. The mixing ratios generally peak between 55 and 65 km, at which point the mixing ratios are 6 - 7 parts per million by volume. The highest peaks occur in January, May, and October.

Nedoluha, Gerald E.↗

Evaluation of phase-diversity techniques for solar-image restoration

Phase-diversity techniques provide a novel observational method for overcomming the effects of turbulence and instrument-induced aberrations in ground-based astronomy. Two implementations of phase-diversity techniques that differ with regard to noise model, estimator, optimization algorithm, method of regularization, and treatment of edge effects are described. Reconstructions of solar granulation derived by applying these two implementations to common data sets are shown to yield nearly identical images. For both implementations, reconstructions from phase-diverse speckle data (involving multiple realizations of turbulence) are shown to be superior to those derived from conventional phase-diversity data (involving a single realization). Phase-diverse speckle reconstructions are shown to achieve near diffraction-limited resolution and are validated by internal and external consistency tests, including a comparison with a reconstruction using a well-accepted speckle-imaging method.

Paxman, Richard G.↗

Ozone vertical distribution retrieval from ground-based high resolution infrared solar spectra

A practical procedure for the retrieval of ozone vertical profiles from ground-based high resolution Fourier transform infrared solar spectra has been developed. The analysis is based on a multilayer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method of Rodgers. The 1002.6-1003.2 cm(exp -1) spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. This interval contains numerous ozone lines covering a range of intensities and providing retrieval sensitivity from ground level to about 35 km. Characterization of the method and an error analysis have been performed. For a spectral resolution of 0.05-0.01 cm(exp -1) and a signal-to-noise ratio greater than or equal to 100 the retrieval is stable with a vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Synthetic spectra studies show that the a priori profile and weak constraints selected for the retrievals introduce no significant biases for a wide range of ozone profiles.

Pougatchev, N. S.↗

Spectroscopy of Solid State Laser Materials

We retrieved the vertical distribution of ozone from a series 0.005-0.013/cm resolution infrared solar spectra recorded with the McMath Fourier Transform spectrometer at the Kitt Peak National Solar Observatory. The analysis is based on a multi-layer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method by Rodgers. The 1002.6-1003.2/cm spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. The characterization and error analysis of the method have been performed. It was shown that for the Kitt Peak spectral resolution and typical signal-to-noise ratio (greater than or equal to 100) the retrieval is stable, with the vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Spectra recorded from 1980 through 1993 have been analyzed. The retrieved total ozone and vertical profiles have been compared with total ozone mapping spectrometer (TOMS) satellite total columns for the location and dates of the Kitt Peak Measurements and about 100 ozone ozonesoundings and Brewer total column measurements from Palestine, Texas, from 1979 to 1985. The total ozone measurements agree to +/- 2%. The retrieved profiles reproduce the seasonally averaged variations with altitude, including the ozone spring maximum and fall minimum measured by Palestine sondes, but up to 15% differences in the absolute values are obtained.

Buoncristiani, A. M.↗