exa7_3.m

matlab应用实例。如：FIR、IIR数字滤波器设计等。

% exa070301_tfrgabor.m ,for example 7.3.1 
% to test tfrgabor.m 

clear;
t=0:0.001:1.024-.001; 
N=1024; 

% 得到两个Chirp 信号，并相加；
y1=chirp(t,0,1,350);
y2=chirp(t,350,1,0);
y=y1+y2;

% 求两个Chirp 信号和的短时傅里叶变换；
ym=y';
tfrgabor(ym);

%tfr=tfrgabor(ym);
%imagesc(abs(tfr(1:64,:)))
%axis('xy')
%mesh(abs(tfr))
%contour(abs(tfr));