%*********************STFT函數***************************
%  輸入變量X表示是信號數據（列向量）；
%  t_interval是窗函數的移動步長（整數,一般為1）；
%  f_interval是做時的頻譜采樣間隔；
%  h為窗函數;

function st_ft=STFT(X,t_interval,fft_number,h)
N=length(X);M=length(h);                       %求信號數據&窗函數長度
n=N/t_interval;                                %由窗的移動步長，確定stft時間軸采樣數目
st_ft=zeros(fft_number,n);                     %申請stft數據空間;fft_number定頻率軸的采樣數目

X=[X;zeros(M,1)];                              %信號數據末端補零，使窗移到末端時可以取值

for i=1:t_interval:N;
    Xi=X(i:(i+M-1),1);
    Yi=Xi.*h;
    Yi=fft(Yi,fft_number);
    st_ft(:,i)=Yi;
    %st_ft(:,i)=fftshift(Yi);                    %將零頻，移到頻率的中心處
end