Search NASA⌕ Search

SEARCH · Search NASA

Results for “periodogram”

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 19 records

A note on the use of coherent integration in periodogram analysis of MST radar signals

The effect of coherent integration on the periodogram method to estimate the power spectra of MST radar signals is examined. The spectrum estimate usually is biased, even when care is taken to reduce the aliasing effects. Due to this bias, the signal power for Doppler shifted signals is underestimated by as much as 4 dB. The use of coherent integration in reducing the effect of aliased power line harmonics is pointed out.

Rastogi, P. K.↗

Periodogram analysis of sunspot numbers and the relation with solar activities

The time series of average monthly sunspot numbers during 1900-1992 is studied by using power spectral analysis. This prediction method is used to study the sunspot periodicities relations between its, and with the other periodicities by solar activities. There are periodicities (between few days and 5 years) overwhelm on the mean solar cycle. ( 11 year cycle). These periodicities have the same relation with variations of solar constant and solar radiation reaching the Earth's atmosphere in the last solar cycle. These periods are related to the solar magnetic activity and to the modulation of solar features due to solar rotation.

Hady, Ahmed A.↗

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data

Detection of a periodic signal hidden in noise is frequently a goal in astronomical data analysis. This paper does not introduce a new detection technique, but instead studies the reliability and efficiency of detection with the most commonly used technique, the periodogram, in the case where the observation times are unevenly spaced. This choice was made because, of the methods in current use, it appears to have the simplest statistical behavior. A modification of the classical definition of the periodogram is necessary in order to retain the simple statistical behavior of the evenly spaced case. With this modification, periodogram analysis and least-squares fitting of sine waves to the data are exactly equivalent. Certain difficulties with the use of the periodogram are less important than commonly believed in the case of detection of strictly periodic signals. In addition, the standard method for mitigating these difficulties (tapering) can be used just as well if the sampling is uneven. An analysis of the statistical significance of signal detections is presented, with examples

Scargle, J. D.↗

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

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.

Komm, R. W.↗

Long-Term X-Ray Variability of Circinus X-1

We present an analysis of long term X-ray monitoring observations of Circinus X-1 (Cir X-1) made with four different instruments: Vela 5B, Ariel V ASM, Ginga ASM, and RXTE ASM, over the course of more than 30 years. We use Lomb-Scargle periodograms to search for the approx. 16.5 day orbital period of Cir X-1 in each of these data sets and from this derive a new orbital ephemeris based solely on X-ray measurements, which we compare to the previous ephemerides obtained from radio observations. We also use the Phase Dispersion Minimization (PDM) technique, as well as FFT analysis, to verify the periods obtained from periodograms. We obtain dynamic periodograms (both Lomb-Scargle and PDM) of Cir X-1 during the RXTE era, showing the period evolution of Cir X-1, and also displaying some unexplained discrete jumps in the location of the peak power.

CIRCINUS X-1↗

On the detection of other planetary systems by astrometric techniques

A quantitative method for astrometrically detecting perturbations induced in a star's motion by the presence of a planetary object is described. A periodogram is defined, wherein signals observed from a star show exactly periodic variations, which can be extracted from observational data using purely statistical methods. A detection threshold is defined for the frequency of occurrence of some detectable signal, e.g., the Nyquist frequency. Possible effects of a stellar orbital eccentricity and multiple companions are discussed, noting that assumption of a circular orbit assures the spectral purity of the signal described. The periodogram technique was applied to 12 yr of astrometric data from the U.S. Naval Observatory for three stars with low mass stellar companions. Periodic perturbations were confirmed. A comparison of the accuracy of different astrometric systems shows that the detection accuracy of a system is determined by the measurement accuracy and the number of observations, although the detection efficiency can be maximized by minimizing the number of data points for the case when observational errors are proportional to the square root of the number of data points. It is suggested that a space-based astrometric telescope is best suited to take advantage of the method.

Black, D. C.↗

Criteria and algorithms for spectrum parameterization of MST radar signals

The power spectra S(f) of MST radar signals contain useful information about the variance of refractivity fluctuations, the mean radial velocity, and the radial velocity variance in the atmosphere. When noise and other contaminating signals are absent, these quantities can be obtained directly from the zeroth, first and second order moments of the spectra. A step-by-step procedure is outlined that can be used effectively to reduce large amounts of MST radar data-averaged periodograms measured in range and time to a parameterized form. The parameters to which a periodogram can be reduced are outlined and the steps in the procedure, that may be followed selectively, to arrive at the final set of reduced parameters are given. Examples of the performance of the procedure are given and its use with other radars are commented on.

Rastogi, P. K.↗

Coherent lidar design and performance verification

The verification of LAWS beam alignment in space can be achieved by a measurement of heterodyne efficiency using the surface return. The crucial element is a direct detection signal that can be identified for each surface return. This should be satisfied for LAWS but will not be satisfied for descoped LAWS. The performance of algorithms for velocity estimation can be described with two basic parameters: the number of coherently detected photo-electrons per estimate and the number of independent signal samples per estimate. The average error of spectral domain velocity estimation algorithms are bounded by a new periodogram Cramer-Rao Bound. Comparison of the periodogram CRB with the exact CRB indicates a factor of two improvement in velocity accuracy is possible using non-spectral domain estimators. This improvement has been demonstrated with a maximum-likelihood estimator. The comparison of velocity estimation algorithms for 2 and 10 micron coherent lidar was performed by assuming all the system design parameters are fixed and the signal statistics are dominated by a 1 m/s rms wind fluctuation over the range gate. The beam alignment requirements for 2 micron are much more severe than for a 10 micron lidar. The effects of the random backscattered field on estimating the alignment error is a major problem for space based lidar operation, especially if the heterodyne efficiency cannot be estimated. For LAWS, the biggest science payoff would result from a short transmitted pulse, on the order of 0.5 microseconds instead of 3 microseconds. The numerically errors for simulation of laser propagation in the atmosphere have been determined as a joint project with the University of California, San Diego. Useful scaling laws were obtained for Kolmogorov atmospheric refractive turbulence and an atmospheric refractive turbulence characterized with an inner scale. This permits verification of the simulation procedure which is essential for the evaluation of the effects of refractive turbulence on coherent Doppler lidar systems. The analysis of 2 micron Doppler lidar data from Coherent Technologies, Inc. (CTI) has demonstrated many of the advantages of doppler lidar measurements of boundary layer winds. The effects of wind shear and wind turbulence over the pulse volume are probably the dominant source of the reduced performance. The effects of wind shear and wind turbulence on the statistical description of doppler lidar data has been derived and calculated.

Frehlich, Rod↗

Investigation of dispersion-relation-preserving scheme and spectral analysis methods for acoustic waves

Important characteristics of the aeroacoustic wave propagation are mostly encoded in their dispersion relations. Hence, a computational aeroacoustic (CAA) algorithm, which reasonably preserves these relations, was investigated. It was derived using an optimization procedure to ensure, that the numerical derivatives preserved the wave number and angular frequency of the differential terms in the linearized, 2-D Euler equations. Then, simulations were performed to validate the scheme and a compatible set of discretized boundary conditions. The computational results were found to agree favorably with the exact solutions. The boundary conditions were transparent to the outgoing waves, except when the disturbance source was close to a boundary. The time-domain data generated by such CAA solutions were often intractable until their spectra was analyzed. Therefore, the relative merits of three different methods were included in the study. For simple, periodic waves, the periodogram method produced better estimates of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window was more effective when used with the weighted-overlapped-segment-averaging and Blackman-Tukey methods gave better results than the periodogram method. Finally, it was demonstrated that the representation of time domain-data was significantly dependent on the particular spectral analysis method employed.

Vanel, Florence O.↗

Solar wind plasma periodicities observed at 1 AU by IMP 8

The IMP 8 spacecraft has been in Earth orbit since 1973, gathering plasma data over one complete 22-year solar cycle. These data are being examined to look for periodicities at time scales ranging from several hours to the entire span of the data set. A 1.3-year periodicity in the radial speed observed by IMP 8 and Voyager 2 has already been reported for the years from 1987 to 1993. The periodogram method, useful for unevenly sampled data such as the IMP 8 plasma data, has been used to search for other periods. It is interesting to note that the 13-year period is not present in the out-of-the-ecliptic component of the velocity (Vz), although a 1-year period is very obvious both visually and on the periodogram. Both components show a very strong peak associated with the 11-year solar cycle variation. This work will be extended to the thermal speed (a measure of the wind's temperature) and density, although the frequent correlations between these parameters and the velocity are expected to cause similar results. Additionally, the fine resolution data will be examined for shorter time periods than are visible using the hourly average data which are appropriate for longer periods. A comparison with periods observed at other spacecraft may also be made.

Paularena, K. I.↗

Coincident 1.3-year Periodicities in the ap Geomagnetic Index and the Solar Wind

Recent observations show an approximately 1.3-year period in the speed of the solar wind detected by the IMP 8 and Voyager 2 spacecraft. A similar period is also seen in the north-south (GSE) component of the magnetic field observed by IMP 8. Since both parameters are commonly used as input to models of geomagnetic activity, the 'ap' index (a measure of geomagnetic disturbance) is examined to look for this periodicity. The Lomb-Scargle periodogram method is used on the ap, plasma, and magnetic field data during the 1973-1994 time range. A dynamic FFT periodogram method is also used to analyze the ap data during this time, as well as to look for periods present between 1932 and 1972. A clear 1.3-year periodicity is present in the post-1986 data when the same period is observed in the plasma and field data. The V(2)B(zsm) and V(2)B(s) proxies for geomagnetic activity also show this periodicity. However, the southward (GSM) component of the magnetic field does not have a 1.3-year period, and neither do solar wind or ap data from 1973-1985. This demonstrates that the ap geomagnetic index can act as a proxy for solar wind periodicities at this time scale. Historic ap data are examined, and show that a similar periodicity in ap exists around 1942. Since auroral data show a 1.4-year periodicity, all these similar periods may result from a common underlying solar mechanism.

Paularena, K. I.↗

A procedure of multiple period searching in unequally spaced time-series with the Lomb-Scargle method

Periodogram analysis of unequally spaced time-series, as part of many biological rhythm investigations, is complicated. The mathematical framework is scattered over the literature, and the interpretation of results is often debatable. In this paper, we show that the Lomb-Scargle method is the appropriate tool for periodogram analysis of unequally spaced data. A unique procedure of multiple period searching is derived, facilitating the assessment of the various rhythms that may be present in a time-series. All relevant mathematical and statistical aspects are considered in detail, and much attention is given to the correct interpretation of results. The use of the procedure is illustrated by examples, and problems that may be encountered are discussed. It is argued that, when following the procedure of multiple period searching, we can even benefit from the unequal spacing of a time-series in biological rhythm research.

NASA Discipline Space Human Factors↗

Discovery of a 115 Day Orbital Period in the Ultraluminous X-ray Source NGC 5408 X-1

We report the detection of a 115 day periodicity in SWIFT/XRT monitoring data from the ultraluminous X-ray source (ULX) NGC 5408 X-1. Our o ngoing campaign samples its X-ray flux approximately twice weekly and has now achieved a temporal baseline of ti 485 days. Periodogram ana lysis reveals a significant periodicity with a period of 115.5 +/- 4 days. The modulation is detected with a significance of 3.2 x 10(exp -4) . The fractional modulation amplitude decreases with increasing e nergy, ranging from 0.13 +/- 0.02 above 1 keV to 0.24 +/- 0.02 below 1 keV. The shape of the profile evolves as well, becoming less sharply peaked at higher energies. The periodogram analysis is consistent wi th a periodic process, however, continued monitoring is required to c onfirm the coherent nature of the modulation. Spectral analysis indic ates that NGC 5408 X-1 can reach 0.3 - 10 keV luminosities of approxi mately 2 x 10 40 ergs/s . We suggest that, like the 62 day period of the ULX in M82 (X41.4-1-60), the periodicity detected in NGC 5408 X-1 represents the orbital period of the black hole binary containing the ULX. If this is true then the secondary can only be a giant or super giant star.

Strohmayer, Tod E.↗

Kepler Observations of Rapid Optical Variability in the BL Lac Object W2r192+42

We present the first Kepler monitoring of a strongly variable BL Lac, W2R1926+42. The light curve covers 181 days with approx. 0.2% errors, 30 minute sampling and >90% duty cycle, showing numerous delta-I/I > 25% flares over timescales as short as a day. The flux distribution is highly skewed and non-Gaussian. The variability shows a strong rms-flux correlation with the clearest evidence to date for non-linearity in this relation. We introduce a method to measure periodograms from the discrete autocorrelation function, an approach that may be well-suited to a wide range of Kepler data. The periodogram is not consistent with a simple power-law, but shows a flattening at frequencies below 7x10(exp -5) Hz. Simple models of the power spectrum, such as a broken power law, do not produce acceptable fits, indicating that the Kepler blazar light curve requires more sophisticated mathematical and physical descriptions than currently in use.

galaxies↗

Searching for Quasi-periodic Oscillations in Astrophysical Transients Using Gaussian Processes

Analyses of quasi-periodic oscillations(QPOs)are important to understanding the dynamic behavior in manyastrophysical objects during transient events like gamma-ray bursts, solarflares, magnetarflares, and fast radiobursts. Astrophysicists often search for QPOs with frequency-domain methods such as(Lomb–Scargle)periodograms, which generally assume power-law models plus some excess around the QPO frequency. Time-series data can alternatively be investigated directly in the time domain using Gaussian process(GP)regression.While GP regression is computationally expensive in the general case, the properties of astrophysical data andmodels allow fast likelihood strategies. Heteroscedasticity and nonstationarity in data have been shown to causebias in periodogram-based analyses. GPs can take account of these properties. Using GPs, we model QPOs as astochastic process on top of a deterministicflare shape. Using Bayesian inference, we demonstrate how to infer GPhyperparameters and assign them physical meaning, such as the QPO frequency. We also perform model selectionbetween QPOs and alternative models such as red noise and show that this can be used to reliablyfind QPOs. Thismethod is easily applicable to a variety of different astrophysical data sets. We demonstrate the use of this methodon a range of short transients: a gamma-ray burst, a magnetarflare, a magnetar giantflare, and simulated solarflare data.

Moritz Hubner↗

Optical observation of WRA 977 (3U 1233-62?)

UBV photometry of WRA 977 on 36 nights between January and July of 1974 shows that this object is active on a time scale of days at the 0.1 mag level, but that it remains quite constant during monitoring intervals lasting up to 1 hour. Periodogram analysis reveals no significant periodic variation in the brightness of this object.

Bord, D. J.↗