# 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.
Thanks.
%Code Fragment:
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.

サインインしてコメントする。

### カテゴリ

Find more on Surface and Mesh Plots in Help Center and File Exchange

R2019b

### Community Treasure Hunt

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

Start Hunting!

Translated by