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

 
%%%%ÓïÆ×ͼ 
 
clear; 
close all; 
clc; 
 
[x,fs,p]=wavread('1y.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]); 
 
tm=len_x/fs; 
axis([0 tm]);