1 / 33

Introduction To Signal Processing & Data Analysis

The L aser I nterferometer G ravitational-Wave O bservatory. http://www.ligo.caltech.edu. Introduction To Signal Processing & Data Analysis. Gregory Mendell LIGO Hanford Observatory. Supported by the United States National Science Foundation. Fourier Transforms. Digital Data. x. ….

minh
Download Presentation

Introduction To Signal Processing & Data Analysis

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. The Laser Interferometer Gravitational-Wave Observatory http://www.ligo.caltech.edu Introduction To Signal Processing & Data Analysis Gregory Mendell LIGO Hanford Observatory Supported by the United States National Science Foundation LIGO-G1200759

  2. Fourier Transforms

  3. . . . . Digital Data . x . . … t j=0 j=1 j=2 j=3 j=4 j=N-1 T Sample rate, fs: Duration T, number of samples N, Nyquist frequency: Time index j: Frequency resolution: Frequency index k:

  4. Discrete Fourier Transforms(DFT) The DFT is of order N2. The Fast Fourier Transform (FFT) is a fast way of doing the DFT, of order Nlog2N.

  5. DFT Aliasing For integer m: Power outside this band is aliased into this band. Thus need to filter data before digitizing, or when changing fs, to prevent aliasing of unwanted power into this band.

  6. DFT Aliasing

  7. DFT Orthogonality

  8. Parseval’s Theorem

  9. Correlation Theorem

  10. Convolution Theorem

  11. Example: Calibration DARM_ERR • Measure open loop gain G, input unity gain to model • Extract gain of sensing function C=G/AD from model • Produce response function at time of the calibration, R=(1+G)/C • Now, to extrapolate for future times, monitor calibration lines in DARM_ERR error signal, plus any changes in gains. • Can then produce R at any later time t. • A photon calibrator, using radiation pressure, gives results consistent with the standard calibration.

  12. One-sided Power Spectral Density (PSD) Estimation Note that the absolute square of a Fourier Transform gives what we call “power”. A one-sided PSD is defined for positive frequencies (the factor of 2 counts the power from negative frequencies). The angle brackets, < >, indicate “average value”. Without the angle brackets, the above is called a periodogram. Thus, the PSD estimate is found by averaging periodograms. The other factors normalize the PSD so that the area under the PDS curves gives the RMS2 of the time domain data.

  13. Gaussian White Noise Parseval’s Theorem is used above.

  14. Power Spectral Density, Gaussian White Noise This ratio is constant, independent of fs. The square root of the PSD is the Amplitude Spectral Density. It has the units of sigma per root Hz. For gaussian white noise, the square root of the area under the PSD gives the RMS of time domain data.

  15. Amplitude and Phase of a spectral line If the product of f and T is an integer k, we call say this frequency is “bin centered”. For a sinusoidal signal with a “bin centered” frequency, all the power lies in one bin.

  16. DFT Bin Centered Frequency

  17. DFT Leakage For kappa not an integer get sinc function response: For a signal (or spectral disturbance) with a non- bin-centered frequency, power leaks out into the neighboring bins.

  18. DFT Leakage

  19. Windowing reduces leakage

  20. Corrections to PSD

  21. Example PSD Estimation

  22. For iLIGO and aLIGO

  23. PSD Statistics Rayleigh Distribution: Chi-squared for 2 degrees of freedom:

  24. Histograms, Real and Imaginary Part of DFT For gaussian noise, the above are also gaussian distributions.

  25. Rayleigh Distribution Histogram of the square root of the DFT power.

  26. Histogram DFT Power Result is a chi-squared distribution for 2 degrees of freedom.

  27. Chi-Squared Statistics A chi-squared variable with ν degrees of freedom is the sum of the squares of ν gaussian distributed variables with zero mean and unit variance.

  28. Maximum Likelihood Likelihood of getting data x for model h for Gaussian Noise: Chi-squared Minimize Chi-squared = Maximize the Likelihood

  29. Matched Filtering This is called the log likelihood; x is the data, h is a model of the signal. Maximizing over a single amplitude gives: This can be generalized to maximize over more parameters, giving the formulas for matched-filtering in terms of dimensionless templates.

  30. Frequentist Confidence Region Matlab simulation for signal: • Estimate parameters by maximizing likelihood or minimizing chi-squared • Injected signal with estimated parameters into many synthetic sets of noise. • Re-estimate the parameters from each injection LIGO-G050492-00-W Figure courtesy Anah Mourant, Caltech SURF student, 2003.

  31. Bayesian Analysis Bayes’ Theorem: Likelihood: Posterior Probability: Confidence Interval: Prior

  32. Bayesian Confidence Region Matlab simulation for signal: Uniform Prior: x Figures courtesy Anah Mourant, Caltech SURF student, 2003. LIGO-G050492-00-W

  33. The End

More Related