exa110702_specgram.m

来自「这是一个对时间序列进行短时傅立叶变换的源程序,测试过,好用的.」· M 代码 · 共 26 行

%-------------------------------------------------------------------------
% exa110702_specgram.m ,for example 11.7.2  and fig11.7.1
% to test specgram.m 
%-------------------------------------------------------------------------
clear;
load ttt1;
subplot(211);plot(ttt1);
ylabel('GR')
Xlabel('DEP')

% 求信号的短时傅里叶变换；
[S,F,T]=specgram(ttt1,256,1,hanning(200),100); 
save S.dat S -ascii 
save F.dat F -ascii
save T.dat T -ascii 
subplot(212);
surf(T,F,abs(S).^2)
view(-1,1);
shading facted;
colormap(red);
xlabel('Time')
ylabel('Frequency')
zlabel('spectrogram')