test_0db.m

用MATLAB实现了该算法

% 画原始音频和录制音频频谱图
[y,Fs,bits]=wavread('record_35db.wav'); % 录制音频
[y1,Fs1,bits1]=wavread('h_35b.wav');  %原始音频  

y=y(:,1); 
y1=y1(:,1);

sigLength=length(y);
Y=fft(y,sigLength);   
Pyy=Y.*conj(Y)/sigLength; 
halflength=floor(sigLength/2); 
f=Fs*(0:halflength)/sigLength; 

subplot(2,1,1);
plot(f,Pyy(1:halflength+1));
title('原始音频功率谱');

sigLength1=length(y1);
Y1=fft(y1,sigLength1);   
Pyy1=Y1.*conj(Y1)/sigLength1; 
halflength1=floor(sigLength1/2); 
f1=Fs1*(0:halflength1)/sigLength1; 

subplot(2,1,2);
plot(f1,Pyy(1:halflength1+1));
title('录制音频功率谱');