1

我需要帮助来找到 FFT 信号的最大幅度。

假设我对音频文件执行 FFT 并获得一列复数，我如何从 FFT 信号中提取最大幅度及其索引？我尝试使用“max”语法，但出现错误：??? 下标索引必须是实数正整数或逻辑数。

会很感激一些帮助..谢谢

这是我用过的代码

[wave,fs]=wavread('c scale fast.wav'); % read file into memory */
%sound(wave,fs); % see what it sounds like */
t=0:1/fs:(length(wave)-1)/fs; % and get sampling frequency */


figure(90);
          subplot(2,1,1)
          %plot(t,wave)
          plot(t,abs(wave))
          title('Wave File')
          ylabel('Amplitude')
          xlabel('Length (in seconds)')


L = length(wave);
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(wave,NFFT)/L;
f = fs/2*linspace(0,1,NFFT/2+1);


% Plot single-sided amplitude spectrum.
        subplot(2,1,2)
        plot(f,2*abs(Y(1:NFFT/2+1))) 
        title('Single-Sided Amplitude Spectrum of y(t)')
        xlabel('Frequency (Hz)')
        ylabel('|Y(f)|')

 A = max(Y)
4

1 回答 1

1

fft 返回一个复数。您应该使用绝对值来找到最大值：

[maxY,indexOfMaxY] = max(abs(Y));
于 2013-11-28T09:22:29.840 回答