在matlab中只绘制FFT的一半

Question

我想在一个图中显示幅度和频率。它只显示两半。 （我的声音长度是 1 秒！）

fs_orig = 44100;
nbits = 16;
N=44100; %amount of samples???
[x1,fs_orig,nbits]=wavread('Mysong.wav');
[B2,A2] = butter(4,[5875 6500]/(fs_orig/2),'bandpass');

dataIn = x1; %your music
dataOut2 = filter(B2, A2, dataIn);
data_pow2 = dataOut2.^2;

X_mags2 = abs(fft(data_pow2));

bin_vals = [0 : N-1];
fax_Hz = bin_vals*fs_orig/N;
N_2 = ceil(N/2);

plot(fax_Hz(1:N_2), 20*log10(X_mags2(1:N_2)))
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)');

现在它仍然显示两边，有人知道为什么吗？

Answer 1

您没有得到正确结果的原因是您混淆了 FFT 产生的总点数与您的采样频率（44.1 kHz）相比。

FFT 计算的默认点数等于输入信号的长度。因此，如果您的信号data_pow2 有 128 个点，那么输出也将有 128 个点。您可以通过为其指定附加参数来手动指定 FFT 中所需的点数：

B = fft(A, N);

这里，N 将是您希望 FFT 具有的点数。不指定这将假定您希望点的总数是输入的长度。因此，您实际需要做的是将 N 更改为信号的长度。因此，您需要这样做：

[x1,fs_orig,nbits]=wavread('Mysong.wav');
N = numel(x1);

声明N 之后你做了wavread。因此，如果你说你的信号是 1 秒，那么这应该意味着元素的总数应该是 44100。当你用 wavread 读入你的信号时，验证 N 是否等于 44100。