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))));

如何从双面光谱绘制单侧光谱？

1 个答案: