6-3.m

来自「Matlab数字信号处理与应用源代码」· M 代码 · 共 20 行

%例程6－3  计算非平稳随机信号的STFT时频谱
% e.g.6-3.m for example6-3;

clear all;
t=0:0.001:1.024-0.001;
N=1024;
% 得到两个Chirp信号，并相加；
y1=chirp (t, 0, 1,350);
y2=chirp (t, 350, 1, 0);
y=y1+y2;
subplot (211);plot(t,y);
title ('Chirp Signal y')
%求两个Chirp信号和的短时傅立叶变换；
ym=y';
tfrgabor(ym);
subplot(212)
xlabel('Time[ms]');
ylabel('Frequency[kHz]')
title('GABOR时频谱,Lh=16,Nf=64,N=1024,Q=4')