test0807gabor1.m

来自「信号处理实例」· M 代码 · 共 29 行

%---------------------------------------------------------------------------------------
% test0807gabor1.m,  for  example 08.07
% chirp signal and its gabor transform
% 2006.12.19 LYB
%----------------------------------------------------------------------------------------

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);
tfrgabor(sig,32,q,h);

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

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