信噪比matlab代码

function snr1=SNR1（ys3）
% 去躁后的信噪比
% I: original signal
% In: noisy signal（ie. original signal + noise signal）
% snr=10*log（sigma2（I2）/sigma2（I2-I1））
Ps=sum（sum（y.^2））;%signal power
Pn=sum（sum（（s3-y）.^2））;%noise power
snr1=10*log10（Ps/Pn）;
disp（[‘SNR1 =‘num2str（snr1）]）;