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）]）;