短时傅里叶变换，画出时间频率图LOFAR

load F:\frogman.wav  %载入数据（声音）  
signalOfChirp = chirpData‘;  
signalOfChirp_5 = signalOfChirp/5;  
lengthOfSignal = length（signalOfChirp）;  
Fs = lengthOfSignal;  %采样频率  
T = 1/Fs;  
t = （0:lengthOfSignal-1）*T;  
  
% 绘原始波形图  
figure;  
subplot（311）;  
plot（t signalOfChirp）;  
grid on  
title（‘Original signal（chirpData）‘）;  
xlabel（‘Time‘）;  
ylabel（‘Amp‘）;  
  
%% 计算信号的fft  
NFFT = 2^nextpow2（lengthOfSignal）; %求最接近2的幂的数作为FFT计算的点数  
fMax = NFFT/2 + 1;