snr.m

用MATLAB实现了该算法

%求经伸缩归一化处理后的录制音频信噪比
F=wavread('h_35b');
F=F(:,1);
F1=wavread('record_35db');
F1=F1(:,1);  %第一声道
%归一化处理
N=length(F);
N1=length(F1);
if N>N1
    min=N1;
else
    min=N;
end
sum=0;
sum1=0;
for i=1:min
    sum=sum+abs( F(i) );
    sum1=sum1+abs( F1(i) );
end
fenmu=sum1/sum;
F2=F1/fenmu;%%%
F2_renew=renew(F2,N,N1); 
%线性伸缩处理
F=F';
%test=F-F2_renew;
%a=length(test)

a=norm(F-F2_renew,1).^2;    %根据公式求信噪比
b=norm(F,1).^2;
SNR=-10*log10(a/b)