matlab心电信号滤波R波提取

hao=ecg;
fs=1000;
N=length（ecg）;
%低通滤波器
B=[1 0 0 0 0 0 -2 0 0 0 0 0 1];
A=[1 -2 1];
shi1=filter（BAhao）;
%高通滤波器
B1=[-1 zeros（115） 32 -32 zeros（114） 1];
A1=[32 -32];
shi2=filter（B1A1shi1）;
%50Hz陷波
B3=[1 -2*cos（2*pi*50/fs） 1];%系统函数的分子系数矢量
A3=[1 -2*0.8*cos（2*pi*50/fs） 0.64];%系统函数的分母系数矢量
shi3=filter（B3A3shi2）;
%300Hz带阻滤波器
B4=[1 -2*cos（2*pi*300/fs） 1];%系统函数的分子系数矢量
A4=[1 -2*0.8*cos（2*pi*300/fs） 0.64];%系统函数的分母系数矢量
shi4=filter（B4A4shi3）;
figure
subplot（511）;plot（hao）;title（‘原始信号‘）;
subplot（512）;plot（shi1）;title（‘低通滤波后‘）;
subplot（513）;plot（shi2）;title（‘高通滤波后‘）;
subplot（514）;plot（shi3）;title（‘50HZ陷波后‘）;
subplot（515）;plot（shi4）;title（‘300HZ带阻滤波后‘）;
%使用最小均方多项式平均（5点）
for i = 1:N-4
    zuixiaojunfang（i） = （（-3）*shi4（i）+12*shi4（i+1）+17*shi4（i+2）+12*shi4（i+3）+（-3）*shi4（i+4））/35;
end
%中值滤波
m=10;
sig_mfil=medfilt1（shi4m）;
sig_smooth=smooth（shi420）;
figure;
subplot（411）;plot（shi4）;title（‘滤波处理后的心电信号‘）;
subplot（412）;plot（zuixiaojunfang）;title（‘最小均方多