Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “singular value decomposition”

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

East Asian Summer Monsoon Rainfall: A Historical Perspective of the 1998 Flood over Yangtze River

One of the main factors that might have caused the disastrous flood in China during 1998 summer is long-term variations that include a trend indicating increasing monsoon rainfall over the Yangtze River Valley. China's 160-station monthly rainfall anomaly for the summers of 1955-98 is analyzed for exploring such long-term variations. Singular value decomposition (SVD) between the summer rainfall and the global sea surface temperature (SST) anomalies reveals that the rainfall over Yangtze River Valley is closely related to global and regional SST variabilities at both interannual and interdecadal timescales. SVD1 mode links the above normal rainfall condition in central China to an El Nino-like SSTA distribution, varying on interannual timescale modified by a trend during the period. SVD3 mode links positive rainfall anomaly in Yangtze River Valley to the warm SST anomaly in the subtropical western Pacific, varying on interannual timescales modified by interdecadal timescales. This link tends to be stronger when the Nino3 area becomes colder and the western subtropical Pacific becomes warmer. The 1998 summer is a transition season when the 1997/98 El Nino event was in its decaying phase, and the SST in the Nino3 area emerged below normal anomaly while the subtropical western Pacific SST above normal. Thus, the first and third SVD modes become dominant in 1998 summer, favoring more Asian summer monsoon rainfall over the Yangtze River Valley.

Weng, H.-Y.↗
NASA NTRS1999

How to Estimate Attitude from Vector Observations

The most robust estimators minimizing Wahba's loss function are Davenport's q method and the Singular Value Decomposition (SVD) method. The q method is faster than the SVD method with three or more measurements. The other algorithms are less robust since they solve the characteristic polynomial equation to find the maximum eigenvalue of Davenport's K matrix. They are only preferable when speed or processor power is an important consideration. Of these, Fast Optimal Attitude Matrix (FOAM) is the most robust and faster than the q method. Robustness is only an issue for measurements with widely differing accuracies, so the fastest algorithms, Quaternion ESTimator (QUEST), EStimator of the Optimal Quaternion (ESOQ), and ESOQ2, are well suited to star sensor applications.

Markley, F. Landis↗
NASA NTRS2000

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↗
NASA NTRS1999

A Method to Solve Interior and Exterior Camera Calibration Parameters for Image Resection

An iterative method is presented to solve the internal and external camera calibration parameters, given model target points and their images from one or more camera locations. The direct linear transform formulation was used to obtain a guess for the iterative method, and herein lies one of the strengths of the present method. In all test cases, the method converged to the correct solution. In general, an overdetermined system of nonlinear equations is solved in the least-squares sense. The iterative method presented is based on Newton-Raphson for solving systems of nonlinear algebraic equations. The Jacobian is analytically derived and the pseudo-inverse of the Jacobian is obtained by singular value decomposition.

Samtaney, Ravi↗
NASA NTRS2001

Nonstationary Dynamics Data Analysis with Wavelet-SVD Filtering

Nonstationary time-frequency analysis is used for identification and classification of aeroelastic and aeroservoelastic dynamics. Time-frequency multiscale wavelet processing generates discrete energy density distributions. The distributions are processed using the singular value decomposition (SVD). Discrete density functions derived from the SVD generate moments that detect the principal features in the data. The SVD standard basis vectors are applied and then compared with a transformed-SVD, or TSVD, which reduces the number of features into more compact energy density concentrations. Finally, from the feature extraction, wavelet-based modal parameter estimation is applied.

Brenner, Marty↗
NASA NTRS1995

Model-Based Obstacle Detention From Image Sequences

This paper proposes a method for detecting obstacles on a runway from monocular image sequences. The surface of the runway is modeled as a plane and the model flow field corresponding to the runway is described by 8 coefficients. The set of 8 coefficients describing the initial model flow field is given by the data from Inertial Navigation System (INS). The uncertainties of this initial model flow field are estimated and used to obtain the accurate model flow field. The residual flow field after warping images by the accurate model flow field is computed by solving over determined gradient-based optical flow equations using a singular value decomposition (SVD). This SVD gives us a new insight of the uncertainties inherent in the optical flow computation. Those pixels with large residual flow vectors are considered as obstacles. Experimental results for two real image sequences are presented.

Sull, Sanghoon↗
NASA NTRS2004

Trend and Variability of China Precipitation in Spring and Summer: Linkage to Sea Surface Temperatures

Observational records in the past 50 years show an upward trend of boreal-summer precipitation over central eastern China and a downward trend over northern China. During boreal spring, the trend is upward over southeastern China and downward over central eastern China. This study explores the forcing mechanism of these trends in association with the global sea-surface temperature (SST) variations on the interannual and inter-decadal timescales. Results based on Singular Value Decomposition analyses (SVD) show that the interannual variability of China precipitation in boreal spring and summer can be well defined by two centers of actions for each season, which are co-varying with two interannual modes of SSTs. The first SVD modes of precipitation in spring and summer, which are centered in southeastern China and northern China, respectively, are linked to an ENSO-like mode of SSTs. The second SVD modes of precipitation in both seasons are confined to central eastern China, and are primarily linked to SST variations over the warm pool and Indian Ocean. Features of the anomalous 850-hPa winds and 700-Wa geopotential height corresponding to these modes support a physical mechanism that explains the causal links between the modal variations of precipitation and SSTs. On the decadal and longer timescale, similar causal links are found between the same modes of precipitation and SSTs, except for the case of springtime precipitation over central eastern China. For this case, while the interannual mode of precipitation is positively correlated with the interannual variations of SSTs over the warm pool and Indian Ocean; the inter-decadal mode is negatively correlated with a different SST mode, the North Pacific mode. The later is responsible for the observed downward trend of springtime precipitation over central eastern China. For all other cases, both the interannual and inter-decadal variations of precipitation can be explained by the same mode of SSTs. The upward trend of springtime precipitation over southeastern China and downward trend of summertime precipitation over northern China are attributable to the warming trend of the ENSO-like mode. The recent frequent summertime floods over central eastern China are linked to the warming trend of SSTs over the warm pool and Indian Ocean.

Yang, Fanglin↗
NASA NTRS2001

Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

NASA Discipline Plant Biology↗
NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗
NASA NTRS2005

SVD analysis of Aura TES spectral residuals

Singular Value Decomposition (SVD) analysis is both a powerful diagnostic tool and an effective method of noise filtering. We present the results of an SVD analysis of an ensemble of spectral residuals acquired in September 2004 from a 16-orbit Aura Tropospheric Emission Spectrometer (TES) Global Survey and compare them to alternative methods such as zonal averages. In particular, the technique highlights issues such as the orbital variation of instrument response and incompletely modeled effects of surface emissivity and atmospheric composition.

spectra↗
NASA NTRS2007

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↗
NASA NTRS2007

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

Multi-channel System for Beat to Beat QT Interval Variability and its Use in Screening for Coronary Artery Disease and Cardiomyopathy

We investigated the ability of beat-to-beat QT interval variability (QTV) and related parameters to differentiate healthy individuals from patients with obstructive coronary artery disease (CAD) and cardiomyopathy (CM). For this purpose we developed a PC-based ECG software program that in real time, acquires, analyzes and displays QTV in each of the eight independent channels that constitute the 12-lead conventional ECG. The system also analyzes and displays the QTV from QT interval signals that are derived from multiple channels and from singular value decomposition (SVD) to substantially reduce the effect of noise and other artifacts on the QTV results. It also provides other useful SVD-related parameters such as the normalized 3-dimensional volume of the T wave (nTV) = 100*(rho(sub 2)*rho(sub 3)rho(sub 1^2). 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. All QTV parameters that were studied for their accuracy in detecting CM and CAD significantly differentiated both CM and CAD from controls (p less than 0.0001). Retrospective areas under the ROC curve (AUC) of SDNN-QTV, rmsSD-QTV, and QTV Index (QTVI) for CM vs. controls in the lead V5 were 0.85, 0.90, and 0.99, respectively, and those for CAD vs. controls in the lead II were 0.82, 0.82, and 0.89. Other advanced ECG parameters, such as HFQRS RAZ score, LF Lomb of RRV or QRS-T angle, differentiated both CM and CAD from controls less significantly, with the respective AUC values of 0.89, 0.88 and 0.98 for CM vs. controls, and 0.73, 0.71 and 0.80 for CAD vs. controls. QTV parameters (especially QTVI, which is QTV as indexed to RRV) were, diagnostically speaking, amongst the best performing of the advanced ECG techniques studied thus far.

Starc, V.↗
NASA NTRS2007

Construction and Use of Resting 12-Lead High Fidelity ECG "SuperScores" in Screening for Heart Disease

We investigated the accuracy of several conventional and advanced resting ECG parameters for identifying 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. Multiple conventional and advanced ECG parameters were studied for their individual and combined retrospective accuracies in detecting underlying disease, the advanced parameters falling within the following categories: 1) Signal averaged ECG, including 12-lead high frequency QRS (150-250 Hz) plus multiple filtered and unfiltered parameters from the derived Frank leads; 2) 12-lead P, QRS and T-wave morphology via singular value decomposition (SVD) plus signal averaging; 3) Multichannel (12-lead, derived Frank lead, SVD lead) beat-to-beat QT interval variability; 4) Spatial ventricular gradient (and gradient component) variability; and 5) Heart rate variability. Several multiparameter ECG SuperScores were derivable, using stepwise and then generalized additive logistic modeling, that each had 100% retrospective accuracy in detecting underlying CM or CAD. The performance of these same SuperScores was then prospectively evaluated using a test set of another 120 individuals (40 new individuals in each of the CM, CAD and control groups, respectively). All 12-lead ECG SuperScores retrospectively generated for CM continued to perform well in prospectively identifying CM (i.e., areas under the ROC curve greater than 0.95), with one such score (containing just 4 components) maintaining 100% prospective accuracy. SuperScores retrospectively generated for CAD performed somewhat less accurately, with prospective areas under the ROC curve typically in the 0.90-0.95 range. We conclude that resting 12-lead high-fidelity ECG employing and combining the results of several advanced ECG software techniques shows great promise as a rapid and inexpensive tool for screening of heart disease.

Schlegel, T. T.↗
NASA NTRS2007

Wavefront Control Toolbox for James Webb Space Telescope Testbed

We have developed a Matlab toolbox for wavefront control of optical systems. We have applied this toolbox to the optical models of James Webb Space Telescope (JWST) in general and to the JWST Testbed Telescope (TBT) in particular, implementing both unconstrained and constrained wavefront optimization to correct for possible misalignments present on the segmented primary mirror or the monolithic secondary mirror. The optical models implemented in Zemax optical design program and information is exchanged between Matlab and Zemax via the Dynamic Data Exchange (DDE) interface. The model configuration is managed using the XML protocol. The optimization algorithm uses influence functions for each adjustable degree of freedom of the optical mode. The iterative and non-iterative algorithms have been developed to converge to a local minimum of the root-mean-square (rms) of wavefront error using singular value decomposition technique of the control matrix of influence functions. The toolkit is highly modular and allows the user to choose control strategies for the degrees of freedom to be adjusted on a given iteration and wavefront convergence criterion. As the influence functions are nonlinear over the control parameter space, the toolkit also allows for trade-offs between frequency of updating the local influence functions and execution speed. The functionality of the toolbox and the validity of the underlying algorithms have been verified through extensive simulations.

Shiri, Ron↗
NASA NTRS2006

Realtime Multichannel System for Beat to Beat QT Interval Variability

The measurement of beat-to-beat QT interval variability (QTV) shows clinical promise for identifying several types of cardiac pathology. However, until now, there has been no device capable of displaying, in real time on a beattobeat basis, changes in QTV in all 12 conventional leads in a continuously monitored patient. While several software programs have been designed to analyze QTV, heretofore, such programs have all involved only a few channels (at most) and/or have required laborious user interaction or offline calculations and postprocessing, limiting their clinical utility. This paper describes a PC-based ECG software program that in real time, acquires, analyzes and displays QTV and also PQ interval variability (PQV) in each of the eight independent channels that constitute the 12lead conventional ECG. The system also processes certain related signals that are derived from singular value decomposition and that help to reduce the overall effects of noise on the realtime QTV and PQV results.

Starc, Vito↗
NASA NTRS2008

An Automated SVD for Alignment and Control of James Webb Space Telescope Mirrors

The James Webb Space Telescope (JWST) is a three-mirror anastigmatic telescope. The alignment of the segmented primary and secondary mirrors in the wavefront sensing and control process involves a series of actuators to control the six degrees-of-freedom motion on each surface in addition to the radius of curvature. The control matrix developed from the alignment parameters is over-determined and singular value decomposition (SVD) method is used to solve it in the least square sense. An automated SVD scheme has been developed to identify the most contributing modes in a typical alignment process and reduce the impact of error-prone modes from the control process.

Shiri, Sharam↗
NASA NTRS2008

Fast Plasma Instrument for MMS: Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. The Dual Electron Spectrometer (DES) of the Fast Plasma Instrument (FPI) for MMS meets these demanding requirements by acquiring the electron velocity distribution functions (VDFs) for the full sky with high-resolution angular measurements every 30 ms. This will provide unprecedented access to electron scale dynamics within the reconnection diffusion region. The DES consists of eight half-top-hat energy analyzers. Each analyzer has a 6 deg. x 11.25 deg. Full-sky coverage is achieved by electrostatically stepping the FOV of each of the eight sensors through four discrete deflection look directions. Data compression and burst memory management will provide approximately 30 minutes of high time resolution data during each orbit of the four MMS spacecraft. Each spacecraft will intelligently downlink the data sequences that contain the greatest amount of temporal structure. Here we present the results of a simulation of the DES analyzer measurements, data compression and decompression, as well as ground-based analysis using as a seed re-processed Cluster/PEACE electron measurements. The Cluster/PEACE electron measurements have been reprocessed through virtual DES analyzers with their proper geometrical, energy, and timing scale factors and re-mapped via interpolation to the DES angular and energy phase-space sampling measurements. The results of the simulated DES measurements are analyzed and the full moments of the simulated VDFs are compared with those obtained from the Cluster/PEACE spectrometer using a standard quadrature moment, a newly implemented spectral spherical harmonic method, and a singular value decomposition method. Our preliminary moment calculations show a remarkable agreement within the uncertainties of the measurements, with the results obtained by the Cluster/PEACE electron spectrometers. The data analyzed was selected because it represented a potential reconnection event as currently published.

Figueroa-Vinas, Adolfo↗