calculate the frequency where psd is maximun
8 views (last 30 days)
Wayne King on 28 Sep 2011
Hi, one thing you can do is to obtain a periodogram and then calculate the average power in an interval around the maximum frequency and then compare that to the average power from 0 to the Nyquist.
Fs = 1e3;
t = 0:1/Fs:1-(1/Fs);
x = cos(2*pi*100*t)+randn(size(t));
psdest = psd(spectrum.periodogram,x,'Fs',1e3,'NFFT',length(x));
[mx,I] = max(psdest.Data);
relperc = ...
100*avgpower(psdest,[psdest.Frequencies(I-2) psdest.Frequencies(I+2)])/avgpower(psdest,[0 Fs/2])