CN103007431B - 一种可电击信号的检测方法和装置 - Google Patents

一种可电击信号的检测方法和装置 Download PDF

Info

Publication number
CN103007431B
CN103007431B CN201210593464.9A CN201210593464A CN103007431B CN 103007431 B CN103007431 B CN 103007431B CN 201210593464 A CN201210593464 A CN 201210593464A CN 103007431 B CN103007431 B CN 103007431B
Authority
CN
China
Prior art keywords
signal
window
electric shock
parameter
detected
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201210593464.9A
Other languages
English (en)
Other versions
CN103007431A (zh
Inventor
李健
强媛
李德东
洪洁新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SHENZHEN BIOCARE BIO-MEDICAL EQUIPMENT Co Ltd
Original Assignee
SHENZHEN BIOCARE BIO-MEDICAL EQUIPMENT Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by SHENZHEN BIOCARE BIO-MEDICAL EQUIPMENT Co Ltd filed Critical SHENZHEN BIOCARE BIO-MEDICAL EQUIPMENT Co Ltd
Priority to CN201210593464.9A priority Critical patent/CN103007431B/zh
Publication of CN103007431A publication Critical patent/CN103007431A/zh
Application granted granted Critical
Publication of CN103007431B publication Critical patent/CN103007431B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种可电击信号的检测方法及装置,该方法包括:1)获得待检测心电信号;2)将时间长度为m的窗滑动截取所述待检测心电信号得到若干第一窗信号,对每个所述第一窗信号分别进行如下处理:2.1)对所述第一窗信号进行延时重构获得第一重构点,获取所述第一重构点的离散度参数Di;2.2)对所述第一窗信号进行归一化处理得到二进制数据序列,获取所述二进制数据序列中上升沿和下降沿的变化度参数Ni;2.3)以所述第一窗信号中相邻的波峰间或者相邻的波谷间的时间差值进行重构获得第二重构点,获取所述第二重构点的偏离度参数Si;3)根据所有第一窗信号的离散度参数Di、变化度参数Ni及偏离度参数Si,判断所述待检测心电信号是否为可电击信号。

Description

一种可电击信号的检测方法和装置
技术领域
本发明涉及医疗电子技术领域,具体涉及一种可电击信号的检测方法和装置。
背景技术
自从1970年后期由Dyack和Wellborn发明第一台AED以来,可电击信号的检测方法也开始研究。但是由于没有得到足够的重视,并没有达到理想的效果。今年来由于心室纤颤出现越来越多,但是通过抢救存活的却没有提高,人们开始逐渐重视可电击信号的检测方法,以便满足心室纤颤除颤需求。
通过对目前出现的检测方法的研究后发现,对于可电击信号的检测方法主要分为时域分析方法和频域分析方法等分析方法。
上述方法对部分心率偏小、部分心室扑动的心室纤颤漏检率高。
发明内容
本发明为了解决现有技术中可电击信号的检测方法对部分心率偏小、部分心室扑动的心室纤颤漏检率高问题,提供一种可电击信号的检测方法和装置。
本发明提供的技术方案如下:
本发明的一个目的是提供一种可电击信号的检测方法,包括以下步骤:
1)获得待检测心电信号y(n);
2)将时间长度为m的窗滑动截取所述待检测心电信号y(n)得到若干第一窗信号yi(m),对每个所述第一窗信号yi(m)分别进行如下处理:
2.1)对所述第一窗信号yi(m)进行延时重构获得第一重构点,获取所述第一重构点的离散度参数Di
2.2)对所述第一窗信号yi(m)进行归一化处理得到二进制数据序列,获取所述二进制数据序列中上升沿和下降沿的变化度参数Ni
2.3)以所述第一窗信号yi(m)中相邻的波峰间或者相邻的波谷间的时间差值进行重构获得第二重构点,获取所述第二重构点的偏离度参数Si
3)根据所有第一窗信号yi(m)的离散度参数Di、变化度参数Ni及偏离度参数Si,判断所述待检测心电信号是否为可电击信号。
作为一种优选,在所述步骤1)和2)之间还包括幅值判断步骤:将所述待检测心电信号y(n)幅值A与预设的第二幅度阈值A1进行比较,如果所述待检测心电信号y(n)幅值A小于第二幅度阈值A1,所述待检测心电信号y(n)归为不可可电击信号,结束;否则,进入下一步。
作为一种优选,所述待检测心电信号y(n)幅值A具体计算步骤如下:
3.1)将时间长度为T0的窗滑动截取所述待检测心电信号y(n)得到若干第二窗信号;
3.2)得到每个所述第二窗信号内的幅度最大值和最小值的差值;
3.3)所述待检测心电信号y(n)幅值A为所有所述第二窗信号幅度最大值和最小值的差值的平均值。
作为一种优选,所述步骤2.1)中的离散度参数Di具体计算如下:
4.1)将所述第一窗信号yi(m)重新采样获得重采样信号zi(m1);
4.2)对所述重采样信号zi(m1)采用下面的方法重构:以zi(m1)作为x轴坐标,以zi(m1+τ)作为y轴纵坐标,绘出第一重构点,其中,0<m1<m-τ;
4.3)用一个40×40的网格覆盖所有所述第一重构点,被所述第一重构点占据的格子数与所有格子数的比值即为离散度参数Di
作为一种优选,所述τ为0.8~0.95。
作为一种优选,所述步骤2.2)中变化度参数Ni具体计算如下:
6.1)将所述第一窗信号yi(m)分为时间长度为m2的第一滑动窗,得到第一滑动窗信号hii(m2);
6.2)对第一滑动窗信号hii(m2)进行归一化得到第一滑动窗归一化信号hi1(m2);
6.3)通过将所述第一滑动窗归一化信号hi1(m2)与给定第一幅度阈值V比较,得到每个所述第一滑动窗归一化信号hi1(m2)的二进制数据序列;
6.4)计算所有所述第一滑动窗归一化信号hi1(m2)的二进制数据序列中上升沿和下降沿的个数,进一步计算单个第一滑动窗归一化信号hi1(m2)的二进制数据序列中上升沿和下降沿的个数平均值即为变化度参数Ni
作为一种优选,所述第一幅度阈值V为0.15~0.25。
作为一种优选,如果所述第一滑动窗的归一化信号hi1(m2)的二进制数据序列中的上升沿或下降沿与相邻的上一个上升沿或下降沿间隔小于n1,则删除所述上升沿或下降沿,并且计一个下降沿后间隔n2点后再进行计算。
作为一种优选,所述步骤2.3)中偏离度参数Si具体计算如下:
9.1)将所述第一窗信号yi(m)每t2分一个小窗信号;
9.2)求得所述每个小窗信号波峰以及波峰所在的位置和波谷以及波谷所在的位置;
9.3)将所有所述小窗信号波峰和波谷按时间顺序分别连续排列,如果所有所述小窗信号波峰的均值大于波谷的均值,则对波峰对应的时间间隔进行重构获得第二重构点,如果所有所述小窗信号波峰的均值小于波谷的均值,则对波谷对应的时间间隔进行重构获得第二重构点;第二重构点的x轴横坐标为ti-ti-1(i=2,Λ,n),第二重构点的y轴纵坐标为ti+1-ti(i=2,Λ,n),ti代表连续排列中第i个波峰或波谷的位置;
9.4)计算每个所述第二重构点与原点连线之直线和直线y=x的直线夹角,其所有直线夹角相加总和即为偏离度参数Si
作为一种优选,所述t2为1.5~2.5秒。
作为一种优选,所述步骤3)中根据所有所述第一窗信号yi(m)的离散度参数Di、变化度参数Ni及偏离度参数Si,判断所述心电信号是否为可电击信号具体为:
对于每个第一窗信号yi(m),
如果P1≤Di≤P2,T1≤Ni≤T2,则
如果P1≤Di≤P2,T2<Ni≤T3,R1<Si,则
如果P1≤Di≤P2,T2<Ni≤T3,R1>Si,则
如果P1≤Di≤P2,T3<Ni≤T4,则
如果P1≤Di≤P2,T4<Ni<T5,R2<Si,则
如果P1≤Di≤P2,T4<Ni<T5,但不满足R2<Si,则
如果P1≤Di≤P2,Ni≥T5,则
如果P1≤Di≤P2,但不满足T1≤Ni≤T2和T2<Ni≤T3和T3<Ni≤T4和T4<Ni<T5和Ni≥T5,则
如果P2<Di,则
如果Di<P1,则
进一步计算所有所述第一窗信号yi(m)的的和sr1;
如果sr1大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号,结束;如果sr1小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号,结束。
作为一种优选,如果所述待检测心电信号为不可电击信号,则根据所有窗信号的离散度参数Di、变化度参数Ni及偏离度参数Si,进一步判断所述心电信号是否为可电击信号具体为:
对于每个第一窗信号yi(m),
如果R3≤Si≤R4,T6<Ni,则
如果R3≤Si≤R4,但不满足T6<Ni,则
如果Si<R4,但不满足R3≤Si≤R4,则
如果不满足R3≤Si≤R4和Si<R4,则
进一步计算所有第一窗信号yi(m)的的和sr2;
如果sr2大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号,结束;如果sr2小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号,结束。
作为一种优选,如果所述待检测心电信号再一次判断为不可电击信号,则根据所有窗信号的离散度参数Di、变化度参数Ni及偏离度参数Si,进一步判断所述待检测心电信号是否为可电击信号具体为:
对于每个第一窗信号yi(m),
如果T7<Ni≤T8,T8<Ni,则
如果T7<Ni≤T8,但不满足T8<Ni,则
如果R5<Si,但不满足T7<Ni≤T8,则
如果不满足T7<Ni≤T8和R5<Si,则
进一步计算所有第一窗信号yi(m)的的和sr3;
如果sr3大于等于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号;如果sr3小于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号,结束。
本发明实施例还公开了一种可电击信号的检测装置,包括数据获取模块、窗截取模块、离散度参数计算模块、变化度参数计算模块、偏离度参数计算模块和可电击信号判断模块;
所述数据获取模块用于获得待检测心电信号y(n);
窗截取模块用于将时间长度为m的窗滑动截取所述待检测心电信号y(n)得到若干第一窗信号yi(m);
离散度参数计算模块用于对所述第一窗信号yi(m)进行延时重构获得第一重构点,获取第一重构点的离散度参数Di
变化度参数计算模块用于对所述第一窗信号yi(m)进行归一化处理得到二进制数据序列,获取二进制数据序列中上升沿和下降沿的变化度参数Ni
偏离度参数计算模块用于以所述第一窗信号yi(m)中相邻的波峰间或者相邻的波谷间的时间差值进行重构获得第二重构点,获取第二重构点的偏离度参数Si
可电击信号判断模块用于根据所有第一窗信号yi(m)的离散度参数Di、变化度参数Ni及偏离度参数Si,判断所述待检测心电信号是否为可电击信号。
作为一种优选,在数据获取模块、窗截取模块之间还包括幅值判断模块,用于将所述预待检测心电信号y(n)幅值A与预设的第二幅度阈值A1进行比较,如果所述待检测心电信号y(n)幅值A小于第二幅度阈值A1,所述待检测心电信号y(n)归为不可可电击信号,结束;否则,所述待检测心电信号y(n)进入窗截取模块进行下一步判断。
作为一种优选,所述离散度参数计算模块具体为:将所述第一窗信号yi(m)重新采样获得重采样信号zi(m1);对所述重采样信号zi(m1)采用下面的方法重构:以zi(m1)作为x轴坐标,以zi(m1+τ)作为y轴纵坐标,绘出第一重构点,其中,0<m1<m-τ;用一个40×40的网格覆盖所有第一重构点,被所述第一重构点占据的格子数与所有格子数的比值即为离散度参数Di
作为一种优选,所述变化度参数计算模块具体为:将所述第一窗信号yi(m)分为时间长度为m2的第一滑动窗,得到第一滑动窗信号hii(m2);对第一滑动窗信号hii(m2)进行归一化得到第一滑动窗归一化信号hi1(m2);通过将所述第一滑动窗归一化信号hi1(m2)与给定第一幅度阈值V比较,得到每个所述第一滑动窗归一化信号hi1(m2)的二进制数据序列;计算所有所述第一滑动窗归一化信号hi1(m2)的二进制数据序列中上升沿和下降沿的个数,进一步计算单个第一滑动窗的二进制数据序列b1b2b3Λbn中上升沿和下降沿的个数平均值即为变化度参数Ni
作为一种优选,所述偏离度参数计算模块具体包括:将所述第一窗信号yi(m)每t2分一个小窗信号;求得所述每个小窗信号波峰以及波峰所在的位置和波谷以及波谷所在的位置;将所有所述小窗信号波峰和波谷按时间顺序分别连续排列,如果所有所述小窗信号波峰的均值大于波谷的均值,则对波峰对应的时间间隔进行重构获得第二重构点,如果所有所述小窗信号波峰的均值小于波谷的均值,则对波谷对应的时间间隔进行重构获得第二重构点,即第二重构点的x轴横坐标为ti-ti-1(i=2,Λ,n),第二重构点的y轴纵坐标为ti+1-ti(i=2,Λ,n),ti代表连续排列中第i个波峰或波谷的位置;计算每个所述第二重构点重构点与原点连线之直线和直线y=x的直线夹角,其所有直线夹角相加总和即为偏离度参数Si
作为一种优选,所述可电击信号判断模块具体为:
对于每个第一窗信号yi(m),
如果P1≤Di≤P2,T1≤Ni≤T2,则
如果P1≤Di≤P2,T2<Ni≤T3,R1<Si,则
如果P1≤Di≤P2,T2<Ni≤T3,R1>Si,则
如果P1≤Di≤P2,T3<Ni≤T4,则
如果P1≤Di≤P2,T4<Ni<T5,R2<Si,则
如果P1≤Di≤P2,T4<Ni<T5,但不满足R2<Si,则
如果P1≤Di≤P2,Ni≥T5,则
如果P1≤Di≤P2,但不满足T1≤Ni≤T2和T2<Ni≤T3和T3<Ni≤T4和T4<Ni<T5和Ni≥T5,则
如果P2<Di,则
如果Di<P1,则
进一步计算所有所述第一窗信号yi(m)的的和sr1;
如果sr1大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号;如果sr1小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号。
作为一种优选,如果所述待检测心电信号为不可电击信号,所述可电击信号判断模块还包括二阶判断单元,在所述待检测心电信号为不可电击信号时,所述二阶判断单元进行如下判断::
对于每个第一窗信号yi(m),
如果R3≤Si≤R4,T6<Ni,则
如果R3≤Si≤R4,但不满足T6<Ni,则
如果Si<R4,但不满足R3≤Si≤R4,则
如果不满足R3≤Si≤R4和Si<R4,则
进一步计算所有第一窗信号yi(m)的的和sr2;
如果sr2大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号,结束;如果sr2小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号。
作为一种优选,所述可电击信号判断模块还包括三阶判断单元,如果所述二阶判断单元判断所述待检测心电信号为不可电击信号时,所述三阶判断单元还进行如下判断:
对于每个第一窗信号yi(m),
如果T7<Ni≤T8,T8<Ni,则
如果T7<Ni≤T8,但不满足T8<Ni,则
如果R5<Si,但不满足T7<Ni≤T8,则
如果不满足T7<Ni≤T8和R5<Si,则
进一步计算所有第一窗信号yi(m)的的和sr3;
如果sr3大于等于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号;如果sr3小于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号,结束。
本发明检测方法及装置通过采用离散度参数、变化度参数和偏离度参数检测,提高了可电击信号检测的灵敏度和准确率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明实施例整体检测方法流程图。
图2是本发明实施例检测方法流程图。
图3是本发明又一实施例检测方法部分流程图。
图4是本发明再一实施例检测方法部分流程图。
图5是本发明实施例检测装置框图。
具体实施方式
下面对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本检测方法只是对心电信号处理的中间结果可电击信号的检测。
如图1所示,一种可电击信号的检测方法,包括以下步骤:
S1、按固定的采样率,获得采集点数为n的初始心电信号;
S2、用滤波方法对初始心电信号进行预处理得到待检测心电信号y(n);
S3、用信号波峰与波谷差值的方法计算出待检测心电信号y(n)的幅值A,与设定的第二幅度阈值A1比较;
将时间长度为T0的窗滑动截取所述待检测心电信号y(n)得到若干第二窗信号。窗的宽度为T0,在时间轴上以1S的时间间隔滑动,得到若干个第二窗信号,取每个第二窗信号的波峰与波谷差,然后对所有第二窗信号的波峰与波谷差取平均,或者先舍弃所有波峰与波谷差中最大值和最小值再取平均,该平均即为待检测心电信号y(n)的幅值A。本实施例优选T0为2S。
S4、如果待检测心电信号y(n)的幅值A小于设定的第二幅度阈值A1,第二幅度阈值A1可以设定为0.2mv,则检测为停搏或幅值小于0.2mv的心电信号,进入S6、归为不可电击信号;如果A大于设定的阈值A1,则该信号进行步骤S5检测;
也可跳过S3、S4直接进入S5;
S5、将时间长度为m的窗,在时间轴上以t1的时间间隔滑动截取所述待检测心电信号y(n)得到若干第一窗信号yi(m),其中m、t1分别可以取4秒和1秒;例如第一个第一窗信号y1(m)的时间为范围(0,4),第二个第一窗信号y2(m)的时间为范围(1,5),以此类推。
S7、对所述第一窗信号yi(m)进行延时重构获得第一重构点,获取所述第一重构点的离散度参数Di;对所述第一窗信号yi(m)进行归一化处理得到二进制数据序列,获取所述二进制数据序列中上升沿和下降沿的变化度参数Ni;以所述第一窗信号yi(m)中相邻的波峰间或者相邻的波谷间的时间差值进行重构获得第二重构点,获取所述第二重构点的偏离度参数Si;所述获取离散度参数Di、变化度参数Ni和偏离度参数Si并无必然顺序关系;
将第一窗信号yi(m)降低采样率到50Hz,重新采样得到信号zi(m1);对信号zi(m1)进行下面的重构方法:x轴横坐标为zi(m1),y轴纵坐标为zi(m1+τ),其中τ=0.9。绘出相空间重构点;在图中产生一个40×40的网格延伸,这个网格覆盖了所有重构点,并计算被这些点所填充的格子。然后我们计算离散度参数Di
将第一窗信号yi(m)分为长度为m2(小窗)的滑动窗,得到信号hii(m2);通过绝对最大值的方法对hii(m2)进行归一化(除以最大值)得到hi1(m2);通过hi1(m2)与给定阈值V(0.2)的比较,得到0-1序列b1b2b3Λbn;计算0-1序列中上升沿和下降沿的个数Ni(其中上升沿和下降沿间隔小于n1的个数的不计,并且计一个下降沿后间隔n2点后再进行计算);进而计算所有滑动窗信号hii(m2)的平均上升沿和下降沿的个数,记为变化度参数Ni
将第一窗信号yi(m)每t(2s)分一个小窗,应用求局部极大值和局部极小值进而求得波峰以及波峰所在的位置和波谷以及波谷所在的位置。将求得的所有小窗的结果合并,进而求得第一窗信号yi(m)的波峰x(t1),x(t2),x(t3),Λ,x(ti-1),x(ti),x(ti-1),Λx(tn)和波谷x(t′1),(t′2),x(t′3)Λ,x(t′i-1),x(t′i),x(t′i+1),Λx(t′n),如果波峰绝对值的均值大于波谷绝对值的均值,则对波峰对应的时间隔进行重构,即x轴横坐标为ti-ti-1(i=2,Λ,n)。y轴纵坐标为ti+1-ti(i=2,Λ,n);否则对波谷对应的时间间隔做同正峰值的重构;计算所有重构点与原点连线和直线y=x的夹角的和即为偏离度参数Si
S8、计算sr1:
对于每个第一窗信号yi(m),
如果P1≤Di≤P2,T1≤Ni≤T2,则
如果P1≤Di≤P2,T2<Ni≤T3,R1<Si,则
如果P1≤Di≤P2,T2<Ni≤T3,R1>Si,则
如果P1≤Di≤P2,T3<Ni≤T4,则
如果P1≤Di≤P2,T4<Ni<T5,R2<Si,则
如果P1≤Di≤P2,T4<Ni<T5,但不满足R2<Si,则
如果P1≤Di≤P2,Ni≥T5,则
如果P1≤Di≤P2,但不满足T1≤Ni≤T2和T2<Ni≤T3和T3<Ni≤T4和T4<Ni<T5和Ni≥T5,则
如果P2<Di,则
如果Di<P1,则
进一步计算所有所述第一窗信号yi(m)的的和sr1;
S9、如果sr1大于等于三分之二所述第一窗信号yi(m)的个数,则进入S11、判定所述待检测心电信号为可电击信号,结束;否则;进入步骤S10;当然,在某些实施例中,如果sr1小于三分之二所述第一窗信号yi(m)的个数,也可以将所述待检测心电信号为不可电击信号,结束;
S10、计算sr2:
对于每个第一窗信号yi(m),
如果R3≤Si≤R4,T6<Ni,则
如果R3≤Si≤R4,但不满足T6<Ni,则
如果Si<R4,但不满足R3≤Si≤R4,则
如果不满足R3≤Si≤R4和Si<R4,则
进一步计算所有第一窗信号yi(m)的的和sr2;
S12、如果sr2大于等于三分之二所述第一窗信号yi(m)的个数,则进入S14、判定所述待检测心电信号为可电击信号,结束;否则,进入步骤S13;当然,在某些实施例中,如果sr2小于三分之二所述第一窗信号yi(m)的个数,也可以将所述待检测心电信号为不可电击信号,结束;
S13、计算sr3(参见图4):
对于每个第一窗信号yi(m),
如果T7<Ni≤T8,T8<Ni,则
如果T7<Ni≤T8,但不满足T8<Ni,则
如果R5<Si,但不满足T7<Ni≤T8,则
如果不满足T7<Ni≤T8和R5<Si,则
进一步计算所有第一窗信号yi(m)的的和sr3;
S15、如果sr3大于等于三分之二所述第一窗信号yi(m)的个数,则进入S17、判定所述待检测心电信号为可电击信号;否则进入S16、判定不可电击信号,结束;
S18、取n点新的数据。
通过变化度参数Ni可以检测出心室扑动心电信号,并且本变化度参数Ni和现有技术相比无需乘以W(t)函数,具有计算量小优点;离散度参数Di可以检测出心率偏低的心室纤颤信号。偏离度参数Si检测可以在中间平衡变化度参数Ni检测和离散度参数Di检测,可以保证算法检测的稳定。
实施例:
1.数据采集:按照采样率1000Hz,采集6s的心电信号存入预设的固定长度的存储空间获得原始心电信号x(n),对该6s原始心电信号进行分析,其中,n=6000。
2.对采集的原始心电信号x(n)进行预处理得到待检测心电信号y(n):使用一个滑动平均滤波器;再使用一个截止频率1Hz的高通滤波器;再使用一个截止频率为30Hz的巴特沃斯低通滤波器,得到待检测心电信号y(n)。
3.用信号波峰与波谷差值的方法计算出待检测心电信号y(n)的相对幅值A,与设定的阈值0.2比较,具体过程:
1)用时间长度为两秒的滑动窗截取所述待检测心电信号y(n)得到若干第二窗信号,得到每两秒第二窗信号的最大值maxi和最小值mini;
2)根据 A = ( Σ ( max i - min i ) - max i ( max i - min i ) - min i ( max i - min i ) ) / ( T - 2 ) , 其中T为第二窗个数,计算A;
4.设定阈值A1为0.2;
5.将时间长度为4s的窗滑动截取所述待检测心电信号y(n)得到3第一窗信号yi(m),计算出Di,Ni,Si
6.设定算法中的各个阈值:P1在0.07-0.08之间,P2在0.08-0.10之间,T1在6-8之间,T2在8-10之间,T3在9-11之间,T4在11-13之间,T5在12-15之间,R1在3-5之间,R2在2-5之间,R3在4-7之间,R4在8-10之间,T6在9-18之间,T7在9-13之间,T8在16-22之间,R5在5-8之间优选。
7.对于第一窗信号yi(m),如果P1≤Di≤P2,T1≤Ni≤T2,则ri=1;
如果P1≤Di≤P2,T2<Ni≤T3,R1<Si,则ri=1;
如果P1≤Di≤P2,T2<Ni≤T3,R1>Si,则ri=0;
如果P1≤Di≤P2,T3<Ni≤T4,则ri=1;
如果P1≤Di≤P2,T4<Ni<T5,R2<Si,则ri=1;
如果P1≤Di≤P2,T4<Ni<T5,但不满足R2<Si,则ri=0;
如果P1≤Di≤P2,Ni≥T5,则ri=1;
如果P1≤Di≤P2,但不满足T1≤Ni≤T2和T2<Ni≤T3和T3<Ni≤T4和T4<Ni<T5和Ni≥T5,则ri=0;
如果P2<Di,则ri=1;
如果Di<P1,则ri=0;
进一步计算所有所述第一窗信号yi(m)的ri的和sr1。
8.如果sr1大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号,结束;如果sr1小于三分之二所述第一窗信号yi(m)的个数,则进入下一步。
9.对于每个第一窗信号yi(m),
如果R3≤Si≤R4,T6<Ni,则ri=1;
如果R3≤Si≤R4,但不满足T6<Ni,则ri=0;
如果Si<R4,但不满足R3≤Si≤R4,则ri=1;
如果不满足R3≤Si≤R4和Si<R4,则ri=0;
进一步计算所有第一窗信号yi(m)的ri的和sr2。
10.如果sr2大于等于三分之二第一窗信号yi(m)的个数,则第一窗信号yi(m)心电信号为可电击信号,结束;如果sr2小于三分之二所述第一窗信号yi(m)的个数,则进入下一步。
11.对于每个第一窗信号yi(m),
如果T7<Ni≤T8,T8<Ni,则ri=1;
如果T7<Ni≤T8,但不满足T8<Ni,则ri=0;
如果R5<Si,但不满足T7<Ni≤T8,则ri=1;
如果不满足T7<Ni≤T8和R5<Si,则ri=0;
进一步计算所有第一窗信号yi(m)的ri的和sr3;
12.如果sr3大于等于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号;如果sr3小于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号,结束。
如图5所示,本发明还公开了一种可电击信号检测装置,包括数据获取模块1、窗截取模块2、离散度参数计算模块3、变化度参数计算模块4、偏离度参数计算模块5和可电击信号判断模块6;
所述数据获取模块1用于获得待检测心电信号y(n);
窗截取模块2用于将时间长度为m的窗滑动截取所述待检测心电信号y(n)得到若干第一窗信号yi(m);
离散度参数计算模块3用于对所述第一窗信号yi(m)进行延时重构获得第一重构点,获取第一重构点的离散度参数Di
变化度参数计算模块4用于对所述第一窗信号yi(m)进行归一化处理得到二进制数据序列,获取二进制数据序列中上升沿和下降沿的变化度参数Ni
偏离度参数计算模块5用于以所述第一窗信号yi(m)中相邻的波峰间或者相邻的波谷间的时间差值进行重构获得第二重构点,获取第二重构点的偏离度参数Si
可电击信号判断模块6用于根据所有窗信号的离散度参数Di、变化度参数Ni及偏离度参数Si,判断所述待检测心电信号是否为可电击信号。
在数据获取模块、窗截取模块之间还包括幅值判断模块,用于将所述待检测心电信号y(n)幅值A与预设的第二幅度阈值A1进行比较,如果所述待检测心电信号y(n)幅值A小于第二幅度阈值A1,所述待检测心电信号y(n)归为不可电击信号,结束;否则,所述待检测心电信号y(n)进入窗截取模块进行下一步判断。
所述离散度参数计算模块具体为:将所述第一窗信号yi(m)重新采样获得重采样信号zi(m1);对所述重采样信号zi(m1)采用下面的方法重构:以zi(m1)作为x轴坐标,以zi(m1+τ)作为y轴纵坐标,绘出第一重构点,其中,0<m1<m-τ;用一个40×40的网格覆盖所有第一重构点,被所述第一重构点占据的格子数与所有格子数的比值即为离散度参数Di
所述变化度参数计算模块具体为:将所述第一窗信号yi(m)分为时间长度为m2的第一滑动窗,得到第一滑动窗信号hii(m2);对第一滑动窗信号hii(m2)进行归一化得到第一滑动窗归一化信号hi1(m2);通过将所述第一滑动窗归一化信号hi1(m2)与给定第一幅度阈值V比较,得到每个所述第一滑动窗归一化信号hi1(m2)的二进制数据序列;计算所有所述第一滑动窗归一化信号hi1(m2)的二进制数据序列中上升沿和下降沿的个数,进一步计算单个第一滑动窗的二进制数据序列b1b2b3Λbn中上升沿和下降沿的个数平均值即为变化度参数Ni
所述偏离度参数计算模块具体包括:将所述第一窗信号yi(m)每t2分一个小窗信号;求得所述每个小窗信号波峰以及波峰所在的位置和波谷以及波谷所在的位置;将所有所述小窗信号波峰和波谷按时间顺序分别连续排列,如果所有所述小窗信号波峰的均值大于波谷的均值,则对波峰对应的时间间隔进行重构获得第二重构点,如果所有所述小窗信号波峰的均值小于波谷的均值,则对波谷对应的时间间隔进行重构获得第二重构点,即第二重构点的x轴横坐标为ti-ti-1(i=2,Λ,n),第二重构点的y轴纵坐标为ti+1-ti(i=2,Λ,n),ti代表连续排列中第i个波峰或波谷的位置;计算每个所述第二重构点与原点连线之直线和直线y=x的直线夹角,其所有直线夹角相加总和即为偏离度参数Si
所述可电击信号判断模块具体为:
对于每个第一窗信号yi(m),
如果P1≤Di≤P2,T1≤Ni≤T2,则
如果P1≤Di≤P2,T2<Ni≤T3,R1<Si,则
如果P1≤Di≤P2,T2<Ni≤T3,R1>Si,则
如果P1≤Di≤P2,T3<Ni≤T4,则
如果P1≤Di≤P2,T4<Ni<T5,R2<Si,则
如果P1≤Di≤P2,T4<Ni<T5,但不满足R2<Si,则
如果P1≤Di≤P2,Ni≥T5,则
如果P1≤Di≤P2,但不满足T1≤Ni≤T2和T2<Ni≤T3和T3<Ni≤T4和T4<Ni<T5和Ni≥T5,则
如果P2<Di,则
如果Di<P1,则
进一步计算所有所述第一窗信号yi(m)的的和sr1;
如果sr1大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号;如果sr1小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号。
所述可电击信号判断模块还包括二阶判断单元,在所述待检测心电信号为不可电击信号时,所述二阶判断单元进行如下判断:
对于每个第一窗信号yi(m),
如果R3≤Si≤R4,T6<Ni,则
如果R3≤Si≤R4,但不满足T6<Ni,则
如果Si<R4,但不满足R3≤Si≤R4,则
如果不满足R3≤Si≤R4和Si<R4,则
进一步计算所有第一窗信号yi(m)的的和sr2;
如果sr2大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号;如果sr2小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号。
所述可电击信号判断模块还包括三阶判断单元,如果所述二阶判断单元判断所述待检测心电信号为不可电击信号时,所述三阶判断单元还进行如下判断:
对于每个第一窗信号yi(m),
如果T7<Ni≤T8,T8<Ni,则
如果T7<Ni≤T8,但不满足T8<Ni,则
如果R5<Si,但不满足T7<Ni≤T8,则
如果不满足T7<Ni≤T8和R5<Si,则
进一步计算所有第一窗信号yi(m)的的和sr3;
如果sr3大于等于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为可电击信号;如果sr3小于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号为不可电击信号。
以上对本发明实施例所提供的方案,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (19)

1.一种可电击信号的检测方法,其特征在于,包括以下步骤:
1)获得待检测心电信号y(n);
2)将时间长度为m的窗滑动截取所述待检测心电信号y(n)得到若干第一窗信号yi(m),对每个所述第一窗信号yi(m)分别进行如下处理:
2.1)对所述第一窗信号yi(m)进行延时重构获得第一重构点,获取所述第一重构点的离散度参数Di
2.2)对所述第一窗信号yi(m)进行归一化处理得到二进制数据序列,获取所述二进制数据序列中上升沿和下降沿的变化度参数Ni
2.3)以所述第一窗信号yi(m)中相邻的波峰间或者相邻的波谷间的时间差值进行重构获得第二重构点,获取所述第二重构点的偏离度参数Si
3)根据所有第一窗信号yi(m)的离散度参数Di、变化度参数Ni及偏离度参数Si,判断所述待检测心电信号y(n)是否为可电击信号,具体步骤如下:
对于每个第一窗信号yi(m),
如果P1≤Di≤P2,T1≤Ni≤T2,则
如果P1≤Di≤P2,T2<Ni≤T3,R1<Si,则
如果P1≤Di≤P2,T2<Ni≤T3,R1>Si,则
如果P1≤Di≤P2,T3<Ni≤T4,则
如果P1≤Di≤P2,T4<Ni<T5,R2<Si,则
如果P1≤Di≤P2,T4<Ni<T5,但不满足R2<Si,则
如果P1≤Di≤P2,Ni≥T5,则
如果P1≤Di≤P2,但不满足T1≤Ni≤T2和T2<Ni≤T3和T3<Ni≤T4和T4<Ni<T5和Ni≥T5,则
如果P2<Di,则
如果Di<P1,则
进一步计算所有所述第一窗信号yi(m)的的和sr1;
如果sr1大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为可电击信号,结束;如果sr1小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为不可电击信号,结束,P1和P2为离散度参数的预设值,T1-T5为变化度参数的预设值,R1-R2为偏离度参数预设值,sr1为所有所述第一窗信号yi(m)的的和。
2.根据权利要求1所述可电击信号的检测方法,其特征在于,在所述步骤1)和2)之间还包括幅值判断步骤:将所述待检测心电信号y(n)幅值A与预设的第二幅度阈值A1进行比较,如果所述待检测心电信号y(n)幅值A小于第二幅度阈值A1,所述待检测心电信号y(n)归为不可电击信号,结束;否则,进入下一步。
3.根据权利要求2所述可电击信号的检测方法,其特征在于,所述待检测心电信号y(n)幅值A具体计算步骤如下:
3.1)将时间长度为T0的窗滑动截取所述待检测心电信号y(n)得到若干第二窗信号;
3.2)得到每个所述第二窗信号内的幅度最大值和最小值的差值;
3.3)所述待检测心电信号y(n)幅值A为所有所述第二窗信号幅度最大值和最小值的差值的平均值。
4.根据权利要求1所述可电击信号的检测方法,其特征在于,所述步骤2.1)中的离散度参数Di具体计算如下:
4.1)将所述第一窗信号yi(m)重新采样获得重采样信号zi(m1);
4.2)对所述重采样信号zi(m1)采用下面的方法重构:以zi(m1)作为x轴坐标,以zi(m1+τ)作为y轴纵坐标,绘出第一重构点,τ为延迟时间,其中,0<m1<m-τ;
4.3)用一个40×40的网格覆盖所有所述第一重构点,被所述第一重构点占据的格子数与所有格子数的比值即为离散度参数Di
5.根据权利要求4所述可电击信号的检测方法,其特征在于,所述τ为0.8~0.95,τ为延迟时间。
6.根据权利要求1所述的可电击信号的检测方法,其特征在于,所述步骤2.2)中变化度参数Ni具体计算如下:
6.1)将所述第一窗信号yi(m)分为时间长度为m2的第一滑动窗,得到第一滑动窗信号hii(m2);
6.2)对第一滑动窗信号hii(m2)进行归一化得到第一滑动窗归一化信号hi1(m2);
6.3)通过将所述第一滑动窗归一化信号hi1(m2)与给定第一幅度阈值V比较,得到每个所述第一滑动窗归一化信号hi1(m2)的二进制数据序列;
6.4)计算所有所述第一滑动窗归一化信号hi1(m2)的二进制数据序列中上升沿和下降沿的个数,进一步计算单个第一滑动窗归一化信号hi1(m2)的二进制数据序列中上升沿和下降沿的个数平均值即为变化度参数Ni
7.根据权利要求6所述可电击信号的检测方法,其特征在于,所述第一幅度阈值V为0.15~0.25,V为第一幅度阈值,具有归一化幅度单位。
8.根据权利要求6所述可电击信号的检测方法,其特征在于,如果所述第一滑动窗的归一化信号hi1(m2)的二进制数据序列中的上升沿或下降沿与相邻的上一个上升沿或下降沿间隔小于n1,则删除所述上升沿或下降沿,并且在一个下降沿后间隔n2点后再进行计算,n1和n2分别为预设的时间值。
9.根据权利要求1所述可电击信号的检测方法,其特征在于,所述步骤2.3)中偏离度参数Si具体计算如下:
9.1)将所述第一窗信号yi(m)每t2分一个小窗信号;
9.2)求得所述每个小窗信号波峰以及波峰所在的位置和波谷以及波谷所在的位置;
9.3)将所有所述小窗信号波峰和波谷按时间顺序分别连续排列,如果所有所述小窗信号波峰的均值大于波谷的均值,则对波峰对应的时间间隔进行重构获得第二重构点,如果所有所述小窗信号波峰的均值小于波谷的均值,则对波谷对应的时间间隔进行重构获得第二重构点;第二重构点的x轴横坐标为ti-ti-1(i=2,…,n),第二重构点的y轴纵坐标为ti+1-ti(i=2,…,n),ti代表连续排列中第i个波峰或波谷的位置;
9.4)计算每个所述第二重构点与原点连线之直线和直线y=x的直线夹角,其所有直线夹角相加总和即为偏离度参数Si
10.根据权利要求9所述可电击信号的检测方法,其特征在于,所述t2为1.5~2.5秒。
11.根据权利要求1所述可电击信号的检测方法,其特征在于,如果所述待检测心电信号y(n)为不可电击信号,则根据所有窗信号的离散度参数Di、变化度参数Ni及偏离度参数Si,进一步判断所述待检测心电信号y(n)是否为可电击信号,具体为:
对于每个第一窗信号yi(m),
如果R3≤Si≤R4,T6<Ni,则
如果R3≤Si≤R4,但不满足T6<Ni,则
如果Si<R4,但不满足R3≤Si≤R4,则
如果不满足R3≤Si≤R4和Si<R4,则
进一步计算所有第一窗信号yi(m)的的和sr2;
如果sr2大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为可电击信号,结束;如果sr2小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为不可电击信号,结束,T6为变化度参数的预设值,R3-R4为偏离度参数预设值,sr2为所有第一窗信号yi(m)的的和。
12.根据权利要求11所述可电击信号的检测方法,其特征在于,如果所述待检测心电信号y(n)再一次判断为不可电击信号,则根据所有窗信号的离散度参数Di、变化度参数Ni及偏离度参数Si,进一步判断所述待检测心电信号y(n)是否为可电击信号具体为:
对于每个第一窗信号yi(m),
如果T7<Ni≤T8,T8<Ni,则
如果T7<Ni≤T8,但不满足T8<Ni,则
如果R5<Si,但不满足T7<Ni≤T8,则
如果不满足T7<Ni≤T8和R5<Si,则
进一步计算所有第一窗信号yi(m)的的和sr3;
如果sr3大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为可电击信号;如果sr3小于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为不可电击信号,结束,T7-T8为变化度参数的预设值,R5为偏离度参数预设值,sr3为所有第一窗信号yi(m)的的和。
13.一种可电击信号的检测装置,其特征在于,包括数据获取模块、窗截取模块、离散度参数计算模块、变化度参数计算模块、偏离度参数计算模块和可电击信号判断模块;
所述数据获取模块用于获得待检测心电信号y(n);
窗截取模块用于将时间长度为m的窗滑动截取所述待检测心电信号y(n)得到若干第一窗信号yi(m);
离散度参数计算模块用于对所述第一窗信号yi(m)进行延时重构获得第一重构点,获取第一重构点的离散度参数Di
变化度参数计算模块用于对所述第一窗信号yi(m)进行归一化处理得到二进制数据序列,获取二进制数据序列中上升沿和下降沿的变化度参数Ni
偏离度参数计算模块用于以所述第一窗信号yi(m)中相邻的波峰间或者相邻的波谷间的时间差值进行重构获得第二重构点,获取第二重构点的偏离度参数Si
可电击信号判断模块用于根据所有第一窗信号yi(m)的离散度参数Di、变化度参数Ni及偏离度参数Si,判断所述待检测心电信号y(n)是否为可电击信号,具体步骤如下:
对于每个第一窗信号yi(m),
如果P1≤Di≤P2,T1≤Ni≤T2,则
如果P1≤Di≤P2,T2<Ni≤T3,R1<Si,则
如果P1≤Di≤P2,T2<Ni≤T3,R1>Si,则
如果P1≤Di≤P2,T3<Ni≤T4,则
如果P1≤Di≤P2,T4<Ni<T5,R2<Si,则
如果P1≤Di≤P2,T4<Ni<T5,但不满足R2<Si,则
如果P1≤Di≤P2,Ni≥T5,则
如果P1≤Di≤P2,但不满足T1≤Ni≤T2和T2<Ni≤T3和T3<Ni≤T4和T4<Ni<T5和Ni≥T5,则
如果P2<Di,则
如果Di<P1,则
进一步计算所有所述第一窗信号yi(m)的的和sr1;
如果sr1大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为可电击信号,结束;如果sr1小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为不可电击信号,结束,P1和P2为离散度参数的预设值,T1-T5为变化度参数的预设值,R1-R2为偏离度参数预设值,sr1为所有所述第一窗信号yi(m)的的和。
14.根据权利要求13所述可电击信号的检测装置,其特征在于,在数据获取模块、窗截取模块之间还包括幅值判断模块,用于将所述待检测心电信号y(n)幅值A与预设的第二幅度阈值A1进行比较,如果所述待检测心电信号y(n)幅值A小于第二幅度阈值A1,所述待检测心电信号y(n)归为不可电击信号,结束;否则,所述待检测心电信号y(n)进入窗截取模块进行下一步判断。
15.根据权利要求13所述可电击信号的检测装置,其特征在于,所述离散度参数计算模块具体为:将所述第一窗信号yi(m)重新采样获得重采样信号zi(m1);对所述重采样信号zi(m1)采用下面的方法重构:以zi(m1)作为x轴坐标,以zi(m1+τ)作为y轴纵坐标,绘出第一重构点,其中,0<m1<m-τ;用一个40×40的网格覆盖所有第一重构点,被所述第一重构点占据的格子数与所有格子数的比值即为离散度参数Di
16.根据权利要求13所述可电击信号的检测装置,其特征在于,所述变化度参数计算模块具体为:将所述第一窗信号yi(m)分为时间长度为m2的第一滑动窗,得到第一滑动窗信号hii(m2);对第一滑动窗信号hii(m2)进行归一化得到第一滑动窗归一化信号hi1(m2);通过将所述第一滑动窗归一化信号hi1(m2)与给定第一幅度阈值V比较,得到每个所述第一滑动窗归一化信号hi1(m2)的二进制数据序列;计算所有所述第一滑动窗归一化信号hi1(m2)的二进制数据序列中上升沿和下降沿的个数,进一步计算单个第一滑动窗的二进制数据序列b1b2b3…bn中上升沿和下降沿的个数平均值即为变化度参数Ni
17.根据权利要求13所述可电击信号的检测装置,其特征在于,所述偏离度参数计算模块具体包括:将所述第一窗信号yi(m)每t2分一个小窗信号;求得所述每个小窗信号波峰以及波峰所在的位置和波谷以及波谷所在的位置;将所有所述小窗信号波峰和波谷按时间顺序分别连续排列,如果所有所述小窗信号波峰的均值大于波谷的均值,则对波峰对应的时间间隔进行重构获得第二重构点,如果所有所述小窗信号波峰的均值小于波谷的均值,则对波谷对应的时间间隔进行重构获得第二重构点,即第二重构点的x轴横坐标为ti-ti-1(i=2,…,n),第二重构点的y轴纵坐标为ti+1-ti(i=2,…,n),ti代表连续排列中第i个波峰或波谷的位置;计算每个所述第二重构点与原点连线之直线和直线y=x的直线夹角,其所有直线夹角相加总和即为偏离度参数Si
18.根据权利要求13所述可电击信号的检测装置,其特征在于,所述可电击信号判断模块还包括二阶判断单元,在所述待检测心电信号y(n)为不可电击信号时,所述二阶判断单元进行如下判断:
对于每个第一窗信号yi(m),
如果R3≤Si≤R4,T6<Ni,则
如果R3≤Si≤R4,但不满足T6<Ni,则
如果Si<R4,但不满足R3≤Si≤R4,则
如果不满足R3≤Si≤R4和Si<R4,则
进一步计算所有第一窗信号yi(m)的的和sr2;
如果sr2大于等于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为可电击信号;如果sr2小于三分之二所述第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为不可电击信号,T6为变化度参数的预设值,R3-R4为偏离度参数预设值,sr2为所有第一窗信号yi(m)的的和。
19.根据权利要求18所述可电击信号的检测装置,其特征在于,所述可电击信号判断模块还包括三阶判断单元,如果所述二阶判断单元判断所述待检测心电信号y(n)为不可电击信号时,所述三阶判断单元还进行如下判断:
对于每个第一窗信号yi(m),
如果T7<Ni≤T8,T8<Ni,则
如果T7<Ni≤T8,但不满足T8<Ni,则
如果R5<Si,但不满足T7<Ni≤T8,则
如果不满足T7<Ni≤T8和R5<Si,则
进一步计算所有第一窗信号yi(m)的的和sr3;
如果sr3大于等于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为可电击信号;如果sr3小于三分之二第一窗信号yi(m)的个数,则所述待检测心电信号y(n)为不可电击信号,T7-T8为变化度参数的预设值,R5为偏离度参数预设值,sr3为所有第一窗信号yi(m)的的和。
CN201210593464.9A 2012-12-31 2012-12-31 一种可电击信号的检测方法和装置 Active CN103007431B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210593464.9A CN103007431B (zh) 2012-12-31 2012-12-31 一种可电击信号的检测方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210593464.9A CN103007431B (zh) 2012-12-31 2012-12-31 一种可电击信号的检测方法和装置

Publications (2)

Publication Number Publication Date
CN103007431A CN103007431A (zh) 2013-04-03
CN103007431B true CN103007431B (zh) 2015-05-20

Family

ID=47956884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210593464.9A Active CN103007431B (zh) 2012-12-31 2012-12-31 一种可电击信号的检测方法和装置

Country Status (1)

Country Link
CN (1) CN103007431B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104274907B (zh) * 2014-09-19 2016-08-31 中国人民解放军第三军医大学 一种除颤仪
CN104382586A (zh) * 2014-11-28 2015-03-04 深圳邦健生物医疗设备股份有限公司 一种电击信号检测方法和装置
CN108670245A (zh) * 2018-05-31 2018-10-19 哈尔滨工业大学深圳研究生院 一种心电信号检测方法及系统
CN114138454B (zh) * 2021-10-19 2022-07-15 中标慧安信息技术股份有限公司 一种边缘计算平台算力的公平分配方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000047279A1 (en) * 1999-02-08 2000-08-17 Heartsine Technologies, Inc. Apparatus for determining when a patient is susceptible to defibrillation
US6480734B1 (en) * 2000-06-30 2002-11-12 Cardiac Science Inc. Cardiac arrhythmia detector using ECG waveform-factor and its irregularity
US20020165585A1 (en) * 2001-05-01 2002-11-07 Dupelle Michael R. Pulse sensors
US7894894B2 (en) * 2006-03-29 2011-02-22 Medtronic, Inc. Method and apparatus for detecting arrhythmias in a subcutaneous medical device
US20080147133A1 (en) * 2006-12-15 2008-06-19 Garth Garner Implantable medical device
CN102389302B (zh) * 2011-07-20 2013-07-10 哈尔滨工业大学深圳研究生院 一种心电信号动力特征分析方法

Also Published As

Publication number Publication date
CN103007431A (zh) 2013-04-03

Similar Documents

Publication Publication Date Title
CN103007431B (zh) 一种可电击信号的检测方法和装置
CN103869162B (zh) 一种基于时域准同步的动态信号相量测量方法
CN104658100B (zh) 卷钞轮式暂存部纸币残留检测方法、装置及atm机
CN105078444B (zh) 一种噪声检测方法、装置及医疗检测设备
CN105105737A (zh) 基于光电容积描记和谱分析的运动状态心率监测方法
CN107329049B (zh) 一种基于卡尔曼滤波器的输电线路早期故障检测方法
CN110236573A (zh) 心理压力状态的检测方法及相关装置
CN106405230B (zh) 频率测量方法和装置
CN104545887A (zh) 伪差心电波形识别方法和装置
CN114142464B (zh) 基于特征电流信号检测的低压台区拓扑识别方法和装置
CN104720798A (zh) 基于脑电频谱特征的疲劳检测方法及系统
CN109350031A (zh) 一种光电容积脉搏波信号特征点检测方法及装置
CN102809687A (zh) 一种交流电频率的数字化测量方法
CN104219008B (zh) 一种宽带频谱检测方法和装置
CN105559771B (zh) 一种心电信号质量检测的方法及装置
CN103823177A (zh) 基于窗函数设计的滤波器性能检测方法和系统
CN103559714A (zh) 一种异常检测方法及装置
CN106073800A (zh) 基于绝对差值和提取的动态光谱数据处理方法及其装置
CN103575979A (zh) 一种数字化测量交流电频率的方法
CN106970265B (zh) 一种采用多时间尺度不完全s变换估计谐波参数的方法
CN104617914B (zh) 一种电器设备的信号滤波方法及系统
CN103575981A (zh) 一种交流电频率的精确测量方法
CN104048680B (zh) 基于donoho阈值的自主式水下机器人外部干扰抑制方法
CN103071236B (zh) 一种电击信号检测方法及装置
CN106357300B (zh) Bpsk扩频系统抗脉冲干扰方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant