检测ECG信号的R波

在matlab环境下，实现ECG信号R波的检测，方法是小波变换。

clc
clear
% read the ECG singal data and translate original binary data into decimal data 
% picture the ECG data
% parameter settings： 
 PATH= ‘E:\02--ECG\Database\PTBdb\patient001‘; % set data storage path 
 HEADERFILE= ‘s0010_re.hea‘;      % .hea format，header file 
 DATAFILE=‘s0010_re.dat‘;         % .dat format，ECG data 
% output parameter：M —a matrix with SAMPLES2READ rows and 15 columnseach column represents a channel of data. 

% read the headfile
signalh= fullfile（PATH HEADERFILE）;                  % through the function fullfile for the full path of header file 
fid1=fopen（signalh‘r‘）;                              % open the header file with attribute ‘r-‘ 
z= fgetl（fid1）;                                       % read the first line data of the header file，strings format
A= sscanf（z ‘%*s %d %d %d‘[13]）;                   % according to the format ‘%*s %d %d %d‘translate the data and store them into matrix A ‘*s‘means ignore the strings 
nosig=A（1）;                                           % signal channel number
sfreq=A（2）;                                           % data frequency 
allnum=A（3）;                                          % data number   
clear A;                                              % clear matrix A and prepare to read the next line data 
for k=1:nosig                                         % read the data information of each channel
    z= fgetl（fid1）; 
    A= sscanf（z ‘%*s %d %d %d %d %d‘[15]）; 
    dformat（k）= A（1）;                                 % signal format 16bit 
    gain（k）= A（2）;                                    % the number of integer in each channel
    bitres（k）= A（3）;                                  % sampling accuracy 
    zerovalue（k）= A（4）;                               % ECG zero corresponding integer value
end; 
fclose（fid1）; 
clear A; 
 
% load binary data   
% .dat format file-datafile.dat   
%  two Byte represent a integer m
if dformat~= 16
    error（‘this script does not apply binary formats different to 16.‘）;  
end; 
SAMPLES2READ=allnum/nosig;                            % samples need to load  
signald= fullfile（PATH DATAFILE）;                    
fid2=fopen（signald‘r‘）;                              % open the header file with attribute ‘r-‘ 
M= fread（fid2 [nosig SAMPLES2READ] ‘int16‘）‘;      % load the data into the first columnsecond column......

i=input（‘input the channel number（1=X=0:1/sfreq:（SAMPLES2READ-1）/sfreq;
B=（M（:i）-zerovalue（i））/gain（i）; 

% filter the original sianal 
% Butterworth filter 
% the frequency response curve within the passband maintains maxinum flat 
% there are no ups and downs and within stopband is gradually reduced to zero
%wp：Passband cut-off frequency（通带阻止频率），ws：Stopband cut-off frequency（阻带截止频率）
%rp：Passband maximum attenuation（通带最大衰减），rs：Stopband minimum attenuation（阻带最小衰减） 
wp=150;ws=220;rp=0.5;rs=40;Fs=sfreq;             
[NWn]=buttord（wp/（Fs/2）ws/（Fs/2）rprs）;

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

     文件       5052  2013-03-15 19:26  detect_R_Wavelet.m

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

                 5052                    1