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

clear
N1=256;
Nh=33;
q=4;
%q=4，表示現在是過抽樣；

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


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

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