段数据的均方差和峰峰值，若某段均方差和峰峰值大于方差阈值和峰峰值阈值，则将该段数据判断为不合理的数据段；3)将不合理的数据段用合理数据进行替换，合理数据的选择使得替换后的数据段与它的前一段数据段之间没有跳变，同时调整后面所有数据段确保没有跳变，依此方法获得重新构建的数据；4)对重新构建的数据做快速傅立叶变换，取模得到频域空间的频谱值数据；5)对于频域空间的频谱值数据，频率小于最低频率值的频谱值全部取0，其中，所述最低频率值为人体最低心率对应的频率值；6)搜索频谱峰，按照峰的高低降序排列，排列后的峰，从第一个峰开始向后搜索，出现在该峰左右最低频率值内的峰直接删除；7)保留最高的三个峰，并修正每个峰的位置；8)对于最高三个峰的频率位置，判断其最优频率位置关系，如果最优频率位置关系判据小于判断阈值，则转到步骤a)对心率原始数据做中心差分获得新的原始心率数据，其中，新的原始心率数据中第η个数据为旧的原始数据中第n+1个数据减去第η-1个数据，并且，将新的原始心率数据作为后续步骤的原始心率数据；其中，步骤a)执行在步骤I)与步骤2)之间；9)按照最优频率位置关系做频率的直线拟合，并计算拟合后的各个峰的频率值，若拟合优度低于拟合阈值，或者拟合直线的截距绝对值大于截距阈值，则转到步骤a)对心率原始数据做中心差分获得新的原始心率数据，其中，新的原始心率数据中第η个数据为旧的原始数据中第n+1个数据减去第n-1个数据，并且，将新的原始心率数据作为后续步骤的原始心率数据；其中，步骤a)执行在步骤I)与步骤2)之间；10)由采用的频率位置关系拟合得到的三个峰所在的频率，分别计算拟合数据，取均值作为最终心率结果。本发明用统计方法判断时域数据的取舍；频域数据，低于最小可能心率的频率舍弃，各个主峰的最小可能心率邻域内的峰舍弃；采用的频谱峰按照可能的几种倍频关系逐一判断决定；在原始数据中倍频分量较小的情况，对原始数据差分以突出倍频；最后直线拟合的优劣和直线截距决定心率结果是否采用。本发明可以获得比较精确的心率值，在心率测量领域具有广泛及深远的应用前景。所以，本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
[0078]上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
【主权项】
1.一种心率测量中的信号处理方法，其特征在于，所述信号处理方法包括步骤: 1)以一定采样频率采集一段时间的心率原始数据； 2)将心率原始数据分段，计算每段数据的均方差和峰峰值，若某段均方差和峰峰值大于方差阈值和峰峰值阈值，则将该段数据判断为不合理的数据段； 3)将不合理的数据段用合理数据进行替换，合理数据的选择使得替换后的数据段与它的前一段数据段之间没有跳变，同时调整后面所有数据段确保没有跳变，依此方法获得重新构建的数据； 4)对重新构建的数据做快速傅立叶变换，取模得到频域空间的频谱值数据； 5)对于频域空间的频谱值数据，频率小于最低频率值的频谱值全部取O，其中，所述最低频率值为人体最低心率对应的频率值； 6)搜索频谱峰，按照峰的高低降序排列，排列后的峰，从第一个峰开始向后搜索，出现在该峰左右最低频率值内的峰直接删除； 7)保留最高的三个峰，并修正每个峰的位置； 8)对于最高三个峰的频率位置，判断其最优频率位置关系，如果最优频率位置关系判据小于判断阈值，则退出计算，否则继续计算； 9)按照最优频率位置关系做频率的直线拟合，并计算拟合后的各个峰的频率值，拟合优度低于拟合阈值，或者拟合直线的截距绝对值大于截距阈值，则退出计算，否则继续计算; 10)由采用的频率位置关系拟合得到的三个峰所在的频率，分别计算拟合数据，取均值作为最终心率结果。2.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤8)中，如果最优频率位置关系判据小于判断阈值，则转到步骤a)对心率原始数据做中心差分获得新的原始心率数据，其中，新的原始心率数据中第η个数据为旧的原始数据中第n+l个数据减去第η-1个数据，并且，将新的原始心率数据作为后续步骤的原始心率数据；其中，步骤a)执行在步骤I)与步骤2)之间 若执行过步骤a)?步骤8)之后，最优频率位置关系判据依然小于判断阈值，则退出计算，否则继续计算。3.根据权利要求1或2所述的心率测量中的信号处理方法，其特征在于:步骤9)中，若拟合优度低于拟合阈值，或者拟合直线的截距绝对值大于截距阈值，则转到步骤a)对心率原始数据做中心差分获得新的原始心率数据，其中，新的原始心率数据中第η个数据为旧的原始数据中第n+l个数据减去第η-1个数据，并且，将新的原始心率数据作为后续步骤的原始心率数据；其中，步骤a)执行在步骤I)与步骤2)之间。 若执行过步骤a)?步骤9)之后，拟合优度依然低于拟合阈值，或者拟合直线的截距绝对值依然大于截距阈值，则退出计算，否则继续计算。4.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤I)中，最低采样频率满足:需要采用的最大频率是人最大可能心率对应的频率的三倍，按照采样定理决定的采样频率至少为该最大频率的两倍。5.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤2)中，心率原始数据分段依据是:以人最小可能心率考虑，每段数据都包含至少一个完整心率周期。6.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤2)中，方差阈值和峰峰值阈值可以设定一个固定值，也可以在程序运行中实时产生，或者两种方法的综合，包括: 最初设定一个比较大的值，以后每次计算出的判断认为可以接受的均方差分别乘系数后与原有的阈值加权平均后分别成为新的方差阈值及峰峰值阈值；或者 每次数据分段后，从各段中取最小的均方差分别乘系数后作为阈值，同时，该阈值不能超过一个数值，即最大的阈值，超过该最大阈值，则所有数据段都是不能接受的，无法将进行后继计算。7.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤2)中，方差阈值和峰峰值阈值的确定方法为:每次数据分段后，从各段中选出最小均方差的段，以它的均方差的三倍为方差阈值，九倍为峰峰值阈值。8.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤3)中，将不合理的数据段进行替换包括: 直接用一个固定值替换不合理数据段，替换后的直线段与前后数据段没有跳变；或者 用合理的数据段重复填充错误数据段，替换后的数据段与前后数据段没有跳变，如果所有数据段都是错误的，则停止计算，等待下一批数据；如果所有数据段合理则数据不变。9.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤5)中，所述人体最低心率为不小于40次/分钟，所述最低频率值为不大于0.67Hz。10.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤7)中，采用抛物线拟合、高斯函数形状拟合、或柯西-洛仑兹函数形状拟合的方法修正每个峰的位置。11.根据权利要求1所述的心率测量中的信号处理方法，其特征在于:步骤8)中，对于最高三个峰的频率位置，判断其位置关系包括:1:2:3、1:3:4、1:2:4，优先以上次计算的位置关系判断，第一次默认判断1:2: 3，判断后取最优，作为最优频率位置关系。
【专利摘要】本发明提供一种心率测量中的信号处理方法，所述信号处理方法包括步骤：用统计方法判断时域数据的取舍；频域数据，低于最小可能心率的频率舍弃，各个主峰的最小可能心率邻域内的峰舍弃；采用的频谱峰按照可能的几种倍频关系逐一判断决定；在原始数据中倍频分量较小的情况，对原始数据差分以突出倍频；最后直线拟合的优劣和直线截距决定心率结果是否采用。本发明可以获得比较精确的心率值，在心率测量领域具有广泛及深远的应用前景。
【IPC分类】A61B5/0245
【公开号】CN104887215
【申请号】CN201510323890
【发明人】黄建华, 李宏, 彭铁刚, 张廷禄
【申请人】上海新微技术研发中心有限公司
【公开日】2015年9月9日
【申请日】2015年6月12日