Search NASA⌕ Search

SEARCH · Search NASA

Results for “likelihood 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 451 records · Page 25

Likelihood Analysis for Mega Pixel Maps

The derivation of cosmological parameters from astrophysical data sets routinely involves operations counts which scale as O(N(exp 3) where N is the number of data points. Currently planned missions, including MAP and Planck, will generate sky maps with N(sub d) = 10(exp 6) or more pixels. Simple "brute force" analysis, applied to such mega-pixel data, would require years of computing even on the fastest computers. We describe an algorithm which allows estimation of the likelihood function in the direct pixel basis. The algorithm uses a conjugate gradient approach to evaluate X2 and a geometric approximation to evaluate the determinant. Monte Carlo simulations provide a correction to the determinant, yielding an unbiased estimate of the likelihood surface in an arbitrary region surrounding the likelihood peak. The algorithm requires O(N(sub d)(exp 3/2) operations and O(Nd) storage for each likelihood evaluation, and allows for significant parallel computation.

Kogut, Alan J.↗

Noncoherent DTTLs for Symbol Synchronization

Noncoherent data-transition tracking loops (DTTLs) have been proposed for use as symbol synchronizers in digital communication receivers. [Communication- receiver subsystems that can perform their assigned functions in the absence of synchronization with the phases of their carrier signals ( carrier synchronization ) are denoted by the term noncoherent, while receiver subsystems that cannot function without carrier synchronization are said to be coherent. ] The proposal applies, more specifically, to receivers of binary phase-shift-keying (BPSK) signals generated by directly phase-modulating binary non-return-to-zero (NRZ) data streams onto carrier signals having known frequencies but unknown phases. The proposed noncoherent DTTLs would be modified versions of traditional DTTLs, which are coherent. The symbol-synchronization problem is essentially the problem of recovering symbol timing from a received signal. In the traditional, coherent approach to symbol synchronization, it is necessary to establish carrier synchronization in order to recover symbol timing. A traditional DTTL effects an iterative process in which it first generates an estimate of the carrier phase in the absence of symbol-synchronization information, then uses the carrier-phase estimate to obtain an estimate of the symbol-synchronization information, then feeds the symbol-synchronization estimate back to the carrier-phase-estimation subprocess. In a noncoherent symbol-synchronization process, there is no need for carrier synchronization and, hence, no need for iteration between carrier-synchronization and symbol- synchronization subprocesses. The proposed noncoherent symbolsynchronization process is justified theoretically by a mathematical derivation that starts from a maximum a posteriori (MAP) method of estimation of symbol timing utilized in traditional, coherent DTTLs. In that MAP method, one chooses the value of a variable of interest (in this case, the offset in the estimated symbol timing) that causes a likelihood function of symbol estimates over some number of symbol periods to assume a maximum value. In terms that are necessarily oversimplified to fit within the space available for this article, it can be said that the mathematical derivation involves a modified interpretation of the likelihood function that lends itself to noncoherent DTTLs. The proposal encompasses both linear and nonlinear noncoherent DTTLs. The performances of both have been computationally simulated; for comparison, the performances of linear and nonlinear coherent DTTLs have also been computationally simulated. The results of these simulations show that, among other things, the expected mean-square timing errors of coherent and noncoherent DTTLs are relatively insensitive to window width. The results also show that at high signal-to-noise ratios (SNRs), the performances of the noncoherent DTTLs approach those of their coherent counterparts at, while at low SNRs, the noncoherent DTTLs incur penalties of the order of 1.5 to 2 dB.

Simon, Marvin↗

Maximum likelihood algorithm using an efficient scheme for computing sensitivities and parameter confidence intervals

Improved techniques for estimating airplane stability and control derivatives and their standard errors are presented. A maximum likelihood estimation algorithm is developed which relies on an optimization scheme referred to as a modified Newton-Raphson scheme with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort compared to integrating the analytically-determined sensitivity equations or using a finite difference scheme. An aircraft estimation problem is solved using real flight data to compare MNRES with the commonly used modified Newton-Raphson technique; MNRES is found to be faster and more generally applicable. Parameter standard errors are determined using a random search technique. The confidence intervals obtained are compared with Cramer-Rao lower bounds at the same confidence level. It is observed that the nonlinearity of the cost function is an important factor in the relationship between Cramer-Rao bounds and the error bounds determined by the search technique.

Murphy, P. C.↗

Regularized Differentiation for Bioburden Density Estimation in Planetary Protection

In this paper, we propose and investigate the performance of two novel shrinkage estimators for bioburden density estimation in planetary protection. The estimators are based on the regularized differentiation of a cumulative count of colony forming units collected throughout the data collecting session or the life cycle of the entire mission. The regularized differentiation recasts the problem of bioburden density estimation as a linear least squares problem. The least squares problem is then solved through regularization techniques, such as truncated singular value decomposition and penalized least squares. The regularization is necessary to avoid noise amplification during the differentiation of noisy data. The two regularization estimators are compared with four other commonly used estimators to simultaneously evaluate the means of multivariable independent Poisson distributions: the maximum likelihood, noninformative Bayes estimator with Jeffreys prior, Empirical Bayes using conjugate gamma-Poisson model with gamma parameters selected by method of moments, and the Clevenson-Zidek estimator. It is shown through computer-simulated data that the regularized differentiation based on ridge regression has the smallest mean-squared error among all estimators. The analysis of shrinkage mechanism implemented by regularized differentiation is performed, and it is shown that the regularized differentiation amounts to performing a weighted averaging of all the samples. The weights are determined by the regularization parameter automatically selected by the L-curve technique. Since the method of least squares makes no distributional assumptions about the data, it presents an attractive technique for bioburden density estimation when there are concerns about the misspecification of the distributional model. The paper concludes with the analysis of the bioburden data collected during InSight mission and directions for future work.

97 - MATHEMATICS AND COMPUTING↗

Proportion estimation using prior cluster purities

The prior distribution of CLASSY component purities is studied, and this information incorporated into maximum likelihood crop proportion estimators. The method is tested on Transition Year spring small grain segments.

Terrell, G. R.↗

The atmospheric lifetime experiment. III - Lifetime methodology and application to three years of CFCL3 data

Observations of the chlorofluorocarbon CFCl3 obtained several times daily over the period July 1978 to June 1981 at Adrigole, Ireland; Ragged Point, Barbados; Point Matatula, American Samoa; and Cape Grim, Tasmania are reported. In addition, observations at Cape Meares, Oregon are given for the period January 1980 to June 1981. On January 1, 1980, the average mixing ratio of CFCl3 in the lower troposphere is esimated to have been 168 pptv, and this is calculated to have been increasing 5.7 percent annually. Assuming that the only destruction of CFCl3 occurs in the stratosphere, the lifetime, on January 1, 1980, estimated by a trend technique is 83 + 73, or -27 years; the lifetime estimated from the global inventory of CFCl3 is to + 89 or -25 years. The maximum likelihood current lifetime estimate obtained by combining the estimates from both analysis techniques is 78 years.

Cunnold, D. M.↗

Identification of dynamic systems, theory and formulation

The problem of estimating parameters of dynamic systems is addressed in order to present the theoretical basis of system identification and parameter estimation in a manner that is complete and rigorous, yet understandable with minimal prerequisites. Maximum likelihood and related estimators are highlighted. The approach used requires familiarity with calculus, linear algebra, and probability, but does not require knowledge of stochastic processes or functional analysis. The treatment emphasizes unification of the various areas in estimation in dynamic systems is treated as a direct outgrowth of the static system theory. Topics covered include basic concepts and definitions; numerical optimization methods; probability; statistical estimators; estimation in static systems; stochastic processes; state estimation in dynamic systems; output error, filter error, and equation error methods of parameter estimation in dynamic systems, and the accuracy of the estimates.

Maine, R. E.↗

Probabilities and statistics for backscatter estimates obtained by a scatterometer

Methods for the recovery of winds near the surface of the ocean from measurements of the normalized radar backscattering cross section must recognize and make use of the statistics (i.e., the sampling variability) of the backscatter measurements. Radar backscatter values from a scatterometer are random variables with expected values given by a model. A model relates backscatter to properties of the waves on the ocean, which are in turn generated by the winds in the atmospheric marine boundary layer. The effective wind speed and direction at a known height for a neutrally stratified atmosphere are the values to be recovered from the model. The probability density function for the backscatter values is a normal probability distribution with the notable feature that the variance is a known function of the expected value. The sources of signal variability, the effects of this variability on the wind speed estimation, and criteria for the acceptance or rejection of models are discussed. A modified maximum likelihood method for estimating wind vectors is described. Ways to make corrections for the kinds of errors found for the Seasat SASS model function are described, and applications to a new scatterometer are given.

Pierson, Willard J., Jr.↗

Path Following with Slip Compensation for a Mars Rover

A software system for autonomous operation of a Mars rover is composed of several key algorithms that enable the rover to accurately follow a designated path, compensate for slippage of its wheels on terrain, and reach intended goals. The techniques implemented by the algorithms are visual odometry, full vehicle kinematics, a Kalman filter, and path following with slip compensation. The visual-odometry algorithm tracks distinctive scene features in stereo imagery to estimate rover motion between successively acquired stereo image pairs, by use of a maximum-likelihood motion-estimation algorithm. The full-vehicle kinematics algorithm estimates motion, with a no-slip assumption, from measured wheel rates, steering angles, and angles of rockers and bogies in the rover suspension system. The Kalman filter merges data from an inertial measurement unit (IMU) and the visual-odometry algorithm. The merged estimate is then compared to the kinematic estimate to determine whether and how much slippage has occurred. The kinematic estimate is used to complement the Kalman-filter estimate if no statistically significant slippage has occurred. If slippage has occurred, then a slip vector is calculated by subtracting the current Kalman filter estimate from the kinematic estimate. This slip vector is then used, in conjunction with the inverse kinematics, to determine the wheel velocities and steering angles needed to compensate for slip and follow the desired path.

Helmick, Daniel↗

On the estimation of climatological Z-R relationships

A statistical framework for climatological Z-R parameter estimation is developed and simulation experiments are conducted to examine sampling properties of the estimators. Both parametric and nonparametric models are considered. For parametric models, it is shown that Z-R parameters can be estimated by maximum likelihood, a procedure with optimal large sample properties. A general nonparametric framework for climatological Z-R estimation is also developed. Nonparametric procedures are attractive because of their flexibility in dealing with certain types of measurement errors common to radar data. Simulation experiments show that even under favorable assumptions on error characteristics of radar and raingages, large datasets are required to obtain accurate Z-R parameter estimates. Another important conclusion is that estimation results are generally quite sensitive to radar and raingage measurement thresholds. For fixed sample size, the simulation results can be used to provide quantitative assessments of the accuracy of Z-R model parameter estimates. These results are particularly useful for error analysis of precipitation products that are derived using climatological Z-R relations. One example is the large-area rainfall estimates derived using the height-area rainfall threshold (HART) technique.

Krajewski, Witold F.↗

Aerodynamic derivatives for an oblique wing aircraft estimated from flight data by using a maximum likelihood technique

There are several practical problems in using current techniques with five degree of freedom equations to estimate the stability and control derivatives of oblique wing aircraft from flight data. A technique was developed to estimate these derivatives by separating the analysis of the longitudinal and lateral directional motion without neglecting cross coupling effects. Although previously applied to symmetrical aircraft, the technique was not expected to be adequate for oblique wing vehicles. The application of the technique to flight data from a remotely piloted oblique wing aircraft is described. The aircraft instrumentation and data processing were reviewed, with particular emphasis on the digital filtering of the data. A complete set of flight determined stability and control derivative estimates is presented and compared with predictions. The results demonstrated that the relatively simple approach developed was adequate to obtain high quality estimates of the aerodynamic derivatives of such aircraft.

Maine, R. E.↗

Toward The Development of Hailstorm Climatologies Derived From Reanalyses and Infared/Passive Microwave Satellite Imagers

Geostationary satellite imagers, such as those of the Geostationary Operational Environmental Satellite (GOES) and Meteosat series, provide both historical and near-real-time observations of cloud top patterns that are commonly associated with severe convection. Environmental conditions favorable for severe weather are thought to be represented well by reanalyses. Predicting exactly where convection and costly storm hazards like hail will occur using models or satellite imagery alone, however, is extremely challenging. The multivariate combination of satellite-observed cloud patterns with reanalysis environmental parameters, linked to United States Next Generation Weather Radar- (NEXRAD-) estimated Maximum Expected Size of Hail (MESH) using a deep neural network (DNN), enables estimation of potentially severe hail likelihood for any observed storm cell. These estimates are specifically designed to make hail likelihood distinctions based on satellite-indicated points of deep convection within environments favorable for storm development. We seek an approach that can be used to estimate climatological hailstorm frequency and risk throughout the historical satellite data record. This presentation demonstrates that statistical distributions of convective parameters from satellite and reanalysis show separation between non-severe/severe hailstorm classes for predictors including overshooting cloud top temperature and area characteristics, convective available potential energy, vertical wind shear, 500 hPa temperature, mid-level lapse rate, precipitable water, and convective inhibition. These complex, multivariate predictor relationships are exploited within a DNN to produce a hail likelihood metric with a critical success index of 0.504 and Heidke skill score of 0.403, which is exceptional among recent analogous hail studies. Furthermore, applications of the DNN to select case studies demonstrate good qualitative agreement between hail likelihood and MESH. These hail classifications are aggregated across an 11-year GOES-12/13 image database to derive a hail frequency and severity climatology, which denotes the Central Plains, the Midwest, and northwestern Mexico as being the most hail-prone regions within the domain studied. Opportunities for training and applying DNN-based hailstorm predictions to recently developed GOES-8/10/12/13/16 and Meteosat Second Generation convective storm detection and characterization climatologies over South America and South Africa, respectively, will also be presented.

Kristopher Bedka↗

On-Orbit Multi-Field Wavefront Control with a Kalman Filter

A document describes a multi-field wavefront control (WFC) procedure for the James Webb Space Telescope (JWST) on-orbit optical telescope element (OTE) fine-phasing using wavefront measurements at the NIRCam pupil. The control is applied to JWST primary mirror (PM) segments and secondary mirror (SM) simultaneously with a carefully selected ordering. Through computer simulations, the multi-field WFC procedure shows that it can reduce the initial system wavefront error (WFE), as caused by random initial system misalignments within the JWST fine-phasing error budget, from a few dozen micrometers to below 50 nm across the entire NIRCam Field of View, and the WFC procedure is also computationally stable as the Monte-Carlo simulations indicate. With the incorporation of a Kalman Filter (KF) as an optical state estimator into the WFC process, the robustness of the JWST OTE alignment process can be further improved. In the presence of some large optical misalignments, the Kalman state estimator can provide a reasonable estimate of the optical state, especially for those degrees of freedom that have a significant impact on the system WFE. The state estimate allows for a few corrections to the optical state to push the system towards its nominal state, and the result is that a large part of the WFE can be eliminated in this step. When the multi-field WFC procedure is applied after Kalman state estimate and correction, the stability of fine-phasing control is much more certain. Kalman Filter has been successfully applied to diverse applications as a robust and optimal state estimator. In the context of space-based optical system alignment based on wavefront measurements, a KF state estimator can combine all available wavefront measurements, past and present, as well as measurement and actuation error statistics to generate a Maximum-Likelihood optimal state estimator. The strength and flexibility of the KF algorithm make it attractive for use in real-time optical system alignment when WFC alone cannot effectively align the system.

Lou, John↗

Aerodynamic coefficient estimation for dynamic wind tunnel models

Maximum likelihood parameter identification is used to estimate the aerodynamic coefficients of a 'flying' aircraft model which is mounted in a transonic wind tunnel by a system of cables and pulleys. The model's motion is governed by both cable and aerodynamic forces, where the parameters to be identified are functions of the aerodynamic forces but cable forces greatly predominate. The test data were successfully processed assuming a linear equation model in the identification procedure. In one instance, the procedure was able to detect a failing measurement instrument. The final parameter estimates are compared to estimates obtained in independent tests of a similar aircraft model.

Mohr, R. L.↗

Guest Investigator Studies with the Compton Gamma Ray Observatory

The cumulative all-sky survey by the Compton Gamma Ray Observatory (CGRO), composed of data acquired during the first three years of the mission, included a number of regions of very limited exposure. The most glaring deficiency in coverage was toward the region of the South Galactic Pole (SGP), which received significantly less exposure than other directions- by a factor of at least 2 to 3. Furthermore, nearly all of the SGP exposure was accumulated during the first year of the mission. Since blazars are known to be time-variable, and of unknown duty cycle, a pointing of the CCRO in that direction was considered highly desirable, and long overdue. In addition, data gathered from a pointing toward the SGP and its comparison with comprehensive data available for the North Galactic Pole would be extremely valuable to investigators studying the extragalactic diffuse emission. The reasons outlined above prompted our initiation of a Cycle 4 campaign to systematically search with EGRET and COMPTEL for gamma-ray emission from sources near the South Galactic Pole. The Cycle 4 SGP campaign consisted of tnvo 14-day observations separated in in time by approximately 10 months. The temporal separation of the observations was requested to allow a test for possible variations in the detected sources. Our primary targets were 38 FSRQs which lie within 30 degrees of the SGP, and which satisfy the basic criteria for candidate gamma-ray AGNs,flat-spectrum radio sources, many of which exhibit blazar-type properties). These targets were selected from the standard references, and from the available on-line databases (e.g., the NASA Extragalactic Database, NED), as the most promising AGN targets in the vicinity of the SGP. A 30 radius from the SGP was chosen as the boundary of our survey, since the selected targets would then fall within the most sensitive portion of the fields of view of EGRET and COMPTEL (i.e., within a 30 zenith angle), for a CGRO pointing directed exactly at the SGP. Our South Galactic Pole Survey yielded a number of exciting results. The EGRET data were analyzed using the maximum likelihood techniques to estimate the intensity, spectrum, and position of gamma-ray sources in the field of view. Our analysis revealed four sources at energies greater than 100 MeV with likelihood ratios corresponding to greater than 30 detections (Vestrand et al. 1996). One of the sources is associated with the well known gamma-ray blazar PKS 0208-512, but the other three were previously unknown. Among the new detections was PKS 2155-304 which is often considered a prototype of the x-ray selected BL Lacs. PKS 2155-304, which was also detected at hard x-ray energies by CGRO/OSSE, is one of the brightest BL Lac objects in the sky at optical through x-ray energies and has a history of rapid, strong multiwavelength variability. As such, it has been the subject of intensive, contemporaneous, multiwavelength monitoring covering radio frequencies to x-ray energies.

Vestrand, W. T.↗

Real-time combining of residual carrier array signals using ML weight estimates

A real-time digital signal combining system for use with array feeds is proposed. The combining system attempts to compensate for signal-to-noise ratio (SNR) loss resulting from antenna deformations induced by gravitational and atmospheric effects. The combining weights are obtained directly from the observed residual carrier samples in each channel using a 'sliding-window' implementation of a maximum-likelihood (ML) parameter estimator. It is shown that with averaging times of about 0.1 s, combining loss for a seven-element array can be limited to about 0.1 dB in a realistic operational environment. This result suggests that the real-time combining system proposed here is capable of recovering virtually all of the signal power captured by the array feed, even in the presence of severe wind gusts and similar disturbances.

Vilnrotter, Victor A.↗

Maximum likelihood failure detection techniques applied to the shuttle orbiter reaction control subsystem

A technique for on-board detection and identification of hard failures and leaks of the shuttle orbiter reaction control subsystem jets, during the orbital flight phase, is presented. The method uses gimbal angle and linear accelerometer measurements from the orbiter inertial measurement unit and requires no additional hardware. Extended Kalman filters with residual traps are employed for state estimation, and generalized likelihood ratio tests for jet failure identification. Rigid body simulation results indicate identification times of less than 2 seconds for hard jet failures and less than 70 seconds for jet leaks.

Deckert, J. C.↗