%---------------------------------------------------------------------------------------
%  exa020503.m,  for  example 2.5.3 and fig 2.5.1,同時也給出圖 2.6.1(a)
%  說明在Gabor變換中，在臨界抽樣的情況下，若窗函數(shù)為高斯窗，則對偶窗是非連續(xù)的，
%                     不再具有高斯窗的形狀。該程序同時給出了一個信號的Gabor變換；
%  注：在該程序中，用到了子程序  fmlin.m 及tfgabor.m,
%      還要用到顯示程序 tfrqview及其他程序，
%      所有這些程序請讀者在如下的網(wǎng)站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因為該網(wǎng)站上的時頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識產權問題，
%      因此不能將其列入國內公開出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------

clear
N1=256;
Nh=33;
q=1;
% q=1 表示臨界抽樣。

sig=fmlin(N1);
h=window(@gausswin,Nh);
h=h/norm(h);
figure(1)
[tfr,dgr,g]=tfrgabor(sig,16,q,h);


figure(2)
subplot(211)
plot(h);grid;
ylabel('h(t)')

subplot(212)
plot(g);grid;
ylabel('g(t)')