plot fails with -Inf

23 ビュー (過去 30 日間)
Jonathan 2020 年 3 月 23 日
コメント済み: Jonathan 2020 年 3 月 23 日
I need to plot a fft of an audio file. It appears that if any of the fft bins are zero, the conversion to dB produces an -Inf. Then the plot function does draw the figure but does not plot any data.
Is there an efficent way to trap the 0 value bins in the conversion to dB? For my work -Inf PSD would be equiv to zero.
The test audio file was a code generated audio file, 100Hz Sine for 10 seconds. Since only one fft bin will have a value, all the other bins are zero.
%Code Fragment:
[yS, Fs] = audioread(PathFN,'double');
L = length(yS);
Y = fft(yS,L);
% Version One
Px = Y.*conj(Y); % PSD Power of each freq components
Px = Px(1:L/2);
Px = 10*log10(Px); % Convert to dB
fVals = Fs * (0:(L/2)-1)/L;
plot(fVals, Px); % Plot
  3 件のコメント
Peng Li
Peng Li 2020 年 3 月 23 日
did you try my solution? If you add a small number like eps to Px, there won't be any -inf.


回答 (1 件)

Stijn Haenen
Stijn Haenen 2020 年 3 月 23 日
you can use Px(isinf(abs(Px)))=0; to set all -Inf values to 0
  3 件のコメント
Jonathan 2020 年 3 月 23 日
Yes, it was placed there and no just real values. aside from a typo, it does plot now. The single positive is at 100Hz as expected. However, it plots also at every 40Hz a spectra that is at/about -250db for the full range of Fs/2 except where the -Inf were. I didn't think it possible to have a negative PSD. Round off/precision errors from the program that generated the audio file? Thanks.





Community Treasure Hunt

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

Start Hunting!

Translated by