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

clear
N1=256;
Nh=33;
q=4;
%q=4，表示現(xiàn)在是過(guò)抽樣；

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)')