一种心电信号R波定位方法、装置、终端及存储介质

一种心电信号r波定位方法、装置、终端及存储介质
技术领域
1.本发明涉及心电信号r波识别技术领域，尤其涉及一种心电信号r波定位方法、装置、终端及存储介质。


背景技术：

2.正常健康的人体心电波形图由p波、qrs波群以及t波组成，其中，qrs波群由向下的q波、尖高向上的r波以及与r波相连并向下的s波组成。在ecg中，r波是最明显的特征，确定r波的位置之后，可以以此为基准，寻找其他波形的位置，因此，确定r波峰值位置是ecg检测的关键。
3.目前实现r波峰值位置识别的方法采用小波变换法，其中，小波变换法主要是在分解后得到的高频分量上通过找到小波变换模极大值对，通过该极值对来定位r波的位置。
4.但是小波变换法实现时，通常需要对心电信号中的高频心电信号分量进行数据扩展，使其数据长度与原始心电信号数据长度一致，这样就会增加检测时间，对r波峰值的识别带来困难。


技术实现要素：

5.本发明提供了一种心电信号r波定位方法、装置、终端及存储介质，以解决现有技术中r波峰值定位速度慢的问题。
6.第一方面，本发明提供了一种心电信号r波定位方法，包括：
7.获取待识别的心电信号中高频心电信号分量的待检测高频数据序列的极大值序列和极小值序列，其中，所述极大值序列包括多个极大值，所述极小值序列包括多个极小值，所述极大值和所述极小值均包括位置和峰值；
8.根据所述极大值序列计算极大值序列阈值，根据所述极小值序列计算极小值序列阈值，并根据所述极大值序列阈值和所述极大值序列的峰值，获得真极大值序列，根据所述极小值序列阈值和所述极小值序列的峰值，获得真极小值序列；
9.对所述真极大值序列和所述真极小值序列进行极值对匹配，定位r波峰值范围；
10.计算所述r波峰值范围内每个极值对的极大值和极小值之间的过零点，组成过零点集合，并将所述过零点集合映射到所述心电信号，定位r波峰值位置。
11.第二方面，本发明提供了一种心电信号r波定位装置，包括：
12.获取模块，用于获取待识别的心电信号中的高频心电信号分量的待检测高频数据序列的极大值序列和极小值序列，其中，所述极大值序列包括多个极大值，所述极小值序列包括多个极小值，所述极大值和所述极小值均包括位置和峰值；
13.计算模块，用于根据所述极大值序列计算极大值序列阈值，根据所述极小值序列计算极小值序列阈值，并根据所述极大值序列阈值和所述极大值序列的峰值，获得真极大值序列，根据所述极小值序列阈值和所述极小值序列的峰值，获得真极小值序列；
14.第一定位模块，用于对所述真极大值序列和所述真极小值序列进行极值对匹配，
定位r波峰值范围；
15.第二定位模块，用于计算所述r波峰值范围内每个极值对的极大值和极小值之间的过零点，组成零点集合，并将所述零点集合映射到所述心电信号，定位r波峰值位置。
16.第三方面，本发明提供了一种终端，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上第一方面或第一方面的任一种可能的实现方式所述方法的步骤。
17.第四方面，本发明提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上第一方面或第一方面的任一种可能的实现方式所述方法的步骤。
18.本发明提供一种心电信号r波定位方法、装置、终端及存储介质，与现有技术相比，本发明所采用的高频心电信号分量不需要进行数据扩展，可以降低计算量，降低检测时间，同时，通过获得真极大值序列和真极小值序列，并进行极值对匹配，粗定位r波峰值位置，再通过计算过零点，映射到心电信号上，精准定位r波峰值位置，进一步提高r波峰值定位的正确度。
附图说明
19.为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
20.图1是本发明实施例提供的心电信号r波定位方法的实现流程图；
21.图2是本发明实施例提供的心电信号r波定位方法的极值序列图；
22.图3是本发明实施例提供的心电信号r波定位方法的极值序列图；
23.图4是本发明实施例提供的心电信号r波定位方法的极值序列图；
24.图5是本发明实施例提供的心电信号r波定位方法的极值序列图；
25.图6是本发明实施例提供的心电信号r波定位方法的过零点计算图；
26.图7是本发明实施例提供的心电信号r波定位装置的结构示意图；
27.图8是本发明实施例提供的终端的示意图。
具体实施方式
28.以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本发明实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本发明的描述。
29.为使本发明的目的、技术方案和优点更加清楚，下面将结合附图通过具体实施例来进行说明。
30.图1为本发明实施例提供的心电信号r波定位方法的实现流程图，详述如下：
31.在s101中，获取待识别的心电信号中高频心电信号分量的待检测高频数据序列的极大值序列和极小值序列，其中，极大值序列包括多个极大值，极小值序列包括多个极小
值，极大值和极小值均包括位置和峰值。
32.在一种可能的实现方式中，s101的具体实现流程包括：
33.根据预设极限漂移频率和心电信号的采样频率，确定心电信号的分解级数；
34.将心电信号分为n个心电信号分组，获取每个心电信号分组的预设噪声值，其中，n为大于等于2的正整数；
35.根据分解级数和每个心电信号分组的预设噪声值，对各个心电信号分组的心电信号进行分解，得到高频心电信号分量；
36.获取高频心电信号分量的待检测高频数据序列的极大值序列和极小值序列。
37.其中，基线漂移(baseline drift)是指检测器基线随时间的增加朝单一方向的偏离。
38.本发明实施例中，需要先对待识别的心电信号进行预处理，得到高频心电信号分量，再获取高频心电信号分量的待检测高频数据序列的极大值序列和极小值序列，其中，对待识别的心电信号进行预处理的实现过程包括：
39.第一步，对心电信号的采样频率f进行第一次分解，得到第一频率分量对第一频率分量进行第二次分解，得到第二频率分量
……
对第i频率分量进行第i 1次分解，得到第i 1频率分量其中，i为大于等于2的正整数；当第i 1频率分量小于预设基线漂移频率，对第i 1频率分量停止分解；将心电信号的分解级数l确定为i 1。
40.举例说明，预设基线漂移频率为2hz，预先采集的心电信号的采集频率范围为0～360hz，对心电信号进行分解。当分解级数为8时，此时的频率分量为1.4hz，此时的频率分量1.4hz小于预设基线漂移频率2hz，从而确定信号的分解级数为8级为最合适。此时的频率信号范围为0～1.4hz。
41.第二步，获取第一预设值、第二预设值和分组对应的标准差，其中，第一预设值为一固定值，且所有心电信号分组的第一预设值可以是相同的；当第一预设值大于等于该标准差时，将第一预设值确定为每个心电信号分组的预设噪声值；当第一预设值小于该标准差时，将第二预设值确定为每个心电信号分组的预设噪声值。
42.需要说明的是，第一预设值和第二预设值分别为不同的固定值。
43.例如：假设心电信号的数据为，d＝{d1,d2,
…
，dn}，d＝2k(k∈n

)，将心电信号的数据分为多个分组，其中，每t个心电信号的数据分为一组，需要说明的是为t个不重叠的心电信号的数据，其中，t∈{21,22,
…
,2k}，t可以选为4。从而将整个心电信号的数据分为n/t组。
44.首先，依次计算目标分组对应的标准差s，首先，依次计算目标分组对应的标准差s，其中，目标分组为多个分组中的任意一组。
45.然后，计算所有分组对应的标准差std的平均值std_mean，然后，计算所有分组对应的标准差std的平均值std_mean，其中，n为心电信号的数据的个数，j和i为自然数。
46.之后，设心电信号的所有数据的第一预设值为thr＝β
×
std_mean，其中，β为调节因子，β的范围1.0-2.0。需要说明的是，本发明实施例中，心电信号的所有数据具有相同的第一预设值thr。
47.最后，根据第一预设值thr和目标分组对应的标准差的大小，确定目标分组对应的预设噪声值δ。具体的，当第一预设值thr大于或等于目标分组对应的标准差的情况下，则将第一预设值thr设为目标分组对应的预设噪声值δ。当第一预设值thr小于目标分组对应的标准差的情况下，将第二预设值设为分组对应的预设噪声值δ。具体的，第二预设值可以根据具体情况设定，本发明实施例中，第二预设值可以设为0。
48.第三步，具体的，本发明中采用归一化一维转换压缩的方法，对心电信号的所有分组中的数据进行分解。
49.首先，根据每个心电信号分组对应的预设噪声值和预设数据扩展模型，对每个心电信号分组内的目标数据进行数据扩展，得到目标扩展数据，其中，目标数据为每个心电信号分组中的任意一个数据。
50.例如：给定的心电信号的初始行向量d＝{d1,d2,
…
,dn}，d＝2k(k∈n

)为心电信号初始行向量的数据个数，根据上述步骤中确定的预设噪声值δ，按照预设数据扩展模型，对目标数据进行数据扩展，得到目标扩展数据。
51.其中，预设数据扩展模型为：
[0052][0053]
其中，si为心电信号中的任意一个信号数据，δ为预设噪声值，为目标低频信号分量。
[0054]
其次，根据预设数据更新模型，对目标扩展数据进行数据更新，得到目标低频信号分量和目标高频信号分量。具体为：
[0055]
对相邻两区间的目标数据和进行处理，其中，i为奇数。依次对相邻两区间使用预设数据更新模型，对目标扩展数据进行数据更新。
[0056]
具体的，预设数据更新模型为：
[0057][0058]
当时，
[0059]
当时，y＝0。
[0060]
其中，y为目标高频信号分量。
[0061]
将目标数据通过预设数据扩展模型和预设数据更新模型后，得到相应的目标低频信号分量和目标高频信号分量。得到的相应的目标低频信号分量和目标高频信号分量都分别有固定的存放位置，必须按照相应的位置存放。
[0062]
其中存储原则为：目标数据si的目标低频信号分量存储于初始行向量的
处，目标数据si的目标高频信号分量y存储于初始行向量处，l1表示当前分解的级数，1≤l1≤l，l为上述步骤的心电信号的分解级数，l≤k。
[0063]
对每个心电信号分组中的所有数据全部进行数据扩展和数据更新后，得到第一低频信号分量和第一高频信号分量，其中，第一低频信号分量为所有目标低频信号分量的集合，第一高频信号分量为所有目标高频信号分量的集合。
[0064]
再次，根据预设数据更新模型，对第一低频信号分量中的所有数据进行数据更新，得到第二低频信号分量和第二高频信号分量。
[0065]
直至对第l-1低频信号分量和第l-1高频信号分量中的所有数据进行第l级数据数据更新后，得到第l低频信号分量和第l高频信号分量。对第l低频信号分量中的所有信号分量按照预设规则进行变换，得到第l低频信号变换分量，将第l低频信号变换分量确定为低频心电信号分量。具体的，经过数据扩展和数据更新后得到的第l低频信号分量中的所有信号分量都是区间值，需要对每个区间值变换，才能得到具体的心电数据。对于第l低频信号分量中的任意一个信号分量，由于信号分量是一个区间的值，需要对区间的值进行变换，可以采用将信号分量的区间2个端点的平均值作为所需的信号数据，并将得到的平均值存储在相应的位置，对第l低频信号分量中的所有信号分量均采用上述方法变换及存储后，得到第l低频信号变换分量。
[0066]
基于对心电信号的l级分解，得到低频心电信号分量c
low_l
和高频心电信号分量c
high_l
；其中，低频心电信号分量c
low_l
为第l低频信号分量，高频心电信号分量c
high_l
为第l高频信号分量至第一高频信号分量的集合。
[0067]
在本发明实施例中，通过设置预设噪声值，通过对心电信号的数据多次进行数据扩展和多次数据更新，当时，y＝0，即当前级分解得到的目标高频信号分量为0，通过l级分解，从而得到了去除小于预设噪声值的高频心电信号分量。
[0068]
例如：给定的心电信号的初始行向量中有8192个数据，如果对心电信号进行8次分解，则得到的低频心电信号分量c
low_l
的数据个数为8192/28，即32个低频心电信号分量。得到的高频心电信号分量c
high_l
的个数为8160个。其中，8160个高频心电信号分量c
high_l
中包含了数值为0的高频心电信号分量。此时得到的高频分量分别由第一级高频分量、第二级高频分量、第三级高频分量、第四级高频分量、第五级高频分量、第六级高频分量、第七级高频分量、第八级高频分量组成。一至八级高频分量数据量分别为：4096、2048、1024、512、256、128、64、32。
[0069]
根据上述内容，进行预处理操作可以根据获取到的不同心电信号的数据，从而确定不同的预设噪声值，进而通过对心电信号的数据的分解和扩展的处理，即可同时去除心电信号中的基线漂移和高频噪声，缩短计算时间。
[0070]
在一种可能的实现方式中，获取高频心电信号分量的待检测高频数据序列的极大值序列，包括：
[0071]
根据高频心电信号分量，获取高频心电信号分量中的预设级高频心电信号分量，并将预设级高频心电信号分量作为待检测高频数据序列；
[0072]
获取待检测高频数据序列的极大值；
[0073]
从待检测高频序列中筛选相邻两个极大值的位置距离大于预设距离的极大值构
成极大值序列。
[0074]
在本发明实施例中，对于获取的高频心电信号分量的待检测高频数据序列，预设级高频心电信号分量定义为第l级高频心电信号分量，其中，l的取值范围为1≤l≤l/2。获取待检测高频数据序列的过程包括：
[0075]
取出高频心电信号分量中的第l级(1≤l≤l/2高频心电信号分量h
l
作为待检测高频数据序列，其中，该级分量为原始数据长度的因此，在该数据上进行检测，可以降低检测时间。
[0076]
例如：给定的心电信号的初始行向量中有8192个数据，如果对心电信号进行8次分解，可以得到由低频分量和高频分量的组合分量w。此时得到的高频分量分别由第一级高频分量、第二级高频分量、第三级高频分量、第四级高频分量、第五级高频分量、第六级高频分量、第七级高频分量、第八级高频分量组成。一至八级高频分量数据量分别为：4096、2048、1024、512、256、128、64、32。根据数据存储规则，可以在组合分量w中找出每一级高频分量和低频分量的起始和结束位置。如第一级高频分量的起始和结束位置位于组合分量w的第4097个位置处和第8192个位置处，第二级高频分量的起始和结束位置位于组合分量w的第2049个位置处和第4096个位置处。第三级高频分量的起始和结束位置位于组合分量w的第1025个位置处和第2048个位置处。第四级高频分量的起始和结束位置位于组合分量w的第513个位置处和第1024个位置处。第五级高频分量的起始和结束位置位于组合分量w的第257个位置处和第512个位置处。第六级高频分量的起始和结束位置位于组合分量w的第129个位置处和第256个位置处。第七级高频分量的起始和结束位置位于组合分量w的第65个位置处和第128个位置处。第八级高频分量的起始和结束位置位于组合分量w的第33个位置处和第64个位置处。第八级低频分量的起始和结束位置位于组合分量w的第1个位置处和第32个位置处。
[0077]
对于从待检测高频序列中筛选相邻两个极大值的位置距离大于预设距离的极大值构成极大值序列的实现过程包括：
[0078]
设定相邻两个极大值的最小峰值间的位置间隔数为预设间隔g，在待检测高频数据序列上找到位置间隔大于g的极大值，找到的极大值构成极大值序列，由ma表示极大值序列。
[0079]
其中，预设间隔的表达式为：
[0080][0081]
其中，g为预设间隔，f为心电信号的采样频率，l为预设级，β为系数。优选的，l＝1，β＝0.15。
[0082]
在本发明实施例中，获取高频心电信号分量的待检测高频数据序列的极小值序列的实现方法与获取极大值序列的方法一致，具体如下：
[0083]
获取待检测高频数据序列的极小值，从待检测高频数据序列中筛选相邻两个极小值的位置间隔大于预设间隔的极小值构成极小值序列，其中，由mi表示极大值序列。
[0084]
举例说明，参照图2，图中所有的极大值由波峰
①
、
②
、
③
、
④
、
⑤
、
⑥
、
⑦
、
⑧
、
⑨
组成，坐标位置分别如图所示。所有极小值由波谷ⅰ、ⅱ、ⅲ、ⅳ、
ⅴ
组成，坐标位置分别如图所
示。由l＝1，β＝0.15，f＝360,得g＝27。
[0085]
以极大值为例，由上述论述可知，对于波峰
①
，波峰
②
距离波峰
①
大于g＝27，故波峰
①
需要保留。然后对于波峰
②
，在其后的波峰
③
在范围g＝27之内，由于波峰
③
峰值大于波峰
②
，故应保留波峰
③
而删除波峰
②
，以此保证h＝27范围内只能存在一个最大的波峰。同理可以得到波峰
④
、
⑤
、
⑧
需要被保留。由此预设规则可以获得极大值序列ma，ma由波峰
①
、
③
、
④
、
⑤
、
⑧
组成。
[0086]
同理可以得到极小值序列mi，mi由波谷ⅰ、ⅱ、ⅲ、ⅳ、
ⅴ
组成。
[0087]
在s102中，根据极大值序列计算极大值序列阈值，根据极小值序列计算极小值序列阈值，并根据极大值序列阈值和极大值序列的峰值，获得真极大值序列，根据极小值序列阈值和极小值序列的峰值，获得真极小值序列。
[0088]
在一种可能的实现方式中，根据极大值序列计算极大值序列阈值，包括：
[0089]
去除极大值序列中峰值小于等于零的极大值，并对处理后的极大值序列按照峰值从小到大的顺序进行排序，得到第一极大值序列，去掉第一极大值序列中后n个极大值，得到第二极大值序列；
[0090]
对第二极大值序列中后n个极大值求平均，得到第一信号值，对第二极大值序列中前m个极大值求平均，得到第二信号值；
[0091]
根据第一信号值和第二信号值，计算得到极大值序列阈值。
[0092]
具体的，根据s102中得到的极大值序列ma，首先去除ma中峰值小于等于0的极大值，然后按照峰值从小到大的顺序对处理后的极大值序列进行排序，得到第一极大值序列，再去掉第一极大值序列中的峰值最大的n个极大值，得到第二极大值序列，记为ma′
，每一个极大值ma′
＝(xai,yai)，其中，xai为该极大值的位置，yai为该极大值的峰值。这样做的目的是去除一些离散的极大值。
[0093]
取第二极大值序列ma′
中峰值最大的前n个极大值求平均，得到第一信号值，记为sig，取第二极大值序列ma′
中峰值最小的前m个极大值求平均，得到第二信号值，记为noi，其中，第二信号值也为噪声信号值。根据上述求取的第一信号值和第二信号值，计算极大值序列阈值。
[0094]
在本发明实施例中，可以将n设为10，可以将m设为5。
[0095]
其中，极大值序列阈值表达式为：
[0096]
thr1＝(sig-noi)
×
α noi
[0097]
其中，thr1为极大值序列阈值，sig为第一信号值，noi为第二信号值，α为系数，优选的，α＝0.05。
[0098]
在本发明实施例中，根据极小值序列计算极小值序列阈值的实现方式和计算极大值序列阈值的方法一致，具体如下：
[0099]
首先，将极小值序列中的极小值取绝对值，获得新极小值序列。
[0100]
去除新极小值序列中峰值小于等于零的新极小值，并对处理后的新极小值序列按照峰值从小到大的顺序进行排序，得到第一极小值序列，去掉第一极小值序列中后n个的极小值，得到第二极小值序列；对第二极小值序列中后n个极小值求平均，得到第三信号值，对第二极小值序列中前m个极小值求平均，得到第四信号值；根据第三信号值和第四信号值，
计算得到极小值序列阈值。
[0101]
在一种可能的实现方式中，根据极大值序列阈值和极大值序列的峰值，获得真极大值序列，包括：
[0102]
针对极大值序列中的任一个极大值，若该极大值的峰值大于第一阈值，则将该极大值作为初始真极大值，并将该初始真极大值存入初始真极大值序列；
[0103]
当初始真极大值序列包括至少三个初始真极大值时，将初始真极大值序列中前两位初始真极大值作为真极大值存入真极大值序列中；
[0104]
依次针对初始真极大值序列中第二位初始真极大值后的每个初始真极大值，计算该初始真极大值与第一真极大值之间的位置距离，并采用位置距离更新当前平均值，其中，初始真极大值序列中的初始真极大值按照进入初始真极大值序列的时间先后顺序排序，真极大值序列中的真极大值按照进入真极大值序列的时间先后顺序排序，第一真极大值为当前真极大值序列中最后一位真极大值；
[0105]
若该初始真极大值与第一真极大值之间的位置距离满足第一条件，则确定该初始真极大值为真极大值，并将该初始真极大值存入真极大值序列，其中，第一条件为：
[0106]
a*avg《d《b*avg；
[0107]
其中，avg为更新后的当前平均值，d为该初始真极大值与第一真极大值之间的位置距离，a为第一系数，b为第二系数，且a≠b≠0；
[0108]
若该初始真极大值与第一真极大值之间的位置距离满足第二条件，且该初始真极大值的峰值大于第一真极大值的峰值，且该初始真极大值与第二真极大值之间的位置距离满足第三条件，则确定该初始真极大值为真极大值，并去除第一真极大值，将该初始真极大值存入真极大值序列，其中，第二真极大值为第一真极大值前一顺位的真极大值，第二条件为：
[0109]
d《a*avg；
[0110]
第三条件为：
[0111]
d1《b*avg；
[0112]
其中，d1为该初始真极大值与第二真极大值之间的位置距离；
[0113]
若该初始真极大值与第一真极大值之间的位置距离满足第四条件，则在极大值序列中查找位于该初始真极大值与第一真极大值之间的极大值作为目标真极大值，若目标真极大值大于第二阈值，则确定目标真极大值为真极大值，将目标真极大值存入真极大值序列，其中，第二阈值为第一阈值的二分之一，第四条件为：
[0114]
d≥b*avg；
[0115]
若该初始真极大值与第一真极大值之间的位置距离满足第五条件，则确定该初始真极大值为真极大值，将该初始真极大值存入真极大值序列，其中，第五条件为：
[0116]
d》a*avg。
[0117]
其中，当前平均值为每次更新后的真极大值序列中所有真极大值与其前一顺位真极大值的位置距离之和的平均值。
[0118]
本发明实施例中，第一阈值为极大值序列阈值，则第二阈值的表达式为：
[0119]
[0120]
其中，thr2为第二阈值。
[0121]
具体的，如果只采用固定的阈值进行峰值筛选，有可能导致一部分峰值较小的极大值被漏检，或者有可能导致一部分峰值较大的噪声信号值被误认为为极大值。因此需要通过降低阈值，采用漏检机制进行补偿，同时采用位置距离d作为进一步判断极大值的依据。具体步骤如下：
[0122]
针对极大值序列ma中的任一个极大值，若该极大值大于第一阈值thr1，则将该极大值作为初始真极大值，存入初始真极大值序列m0，m0由位置序列和峰值序列组成，否则不存储。
[0123]
当初始真极大值序列m0中有3个及以上的初始真极大值时，将初始真极大值序列中前两位初始真极大值作为真极大值存入真极大值序列m中。
[0124]
依次针对初始真极大值序列中第二位初始真极大值后的每个初始真极大值，计算该初始真极大值与第一极大值之间的位置距离，并求出该位置距离的平均值avg。
[0125]
若a*avg《d《b*avg，则确定该初始真极大值为真极大值，将该初始真极大值存入极大值序列m，并更新当前平均值avg，否则不存储到m中；
[0126]
若d《a*avg，且该初始真极大值的峰值大于第一真极大值的峰值，且d1《b*avg，则确定该初始真极大值为真极大值，将该初始真极大值存入真极大值序列m，并更新当前平均值avg，否则不存储到m中，该步骤可以降低误检率；
[0127]
若d≥b*avg时，说明有漏检，则需要在极大值序列中查找位于该初始真极大值与第一真极大值之间的目标真极大值，若目标真极大值大于第二阈值，则确定目标真极大值为真极大值，将目标真极大值存入真极大值序列m，并更新当前平均值avg，否则不存储到m中；其中，若时，同样适用于该步骤；
[0128]
若d》a*avg，则确定该初始真极大值为真极大值，将该初始真极大值存入真极大值序列m，并更新当前平均值avg，否则不存储到m中，该步骤也可以降低漏检。
[0129]
在本发明实施例中，根据极小值序列阈值和极小值序列的峰值，获得真极小值序列的实现方式与真极大值序列的获取方式一致，但是，前提需要对极小值序列中的极小值取绝对值，获得新极小值序列，就不在详细赘述。
[0130]
举例说明，以极大值序列为例，参照图3，每一个极大值的坐标如图所示。首先根据阈值选取方法，可得到极大值序列ma的信号大小sig＝0.3030，噪声大小noi＝0.0283，由第一阈值公式thr1＝(sig-noi)
×
α noi，α＝0.05，得，thr1＝0.042，
[0131]
然后，根据获得初始真极大值序列，可知波峰
①
峰值大于第一阈值th1，说明当前波峰
①
是初始真极大值，需要存储到真极大值序列m中；然后判断波峰
②
，也是初始真极大值，需要将其存储到真极大值序列m中，此时m中已经有两个真极大值，因此需要计算出其位置距离的平均值，可得avg＝(226-60)/1＝166；然后判断波峰
③
：此时波峰
③
大于第一阈值thr1，为初始真极大值，其与m序列最后一个真极大值的位置距离d＝390-226＝164，满足条件0.4
×
avg《d《1.5
×
avg，因此判定当前波峰
③
是真极大值，需要将其存储到m序列中，此时m序列中有3个极大值，也就是两个位置距离差值，因此需要更新然后判断波峰
④
，由于其数值小于第一阈值thr1，因此波峰
④
不是初始真极大值，不将其存
储到m序列中，不更新avg；然后判断波峰
⑤
，其数值大于第一阈值thr1，为初始真极大值，并且其位置距离d＝695-390＝305，满足条件或者即满足第四条件，说明有漏检，需要在波峰
⑤
之前，m序列最后一个波峰
③
之后重新查找，此时需将阈值降低一半，即采用第二阈值thr2＝0.021进行检测，此时检测到波峰
④
，由于波峰
④
的数值大于第二阈值thr2，因此需将波峰
④
存储到m序列中，更新avg＝161.33，接着判断波峰
⑤
，其位置距离d＝695-544＝151，满足d》0.4
×
avg，因此波峰
⑤
为真极大值，将其存储到m序列中，更新avg＝158.75；然后判断波峰
⑥
，其数值小于第一阈值thr1，因此其不是初始真极大值，不保存到m序列中，不更新avg；然后判断波峰
⑦
，其满足0.4
×
avg《d《1.5
×
avg，因此波峰
⑦
为真极大值，将其存储到m序列中，更新avg＝156.4；然后判断波峰
⑧
，其满足0.4
×
avg《d《1.5
×
avg，因此波峰
⑧
为真极大值，将其存储到m序列中，更新avg＝149；然后判断波峰
⑨
，其阈值大于th1，为初始真极大值，其与m序列最后一个波峰
⑧
之间的位置距离d＝984-954＝30，满足d≤0.4
×
avg，且波峰
⑨
大于m序列最后的波峰
⑧
，且波峰
⑨
与m序列最后第2个波峰
⑦
的位置距离为142，其位置距离小于1.5
×
avg，因此满足d《0.4
×
avg，且d1《1.5
×
avg，因此需要用波峰
⑨
替换掉m序列中的最后一个波峰
⑧
，并更新avg＝154；然后判断波峰
⑩
，其数值大于第一阈值thr1，为初始真极大值，其位置距离d＝1042-984＝58，满足d≤0.4
×
avg，但是其数值小于m序列最后的波峰
⑨
，故当前波峰
⑩
不能替换波峰
⑨
，不用存储到m中。
[0132]
由以上步骤可以得到极大值序列ma的真极大值序列m，其由波峰
①
、
②
、
③
、
④
、
⑤
、
⑦
、
⑨
组成。同理也可以对极小值序列mi做同样处理得到其真极小值序列v。
[0133]
在s103中，对真极大值序列和真极小值序列进行极值对匹配，定位r波峰值范围。
[0134]
在一种可能的实现方式中，s103的具体实现流程包括：
[0135]
根据真极大值序列中极大值的位置和真极小值序列中极小值的位置，判断相邻的极大值和极小值之间的位置距离是否小于给定阈值；
[0136]
若相邻的极大值和极小值之间的位置距离小于给定阈值，则将相邻的极大值和极小值作为一对极值对；
[0137]
分别对每一对极值对进行斜率计算，若存在一对极值对的斜率大于预设阈值，则确定该极值对为真极值对；
[0138]
根据真极值对，定位r波峰值范围。
[0139]
在本发明实施例中，定位r波峰值范围包括两个步骤：
[0140]
第一步骤，根据真极大值序列和真极小值序列，若相邻的极大值和极小值之间的距离小于给定阈值thr3，则确定为一对可匹配的极值对，这样可以去掉孤立的极值点。据此可获得匹配的极值对组合[p
l
,pr]，p
l
＝[p
l
(1),p
l
(2),
…
,p
l
(i),
…
,p
l
(k)]，pr＝[pr(1),pr(2),
…
,pr(i),
…
,pr(k)]，每一个极值对p
l
(i)＝(loc
l
(i),pks
l
(i)),pr(i)＝(locr(i),pksr(i))，其中坐标(loc
l
(i),pks
l
(i))和(locr(i),pksr(i))分别代表第i个极值对的极大值和极小值的坐标，loc
l
(i)和locr(i)是横坐标，代表峰值的位置，pks
l
(i)和pksr(i)是纵坐标，代表峰值的大小，k为极值对的个数。
[0141]
其中，给定阈值的表达式为：
[0142]
[0143]
其中，thr3为给定阈值。
[0144]
举例说明，参照图4，真极大值序列m(真极大值点分别为
①
、
②
、
③
)和真极小值序列v(真极小值点分别为ⅰ、ⅱ)构成的极值序列图。
[0145]
极小值ⅰ和极大值
①
，并且二者之间的位置距离为8，小于给定阈值因此极小值ⅰ和极大值
①
可以匹配，由此可以获得极值对组合[p
l
,pr]，p
l
对应极小值ⅰ，pr对应极大值
①
；对于极大值
①
和极大值
②
，因为二者都是极大值，因此不能匹配。对于极大值
②
和极小值ⅱ，二者分别是极大值和极小值，但是二者的位置距离为49，其大于给定阈值thr3，二者不能匹配，因此可以去除孤立的极大值
②
；对于极小值ⅱ和极大值
③
，二者分别是极小值和极大值，并且二者之间的位置距离为8，小于给定阈值thr3，可以匹配，由此可以获得极值对组合[p
l
,pr]，p
l
对应极小值ⅱ，pr对应极大值
③
。
[0146]
第二步，以斜率作为判据，去掉伪极值对。其中，每一个极值对的斜率的计算公式为：
[0147][0148]
其中，slope(i)为每一个极值对的斜率。
[0149]
对于前a个极值对，预设阈值设置为第一预设阈值，对于剩下的极值对，预设阈值设置为第二预设阈值，则判定该极值对为真极值对，保留真极值对组合[p
′
l
,p
′r]。
[0150]
例如，a＝3，则第一预设阈值即等于前三个极值对的斜率平均值的二分之一，第二预设阈值即等于前三个极值对的斜率平均值的二分之一，第二预设阈值即等于真极值对组合[p
′
l
,p
′r]的最新的三个极值对斜率平均值的三分之一。
[0151]
根据上述步骤可得到最终的极值对组合[p
′
l
,p
′r]，p
′
l
＝[p
′
l
(1),p
′
l
(2),
…
,p
′
l
(i),
…
,p
′
l
(g)]，p
′r＝[p
′r(1),p
′r(2),
…
,p
′r(i),
…
,p
′r(g)]，每一个极值对p
′
l
(i)＝(loc
′
l
(i),pks
′
l
(i))，p
′r(i)＝(loc
′r(i),pks
′r(i))，其中坐标(loc
′
l
(i),pks
′
l
(i))和(loc
′r(i),pks
′r(i))分别代表第i个极值对的极大值和极小值的坐标，loc
′
l
(i)和loc
′r(i)是横坐标，代表峰值的位置，pks
′
l
(i)和pks
′r(i)是纵坐标，代表峰值的大小，g为极值对的个数。
[0152]
举例说明，参照图5，图为根据位置距离匹配后获得的极值对，图中极大值
①
和极小值ⅰ匹配、极大值
②
和极小值ⅱ匹配、极大值
③
和极小值ⅲ匹配，极大值
④
和极小值ⅳ匹配、极大值
⑤
和极小值
ⅴ
匹配、极大值
⑥
和极小值ⅵ匹配，极大值
⑦
和极小值ⅶ匹配。
[0153]
接下来需要根据斜率去掉伪极值对。根据斜率计算公式可求得初始阈值：
[0154][0155]
[0156][0157][0158]
逐一判断极值对是否是真极值对。首先判断第一对极值对：极大值
①
和极小值ⅰ，其斜率为0.009425，其数值大于第一预设阈值thr4，故第一对极值对是真极值对，要存储到真极值对[p
′
l
,p
′r]中；接着判断第二对极值对：极大值
②
和极小值ⅱ，其斜率为0.008154，其数值大于第一预设阈值thr4，故第二对极值对是真极值对，要存储到真极值对[p
′
l
,p
′r]中；接着判断第三对极值对：极大值
③
和极小值ⅲ，其斜率为0.009717，其数值大于第一预设阈值thr4，故第三对极值对是真极值对，要存储到真极值对[p
′
l
,p
′r]中；接着判断第四对极值对：极大值
④
和极小值ⅳ，其斜率为0.011492，此时真极值[p
′
l
,p
′r]中已经有三个极值对，因此需要根据最新的三个极值对计算出第二预设阈值对，因此需要根据最新的三个极值对计算出第二预设阈值此时斜率大于第二预设阈值thr5，故第四对极值对是真极值对，要存储到真极值对[p
′
l
,p
′r]中；接着判断第五对极值对：极大值
⑤
和极小值
ⅴ
，其斜率为0.02404，此时第二预设阈值thr5＝0.003263，此时斜率大于第二预设阈值thr5，故第五对极值对是真极值对，要存储到真极值对[p
′
l
,p
′r]中；接着判断第六对极值对：极大值
⑥
和极小值ⅵ，其斜率为0.002828，此时第二预设阈值thr5＝0.004453，此时斜率小于第二预设阈值thr5，故第六对极值对不是真极值对，不需要存储到真极值对[p
′
l
,p
′r]中；接着判断第七对极值对：极大值
⑦
和极小值ⅶ，其斜率为0.006483，此时第二预设阈值thr5＝0.004453，此时斜率大于第二预设阈值thr5，故第七对极值对是真极值对，需要存储到真极值对[p
′
l
,p
′r]中。故经过此步骤可以去除伪极值对，如极大值
⑥
和极小值ⅵ，得到真极值对：极大值
①
和极小值ⅰ、极大值
②
和极小值ⅱ、极大值
③
和极小值ⅲ，极大值
④
和极小值ⅳ、极大值
⑤
和极小值
ⅴ
、极大值
⑦
和极小值ⅶ。
[0159]
在s104中，计算r波峰值范围内每个极值对的极大值和极小值之间的过零点，组成过零点集合，并将过零点集合映射到心电信号，定位r波峰值位置。
[0160]
在一种可能的实现方式中，s104的具体实现流程包括：
[0161]
针对r波峰值范围内的每个极值对，基于该极值对的极大值和极小值进行直线拟合，得到直线拟合公式，并基于直线拟合公式确定过零点；将该极值对的过零点存入过零点集合；
[0162]
将过零点集合中过零点按照映射方法映射到心电信号，定位r波峰值位置；
[0163]
其中，映射方法的计算公式为：
[0164][0165]
其中，pk为r波峰值位置，round()为对r波峰值位置进行四舍五入的函数，pos为过零点集合中的过零点，l为高频心电信号分量中的第l级高频心电信号分量。
[0166]
在本发明实施例中，可以根据每个极值对的极大值和极小值的位置和大小，计算极大值和极小值之间的过零点。在这里将极大值和极小值按照直线相连接，通过直线拟合公式可以计算出直线的过零点，将计算的过零点组成过零点集合pos，其中，pos＝{pos(1),
pos(2),
…
,pos(g)}。然后将计算后的过零点集合映射到心电信号中，可以获得心电信号的r波峰值位置。
[0167]
其中，参照图6所示，由极大值和极小值组成的直线，其过零点o的横坐标x的计算公式如下：由推导得：
[0168]
例如，参照图5，由极大值
①
和极小值ⅰ组成的极值对，其坐标分别为(35,0.0707)和(47,-0.0424)，由上述公式可以求得过零点位置为36.13，映射到心电信号中r波峰值位置：此时假设l＝1，那么置：此时假设l＝1，那么
[0169]
由此可以求出每一个极值对连线的过零点，从而定位每一个r波峰值的位置。
[0170]
本发明提供一种心电信号r波定位方法，与现有技术相比，本发明所采用的高频心电信号分量不需要进行数据扩展，可以降低计算量，降低检测时间，同时，通过获得真极大值序列和真极小值序列，并进行极值对匹配，粗定位r波峰值位置，再通过计算过零点，映射到心电信号上，精准定位r波峰值位置，进一步提高r波峰值定位的正确度。
[0171]
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。
[0172]
以下为本发明的装置实施例，对于其中未详尽描述的细节，可以参考上述对应的方法实施例。
[0173]
图7示出了本发明实施例提供的心电信号r波定位装置的结构示意图，为了便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0174]
如图7所示，心电信号r波定位装置7包括：
[0175]
获取模块71，用于获取待识别的心电信号中的高频心电信号分量的待检测高频数据序列的极大值序列和极小值序列，其中，极大值序列包括多个极大值，极小值序列包括多个极小值，极大值和极小值均包括位置和峰值；
[0176]
计算模块72，用于根据极大值序列计算极大值序列阈值，根据极小值序列计算极小值序列阈值，并根据极大值序列阈值和极大值序列的峰值，获得真极大值序列，根据极小值序列阈值和极小值序列的峰值，获得真极小值序列；
[0177]
第一定位模块73，用于对真极大值序列和真极小值序列进行极值对匹配，定位r波峰值范围；
[0178]
第二定位模块74，用于计算r波峰值范围内每个极值对的极大值和极小值之间的过零点，组成零点集合，并将零点集合映射到心电信号，定位r波峰值位置。
[0179]
本发明提供一种心电信号r波定位装置，与现有技术相比，本发明所采用的高频心电信号分量不需要进行数据扩展，可以降低计算量，降低检测时间，同时，通过获得真极大值序列和真极小值序列，并进行极值对匹配，粗定位r波峰值位置，再通过计算过零点，映射到心电信号上，精准定位r波峰值位置，进一步提高r波峰值定位的正确度。
[0180]
在一种可能的实现方式中，获取模块可以用于：
[0181]
根据预设极限漂移频率和心电信号的采样频率，确定心电信号的分解级数；
[0182]
将心电信号分为n个心电信号分组，获取每个心电信号分组的预设噪声值，其中，n
为大于等于2的正整数；
[0183]
根据分解级数和每个心电信号分组的预设噪声值，对各个心电信号分组的心电信号进行分解，得到高频心电信号分量；
[0184]
获取高频心电信号分量的待检测高频数据序列的极大值序列和极小值序列。
[0185]
在一种可能的实现方式中，获取模块还可以用于：
[0186]
根据高频心电信号分量，获取高频心电信号分量中的预设级高频心电信号分量，并将预设级高频心电信号分量作为待检测高频数据序列；
[0187]
获取待检测高频数据序列的极大值；
[0188]
从待检测高频序列中筛选相邻两个极大值的位置距离大于预设距离的极大值构成极大值序列。
[0189]
在一种可能的实现方式中，计算模块可以用于：
[0190]
去除极大值序列中峰值小于等于零的极大值，并对处理后的极大值序列按照峰值从小到大的顺序进行排序，得到第一极大值序列，去掉第一极大值序列中后n个极大值，得到第二极大值序列；
[0191]
对第二极大值序列中后n个极大值求平均，得到第一信号值，对第二极大值序列中前m个极大值求平均，得到第二信号值；
[0192]
根据第一信号值和第二信号值，计算得到极大值序列阈值。
[0193]
在一种可能的实现方式中，计算模块还可以用于：
[0194]
针对极大值序列中的任一个极大值，若该极大值的峰值大于第一阈值，则将该极大值作为初始真极大值，并将该初始真极大值存入初始真极大值序列；
[0195]
当初始真极大值序列包括至少三个初始真极大值时，将初始真极大值序列中前两位初始真极大值作为真极大值存入真极大值序列中；
[0196]
依次针对初始真极大值序列中第二位初始真极大值后的每个初始真极大值，计算该初始真极大值与第一真极大值之间的位置距离，并采用位置距离更新当前平均值，其中，初始真极大值序列中的初始真极大值按照进入初始真极大值序列的时间先后顺序排序，真极大值序列中的真极大值按照进入真极大值序列的时间先后顺序排序，第一真极大值为当前真极大值序列中最后一位真极大值；
[0197]
若该初始真极大值与第一真极大值之间的位置距离满足第一条件，则确定该初始真极大值为真极大值，并将该初始真极大值存入真极大值序列，其中，第一条件为：
[0198]
a*avg《d《b*avg；
[0199]
其中，avg为更新后的当前平均值，d为该初始真极大值与第一真极大值之间的位置距离，a为第一系数，b为第二系数，且a≠b≠0；
[0200]
若该初始真极大值与第一真极大值之间的位置距离满足第二条件，且该初始真极大值的峰值大于第一真极大值的峰值，且该初始真极大值与第二真极大值之间的位置距离满足第三条件，则确定该初始真极大值为真极大值，并去除第一真极大值，将该初始真极大值存入真极大值序列，其中，第二真极大值为第一真极大值前一顺位的真极大值，第二条件为：
[0201]
d《a*avg；
[0202]
第三条件为：
[0203]
d1《b*avg；
[0204]
其中，d1为该初始真极大值与第二真极大值之间的位置距离；
[0205]
若该初始真极大值与第一真极大值之间的位置距离满足第四条件，则在极大值序列中查找位于该初始真极大值与第一真极大值之间的极大值作为目标真极大值，若目标真极大值大于第二阈值，则确定目标真极大值为真极大值，将目标真极大值存入真极大值序列，其中，第二阈值为第一阈值的二分之一，第四条件为：
[0206]
d≥b*avg；
[0207]
若该初始真极大值与第一真极大值之间的位置距离满足第五条件，则确定该初始真极大值为真极大值，将该初始真极大值存入真极大值序列，其中，第五条件为：
[0208]
d》a*avg。
[0209]
在一种可能的实现方式中，第一定位模块具体用于：
[0210]
根据真极大值序列中极大值的位置和真极小值序列中极小值的位置，判断相邻的极大值和极小值之间的位置距离是否小于给定阈值；
[0211]
若相邻的极大值和极小值之间的位置距离小于给定阈值，则将相邻的极大值和极小值作为一对极值对；
[0212]
分别对每一对极值对进行斜率计算，若存在一对极值对的斜率大于预设阈值，则确定该极值对为真极值对；
[0213]
根据真极值对，定位r波峰值范围。
[0214]
在一种可能的实现方式中，第二定位模块具体用于：
[0215]
针对r波峰值范围内的每个极值对，基于该极值对的极大值和极小值进行直线拟合，得到直线拟合公式，并基于直线拟合公式确定过零点；将该极值对的过零点存入过零点集合；
[0216]
将过零点集合中过零点按照映射方法映射到心电信号，定位r波峰值位置；其中，映射方法的计算公式为：
[0217][0218]
其中，pk为r波峰值位置，round()为对r波峰值位置进行四舍五入的函数，pos为过零点集合中的过零点，l为高频心电信号分量中的第l级高频心电信号分量。
[0219]
图8是本发明实施例提供的终端的示意图。如图8所示，该实施例的终端8包括：处理器80、存储器81以及存储在所述存储器81中并可在所述处理器80上运行的计算机程序82。所述处理器80执行所述计算机程序82时实现上述各个心电信号r波定位方法实施例中的步骤，例如图1所示的s101至s104。或者，所述处理器80执行所述计算机程序82时实现上述各装置实施例中各模块/单元的功能，例如图7所示模块71至74的功能。
[0220]
示例性的，所述计算机程序82可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器81中，并由所述处理器80执行，以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序82在所述终端8中的执行过程。例如，所述计算机程序82可以被分割成图7所示的模块71至74。
[0221]
所述终端8可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所
述终端8可包括，但不仅限于，处理器80、存储器81。本领域技术人员可以理解，图8仅仅是终端8的示例，并不构成对终端8的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述终端还可以包括输入输出设备、网络接入设备、总线等。
[0222]
所称处理器80可以是中央处理单元(central processing unit，cpu)，还可以是其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
[0223]
所述存储器81可以是所述终端8的内部存储单元，例如终端8的硬盘或内存。所述存储器81也可以是所述终端8的外部存储设备，例如所述终端8上配备的插接式硬盘，智能存储卡(smart media card，smc)，安全数字(secure digital，sd)卡，闪存卡(flash card)等。进一步地，所述存储器81还可以既包括所述终端8的内部存储单元也包括外部存储设备。所述存储器81用于存储所述计算机程序以及所述终端所需的其他程序和数据。所述存储器81还可以用于暂时地存储已经输出或者将要输出的数据。
[0224]
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0225]
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
[0226]
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
[0227]
在本发明所提供的实施例中，应该理解到，所揭露的装置/终端和方法，可以通过其它的方式实现。例如，以上所描述的装置/终端实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。
[0228]
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目
的。
[0229]
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0230]
所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个心电信号r波定位方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、u盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括是电载波信号和电信信号。
[0231]
以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。