完全实施Pan Tompkins;
如果您发现此脚本有用，请引用以下引用;
％％ 参考 ：
％[1] Sedghamiz。H，“Pan Tompkins ECG QRS探测器的Matlab实施。”，2014年3月。https://www.researchgate.net/publication/313673153_matlab_implementation_of_pan_tompkins_ecg_qrs_detect.
和
%[2]。J,汤普金斯。王志强，“一种实时QRS检测算法”
《生物医学工程学报》，vol. bme-32, no . 33、1985年3月。
作者:胡曼·塞德哈密兹
%林雪平大学
邮箱:hoose792@student.liu.se
版权所有2014年3月
-----------------
％％ 方法 ：
% %预处理
％1）带通滤波器（5-15 Hz）
％2）终止滤波器到高光QRS复合物。
% 3)信号平方。
% 4)信号是噪声的平均值(0.150秒长度)。
％5）根据信号的采样频率滤波
％选项更改为最佳匹配ECG信号的特性
%%决定规则
在算法的这一点上，前面的步骤已经产生了一个大致的脉冲形状
在MWI输出的%波形。关于这个脉冲是否
％对应于QRS复合物（而不是高倾斜的T波或噪声假物）是
％以自适应阈值操作和其他决定执行
%规则概述如下;
% a)基准标记-首先对波形进行处理以产生一组加权单位
MWI Maxima位置处的％样本。这是完成的，以便本地化QRS
％复合物到单个时间的时间。W [k]加权是最大值。
% b)阈值-当分析MWI输出的幅值时，该算法使用
%两个阈值(THR_SIG和THR_NOISE，适当初始化在简短
% 2第二训练阶段)，不断适应变化的心电信号质量。的
％首先通过y [n]使用这些阈值来分类每个非零样本
% (CURRENTPEAK)作为信号或噪声:
％if currentpeak> thr_sig，该位置被标识为“qrs复杂”
% candidate”和信号电平(SIG_LEV)被更新:
% sig _ Lev = 0.125 ×currentpeak + 0.875× sig _ Lev
％如果thr_noise % " noise peak "和noise level (NOISE_LEV)被更新:
% noise _ Lev = 0.125×currentpeak + 0.875× noise _ Lev
%根据信号和噪声水平(SIG_LEV和NOISE_LEV)的新估计，
％分别）在ECG中的该点处，阈值调整如下：
％thr _ sig =噪声_ lev + 0.25×（sig _ lev？噪音_ lev）
% thr_ noise = 0.5× (thr_ sig)
％这些调整在被认为的信号段中逐渐降低阈值
质量较差。
％c）对错过的QRS复合物的搜查 - 在上面的阈值步骤中，如果
% CURRENTPEAK < THR_SIG，峰值被认为不是由QRS引起的
％ 复杂的。但是，如果没有禁止的禁止，那么不合理的长期已经过期
%峰值时，算法将假定QRS丢失并执行
％申请。这限制了假底片的数量。用于触发的最短时间
％a的验证是r峰值时间段的1.66倍（称为rr
%间隔)。这个值有一个生理来源-相邻时间值
％心跳不能比这更快地改变。假设错过的QRS复合体
％要在最高峰值的位置处于位于thr_sig和
％thr_noise。在该算法中，存储了两个平均RR间隔，第一个RR间隔是
%计算为最后八个QRS位置的平均值，以适应不断变化的心脏
%率和第二个RR区间平均值为平均值
最常常规的rr间隔的百分比。如果心率不规律，则阈值降低
％改善检测。
% d)在不应期内消除多次检测-是
如果在以前在200ms之后，则不可能发生合法的QRS复合物
%检测。这是一个生理上的限制-由于不应期期间
尽管刺激措施，患心室去极化的百分比可以[1]。作为QRS复合体
%候选，算法消除了这些物理上不可能的事件，
%从而减少误报。
％e）T波歧视 - 最后，如果QRS候选人发生在200ms之后
%不应期，但在之前QRS的360ms内，算法确定
％这是一个真正的QRS复合物的下一个心跳还是异常突出的
％t波。该决定基于该位置的波形的平均斜率。一个斜坡
%小于以前QRS复合体的一半则是一致的较慢
%改变T波的行为-否则，它成为QRS检测。
％额外的概念：在论文中提到的点旁边，这段代码也是如此
％检查是否发生少于360毫秒延迟的峰值也有一个
％延迟小于0.5 * mean_rr如果是，则这被视为噪声