买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：河北大学

摘要：本发明提供了一种心电信号心率变异性特征提取方法。该方法首先获取原始心电信号，采用Pan‑Tompkins方法对心电信号进行R特征波检测并提取HRV信号，在HRV信号中顺序选取N个采样点构成一组HRV时间序列片段，然后对HRV时间序列片段进行50尺度庞加莱散点图分析，得到HRV时间序列片段在50个尺度上的庞加莱分布散点图，并提取每一尺度的特征参数SD1，然后提取HRV时间序列片段的阈值特征w和尺度变化特征Sv；通过w或Sv识别健康个体和心脏疾病个体。本发明通过多尺度庞加莱散点图分析建立了HRV时间序列片段的多尺度特征，通过研究不同尺度下的特征变化特性，采用w和Sv均可准确对心电信号心率变异性特征进行提取。

主权项：1.一种心电信号心率变异性特征提取方法，其特征是，包括如下步骤：S1、获取原始心电信号，采用Pan-Tompkins方法对原始心电信号进行R特征波检测，提取心率变异性信号RRm，1≤m≤K，K为正整数；S2、在心率变异性信号RRm中顺序选取N个采样点构成一组HRV时间序列片段{RRi：1≤i≤N}；N≤K；S3、对步骤S2中所形成的HRV时间序列片段{RRi：1≤i≤N}进行50尺度庞加莱散点图分析，得到HRV时间序列片段在50个尺度上的庞加莱散点分布图，并提取每一尺度的特征参数SD1s，1≤s≤50；S4、提取HRV时间序列片段的阈值特征w，阈值特征w的计算公式为： w=Z－y1其中，Z为第sx尺度下的SD1值，y1为阈值；阈值y1的计算公式为： 其中，H和D分别表示第sx尺度下健康个体和心脏疾病个体SD1的均值； s x 的求解过程如下：①选取若干健康个体和若干心脏疾病个体的原始心电信号，并依据步骤S1获得心率变异性信号，依据步骤S2构建HRV时间序列片段数据集，依据步骤S3对HRV时间序列片段数据集中的每一组HRV时间序列片段均在50个尺度上提取特征参数SD1s；②分别绘制健康个体和心脏疾病个体HRV时间序列片段各尺度下SD1的误差棒图；③对健康个体和心脏疾病个体的误差棒图进行多尺度分析，寻找使得健康个体和心脏疾病个体误差棒图重合部分最小的尺度因子，该尺度因子即为sx；上述H和D即可在求解sx的过程中获得；S5：若阈值特征w大于0，则判定该数据为健康个体对应的数据，若小于0，则该数据为心脏疾病个体对应的数据；步骤S3具体如下：首先，对HRV时间序列片段进行粗粒化计算，各尺度的粗粒化时间序列计算公式为： 其中，s为尺度因子，1≤s≤50；则每个尺度下的时间序列长度等于原始时间序列长度N除以尺度因子s，即Ns=Ns；则第s尺度下的时间序列为{RR1，RR2，RR3，……，RRj，……，RRNs}；接着，依次求出1-50个尺度下时间序列的庞加莱散点图的短轴SD1，即{SD11，SD12，……，SD150}；各尺度下的庞加莱散点图的短轴SD1计算公式为： 其中，，为△RR的均值，Ns为第s尺度时间序列的长度，即Ns=Ns。

全文数据：

权利要求：

百度查询： 河北大学 一种心电信号心率变异性特征提取方法