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

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

Info

Publication number
CN105939657A
CN105939657A CN201580004021.3A CN201580004021A CN105939657A CN 105939657 A CN105939657 A CN 105939657A CN 201580004021 A CN201580004021 A CN 201580004021A CN 105939657 A CN105939657 A CN 105939657A
Authority
CN
China
Prior art keywords
heart rate
frequency
module
probability density
sub
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201580004021.3A
Other languages
English (en)
Other versions
CN105939657B (zh
Inventor
杜凯萌
傅荃炜
邓志强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sky Light Industrial Ltd
Original Assignee
Sky Light Industrial Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sky Light Industrial Ltd filed Critical Sky Light Industrial Ltd
Publication of CN105939657A publication Critical patent/CN105939657A/zh
Application granted granted Critical
Publication of CN105939657B publication Critical patent/CN105939657B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802—Sensor mounted on worn items
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Power Engineering (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)

Abstract

一种运动心率测量方法和穿戴式设备。该方法包括：预设心率区间包括多个子心率区间Bi。对采集的心率信号、运动信号分别进行峰值检测得到心率频率和第一频谱峰值、运动频率和第二频谱峰值。第一频率峰值、第二频率峰值分别经预设规则计算得到心率权重C1、运动权重C2。根据预设心率区间对预设时间段内获取的心率频率、运动频率分别经频率分布处理得到每一个子心率区间Bi的心率频率概率密度P(Qi)、运动频率概率密度P(Gi)。以每一个心率子区间Bi的心率频率概率密度P(Qi)、运动频率概率密度P(Gi)和平均值，以及心率权重C1、运动权重C2作为加权平均算法的输入量计算心率，该方法减少了计算量，且提升了精确度。

Description

一种运动心率测量方法及其穿戴式设备
技术领域
本发明涉及移动健康技术领域，尤其涉及一种运动心率测量方法及其穿戴式设备。
背景技术
现有的用于测量心率的穿戴式设备主要包括血氧计、运动胸带、运动耳塞以及运动腕带。血氧计主要适用于静态测量心率，因此，运动时测量的心率的精度度不高。此外，运动胸带采用的是心电测量原理，由于体表心电信号较微弱，因此，运动胸带采用的电极较大，且需要运动胸带紧贴皮肤，以致穿戴不适，且检测到心电信号易受干扰，测量得到的心率精确度不高。此外，运动耳塞的体积小，因此，信号强度弱，以致测量得到的心率精确度不高。此外，运动腕带在运动时受到的噪声干扰比较大，为此，采用傅里叶变换抑制噪声干扰，但是，傅里叶变换计算量较大，以致不适宜应用到运动腕带的嵌入式系统中。
发明内容
有鉴于此，实有必要提供一种计算量少，且测量精确度高的穿戴式设备，以及该穿戴式设备采用的运动心率测量方法。
一种运动心率测量方法，其包括如下步骤：
采集人体的心率信号，心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。
采集人体的运动信号，运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。
第一频率峰值经第一预设规则计算得到心率权重C1，第二频率峰值经第二预设规则计算得到运动权重C2。
对预设时间段内获取的心率频率在包括多个子心率区间Bi的预设心率区间内进行频率分布处理，得到每一个子心率区间Bi的心率频率概率密度P(Qi)。
对预设时间段内获取的运动频率在预设心率区间内进行频率分布处理，得到每一个子心率区间Bi的运动频率概率密度P(Gi)。
计算每一个子心率区间Bi内的所有的心率频率的平均值
根据公式(1)计算得到心率HR:
H R = Σ HR i ‾ * P ( Q i ) * C 1 P ( G i ) * C 2 Σ P ( Q i ) * C 1 P ( G i ) * C 2 - - - ( 1 ) .
优选地，对预设时间段内获取的心率频率在包括多个子心率区间Bi的预设心率区间内进行频率分布处理，得到每一个子心率区间Bi的心率频率概率密度P(Qi)的步骤，包括：
获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。
根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：
其中c∈(0，1)，若则P(Qi)＝0。
优选地，获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi的步骤之后，包括：
根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：
其中，若∑Gi＝0，则P(Gi)＝0.0000001。
优选地，采集人体的心率信号，心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值的步骤，包括：
接收人体反射的激发光源发射出的光信号，并将光信号转换为电压信号。
将电压信号进行滤波处理、放大处理后得到放大信号。
将放大信号转换为数字信号。
数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。
根据第一频谱峰值和第三预设规则调整激发光源的亮度，且根据第一频谱峰值和第四预设规则调整电压信号的放大倍数。
优选地，根据公式(1)计算得到心率HR的步骤之后，包括：
输出显示心率HR和/或发送心率HR至外部。
一种实现上述的运动心率测量方法的穿戴式设备，其包括第一采集模块、第二采集模块、第一峰值检测模块、第二峰值检测模块、心率权重计算模块、运动权重计算模块、第一概率密度处理模块、第二概率密度处理模块、平均值计算模块和心率计算模块。预设心率区间包括多个子心率区间Bi。第一采集模块，用于采集人体的心率信号。第二采集模块，用于采用人体的运动信号。第一峰值检测模块，用于对心率信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。第二峰值检测模块，用于对运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。心率权重计算模块，用于第一频率峰值经第一预设规则计算得到心率权重C1。运动权重计算模块，用于第二频率峰值经第二预设规则计算得到心运动权重C2。第一概率密度处理模块，用于对预设时间段内获取的心率频率在预设心率区间进行频率分布处理，得到每一个子心率区间Bi的心率频率概率密度P(Qi)。第二概率密度处理模块，用于对预设时间段内获取的运动频率在预设心率区间进行频率分布处理，得到每一个子心率区间Bi的运动频率概率密度P(Gi)。平均值计算模块，用于计算每一个子心率区间Bi内的所有的心率频率的平均值心率计算模块，用于根据公式(1)计算得到心率HR：
H R = Σ HR i ‾ * P ( Q i ) * C 1 P ( G i ) * C 2 Σ P ( Q i ) * C 1 P ( G i ) * C 2 - - - ( 1 ) .
优选地，穿戴式设备还包括获取模块。获取模块，用于获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。第一概率密度处理模块，根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：
其中c∈(0，1)，若则P(Qi)＝0。
优选地，第二概率密度处理模块，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：
若∑Gi＝0，则P(Gi)＝0.0000001。
优选地，穿戴式设备还包括调整模块。第一采集模块包括光电转换单元、放大单元和模数转换单元。光电转换单元，用于接收人体反射的激发光源发射出的光信号，并将光信号转换为电压信号。放大单元，用于将电压信号进行滤波处理、放大处理后得到放大信号。模数转换单元，用于将放大信号转换为数字信号。第一峰值检测模块，用于对数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。调整模块，用于根据第一频谱峰值和第三预设规则调整激发光源的亮度，且根据第一频谱峰值和第四预设规则调整所述放大单元的放大倍数。
优选地，穿戴式设备还包括显示模块和/或输出模块，显示模块，用于输出显示心率HR，输出模块，用于发送心率HR至外部。
本发明中的穿戴式设备，以每一个心率子区间Bi的心率频率概率密度P(Qi)、运动频率概率密度P(Gi)和平均值以及心率权重C1、运动权重C2作为加权平均算法的输入量计算心率，减少了计算量，且提升了计算精确度。
附图说明
图1为一种实施方式穿戴式设备的功能模块示意图。
图2为另一种实施方式穿戴式设备的功能模块示意图。
图3为另一种实施方式穿戴式设备的功能模块示意图。
图4为一种实施方式运动心率测量方法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用来限定本发明。
如图1所示，其为一种实施方式穿戴式设备的功能模块示意图。该穿戴式设备包括第一采集模块1、第二采集模块2、第一峰值检测模块3、第二峰值检测模块4、获取模块5、心率权重计算模块6、运动权重计算模块7、第一概率密度处理模块8、第二概率密度处理模块9、平均值计算模块10、心率计算模块11、调整模块13和显示模块12。该第一采集模块1包括激发光源101、光电转换单元102、放大单元103和模数转换单元104。该光电转换单元102包括光电传感器。放大单元103包括放大器。模数转换单元104包括模数转换器。第一峰值检测模块3包括数字滤波器和峰值检测单元。第二峰值检测模块4包括数字滤波器和峰值检测单元。第二采集模块2包括加速度计。
光电转换单元102，用于接收人体反射的激发光源101发射出的光信号，并将光信号转换为电压信号。放大单元103，用于将电压信号进行滤波处理、放大处理后得到放大信号。模数转换单元104，用于将放大信号转换为数字信号。第一峰值检测模块3，用于对数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。第二采集模块2，用于采用人体的运动信号。(例如：第二采集模块2包括加速度计，该加速度计用于采集人体运动时的加速度，即为运动信号)第二峰值检测模块4，用于对运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。调整模块13，用于根据第一频谱峰值和第三预设规则调整激发光源101的亮度(通过该第三预设规则和第一频谱峰值负反馈调节激发光源101亮度，即若第一频谱峰值大时，降低激发光源101的亮度，以减小后续的第一频谱峰值，若第一频谱峰值小时，增强激发光源101的亮度，以增大后续的第一频谱峰值)，且根据第一频谱峰值和第四预设规则调整放大单元103的放大倍数(通过该第四预设规则和第二频谱峰值负反馈调节放大单元103的放大倍数，即若第二频谱峰值大时，降低放大单元103的放大倍数，以减小后续的第二频谱峰值，若第二频谱峰值小时，增大放大单元103的放大倍数，以增大后续的第二频谱峰值)。心率权重计算模块6，用于第一频率峰值经第一预设规则计算得到心率权重C1。(通过该第一预设规则计算得到的心率权重C1与第一频谱峰值成正比，即得到的第一频谱峰值越大，则得到的心率权重C1越大)。运动权重计算模块7，用于第二频率峰值经第二预设规则计算得到心运动权重C2。(通过该第二预设规则计算得到的运动权重C2与第二频谱峰值成正比，即得到的第二频谱峰值越大，则得到的运动权重C2越大)获取模块5，用于获取预设时间段内，每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。第一概率密度处理模块8，根据公式(2)计算每一个子心率区间Bi的心率频率概率密度P(Qi)：
其中c∈(0，1)，其中，若
则P(Qi)＝0。
第二概率密度处理模块9，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：
若∑Gi＝0，则P(Gi)＝0.0000001。
平均值计算模块10，用于计算每一个子心率区间Bi内的所有的心率频率的平均值心率计算模块11，用于根据公式(1)计算得到心率HR：
显示模块12，用于输出显示心率HR。
如图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包括加速度计。
光电转换单元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)：
其中c∈(0，1)，其中，若
则P(Qi)＝0。
第二概率密度处理模块9，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：
若∑Gi＝0，则P(Gi)＝0.0000001。
平均值计算模块10，用于计算每一个子心率区间Bi内的所有的心率频率的平均值心率计算模块11，用于根据公式(1)计算得到心率HR：
输出模块12，用于发送心率HR至外部。
如图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包括加速度计。
光电转换单元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)：
其中c∈(0，1)，其中，若
则P(Qi)＝0。
第二概率密度处理模块9，还用于根据公式(3)计算每一个子心率区间Bi的运动频率概率密度P(Gi)：
若∑Gi＝0，则P(Gi)＝0.0000001。
平均值计算模块10，用于计算每一个子心率区间Bi内的所有的心率频率的平均值心率计算模块11，用于根据公式(1)计算得到心率HR：
输出模块12，用于发送心率HR至外部，且显示模块14，用于输出显示心率HR。
如图4所示，其为一种实施方式运动心率测量方法的流程示意图。该运动心率测量方法包括如下步骤：
步骤S1，接收人体反射的激发光源发射出的光信号，并将光信号转换为电压信号。将电压信号进行滤波处理、放大处理后得到放大信号。将放大信号转换为数字信号。数字信号进行数字滤波处理后，进行峰值检测得到心率频率和第一频谱峰值。
步骤S2，采集人体的运动信号，运动信号进行数字滤波处理后，进行峰值检测得到运动频率和第二频谱峰值。
步骤S3，第一频率峰值经第一预设规则计算得到心率权重C1，第二频率峰值经第二预设规则计算得到运动权重C2。
步骤S4，根据第一频谱峰值和第三预设规则调整激发光源的亮度，且根据第一频谱峰值和第四预设规则调整电压信号的放大倍数。
步骤S5，获取预设时间段内，在预设心率区间内的每一个子心率区间Bi内出现的心率频率的第一数量Qi和出现的运动频率的第二数量Gi。
步骤S6，根据心率频率概率密度计算公式计算每一个子心率区间Bi的心率频率概率密度P(Qi)。心率频率概率密度计算公式为：
其中c∈(0，1)，其中，若则P(Qi)＝0。
步骤S7，根据运动频率概率密度计算公式计算每一个子心率区间Bi的运动频率概率密度P(Gi)。运动频率概率密度计算公式为：
其中，若∑Gi＝0，则P(Gi)＝0.0000001。
步骤S8，计算每一个子心率区间Bi内的所有的心率频率的平均值
步骤S9，根据心率计算公式计算得到心率HR。该心率计算公式为：
H R = Σ HR i ‾ * P ( Q i ) * C 1 P ( G i ) * C 2 Σ P ( Q i ) * C 1 P ( G i ) * C 2 .
步骤S10，输出显示心率HR和/或发送心率HR至外部。
以上对发明的具体实施方式进行了详细说明，但其只作为范例，本发明并不限制与以上描述的具体实施方式。对于本领域的技术人员而言，任何对该发明进行的等同修改或替代也都在本发明的范畴之中，因此，在不脱离本发明的精神和原则范围下所作的均等变换和修改、改进等，都应涵盖在本发明的范围内。

Claims (10)

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

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2015/087200 WO2017028149A1 (zh) 2015-08-17 2015-08-17 一种运动心率测量方法及其穿戴式设备

Publications (2)

Publication Number Publication Date
CN105939657A true CN105939657A (zh) 2016-09-14
CN105939657B CN105939657B (zh) 2019-01-22

Family

ID=57153132

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580004021.3A Expired - Fee Related CN105939657B (zh) 2015-08-17 2015-08-17 一种运动心率测量方法及其穿戴式设备

Country Status (2)

Country Link
CN (1) CN105939657B (zh)
WO (1) WO2017028149A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106913325A (zh) * 2017-02-17 2017-07-04 北京小米移动软件有限公司 计步方法及装置
CN109330580A (zh) * 2018-11-22 2019-02-15 深圳市元征科技股份有限公司 一种心率监测方法、系统及相关组件
CN109480815A (zh) * 2018-12-29 2019-03-19 出门问问信息科技有限公司 心率测量方法及装置
CN110494079A (zh) * 2018-08-03 2019-11-22 高驰运动科技(深圳)有限公司 心率检测方法、装置、检测设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030198985A1 (en) * 2002-02-22 2003-10-23 Hogan Kirk J. Assay for nitrous oxide neurologic syndrome
CN102670190A (zh) * 2012-05-21 2012-09-19 南京大学 基于心率变异性非线性特性的充盈性心衰自动诊断方法
CN104224132A (zh) * 2014-09-26 2014-12-24 天彩电子（深圳）有限公司 睡眠监测装置及其监测方法
CN204049629U (zh) * 2014-06-27 2014-12-31 天彩电子（深圳）有限公司 一种心率监测仪
CN104462858A (zh) * 2014-12-30 2015-03-25 天津迈沃医药技术有限公司 基于多阶隐马尔可夫模型的健康预警方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5448515B2 (ja) * 2009-03-25 2014-03-19 シチズンホールディングス株式会社 生体信号測定装置
JP5682383B2 (ja) * 2011-03-09 2015-03-11 セイコーエプソン株式会社 拍動検出装置
JP6149037B2 (ja) * 2011-09-16 2017-06-14 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． ポータブル装置、心拍を決定する方法、システム及びコンピュータプログラム
CN103315728B (zh) * 2012-03-20 2015-12-09 深圳市飘浮互动科技有限公司 心率检测和显示方法及其装置
CN104027103A (zh) * 2013-03-06 2014-09-10 精工爱普生株式会社 生物体信息检测装置以及脉搏计
JP5794344B2 (ja) * 2014-04-22 2015-10-14 セイコーエプソン株式会社 生体情報測定装置及び方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030198985A1 (en) * 2002-02-22 2003-10-23 Hogan Kirk J. Assay for nitrous oxide neurologic syndrome
CN102670190A (zh) * 2012-05-21 2012-09-19 南京大学 基于心率变异性非线性特性的充盈性心衰自动诊断方法
CN204049629U (zh) * 2014-06-27 2014-12-31 天彩电子（深圳）有限公司 一种心率监测仪
CN104224132A (zh) * 2014-09-26 2014-12-24 天彩电子（深圳）有限公司 睡眠监测装置及其监测方法
CN104462858A (zh) * 2014-12-30 2015-03-25 天津迈沃医药技术有限公司 基于多阶隐马尔可夫模型的健康预警方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
涂岳文: "Holter系统中运动伪差自动识别的关键技术及算法研究", 《中国博士学位论文全文数据库（电子期刊）信息科技辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106913325A (zh) * 2017-02-17 2017-07-04 北京小米移动软件有限公司 计步方法及装置
CN110494079A (zh) * 2018-08-03 2019-11-22 高驰运动科技(深圳)有限公司 心率检测方法、装置、检测设备及存储介质
CN110494079B (zh) * 2018-08-03 2022-09-02 广东高驰运动科技有限公司 心率检测方法、装置、检测设备及存储介质
CN109330580A (zh) * 2018-11-22 2019-02-15 深圳市元征科技股份有限公司 一种心率监测方法、系统及相关组件
CN109480815A (zh) * 2018-12-29 2019-03-19 出门问问信息科技有限公司 心率测量方法及装置
CN109480815B (zh) * 2018-12-29 2021-10-29 出门问问信息科技有限公司 心率测量方法及装置

Also Published As

Publication number Publication date
CN105939657B (zh) 2019-01-22
WO2017028149A1 (zh) 2017-02-23

Similar Documents

Publication Publication Date Title
CN109907729B (zh) 一种睡眠时生命特征的检测方法
CN105939657A (zh) 一种运动心率测量方法及其穿戴式设备
CN103596492B (zh) 耳戴式生命体征监视器
CN103099615B (zh) 一种消除运动心电信号干扰的方法和装置
EP3106086B1 (en) Photoelectric-type pulse signal measurement method and measurement device
CN108009519B (zh) 一种光辐照信息监测方法及装置
CN108294737B (zh) 心率测量方法、装置及智能穿戴设备
CN102551727A (zh) 呼吸信息检测方法及装置
CN109091142B (zh) 一种检测肌肉内源性疲劳度的方法、装置及存储介质
KR101263759B1 (ko) 3축 가속도 센서를 이용한 칼로리 소모량 산출 장치 및 그 방법
CN106344023A (zh) 一种基于气压和加速度的非稳态呼吸波检测装置
US20160198964A1 (en) Vital signs measurement system, detecting method of the vital signs measurement system, and vital signs measurement earphone
CN205120130U (zh) 一种基于改进的fft计算步数和人体运动消耗卡路里的设备
CN105147293A (zh) 实现呼吸频率测量的系统及方法
CN201847678U (zh) 腰带式无线心率检测装置
CN106805943B (zh) 基于分数阶微分方程的血糖数据处理方法及装置
Prawiro et al. A heart rate detection method for low power exercise intensity monitoring device
WO2017113152A1 (zh) 基于多组光源的心率测量方法及实现该方法的穿戴式设备
CN106618570A (zh) 一种基于生物介电谱的皮肤生化指标检测方法及系统
CN202604793U (zh) 一种人体平衡评定及训练系统
CN115153580A (zh) 一种用于胎儿心电监测的皮肤阻抗测量方法和系统
KR100995318B1 (ko) 신호 검출 장치 및 방법
CN204394493U (zh) 一种智能血压计
CN105982665A (zh) 一种从人耳耳道采集脑电波信号的方法
Shi et al. Design, Development and Testing of a Wearable sEMG Acquisition System

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190122

Termination date: 20190817

CF01 Termination of patent right due to non-payment of annual fee