exa110702_specgram.asv

来自「这是一个对时间序列进行短时傅立叶变换的源程序,测试过,好用的.」· ASV 代码 · 共 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,512,1,hanning(200),64); 
save S.dat S -ascii 
save F.dat F -ascii
save T.dat T -ascii 
subplot(212);
surf(T,F,abs(S).^2)
view(0,50);
shading facted;
colormap(red);
xlabel('Time')
ylabel('Frequency')
zlabel('spectrogram')