Matlab计算语音信噪比

function snr = calcu_snr(clean_file,noisy_file)
%读取干净语音文件
[cl,~]=audioread(clean_file);  %读取干净语音信号
[ny,~]=audioread(noisy_file);  %读取带噪语音信号
len=size(ny,1);   %计算带噪语音信号长度
clean=cl(1:len);   %截取与带噪语音信号长度相同的干净语音信号
Ps=sum(sum((clean-mean(mean(clean))).^2));%干净语音信号功率
Pn=sum(sum((clean-ny).^2));           %噪声信号功率
snr=10*log10(Ps/Pn);
end