exa020105.m

信号的短时傅立叶分析

%---------------------------------------------------------------------------------------
%  exa020105.m,  for  example 2.1.5 and fig 2.1.5（a），（b）
%  说明窗函数宽度的选择对 STFT 时－频分辨率的影响；
%  注：在该程序中，用到了子程序 amagauss.m, fmlin.m 及tfstft.m, 
%      还要用到显示程序 tfrqview及其他程序，
%      所有这些程序请读者在如下的网站上下载:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因为该网站上的时频分析软件并没有列入MATLAB的工具箱，涉及到知识产权问题，
%      因此不能将其列入国内公开出版物上，但读者可以自由下载。
%----------------------------------------------------------------------------------------
clear;

sig=atoms(128,[50,.25,32,1;90,.25,32,1],0);

figure(1)
h=window(@hamming,55);
tfrstft(sig,1:128,128,h);

figure(2)
h=window(@hanning,13);
tfrstft(sig,1:128,128,h);