Power Spectral Density Estimates Using FFT

3 views (last 30 days)
Yigao Shao
Yigao Shao on 21 May 2015
Answered: Jeremy on 18 Jun 2015
In the code example, psdx = (1/(Fs*N)) * abs(xdft).^2;
why isn't it instead (Fs/N) such that psdx = (Fs/N) * abs(xdft).^2? Isn't the period equal to the length of the signal, which is N, multiplied by the sampling interval, which is 1/Fs ?

Answers (1)

Jeremy
Jeremy on 18 Jun 2015
There are different ways of looking at it but the PSD is often written as:
(fft*conj(fft))/(bin_width*NFFT^2)
(fft*conj(fft)) = abs(fft) and bin_width = Fs/NFFT, so you can see how it reduces down.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!