利用matlab对音频做FFT进行频谱分析

利用matlab对音频的左右声道做频谱分析。

syms y fs;

syms in_put;
in_put= 'C:\Users\getao\Desktop\立体声.wav'

%in_put 为音频文件所在的位置，.wav文件为文件的名称。

[y,fs] = audioread(in_put);

info=audioinfo(in_put)  ;
sound(y,fs);
T=1/fs; 
t=(0:length(y)-1)*T;
f=(0:length(y)-1)*fs/length(y);
figure(1);
yz=y(:,1);

     %提取音频信号的左声道数据
subplot(2,2,1);
plot(t,yz);
title('原始信号时域（左声道）');
xlabel('时间');
ylabel('振幅');

yz2=y(:,2);

   %提取音频信号的右声道数据
subplot(2,2,2);
plot(t,yz2);
title('原始信号时域（右声道）');
xlabel('时间');
ylabel('振幅');

subplot(2,2,3);
n=length(yz);
y1=fft(yz,n); 
F=fs/length(yz);

plot(f,abs(y1));
title('原始信号频谱（左声道）');
xlabel('F(Hz)');
ylabel('H(jw)');

subplot(2,2