基于matlab去除心电信号工频干扰的方法,小波变换消除心电信号中人为混入的50HZ工频干扰！...

fs=1000;                 %采样频率，采样定理2Fmax

N=3600;                  %总采样长度

t=0.001:1/fs:N/fs;       %时间的变化范围

data=load('rec_1_2.txt');

s=data(:,2);             %输入心电信号

A=0.2;                   %干扰信号的幅值

fai=pi/3;                %干扰信号的相移

noise=A*cos(2*pi*50*t+fai);  %50Hz工频干扰信号

x=s'+noise;

%N=length(t)

%ns=[zeros(1,100),0.2*randn(1,200),zeros(1,101)];%区间为101-201的白噪声

%ns=zeros(1,401);

sig=x;%叠加了噪声的信号

[c,l]=wavedec(sig,4,'db4');%对有噪信号用db4小波函数进行4层分解

a1=appcoef(c,l,'db4',1);  %提取第一层的近似分量

d1=detcoef(c,l,1);        %提取第一层的细节分量

a2=appcoef(c,l,'db4',2);  %提取第二层的近似分量

d2=detcoef(c,l,2);        %提取第二层的细节分量

a3=appcoef(c,l,'db4',3);  %提取第三层的近似分量

d3=detcoef(c,l,3);