Nothing Special   »   [go: up one dir, main page]

CN105997043B - A kind of pulse frequency extracting method based on wrist wearable device - Google Patents

A kind of pulse frequency extracting method based on wrist wearable device Download PDF

Info

Publication number
CN105997043B
CN105997043B CN201610471150.XA CN201610471150A CN105997043B CN 105997043 B CN105997043 B CN 105997043B CN 201610471150 A CN201610471150 A CN 201610471150A CN 105997043 B CN105997043 B CN 105997043B
Authority
CN
China
Prior art keywords
pulse rate
value
acceleration
motion
state
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.)
Expired - Fee Related
Application number
CN201610471150.XA
Other languages
Chinese (zh)
Other versions
CN105997043A (en
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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201610471150.XA priority Critical patent/CN105997043B/en
Publication of CN105997043A publication Critical patent/CN105997043A/en
Application granted granted Critical
Publication of CN105997043B publication Critical patent/CN105997043B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring pulse rate or heart rate
    • A61B5/02444Details of sensor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring pulse rate or heart rate
    • A61B5/02438Measuring pulse rate or heart rate with portable devices, e.g. worn by the patient
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring pulse rate or heart rate
    • A61B5/0245Measuring pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements 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/6813Specially adapted to be attached to a specific body part
    • A61B5/6824Arm or wrist
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7207Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02Details of sensors specially adapted for in-vivo measurements
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明公开了一种基于腕式可穿戴设备的脉率提取方法。本发明首先通过对加速度传感器采集的加速度信息进行分析,将当前的运动状态进行了分类。对于不同运动状态,采用不同的算法获得脉率信息。对于运动下的测试,还分析了其运动是否规律,对于规律的运动和不规律的运动,分别采用不同的算法获得脉率信息。这种分类方式使得脉率的获取更为准确。同时,也解决了对于存在噪声信号无法获取脉率的问题。

The invention discloses a pulse rate extraction method based on a wrist wearable device. The present invention first classifies the current motion state by analyzing the acceleration information collected by the acceleration sensor. For different exercise states, different algorithms are used to obtain pulse rate information. For the test under exercise, it is also analyzed whether its movement is regular. For regular exercise and irregular exercise, different algorithms are used to obtain pulse rate information. This classification method makes the acquisition of pulse rate more accurate. At the same time, the problem that the pulse rate cannot be obtained for the noise signal is also solved.

Description

一种基于腕式可穿戴设备的脉率提取方法A pulse rate extraction method based on wrist wearable device

技术领域technical field

本发明涉及腕式可穿戴设备领域,具体涉及一种基于腕式可穿戴设备的脉率提取方法。The invention relates to the field of wrist wearable devices, in particular to a pulse rate extraction method based on wrist wearable devices.

背景技术Background technique

随着低功耗技术和生理参数监测技术的发展,可穿戴设备在医疗设备领域掀起了一场革命。在可穿戴设备监测的各种参数中,脉率尤其重要。因为脉率监测不仅可以预防运动中心脏病的突然发作,而且在日常训练中可以提供指导作用。通常,人们通过捆绑在胸前的脉率带获取ECG来实现脉率监测。这种脉率带需要勒紧胸口来获取可靠的信号,舒适性较差。With the development of low power consumption technology and physiological parameter monitoring technology, wearable devices have started a revolution in the field of medical equipment. Among the various parameters monitored by wearable devices, pulse rate is particularly important. Because pulse rate monitoring can not only prevent sudden heart attacks during exercise, but also provide guidance in daily training. Usually, people realize pulse rate monitoring by obtaining ECG through the pulse rate belt that is tied on the chest. This kind of pulse rate belt needs to tighten the chest to obtain a reliable signal, and the comfort is poor.

PPG(光电溶剂脉搏波)信号是表征毛细血管中血液容量变化的物理量。伴随着心脏的搏动,血液流向毛细血管,从而导致血管中血液容量发生变化。因此PPG作为一种生理信号,与ECG(心电)信号具有同样的生物学意义。人类血管中血红蛋白吸收绿光的能力比吸收其它光的能力强,根据这一特性,可以通过向表皮发射绿光然后检测毛细血管反射的绿光强度的变化来获取PPG信号。绿光的产生和检测可以通过光电传感器来实现。现有技术中已出现通过腕表采集的PPG信号,替代脉率带采集的ECG信号提取脉率的方法,该方法较大程度提升了舒适性。The PPG (photoelectric solvent pulse wave) signal is a physical quantity that characterizes changes in blood volume in capillaries. As the heart beats, blood flows into the capillaries, causing changes in the volume of blood in the vessels. Therefore, as a physiological signal, PPG has the same biological significance as the ECG (cardiac electricity) signal. The ability of hemoglobin in human blood vessels to absorb green light is stronger than that of other light. According to this characteristic, PPG signals can be obtained by emitting green light to the epidermis and then detecting changes in the intensity of green light reflected by capillaries. The generation and detection of green light can be realized by photoelectric sensors. In the prior art, there has been a method of extracting the pulse rate by replacing the ECG signal collected by the pulse rate belt with the PPG signal collected by the watch, which greatly improves the comfort.

但是,PPG信号很容易被运动产生的噪声所干扰,这些噪声大部分是由于运动过程中传感器相对于皮肤的滑动所产生的,具有周期性而且频率与脉率接近。因此在运动状态下从PPG信号中提取脉率非常困难。常见的腕式设备很难对运动状态下脉率实现准确测量。However, the PPG signal is easily disturbed by motion-generated noise, most of which are generated by the sliding of the sensor relative to the skin during motion, which is periodic and has a frequency close to the pulse rate. Therefore, it is very difficult to extract the pulse rate from the PPG signal in the state of exercise. It is difficult for common wrist devices to accurately measure the pulse rate during exercise.

发明内容Contents of the invention

有鉴于此,本发明提供了一种基于腕式可穿戴设备的脉率提取方法,能够在运动状态和非运动状态下,获得较为精准的脉率。In view of this, the present invention provides a pulse rate extraction method based on a wrist-type wearable device, which can obtain a more accurate pulse rate in an exercise state and a non-exercise state.

一种基于腕式可穿戴设备的脉率提取方法,所述腕式可穿戴设备主要由加速度传感器和光电传感器构成;提取方法具体包括如下步骤:A pulse rate extraction method based on a wrist-type wearable device, the wrist-type wearable device is mainly composed of an acceleration sensor and a photoelectric sensor; the extraction method specifically includes the following steps:

步骤一、开启腕式可穿戴设备,加速度传感器和光电传感器各自进行采集工作;待加速度传感器采集n个加速度值后,将采样的n个加速度数据储存在数组d中;执行步骤二;Step 1. Turn on the wrist wearable device, and the acceleration sensor and the photoelectric sensor perform collection work respectively; after the acceleration sensor collects n acceleration values, store the sampled n acceleration data in the array d; perform step 2;

步骤二、将步骤一中采集的n个加速度值求模,并依次与阈值q进行比较;若模值大于阈值q的个数超过设定值x时,且x=n×10%,则说明当前处于运动状态,执行步骤四;若模值大于阈值q的个数未超过设定值x时,则说明当前处于静止状态,执行步骤三;Step 2. Calculate the modulus of the n acceleration values collected in step 1, and compare them with the threshold q in turn; if the number of modulus values greater than the threshold q exceeds the set value x, and x=n×10%, then it means Currently in a state of motion, go to step 4; if the number of modulus values greater than the threshold q does not exceed the set value x, it means that you are currently in a static state, go to step 3;

步骤三、采用光电传感器进行脉率提取,得到光电容积脉搏波PPG信号,通过对PPG信号中相邻R点的间期算法获得静止状态下的脉率;Step 3, using the photoelectric sensor to extract the pulse rate to obtain the photoplethysmography PPG signal, and obtain the pulse rate in the static state by the interval algorithm of the adjacent R points in the PPG signal;

步骤四、根据运动分类算法将运动状态分为规律运动和非规律运动,其确定方法如下;Step 4, according to the motion classification algorithm, the motion state is divided into regular motion and irregular motion, and its determination method is as follows;

S41、从加速度模值中选出最大值amax S41, select the maximum value a max from the acceleration modulus

S42、对数组d中的加速度数据进行归一化处理;S42. Normalize the acceleration data in the array d;

S43、选取归一化后的数组d中的最大值pmaxS43. Select the maximum value p max in the normalized array d;

S44、设定阈值P,以pmax为基准值,在归一化后的数组d中遍历基准值左侧的其他数值,获取首个超过阈值P的数值,定义为pb;同样地,遍历基准值右侧的其他数值,获取首个超过阈值P的数值,定义为pc;当遍历完所有数值后,一旦有一侧没有获取到数值,则判断为非规律运动状态;执行步骤六;若pb和pc均获得,则执行步骤S45;S44. Set the threshold value P, take p max as the reference value, traverse other values on the left side of the reference value in the normalized array d, obtain the first value exceeding the threshold value P, and define it as p b ; similarly, traverse For other values on the right side of the reference value, obtain the first value that exceeds the threshold P, which is defined as p c ; after traversing all the values, once the value is not obtained on one side, it is judged as an irregular motion state; perform step 6; if Both p b and p c are obtained, then perform step S45;

S45、定义pb与pmax之间的距离为L1,pmax与pc之间的距离为L2;若L1和L2相等,则说明做规律运动,执行步骤五;否则,判定为非规律运动,执行步骤六;S45. Define the distance between p b and p max as L 1 , and the distance between p max and p c as L 2 ; if L 1 and L 2 are equal, it means that you are doing regular exercise and go to step 5; otherwise, determine For irregular exercise, perform step 6;

步骤五、对于规律运动状态下的脉率提取方法,采用光电传感器采集PPG信号,并提取脉率;Step 5. For the pulse rate extraction method under the regular exercise state, the photoelectric sensor is used to collect the PPG signal, and the pulse rate is extracted;

步骤六、对于非规律运动状态下的脉率提取方法:Step 6. For the pulse rate extraction method in the state of irregular exercise:

先采用光电传感器采集PPG信号,之后采用ANC算法对PPG信号进行去噪处理,通过提取PPG信号中相邻P点的间期,获得脉率。First, the photoelectric sensor is used to collect the PPG signal, and then the ANC algorithm is used to denoise the PPG signal, and the pulse rate is obtained by extracting the interval between adjacent P points in the PPG signal.

有益效果:Beneficial effect:

本发明首先通过对加速度传感器采集的加速度信息进行分析,将当前的运动状态进行了分类。对于不同运动状态,采用不同的算法获得脉率信息。对于运动下的测试,还分析了其运动是否规律,对于规律的运动和不规律的运动,分别采用不同的算法获得脉率信息。这种分类方式使得脉率的获取更为准确。同时,也解决了对于存在噪声信号无法获取脉率的问题。The present invention first classifies the current motion state by analyzing the acceleration information collected by the acceleration sensor. For different exercise states, different algorithms are used to obtain pulse rate information. For the test under exercise, it is also analyzed whether its movement is regular. For regular exercise and irregular exercise, different algorithms are used to obtain pulse rate information. This classification method makes the acquisition of pulse rate more accurate. At the same time, the problem that the pulse rate cannot be obtained for the noise signal is also solved.

附图说明Description of drawings

图1为本发明示意图。Fig. 1 is a schematic diagram of the present invention.

图2为静止状态心率提取。Figure 2 is the resting state heart rate extraction.

图3为运动分类算法示意图。Figure 3 is a schematic diagram of the motion classification algorithm.

图4为ANC算法原理图。Figure 4 is a schematic diagram of the ANC algorithm.

图5为频域处理示意图。Fig. 5 is a schematic diagram of frequency domain processing.

具体实施方式Detailed ways

下面结合附图并举实施例,对本发明进行详细描述。The present invention will be described in detail below with reference to the accompanying drawings and examples.

本发明提供了一种基于腕式可穿戴设备的脉率提取方法,其主要思想在于:The present invention provides a pulse rate extraction method based on a wrist wearable device, the main idea of which is:

由于在非规律运动状态下,光电传感器和加速度传感器所采集的图像存在噪声,无法实现对非规律运动状态下脉率的提取。所以,本发明基于此,提出了一种新的脉率提取方法。首先,本发明先对加速度传感器提取的数据进行分析,判断其运动状态和运动规律。对于非规律运动状态下所采集的信号,进行去噪处理,进而获得一个较为准确的脉率值。即:本发明实现了对光电传感器和加速度传感器提取的信号得分类处理,进而能够针对不同状态下的信号,采用不同的处理方法(处理方法为现有技术),进而获得较为准确的脉率。Due to the presence of noise in the image collected by the photoelectric sensor and the acceleration sensor in the state of irregular motion, it is impossible to extract the pulse rate in the state of irregular motion. Therefore, based on this, the present invention proposes a new pulse rate extraction method. Firstly, the present invention firstly analyzes the data extracted by the acceleration sensor, and judges its motion state and motion law. For the signal collected under the state of irregular motion, denoise processing is performed to obtain a more accurate pulse rate value. That is: the present invention realizes the classification processing of the signals extracted by the photoelectric sensor and the acceleration sensor, and then can adopt different processing methods (the processing method is the prior art) for the signals in different states, and then obtain a more accurate pulse rate.

其中,本发明的腕式可穿戴设备由加速度传感器模块、光电传感器和处理模块构成。Wherein, the wrist wearable device of the present invention is composed of an acceleration sensor module, a photoelectric sensor and a processing module.

具体实施方式如下,如图1所示:The specific implementation is as follows, as shown in Figure 1:

步骤一、开启腕式可穿戴设备,此时,加速度传感器和光电传感器将开始工作。光电传感器实时采集当前待测者的光电容积脉搏波(PPG)信号,加速度传感器模块对加速度信号进行采样,共采集n次,并将采样的n个数据储存在数组d中。执行步骤二。为了保证数据的准确性,本发明设定了一个时间阈值,从第一次采集时刻开始计时,一旦达到该时间阈值,则将数据d中的数据清除,重新进行采集,数组始终存储最新一段时间的加速度信息。Step 1. Turn on the wrist wearable device. At this time, the acceleration sensor and photoelectric sensor will start to work. The photoelectric sensor collects the photoplethysmography (PPG) signal of the current subject in real time, and the acceleration sensor module samples the acceleration signal for a total of n times, and stores the sampled n data in the array d. Go to step 2. In order to ensure the accuracy of the data, the present invention sets a time threshold, counting from the first collection time, once the time threshold is reached, the data in the data d will be cleared, and the collection will be carried out again, and the array will always store the latest period of time acceleration information.

步骤二、根据步骤一中采集的加速度信号进行分析,判断当前运动状态:Step 2. Analyze the acceleration signal collected in step 1 to determine the current motion state:

首先,根据之前的实验数据及经验,设定一个用于判断当前运动状态的阈值q;First, according to the previous experimental data and experience, set a threshold q for judging the current motion state;

之后,对数组d中的n个加速度值求模,并依次与阈值q进行比较。若模值大于阈值q的个数超过设定值x时,x<n,且x=n×10%,则说明当前处于运动状态,执行步骤四;若模值大于阈值q的个数未超过设定值x时,则说明当前处于静止状态,执行步骤三。After that, the n acceleration values in the array d are calculated modulo, and compared with the threshold q in turn. If the number of modulus values greater than the threshold q exceeds the set value x, x<n, and x=n×10%, it means that the current state is in motion, and step 4 is performed; if the number of modulus values greater than the threshold q does not exceed When the value x is set, it means that it is currently in a static state, and step 3 is performed.

步骤三、对处于静止状态下的脉率提取方法:Step 3. The pulse rate extraction method in a static state:

由于PPG信号是表征毛细血管中血液容量变化的物理量,伴随着心脏的搏动,血液流向毛细血管,从而导致血管中血液容量发生变化。所以在静止状态下PPG的波形如同ECG一样稳定,如图2所示,故对处于静止状态下的脉率提取方法为:采用光电传感器进行脉率提取,得到PPG信号,提取PPG信号中相邻R点,采用RR间期算法计算其重复周期TRR(RR间期),TRR倒数即为当前脉率。该方法为现有技术,再此不做赘述。Since the PPG signal is a physical quantity representing the change of blood volume in the capillary, along with the beating of the heart, the blood flows to the capillary, resulting in the change of the blood volume in the blood vessel. Therefore, the waveform of PPG in the static state is as stable as ECG, as shown in Figure 2, so the pulse rate extraction method in the static state is: use a photoelectric sensor to extract the pulse rate, obtain the PPG signal, and extract the adjacent pulse rate in the PPG signal. At point R, use the RR interval algorithm to calculate its repetition period T RR (RR interval), and the reciprocal of T RR is the current pulse rate. This method is a prior art, and will not be repeated here.

步骤四、当处于运动状态时,根据运动分类算法将运动状态分为规律运动和非规律运动。其具体判别方法如下:Step 4. When in a motion state, classify the motion state into regular motion and irregular motion according to a motion classification algorithm. The specific identification method is as follows:

S41、从加速度模值中选出最大值amax S41, select the maximum value a max from the acceleration modulus

S42、根据公式(1),对数组d进行归一化处理,获得新的数据p={p1,p2,...,pi,...,pn};S42. According to the formula (1), normalize the array d to obtain new data p={p 1 ,p 2 ,...,p i ,...,p n };

S43、在数据p中,选取数据p中的最大值pmaxS43. In the data p, select the maximum value p max in the data p;

S44、设定阈值P,以pmax为基准值,依次遍历数据p中基准值左侧的其他数值,获取首个超过阈值P的数值,分别定义为pb;同样地,遍历数据p中基准值右侧的其他数值,获取首个超过阈值P的数值,分别定义为pc;当遍历完所有数值后,一旦有一侧没有获取到数值,则判断为非规律运动状态;执行步骤六;若pb和pc均获得,则执行步骤S45;S44. Set the threshold value P, take p max as the reference value, sequentially traverse other values on the left side of the reference value in the data p, and obtain the first value exceeding the threshold value P, which is respectively defined as p b ; similarly, traverse the reference values in the data p For other values on the right side of the value, get the first value that exceeds the threshold P, which is defined as p c respectively; after traversing all the values, once the value is not obtained on one side, it is judged as an irregular motion state; perform step 6; if Both p b and p c are obtained, then perform step S45;

S45、通过pb、pmax和pc,进一步判断当前运动状态:S45. Further judge the current motion state through p b , p max and p c :

定义pb与pmax之间的距离为L1,pmax与pc之间的距离为L2;如图3所示,L1和L2实际代表了pb与pmax,pmax与pc之间的时间差。若L1和L2相等,则说明做规律运动,执行步骤五。否则,判定为非规律运动,执行步骤六。Define the distance between p b and p max as L 1 , and the distance between p max and p c as L 2 ; as shown in Figure 3, L 1 and L 2 actually represent p b and p max , p max and p max The time difference between pc . If L 1 and L 2 are equal, it means doing regular exercise, go to step 5. Otherwise, it is judged as irregular motion, and go to step 6.

步骤五、对于规律运动状态下的脉率提取方法:Step 5. For the pulse rate extraction method in the state of regular exercise:

采用现有常规技术,采用光电传感器采集PPG信号,对获得的PPG信号的加速度信号按如下方法进行处理后,便能够获得脉率。Using the existing conventional technology, the photoelectric sensor is used to collect the PPG signal, and the pulse rate can be obtained after the acceleration signal of the obtained PPG signal is processed according to the following method.

在规律运动情况下通常如图5所示,加速度信号所对应的功率谱曲线通常会出现两个波峰,这两个波峰分别对应规律运动产生的基波和谐波,而这两个波峰也会同样出现在PPG信号的功率谱曲线中,因此只要从PPG信号功率谱中滤去加速度信号功率谱就可以得到理想的PPG信号。故:In the case of regular motion, usually as shown in Figure 5, the power spectrum curve corresponding to the acceleration signal usually has two peaks, which correspond to the fundamental wave and harmonics generated by regular motion, and these two peaks will also It also appears in the power spectrum curve of the PPG signal, so as long as the acceleration signal power spectrum is filtered out from the PPG signal power spectrum, an ideal PPG signal can be obtained. Therefore:

第一步、对加速度信号和光电容积脉搏波信号分别做快速傅里叶变换得到各自的频谱Facc,FppgIn the first step, fast Fourier transform is performed on the acceleration signal and the photoplethysmography signal respectively to obtain respective frequency spectra F acc , F ppg .

第二步、将Facc和Fppg转化为功率谱Xacc和Xppg,并分别找出对应的最大幅值Amax和BmaxThe second step is to convert F acc and F ppg into power spectra X acc and X ppg , and find out the corresponding maximum amplitudes A max and B max respectively.

第三步、归一化处理,获之后,进行差异化处理,获得频谱中幅值最大的点。The third step, normalization processing, obtained and Afterwards, differential processing is performed to obtain the point with the largest amplitude in the spectrum.

第四步、将频谱中幅值最大的点所对应的频率就是脉率。The fourth step, the frequency corresponding to the point with the largest amplitude in the frequency spectrum is the pulse rate.

由于方法为常规技术,本发明的重点也不在于此,本发明仅做简要介绍。Since the method is a conventional technology, the focus of the present invention is not here, and the present invention is only briefly introduced.

步骤六、对于非规律运动状态下的脉率提取方法:Step 6. For the pulse rate extraction method in the state of irregular exercise:

由于在非规律运动状态下,光电传感器除了采集到PPG信号之外,还会采集到噪声干扰信号,如图4所示。为此,本发明采用ANC算法进行去噪处理。由于ANC算法为现有常规技术手段,不在过多赘述。ANC算法的基本原理为:通过加速度信号产生一个噪声干扰信号的最优估计信号,从光电传感器获取的信号中滤去该噪声估计信号就可以得到理想的PPG。Because in the state of irregular motion, the photoelectric sensor will also collect noise interference signals in addition to collecting PPG signals, as shown in Figure 4. For this reason, the present invention adopts ANC algorithm to carry out denoising processing. Since the ANC algorithm is an existing conventional technical means, details are not repeated here. The basic principle of the ANC algorithm is: the acceleration signal generates an optimal estimation signal of the noise interference signal, and the ideal PPG can be obtained by filtering the noise estimation signal from the signal obtained by the photoelectric sensor.

在获得理想的PPG信号后,通过计算PP间期TPP,TPP倒数即为当前脉率。该方法为现有技术,再此不做赘述。After obtaining an ideal PPG signal, the PP interval T PP is calculated, and the reciprocal of T PP is the current pulse rate. This method is a prior art, and will not be repeated here.

综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。To sum up, the above are only preferred embodiments of the present invention, and are not intended to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.

Claims (1)

1.一种基于腕式可穿戴设备的脉率提取方法,所述腕式可穿戴设备主要由加速度传感器和光电传感器构成;其特征在于,提取方法具体包括如下步骤:1. a pulse rate extraction method based on wrist-type wearable equipment, described wrist-type wearable equipment is mainly made of acceleration sensor and photoelectric sensor; It is characterized in that, extraction method specifically comprises the steps: 步骤一、开启腕式可穿戴设备,加速度传感器和光电传感器各自进行采集工作;待加速度传感器采集n个加速度值后,将采样的n个加速度数据储存在数组d中;执行步骤二;Step 1. Turn on the wrist wearable device, and the acceleration sensor and the photoelectric sensor perform collection work respectively; after the acceleration sensor collects n acceleration values, store the sampled n acceleration data in the array d; perform step 2; 步骤二、将步骤一中采集的n个加速度值求模,并依次与阈值q进行比较;若模值大于阈值q的个数超过设定值x时,且x=n×10%,则说明当前处于运动状态,执行步骤四;若模值大于阈值q的个数未超过设定值x时,则说明当前处于静止状态,执行步骤三;Step 2. Calculate the modulus of the n acceleration values collected in step 1, and compare them with the threshold q in turn; if the number of modulus values greater than the threshold q exceeds the set value x, and x=n×10%, then it means Currently in a state of motion, go to step 4; if the number of modulus values greater than the threshold q does not exceed the set value x, it means that you are currently in a static state, go to step 3; 步骤三、采用光电传感器进行脉率提取,得到光电容积脉搏波PPG信号,通过对PPG信号中相邻R点的间期算法获得静止状态下的脉率;Step 3, using the photoelectric sensor to extract the pulse rate to obtain the photoplethysmography PPG signal, and obtain the pulse rate in the static state by the interval algorithm of the adjacent R points in the PPG signal; 步骤四、根据运动分类算法将运动状态分为规律运动和非规律运动,其确定方法如下;Step 4, according to the motion classification algorithm, the motion state is divided into regular motion and irregular motion, and its determination method is as follows; S41、从加速度模值中选出最大值amax S41, select the maximum value a max from the acceleration modulus S42、对数组d中的加速度数据进行归一化处理;S42. Normalize the acceleration data in the array d; S43、选取归一化后的数组d中的最大值pmaxS43. Select the maximum value p max in the normalized array d; S44、设定阈值P,以pmax为基准值,在归一化后的数组d中遍历基准值左侧的其他数值,获取首个超过阈值P的数值,定义为pb;同样地,遍历基准值右侧的其他数值,获取首个超过阈值P的数值,定义为pc;当遍历完所有数值后,一旦有一侧没有获取到数值,则判断为非规律运动状态;执行步骤六;若pb和pc均获得,则执行步骤S45;S44. Set the threshold value P, take p max as the reference value, traverse other values on the left side of the reference value in the normalized array d, obtain the first value exceeding the threshold value P, and define it as p b ; similarly, traverse For other values on the right side of the reference value, obtain the first value that exceeds the threshold P, which is defined as p c ; after traversing all the values, once the value is not obtained on one side, it is judged as an irregular motion state; perform step 6; if Both p b and p c are obtained, then perform step S45; S45、定义pb与pmax之间的距离为L1,pmax与pc之间的距离为L2;若L1和L2相等,则说明做规律运动,执行步骤五;否则,判定为非规律运动,执行步骤六;S45. Define the distance between p b and p max as L 1 , and the distance between p max and p c as L 2 ; if L 1 and L 2 are equal, it means that you are doing regular exercise and go to step 5; otherwise, determine For irregular exercise, perform step 6; 步骤五、对于规律运动状态下的脉率提取方法,采用光电传感器采集PPG信号,并提取脉率;Step 5. For the pulse rate extraction method under the regular exercise state, the photoelectric sensor is used to collect the PPG signal, and the pulse rate is extracted; 步骤六、对于非规律运动状态下的脉率提取方法:Step 6. For the pulse rate extraction method in the state of irregular exercise: 先采用光电传感器采集PPG信号,之后采用ANC算法对PPG信号进行去噪处理,通过提取PPG信号中相邻P点的间期,获得脉率。First, the photoelectric sensor is used to collect the PPG signal, and then the ANC algorithm is used to denoise the PPG signal, and the pulse rate is obtained by extracting the interval between adjacent P points in the PPG signal.
CN201610471150.XA 2016-06-24 2016-06-24 A kind of pulse frequency extracting method based on wrist wearable device Expired - Fee Related CN105997043B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610471150.XA CN105997043B (en) 2016-06-24 2016-06-24 A kind of pulse frequency extracting method based on wrist wearable device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610471150.XA CN105997043B (en) 2016-06-24 2016-06-24 A kind of pulse frequency extracting method based on wrist wearable device

Publications (2)

Publication Number Publication Date
CN105997043A CN105997043A (en) 2016-10-12
CN105997043B true CN105997043B (en) 2018-11-20

Family

ID=57083859

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610471150.XA Expired - Fee Related CN105997043B (en) 2016-06-24 2016-06-24 A kind of pulse frequency extracting method based on wrist wearable device

Country Status (1)

Country Link
CN (1) CN105997043B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106448051A (en) * 2016-11-25 2017-02-22 广东电网有限责任公司电力科学研究院 Wearable physiological sensation equipment applicable to high-altitude operation protection
CN107223036B (en) * 2017-04-10 2019-03-29 深圳市汇顶科技股份有限公司 Object wearing device, the method and device for adaptively filtering out motion artifacts
EP3406189B1 (en) * 2017-05-25 2024-05-29 Tata Consultancy Services Limited System and method for heart rate estimation
CN109044316B (en) * 2018-07-11 2021-02-02 歌尔科技有限公司 Pure dynamic heart rate signal calculation method and device and intelligent wearable equipment
CN109875541A (en) * 2018-12-28 2019-06-14 北京津发科技股份有限公司 Pulses measure method, pulse measurement device and storage medium
CN110801214A (en) * 2019-11-27 2020-02-18 青岛歌尔智能传感器有限公司 Heart rate real-time detection method and system
CN111616695B (en) * 2020-06-29 2022-09-13 歌尔科技有限公司 Heart rate acquisition method, device, system and medium
EP4197436A4 (en) * 2021-11-01 2024-01-24 Beijing Honor Device Co., Ltd. HEART RATE DETECTION METHOD AND ELECTRONIC DEVICE

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1531901A (en) * 2003-03-19 2004-09-29 ������������ʽ���� Information collection device and pulsometer
CN101039617A (en) * 2004-10-15 2007-09-19 普尔塞特拉瑟技术有限公司 Motion cancellation of optical input signals for physiological pulse measurement
CN104605827A (en) * 2015-01-15 2015-05-13 辛勤 Human body vigorous motion state judgment method and device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105530858B (en) * 2013-09-16 2022-06-07 皇家飞利浦有限公司 System and method for estimating cardiovascular health of a person

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1531901A (en) * 2003-03-19 2004-09-29 ������������ʽ���� Information collection device and pulsometer
CN101039617A (en) * 2004-10-15 2007-09-19 普尔塞特拉瑟技术有限公司 Motion cancellation of optical input signals for physiological pulse measurement
CN104605827A (en) * 2015-01-15 2015-05-13 辛勤 Human body vigorous motion state judgment method and device

Also Published As

Publication number Publication date
CN105997043A (en) 2016-10-12

Similar Documents

Publication Publication Date Title
CN105997043B (en) A kind of pulse frequency extracting method based on wrist wearable device
Nallathambi et al. Integrate and fire pulse train automaton for QRS detection
US20160120431A1 (en) Medical device having automated ecg feature extraction
CN109907752A (en) A kind of cardiac diagnosis and monitoring method and system of the interference of removal motion artifacts and ecg characteristics detection
CN105266772A (en) Measuring method for physiological parameter
Satija et al. A simple method for detection and classification of ECG noises for wearable ECG monitoring devices
CN107041743A (en) A kind of real-time R wave detecting methods of electrocardiosignal
Gupta et al. Higher order derivative-based integrated model for cuff-less blood pressure estimation and stratification using PPG signals
US20170027464A1 (en) Detecting Artifacts in a Signal
CN108937957A (en) Detection method, device and detection device
CN104644160A (en) Electrocardiogram pseudo-difference signal identification method and device
CN109330582A (en) Heart rate and its characteristic index detection method based on ECG Signal Analysis
Liu et al. A novel P-QRS-T wave localization method in ECG signals based on hybrid neural networks
Chin et al. Bayesian real-time QRS complex detector for healthcare system
Kumar et al. STSR: Spectro-temporal super-resolution analysis of a reference signal less photoplethysmogram for heart rate estimation during physical activity
Satija et al. Low-complexity detection and classification of ECG noises for automated ECG analysis system
CN109498022A (en) A kind of respiratory rate extracting method based on photoplethysmographic
KR20150009843A (en) System for classification electrocardiogram signals
CN115486849B (en) Electrocardiosignal quality assessment method, device and equipment
Wang et al. Analysis of pulse waveforms preprocessing
CN106175698B (en) Sleep cycle detection device in sleep state analysis
Cherif et al. Detection of artifacts on photoplethysmography signals using random distortion testing
Jin et al. A novel heart rate detection algorithm in ballistocardiogram based on wavelet transform
CN106419884B (en) A kind of rate calculation method and system based on wavelet analysis
TW201440725A (en) Denoising method and apparatus of pulse wave signal and pulse oximetry

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181120