A message pops up in Matlab when I use spectrum.periodogram to find the median frequency of a signal, saying that it is not recommended. Why is this? I also get the same message for spectrum.pwelch
psdest = psd(spectrum.periodogram,x,'Fs',1000,'NFFT',length(x)); normcumsumpsd = cumsum(psdest.Data)./sum(psdest.Data); Ind = find(normcumsumpsd <=0.5,1,'last'); fprintf('Median frequency is %2.3f Hz\n',psdest.Frequencies(Ind));
After lots of research I still don't understand the output of the psd well enough so that I can susbstitute spectrum.periodgram. Normally I find my PSD using:
[Pxx,Fx] = pwelch(s,[],[],[],fs);plot(Fx,10*log10(Pxx))
Ideally I would like to find the median frequency by manipulating the Pxx and Fx values but I am struggling to relate Pxx directly to psdest.Data. I would be grateful if someone could point me in the right direction. Thank you in advance!
Best Answer