Pan Tompkins的完整实现
如果你觉得这个脚本有用，请引用下面的参考资料;
% %引用:
% [1] Sedghamiz。H，“Pan Tompkins心电QRS探测器的Matlab实现”，2014年3月。https://www.researchgate.net/publication/313673153_Matlab_Implementation_of_Pan_Tompkins_ECG_QRS_detect
和
％[2] Pan.j，Tompkins。W.J，“实时QRS检测算法”IEEE
生物医学工程百分比交易，卷。BME-32，没有。1985年3月3日。
%%作者：Hooman Sedghamiz
百分比大学
％电子邮件：hoose792@student.liu.se
％版权所有2014年3月
-----------------
% %方法:
%%预处理
% 1)带通滤波器(5-15 Hz)
% 2)引出滤波器，以高光QRS络合物。
％3）信号是平方的。
％4）信号平均噪声（0.150秒长）。
% 5)根据信号的采样频率进行滤波
%选项被改变，以最佳匹配您的心电图信号的特征
% %决策规则
％在算法中的这一点，前一阶段产生了大致脉冲形状的
MWI输出处的％波形。确定该脉冲的决定
%对应QRS复合体(相对于高倾斜t波或噪声伪影)是
%通过自适应阈值操作和其他决策执行
下面列出的％规则;
％a）基准标记 - 首先处理波形以产生一组加权单元
%样本在MWI最大值的位置。这样做是为了本地化QRS
%复杂到一个瞬间的时间。w[k]权重为最大值。
％b）阈值 - 当分析MWI输出的幅度时，算法使用
％两个阈值（thr_sig和thr_noise，在简短的情况下适当地初始化
％2秒训练阶段），不断适应改变ECG信号质量。这
%第一次通过y[n]使用这些阈值对每个非零样本进行分类
％（currentpeak）作为信号或噪音：
%如果CURRENTPEAK > THR_SIG，则该位置被标识为“QRS复合体”
％候选人“和信号级别（sig_lev）更新：
％SIG _ LEV = 0.125×CurrentPeak + 0.875×SIG _ LEV
%如果THR_NOISE < CURRENTPEAK < THR_SIG，则该位置被标识为
％“噪声峰值”和噪声级别（噪声_LEV）更新：
％噪声_ lev = 0.125×CurrentPeak + 0.875×噪声_ Lev
％基于对信号和噪声水平的新估计（Sig_Lev和keek_lev，
%)时，阈值调整如下:
% THR _ sig = noise _ Lev + 0.25 × (sig _ Lev ?噪音_ Lev)
％thr _噪声= 0.5×（thr _ sig）
这些调整在信号段中逐渐降低阈值
质量较差。
% c) SEARCHBACK FOR MISSED QRS complex -在上面的阈值步骤中，如果
％currentpeak %复杂。然而，如果一个不合理的长期限已经过期而没有超过阈值
％峰值，算法将假设已错过QRS并执行a
% searchback。这就限制了假阴性的数量。用来触发的最小时间
% a搜索返回是当前R峰值到R峰值时间段(称为RR)的1.66倍
％ 间隔）。该值具有生理来源 - 相邻之间的时间值
心跳的变化不能比这更快。假设缺失的QRS复合体
%发生在THR_SIG和之间的间隔中的最高峰值位置
% THR_NOISE。该算法存储两个平均RR间隔，第一个RR间隔为
%计算为最后八个QRS位置的平均值，以适应不断变化的心脏
％率和第二个rr间隔均值是平均值
最规则RR间隔的%。如果心率不规律，阈值就会降低
%提高检测能力。
％d）消除难敏时期内的多种检测 - 它是
%不可能发生一个合法的QRS复合体，如果它位于200毫秒后的前一个
％检测到一个。由于耐火期间，这一约束是一种生理学
%，尽管有刺激[1]也不能发生心室去极化。作为QRS波群
生成％候选物，该算法消除了这种物理上不可能的事件，
从而减少了误报。
% e) T波识别-最后，如果QRS候选出现在200ms之后
％耐火期，但在以前的QRS的360ms内，算法确定
是否为真正的下一个心跳QRS复合体或异常突出
% T波。这个决定是基于那个位置的波形的平均斜率。斜率为
％小于前一QRS复合物的一半与较慢的一致
T波的％改变行为 - 否则，它将成为QRS检测。
附加概念:除文中提到的要点外，此代码也适用
%检查出现的小于360毫秒的峰值延迟是否也有
% latency小于0,5*mean_RR如果是，这被计算为噪声