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（‘最小均方多