基于matlab对ECG信号进行滤波处理

由单片机采集心率数据，截取部分数据处理如下 ：

前提：根据采样定理，心率信号采样频率为50hz  （大于心率的频率两倍），得带数据信号数据源“B.TXT”文件。

1. 原始数据的采集截取有效部分分析。

2.数据滤波（50Hz的工频干扰以及放大器本身的漂移），要根据频谱具体分析

3.FDA数字滤波器设计，相应系数的得到（本设计2.5HZ滤波高通滤波器(设计的不合理，只做实验测试) ）

4.滤波数据的显示

PS:虽然滤波后数据还是有一点不规律，但是已经相对原始数据好看很多，也还有很多不足之处。

clc;
clear all;

A=importdata('b.TXT') ;
L=length(A);
data=A(3900:4200);
figure
% subplot(311)
% plot(A(3900:4200));grid on;
subplot(312)
plot(data);grid on;
title('数据库数据')

FS=50;T=1/FS;               %50HZ  
L=length(data);  
t=(0:L-1)*T;  
  
NFFT = 2^nextpow2(L)  
Y = fft(A,NFFT)/L;  
f = FS/2*linspace(0,1,NFFT/2+1);  
subplot(313);  
plot(f,2*abs(Y(1:NFFT/2+1)));grid on; 


IIR_B = [
    0.7294407226391,   -2.188322167917,    2.188322167917,  -0.7294407226391];

IIR_A= [
                  1,   -2.374094743709,    1.929355669091,  -0.5320753683121
];
         w01=0;
         w02=0;
         w03=0;
         w04=0;
   y0=zeros(1,300);      
         for i=1:300
                w01=data(i)-IIR_A(2)*w02-IIR_A(3)*w03-IIR_A(4)*w04;
                y0(i)=IIR_B(1)*w01+IIR_B(2)*w02+IIR_B(3)*w03+IIR_B(4)*w04;
               w04=w03;
               w03=w02;
               w02=w01;
         end
   subplot(311)
   plot(y0);