www.pudn.com > spectrogram_analysis.rar > yuputu.m

 
%%%%ÓïÆ×ͼ 
%%%%tm=len_x/fs; 
%%%%axis([0 tm 0 fs/2]); 
 
clear; 
close all; 
clc; 
 
[x,fs,p]=wavread('0.wav'); 
len_x=length(x); 
x=filter([1-0.9375],1,x); 
 
len=150; 
lenshift=30; 
y=enframe(x,len,lenshift);   %% 
 
[m,n]=size(y); 
 
for i=1:m 
    fft_mag=abs(fft(y(i,:),n)); 
    yupu(i,:)=fft_mag.^2; 
end 
 
s=sum(yupu); 
ss=sum(s); 
ave=ss/(m*n); 
yupu2=yupu/ave; 
yupu3=yupu2(:,1:len/2); 
 
figure; 
imshow(yupu3',[0 1]); 
 
title('ÓïÆ×ͼ'); 
 
xlabel('t/s'); 
ylabel('f/Hz');