# plot fails with -Inf

조회 수: 14(최근 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.
Thanks.
%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표시숨기기 이전 댓글 수: 2
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 2020년 3월 23일
you can use Px(isinf(abs(Px)))=0; to set all -Inf values to 0
##### 댓글 수: 3표시숨기기 이전 댓글 수: 2
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.

댓글을 달려면 로그인하십시오.

R2019b

### Community Treasure Hunt

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

Start Hunting!