用MATLAB对WAV文件做FFT、IFFT及短时傅里叶变换code

%%%%%%傅里叶变换/逆变换/短时傅里叶变换%%%%%%

%[y,Fs]=wavread('C:\Users\HSF\Desktop\AgilentTechnologies.wav'); %读出信号，采样率和采样位数。

[y,Fs]=wavread('C:\Users\HSF\Desktop\sound\beiyou.wav');%读出信号，采样率和采样位数。

sound(y,Fs);

y=y(:,1);%我这里假设你的声音是双声道，我只取单声道作分析，如果你想分析另外一个声道，请改成y=y(:,2)

figure(1)

sigLength=length(y);

t=(0:length(y)-1)/Fs; %计算时间轴

subplot(3,1,1);plot(t,y);xlabel('Time(s)');%在第一个窗口画波形

grid on;

n=0:sigLength-1;

Y= fft(y,length(y));%在sigLength这个有限区间内做快速傅立叶变换

mag=abs(Y);

f=Fs*n/sigLength;

subplot(3,1,2);plot(f,abs(Y));

xlabel('Frequency(Hz)');

%在第一个窗口画率谱,

grid on;

subplot(3,1,3);

xifft=ifft(Y);

sound(xifft,Fs);

plot(t,xifft);xlabel('Time(s)')

grid on;

figure(2)

Nw=20;         %窗函数长window length

L=Nw/2;        %窗函数每次移动的样点数

Ts=round((sigLength-Nw)/L)+1;    %计算把数据y共分成多少段

nfft=512;       % FFT的长度

TF=zeros(Ts,nfft); %将存放三维图谱，先清零

for i=1:Ts

xw=y((i-1)*L+1:i*L+L);   %取一段数据

temp=fft(xw,nfft);      %FFT变换

%temp=fftshift(temp);    %频谱以0频为中心

TF(i,:)=temp;          %把谱图存放在TF中代表TF矩阵中的第i行

end

fn=(1:nfft)*Fs/nfft;

tn=(1:Ts)*Nw/2/Fs;

[T,F]=meshgrid(tn,fn);

mesh(T,F,abs(TF.'));  %三维绘图

axis tight;

title('三维时频图');

ylabel('Frequency(Hz)');

xlabel('Time(s)');

grid on;

仿真结果图如下：?