心电信号的预处理滤波器程序

心电信号的预处理滤波器的设计（包括数据文件和程序） 包括了FIR滤波器和IIR滤波器以及整系数滤波器，可以有效滤除基线漂移与工频干扰

%%%%%%%%%%  数据读取
ECG=load（‘1record.dat‘）;
l=length（ECG）;
t=[150/l:150/l:30]
x=ECG（4:）;    %可取不同的组
y=x（1:5:end）
figure（1）
plot（t 20*log（abs（y）））;
title（‘含噪心电信号‘）;
xlabel（‘时间（s）‘）;
ylabel（‘幅度（dB）‘）;
axis（[0 15 179 183]）

%%%%%%%%%%  频域分析
n=4000
m=abs（fft（yn））;
fs=200;
f=fs/n*（0:n-1）;
figure（2）
plot（f m）;
title（‘心电信号的频谱图‘）;
xlabel（‘频率 f/Hz‘）;ylabel（‘幅值/db‘）;
axis（[0 100 -100 40000]）


%%%%%%%%%%  加入工频噪声  时域及频域分析
x=20*sin（2*50*pi*t）;
y1=y+x;
figure（3）;
subplot（211）;
plot（t 20*log（abs（y1）））;
axis（[0 15 179 183]）
title（‘加工频噪声时域波形‘）;
xlabel（‘时间t/s‘）;ylabel（‘幅值/db‘）;
k=abs（fft（y1n））;	
subplot（212）;
plot（fk）;
axis（[0 100 0 40000]）
title（‘加工频噪声频谱图‘）;
xlabel（‘频率 /Hz‘）;ylabel（‘幅值‘）

%%%%%%%%%%  工频噪声的滤除   巴特沃思
fs=200;
wp=[47 53]*2/fs;
ws=[42 58]*2/fs;
Rp=3;Rs=45;
[Nwn]=buttord（wpwsRpRs）           %  N=6
[ba]=butter（Nwn‘stop‘）
figure（4）
[HW]=freqz（ba1024）;%生成滤波器的幅频响应
k=0:1023;
plot（（fs/2）/1024 *kabs（H））;%输出滤波器的频率响应
axis（[0 100 0 1

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

     文件       7814  2013-05-30 20:14  ECG.m

     文件    1800012  2013-04-19 11:02  1record.dat

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

              1807826                    2