%加噪声snr=10;sp=mean(yy.^2);attn=sp./ 10^(snr/10);attn = sqrt(attn);noise=randn(1,length(yy)).*attn;np=mean(noise.^2);snr1=10*log10(sp/np);data=yy+noise;