matlab - 如何从双面频谱绘制出单边频谱？

Question

我的代码中有一个wav文件。我在音频信号上使用了FFT来进行STFT。基本上我正在做的是“实时频谱分析”​​。我从信号中得到的是双向频谱。

这里提供的代码可以帮助您理解我的意思。我建议您尝试一下。只需提取任何wav文件并将其放在Matlab / Octave目录中即可。变量Yres包含两侧频谱。

任何想法，如何使它单面？我敢肯定，这很简单，我很想念。

[y,fs]=wavread('UnchainMyHeart.wav');
t=linspace(0,length(y)/fs,length(y));
plot(t,y)

%Plotting my signal in the time domain
fftlen = 4096; 
segl =floor(0.05*fs); 
windowshift=segl/2; 
window=hann(segl); 
window=window.'; 

si=1; 
ei=segl; 
AOS= length(y)/windowshift - 1;
f1=figure;
f=0:1:fftlen-1;
f=f/(fftlen-1)*fs;

Ya=zeros(1,fftlen);
plot(f,Ya),axis([0 fs -90 20])
grid on 

n=0;
for m = 1:1:AOS
  y_a = y(si:ei);
  y_a= y_a.*window;
  Ya=fft(y_a, fftlen);
  n=n+1;
  if n==1
    Yres=abs(Ya);
  else
    Yres=Yres+abs(Ya);
  end

  if n=10
    Yres=Yres/10;
    n=0;
    drawnow; 
    figure(f1);
    plot(f, 20*log10(abs(Yres)));
    ylim([-90 20]);
    title("Spectrum of a audio signal");
    xlabel("f(Hz)");
    ylabel("dB");
    grid on;
  end

  si=si+windowshift;   
  ei=ei+windowshift; 
end     

Answer 1

您可以绘制已经拥有的频谱的一半:

plot(f(1:end/2), 20*log10(abs(Yres(1:end/2))));

或者，您也可以选择完整但集中的内容:

plot(fftshift(f), 20*log10(abs(fftshift(Yres))));