testpeakdetection.m

用于心电信号滤波

%
% Test program for R-peak detection of ECG signals.
%
% Dependencies: The baseline wander and ECG filtering toolboxes of the Open Source ECG Toolbox
%
% Open Source ECG Toolbox, version 1.0, November 2006
% Released under the GNU General Public License
% Copyright (C) 2006  Reza Sameni
% Sharif University of Technology, Tehran, Iran -- LIS-INPG, Grenoble, France
% [e-mail]

% This program is free software; you can redistribute it and/or modify it
% under the terms of the GNU General Public License as published by the
% Free Software Foundation; either version 2 of the License, or (at your
% option) any later version.
% This program is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
% Public License for more details. You should have received a copy of the
% GNU General Public License along with this program; if not, write to the
% Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
% MA  02110-1301, USA.

clc
clear all
close all;

% load('SampleECG1.mat');
filetext = '300';
PATH= 'E:\300'; % path, where data are saved
HEADERFILE = [filetext,'.hea'];     % header-file in text format
ATRFILE    = [filetext,'.atr'];         % attributes-file in binary format
DATAFILE   = [filetext,'.dat'];         % data-file
SAMPLES2READ=10000;         % number of samples to be read
[M sfreq]= RDDataFun(PATH,HEADERFILE,ATRFILE,DATAFILE,SAMPLES2READ);
ARRAY=M(:,1);
ARRAY = ARRAY(1:10000,5)';

fs = 10000;
t = [0:length(ARRAY)-1]/fs;

f = 1;                                      % approximate R-peak frequency

bsline = LPFilter(ARRAY,.7/fs);              % baseline wander removal

peaks = PeakDetection(ARRAY-bsline,f/fs);    % peak detection

figure;
plot(t,ARRAY','b');
hold on
plot(t,peaks,'r');
xlabel('time (sec.)');
legend('ECG','ECG Peaks');