clc;clear all;

clc;clear all;
close all;
xlsread D:\work\denoise\lead\TEST0010.xls;
so1=ans;
n=4096/2; %计算点数
dt=1/(1e6); %采样频率
s2=so1(1:n);
%-------add white guassian noise-------
%yy = awgn(s2,-2,'measured'); % Add white Gaussian noise.
xlsread D:\work\denoise\onsite\6_150filter.xls;
yo=ans(1:n);
%yy = awgn(s(1:n),-30,'measured'); % Add white Gaussian noise.
yy=yo+s2;
CC0=corrcoef(s2(1:1024),yy(1:1024));
Es0=0;
En0=0;
for i=1:1024
Es0=Es0+s2(i)^2;
En0=En0+(s2(i)-yy(i))^2;
end
SNR0=10*log10(Es0/En0);
noise=yy-s2;