Search NASA⌕ Search

SEARCH · Search NASA

Results for “randomized algorithms”

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 361 records · Page 20

Synergistic Retrievals of Ice in High Clouds From Elastic Backscatter Lidar, Ku-band Radar and Submillimeter Wave Radiometer Observations

In this study, we investigate the synergy of elastic backscatter lidar, Ku-band radar, and sub-millimeter-wave radiometer measurements in the retrieval of ice from satellite observations. The synergy is analyzed through the generation of a large dataset of IceWater Content (IWC) profiles and simulated lidar, radar and radiometer observations. The characteristics of the instruments e.g. frequencies, sensitivities, etc. are set based on the expected characteristics of instruments of the Atmosphere Observing System (AOS) mission. A hold-out validation methodology is used to assess the accuracy of the IWC profiles retrieved from various combinations of observations from the three instruments. Specifically, the IWC and associated observations are randomly divided into two datasets, one for training and the other for evaluation. The training dataset is used to train the retrieval algorithm, while the evaluation dataset is used to assess the retrieval performance. The dataset of IWC profiles is derived from CloudSat reflectivity and CALIOP lidar observations. The retrieval of the ice water content IWC profiles from the computed observations is achieved in two steps. In the first step, a class, out of 18 potential classes characterized by different vertical distribution of IWC, is estimated from the observations. The 18 classes are predetermined based on the k-Means clustering algorithm. In the second step, the IWC profile is estimated using an Ensemble Kalman Smoother (EKS) algorithm that uses the estimated class as a priori information. The results of the study show that the synergy of lidar, radar, and radiometer observations is significant in the retrieval of the IWC profiles. Nevertheless, it should be mentioned that this synergy was found under idealized conditions, and additional work might be required to materialize it in practice. The inclusion of the lidar backscatter observations in the retrieval process has a larger impact on the retrieval performance than the inclusion of the radar observations. As ice clouds have a significant impact on atmospheric radiative processes, this work is relevant to ongoing efforts to reduce uncertainties in climate analyses and projections.

Mircea Grecu↗

Fast Quantum Algorithm for Predicting Descriptive Statistics of Stochastic Processes

Stochastic processes are used as a modeling tool in several sub-fields of physics, biology, and finance. Analytic understanding of the long term behavior of such processes is only tractable for very simple types of stochastic processes such as Markovian processes. However, in real world applications more complex stochastic processes often arise. In physics, the complicating factor might be nonlinearities; in biology it might be memory effects; and in finance is might be the non-random intentional behavior of participants in a market. In the absence of analytic insight, one is forced to understand these more complex stochastic processes via numerical simulation techniques. In this paper we present a quantum algorithm for performing such simulations. In particular, we show how a quantum algorithm can predict arbitrary descriptive statistics (moments) of N-step stochastic processes in just O(square root of N) time. That is, the quantum complexity is the square root of the classical complexity for performing such simulations. This is a significant speedup in comparison to the current state of the art.

Williams Colin P.↗

Lila: Optimal Dispatching in Probabilistic Temporal Networks using Monte Carlo Tree Search

Executing a Probabilistic Simple Temporal Network (PSTN) amounts at scheduling, i.e. \textit{dispatch}, a set of events under time uncertainty. This constitutes a NP-hard online optimization problem. The right execution time must be dynamically assigned to each event of the PSTN such that the temporal constraints are met, whereas activity durations are progressively observed as the execution unfolds. We propose a dispatching algorithm based on Monte Carlo Tree Search, called Lila, with the following characteristics: (i) it is an anytime algorithm, both offline and online, proven asymptotically optimal; (ii) it returns the current probability of success, either before or at any moment during operations; (iii) it handles any possible continuous or discrete, even non-parametric, probability distributions, as well as inter-dependencies between random variables, exogenous and endogenous uncertainty; and (iv) can be easily extended to handle probabilistic external events, PSTNs with resources, PSTNs with cutoff times and precondition chains, etc. Lila is universal in the sense that it can handle any dispatching protocol, simply by specifying it to the algorithm. It has the unlimited flexibility offered by the simulation paradigm, whilst it asymptotically converges to optimal decisions and/or robustness approximations.

Chien, Steve A.↗

Comparison of genetic algorithms with conjugate gradient methods

Genetic algorithms for mathematical function optimization are modeled on search strategies employed in natural adaptation. Comparisons of genetic algorithms with conjugate gradient methods, which were made on an IBM 1800 digital computer, show that genetic algorithms display superior performance over gradient methods for functions which are poorly behaved mathematically, for multimodal functions, and for functions obscured by additive random noise. Genetic methods offer performance comparable to gradient methods for many of the standard functions.

Bosworth, J. L.↗

A study of digital holographic filters generation. Phase 2: Digital data communication system, volume 1

An empirical study of the performance of the Viterbi decoders in bursty channels was carried out and an improved algebraic decoder for nonsystematic codes was developed. The hybrid algorithm was simulated for the (2,1), k = 7 code on a computer using 20 channels having various error statistics, ranging from pure random error to pure bursty channels. The hybrid system outperformed both the algebraic and the Viterbi decoders in every case, except the 1% random error channel where the Viterbi decoder had one bit less decoding error.

Ingels, F. M.↗

Turbulence control on airborne laser platform

An avctive flow control device to generate large-scale, periodic structures in a turbulent shear flow is developed. Together with adaptive optics, the device may be used on airborne laser platforms to reduce or eliminate optical distortion caused by the turbulence in the aircraft's boundary layer. A flat plate towed in a water channel is used as a test bed. A cyclic jet issuing from a spanwise slot is used to collect the turbulent boundary layer for a finite time during its 'on' period. When the jet is turned 'off', all of the turbulent fluid is released instantaneously in one large eddy that convects downstream. Flow visualization and hot-film probe measurements are used together with pattern recognition algorithms to demonstrate the viability of the flow control method. The instantaneous velocity signal is used to compute important statistical quantities of the random velocity field, such as the mean, the root-mean-square, the spectral distribution, and the probability density function. When optimized for a given boundary layer, the cyclic jet produces periodic structures that are qualitatively similar to the random, naturally occurring ones. These structures seem to trigger the onset of bursting events near the wall. Thus, the present device generates periodic structures in both the outer and inner regions of a turbulent boundary layer.

Gad-El-hak, Mohamed↗

Estimation of optical flow in airborne electro-optical sensors by stochastic approximation

The essence of motion or range estimation by passive electrooptical means is the ability to determine the correspondence of picture elements in pairs of image frames and to estimate their coordinates and their disparity (relative shifts) in the image plane of an electrooptical imaging sensor. The disparity can be in successive frames due to self-motion or in simultaneous frames of a stereo pair. A key issue is to provide these estimates on-line. This paper describes the theoretical background of such an interframe shift estimator. It is based on a stochastic gradient algorithm, specifically implementing a form of stochastic approximation, which can achieve rapid convergence of the shift estimate. Analytical and numerical simulation examples for random texture and isolated features validate the feasibility and the effectiveness of the estimator.

Merhav, S. J.↗

Aerosol Abundances and Optical Characteristics in the Pacific Basin Free Troposphere

During NASA's Global Backscatter Experiment (GLOBE) mission flights in November 1989 and May 1990, a DC-8 research aircraft probed the Pacific Basin free troposphere for about 90 flight hours in each month between +72 and -62 degrees latitude, +130 and -120 degrees longitude, and up to 39,000 feet pressure altitudes. Aerosols were sampled continuously in situ by optical particle counters to measure concentration and particle size, and during 48 10-min intervals during each mission by wire impactors for concentration, size, composition, phase and shape analyses. The optical particle counters cover a particle diameter range between 0.3 and 20 microns; wire impactors extend the range down to 0.03 microns. Results of particle number, size, shape, together with the assumption of a refractive index corresponding to (NH4)2SO4 to account for the prevalence of aerosol sulfur, were utilized in a Mie algorithm to calculate aerosol extinction and backscatter for a range of wavelengths (0.385 less than lambda less than 10.64 microns). Computations for 22 randomly selected size distributions yield coefficients of extinction E(0.525) = (2.03 +/- 1.20) x 10(exp -4) km(exp -1) and backscatter beta(0.525) = (6.45 +/- 3.49) x 10(exp -6) km(exp -1) sr(exp -1) in the visible, and E(10.64) = (8.13 +/- 6.47) x 10(exp -6) km(exp -1) and beta(10.64) = (9.98 +/- 10.69) x 10(exp -8) km(exp -1) sr(exp -1) in the infrared, respectively. Large particles (D greater than 0.3 microns) contribute two-thirds to the total extinction in the visible (lambda = 0.525 microns), and almost 100% in the infrared (lambda = 10.64 microns). These results have been used to define an IR optical aerosol climatology of the Pacific Basin free troposphere, from which it follows that the infrared backscatter coefficient at lambda = 9.25 microns wavelength fluctuates between 5.0 x 10(exp -10) and 2.0 x 10(exp -7) km(exp -1) sr(exp -1) with a modal value 2.0 x 10(exp -8) km(exp -1) sr(exp -1).

Pueschel, R. F.↗

Using MERRA Gridded Innovations for Quantifying Uncertainties in Analysis Fields and Diagnosing Observing System Inhomogeneities

MERRA is a NASA reanalysis for the satellite era using a major new version of the Goddard Earth Observing System Data Assimilation System Version 5 (GEOS-5). The project focuses on historical analyses of the hydrological cycle on a broad range of weather and climate time scales and places the NASA EOS suite of observations in a climate context. The characterization of uncertainty in reanalysis fields is a commonly requested feature by users of such data. While intercomparison with reference data sets is common practice for ascertaining the realism of the datasets, such studies typically are restricted to long term climatological statistics and seldom provide state dependent measures of the uncertainties involved. In principle, variational data assimilation algorithms have the ability of producing error estimates for the analysis variables (typically surface pressure, winds, temperature, moisture and ozone) consistent with the assumed background and observation error statistics. However, these "perceived error estimates" are expensive to obtain and are limited by the somewhat simplistic errors assumed in the algorithm. The observation minus forecast residuals (innovations) by-product of any assimilation system constitutes a powerful tool for estimating the systematic and random errors in the analysis fields. Unfortunately, such data is usually not readily available with reanalysis products, often requiring the tedious decoding of large datasets and not so-user friendly file formats. With MERRA we have introduced a gridded version of the observations/innovations used in the assimilation process, using the same grid and data formats as the regular datasets. Such dataset empowers the user with the ability of conveniently performing observing system related analysis and error estimates. The scope of this dataset will be briefly described. We will present a systematic analysis of MERRA innovation time series for the conventional observing system, including maximum-likelihood estimates of background and observation errors, as well as global bias estimates. Starting with the joint PDF of innovations and analysis increments at observation locations we propose a technique for diagnosing bias among the observing systems, and document how these contextual biases have evolved during the satellite era covered by MERRA.

da Silva, Arlindo↗

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

Mapping Surface Vapor Pressure Deficits From Geostationary Satellites for Fire Weather Monitoring

The increase in the wildfires were observed globally in accordance with global warming, and to real- time monitoring of wildfire risk in broad scale is demanded for wildfire management to prevent the spread of wildfires. Scientists invented a lot of indices to assess the wildfire risk. Vapor Pressure Deficit (VPD) is one of the most important meteorological components for those indices. Compared to other components of fire weather indices, VPD can change quickly from lower risk to higher risk even in sub-hourly. Therefore, real-time fire risk monitoring requires high-resolution and high- temporal VPD spatial map. Here, we developed VPD estimation method using the GOES Advanced Baseline Imager (ABI) data. Unlike the polar-orbital satellite data, the ABI can observe target region every 10 minutes, so that we can estimate VPD for fire weather in real-time. The method used to estimate VPD is same with the algorithm of NASA Earth Exchange Gridded Daily Meteorology (NEX- GDM), which estimate meteorological variables from ground weather observation and spatial variables based on random forest (RF). We calculated RF importance to select bands of ABI as input of the model. To validate our results, we compared the spatial pattern of our VPD data with the Real- Time Mesoscale Analysis (RTMA) data over the conterminous USA. We sought possibility of applying our method to the region where no real-time high-resolution weather data is available, such as South America. The developed method can produce real-time high-resolution high-frequent VPD data in the continental scale. The derived data from GOES ABI could contribute to improve the fire weather monitoring and lead to prevent wildfires.

Hirofumi Hashimoto↗

Contextual classification of multispectral image data

A general method is presented for exploiting both spatial and spectral information when classifying multispectral image data. This statistical classification algorithm utilizes the tendency of certain ground cover classes to be more likely to occur in some contexts than others. The theoretical model assumes the two-dimensional array of random observations and a 0-1 loss function, a distribution of the p-context array that is spatially invariant, and class-conditional independence for the observations. The problems that prevent the immediate use of this context classifier are the need for a generally applicable method for making adequate estimates of the context distribution and a reduction in the computational intensivity of the classifier. The former problem is being approached by a method that raises the relative frequency value for each class configuration to a power and uses the result as the context distribution estimate. The second is being approached by searching for a less computationally intensive algorithm.

Tilton, J. C.↗

Designing laminated composites using random search techniques

A computer program called UWCODA is presented. UWCODA is intended to assist in the design, analysis and optimization of composite plates. UWCODA combines a state-of-the-art global optimization algorithm (Improving Hit and Run) with classical lamination theory. Optimization results are presented for simple loading conditions as well as for complex, biaxial load conditions. The computer code proved to be very effective in the design of composite plates.

Graesser, D. L.↗

Galileo spacecraft modal identification using an eigensystem realization algorithm

A modal parameter identification technique referred to as the Eigensystem Realization Algorithm (ERA) was applied to free-response measurements from the Galileo spacecraft modal survey test. The data were recorded following single-point random excitation of the structure. This work is one phase in a research project coordinated by the Jet Propulsion Laboratory to compare the performance of various contemporary identification techniques using Galileo data. Principal emphasis is placed on estimating the accuracy of the ERA-identified modal parameters. Various accuracy indicators, such as Modal Amplitude Coherence and Modal Phase Collinearity, are discussed. More than 20 modes of the spacecraft were identified, demonstrating the ability of the ERA method to determine the dynamics of such complex structures using only a few seconds of test data.

Pappa, R. S.↗

Using Markov Models of Fault Growth Physics and Environmental Stresses to Optimize Control Actions

A generalized Markov chain representation of fault dynamics is presented for the case that available modeling of fault growth physics and future environmental stresses can be represented by two independent stochastic process models. A contrived but representatively challenging example will be presented and analyzed, in which uncertainty in the modeling of fault growth physics is represented by a uniformly distributed dice throwing process, and a discrete random walk is used to represent uncertain modeling of future exogenous loading demands to be placed on the system. A finite horizon dynamic programming algorithm is used to solve for an optimal control policy over a finite time window for the case that stochastic models representing physics of failure and future environmental stresses are known, and the states of both stochastic processes are observable by implemented control routines. The fundamental limitations of optimization performed in the presence of uncertain modeling information are examined by comparing the outcomes obtained from simulations of an optimizing control policy with the outcomes that would be achievable if all modeling uncertainties were removed from the system.

Bole, Brian↗

Design and Calibration of Autonomous Coherent Doppler Lidar for Space Missions

Developed a new algorithm for the simulation of three dimensional homogeneous turbulent velocity fields. For typical atmospheric conditions it is impossible to produce a simulated velocity field that simultaneously satisfy a given spatial correlation and the corresponding spatial spectrum because of spectral aliasing. The new algorithms produce a turbulent velocity field which has accurate spatial correlations which is required for performance predictions from space-based systems. Developed a new algorithm for extracting the spatial statistics of the atmospheric velocity field using coherent Doppler lidar. The performance of the algorithm was compared with past methods and the new algorithm produces useful results for space-based data, which was not possible before. Developed new methods for verification of the errors in ground-based and space-based Doppler lidar wind measurements. These new methods do not require independent in situ data. This is an important issue for the verification of space-based Doppler lidar measurements of the global wind field. The performance of the new algorithm was compared with past results for both space-based and ground-based operation. The new algorithm has the best performance and is the only algorithm that performed satisfactory for spacebased operation. The performance of coherent Doppler lidar for a space missions with various scanning geometries was determined using computer simulation which contained the effects of random instrumental velocity errors, wind shear, wind variability along the range-gate and from shot-to-shot, and random variations in atmospheric aerosol backscatter over the measurement volume. The bias in the velocity estimates was small and the accuracy in the is typically less than 0.5 m/s for high signal conditions. For a large number of shot per velocity estimate, the threshold signal level for acceptable estimates is proportional to the number of shots to the minus one half power. This agrees with previous results determined for ground-based measurements. The use of multi-element optical detectors for autonomous operation of coherent Doppler lidar was shown to be a very promising technique. Optimal detector geometries were determined by computer simulation of performance: for ground-based testing with a fixed calibration target and for space-based operation using the random surface returns. The effects of refractive turbulence on ground-based calibration of coherent Doppler lidar was determined by computer simulations and compared with theoretical predictions. New techniques were required to correctly predict performance for the focused beam geometry commonly used for verification of space-based operation. An improved velocity estimator was evaluated for space-based applications were signal shot measurements are used to produce vector wind measurements. This permits more accurate measurements when the signal level is not known a priori or not available from multiple shot measurements. The average Doppler lidar signal spectrum including the effects of velocity turbulence was derived and calculated. This permits new estimation algorithms for turbulence based on spectral estimates. In situ atmospheric measurements were conducted and analyzed using an instrumented kite-platform. This work helps provide the required in situ data for verification of Doppler lidar velocity statistics.

Frehlich, Rod G.↗

Vision Algorithms to Determine Shape and Distance for Manipulation of Unmodeled Objects

This paper discusses the development of a robotic system for general use in an unstructured environment. This is illustrated through pick and place of randomly positioned, un-modeled objects. There are many applications for this project, including rock collection for the Mars Surveyor Program. This system is demonstrated with a Puma560 robot, Barrett hand, Cognex vision system, and Cimetrix simulation and control, all running on a PC. The demonstration consists of two processes: vision system and robotics. The vision system determines the size and location of the unknown objects. The robotics part consists of moving the robot to the object, configuring the hand based on the information from the vision system, then performing the pick/place operation. This work enhances and is a part of the Low Cost Virtual Collaborative Environment which provides remote simulation and control of equipment.

Montes, Leticia↗

Separation of man-made and natural patterns in high-altitude imagery of agricultural areas

A nonstationary linear digital filter is designed and implemented which extracts the natural features from high-altitude imagery of agricultural areas. Essentially, from an original image a new image is created which displays information related to soil properties, drainage patterns, crop disease, and other natural phenomena, and contains no information about crop type or row spacing. A model is developed to express the recorded brightness in a narrow-band image in terms of man-made and natural contributions and which describes statistically the spatial properties of each. The form of the minimum mean-square error linear filter for estimation of the natural component of the scene is derived and a suboptimal filter is implemented. Nonstationarity of the two-dimensional random processes contained in the model requires a unique technique for deriving the optimum filter. Finally, the filter depends on knowledge of field boundaries. An algorithm for boundary location is proposed, discussed, and implemented.

Samulon, A. S.↗