短时傅里叶变换的MATLAB代码

function [Y F] = stft（X NFFT WINDOW NOVERLAP Fs）

% [Y F] = STFT（X NFFT WINDOW NOVERLAP Fs）
% Short-time Fourier Transform Analysis
%  
%     Y : stFT doamin output matrix （NFFT x L）
%     F : real frequencies corresponding to each FFT points
%     Fs : sampling freguency
%     NFFT : # of FFT points
%     WINDOW : window function
%     NOVERLAP : # of overlaped samples
%
%  by Taesu Kim     2003. 1. 26.

if nargin<5 Fs = 8000; end

WLEN = size（WINDOW1）;