%《數字信號處理教程——MATLAB釋義與實現》第九章例9.5.3程序hc953
% 短時間傅立葉分析及時頻分析
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
t=0:0.001:2;                 % 1kHz 采樣頻率，持續2秒
x=chirp(t,0,1,150);          % 從直流開始，在1秒時達到150Hz的線性調頻 
subplot(2,1,1),plot(t,x)     % 畫出信號曲線  
subplot(2,1,2),specgram(x);  % 顯示頻譜圖
set(gcf,'color','w')         % 置圖形背景色為白
colormap('gray')