0

Fourier transform is an estimator of frequency. Because its an estimator, there is always some uncertainty associated with my coefficients as described by the Fourier/Gabor limit. I'm wondering how to generate a probability density of frequency from the transformed data. My intuition is that I like what they do in physics with wave function probability amplitudes: modulus square the coefficients and normalize.

Edit: Fourier transform is NOT an estimator of frequency when integrating across time from -infinity to infinity. When you do not do this, you do not get distinct poles, but gaussian-like curves. First question, does it make sense or can you think of a time-limited fourier-transform in terms of probability? Given my time-limited fourier transform, what is the probability that frequency f is precsent?

Davey
  • 274
  • 1
  • 7
  • 3
    The Fourier transform is not an estimator of frequency. It is a transform from one domain to the related frequency domain (e.g. from time to frequency). There is no estimation. – Peter K. Oct 01 '22 at 19:44
  • 1
    As Peter said, we need to be a bit more nit-picky here (sorry!) as it's really important to understand what you want to estimate when discussing the limits of accuracy/precision of that estimate. Is it the frequency of a tone (or a discrete set of tones) or is it a spectrum? The PSD (what you would be most tempted to call "spectrum") inherently is a stochastic description of a random process in frequency domain. But it's not clear what you would mean with "probability density of frequency". For that you would first need to define "frequency of what", specifically! – Marcus Müller Oct 02 '22 at 00:49
  • 1
    If this concerns finite measurements, staring at DFT won't do much except for pure tones. Prefer time-frequency analysis. – OverLordGoldDragon Oct 02 '22 at 11:17
  • Appreciate the nit-pickiness actually. I am looking for probability density of tone frequencies with time-limited discrete measurements. If I normalized PSD so that area under curve = 1, would that be what I'm looking for? – Davey Oct 02 '22 at 15:54
  • I guess since my signal is transient I would do ESD instead? – Davey Oct 02 '22 at 17:40
  • So, let's formalize a signal model then. Can you write down your signal (for example, the sum of a sinusoid with white noise)? And write down which of the parameters of the model are random (and how they are distributed), and which of these you want to know the distribution of? – Marcus Müller Oct 04 '22 at 00:58
  • Instead of leading with incorrect statements of fact ("Fourier transform is an estimator of frequency"), please edit and revise your question to explain what you really want to do, and ask how you might do it? I.e., tell us the characteristics of the signal -- which, presumably, includes strong sinusoidal components -- and ask how you might estimate the frequencies of those components. – TimWescott Oct 07 '22 at 22:25

1 Answers1

1

The Power Spectral Density can be estimated using a number of methods.

  • The most straightforward, like you said, is taking the magnitude squared of the frequency coefficients: $\mathbf{PSD} = |X(k)|^2$
  • A better estimator can be obtained through Welch's method, which uses frequency averaging to increase the estimate's SNR.
  • Here is a great resource for proper PSD scaling.
Jdip
  • 5,980
  • 3
  • 7
  • 29
  • 1
    I've never known it as Welch's method. But this power spectrum estimation, in general, is not new. Sliding window, pick a good window (I've been using Gaussian of late but have also used Kaiser and Hamming and Hann), apply DFT to the windowed signal, magnitude square, pick out discrete peaks at integer index locations, apply quadratic interpolation to get a precise peak location. Precision is not the same as accuracy, but the expected accuracy can be determined in analysis, and there are different ways of doing that analysis to get an idea of how off some frequency estimate is. – robert bristow-johnson Oct 01 '22 at 21:31
  • 1
    Are you sure you want to "reduce the estimate's SNR"? Most of us have spent a lifetime struggling to increase SNR! :-) – Dilip Sarwate Oct 07 '22 at 18:47
  • Ha, I ALWAYS make that mistake, thank you for catching! Edited – Jdip Oct 07 '22 at 19:01