
about psd estimation by FFT
19 次查看(过去 30 天)
显示 更早的评论
Hi,
In topic "Power Spectral Density Estimates Using FFT (https://www.mathworks.com/help/signal/ug/power-spectral-density-estimates-using-fft.html)", why are squared dft values multipled by (1/(fs*N))?
fs = 1000;
t = 0:1/fs:1-1/fs;
x = cos(2*pi*100*t) + randn(size(t));
N = length(x);
xdft = fft(x);
xdft = xdft(1:N/2+1);
psdx = (1/(fs*N)) * abs(xdft).^2;
psdx(2:end-1) = 2*psdx(2:end-1);
freq = 0:fs/length(x):fs/2;
0 个评论
采纳的回答
Askic V
2022-11-10
It is because of Parseval's theorem, which expresses the energy of a signal in time-domain in terms of the average energy in its frequency components

In this case, n and k are integers, but actually these are samples i.e. discrete time intervals, so they can be considered as n*deltaT, in such case it is
deltaT is a sample time period i.e. deltaT = 1/Fs where Fs is a sampling frequency.
That is why 1/(fs*N)!
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Parametric Spectral Estimation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!