CN105078447B - 一种心电信号r波定位方法 - Google Patents
一种心电信号r波定位方法 Download PDFInfo
- Publication number
- CN105078447B CN105078447B CN201510398092.8A CN201510398092A CN105078447B CN 105078447 B CN105078447 B CN 105078447B CN 201510398092 A CN201510398092 A CN 201510398092A CN 105078447 B CN105078447 B CN 105078447B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- electrocardiosignal
- msup
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000004807 localization Effects 0.000 title claims description 11
- 238000005070 sampling Methods 0.000 claims description 9
- 230000002708 enhancing effect Effects 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000000718 qrs complex Methods 0.000 description 12
- 238000000605 extraction Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 6
- 230000004217 heart function Effects 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000000747 cardiac effect Effects 0.000 description 3
- 210000004413 cardiac myocyte Anatomy 0.000 description 3
- 230000028161 membrane depolarization Effects 0.000 description 3
- 230000002336 repolarization Effects 0.000 description 3
- 230000009885 systemic effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000029052 metamorphosis Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明针对滤波后的心电信号,通过采用插值多项式法对人体心电信号进行拟合,并通过两次曲线微分增强心电信号的峰值,由此确定R波的位置。本发明的方法不仅计算简单、易于实现,而且准确率高,为心电信号快速有效识别提供新途径。
Description
技术领域
该发明是一种心电信号R波定位方法,主要针对人体滤波后的心电信号,采用多项式拟合心电曲线,并对拟合曲线进行两次微分以增强心电信号的峰值,从而确定R波的位置,实现心电信号R波的准确定位,以及R波峰值和RR间期特征的识别。
背景技术
心电信号(Electrocardiogram,ECG)记录了心脏在每个心动周期中,由起搏点、心房、心室相继兴奋,伴随的电能量的变化,其形态变化反映了人体心脏功能,是反应心脏状态的参考依据,准确识别心电信号的变化特征显得至关重要。
心电信号的幅度范围大约在10μV-5mV,频率范围大约在0.05-100Hz。典型的心电信号波形如图1所示:心电信号是由许多特征部分组成,分别为:P波、Q波、R波、S波、T波,其中Q波、R波、S波同属于QRS波群中。QRS波群起点对应于Q波峰之前的第一个斜率突变点,QRS波群的终点对应于S波峰之后的第一个斜率突变点。QRS波群的终点到T波起点的一段称为ST段。在心电信号特征波形的定位与检测过程中,QRS复合波群在心电图中表现为其R波波峰峰值明显大于其他位置幅值,且其在心电信号中占据时间短。通常QRS波群的定位是心电特征提取的关键,而QRS波群的确定首先要先识别R波的波峰,R波是QRS波群中的正向波峰,相较于P波和T波,具有幅值高、上升下降斜率大的特点。R波作为心电信号最显著特征,其定位正确与否直接关系到QRS波群、ST段位移、T波等其它特征的正确提取。
目前常见的R波定位方法差分阈值法、峰值检测法、可变斜率阈值法等,这类方法简单易于实现,但是由于检测时干扰的存在和预先设定阈值的不准确,通常R波错误识别率较高。此外,现在还有神经网络、小波分析、模板匹配等先进的信号处理方法用于R波的定位,但是这些方法都比较复杂,计算量大,不适合心电信号的实时处理。对于实际非标准的动态心电信号R波识别效果并不理想,所以,随着心电监护社会的需求日益增加,有必要研究一种新的计算简单、适应性强的R波定位方法,实现心电信号特征的快速准确识别。
发明内容
人体心电信号反映了心肌细胞有规律的除极复极过程,准确识别心电信号特征对了解心脏功能具有重要意义。采用心电检测仪采集出的人体心脏信号S(t),通常经过滤波方法去除基线漂移、工频干扰、肌电干扰等检测时参杂的干扰信号,可以获得较为清晰的心电信号然后采用特征提取方法提取心电信号特征,以反映心脏功能的强弱。
本发明针对滤波后的心电信号通过采用多项式对人体心电信号进行拟合,并通过对拟合曲线的两次微分增强心电信号的峰值,由此确定R波的位置。最后实现R波的准确定位,R波峰值和RR间期特征的提取,为心电信号其它特征的识别奠定基础。该方法不仅计算简单、易于实现,而且准确率高,也为心电信号快速有效识别提供了新途径。
本发明提供一种基于多项式拟合曲线二次微分的心电信号R波定位方法,其中:
设函数s(x)为滤波后的心电信号,xi为采样点时间,fi为采样点xi心电信号幅值,其中i=0,1,2...,N,
基于多项式拟合曲线二次微分的心电信号R波定位方法具体步骤如下:
设x为心电信号当前采样点时间,心电信号采样时间间隔为Δt,采样时间点x1、x2分别表示x+Δt、x+2Δt,x-1、x-2为分别表示x-Δt、x-2Δt,f-1、f-2、f1、f2分别为时间x-1、x-2、x1、x2对应心电信号的幅值,采用插值多项式s(x)拟合心电信号s(x)可得:
s(x)=L-2(x)f-2+L-1(x)f-1+L1(x)f1+L2(x)f2,
其中
为增强心电信号峰值,计算多项式近似曲线s(x)的斜率s'(x):
其中
设s'(x)在时间x-1、x-2、x1、x2对应心电信号的斜率为f′-1、f′-2、f1′、f2′,对心电信号斜率曲线s'(x)进行多项式插值计算并求取斜率,从而获得心电信号s(x)经过二次插值计算和微分后的峰值增强曲线s”(x):
s”(x)=λ-2f’-2+λ-1f’-1+λ1f1'+λ2f2',
设定R波峰值阈值为θ,当s”(x)曲线采样点xi的幅值s″(xi)>θ时,xi便可确定为R峰的位置xR。
此外,为了简化微分计算,令实验采样点时间为:
xi=i·Δt,
则
所以,获得采样点x0处心电拟合曲线的斜率为:
采样点x0处二次峰值增强后幅值为:
本发明的方法不仅计算简单、易于实现,而且准确率高,为心电信号快速有效识别提供新途径。
附图说明
附图1为典型心电信号波形。
具体实施方式
人体心电信号反映了心肌细胞有规律的除极复极过程,采用心电检测仪采集出人体心脏信号S(t),并经过滤波方法去除基线漂移、工频干扰、肌电干扰等检测时参杂的干扰信号,获得较为清晰的心电信号然后采用特征提取方法提取心电信号QRS波、RR间期、S波、T波、ST段位移、ST段斜率等特征,以反映心脏功能的强弱。
人体心电信号反映了心肌细胞有规律的除极复极过程,采用心电检测仪采集出人体心脏信号S(t),并经过滤波方法去除基线漂移、工频干扰、肌电干扰等检测时参杂的干扰信号,获得较为清晰的心电信号然后采用特征提取方法提取心电信号QRS波、RR间期、S波、T波、ST段位移、ST段斜率等特征,以反映心脏功能的强弱。
本发明主要针对滤波后心电信号的R波位置、R波波峰的特征提取提出了新方法,采用对多项式拟合心电曲线,并对拟合曲线进行二次微分,以增强QRS波群,并减弱P波和T波,然后根据斜率大点确定QRS波群位置,最后实现R波定位、R波波峰和RR间期特征识别。这种方法既没有使用复杂的数学计算,又容易实现。具体过程如下:
设s(x)为滤波后的人体心电信号,xi为采样点时间,fi为采样点xi心电信号幅值,i=0,1,2...,N。基于多项式拟合曲线二次微分的心电信号R波定位方法具体步骤如下:①设x为心电信号当前采样点时间,心电信号采样时间间隔为Δt,采样时间点x1、
x2分别表示x+Δt、x+2Δt,同样地,x-1、x-2为分别表示x-Δt、x-2Δt;再设f-1、
f-2、f1、f2分别为时间x-1、x-2、x1、x2对应心电信号的幅值。采用插值多项式s(x)
拟合心电信号s(x)可得:
s(x)=L-2(x)f-2+L-1(x)f-1+L1(x)f1+L2(x)f2 (1)
其中
为增强心电信号峰值,根据公式(1)计算多项式近似曲线s(x)的斜率s'(x):
其中
为了简化微分计算,令实验采样点时间为:
xi=i·Δt
故有:
所以,获得采样点x0处心电拟合曲线的斜率为:
②由于R峰是心电信号QRS波群的斜率变化最剧烈的地方,为对R峰进行准确定位,设s'(x)在时间x-1、x-2、x1、x2对应心电信号的斜率为f′-1、f′-2、f1′、f2′,再根据公式(2)对心电信号斜率曲线s'(x)进行多项式插值计算并求取斜率,从而获得心电信号s(x)经过二次插值计算和微分后的峰值增强曲线s”(x):
s”(x)=λ-2f’-2+λ-1f’-1+λ1f1'+λ2f2' (4)
并由此获得采样点x0处二次峰值增强后幅值为:
③设定R波峰值阈值为θ,当s”(x)曲线采样点xi的幅值s″(xi)>θ时,xi便可确定为R峰的位置xR,从而在心电信号s(x)中找到其对应的R峰峰值fR,以及两个R峰之间的时间间隔,即RR间期TRR。
上述说明示出并描述了本发明的若干优选实施例,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (2)
1.一种基于多项式拟合曲线二次微分的心电信号R波定位方法,其特征在于:
设函数s(x)为滤波后的心电信号,xi为采样点时间,fi为采样点xi心电信号幅值,其中i=0,1,2...,N,
基于多项式拟合曲线二次微分的心电信号R波定位方法具体步骤如下:
设x为心电信号当前采样点时间,心电信号采样时间间隔为Δt,采样时间点x1、x2分别表示x+Δt、x+2Δt,x-1、x-2为分别表示x-Δt、x-2Δt,f-1、f-2、f1、f2分别为时间x-1、x-2、x1、x2对应心电信号的幅值,采用插值多项式s(x)拟合心电信号s(x)可得:
s(x)=L-2(x)f-2+L-1(x)f-1+L1(x)f1+L2(x)f2,
其中
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
为增强心电信号峰值,计算多项式近似曲线s(x)的斜率s'(x):
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
</mrow>
其中
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>dL</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>dL</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>dL</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>dL</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
设s'(x)在时间x-1、x-2、x1、x2对应心电信号的斜率为f′-1、f′-2、f′1、f′2,对心电信号斜率曲线s'(x)进行多项式插值计算并求取斜率,从而获得心电信号s(x)经过二次插值计算和微分后的峰值增强曲线s”(x):
s”(x)=λ-2f′-2+λ-1f′-1+λ1f′1+λ2f′2,
设定R波峰值阈值为θ,当s”(x)曲线采样点xi的幅值s″(xi)>θ时,xi便可确定为R峰的位置xR。
2.如权利要求1所述的心电信号R波定位方法,其特征在于:
令实验采样点时间为:
xi=i·Δt,
则
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>12</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mn>3</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mn>3</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>12</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
所以,获得采样点x0处心电拟合曲线的斜率为:
<mrow>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>12</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>8</mn>
<msub>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>8</mn>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
采样点x0处二次峰值增强后幅值为:
<mrow>
<msup>
<mi>s</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>12</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<mn>8</mn>
<msubsup>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>+</mo>
<mn>8</mn>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510398092.8A CN105078447B (zh) | 2015-07-08 | 2015-07-08 | 一种心电信号r波定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510398092.8A CN105078447B (zh) | 2015-07-08 | 2015-07-08 | 一种心电信号r波定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105078447A CN105078447A (zh) | 2015-11-25 |
CN105078447B true CN105078447B (zh) | 2017-12-05 |
Family
ID=54560438
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510398092.8A Expired - Fee Related CN105078447B (zh) | 2015-07-08 | 2015-07-08 | 一种心电信号r波定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105078447B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113499082B (zh) * | 2020-03-23 | 2023-12-19 | 疆域康健创新医疗科技成都有限公司 | Qrs波群检测方法、心电检测装置及可读存储介质 |
CN113100778B (zh) * | 2021-04-14 | 2022-03-08 | 福州大学 | 一种心电信号的r峰定位方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101856225A (zh) * | 2010-06-30 | 2010-10-13 | 重庆大学 | 一种心电信号r波峰检测方法 |
CN101947112A (zh) * | 2010-09-20 | 2011-01-19 | 天津工业大学 | 心电监护物联网应用层心电信号综合检测算法实现 |
CN102379694A (zh) * | 2011-10-12 | 2012-03-21 | 中国科学院苏州纳米技术与纳米仿生研究所 | 心电图r波检测方法 |
CN103190901A (zh) * | 2013-04-01 | 2013-07-10 | 天津工业大学 | 基于极值域均值模式分解和改进Hilbert包络的R波检测算法 |
CN103750835A (zh) * | 2013-11-12 | 2014-04-30 | 天津工业大学 | 一种心电信号特征检测算法 |
CN104161509A (zh) * | 2014-08-08 | 2014-11-26 | 申岱 | 一种基于幅值谱的心率变异性分析方法及仪器 |
CN104173043A (zh) * | 2014-09-04 | 2014-12-03 | 东莞理工学院 | 一种适合于移动平台的心电数据分析方法 |
CN104586384A (zh) * | 2013-10-30 | 2015-05-06 | 深圳华清心仪医疗电子有限公司 | 一种心电信号实时心率检测方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03202048A (ja) * | 1989-12-28 | 1991-09-03 | Sony Corp | 信号波形のピークのサーチ方法 |
JP2011147583A (ja) * | 2010-01-21 | 2011-08-04 | Seiko Epson Corp | 心電計測装置 |
US8855752B2 (en) * | 2010-02-09 | 2014-10-07 | Angel Medical Systems, Inc. | QRS onset and offset detection with adaptive temporal windowing |
-
2015
- 2015-07-08 CN CN201510398092.8A patent/CN105078447B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101856225A (zh) * | 2010-06-30 | 2010-10-13 | 重庆大学 | 一种心电信号r波峰检测方法 |
CN101947112A (zh) * | 2010-09-20 | 2011-01-19 | 天津工业大学 | 心电监护物联网应用层心电信号综合检测算法实现 |
CN102379694A (zh) * | 2011-10-12 | 2012-03-21 | 中国科学院苏州纳米技术与纳米仿生研究所 | 心电图r波检测方法 |
CN103190901A (zh) * | 2013-04-01 | 2013-07-10 | 天津工业大学 | 基于极值域均值模式分解和改进Hilbert包络的R波检测算法 |
CN104586384A (zh) * | 2013-10-30 | 2015-05-06 | 深圳华清心仪医疗电子有限公司 | 一种心电信号实时心率检测方法及系统 |
CN103750835A (zh) * | 2013-11-12 | 2014-04-30 | 天津工业大学 | 一种心电信号特征检测算法 |
CN104161509A (zh) * | 2014-08-08 | 2014-11-26 | 申岱 | 一种基于幅值谱的心率变异性分析方法及仪器 |
CN104173043A (zh) * | 2014-09-04 | 2014-12-03 | 东莞理工学院 | 一种适合于移动平台的心电数据分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105078447A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105030233B (zh) | 一种心电信号st段的识别方法 | |
CN101919695B (zh) | 一种基于小波变换的心电信号qrs波检测方法 | |
CN107041743A (zh) | 一种心电信号实时r波检测方法 | |
CN105640545B (zh) | 一种胎儿心电信号提取方法及装置 | |
CN104545887B (zh) | 伪差心电波形识别方法和装置 | |
CN103393435B (zh) | 一种胎心音信号包络的取得方法及装置 | |
CN110477906B (zh) | 一种心电信号qrs波起止点定位方法 | |
CN109381181A (zh) | 心电信号特征波形的端点检测方法 | |
CN104706349A (zh) | 一种基于脉搏波信号的心电信号构建方法 | |
CN105054926A (zh) | 心电信号特征信息的提取方法及装置 | |
CN103584854A (zh) | 心电信号r波的提取方法 | |
CN108427031B (zh) | 基于多项式拟合及非干扰区域划分的间谐波检测方法 | |
CN103829944B (zh) | 基于模式识别的胸阻抗信号处理方法 | |
Manikandan et al. | Straightforward and robust QRS detection algorithm for wearable cardiac monitor | |
CN106154035A (zh) | 一种快速谐波及间谐波检测方法 | |
CN104204823B (zh) | 用于估算电池与电接地极之间的绝缘电阻的方法和系统 | |
CN105078447B (zh) | 一种心电信号r波定位方法 | |
CN102323977B (zh) | 基于心电特征点的心电数据存储方法 | |
US20200100697A1 (en) | Methods and system for processing an emg signal | |
CN103190901A (zh) | 基于极值域均值模式分解和改进Hilbert包络的R波检测算法 | |
Pangerc et al. | Robust detection of heart beats in multimodal records using slope-and peak-sensitive band-pass filters | |
CN104367315A (zh) | 一种提取心电图中j波信号的盲源分离方法 | |
CN110596458B (zh) | Demon谱谐波线谱和基频自动估计方法 | |
CN113499082A (zh) | Qrs波群检测方法、心电检测装置及可读存储介质 | |
Debbabi et al. | Correction of ECG baseline wander application to the Pan & Tompkins QRS detection algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171205 |
|
CF01 | Termination of patent right due to non-payment of annual fee |