matlab计算心率 QRS波群定位

clc
clear
%%
fs=500;%44.1e3;%采样率
samptime=8;%采样时间ECG调制用4秒，其他用0.5秒，便于观察
fm=150;%调制信号频率
am1=1;%调制信号幅度
fc=1e3;%载波频率
[z0t0]=ReadECGData（samptimefs）;%读取心电信号txt文件
Amplify=800;%电路放大倍数
[~z1t]=ReadMIT_BIH（‘117‘‘.\MIT-ECG data‘samptimefs）;
z=z1/1000*Amplify;%单位换算为V
figure（1）
plot（tz）
title（‘ECG原信号‘）;
data=z‘;
differ=zeros（40001）;
dif_max=0;
for i=3:1:3998
    differ（i1）=data（i+11）-data（i-11）+2*data（i+21）-2*data（i-21）;%求导数
end   
dif_max=max（differ）;
thresh=dif_max*0.375;          %阈值为导数最大值的0.375倍
NumOfR=0;
for i=1:4000
    if  （data（i1）>thresh）&&（differ（i1）*differ（i+1）<0） 
        NumOfR=NumOfR+1;
        R_position（NumOfR1）=i;
    end
end

for i=1:NumOfR-1
    DR_R（i1）=60*500/（R_position（i+11）-R_position（i1））;
end
figure（2）
 plot（（1:4000）‘dataR_positiondata（R_position）‘rs‘）;
 title（‘117号ECG信号R波定位‘）
 %%
%Q波定位
NumOfQ=0;
for i=1:NumOfR
    minQvalue=R_position（i）-40;
    for m=1:40
        if （data（minQvalue）>data（R_position（i）-40+m））
        minQvalue=R_position（i）-80+m;
       end
    end
    NumOfQ=NumOfQ+1;
    Q_position（NumOfQ1）=minQvalue;
end
figure（3）
plot（（1:4000）‘dataQ_posit