专利名称：基于二阶导数编码的自动判别室速室颤的体外除颤器的制作方法
技术领域：
本发明属于医用设备技术领域，具体涉及一种自动体外除颤器。
背景技术：
心脏性猝死(S⑶)是指由于心脏原因引起的突然发生的自然死亡。造成心脏性猝死的原因大部分是在各类心血管病变基础上发生的一时性功能障碍和电生理改变，并引起恶性室性心律失常如室性心动过速(室速，VT)、心室纤颤(室颤，VF)等。恶性室性心律失常发作后，抢救成功率会随时间的流逝不断降低，在7-10分钟后，抢救成功率一般会低至 10%，因此迅速地电击除颤是终止多数快速性恶性室性心律失常、挽救患者生命的首选有效方法。自动体外除颤器(AED)是实现早期除颤的主要工具。由于恶性室性心律失常的突发性，医院外的患者存活率很可能依赖于对公众培训的肺复苏抢救(cardiopulmonary resuscitation, CPR)和组织完备的公众除颤(public access defibrillation, PAD) 工程。为了推广公众除颤，1997年美国心脏病协会(AHA)推出了公众场所用AED (AED for Public Access Defibrillation)指南。根据除颤的风险和可能的收益，该指南将心律分为以下三大类可电击复律心律(shockable rhythms, ShR)，不可电击复律心律 (nonshockable rhythms, NShR)禾口中间心律(Intermediate rhythms)。为协、助非专业人士正确使用AED，快速准确地自动检测出可电击复律心律(SiR)是AED实现其功能的关键环节。AHA建议对患有心室纤颤的成人进行初次能量为120J-200J双相波除颤或360J的单相波除颤，对稳定型单型性室性心动过速进行首剂量100J的单相波或双相波电复律，对复杂的多型性心动过速或无脉搏室性心动过速需类同心室纤颤处理。另一方面，植入式自动复律除颤器(ICD)在临床上已取得广泛应用。已有研究表明80%的除颤电击实际上是由室速误触发，这些误电击不仅对患者的心脏造成一定的损伤，给意识尚清楚的患者带来强烈的疼痛感，而且会影响ICD的使用寿命。因此，有必要研究出可靠的室速室颤判别算法。除颤能量过小很可能导致除颤失败；除颤能量过大则会大大增加对患者机体的热损伤，同时增加患者的疼痛感。除颤极板间的经胸阻抗决定了相同大小的电压下施加在患者上的能量。本发明即是在根据二阶导数编码的高波动性矢量概率判别出室速/室颤的基础上，决策合适的除颤方案，并由经胸阻抗调整合适的除颤电压。本发明不仅可提高除颤的针对性和成功率，还可降低对患者心肌的损伤率和疼痛感。
发明内容
本发明的目的在于提供一种不仅可提高除颤的针对性和成功率，还可降低对患者心肌的损伤率和疼痛感的自动体外除颤器。本发明提供的自动体外除颤器，是基于一种本发明提出的自动判别可电击复律心律为室速和室颤的方法的。
本发明提出的自动判别可电击复律心律是否为室速或室颤的方法，具体步骤如下
1.对采到的心电信号进行滤波预处理；具体步骤包括
(1)通过一个5阶的滑动平均滤波器滤除高频噪声；
(2)通过一个截止频率为IHz的高通滤波器抑制基线漂移；
(3)通过一个截止频率为30Hz的低通滤波器进一步滤除无关的高频成分；
2.识别预处理后的心电信号是否为可电击复律心律，根据可电击复律心律识别方法， 判别心电信号是否为可电击复律心律。(此步引用专利《一种改进的可电击复律识别仪器》 (专利号 ZL 2009 1 0045152. 2、ZL 2009 1 0045153. 7、ZL 2009 1 0045151. 8))
3.计算可电击复律心律的二阶导数编码的高波动性矢量概率；具体步骤包括
(1)对可电击复律心律进行降采样处理，以降低计算量；
(2)对长度为N的时间序列{u(i)，i=l，2，"《，N}，求其二阶导数绝对值和均值i&，
就是在判定一段心电图为可电击复律心律之后，再求这段心电信号的二阶导数绝对值和均值；N —般取600—1000，在一个实施例中取N=800 ；
⑶构建m维二阶导数绝对值矢量IO = Cs⑴，s ( +1), -",S (i+m-1) ]，i = 1，2，…，N-m+1，
数值m由实验验证确定，一般m为15 — 30，在一个实施例中m=22 ；
(4)将二阶导数绝对值矢量^TU转化为符号序列SX(i)
Of iJj·) < a* Ms 1，Ij(Z) > a* Ms
其中Z/i)和SX/i)分别表示iTi)和SX(i)的第j个元素，a是调整系数，a值由实验测得；a —般为0. 8—2. 5，在一个实施例中，a=l. 4 ；
(5)统计符号序列中高波动性矢量对应的符号序列的数目L；其中，高波动性矢量是指符号序列中“1”的比例大于1/3的序列对应的矢量；
(β)计算高波动性矢量概率IKK
PHVV =——-——； N - m - I
4.将由上述步骤计算得到的高波动性矢量概率与阈值进行比较，如果高波动性矢量概率小于阈值则为室速，如果大于阈值则为室颤。该阈值根据实际经验和试验确定，一般为0. 1—0. 5，在一个实施例中，取为0. 22 ；
最后根据判别出的室速或室颤，决策除颤方案对于室速患者，施行低能级同步电复律；对于室颤患者，施行高能级非同步电除颤。基于上述判别方法，本发明提出了相应的判别系统。该系统包括如下计算模块 Sl预处理模块对采集到的心电信号进行滤波；
S2识别模块用于识别心电信号是否为可电击复律心律；
S3高波动性矢量概率计算模块，用于计算可电击复律心律的二阶导数编码的高波动性矢量概率；
S4判别模块，通过高波动性矢量概率与阈值进行比较，判别可电击复律心律是室速还
是室颤。上述4个模块分别对应于自动判别可电击复律心律是否为室速或室颤的方法的4 个步骤的计算。本发明中，所述的预处理模块，包括如下子模块
Sll子模块，一个5阶的滑动平均滤波器，用于滤除高频噪声；
S12子模块，一个截止频率为IHz的高通滤波器，用于抑制基线漂移；
S13子模块，一个截止频率为30Hz的低通滤波器，用于进一步滤除无关的高频成分。本发明中，所述的可电击复律心律的识别模块，根据可电击复律心律识别方法，判别心电信号是否为可电击复律心律。本发明中，所述的二阶导数编码的高波动性矢量概率计算模块，包括如下子模块
S31子模块，用于将心电信号降采样至IOOHz ；
S32子模块，用于求心电信号序列的二阶导数绝对值和均值；
S33子模块，用于构建二阶导数绝对值矢量；
S34子模块，用于将绝对值矢量转化为符号序列；
S35子模块，用于统计符号序列中高波动性矢量对应的符号序列的数目；
S36子模块，用于计算高波动性矢量概率。这里的6个子模块对应于可电击复律心律的二阶导数编码的高波动性矢量概率计算的6个步骤。本发明还提出一种自动判别可电击复律心律是否为室速或室颤的自动体外除颤器，其系统框图如图1所示，包括心电采集及分析单元20、经胸阻抗测量单元30、充电/放电电路单元40、人机交互界面单元50、总控单元60和存储打印功能单元70。心电采集及分析单元20和经胸阻抗测量单元30将采集到的心电信号和经胸阻抗值发送给总控单元60 ； 总控单元60自动判别可电击复律心律为室速或室颤后，决策合适的除颤方案，向充电/放电电路单元40发送除颤指令；充电/放电电路单元40接收到除颤指令后对患者实施除颤， 同时将除颤器自身状况参数发送给总控单元60。在这种全自动工作模式之外，通过可视化的人机交互界面50，操作者还可对AEDlO的各项参数进行设置或直接手动操作AEDlO进行除颤。现场的有关数据存储和打印功能则通过存储打印功能单元70完成。上述自动体外除颤器中的总控单元60包括本发明所述的上述判别系统。本发明提供了一种识别准确、计算简单、能满足应用要求的基于二阶导数编码的高波动性矢量概率的心律识别算法，可提高判别可电击复律心律中室速和室颤的敏感性和特异性，简化算法的计算复杂度，缩短AED的判决时间。本发明还提供了基于该方法的判别系统和自动体外除颤器。
图1是本发明的系统框图。图2是本发明的工作流程。图3是本发明的预处理Sl步骤。图4是本发明的计算二阶导数编码的高波动性矢量概率S3的步骤。图5是本发明的判别室速室颤S4和决策除颤方案S5的步骤。
具体实施例方式下面结合附图对本发明作进一步的描述。图1是本发明的系统框图。本发明的主体为AED10，其硬件系统构成为心电采集及分析单元20、经胸阻抗测量单元30、充电/放电电路单元40、人机交互界面单元50、总控单元60和存储打印功能单元70。心电采集及分析单元20和经胸阻抗测量单元30将采集到的心电信号和经胸阻抗值发送给总控单元60 ；总控单元60自动判别可电击复律心律为室速或室颤后，决策合适的除颤方案，向充电/放电电路单元40发送除颤指令；充电/放电电路单元40接收到除颤指令后对患者实施除颤，同时将除颤器自身状况参数发送给总控单元60。在这种全自动工作模式之外，通过可视化的人机交互界面50，操作者还可对AEDlO 的各项参数进行设置或直接手动操作AEDlO进行除颤。现场的有关数据存储和打印功能则通过存储打印功能单元70完成。图2为本发明的工作流程。首先进行预处理(模块Si)，对采集到的心电进行滤波；然后根据电击复律心律识别方法识别预处理后的心电信号是否为可电击复律心律(模块S2)，如果是可电击复律心律则继续下面步骤(模块S3)，如果不是可电击复律心律则返回继续采集心电，此步引用专利《一种改进的可电击复律识别仪器》(专利号ZL 2009 1 0045152. 2, ZL 2009 1 0045153. 7, ZL 2009 1 0045151. 8)；计算可电击复律心律的二阶导数编码的高波动性矢量概率(模块S3)；判别室速/室颤(模块S4)将由上述步骤计算得到的高波动性矢量概率与阈值进行比较，如果高波动性矢量概率小于阈值则为室速，如果大于阈值则为室颤；最后根据判别出的室速或室颤，决策除颤方案S5 对于室速患者，施行低能级同步电复律；对于室颤患者，施行高能级非同步电除颤。图3是本发明的预处理(模块Si)的步骤。依次通过一个5阶的滑动平均滤波器用于滤除高频噪声(子模块S11)，一个截止频率为IHz的高通滤波器用于抑制基线漂移(子模块S12)，一个截止频率为30Hz的低通滤波器进一步滤除无关的高频成分(子模块S13)。图4是本发明的计算二阶导数编码的高波动性矢量概率(模块S3)的步骤。先进行降采样处理(子模块S31)，降采样率为100Hz，以降低计算量；然后由子模块S32对长度为N (本实施例中取800)的时间序列{u(i), i=l,2,-,N},求其二阶导数绝对值 s (V及其均值；再由子模块S33构建m维(本实施例中取22) 二阶导数绝对值矢量 X(i) = Cs⑴，s (i+1),…，s (i+ni-1) ]，i=l, 2，…，N-m+1 ；由子模块SM将二阶导数绝对值矢量KU转化为符号序列SX(i)
权利要求
1. 一种自动判别可电击复律心律是室速还是室颤的方法，其特征在于具体步骤为一、对采到的心电信号进行滤波预处理；二、识别预处理后的心电信号是否为可电击复律心律，即根据可电击复律心律识别方法，判别心电信号是否为可电击复律心律；三、计算可电击复律心律的二阶导数编码的高波动性矢量概率；具体步骤包括(1)对可电击复律心律进行降采样处理，以降低计算量；(2)对长度为N的时间序列|u(i)，i=l，2，…，N}，求其二阶导数绝对值s^和均值；⑶构建m维二阶导数绝对值矢量ZO Ds⑴，s (i+1),…，s (i+m-1)]， i = l，2，’"，N-I^l,数值m由实验验证确定；(4)将二阶导数绝对值矢量KU转化为符号序列SX(i)
2. 一种基于权利要求1所述方法的用于自动体外除颤器的自动判别可电击复律心律是室速还是室颤的系统，其特征在于，包括Sl预处理模块，对采集到的心电信号进行滤波； S2识别模块，用于识别心电信号是否为可电击复律心律；S3高波动性矢量概率计算模块，用于计算可电击复律心律的二阶导数编码的高波动性矢量概率；S4判别模块，通过高波动性矢量概率与阈值进行比较，判别可电击复律心律是室速还是室颤；上述4个模块分别对应于权利要求1中步骤一 步骤四； 其中，S3高波动性矢量概率计算模块，包括如下子模块PHVVS31子模块，用于将心电信号降采样至IOOHz ；S32子模块，用于求心电信号序列的二阶导数绝对值和均值；S33子模块，用于构建二阶导数绝对值矢量；S34子模块，用于将绝对值矢量转化为符号序列；S35子模块，用于统计符号序列中高波动性矢量对应的符号序列的数目；S36子模块，用于计算高波动性矢量概率；这里的6个子模块对应于权利要求1中可电击复律心律的二阶导数编码的高波动性矢量概率计算的6个步骤。
3.根据权利要求2所述的系统，其特征在于所述的预处理模块，包括包括如下子模块 Sll子模块，一个5阶的滑动平均滤波器，用于滤除高频噪声；S12子模块，一个截止频率为IHz的高通滤波器，用于抑制基线漂移；S13子模块，一个截止频率为30Hz的低通滤波器，用于进一步滤除无关的高频成分。
4.一种自动体外除颤器，其特征在于包括心电采集及分析单元(20)、经胸阻抗测量单元(30)、充电/放电电路单元(40)、人机交互界面单元(50)、总控单元(60)和存储打印功能单元(70)；所述心电采集及分析单元(20)和经胸阻抗测量单元(30)分别将采集到的心电信号和经胸阻抗值发送给总控单元(60);总控单元(60)自动判别可电击复律心律为室速或室颤后，决策合适的除颤方案，向充电/放电电路单元(40)发送除颤指令；充电/放电电路单元(40)接收到除颤指令后对患者实施除颤，同时将除颤器自身状况参数发送给总控单元(60)；此外，通过可视化的人机交互界面(50)，操作者对自动体外除颤器的各项参数进行设置或直接手动操作自动体外除颤器进行除颤；存储打印功能单元(70)用于存储和打印现场的有关数据；所述总控单元(60)包括权利要求3所述的系统。
全文摘要
本发明属于医用设备技术领域，具体为一种基于二阶导数编码的自动判别室速室颤的体外除颤器。本发明通过计算信号的二阶导数编码的高波动性矢量概率，自动判别可电击复律心律是室速或室颤，在此基础上决策合适的除颤方案，具体包括预处理对采集到的心电信号进行滤波；识别心电信号是否为可电击复律心律；计算可电击复律心律的二阶导数编码的高波动性矢量概率；根据高波动性矢量概率判别室速/室颤；根据室速/室颤决策除颤方案。本发明可降低对患者的身心损伤，提高除颤的针对性和成功率。
文档编号A61N1/39GK102380164SQ20111022102
公开日2012年3月21日 申请日期2011年8月3日 优先权日2011年8月3日
发明者刘聪, 方祖祥, 杨圣均, 邬小玫 申请人:复旦大学