Removal of the phase noise in the autocorrelation estimates with data windowing

13th European Signal Processing Conference, EUSIPCO 2005; Antalya; Turkey; 4 September 2005 through 8 September 2005The sinusoidal frequency estimation from short data records based on Toeplitz autocorrelation (AC) matrix estimates suffer from phase noise. This effect becomes prominent especially when additive noise vanishes becoming a nuisance, that is at high signal-to-noise ratios (SNR). Based on both analytic derivation of the AC lag terms and simulation experiments, we show that data windowing can mitigate the limitations caused by the phase noise. Thus with proper windowing, the variance of the frequency estimate is no more limited by phase noise, but it continues to decrease linearly with the SNR. The cases of the Pisarenko frequency estimator and of MUSIC, both for the single sinusoid case, are analyzed in detail

First Results from COPSS: The CO Power Spectrum Survey

We present constraints on the abundance of carbon monoxide in the early universe from the CO Power Spectrum Survey. We utilize a data set collected between 2005 and 2008 using the Sunyaev–Zel'dovich Array (SZA), which was previously used to measure arcminute-scale fluctuations of the cosmic microwave background. This data set features observations of 44 fields, covering an effective area of 1.7 square degrees, over a frequency range of 27–35 GHz. Using the technique of intensity mapping, we are able to probe the CO(1–0) transition, with sensitivity to spatial modes between k = 0.5–2 h Mpc^(−1) over a range in redshift of z = 2.3–3.3, spanning a comoving volume of 3.6 × 10^6 h^(−3) Mpc^3. We demonstrate our ability to mitigate foregrounds, and present estimates of the impact of continuum sources on our measurement. We constrain the CO power spectrum to P_(CO) < 2.6 × 10^4 μK^2 (h^(−1) Mpc)^3, or Δ^2_(CO)(k = 1 h Mpc^(−1)) < 1.3 × 10^3 μK^2, at 95% confidence. This limit resides near optimistic predictions for the CO power spectrum. Under the assumption that CO emission is proportional to halo mass during bursts of active star formation, this corresponds to a limit on the ratio of CO(1–0) luminosity to host halo mass of A_(CO) < 1.2 × 10^(−5) L⊙_ M_⊙^(−1). Further assuming a Milky Way-like conversion factor between CO luminosity and molecular gas mass (α_(CO) = 4.3 M_⊙ (K km s^(−1) pc^(−2))^(−1)), we constrain the global density of molecular gas to ρ_(z~3) (M_H_2) ⩽ 2.8 x 10^8 M_☉ Mpc^(-3)

Characterization of the Crab Pulsar's Timing Noise

We present a power spectral analysis of the Crab pulsar's timing noise, mainly using radio measurements from Jodrell Bank taken over the period 1982-1989. The power spectral analysis is complicated by nonuniform data sampling and the presence of a steep red power spectrum that can distort power spectra measurement by causing severe power ``leakage''. We develop a simple windowing method for computing red noise power spectra of uniformly sampled data sets and test it on Monte Carlo generated sample realizations of red power-law noise. We generalize time-domain methods of generating power-law red noise with even integer spectral indices to the case of noninteger spectral indices. The Jodrell Bank pulse phase residuals are dense and smooth enough that an interpolation onto a uniform time series is possible. A windowed power spectrum is computed revealing a periodic or nearly periodic component with a period of about 568 days and a 1/f^3 power-law noise component with a noise strength of 1.24 +/- 0.067 10^{-16} cycles^2/sec^2 over the analysis frequency range 0.003 - 0.1 cycles/day. This result deviates from past analyses which characterized the pulse phase timing residuals as either 1/f^4 power-law noise or a quasiperiodic process. The analysis was checked using the Deeter polynomial method of power spectrum estimation that was developed for the case of nonuniform sampling, but has lower spectral resolution. The timing noise is consistent with a torque noise spectrum rising with analysis frequency as f implying blue torque noise, a result not predicted by current models of pulsar timing noise. If the periodic or nearly periodic component is due to a binary companion, we find a companion mass > 3.2 Earth masses.Comment: 53 pages, 9 figures, submitted to MNRAS, abstract condense

Spectral filtering as a method of visualising and removing striped artefacts in digital elevation data

Spectral filtering was compared with traditional mean spatial filters to assess their ability to identify and remove striped artefacts in digital elevation data. The techniques were applied to two datasets: a 100 m contour derived digital elevation model (DEM) of southern Norway and a 2 m LiDAR DSM of the Lake District, UK. Both datasets contained diagonal data artefacts that were found to propagate into subsequent terrain analysis. Spectral filtering used fast Fourier transformation (FFT) frequency data to identify these data artefacts in both datasets. These were removed from the data by applying a cut filter, prior to the inverse transform. Spectral filtering showed considerable advantages over mean spatial filters, when both the absolute and spatial distribution of elevation changes made were examined. Elevation changes from the spectral filtering were restricted to frequencies removed by the cut filter, were small in magnitude and consequently avoided any global smoothing. Spectral filtering was found to avoid the smoothing of kernel based data editing, and provided a more informative measure of data artefacts present in the FFT frequency domain. Artefacts were found to be heterogeneous through the surfaces, a result of their strong correlations with spatially autocorrelated variables: landcover and landsurface geometry. Spectral filtering performed better on the 100 m DEM, where signal and artefact were clearly distinguishable in the frequency data. Spectrally filtered digital elevation datasets were found to provide a superior and more precise representation of the landsurface and be a more appropriate dataset for any subsequent geomorphological applications

NIMBUS-7 SBUV (Solar Backscatter Ultraviolet) observations of solar UV spectral irradiance variations caused by solar rotation and active-region evolution for the period November 7, 1978 - November 1, 1980

Observations of temporal variations of the solar UV spectral irradiance over several days to a few weeks in the 160-400 nm wavelength range are presented. Larger 28-day variations and a second episode of 13-day variations occurred during the second year of measurements. The thirteen day periodicity is not a harmonic of the 28-day periodicity. The 13-day periodicity dominates certain episodes of solar activity while others are dominated by 28-day periods accompanied by a week 14-day harmonic. Techniques for removing noise and long-term trends are described. Time series analysis results are presented for the Si II lines near 182 nm, the Al I continuum in the 190 nm to 205 nm range, the Mg I continuum in the 210 nm to 250 nm range, the MgII H & K lines at 280 nm, the Mg I line at 285 nm, and the Ca II K & H lines at 393 and 397 nm

The estimation of geoacoustic properties from broadband acoustic data, focusing on instantaneous frequency techniques

The compressional wave velocity and attenuation of marine sediments are fundamental to marine science. In order to obtain reliable estimates of these parameters it is necessary to examine in situ acoustic data, which is generally broadband. A variety of techniques for estimating the compressional wave velocity and attenuation from broadband acoustic data are reviewed. The application of Instantaneous Frequency (IF) techniques to data collected from a normal-incidence chirp profiler is examined. For the datasets examined the best estimates of IF are obtained by dividing the chirp profile into a series of sections, estimating the IF of each trace in the section using the first moments of the Wigner Ville distribution, and stacking the resulting IF to obtain a composite IF for the section. As the datasets examined cover both gassy and saturated sediments, this is likely to be the optimum technique for chirp datasets collected from all sediment environments

Mitigating Internal Instrument Coupling for 21 cm Cosmology. II. A Method Demonstration with the Hydrogen Epoch of Reionization Array

We present a study of internal reflection and cross-coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted Epoch of Reionization (EoR) line-of-sight modes in the range 0.2 h −1 Mpc−1 < ${k}_{\parallel }$ < 0.5 h −1 Mpc−1. In particular, we find evidence for nonnegligible amounts of spectral structure in the raw autocorrelations that overlaps with the EoR window and is suggestive of complex instrumental effects. Through systematic modeling on a single night of data, we find we can recover these modes in the power spectrum down to the integrated noise floor, achieving a dynamic range in the EoR window of 106 in power (mK2 units) with respect to the bright galactic foreground signal. Future work with deeper integrations will help determine whether these systematics can continue to be mitigated down to EoR levels. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field

Modelling and extraction of fundamental frequency in speech signals

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University.One of the most important parameters of speech is the fundamental frequency of vibration of voiced sounds. The audio sensation of the fundamental frequency is known as the pitch. Depending on the tonal/non-tonal category of language, the fundamental frequency conveys intonation, pragmatics and meaning. In addition the fundamental frequency and intonation carry speaker gender, age, identity, speaking style and emotional state. Accurate estimation of the fundamental frequency is critically important for functioning of speech processing applications such as speech coding, speech recognition, speech synthesis and voice morphing. This thesis makes contributions to the development of accurate pitch estimation research in three distinct ways: (1) an investigation of the impact of the window length on pitch estimation error, (2) an investigation of the use of the higher order moments and (3) an investigation of an analysis-synthesis method for selection of the best pitch value among N proposed candidates. Experimental evaluations show that the length of the speech window has a major impact on the accuracy of pitch estimation. Depending on the similarity criteria and the order of the statistical moment a window length of 37 to 80 ms gives the least error. In order to avoid excessive delay as a consequence of using a longer window, a method is proposed ii where the current short window is concatenated with the previous frames to form a longer signal window for pitch extraction. The use of second order and higher order moments, and the magnitude difference function, as the similarity criteria were explored and compared. A novel method of calculation of moments is introduced where the signal is split, i.e. rectified, into positive and negative valued samples. The moments for the positive and negative parts of the signal are computed separately and combined. The new method of calculation of moments from positive and negative parts and the higher order criteria provide competitive results. A challenging issue in pitch estimation is the determination of the best candidate from N extrema of the similarity criteria. The analysis-synthesis method proposed in this thesis selects the pitch candidate that provides the best reproduction (synthesis) of the harmonic spectrum of the original speech. The synthesis method must be such that the distortion increases with the increasing error in the estimate of the fundamental frequency. To this end a new method of spectral synthesis is proposed using an estimate of the spectral envelop and harmonically spaced asymmetric Gaussian pulses as excitation. The N-best method provides consistent reduction in pitch estimation error. The methods described in this thesis result in a significant improvement in the pitch accuracy and outperform the benchmark YIN method

Markov speckle for efficient random bit generation

Optical speckle is commonly observed in measurements using coherent radiation. While lacking experimental validation, previous work has often assumed that speckle’s random spatial pattern follows a Markov process. Here, we present a derivation and experimental confirmation of conditions under which this assumption holds true. We demonstrate that a detected speckle field can be designed to obey the first-order Markov property by using a Cauchy attenuation mask to modulate scattered light. Creating Markov speckle enables the development of more accurate and efficient image post-processing algorithms, with applications including improved de-noising, segmentation and super-resolution. To show its versatility, we use the Cauchy mask to maximize the entropy of a detected speckle field with fixed average speckle size, allowing cryptographic applications to extract a maximum number of useful random bits from speckle images