[发明专利]一种抗运动干扰的脉率提取方法

权利要求书

1.一种抗运动干扰的脉率提取方法，其特征在于，包括如下步骤：

步骤一、采集脉搏信号并获取其交流成分后转换为频谱信号；

步骤二、从频谱信号的峰值中搜索具有近似二倍关系的频率对；

步骤三、从搜索得到的频率对中寻找与脉率预测值之间的关系符合设定近似条件的频率值，并从中确定出最佳脉率值；

步骤四、将选择的最佳脉率值作为下一轮脉率提取所使用的脉率预测值。

2.如权利要求1所述的方法，其特征在于，所述步骤三具体包括：

将步骤二搜索得到的频率对记为[A,B]，其中，B近似为A的二倍，将A与B的一半即C构成近似频率对[A,C]储存到数组TA；

当数组TA中存在且仅存在1对[A,C]符合A与C均近似于预测脉率值，且A和C位置对应的幅度值均比其他位置的幅值大，则当前为非运动状态，从A和C二者选取其一为最佳脉率值；否则，当前有运动干扰，找到符合A与C均近似于预测脉率值的[A,C]，以及符合A或C更近似于预测脉率值的[A,C]，从找到的近似于预测脉率值的A和C以及找到的更近似于预测脉率值的A或C中选取最佳脉率值；

所述近似于预测脉率值的条件是：数据在预测脉率值的一临近范围之内；

所述更加近似于预测脉率值的条件是：数据在预测脉率值的更小临近范围之内。

3.如权利要求2所述的方法，其特征在于，

在所述步骤二中，将频谱信号中的峰值对应的频率值存入峰值数组Pm；

在所述步骤三中，当数组TA中的所有[A,C]中的A和C均不近似于预测脉率值或TA数组的长度为0时，则遍历数组Pm：若存在且仅存在一个值X或一个值Y的二分之一即Y/2符合所述更加近似于预测脉率值的条件，则选取符合更加近似于预测脉率值条件的X或Y/2为最佳脉率值，否则采取脉率保护措施，即最佳脉率值沿用上一个循环的值。

4.如权利要求2或3所述的方法，其特征在于，步骤三中，进一步从找到的近似于预测脉率值的A和C以及找到的更近似于预测脉率值的A或C中选取备选脉率值；

当在本次脉率值提取过程中无法找到最佳脉率值时，进一步采用上一个循环中确定的备选脉率值替代本次使用的预测脉率值重复步骤三；如果仍不能获得最佳脉率值，则采取脉率保护措施，即最佳脉率值沿用上一个循环的值。

5.如权利要求1所述的方法，其特征在于，在所述步骤二中，将频谱信号中的峰值对应的频率值存入峰值数组Pm；将搜索得到的频率对记为[A,B]，其中，B近似为A的二倍，计算B/2=C，将A与C构成近似频率对[A,C]储存到数组TA；

所述步骤三包括：

将TA中的近似频率对[A,C]逐一与预测脉率值进行比较，分情况确定最佳脉率值：

a）当数组TA中存在且存在1对[A,C]符合A与C均近似于预测脉率值，且A和C位置对应的幅度值均比其他位置的幅值大，则判定当前为非运动状态，从A和C二者选取其一为最佳脉率值；以下情况b）～g）均为运动状态；

b）当在数组TA中存在2对[A,C]符合A与C均近似于预测脉率值时，从两个A和两个C中选取最接近预测脉率值的一个为最佳脉率值；选取次接近预测脉率值的一个为备选脉率值；

c）当在数组TA中存在且仅存在1对[A,C]符合A与C均近似于预测脉率值，但A和C不符合二者对应的幅度值均比其他位置的幅值大的条件，则从A和C二者选取最接近预测脉率值的一个为最佳脉率值，选取另一个为备选脉率值；

d）当在数组TA中存在且仅存在1对[A,C]符合A与C均近似于预测脉率值，同时存在1对[A,C]中的A或C亦近似于预测脉率值且符合更加近似于预测脉率值的条件时，从前一[A,C]中选取A、C其一为近似值1，从后一[A,C]中选取符合更加近似于预测脉率值条件的A或C确定为近似值2，从近似值1和近似值2中选取其一为最佳脉率值，选取另一个为备选脉率值；

e）当在数组TA中有且仅有1对[A,C]中的A或C符合更加近似于预测脉率值的条件，则选取符合更加近似于预测脉率值的A或C为最佳脉率值；备选脉率值沿用上一次循环的值；

f）当数组TA中的所有[A,C]中的A和C均不近似于预测脉率值或TA数组的长度为0时，则遍历数组Pm：若存在且仅存在一个值X或一个值Y的二分之一即Y/2符合所述更加近似于预测脉率值的条件，则选取符合更加近似于预测脉率值条件的X或Y/2为最佳脉率值，备选脉率值沿用上一次循环的值；否则属于情况g）；

g）先采用上一个循环中确定的备选脉率值替代预测脉率值重复步骤三；如果仍不能获得最佳脉率值，则采取脉率保护措施，即最佳脉率值沿用上一个循环的值，且备选脉率值沿用上一次循环的值。