ECG心率计算

fs=200; 
sigtemp=y;  
siglen=length（y）; 
sigmax=[];  
yabs=abs（y）;
rate=1/fs;
for i=1:siglen-2     
    if  （y（i+1）>y（i）&y（i+1）>y（i+2））|（y（i+1）        sigmax=[sigmax;abs（sigtemp（i+1））i+1];     
    end; 
end; 

%取阈值阈值为相对幅值的差的60%   
%最大幅度平均值，8个最大幅值点的平均值    
thrtemp=sort（sigmax）; 
thrlen=length（sigmax）; 
thr=0;  
for i=（thrlen-7）:thrlen     
    thr=thr+thrtemp（i）; 
end;  
thrmax=thr/8;               
 %最小幅度平均值，对消幅度，100个最小幅值点的平均值  
zerotemp=sort（y）; 
zerovalue=0; 
for i=1:100      
    zerovalue=zerovalue+zerotemp（i）; 
end;  
zerovalue=zerovalue/100;   
  
thr=（thrmax-zerovalue）*0.3; %最大、最小幅度的差值的30%为判别R波的阈值  

%定位R波  
rvalue=[];  
for i=1:thrlen      
    if sigmax（i1）>thr         
        rvalue=[rvalue; sigmax（i2）];     
    end;