Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix inversion”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Polarimetric radar data decomposition and interpretation

Significant efforts have been made to decompose polarimetric radar data into several simple scattering components. The components which are selected because of their physical significance can be used to classify SAR (Synthetic Aperture Radar) image data. If particular components can be related to forest parameters, inversion procedures may be developed to estimate these parameters from the scattering components. Several methods have been used to decompose an averaged Stoke's matrix or covariance matrix into three components representing odd (surface), even (double-bounce) and diffuse (volume) scatterings. With these decomposition techniques, phenomena, such as canopy-ground interactions, randomness of orientation, and size of scatters can be examined from SAR data. In this study we applied the method recently reported by van Zyl (1992) to decompose averaged backscattering covariance matrices extracted from JPL SAR images over forest stands in Maine, USA. These stands are mostly mixed stands of coniferous and deciduous trees. Biomass data have been derived from field measurements of DBH and tree density using allometric equations. The interpretation of the decompositions and relationships with measured stand biomass are presented in this paper.

Sun, Guoqing↗
NASA NTRS2018

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We describe a variation of the Optimal Estimation (OE) method for greenhouse gas remote sensing retrievals using a singular value decomposition (SVD) and an uninformative prior. The SVD method is capable of discerning vertical information in column integrated absorption measurements. While traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We discuss the method, illustrate its use on an integrated path differential absorption CO2 lidar measurement model, and compare it to traditional optimal estimation using numerical simulations. In the absence of forward model error, the mathematics behind the SVD method guarantee it to be bias-free, which is confirmed by the numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. While the SVD approach can be used for most trace gas retrievals, it is particularly useful for situations where prior knowledge of the trace gas profile is poor. The SVD analysis is also useful in illustrating how vertical information is treated by the traditional OE approach.

Ramanathan, Anand K.↗
NASA NTRS1991

Vibration-dissociation coupling in nonequilibrium flows

The purpose of the study is to evaluate all important physical phenomena contributing to the vibration-dissociation problem by developing a numerical technique for the simulation of the vibrational relaxation and dissociation/recombination of N(2) in an isothermal box. A set of coupled master equations representing each vibrational level of N(2) is solved, and both V-V and V-T exchanges between molecules and atoms are included, along with the effects of dissociation and recombustion from/to each level. It is observed that the inverse range parameter has a strong impact on the vibrational exchange rates and relaxation. In a heating case, the lower vibrational levels relax through a series of equilibrium distributions, while in a cooling case, a population inversion is formed in the upper levels. Other cases of heating and cooling will be simulated, and together with this work, used to develop a matrix of nonequilibrium relaxation runs.

Landrum, D. Brian↗
NASA NTRS1987

Improved VAS regression soundings of mesoscale temperature structure observed during the 1982 atmospheric variability experiment

An Atmospheric Variability Experiment (AVE) was conducted over the central U.S. in the spring of 1982, collecting radiosonde date to verify mesoscale soundings from the VISSR Atmospheric Sounder (VAS) on the GOES satellite. Previously published VAS/AVE comparisons for the 6 March 1982 case found that the satellite retrievals scarcely detected a low level temperature inversion or a mid-tropospheric cold pool over a special mesoscale radiosonde verification network in north central Texas. The previously published regression and physical retrieval algorithms did not fully utilize VAS' sensitivity to important subsynoptic thermal features. Therefore, the 6 March 1982 case was reprocessed adding two enhancements to the VAS regression retrieval algorithm: (1) the regression matrix was determined using AVE profile data obtained in the region at asynoptic times, and (2) more optimistic signal-to-noise statistical conditioning factors were applied to the VAS temperature sounding channels. The new VAS soundings resolve more of the low level temperature inversion and mid-level cold pool. Most of the improvements stems from the utilization of asynoptic radiosonde observations at NWS sites. This case suggests that VAS regression soundings may require a ground-based asynoptic profiler network to bridge the gap between the synoptic radiosonde network and the high resolution geosynchronous satellite observations during the day.

Chesters, Dennis↗
NASA NTRS1986

Concurrent Cholesky factorization of positive definite banded Hermitian matrices

First, the Cholesky factorization is extended to cover uniformly partitioned banded positive definite matrices of rank n which may be real symmetric or Hermitian. Then, two stratagems are given for the use of the algorithm in concurrent machines where the number of processing elements is less than required to factor the matrix in as few serial steps as possible, and where uniformly high efficiency is expected from all processing elements. Expressions are given for the efficiency factor e appearing in the speed-up expression q = eN, and these are specialized for the N node hypercube machine as a function of partition size s, the number N of processing elements of the hypercube machine, and the cost mu of interelement transmission relative to computation. It is shown that the efficiency factor e is inversely proportional to mu/s, and that e is almost independent of N when N is large and mu/s = 0. The task is completed in n/s serial steps with no limit on n. The half bandwidth b of the matrix is 2 Ns.

Utku, S.↗
NASA NTRS1993

Dynamic Forms. Part 1: Functions

The formalism of dynamic forms is developed as a means for organizing and systematizing the design control systems. The formalism allows the designer to easily compute derivatives to various orders of large composite functions that occur in flight-control design. Such functions involve many function-of-a-function calls that may be nested to many levels. The component functions may be multiaxis, nonlinear, and they may include rotation transformations. A dynamic form is defined as a variable together with its time derivatives up to some fixed but arbitrary order. The variable may be a scalar, a vector, a matrix, a direction cosine matrix, Euler angles, or Euler parameters. Algorithms for standard elementary functions and operations of scalar dynamic forms are developed first. Then vector and matrix operations and transformations between parameterization of rotations are developed in the next level in the hierarchy. Commonly occurring algorithms in control-system design, including inversion of pure feedback systems, are developed in the third level. A large-angle, three-axis attitude servo and other examples are included to illustrate the effectiveness of the developed formalism. All algorithms were implemented in FORTRAN code. Practical experience shows that the proposed formalism may significantly improve the productivity of the design and coding process.

Meyer, George↗
NASA NTRS1997

A Portable MPI Implementation of the SPAI Preconditioner in ISIS++

A parallel MPI implementation of the Sparse Approximate Inverse (SPAI) preconditioner is described. SPAI has proven to be a highly effective preconditioner, and is inherently parallel because it computes columns (or rows) of the preconditioning matrix independently. However, there are several problems that must be addressed for an efficient MPI implementation: load balance, latency hiding, and the need for one-sided communication. The effectiveness, efficiency, and scaling behavior of our implementation will be shown for different platforms.

Barnard, Stephen T.↗
NASA NTRS2001

The Role of Magnetograms in Testing Models of Coronal Mass Ejections

A general review of current and future magnetographic capabilities is presented with emphasis on how they address observational predictions of CME models (e.g., the model/observation matrix developed in the SHINE 2000 workshop). Existing line-of-sight magnetographs can observe certain structural and kinematic properties of CME source regions when they are on the visible disk. These include polarity configuration, location and structure of the magnetic inversion line, proper motions of magnetic elements, and time variation of flux. Additional chromospheric magnetograms provide information on how the field spreads with height and, if properly sequenced with photospheric magnetograms, can determine whether flux is emerging or submerging. Vector magnetographs can provide information about currents and magnetic shear but will only begin to attain the spatial and temporal coverage appropriate to CMEs with the advent of the SOLIS vector spectromagnetograph and, hopefully, its copies at several widely spaced longitudes. Important additions to the SHINE 2000 matrix are required to define an observational program. These include spatial and temporal resolutions and extents as well as needed sensitivity.

Jones, Harrison P.↗
NASA NTRS1987

Recursive inverse kinematics for robot arms via Kalman filtering and Bryson-Frazier smoothing

This paper applies linear filtering and smoothing theory to solve recursively the inverse kinematics problem for serial multilink manipulators. This problem is to find a set of joint angles that achieve a prescribed tip position and/or orientation. A widely applicable numerical search solution is presented. The approach finds the minimum of a generalized distance between the desired and the actual manipulator tip position and/or orientation. Both a first-order steepest-descent gradient search and a second-order Newton-Raphson search are developed. The optimal relaxation factor required for the steepest descent method is computed recursively using an outward/inward procedure similar to those used typically for recursive inverse dynamics calculations. The second-order search requires evaluation of a gradient and an approximate Hessian. A Gauss-Markov approach is used to approximate the Hessian matrix in terms of products of first-order derivatives. This matrix is inverted recursively using a two-stage process of inward Kalman filtering followed by outward smoothing. This two-stage process is analogous to that recently developed by the author to solve by means of spatial filtering and smoothing the forward dynamics problem for serial manipulators.

Rodriguez, G.↗
NASA NTRS1988

Recursive inversion of externally defined linear systems

The approximate inversion of an internally unknown linear system, given by its impulse response sequence, by an inverse system having a finite impulse response, is considered. The recursive least squares procedure is shown to have an exact initialization, based on the triangular Toeplitz structure of the matrix involved. The proposed approach also suggests solutions to the problems of system identification and compensation.

Bach, Ralph E., Jr.↗
NASA NTRS2007

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗
NASA NTRS2010

An Atmospheric Science Observing System Simulation Experiment (OSSE) Environment

An atmospheric sounding mission starts with a wide range of concept designs involving measurement technologies, observing platforms, and observation scenarios. Observing system simulation experiment (OSSE) is a technical approach to evaluate the relative merits of mission and instrument concepts. At Jet Propulsion Laboratory (JPL), the OSSE team has developed an OSSE environment that allows atmospheric scientists to systematically explore a wide range of mission and instrument concepts and formulate a science traceability matrix with a quantitative science impact analysis. The OSSE environment virtually creates a multi-platform atmospheric sounding testbed (MAST) by integrating atmospheric phenomena models, forward modeling methods, and inverse modeling methods. The MAST performs OSSEs in four loosely coupled processes, observation scenario exploration, measurement quality exploration, measurement quality evaluation, and science impact analysis.

mission concepts↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2025

Properties of states near 𝐸 𝑥 = 6 MeV in 18 Ne through 17 F +𝑝 scattering

The rate of energy production in the hot-CNO cycle and breakout to the rapid-proton capture process in Type I X-ray bursts is strongly related to the 14 O⁡(𝛼,𝑝)⁢ 17 F reaction rate. The properties of states in 18 Ne near 𝐸 𝑥 = 6.1–6.3 MeV are important for understanding the 14 O⁡(𝛼,𝑝)⁢ 17 F reaction rate. In order to study 18 Ne resonances around this energy region, the RESOLUT radioactive-ion beam facility at Florida State University was used to perform 17 F⁡(𝑝,𝑝)⁢ 17 F elastic scattering on a polypropylene target under inverse kinematics. Scattered protons were detected in a silicon-strip detector array while recoiling 17 F ions were detected in coincidence in a gas ionization detector. An 𝑅-matrix analysis of measured cross sections was conducted along with a reanalysis of data from previous 17 F +𝑝 measurements. All the data analyzed are well described by a consistent set of parameters with a 1 − assignment for a state at 6.14(1) MeV. A second comparable solution is also found with a 3 − assignment for the 6.14(1) MeV state. There is no conclusive evidence supporting one of the two solutions over the other. The rates of the 14 O⁡(𝛼,p)⁢ 17 F reaction that are determined from the two solutions differ by up to an order of magnitude.

6 ≤ A ≤ 19↗
NASA NTRS1989

Recursive inversion of externally defined linear systems by FIR filters

The approximate inversion of an internally unknown linear system, given by its impulse response sequence, by an inverse system having a finite impulse response, is considered. The recursive least-squares procedure is shown to have an exact initialization, based on the triangular Toeplitz structure of the matrix involved. The proposed approach also suggests solutions to the problem of system identification and compensation.

Bach, Ralph E., Jr.↗
DOE OSTI2026

Gaussian processes for inferring parton distributions

The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.

hadronic spectroscopy↗