一种运动心率测量方法及其穿戴式设备

1.一种运动心率测量方法，其特征在于，其包括如下步骤：
采集人体的心率信号，所述心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值；
采集人体的运动信号，所述运动信号进行数字滤波处理后，进行所述峰值检测得到运动频率和第二频谱峰值；
所述第一频谱峰值经第一预设规则计算得到心率权重C1，所述第二频谱峰值经第二预设规则计算得到运动权重C2；
对预设时间段内获取的所述心率频率在包括多个子心率区间Bi的预设心率区间内进行频率分布处理，得到每一个所述子心率区间Bi的心率频率概率密度P(Qi)；
对预设时间段内获取的所述运动频率在所述预设心率区间内进行频率分布处理，得到每一个所述子心率区间Bi的运动频率概率密度P(Gi)；
计算每一个所述子心率区间Bi内的所有的心率频率的平均值
根据公式(1)计算得到心率HR：
2.根据权利要求1所述的运动心率测量方法，其特征在于，对预设时间段内获取的所述心率频率在包括多个子心率区间Bi的预设心率区间内进行频率分布处理，得到每一个所述子心率区间Bi的心率频率概率密度P(Qi)的步骤，包括：
获取预设时间段内，每一个所述子心率区间Bi内出现的所述心率频率的第一数量Qi和出现的所述运动频率的第二数量Gi；
根据公式(2)计算每一个所述子心率区间Bi的心率频率概率密度P(Qi)：
其中 若 则P
(Qi)＝0。
3.根据权利要求2所述的运动心率测量方法，其特征在于，获取预设时间段内，每一个所述子心率区间Bi内出现的所述心率频率的第一数量Qi和出现的所述运动频率的第二数量Gi的步骤之后，还包括：
根据公式(3)计算每一个所述子心率区间Bi的运动频率概率密度P(Gi)：
其中，若∑Gi＝0，则P(Gi)＝0.0000001。
4.根据权利要求1所述的运动心率测量方法，其特征在于，采集人体的心率信号，所述心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值的步骤，包括：
接收人体反射的激发光源发射出的光信号，并将所述光信号转换为电压信号；
将所述电压信号进行滤波处理、放大处理后得到放大信号；
将所述放大信号转换为数字信号；
所述数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值；
根据所述第一频谱峰值和第三预设规则调整所述激发光源的亮度，且根据所述第一频谱峰值和第四预设规则调整所述电压信号的放大倍数。
5.根据权利要求1所述的运动心率测量方法，其特征在于，根据公式(1)计算得到心率HR的步骤之后，包括：
输出显示所述心率HR和/或发送所述心率HR至外部。
6.一种实现权利要求1～5所述的运动心率测量方法的穿戴式设备，其特征在于，其包括第一采集模块、第二采集模块、第一峰值检测模块、第二峰值检测模块、心率权重计算模块、运动权重计算模块、第一概率密度处理模块、第二概率密度处理模块、平均值计算模块和心率计算模块；预设心率区间包括多个子心率区间Bi；所述第一采集模块，用于采集人体的心率信号；所述第二采集模块，用于采用人体的运动信号；所述第一峰值检测模块，用于对所述心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值；所述第二峰值检测模块，用于对所述运动信号进行数字滤波处理后，进行所述峰值检测得到运动频率和第二频谱峰值；所述心率权重计算模块，用于所述第一频谱峰值经第一预设规则计算得到心率权重C1；所述运动权重计算模块，用于所述第二频谱峰值经第二预设规则计算得到运动权重C2；所述第一概率密度处理模块，用于对预设时间段内获取的所述心率频率在所述预设心率区间进行频率分布处理，得到每一个所述子心率区间Bi的心率频率概率密度P(Qi)；所述第二概率密度处理模块，用于对预设时间段内获取的所述运动频率在所述预设心率区间进行频率分布处理，得到每一个所述子心率区间Bi的运动频率概率密度P(Gi)；所述平均值计算模块，用于计算每一个所述子心率区间Bi内的所有的心率频率的平均值 所述心率计算模块，用于根据公式(1)计算得到心率HR；
7.根据权利要求6所述的穿戴式设备，其特征在于，其还包括获取模块，所述获取模块，用于获取预设时间段内，每一个所述子心率区间Bi内出现的所述心率频率的第一数量Qi和出现的所述运动频率的第二数量Gi；所述第一概率密度处理模块，根据公式(2)计算每一个所述子心率区间Bi的心率频率概率密度P(Qi)：
其中 若
则P(Qi)＝0。
8.根据权利要求7所述的穿戴式设备，其特征在于，所述第二概率密度处理模块，还用于根据公式(3)计算每一个所述子心率区间Bi的运动频率概率密度P(Gi)：
若∑Gi＝0，则P(Gi)＝0.0000001。
9.根据权利要求6所述的穿戴式设备，其特征在于，其还包括调整模块，所述第一采集模块包括光电转换单元、放大单元和模数转换单元；所述光电转换单元，用于接收人体反射的激发光源发射出的光信号，并将所述光信号转换为电压信号；所述放大单元，用于将所述电压信号进行滤波处理、放大处理后得到放大信号；所述模数转换单元，用于将所述放大信号转换为数字信号；所述第一峰值检测模块，用于对所述数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值；所述调整模块，用于根据所述第一频谱峰值和第三预设规则调整所述激发光源的亮度，且根据所述第一频谱峰值和第四预设规则调整所述放大单元的放大倍数。
10.根据权利要求6所述的穿戴式设备，其特征在于，其还包括显示模块和/或输出模块，所述显示模块，用于输出显示所述心率HR，所述输出模块，用于发送所述心率HR至外部。

一种运动心率测量方法及其穿戴式设备\n技术领域\n[0001] 本发明涉及移动健康技术领域，尤其涉及一种运动心率测量方法及其穿戴式设备。\n背景技术\n[0002] 现有的用于测量心率的穿戴式设备主要包括血氧计、运动胸带、运动耳塞以及运动腕带。血氧计主要适用于静态测量心率，因此，运动时测量的心率的精度度不高。此外，运动胸带采用的是心电测量原理，由于体表心电信号较微弱，因此，运动胸带采用的电极较大，且需要运动胸带紧贴皮肤，以致穿戴不适，且检测到心电信号易受干扰，测量得到的心率精确度不高。此外，运动耳塞的体积小，因此，信号强度弱，以致测量得到的心率精确度不高。此外，运动腕带在运动时受到的噪声干扰比较大，为此，采用傅里叶变换抑制噪声干扰，但是，傅里叶变换计算量较大，以致不适宜应用到运动腕带的嵌入式系统中。\n发明内容\n[0003] 有鉴于此，实有必要提供一种计算量少，且测量精确度高的穿戴式设备，以及该穿戴式设备采用的运动心率测量方法。\n[0004] 一种运动心率测量方法，其包括如下步骤：\n[0005] 采集人体的心率信号，心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。\n[0006] 采集人体的运动信号，运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。\n[0007] 第一频谱峰值经第一预设规则计算得到心率权重C1，第二频谱峰值经第二预设规则计算得到运动权重C2。\n[0008] 对预设时间段内获取的心率频率在包括多个子心率区间Bi的预设心率区间内进行频率分布处理，得到每一个子心率区间Bi的心率频率概率密度P(Qi)。\n[0009] 对预设时间段内获取的运动频率在预设心率区间内进行频率分布处理，得到每一个子心率区间Bi的运动频率概率密度P(Gi)。\n[0010] 计算每一个子心率区间Bi内的所有的心率频率的平均值\n[0011] 根据公式(1)计算得到心率HR：\n[0012]\n[0013] 优选地，对预设时间段内获取的心率频率在包括多个子心率区间Bi的预设心率区间内进行频率分布处理，得到每一个子心率区间Bi的心率频率概率密度P(Qi)的步骤，包括：\n[0014] 获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。\n[0015] 根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：\n[0016] 其 中 c ∈ (0 ，1 ) ，若\n则P(Qi)＝0。\n[0017] 优选地，获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi的步骤之后，包括：\n[0018] 根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：\n[0019] 其中，若∑Gi＝0，则P(Gi)＝0.0000001。\n[0020] 优选地，采集人体的心率信号，心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值的步骤，包括：\n[0021] 接收人体反射的激发光源发射出的光信号，并将光信号转换为电压信号。\n[0022] 将电压信号进行滤波处理、放大处理后得到放大信号。\n[0023] 将放大信号转换为数字信号。\n[0024] 数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。\n[0025] 根据第一频谱峰值和第三预设规则调整激发光源的亮度，且根据第一频谱峰值和第四预设规则调整电压信号的放大倍数。\n[0026] 优选地，根据公式(1)计算得到心率HR的步骤之后，包括：\n[0027] 输出显示心率HR和/或发送心率HR至外部。\n[0028] 一种实现上述的运动心率测量方法的穿戴式设备，其包括第一采集模块、第二采集模块、第一峰值检测模块、第二峰值检测模块、心率权重计算模块、运动权重计算模块、第一概率密度处理模块、第二概率密度处理模块、平均值计算模块和心率计算模块。预设心率区间包括多个子心率区间Bi。第一采集模块，用于采集人体的心率信号。第二采集模块，用于采用人体的运动信号。第一峰值检测模块，用于对心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。第二峰值检测模块，用于对运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。心率权重计算模块，用于第一频谱峰值经第一预设规则计算得到心率权重C1。运动权重计算模块，用于第二频谱峰值经第二预设规则计算得到心运动权重C2。第一概率密度处理模块，用于对预设时间段内获取的心率频率在预设心率区间进行频率分布处理，得到每一个子心率区间Bi的心率频率概率密度P(Qi)。第二概率密度处理模块，用于对预设时间段内获取的运动频率在预设心率区间进行频率分布处理，得到每一个子心率区间Bi的运动频率概率密度P(Gi)。平均值计算模块，用于计算每一个子心率区间Bi内的所有的心率频率的平均值 心率计算模块，用于根据公式(1)计算得到心率HR：\n[0029]\n[0030] 优选地，穿戴式设备还包括获取模块。获取模块，用于获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。第一概率密度处理模块，根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：\n[0031] 其 中 c ∈ (0 ，1 ) ，若\n则P(Qi)＝0。\n[0032] 优选地，第二概率密度处理模块，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：\n[0033] 若∑Gi＝0，则P(Gi)＝0.0000001。\n[0034] 优选地，穿戴式设备还包括调整模块。第一采集模块包括光电转换单元、放大单元和模数转换单元。光电转换单元，用于接收人体反射的激发光源发射出的光信号，并将光信号转换为电压信号。放大单元，用于将电压信号进行滤波处理、放大处理后得到放大信号。\n模数转换单元，用于将放大信号转换为数字信号。第一峰值检测模块，用于对数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。调整模块，用于根据第一频谱峰值和第三预设规则调整激发光源的亮度，且根据第一频谱峰值和第四预设规则调整所述放大单元的放大倍数。\n[0035] 优选地，穿戴式设备还包括显示模块和/或输出模块，显示模块，用于输出显示心率HR，输出模块，用于发送心率HR至外部。\n[0036] 本发明中的穿戴式设备，以每一个心率子区间Bi的心率频率概率密度P(Qi)、运动频率概率密度P(Gi)和平均值 以及心率权重C1、运动权重C2作为加权平均算法的输入量计算心率，减少了计算量，且提升了计算精确度。\n附图说明\n[0037] 图1为一种实施方式穿戴式设备的功能模块示意图。\n[0038] 图2为另一种实施方式穿戴式设备的功能模块示意图。\n[0039] 图3为另一种实施方式穿戴式设备的功能模块示意图。\n[0040] 图4为一种实施方式运动心率测量方法的流程示意图。\n具体实施方式\n[0041] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用来限定本发明。\n[0042] 如图1所示，其为一种实施方式穿戴式设备的功能模块示意图。该穿戴式设备包括第一采集模块1、第二采集模块2、第一峰值检测模块3、第二峰值检测模块4、获取模块5、心率权重计算模块6、运动权重计算模块7、第一概率密度处理模块8、第二概率密度处理模块\n9、平均值计算模块10、心率计算模块11、调整模块13和显示模块12。该第一采集模块1包括激发光源101、光电转换单元102、放大单元103和模数转换单元104。该光电转换单元102包括光电传感器。放大单元103包括放大器。模数转换单元104包括模数转换器。第一峰值检测模块3包括数字滤波器和峰值检测单元。第二峰值检测模块4包括数字滤波器和峰值检测单元。第二采集模块2包括加速度计。\n[0043] 光电转换单元102，用于接收人体反射的激发光源101发射出的光信号，并将光信号转换为电压信号。放大单元103，用于将电压信号进行滤波处理、放大处理后得到放大信号。模数转换单元104，用于将放大信号转换为数字信号。第一峰值检测模块3，用于对数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。第二采集模块2，用于采用人体的运动信号。(例如：第二采集模块2包括加速度计，该加速度计用于采集人体运动时的加速度，即为运动信号)第二峰值检测模块4，用于对运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。调整模块13，用于根据第一频谱峰值和第三预设规则调整激发光源101的亮度(通过该第三预设规则和第一频谱峰值负反馈调节激发光源101亮度，即若第一频谱峰值大时，降低激发光源101的亮度，以减小后续的第一频谱峰值，若第一频谱峰值小时，增强激发光源101的亮度，以增大后续的第一频谱峰值)，且根据第一频谱峰值和第四预设规则调整放大单元103的放大倍数(通过该第四预设规则和第二频谱峰值负反馈调节放大单元103的放大倍数，即若第二频谱峰值大时，降低放大单元\n103的放大倍数，以减小后续的第二频谱峰值，若第二频谱峰值小时，增大放大单元103的放大倍数，以增大后续的第二频谱峰值)。心率权重计算模块6，用于第一频谱峰值经第一预设规则计算得到心率权重C1。(通过该第一预设规则计算得到的心率权重C1与第一频谱峰值成正比，即得到的第一频谱峰值越大，则得到的心率权重C1越大)。运动权重计算模块7，用于第二频谱峰值经第二预设规则计算得到心运动权重C2。(通过该第二预设规则计算得到的运动权重C2与第二频谱峰值成正比，即得到的第二频谱峰值越大，则得到的运动权重C2越大)获取模块5，用于获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。第一概率密度处理模块8，根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：\n[0044] 其中 c∈(0，1)，其中，若\n[0045] 则P(Qi)＝0。\n[0046] 第二概率密度处理模块9，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：\n[0047] 若∑Gi＝0，则P(Gi)＝0.0000001。\n[0048] 平均值计算模块10，用于计算每一个子心率区间Bi内的所有的心率频率的平均值心率计算模块11，用于根据公式(1)计算得到心率HR：\n[0049] 显示模块12，用于输出显示心率HR。\n[0050] 如图2所示，其为另一种实施方式穿戴式设备的功能模块示意图。该穿戴式设备包括第一采集模块1、第二采集模块2、第一峰值检测模块3、第二峰值检测模块4、获取模块5、心率权重计算模块6、运动权重计算模块7、第一概率密度处理模块8、第二概率密度处理模块9、平均值计算模块10、心率计算模块11、调整模块13和输出模块12。该第一采集模块1包括激发光源101、光电转换单元102、放大单元103和模数转换单元104。该光电转换单元102包括光电传感器。放大单元103包括放大器。模数转换单元104包括模数转换器。第一峰值检测模块3包括数字滤波器和峰值检测单元。第二峰值检测模块4包括数字滤波器和峰值检测单元。第二采集模块2包括加速度计。\n[0051] 光电转换单元102，用于接收人体反射的激发光源101发射出的光信号，并将光信号转换为电压信号。放大单元103，用于将电压信号进行滤波处理、放大处理后得到放大信号。模数转换单元104，用于将放大信号转换为数字信号。第一峰值检测模块3，用于对数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。第二采集模块2，用于采用人体的运动信号。(例如：第二采集模块2包括加速度计，该加速度计用于采集人体运动时的加速度，即为运动信号)第二峰值检测模块4，用于对运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。调整模块13，用于根据第一频谱峰值和第三预设规则调整激发光源101的亮度(通过该第三预设规则调整的激发光源101亮度与第一频谱峰值成反比，即第一频谱峰值越大，则降低激发光源101的亮度)，且根据第一频谱峰值和第四预设规则调整放大单元103的放大倍数(通过该第四预设规则调整的放大单元103放大倍数与第二频谱峰值成反比，即第二频谱峰值越大，则降低放大单元103的放大倍数)。心率权重计算模块6，用于第一频谱峰值经第一预设规则计算得到心率权重C1(通过该第一预设规则计算得到的心率权重C1与第一频谱峰值成正比，即得到的第一频谱峰值越大，则得到的心率权重C1越大)。运动权重计算模块7，用于第二频谱峰值经第二预设规则计算得到心运动权重C2(通过该第二预设规则计算得到的运动权重C2与第二频谱峰值成正比，即得到的第二频谱峰值越大，则得到的运动权重C2越大)。获取模块5，用于获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。第一概率密度处理模块8，根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：\n[0052] 其中 c∈(0，1)，其中，若\n[0053] 则P(Qi)＝0。\n[0054] 第二概率密度处理模块9，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：\n[0055] 若∑Gi＝0，则P(Gi)＝0.0000001。\n[0056] 平均值计算模块10，用于计算每一个子心率区间Bi内的所有的心率频率的平均值心率计算模块11，用于根据公式(1)计算得到心率HR：\n[0057] 输出模块12，用于发送心率HR至外部。\n[0058] 如图3所示，其为另一种实施方式穿戴式设备的功能模块示意图。该穿戴式设备包括第一采集模块1、第二采集模块2、第一峰值检测模块3、第二峰值检测模块4、获取模块5、心率权重计算模块6、运动权重计算模块7、第一概率密度处理模块8、第二概率密度处理模块9、平均值计算模块10、心率计算模块11、调整模块13、输出模块12和显示模块14。该第一采集模块1包括激发光源101、光电转换单元102、放大单元103和模数转换单元104。该光电转换单元102包括光电传感器。放大单元103包括放大器。模数转换单元104包括模数转换器。第一峰值检测模块3包括数字滤波器和峰值检测单元。第二峰值检测模块4包括数字滤波器和峰值检测单元。第二采集模块2包括加速度计。\n[0059] 光电转换单元102，用于接收人体反射的激发光源101发射出的光信号，并将光信号转换为电压信号。放大单元103，用于将电压信号进行滤波处理、放大处理后得到放大信号。模数转换单元104，用于将放大信号转换为数字信号。第一峰值检测模块3，用于对数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。第二采集模块2，用于采用人体的运动信号。(例如：第二采集模块2包括加速度计，该加速度计用于采集人体运动时的加速度，即为运动信号)第二峰值检测模块4，用于对运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。调整模块13，用于根据第一频谱峰值和第三预设规则调整激发光源101的亮度(通过该第三预设规则调整的激发光源101亮度与第一频谱峰值成反比，即第一频谱峰值越大，则降低激发光源101的亮度)，且根据第一频谱峰值和第四预设规则调整放大单元103的放大倍数(通过该第四预设规则调整的放大单元103放大倍数与第二频谱峰值成反比，即第二频谱峰值越大，则降低放大单元103的放大倍数)。心率权重计算模块6，用于第一频谱峰值经第一预设规则计算得到心率权重C1(通过该第一预设规则计算得到的心率权重C1与第一频谱峰值成正比，即得到的第一频谱峰值越大，则得到的心率权重C1越大)。运动权重计算模块7，用于第二频谱峰值经第二预设规则计算得到心运动权重C2(通过该第二预设规则计算得到的运动权重C2与第二频谱峰值成正比，即得到的第二频谱峰值越大，则得到的运动权重C2越大)。获取模块5，用于获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。第一概率密度处理模块8，根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：\n[0060] 其中 c∈(0，1)，其中，若\n[0061] 则P(Qi)＝0。\n[0062] 第二概率密度处理模块9，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：\n[0063] 若∑Gi＝0，则P(Gi)＝0.0000001。\n[0064] 平均值计算模块10，用于计算每一个子心率区间Bi内的所有的心率频率的平均值心率计算模块11，用于根据公式(1)计算得到心率HR：\n[0065] 输出模块12，用于发送心率HR至外部，\n且显示模块14，用于输出显示心率HR。\n[0066] 如图4所示，其为一种实施方式运动心率测量方法的流程示意图。该运动心率测量方法包括如下步骤：\n[0067] 步骤S1，接收人体反射的激发光源发射出的光信号，并将光信号转换为电压信号。\n将电压信号进行滤波处理、放大处理后得到放大信号。将放大信号转换为数字信号。数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。\n[0068] 步骤S2，采集人体的运动信号，运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。\n[0069] 步骤S3，第一频谱峰值经第一预设规则计算得到心率权重C1，第二频谱峰值经第二预设规则计算得到运动权重C2。\n[0070] 步骤S4，根据第一频谱峰值和第三预设规则调整激发光源的亮度，且根据第一频谱峰值和第四预设规则调整电压信号的放大倍数。\n[0071] 步骤S5，获取预设时间段内，在预设心率区间内的每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。\n[0072] 步骤S6，根据心率频率概率密度计算公式计算每一个子心率区间Bi的心率频率概率密度P(Qi)。心率频率概率密度计算公式为：\n[0073] 其中 c∈ (0 ，1) ，其中 ，若\n则P(Qi)＝0。\n[0074] 步骤S7，根据运动频率概率密度计算公式计算每一个子心率区间Bi的运动频率概率密度P(Gi)。运动频率概率密度计算公式为：\n[0075] 其中，若∑Gi＝0，则P(Gi)＝0.0000001。\n[0076] 步骤S8，计算每一个子心率区间Bi内的所有的心率频率的平均值[0077] 步骤S9，根据心率计算公式计算得到心率HR。该心率计算公式为：\n[0078]\n[0079] 步骤S10，输出显示心率HR和/或发送心率HR至外部。\n[0080] 以上对发明的具体实施方式进行了详细说明，但其只作为范例，本发明并不限制与以上描述的具体实施方式。对于本领域的技术人员而言，任何对该发明进行的等同修改或替代也都在本发明的范畴之中，因此，在不脱离本发明的精神和原则范围下所作的均等变换和修改、改进等，都应涵盖在本发明的范围内。