基于MATLAB的心电信号预处理

%-----------------高通滤波器抑制肌电干扰-------------------  
fp=35;
fs=50;  
ap=3;    %通带最大衰减 
as=20;                             %阻带最小衰减 
fs1=200;                          %抽样频率 
wp=2*pi*fp/fs1;       
ws=2*pi*fs/fs1;                     %数字角频率  
Fs=fs1/fs1;
T=1/Fs;  
Op=2/T*tan（wp/2）;         
Os=2/T*tan（ws/2）;                   %角频率畸变 
[N Wn]=buttord（OpOsapas‘s‘）;  
[z p k]=buttap（N）;                    %G（p）的零点极点增益 
[b a]=zp2tf（zpk）;                    %G（p）的分子、分数系数 
[B A]=lp2lp（baOp）;                  %H（s）的分子、分数系数 
[Bz Az]=bilinear（BAFs*T/2）;          %H（z）的分子、分数系数
[H w]=freqz（BzAz256Fs*fs1）;        %H 
plot（wabs（H）‘r‘）; 
title（‘低通滤波器‘）; 
grid on; 
data=xlsread（‘maibo1.xlsx‘）;
x=data（:1）;
n=0:1/200:3-1/200;
y=filter（BzAzx）;
figure
plot（ny（9001:9600））;


%-----------------带陷滤波器抑制工频干扰-------------------  
%50Hz陷波器：由一个低通滤波器加上一个高通滤波器组成  
%而高通滤波器由一个全通滤波器减去一个低通滤波器构成  
Me=100;               %滤波器阶数  
L=100;                %窗口长度  
beta=100;             %衰减系数  
Fs=200;  
wc1=49/Fs*pi;     %wc1为高通滤波器截止频率，对应51Hz  
wc2=51/Fs*pi     ;%wc2为低通滤波器截止频率，对应4