exa7_1.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 26 行

% exa070201_specgram.m ,for example 7.2.1 
% to test specgram.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;
subplot(211);plot(t,y1);
ylabel(' Chirp signal y1')
axis([0 1.05 -1.1 1.1])

% 求两个Chirp 信号和的短时傅里叶变换；
subplot(212);
specgram(y,127,1,hanning(127),126); 
xlabel('时间(ms)')
ylabel('频率(kHz)')
title('STFT时频图')