A computer program for power spectral analysis of unequally spaced points. Volume 1 - User's manual
Computer program for power spectral analysis of unequally spaced points - Vol. 1
SEARCH · Search NASA
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.
Computer program for power spectral analysis of unequally spaced points - Vol. 1
The mathematical background in spectral analysis as applied to geodetic applications is summarized. The resolution (cut-off frequency) of the GEOS 3 altimeter data is examined by determining the shortest wavelength (corresponding to the cut-off frequency) recoverable. The data from some 18 profiles are used. The total power (variance) in the sea surface topography with respect to the reference ellipsoid as well as with respect to the GEM-9 surface is computed. A fast inversion algorithm for matrices of simple and block Toeplitz matrices and its application to least squares collocation is explained. This algorithm yields a considerable gain in computer time and storage in comparison with conventional least squares collocation. Frequency domain least squares collocation techniques are also introduced and applied to estimating gravity anomalies from GEOS 3 altimeter data. These techniques substantially reduce the computer time and requirements in storage associated with the conventional least squares collocation. Numerical examples given demonstrate the efficiency and speed of these techniques.
A brief description of the newly developed Empirical Mode Decomposition (ENID) and Hilbert Spectral Analysis (HSA) method will be given. The decomposition is adaptive and can be applied to both nonlinear and nonstationary data. Example of the method applied to a sample earthquake record will be given. The results indicate those low frequency components, totally missed by the Fourier analysis, are clearly identified by the new method. Comparisons with Wavelet and window Fourier analysis show the new method offers much better temporal and frequency resolutions.
The topics of thermal anemometry and high-order spectral analyses were the subject of the technical training. Specifically, the objective of the technical training was to study: (i) the recently introduced constant voltage anemometer (CVA) for high-speed boundary layer; and (ii) newly developed high-order spectral analysis techniques (HOSA). Both CVA and HOSA are relevant tools for studies of boundary layer transition and stability.
Spectral mixture analysis was shown to be a powerful, multifaceted tool for analysis of multi- and hyper-spectral data. The essence of the first phase of the approach is to determine a set of image end-members that best account for the spectral variance in an image cube within a constrained, linear least squares mixing model. The selection of the image end-member is usually achieved using a priori knowledge and successive trial and error solutions to refine the total number and physical location of the end-members. However, in many situations a more objective method of determining these essential components is desired. The problem of image end-member determination was approached objectively by using the inherent variance of the data. Unlike purely statistical methods such as factor analysis, this approach derives solutions that conform to a physically realistic model.
We present a systematic spectral analysis with Suzaku of six AGNs detected in the Swift/BAT hard X-ray (15-200 keV) survey, Swift J0138.6-4001, J0255.2-0011, J0350.1-5019, J0505.7-2348, J0601.9-8636, and J1628.1-5145. This is considered to be a representative sample of new AGNs without X-ray spectral information before the BAT survey. We find that the 0.5-200 keV spectra of these sources can be uniformly fit with a base model consisting of heavily absorbed (log NH >23.5/sq cm) transmitted components, scattered lights, a reflection component, and an iron-K emission line. There are two distinct groups, three "new type" AGNs (including the two sources reported by Ueda et al. 2007) with an extremely small scattered fraction (f(sub scat) < 0:5%) and strong reflection component (R = omega/2pi > or equal to 0.8 where omega is the solid angle of the reflector), and three "classical type" ones with f(sub scat > 0.5% and R < or approx. 0.8. The spectral parameters suggest that the new type has an optically thick torus for Thomson scattering (N(sub H) approx. 10(exp 25)/sq cm) with a small opening angle theta approx. 20deg viewed in a rather face-on geometry, while the classical type has a thin torus (N(sub H) approx. 10(exp 23-24)/sq cm) with theta > or approx. 30deg. We infer that a significant number of new type AGNs with an edge-on view is missing in the current all-sky hard X-ray surveys. Subject headings: galaxies: active . gamma rays: observations . X-rays: galaxies . X-rays: general
Cape Kennedy vertical wind speed profiles measured by spectral analysis using fast Fourier transform method
Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.
Equatorial bottomside sinusoidal (BSS) irregularities have been studied by applying techniques of cross-correlation and spectral analysis to the Atmosphere Explorer data set. The phase of the cross-correlations of the plasma number density is discussed and the two drift velocity components observed using the retarding potential analyzer and ion drift meter on the satellite are discussed. Morphology is addressed, presenting the geographical distributions of the occurrence of BSS events for the equinoxes and solstices. Physical processes including the ion Larmor flux, interhemispheric plasma flows, and variations in the lower F region Pedersen conductivity are invoked to explain the findings.
For the spectral analysis of high-resolution and high signal-to-noise (SN) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. Aims. New Krivvii oscillator strengths for a large number of lines enable us to construct more detailed model atoms for our NLTEmodel-atmosphere calculations. This enables us to search for additional Kr lines in observed spectra and to improve Kr abundance determinations. Methods. We calculated Krivvii oscillator strengths to consider radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of Kr lines that are exhibited in high-resolution and high SN ultraviolet (UV)observations of the hot white dwarf RE 0503.
The maximum entropy method has been applied in the spectral analysis of high-energy cosmic-ray intensity during the large Forbush event of July 13, 1982. An oscillation with period of about 2 hours and amplitude of 1 to 3% was found to be present during the decrease phase. This oscillation can be related to a similar periodicity in the magnetospheric field. However, the variation was not observed at all neutron monitor stations. In the beginning of the recovery phase, the intensity oscillated with a period of about 10 hours and amplitude of 3%.
A method of surface acoustic wave (SAW) spectral analysis and digital integration that has been used previously in coherent CW laser work with CO2 lasers at 10.6 microns is described. Expressions are derived for the signal to noise ratio in the measured voltage spectrum with an approximation for the general case and rigorous treatment for the low signal case. The atmospheric backscatter data accumulated by the airborne LATAS (laser true airspeed) coherent laser radar system are analyzed.
Changes in respiratory variance revealed by power spectral analysis of the pulmonary impedance pneumogram can be used to detect and measure stresses directly or indirectly affecting human respiratory function. When gastrointestinal distress occurred during a series of 5 total-body exposures of 30 R at a rate of 1.5 R/min, it was accompanied by typical shifts in pulmonary impedance power spectra. These changes did not occur after protracted exposure of 250 R (30 R daily) at 1.5 R/hr that failed to cause radiation sickness. This system for quantitating respiratory effort can also be used to detect alterations in one's ability to perform under controlled exercise conditions.
Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.
Preliminary results of an investigation of magnetic fluctuations seen upstream of two interplanetary shocks are presented. The spectral analysis includes calculation of the normalized reduced magnetic helicity spectrum, the normalized reduced cross-helicity spectrum, and the Alfven ratio as discussed by Matthaeus and Goldstein (1982). Minimum variance methods are used to compute wave polarization as a function of frequency. The Taylor 'frozen in flow' hypothesis is assumed to convert frequencies to wave vectors. Some of the basic properties of the waves, including the probable mode of propagation in association with both quasi-parallel forward and reverse shocks, are described. A comparison with previous results on the generation of waves at interplanetary and planetary shocks is presented.
Automatic equipment for sampling and preparing gases by dual column gas chromatography for subsequent mass spectral analysis of carbon isotopes during methane metabolism
Prior studies suggest that postprandial hypotension in elderly persons may be due to defective sympathetic nervous system activation. We examined autonomic control of heart rate (HR) after a meal using spectral analysis of HR data in 13 old (89 +/- 6 years) and 7 young (24 +/- 4 years) subjects. Total spectral power, an index of overall HR variability, was calculated for the frequency band between 0.01 and 0.40 Hz. Relatively low-frequency power, associated with sympathetic nervous system and baroreflex activation, was calculated for the 0.01 to 0.15 Hz band. High-frequency power, representing parasympathetic influences on HR, was calculated for the 0.15 to 0.40 Hz band. Mean arterial blood pressure declined 27 +/- 8 mm Hg by 60 minutes after the meal in elderly subjects, compared with 9 +/- 8 mm Hg in young subjects (p less than or equal to 0.0001, young vs old). The mean change in low-frequency HR power from 30 to 50 minutes after the meal was +19.4 +/- 25.3 U in young subjects versus -0.1 +/- 1.5 U in old subjects (p less than or equal to 0.02). Mean change in total power was also greater in young (19.0 +/- 26.6 U) subjects compared with old subjects (0.0 +/- 1.6 U, p greater than or equal to 0.02). Mean ratio of low:high-frequency power increased 3.1 +/- 3.3 U in young subjects vs 0.5 +/- 2.7 U in old subjects (p less than or equal to 0.01). The increase in low-frequency HR power and in the low:high frequency band ratio in young subjects is consistent with sympathetic activation in the postprandial period.(ABSTRACT TRUNCATED AT 250 WORDS).
Charge-coupled devices (CCDs) are ideally suited for performing sampled-data transversal filtering operations in the analog domain. Two algorithms have been identified for performing spectral analysis in which the bulk of the computation can be performed in a CCD transversal filter; the chirp z-transform and the prime transform. CCD implementation of both these transform algorithms is presented together with performance data and applications.