潘汤普金斯的完整实现;
如果你发现这个脚本有用请引用以下引用;
% %引用:
% [1]Sedghamiz。H,”潘汤普金斯心电图QRS探测器的Matlab实现。”,2014年3月。https://www.researchgate.net/publication/313673153_Matlab_Implementation_of_Pan_Tompkins_ECG_QRS_detect
和
% [2]。J,汤普金斯。W。J, IEEE“实时QRS检测算法”
%生物医学工程事务,BME-32卷,没有。3、1985年3月。
% %作者:Hooman Sedghamiz
%林雪平大学
%电子邮件:hoose792@student.liu.se
%版权2014年3月
- - - - - - - - - - - - - - - - - -
% %方法:
% %预处理
% 1)带通滤波器(5 - 15赫兹)
% 2)由此衍生过滤器高亮度QRS波群。
% 3)信号的平方。
% 4)噪声信号平均长度(0.150秒)。
% 5)根据信号的采样频率过滤
%选项更改为最佳匹配你的ecg信号的特点
% %决策规则
%算法在这一点上,前阶段产生了约脉冲形状
%在MWI的输出波形。确定是否这个脉冲
%对应于一个QRS波群(而不是一个high-sloped让或噪音产物)
%执行与一个自适应的阈值操作和其他决定
%规则下面;
%)框标-波形是首先产生一组加权处理单元
%样本MWI最大值的位置。这样做是为了定位QRS
%复杂到一个即时的时间。w [k]权重最大值的值。
% b)的阈值,当MWI输出振幅的分析,该算法使用
%两个阈值(THR_SIG THR_NOISE,适当地在一个简短的初始化
% 2第二个训练阶段),不断适应变化的ECG信号质量。的
%首先通过y [n]使用这些阈值对每个非零的样本进行分类
% (CURRENTPEAK)作为信号或噪声:
%如果CURRENTPEAK > THR_SIG,位置被确定为一个“QRS波群
%的候选人”和信号电平(SIG_LEV)更新:
%团体_列弗= 0.125×CURRENTPEAK + 0.875×团体_列弗
%如果THR_NOISE < CURRENTPEAK < THR_SIG,那位置被确定为一个
%”噪声峰值”和噪音水平(NOISE_LEV)更新:
%噪音_列弗= 0.125×CURRENTPEAK + 0.875×噪音_列弗
%根据新估计的信号和噪音水平(SIG_LEV NOISE_LEV,
心电图分别%)在这一点上,调整阈值如下:
%刺_团体=噪音_列弗+ 0.25×(SIG _列弗吗?噪音_列弗)
%刺_噪音= 0.5×(刺_团体)
%这些调整阈值逐渐下降视为信号段
%是质量差。
% c) SEARCHBACK错过QRS复合物,在上面的阈值的步骤中,如果
% CURRENTPEAK < THR_SIG,峰值QRS被认为是没有了
%复杂。但是,如果不合理的长期没有abovethreshold已经过期
%的高峰,该算法假设QRS已经错过了和执行
% searchback。这限制了数量的假阴性。用于触发的最小时间
% searchback是目前的1.66倍R峰的峰值时间(称为RR
%间隔)。这个值有一个生理起源——相邻之间的时间价值
%心跳不能改变比这更迅速。错过的QRS波群
%在最高峰的位置出现在位于THR_SIG之间的时间间隔
% THR_NOISE。在该算法中,两个平均RR间隔存储,第一个RR间隔
%计算平均过去八QRS位置,以适应不断变化的心
%率和第二RR间隔的意思是意思
%最常规的RR间隔。如果阈值降低心率不正常
%,以提高检测。
% d)消除多个检测在不应期
合法的QRS波群%不可能发生如果谎言后200毫秒内
%检测。这个约束是一个生理期间由于不应期
%心室去极化不能发生尽管刺激[1]。作为QRS波群
%候选人产生,该算法消除了这些物理上不可能的事件,
%从而减少假阳性。
% e) T波歧视——最后,如果QRS候选人发生后200毫秒
%不应期在360之前的QRS女士,但是该算法决定
%这是否是一个真正的QRS波群的下一个心跳或异常突出
% T波。这个决定是基于波形的平均斜率在那个位置。斜率为
%的一半不到一个前面的QRS波群和慢是一致的
% T波改变行为习惯,否则,它就变成了QRS检测。
%额外的概念:在论文中提到的点,这段代码
%检查不到360毫秒延迟的发生高峰也一个
%延迟小于0,5 * mean_RR如果是的这是算作噪音