aaaassss.m

短时傅立叶变换及其反变换---绝对可运行！

clear all;
clc;
fs=8000;%设定采样频率
N=512;
n=1:N;
t=n/fs;
x=wavread('D:\speech.wav');
time=(1:length(x))/fs;	% 时间轴向量
figure(1);
plot(time, x);
title('原始语音信号图');
grid;

N=512;
NFFT=512;
N1=round((1000-N)*2/N+1);
w=rectwin(N);
for i=1:N1;
    xt=x((i-1)*N/2+1:(i-1)*N/2+N);
    xt=xt.*w;                     %加矩形窗
    STFT=fft(xt,NFFT);
end

figure(2);
plot(20*log10(abs(STFT)));
title('STFT');
xlabel('频率');
ylabel('dB');
grid;
y=rectwin(N).*x(1:512);
figure(3);
plot(y);title('加窗语音');
grid;

ln=log(STFT);
figure(5);
plot(ln);
xlabel('频率(Hz)');
ylabel('对数谱');
title('正弦信号y=2*pi*10t对数谱');
%傅立叶反变换
z=ifft(STFT,512);
figure(4);
plot(z);
xlabel('频率');
ylabel('y');
title('通过IFF恢复的语音波形');
grid;