matlab应用于心电信号检测

%可自动检测RSRR间期;如果你仔细研究代码还会发现它能自动检测QRS
%[hrv R_t R_amp R_index S_t S_amp] = ECGpeak（一维心电数据A/D转换率）;   
%代码:  
data=load（‘Ecgdata.txt‘）;
samp_freq=0;
[hrv R_t R_amp R_index S_t S_amp] = ECGpeak（datasamp_freq） 

% [hrv R_t R_amp R_index S_t S_amp] = ECGpeak（data samp_freq）; 
% R_t == RR points in time;
% R_amp == amplitude of R peak in bpf data;    
% S_amp == amplitude of following minmum; 
thresh = 0.2; 
testmode = 1; 
nargin=1
%%%%%%%%%%% make sample frequency default 256 Hz 
if nargin < 2
  samp_freq = 256; 
  if（testmode==1） 
    fprintf（‘Assuming sampling frequency of %iHz\n‘samp_freq）;  
  end 
end
  
%%%%%%%%%%% check format of data %%%%%%%%%% 
[a b] = size（data）; 
if（a>b） 
len =a;   
end   
if（b>a） 
len =b; 
end 

%%%%%%%%%% if there‘s no time axis - make one 
if （a | b == 1）;
% make time axis  
tt = 1/samp_freq:1/samp_freq:ceil（len/samp_freq）; 
t = tt（1:len）; 
x = data;  
end 
%%%%%%%%%% check if data is in columns or rows 
if （a == 2）
x=data（:1）; 
t=data（:2）;
end  
if （b == 2）   
t=data（:1）; 
x=data（:2）; 
end  

%%%%%%%%% bandpass filter data - assume 256hz data %%%%%  y  
% remove mean  
x = x-mean（x）;  
 
% FIR filtering stage 
bpf=x; %Initialise  
if（ （samp_freq==128） & （exist（‘filterECG128Hz‘）~=0） ） 
    bpf = filterECG128Hz（x）;  
end 
if（ （samp_freq==256） & （exist（‘filterECG256Hz‘）~=0） ）   
    bpf = filterECG256Hz（x）;  
end  

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%% differentiate data %%%%%%%%%%%%%%%%%%%%%%%%%%% F 
dff = diff（bpf）; % now it‘s one datum shorter than before

%%%%%%%%%square data   %%%%%%%%%%%%%%%%%%%%%%%%%%% 
sqr = dff.*dff;   
len = len-1; % how long is the new vector now? 
%%%%%%%%% integrate data over window ‘d‘ %%%%%%%%%%%%%%%%%%%%%%%%%  
d=[1 1 1 1 1 1 1]; % window size - intialise   
if （samp_freq>=256） % adapt for higher sampling rates    
  d = [ones（1round（7*samp_freq/256））]; 
end    
% integrate 
mdfint = medfilt1（filter（d1sqr）10）;
% remove filter delay for scanning back through ECG Z“ ~.Y@]@T  
delay = ceil（length（d）/2）;   
mdfint = mdfint（delay:length（mdfint））;  
%%%%%%%%% segment search area %%%%%%%%%%%%%%%%%%%%%%%