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