Search NASA⌕ Search

SEARCH · Search NASA

Results for “Likelihood Function”

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 289 records · Page 16

A Parameter Estimation Scheme for Multiscale Kalman Smoother (MKS) Algorithm Used in Precipitation Data Fusion

A new approach is presented in this paper to effectively obtain parameter estimations for the Multiscale Kalman Smoother (MKS) algorithm. This new approach has demonstrated promising potentials in deriving better data products based on data of different spatial scales and precisions. Our new approach employs a multi-objective (MO) parameter estimation scheme (called MO scheme hereafter), rather than using the conventional maximum likelihood scheme (called ML scheme) to estimate the MKS parameters. Unlike the ML scheme, the MO scheme is not simply built on strict statistical assumptions related to prediction errors and observation errors, rather, it directly associates the fused data of multiple scales with multiple objective functions in searching best parameter estimations for MKS through optimization. In the MO scheme, objective functions are defined to facilitate consistency among the fused data at multiscales and the input data at their original scales in terms of spatial patterns and magnitudes. The new approach is evaluated through a Monte Carlo experiment and a series of comparison analyses using synthetic precipitation data. Our results show that the MKS fused precipitation performs better using the MO scheme than that using the ML scheme. Particularly, improvements are significant compared to that using the ML scheme for the fused precipitation associated with fine spatial resolutions. This is mainly due to having more criteria and constraints involved in the MO scheme than those included in the ML scheme. The weakness of the original ML scheme that blindly puts more weights onto the data associated with finer resolutions is overcome in our new approach.

multiscale↗

Image Discrimination Models Predict Object Detection in Natural Backgrounds

Object detection involves looking for one of a large set of object sub-images in a large set of background images. Image discrimination models only predict the probability that an observer will detect a difference between two images. In a recent study based on only six different images, we found that discrimination models can predict the relative detectability of objects in those images, suggesting that these simpler models may be useful in some object detection applications. Here we replicate this result using a new, larger set of images. Fifteen images of a vehicle in an other-wise natural setting were altered to remove the vehicle and mixed with the original image in a proportion chosen to make the target neither perfectly recognizable nor unrecognizable. The target was also rotated about a vertical axis through its center and mixed with the background. Sixteen observers rated these 30 target images and the 15 background-only images for the presence of a vehicle. The likelihoods of the observer responses were computed from a Thurstone scaling model with the assumption that the detectabilities are proportional to the predictions of an image discrimination model. Three image discrimination models were used: a cortex transform model, a single channel model with a contrast sensitivity function filter, and the Root-Mean-Square (RMS) difference of the digital target and background-only images. As in the previous study, the cortex transform model performed best; the RMS difference predictor was second best; and last, but still a reasonable predictor, was the single channel model. Image discrimination models can predict the relative detectabilities of objects in natural backgrounds.

Ahumada, Albert J., Jr.↗

Spatial Statistics of Large Astronomical Databases: An Algorithmic Approach

In this AISRP, the we have demonstrated that the correlation function i) can be calculated for MAP in minutes (about 45 minutes for Planck) on a modest 500Mhz workstation ii) the corresponding method, although theoretically suboptimal, produces nearly optimal results for realistic noise and cut sky. This trillion fold improvement in speed over the standard maximum likelihood technique opens up tremendous new possibilities, which will be persued in the follow up.

Szapudi, Istvan↗

Defining Nitrogen Kinetics for Air Break in Prebreath

Actual tissue nitrogen (N2) kinetics are complex; the uptake and elimination is often approximated with a single half-time compartment in statistical descriptions of denitrogenation [prebreathe(PB)] protocols. Air breaks during PB complicate N2 kinetics. A comparison of symmetrical versus asymmetrical N2 kinetics was performed using the time to onset of hypobaric decompression sickness (DCS) as a surrogate for actual venous N2 tension. METHODS: Published results of 12 tests involving 179 hypobaric exposures in altitude chambers after PB, with and without airbreaks, provide the complex protocols from which to model N2 kinetics. DCS survival time for combined control and airbreaks were described with an accelerated log logistic model where N2 uptake and elimination before, during, and after the airbreak was computed with a simple exponential function or a function that changed half-time depending on ambient N2 partial pressure. P1N2-P2 = (Delta)P defined decompression dose for each altitude exposure, where P2 was the test altitude and P1N2 was computed N2 pressure at the beginning of the altitude exposure. RESULTS: The log likelihood (LL) without decompression dose (null model) was -155.6, and improved (best-fit) to -97.2 when dose was defined with a 240 min half-time for both N2 elimination and uptake during the PB. The description of DCS survival time was less precise with asymmetrical N2 kinetics, for example, LL was -98.9 with 240 min half-time elimination and 120 min half-time uptake. CONCLUSION: The statistical regression described survival time mechanistically linked to symmetrical N2 kinetics during PBs that also included airbreaks. The results are data-specific, and additional data may change the conclusion. The regression is useful to compute additional PB time to compensate for an airbreak in PB within the narrow range of tested conditions.

Conkin, Johnny↗

A Capillary-Based Static Phase Separator for Highly Variable Wetting Conditions

The invention, a static phase separator (SPS), uses airflow and capillary wetting characteristics to passively separate a two-phase (liquid and air) flow. The device accommodates highly variable liquid wetting characteristics. The resultant design allows for a range of wetting properties from about 0 to over 90 advancing contact angle, with frequent complete separation of liquid from gas observed when using appropriately scaled test conditions. Additionally, the design accommodates a range of air-to-liquid flow-rate ratios from only liquid flow to over 200:1 air-to-liquid flow rate. The SPS uses a helix input section with an ice-cream-cone-shaped constant area cross section (see figure). The wedge portion of the cross section is on the outer edge of the helix, and collects the liquid via centripetal acceleration. The helix then passes into an increasing cross-sectional area vane region. The liquid in the helix wedge is directed into the top of capillary wedges in the liquid containment section. The transition from diffuser to containment section includes a 90 change in capillary pumping direction, while maintaining inertial direction. This serves to impinge the liquid into the two off-center symmetrical vanes by the airflow. Rather than the airflow serving to shear liquid away from the capillary vanes, the design allows for further penetration of the liquid into the vanes by the air shear. This is also assisted by locating the air exit ports downstream of the liquid drain port. Additionally, any droplets not contained in the capillary vanes are re-entrained downstream by a third opposing capillary vane, which directs liquid back toward the liquid drain port. Finally, the dual air exit ports serve to slow the airflow down, and to reduce the likelihood of shear. The ports are stove-piped into the cavity to form an unfriendly capillary surface for a wetting fluid to carryover. The liquid drain port is located at the start of the containment region, allowing for draining the bulk fluid in a continuous circuit. The functional operation of the SPS involves introducing liquid flow (from a human body, a syringe, or other source) to the two-phase inlet while an air fan pulls on the air exit lines. The fan is operated until the liquid is fully introduced. The system is drained by negative pressure on the liquid drain lines when the SPS containment system is full.

Thomas, Evan A.↗

Super Resolution of Remote Sensing Images Using Edge-Directed Radial Basis Functions

Edge-Directed Radial Basis Functions (EDRBF) are used to compute super resolution(SR) image from a given set of low resolution (LR) images differing in subpixel shifts. The algorithm is tested on remote sensing images and compared for accuracy with other well-known algorithms such as Iterative Back Projection (IBP), Maximum Likelihood (ML) algorithm, interpolation of scattered points using Nearest Neighbor (NN) and Inversed Distance Weighted (IDW) interpolation, and Radial Basis Functin(RBF) . The accuracy of SR depends on various factors besides the algorithm (i) number of subpixel shifted LR images (ii) accuracy with which the LR shifts are estimated by registration algorithms (iii) and the targeted spatial resolution of SR. In our studies, the accuracy of EDRBF is compared with other algorithms keeping these factors constant. The algorithm has two steps: i) registration of low resolution images and (ii) estimating the pixels in High Resolution (HR) grid using EDRBF. Experiments are conducted by simulating LR images from a input HR image with different sub-pixel shifts. The reconstructed SR image is compared with input HR image to measure the accuracy of the algorithm using sum of squared errors (SSE). The algorithm has outperformed all of the algorithms mentioned above. The algorithm is robust and is not overly sensitive to the registration inaccuracies.

Mareboyana, Manohar M.↗

Evaluation of errors in prior mean and variance in the estimation of integrated circuit failure rates using Bayesian methods

The critical point of any Bayesian analysis concerns the choice and quantification of the prior information. The effects of prior data on a Bayesian analysis are studied. Comparisons of the maximum likelihood estimator, the Bayesian estimator, and the known failure rate are presented. The results of the many simulated trails are then analyzed to show the region of criticality for prior information being supplied to the Bayesian estimator. In particular, effects of prior mean and variance are determined as a function of the amount of test data available.

Fletcher, B. C.↗

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.↗

Corn and soybean Landsat MSS classification performance as a function of scene characteristics

In order to fully utilize remote sensing to inventory crop production, it is important to identify the factors that affect the accuracy of Landsat classifications. The objective of this study was to investigate the effect of scene characteristics involving crop, soil, and weather variables on the accuracy of Landsat classifications of corn and soybeans. Segments sampling the U.S. Corn Belt were classified using a Gaussian maximum likelihood classifier on multitemporally registered data from two key acquisition periods. Field size had a strong effect on classification accuracy with small fields tending to have low accuracies even when the effect of mixed pixels was eliminated. Other scene characteristics accounting for variability in classification accuracy included proportions of corn and soybeans, crop diversity index, proportion of all field crops, soil drainage, slope, soil order, long-term average soybean yield, maximum yield, relative position of the segment in the Corn Belt, weather, and crop development stage.

Batista, G. T.↗

Vulnerability of manned spacecraft to crew loss from orbital debris penetration

Orbital debris growth threatens the survival of spacecraft systems from impact-induced failures. Whereas the probability of debris impact and spacecraft penetration may currently be calculated, another parameter of great interest to safety engineers is the probability that debris penetration will cause actual spacecraft or crew loss. Quantifying the likelihood of crew loss following a penetration allows spacecraft designers to identify those design features and crew operational protocols that offer the highest improvement in crew safety for available resources. Within this study, a manned spacecraft crew survivability (MSCSurv) computer model is developed that quantifies the conditional probability of losing one or more crew members, P(sub loss/pen), following the remote likelihood of an orbital debris penetration into an eight module space station. Contributions to P(sub loss/pen) are quantified from three significant penetration-induced hazards: pressure wall rupture (explosive decompression), fragment-induced injury, and 'slow' depressurization. Sensitivity analyses are performed using alternate assumptions for hazard-generating functions, crew vulnerability thresholds, and selected spacecraft design and crew operations parameters. These results are then used to recommend modifications to the spacecraft design and expected crew operations that quantitatively increase crew safety from orbital debris impacts.

Williamsen, J. E.↗

Probing the Masses of the PSR JO621+1002 Binary System Through Relativistic Apsidal Motion

Orbital, spin and astrometric parameters of the millisecond pulsar PSR J0621+1002 have been determined through six years of timing observations at three radio telescopes. The chief result is a measurement of the rate of periastron advance, omega=0 deg.0116 +/-0 deg.0008/yr. Interpreted as a general relativistic effect, this implies the sum of the pulsar mass, m(1), and the companion mass, m(2), to be M=m(1)+m(2)= 2.81 +/-0.30 solar mass. The Keplerian parameters rule out certain combinations of m(1) and m(2), as does the non-detection of Shapiro delay in the pulse arrival times. These constraints, together with the assumption that the companion is a white dwarf, lead to the maximum likelihood values m(1)=1.69((sup +0.30)(sub -0.30)) solar mass and m(2)=0.98((sup +0.32)(sub -0.12) solar mass (68% confidence). The other major finding is that the pulsar experiences dramatic variability in its dispersion measure (DM), with gradients as steep as 0.013 pc/cu cm/yr. A structure function analysis of the DM variations uncovers spatial fluctuations in the interstellar electron density that cannot be fit to a single power law, unlike the Kolmogorov turbulent spectrum that has been seen in the direction of other pulsars. Other results from the timing analysis include the first measurements of the pulsar's proper motion, mu=3.5+/-0.3 mas/yr, and of its spin-down rate, dP/dt=4.7 x 10(exp -20), which, when corrected for kinematic biases and combined with the pulse period, P=28.8 ms, gives a characteristic age of 1.1 x 10(exp 10) yr and a surface magnetic field strength of 1.2 x 10 (exp 9) G.

Spaver, Eric M.↗

A statistical model of the human core-temperature circadian rhythm

We formulate a statistical model of the human core-temperature circadian rhythm in which the circadian signal is modeled as a van der Pol oscillator, the thermoregulatory response is represented as a first-order autoregressive process, and the evoked effect of activity is modeled with a function specific for each circadian protocol. The new model directly links differential equation-based simulation models and harmonic regression analysis methods and permits statistical analysis of both static and dynamical properties of the circadian pacemaker from experimental data. We estimate the model parameters by using numerically efficient maximum likelihood algorithms and analyze human core-temperature data from forced desynchrony, free-run, and constant-routine protocols. By representing explicitly the dynamical effects of ambient light input to the human circadian pacemaker, the new model can estimate with high precision the correct intrinsic period of this oscillator ( approximately 24 h) from both free-run and forced desynchrony studies. Although the van der Pol model approximates well the dynamical features of the circadian pacemaker, the optimal dynamical model of the human biological clock may have a harmonic structure different from that of the van der Pol oscillator.

NASA Discipline Regulatory Physiology↗

Modeling X-Ray Binary Evolution in Normal Galaxies: Insights from SINGS

We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) populations and theoretical models of their production. We construct observational X-ray luminosity functions (oXLFs) using Chandra observations of 12 late-type galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining XRB synthetic models, constructed with the population synthesis code StarTrack, with observational star formation histories (SFHs). We identify highest-likelihood models both for individual galaxies and globally, averaged over the full galaxy sample. Individual tXLFs successfully reproduce about half of oXLFs, but for some galaxies we are unable to find underlying source populations, indicating that galaxy SFHs and metallicities are not well matched and/or XRB modeling requires calibration on larger observational samples. Given these limitations, we find that best models are consistent with a product of common envelope ejection efficiency and central donor concentration approx.. = 0.1, and a 50% uniform - 50% "twins" initial mass-ratio distribution. We present and discuss constituent subpopulations of tXLFs according to donor, accretor and stellar population characteristics. The galaxy-wide X-ray luminosity due to low-mass and high-mass XRBs, estimated via our best global model tXLF, follows the general trend expected from the L(sub X) - star formation rate and L(sub X) - stellar mass relations of Lehmer et al. Our best models are also in agreement with modeling of the evolution both of XRBs over cosmic time and of the galaxy X-ray luminosity with redshift.

Normal Galaxies↗

Health management system for rocket engines

The functional framework of a failure detection algorithm for the Space Shuttle Main Engine (SSME) is developed. The basic algorithm is based only on existing SSME measurements. Supplemental measurements, expected to enhance failure detection effectiveness, are identified. To support the algorithm development, a figure of merit is defined to estimate the likelihood of SSME criticality 1 failure modes and the failure modes are ranked in order of likelihood of occurrence. Nine classes of failure detection strategies are evaluated and promising features are extracted as the basis for the failure detection algorithm. The failure detection algorithm provides early warning capabilities for a wide variety of SSME failure modes. Preliminary algorithm evaluation, using data from three SSME failures representing three different failure types, demonstrated indications of imminent catastrophic failure well in advance of redline cutoff in all three cases.

Nemeth, Edward↗

Likelihood estimation for distributed parameter models for NASA Mini-MAST truss

A maximum likelihood estimation for distributed parameter models of large flexible structures was formulated. Distributed parameter models involve far fewer unknown parameters than independent modal characteristics or finite element models. The closed form solutions for the partial differential equations with corresponding boundary conditions were derived. The closed-form expressions of sensitivity functions led to highly efficient algorithms for analyzing ground or on-orbit test results. For an illustration of this approach, experimental data of the NASA Mini-MAST truss was used. The estimations of modal properties involve lateral bending modes and torsional modes. The results show that distributed parameter models are promising in the parameter estimation of large flexible structures.

Shen, Ji-Yao↗

Likelihood estimation for distributed parameter models of large beam-like structures

A maximum likelihood estimation for distributed parameter models of large flexible structures was formulated. Distributed parameter models involve far fewer unknown parameters than independent modal characteristics or finite element models. The closed form solutions for the partial differential equations with corresponding boundary conditions were derived. The closed-form expressions of sensitivity functions led to highly efficient algorithms for analyzing ground or on-orbit test results. For an illustration of this approach, experimental data of the NASA Mini-Mast truss was used. The estimations of modal properties involve lateral bending modes and torsional modes. The results show that distributed parameter models are promising in the parameter estimation of large flexible structures.

Shen, J.-Y.↗

Quaternion Averaging

Many applications require an algorithm that averages quaternions in an optimal manner. For example, when combining the quaternion outputs of multiple star trackers having this output capability, it is desirable to properly average the quaternions without recomputing the attitude from the the raw star tracker data. Other applications requiring some sort of optimal quaternion averaging include particle filtering and multiple-model adaptive estimation, where weighted quaternions are used to determine the quaternion estimate. For spacecraft attitude estimation applications, derives an optimal averaging scheme to compute the average of a set of weighted attitude matrices using the singular value decomposition method. Focusing on a 4-dimensional quaternion Gaussian distribution on the unit hypersphere, provides an approach to computing the average quaternion by minimizing a quaternion cost function that is equivalent to the attitude matrix cost function Motivated by and extending its results, this Note derives an algorithm that deterniines an optimal average quaternion from a set of scalar- or matrix-weighted quaternions. Rirthermore, a sufficient condition for the uniqueness of the average quaternion, and the equivalence of the mininiization problem, stated herein, to maximum likelihood estimation, are shown.

Markley, F. Landis↗

A dual-mode generalized likelihood ratio approach to self-reorganizing digital flight control system design

Analytic techniques have been developed for detecting and identifying abrupt changes in dynamic systems. The GLR technique monitors the output of the Kalman filter and searches for the time that the failure occured, thus allowing it to be sensitive to new data and consequently increasing the chances for fast system recovery following detection of a failure. All failure detections are based on functional redundancy. Performance tests of the F-8 aircraft flight control system and computerized modelling of the technique are presented.

Source record↗