Search NASA⌕ Search

SEARCH · Search NASA

Results for “SVD”

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 109 records · Page 6

Progress Towards A Measurement of Neutrino Induced Charged Current Neutral Pion Production in the MicroBooNE Experiment

An analysis of MicroBooNE data with a signal of one muon, one neutral pion, and no charged pions is presented. Studying neutral pion production in the MicroBooNE detector provides an opportunity to better understand neutrino-argon interactions, and is crucial for future accelerator-based neutrino oscillation experiments. This analysis presents the progress towards the first measurement of the differential cross section for charged current (CC) $π^0$ production in neutrino-argon interactions. Using a dataset corresponding to about 7 × 10 20 protons on target (POT), we present an analysis which aims to measure the single differential cross sections as a function of the $π^0$ kinematic variables such as momentum and scattering angle. The Wiener-SVD technique for unfolding the measurement is presented and demonstrated using multiple generator predictions. A future iteration of this analysis will compare an unfolded data measurement to these models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DIAS: A Data-Informed Active Subspace Regularization Framework for Inverse Problems

This paper presents a regularization framework that aims to improve the fidelity of Tikhonov inverse solutions. At the heart of the framework is the data-informed regularization idea that only data-uninformed parameters need to be regularized, while the data-informed parameters, on which data and forward model are integrated, should remain untouched. We propose to employ the active subspace method to determine the data-informativeness of a parameter. The resulting framework is thus called a data-informed (DI) active subspace (DIAS) regularization. Four proposed DIAS variants are rigorously analyzed, shown to be robust with the regularization parameter and capable of avoiding polluting solution features informed by the data. They are thus well suited for problems with small or reasonably small noise corruptions in the data. Furthermore, the DIAS approaches can effectively reuse any Tikhonov regularization codes/libraries. Though they are readily applicable for nonlinear inverse problems, we focus on linear problems in this paper in order to gain insights into the framework. Various numerical results for linear inverse problems are presented to verify theoretical findings and to demonstrate advantages of the DIAS framework over the Tikhonov, truncated SVD, and the TSVD-based DI approaches.

Nguyen, Hai↗

A Deep Neural Network for Accurate and Robust Prediction of the Glass Transition Temperature of Polyhydroxyalkanoate Homo- and Copolymers

The purpose of this study was to develop a data-driven machine learning model to predict the performance properties of polyhydroxyalkanoates (PHAs), a group of biosourced polyesters featuring excellent performance, to guide future design and synthesis experiments. A deep neural network (DNN) machine learning model was built for predicting the glass transition temperature, Tg, of PHA homo- and copolymers. Molecular fingerprints were used to capture the structural and atomic information of PHA monomers. The other input variables included the molecular weight, the polydispersity index, and the percentage of each monomer in the homo- and copolymers. The results indicate that the DNN model achieves high accuracy in estimation of the glass transition temperature of PHAs. In addition, the symmetry of the DNN model is ensured by incorporating symmetry data in the training process. The DNN model achieved better performance than the support vector machine (SVD), a nonlinear ML model and least absolute shrinkage and selection operator (LASSO), a sparse linear regression model. The relative importance of factors affecting the DNN model prediction were analyzed. Sensitivity of the DNN model, including strategies to deal with missing data, were also investigated. Compared with commonly used machine learning models incorporating quantitative structure–property (QSPR) relationships, it does not require an explicit descriptor selection step but shows a comparable performance. The machine learning model framework can be readily extended to predict other properties.

quantitative structure–property relationship (QSPR↗

Thermal Radiation Transport with Tensor Trains

We present a novel tensor network algorithm to solve the time-dependent, gray thermal radiation transport equation. The method invokes a tensor train (TT) decomposition for the specific intensity. The efficiency of this approach is dictated by the rank of the decomposition. When the solution is “low rank,” the memory footprint of the specific intensity solution vector may be significantly compressed. The algorithm, following a step-then-truncate approach of a traditional discrete ordinates method, operates directly on the compressed state vector, thereby enabling large speedups for low-rank solutions. To achieve these speedups, we rely on a recently developed rounding approach based on the Gram-SVD. We detail how familiar S N algorithms for (gray) thermal transport can be mapped to this TT framework and present several numerical examples testing both the optically thick and thin regimes. The TT framework finds low-rank structure and supplies up to ≃60× speedups and ≃1000× compressions for problems demanding large angle counts, thereby enabling previously intractable SN calculations and supplying a promising avenue to mitigate ray effects.

79 ASTRONOMY AND ASTROPHYSICS↗

Acoustooptic linear algebra processors - Architectures, algorithms, and applications

Architectures, algorithms, and applications for systolic processors are described with attention to the realization of parallel algorithms on various optical systolic array processors. Systolic processors for matrices with special structure and matrices of general structure, and the realization of matrix-vector, matrix-matrix, and triple-matrix products and such architectures are described. Parallel algorithms for direct and indirect solutions to systems of linear algebraic equations and their implementation on optical systolic processors are detailed with attention to the pipelining and flow of data and operations. Parallel algorithms and their optical realization for LU and QR matrix decomposition are specifically detailed. These represent the fundamental operations necessary in the implementation of least squares, eigenvalue, and SVD solutions. Specific applications (e.g., the solution of partial differential equations, adaptive noise cancellation, and optimal control) are described to typify the use of matrix processors in modern advanced signal processing.

Casasent, D.↗

Landsat D Thematic Mapper image dimensionality reduction and geometric correction accuracy

To characterize and quantify the performance of the Landsat thematic mapper (TM), techniques for dimensionality reduction by linear transformation have been studied and evaluated and the accuracy of the correction of geometric errors in TM images analyzed. Theoretical evaluations and comparisons for existing methods for the design of linear transformation for dimensionality reduction are presented. These methods include the discrete Karhunen Loeve (KL) expansion, Multiple Discriminant Analysis (MDA), Thematic Mapper (TM)-Tasseled Cap Linear Transformation and Singular Value Decomposition (SVD). A unified approach to these design problems is presented in which each method involves optimizing an objective function with respect to the linear transformation matrix. From these studies, four modified methods are proposed. They are referred to as the Space Variant Linear Transformation, the KL Transform-MDA hybrid method, and the First and Second Version of the Weighted MDA method. The modifications involve the assignment of weights to classes to achieve improvements in the class conditional probability of error for classes with high weights. Experimental evaluations of the existing and proposed methods have been performed using the six reflective bands of the TM data. It is shown that in terms of probability of classification error and the percentage of the cumulative eigenvalues, the six reflective bands of the TM data require only a three dimensional feature space. It is shown experimentally as well that for the proposed methods, the classes with high weights have improvements in class conditional probability of error estimates as expected.

Ford, G. E.↗

Statistical analysis of effective singular values in matrix rank determination

A major problem in using SVD (singular-value decomposition) as a tool in determining the effective rank of a perturbed matrix is that of distinguishing between significantly small and significantly large singular values to the end, conference regions are derived for the perturbed singular values of matrices with noisy observation data. The analysis is based on the theories of perturbations of singular values and statistical significance test. Threshold bounds for perturbation due to finite-precision and i.i.d. random models are evaluated. In random models, the threshold bounds depend on the dimension of the matrix, the noisy variance, and predefined statistical level of significance. Results applied to the problem of determining the effective order of a linear autoregressive system from the approximate rank of a sample autocorrelation matrix are considered. Various numerical examples illustrating the usefulness of these bounds and comparisons to other previously known approaches are given.

Konstantinides, Konstantinos↗

Constraint elimination in dynamical systems

Large space structures (LSSs) and other dynamical systems of current interest are often extremely complex assemblies of rigid and flexible bodies subjected to kinematical constraints. A formulation is presented for the governing equations of constrained multibody systems via the application of singular value decomposition (SVD). The resulting equations of motion are shown to be of minimum dimension.

Singh, R. P.↗

On recursive least-squares filtering algorithms and implementations

In many real-time signal processing applications, fast and numerically stable algorithms for solving least-squares problems are necessary and important. In particular, under non-stationary conditions, these algorithms must be able to adapt themselves to reflect the changes in the system and take appropriate adjustments to achieve optimum performances. Among existing algorithms, the QR-decomposition (QRD)-based recursive least-squares (RLS) methods have been shown to be useful and effective for adaptive signal processing. In order to increase the speed of processing and achieve high throughput rate, many algorithms are being vectorized and/or pipelined to facilitate high degrees of parallelism. A time-recursive formulation of RLS filtering employing block QRD will be considered first. Several methods, including a new non-continuous windowing scheme based on selectively rejecting contaminated data, were investigated for adaptive processing. Based on systolic triarrays, many other forms of systolic arrays are shown to be capable of implementing different algorithms. Various updating and downdating systolic algorithms and architectures for RLS filtering are examined and compared in details, which include Householder reflector, Gram-Schmidt procedure, and Givens rotation. A unified approach encompassing existing square-root-free algorithms is also proposed. For the sinusoidal spectrum estimation problem, a judicious method of separating the noise from the signal is of great interest. Various truncated QR methods are proposed for this purpose and compared to the truncated SVD method. Computer simulations provided for detailed comparisons show the effectiveness of these methods. This thesis deals with fundamental issues of numerical stability, computational efficiency, adaptivity, and VLSI implementation for the RLS filtering problems. In all, various new and modified algorithms and architectures are proposed and analyzed; the significance of any of the new method depends crucially on specific application.

Hsieh, Shih-Fu↗

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗

An Assessment of the Impact of the 1997-98 El Nino on the Asian-Australian Monsoon

Using state-of-the-art satellite-gauge monthly rainfall estimate and optimally interpolated sea surface temperature (SST) data, we have assessed the 1997-98 AA-monsoon anomalies in terms of three basic causal factors: basin-scale SST, regional coupling, and internal variability. Singular Value Decomposition analyses of rainfall and SST are carried out globally over the entire tropics and regionally over the AA-monsoon domain. Contributions to monsoon rainfall predictability by various factors are evaluated from cumulative anomaly correlation with dominant regional SVD modes. Results reveal a dominant, large-scale monsoon-El Nino coupled mode with well-defined centers of action in the near-equatorial monsoon regions during the boreal summer and winter respectively. The observed 1997-98 AA-monsoon anomalies are found to be very complex with approximately 34% of the anomalies of the Asian (boreal) summer monsoon and 74% of the Australia (austral) monsoon attributable to basin-scale SST influence associated with El Nino. Regional coupled processes contribute an additional 19% and 10%, leaving about 47% and 16% due to internal dynamics for the boreal and austral monsoon respectively. For the boreal summer monsoon, it is noted that the highest monsoon predictability is not necessary associated with major El Nino events (e.g. 1997, 1982) but rather in non-El Nino years (e.g. 1980, 1988) when contributions from the regional coupled modes far exceed those from the basin-scale SST. The results suggest that in order to improve monsoon seasonal-to-interannual predictability, there is a need to exploit not only monsoon-El Nino relationship, but also intrinsic monsoon regional coupled processes.

Lau, K.-M.↗

Assessment of the 1997-1998 Asian Monsoon Anomalies

Using State-of-the-art satellite-gauge monthly rainfall estimate and optimally interpolated sea surface temperature (SST) data, we have assessed the 1997-98 Asian monsoon anomalies in terms of three basic causal factors: basin-scale SST, regional coupling, and internal variability. Singular Value Decomposition analysis of rainfall and SST are carried out globally over the entire tropics and regionally over the Asian monsoon domain. Contributions to monsoon rainfall predictability by various factors are evaluated from cumulative anomaly correlation with dominant regional SVD modes. Results reveal a dominant, large-scale monsoon-El Nino coupled mode with well-defined centers of action in the near-equatorial monsoon regions. it is noted that some subcontinental regions such as all-India, or arbitrarily chosen land regions over East Asia, while important socio-economically, are not near the centers of influence from El Nino, hence are not necessarily representative of the response of the entire monsoon region to El Nino. The observed 1997-98 Asian monsoon anomalies are found to be very complex with approximately 34% of the anomalies attributable to basin- scale SST influence associated with El Nino. Regional coupled processes contribute an additional 19%, leaving about 47% due to internal dynamics. Also noted is that the highest monsoon predictability is not necessary associated with major El Nino events (e.g. 1997, 1982) but rather in non-El Nino years (e.g. 1980, 1988) when contributions from the regional coupled modes far exceed those from the basin-scale SST. The results suggest that in order to improve monsoon seasonal-to-interannual predictability, there is a need to exploit not only monsoon-El Nino relationship, but also monsoon regional coupled processes and their modulation by long-term climate change.

Lau, William K.-M.↗

New Developments in Quaternion Estimation from Vector Observations

This paper contains a critical comparison of estimators minimizing Wahba's loss function. Some new results are presented for the QUaternion ESTimator (QUEST) and EStimators of the Optimal Quaternion (ESOQ and ESOQ2) to avoid the computational burden of sequential rotations in these algorithms. None of these methods is as robust in principle as Davenport's q method or the Singular Value Decomposition (SVD) method, which are significantly slower. Robustness is only an issue for measurements with widely differing accuracies, so the fastest estimators, the modified ESOQ and ESOQ2, are well suited to sensors that track multiple stars with comparable accuracies. More robust forms of ESOQ and ESOQ2 are developed that are intermediate in speed.

Markley, F. Landis↗

Comparison of Different Measurement Technologies for the In-Flight Assessment of Radiated Acoustic Intensity

A series of tests was planned and conducted in the Interior Noise Test Facility at Boeing Field, on the NASA Aries 757 flight research aircraft, and in the Structural Acoustic Loads and Transmission Facility at NASA Langley Research Center. These tests were designed to answer several questions concerning the use of array methods in flight. One focus of the tests was determining whether and to what extent array methods could be used to identify the effects of an acoustical treatment applied to a limited portion of an aircraft fuselage. Another focus of the tests was to verify that the arrays could be used to localize and quantify a known source purposely placed in front of the arrays. Thus the issues related to backside sources and flanking paths present in the complicated sound field were addressed during these tests. These issues were addressed through the use of reference transducers, both accelerometers mounted to the fuselage and microphones in the cabin, that were used to correlate the pressure holograms. measured by the microphone arrays using either SVD methods or partial coherence methods. This correlation analysis accepts only energy that is coherent with the sources sensed by the reference transducers, allowing a noise control engineer to only identify and study those vibratory sources of interest. The remainder of this paper will present a detailed description of the test setups that were used in this test sequence and typical results of the NAH/IBEM analysis used to reconstruct the sound fields. Also, a comparison of data obtained in the laboratory environments and during flights of the 757 aircraft will be made.

Klos, Jacob↗

Tracking Energy Flow Using a Volumetric Acoustic Intensity Imager (VAIM)

A new measurement device has been invented at the Naval Research Laboratory which images instantaneously the intensity vector throughout a three-dimensional volume nearly a meter on a side. The measurement device consists of a nearly transparent spherical array of 50 inexpensive microphones optimally positioned on an imaginary spherical surface of radius 0.2m. Front-end signal processing uses coherence analysis to produce multiple, phase-coherent holograms in the frequency domain each related to references located on suspect sound sources in an aircraft cabin. The analysis uses either SVD or Cholesky decomposition methods using ensemble averages of the cross-spectral density with the fixed references. The holograms are mathematically processed using spherical NAH (nearfield acoustical holography) to convert the measured pressure field into a vector intensity field in the volume of maximum radius 0.4 m centered on the sphere origin. The utility of this probe is evaluated in a detailed analysis of a recent in-flight experiment in cooperation with Boeing and NASA on NASA s Aries 757 aircraft. In this experiment the trim panels and insulation were removed over a section of the aircraft and the bare panels and windows were instrumented with accelerometers to use as references for the VAIM. Results show excellent success at locating and identifying the sources of interior noise in-flight in the frequency range of 0 to 1400 Hz. This work was supported by NASA and the Office of Naval Research.

Klos, Jacob↗

New and Improved T-wave Morphology Parameters to Differentiate Healthy Individuals from those with Cardiomyopathy and Coronary Artery Disease

We investigated the ability of several known as well as new ECG repolarization parameters to differentiate healthy individuals from patients with obstructive coronary artery disease (CAD) and cardiomyopathy (CM). Advanced high-fidelity 12-lead ECG tests (approx. 5-min supine) were first performed on a "training set" of 99 individuals: 33 with ischemic or dilated CM and low ejection fraction (EF less than 40%); 33 with catheterization-proven obstructive CAD but normal EF; and 33 age-/gender-matched healthy controls. The following multiple parameters of T-wave morphology (TWM) were derived via signal averaging and singular value decomposition (SVD, which yields 8 eigenvalues, rho(sub 1) greater than rho(sub 2)...greater than rho(sub 8) and studied for their retrospective accuracy in detecting underlying disease: 1) Principal component analysis ratio of the T wave (T-PCA) = 100*rho(sub 2)/rho(sub 1); 2) Relative T-wave residuum (rTWR) = 100* SIGMA (rho(sub 4)(sup 2) +...+ rho(sub 8)(sup 2)); 3) Modified complexity ratio of the T wave (T-mCR) = 100*SIGMA(rho(sub 3)(sup 2) +...+rho(sb 8) (sup 2)); and 4) Normalized 3-dimensional volume of the T wave (nTV) = 100*(rho(sub 2)*rho(sub 3)/rho(sub 1)(sup 2). All TWM parameters significantly differentiated CAD from controls (p less than 0.0001), and also CM from controls (p less than 0.0001). Retrospective areas under the ROC curve were 0.77, 0.81, 0.82, and 0.83 (CAD vs. controls) and 0.93, 0.89, 0.95 and 0.96 (CM vs. controls) for T-PCA, rTWR, T-mCR and nTV respectively. The newer TWM parameters (T-mCR and nTV) thus outperformed the more established parameters (T-PCA and rTWR), presumably by putting a greater emphasis on the third T-wave eigenvalue, which in most healthy subjects has little energy compared to the first two eigenvalues. Subsequent prospective analyses have also yielded similar results, such that we conclude that diagnostic differentiation of pathology from non-pathology may be especially aided by detecting the transference of energy from the first and second T-wave eigenvalues into the third T-wave eigenvalue.

Greco, E. C.↗

Three dimensional empirical mode decomposition analysis apparatus, method and article manufacture

An apparatus and method of analysis for three-dimensional (3D) physical phenomena. The physical phenomena may include any varying 3D phenomena such as time varying polar ice flows. A repesentation of the 3D phenomena is passed through a Hilbert transform to convert the data into complex form. A spatial variable is separated from the complex representation by producing a time based covariance matrix. The temporal parts of the principal components are produced by applying Singular Value Decomposition (SVD). Based on the rapidity with which the eigenvalues decay, the first 3-10 complex principal components (CPC) are selected for Empirical Mode Decomposition into intrinsic modes. The intrinsic modes produced are filtered in order to reconstruct the spatial part of the CPC. Finally, a filtered time series may be reconstructed from the first 3-10 filtered complex principal components.

Gloersen, Per↗

Characterization of the Inner Knot of the Crab: the Site of the Gamma-ray Flares?

One of the most intriguing recent discoveries has been the detection of powerful gamma-ray flares from the Crab Nebula. Such events, with a recurrence time of about once per year, can be so dramatic to make the system the brightest source in the gamma‐ray sky as occurred, e.g. in April 2011. These flares challenge our understanding of how pulsar wind nebulae work and defy current astrophysical models for particle acceleration. We present here our study of the inner knot located within a fraction of an arcsecond from the pulsar with the aim of characterizing the feature and asking if this might be the site of the origin of the gamma-ray flares. We took data using Keck, HST, and Chandra obtained as part of our multi‐wavelength campaign to identify the source of the enigmatic flares. We set an upper limit as to the gamma-ray flux from the knot. We also find that the dimensions, surface brightness, flux, etc. of the optical and infrared knot are all correlated with distance from the pulsar. This distance, in turn, varies with time. In addition to this most thorough characterization of the inner knot's properties, we examine the hypothesis that the knot may be the site of the flares by examining the knot separation versus the Fermi/LAT gamma-ray flux. Finally, as part of this research, we make use of a new approach employing singular value decomposition (SVD) for analyzing time series of images and compare the approach to more traditional methods. Our conclusions are only refined but not impacted by using the new approach.

Gamma-ray Flares↗