心电信号预处理去噪

针对心电信号中存在肌电干扰、工频干扰以及基线漂移，本代码可以有效地进行去噪。

%-----------------带陷滤波器抑制工频干扰-------------------  
%50Hz陷波器：由一个低通滤波器加上一个高通滤波器组成  
%而高通滤波器由一个全通滤波器减去一个低通滤波器构成  

Me=100;               %滤波器阶数  
L=100;                %窗口长度  
beta=100;             %衰减系数  
Fs=1500;  
wc1=49/Fs*pi;     %wc1为高通滤波器截止频率，对应51Hz  
wc2=51/Fs*pi     ;%wc2为低通滤波器截止频率，对应49Hz  
h=ideal_lp（0.132*piMe）-ideal_lp（wc1Me）+ideal_lp（wc2Me）; %h为陷波器                                                         
w=kaiser（Lbeta）;  
y=h.*rot90（w）;         %y为50Hz陷波器冲击响应序列  
m2=filter（y1m）;  

figure  
subplot（211）;plot（abs（h））;axis（[0 100 0 0.2]）;  
xlabel（‘频率（Hz）‘）;ylabel（‘幅度（mv）‘）;title（‘陷波器幅度谱‘）;grid;  
N=512;  
P=10*log10（abs（fft（y）.^2）/N）;  
f=（0:length（P）-1）;  
subplot（212）;plot（fP）;  
xlabel（‘频率（Hz）‘）;ylabel（‘功率（dB）‘）;title（‘陷波器功率谱‘）;grid;  
     
figure  
subplot （211）; plot（TIMEm）;  
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘原始信号‘）;grid;  
subplot（212）;plot（TIMEm2）;  
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘带阻滤波后信号‘）;grid;  
    
figure  
N=512 ; 
subplot（211）;plot（abs（fft（m））*2/N）;axis（[0 100 0 1]）;  
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘原始信号频谱‘）;grid;  
subplot（212）;plot（abs（fft（m2））*2/N）;axis（[0 100 0 1]）;  
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘带阻滤波后信号频谱‘）;grid;    
  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2451  2018-08-14 21:34  心电信号预处理\100数据以及仿真结果\data.mat

     文件      10917  2018-08-14 21:09  心电信号预处理\100数据以及仿真结果\IIR零相移数字滤波器纠正基线漂移.png

     文件      10669  2018-08-14 21:55  心电信号预处理\100数据以及仿真结果\低通滤波器滤除肌电信号1.png

     文件       9897  2018-08-14 21:56  心电信号预处理\100数据以及仿真结果\带陷滤波器抑制工频干扰.png

     文件      10552  2018-08-14 22:02  心电信号预处理\100数据以及仿真结果\心电信号预处理.png

     文件       1408  2018-08-14 21:56  心电信号预处理\Filter.m

     文件       1027  2018-09-10 19:32  心电信号预处理\Filter1.m

     文件        218  2018-08-14 20:19  心电信号预处理\ideal_lp.m

     文件       1018  2018-08-14 20:35  心电信号预处理\IIR零相移数字滤波器纠正基线漂移.txt

     文件      91660  2018-08-14 22:03  心电信号预处理\result.mat

     文件       1927  2018-09-10 19:29  心电信号预处理\Test.m

     文件       1944  2018-08-14 15:52  心电信号预处理\低通滤波器滤除肌电信号（Test）.txt

     文件       1695  2018-08-14 20:14  心电信号预处理\带陷滤波器抑制工频干扰（Filter+ideal_lp）.txt

     目录          0  2018-08-22 19:24  心电信号预处理\100数据以及仿真结果

     目录          0  2018-08-22 19:24  心电信号预处理

----------- ---------  ---------- -----  ----

               145383                    15