高斯窗的短时傅里叶变换变换程序

function stft=my_STFT（xfs）
% 
% clc;
% clear;
% 
% addpath（‘D:\MATLAB\R2008b\work\mmread\‘）;
% info=mmfileinfo（‘spilths2.wav‘）;
% [video audio] = mmread（‘spilths2.wav‘[][]‘disableVideo‘）;
% d=audio.data;
% fs=audio.rate;
% df=10;
% d=d（1:df:length（d）1）;
% fs=fs/df;                                                           
% [ba]=butter（10882/1102.5‘low‘）;
% x=filter（bad）;
% x=x（50:50+1024）;
% x=x/max（abs（x））+eps;        

% L=1.2*fs;
% x=x（1:L）;

L=length（x）;

w=round（80e-3*fs）;
win=gausswin（w2）;  % sliding window
% win=rectwin（w）;
F=[];c=0;
hw=ceil（w/2）;
flag=0; k=