Power Spectral Density Estimates Using FFT
Show older comments
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
on 18 Jun 2015
1 vote
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.
Categories
Find more on Spectral Estimation in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!