CN108903929B - Heart rate detection correction method, device, storage medium and system - Google Patents

Heart rate detection correction method, device, storage medium and system Download PDF

Info

Publication number
CN108903929B
CN108903929B CN201810297974.9A CN201810297974A CN108903929B CN 108903929 B CN108903929 B CN 108903929B CN 201810297974 A CN201810297974 A CN 201810297974A CN 108903929 B CN108903929 B CN 108903929B
Authority
CN
China
Prior art keywords
heart rate
exercise
intensity
exercise intensity
user
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.)
Active
Application number
CN201810297974.9A
Other languages
Chinese (zh)
Other versions
CN108903929A (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.)
Guangdong Spark Technology Co ltd
Original Assignee
Guangdong Spark Technology Co 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 Guangdong Spark Technology Co ltd filed Critical Guangdong Spark Technology Co ltd
Priority to CN201810297974.9A priority Critical patent/CN108903929B/en
Publication of CN108903929A publication Critical patent/CN108903929A/en
Application granted granted Critical
Publication of CN108903929B publication Critical patent/CN108903929B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Physiology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Cardiology (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (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 relates to the field of heart rate retrieval, in particular to a heart rate detection correction system, wherein a processor and a computer readable storage medium are arranged on the system, a computer program on the computer readable storage medium can be executed by the processor, the current exercise intensity of a user and a PPG signal are subjected to spectrum difference operation, and a spectrum peak tracking mechanism is used for positioning a tested heart rate HR 2; determining an estimated heart rate HR1 through a pre-stored corresponding relation curve; and then the estimated heart rate HR1 and the tested heart rate HR2 are weighted and corrected, so that more accurate corrected heart rate HR3 is obtained, and the problem of inaccurate heart rate estimation caused by irregular motion is solved.

Description

Heart rate detection correction method, device, storage medium and system
Technical Field
The invention relates to the field of heart rate retrieval, in particular to a heart rate detection correction system which is provided with a processor and a computer readable storage medium, wherein a computer program on the computer readable storage medium can be executed by the processor.
Background
The heart rate detection function is common on intelligent wearable products, and has very important significance for sports people or heart health disease patients. The product is dressed to intelligence on the market at present, and most all adopt PPG blood volume method to measure, and PPG signal detection device shines people's skin surface through emission green glow or infrared light, and according to the different absorption degree of blood to light, the reflected light intensity that photoelectric sensor received is also different to according to this waveform, the analysis obtains the rhythm of the heart.
However, relying on PPG signals alone for heart rate analysis under severe exercise is far from adequate and the algorithm accuracy can be poor. Therefore, a method for assisting in eliminating motion interference by using an acceleration sensor in the prior art appears, and specifically includes that a PPG signal detection device and the acceleration sensor acquire a PPG signal and an acceleration signal of a user, then difference operation is performed on the two signals on a frequency spectrum to remove a signal generated by motion noise, and then a heart rate frequency point is accurately positioned according to a spectral peak tracking mechanism. However, the frequency domain transformation interference elimination method is only suitable for comparison of regular motion, because only under the condition of regular motion, the acceleration signal and the heart rate can present ideal synchronism on the frequency spectrum, so that the interference can be eliminated better by subtraction, and if the motion is irregular, because the synchronism of the two signals on the frequency spectrum is poor, the subtraction cannot be completely eliminated, and at the moment, the positioning by using a spectral peak tracking mechanism can have deviation.
Disclosure of Invention
It is an object of the invention to improve the accuracy of heart rate detection.
To this end, a heart rate detection and correction system is provided, which includes an acceleration sensor, a PPG signal detection device, a computer readable storage medium and a processor, wherein the acceleration sensor and the PPG signal detection device respectively collect an acceleration signal and a PPG signal of a user to the processor, the processor runs a program in the computer readable storage medium according to the acceleration signal and the PPG signal, and when the program is run, the following determination steps are performed, including:
and a frequency domain transformation step: respectively carrying out frequency domain transformation on the current motion intensity of the user and the PPG signal;
the testing steps are as follows: performing difference operation on the frequency spectrum of the current motion intensity and the frequency spectrum of the current PPG signal, and taking the frequency corresponding to the highest peak in the frequency spectrum obtained by operation as the tested heart rate HR2 of the user;
pre-estimating: determining an estimated heart rate HR1 according to a pre-stored corresponding relation between the exercise intensity and the heart rate and the current exercise intensity of the user;
a correction step performed after the test step: the estimated heart rate HR1 and the test heart rate HR2 are weighted and summed to obtain a corrected heart rate HR 3.
In the estimating step, the manner of obtaining the corresponding relationship is specifically as follows: and testing the exercise intensity and the heart rate of different people in different exercise states, and fitting by using all the exercise intensities and the heart rates to obtain the corresponding relation between the exercise intensity and the heart rate.
In the estimation step, fitting is performed by using a least square method.
The exercise intensity test mode specifically comprises the following steps: acquiring n-axis data h of user during movement by using n-axis acceleration sensor1(i)、h2(i)……hn(i) Calculating the exercise intensity in the ith second
Figure BDA0001617126010000021
Where F is the sampling rate.
Has the advantages that:
according to the method, the current motion intensity of a user and the PPG signal are subjected to spectrum difference value operation, and a spectrum peak tracking mechanism is used for positioning a test heart rate HR 2; determining an estimated heart rate HR1 through a pre-stored corresponding relation curve; and then the estimated heart rate HR1 and the tested heart rate HR2 are weighted and corrected, so that more accurate corrected heart rate HR3 is obtained, and the problem of inaccurate heart rate estimation caused by irregular motion is solved.
Drawings
The invention is further described by using the drawings, but the embodiments in the drawings do not constitute any limitation to the invention, and for a user of ordinary skill in the art, other drawings can be obtained according to the following drawings without creative efforts.
Fig. 1 is a flow chart of a heart rate detection correction method of the present invention.
Fig. 2 is a graph of discrete data characterizing the relationship between heart rate and exercise intensity.
Fig. 3 is a graph of the correspondence relationship obtained by fitting the discrete points of fig. 2.
Fig. 4 is a graph of the motion intensity of a three-axis acceleration sensor in three directions.
Fig. 5 is a motion intensity diagram obtained by vector modulo of three directions by the triaxial acceleration sensor.
Detailed Description
The invention will be further described with reference to the following examples.
Fig. 1 is a flowchart of the modification of the present embodiment.
In advance, the processor firstly counts the exercise intensity and heart rate of different people in different exercise states (including regular exercise and irregular exercise), and takes all the exercise intensity and heart rate as samples to form a coordinate axis which takes an abscissa as a heart rate value and an ordinate as the exercise intensity as shown in fig. 2. As can be seen from fig. 2, although the exercise intensity of each person is different, the relationship between the heart rate and the exercise intensity is still largely the same as the relationship between the heart rate and the exercise intensity. Therefore, a least square method is used to perform polynomial fitting on the discrete points in fig. 2, thereby obtaining the corresponding relationship curve in fig. 3. Since the least square method has been applied for many years, the implementation principle and process thereof are not expanded in this embodiment to avoid repeated description, but it should be noted that the least square method is used for curve fitting to minimize the sum of squares of the errors between the estimated data and the actual data, so as to improve the accuracy.
In the statistical process, the processor mainly acquires the heart rate of the user to the processor through a PPG signal detection device; regarding the exercise intensity, the processor mainly obtains the exercise intensity through the built-in three-axis acceleration sensors, specifically, x (i), y (i), z (i) represent the three-axis acceleration data of the three-axis acceleration sensors at the same time, where the sampling rate is 20Hz, and the exercise intensity s (i) is calculated every 1s, for example, the exercise intensity in the ith second is expressed by the following formula (1):
Figure BDA0001617126010000031
the three-axis acceleration sensor detects the motion intensity of three-dimensional vertical spaces of X/Y/Z, the three detection directions are three, and the three-axis acceleration sensor cannot ensure that the three directions are always in the maximum amplitude due to different wearing modes and motion amplitudes, for example, the X-axis data in fig. 4 is obviously smaller, so that the finally obtained motion intensity has a certain amplitude, the three-axis acceleration data is required to be substituted into the formula (1) to be combined in a mode of simulating space vector modulo, the curve waveform with larger signal intensity shown in fig. 5 is obtained, and further the required motion intensity is obtained.
After obtaining the corresponding relation curve, the processor stores the corresponding relation curve in the database.
When a user moves, the processor acquires the current movement intensity of the user through the three-axis acceleration sensor, acquires the current PPG signal of the user through the PPG signal detection device, performs FFT (fast Fourier transform) on the current movement intensity and the current PPG signal respectively, performs difference operation on the frequency spectrum of the current movement intensity and the frequency spectrum of the current PPG signal, and takes the frequency corresponding to the highest peak in the frequency spectrum after operation as the test heart rate HR2 of the user. Then, the processor determines the estimated heart rate HR1 corresponding to the current exercise intensity according to the corresponding relation curve, and then performs weighted correction on the estimated heart rate HR1 and the tested heart rate HR2 by using the following formula (2), so that a corrected heart rate HR3 more fitting the reality is obtained, and the problem of inaccurate heart rate estimation caused by irregular exercise is solved.
HR3＝a×HR1+(1-a)×HR2 (2)
In formula (2), a is a variable factor, which may be actually corrected according to different PPG signal detection devices.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting the protection scope of the present invention, although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions can be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention.

Claims (8)

1. A method of heart rate detection modification, comprising:
and a frequency domain transformation step: respectively carrying out frequency domain transformation on the current motion intensity of the user and the PPG signal;
the testing steps are as follows: performing difference operation on the frequency spectrum of the current motion intensity and the frequency spectrum of the current PPG signal, and taking the frequency corresponding to the highest peak in the frequency spectrum obtained by operation as the tested heart rate HR2 of the user;
the method is characterized by further comprising the following steps:
pre-estimating: determining an estimated heart rate HR1 according to a pre-stored corresponding relation between the exercise intensity and the heart rate and the current exercise intensity of the user;
a correction step performed after the test step: carrying out weighted sum operation on the estimated heart rate HR1 and the tested heart rate HR2 to obtain a corrected heart rate HR 3;
the test mode of the exercise intensity is specifically that an n-axis acceleration sensor is used for collecting n-axis data h of a user during exercise1(i)、h2(i)……hn(i) Calculating the exercise intensity in the ith second
Figure FDA0002618921990000011
Where F is the sampling rate.
2. The method for detecting and correcting heart rate according to claim 1, wherein in the estimating step, the corresponding relationship is obtained by: and testing the exercise intensity and the heart rate of different people in different exercise states, and fitting by using all the exercise intensities and the heart rates to obtain the corresponding relation between the exercise intensity and the heart rate.
3. A method of heart rate detection modification as claimed in claim 2, wherein in the estimating step, the fitting is performed by a least squares method.
4. An apparatus for heart rate detection modification, comprising:
the frequency domain transformation device: respectively carrying out frequency domain transformation on the current motion intensity of the user and the PPG signal;
the testing device comprises: performing difference operation on the frequency spectrum of the current motion intensity and the frequency spectrum of the current PPG signal, and taking the frequency corresponding to the highest peak in the frequency spectrum obtained by operation as the tested heart rate HR2 of the user;
it is characterized by also comprising the following devices:
the estimation device: determining an estimated heart rate HR1 according to a pre-stored estimated corresponding relation between the exercise intensity and the heart rate and the current exercise intensity of the user;
correction means executed after the test means: carrying out weighted sum operation on the estimated heart rate HR1 and the tested heart rate HR2 to obtain a corrected heart rate HR 3;
the test mode of the exercise intensity is specifically that an n-axis acceleration sensor is used for collecting n-axis data h of a user during exercise1(i)、h2(i)……hn(i) Calculating the exercise intensity in the ith second
Figure FDA0002618921990000021
Where F is the sampling rate.
5. The apparatus for detecting and correcting heart rate as claimed in claim 4, wherein the estimating means obtains the corresponding relationship by: and testing the exercise intensity and the heart rate of different people in different exercise states, and fitting by using all the exercise intensities and the heart rates to obtain the corresponding relation between the exercise intensity and the heart rate.
6. The apparatus for heart rate detection modification as claimed in claim 5, wherein the estimating means performs the fitting using a least squares method.
7. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the method of heart rate detection correction according to any one of claims 1-3.
8. Heart rate detects correction system, including acceleration sensor, PPG signal detection device and treater, acceleration sensor and PPG signal detection device, characterized by are connected to the treater electricity respectively: further comprising a computer readable storage medium as claimed in claim 7, the program in the medium being executable on the processor.
CN201810297974.9A 2018-03-30 2018-03-30 Heart rate detection correction method, device, storage medium and system Active CN108903929B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810297974.9A CN108903929B (en) 2018-03-30 2018-03-30 Heart rate detection correction method, device, storage medium and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810297974.9A CN108903929B (en) 2018-03-30 2018-03-30 Heart rate detection correction method, device, storage medium and system

Publications (2)

Publication Number Publication Date
CN108903929A CN108903929A (en) 2018-11-30
CN108903929B true CN108903929B (en) 2021-02-26

Family

ID=64402915

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810297974.9A Active CN108903929B (en) 2018-03-30 2018-03-30 Heart rate detection correction method, device, storage medium and system

Country Status (1)

Country Link
CN (1) CN108903929B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110540118B (en) * 2019-06-21 2021-04-09 浙江大学 Intelligent detection device for detecting motion state of car
CN110575153B (en) * 2019-09-27 2022-06-10 歌尔股份有限公司 Heart rate detection method and intelligent wearable device
CN110916639B (en) * 2019-12-23 2022-09-02 深圳市圆周率智能信息科技有限公司 Method, system, wearable device and computer-readable storage medium for acquiring exercise heart rate recovery rate
CN111728603A (en) * 2020-01-09 2020-10-02 成都维客昕微电子有限公司 Sampling rate self-adjusting method for optical heart rate sensor
CN112790752B (en) * 2021-01-22 2022-09-27 维沃移动通信有限公司 Heart rate value correction method and device and electronic equipment
WO2023071501A1 (en) * 2021-11-01 2023-05-04 北京荣耀终端有限公司 Heart rate detection method and electronic device
CN114521880B (en) * 2022-01-21 2023-09-01 中国人民解放军陆军军医大学 Method, system and computer storage medium for calculating heart rate under exercise state

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130231576A1 (en) * 2010-10-12 2013-09-05 Fukuoka University Somatic data-measuring apparatus and somatic data measurement method
CN105286842B (en) * 2015-11-06 2018-04-03 深圳风景网络科技有限公司 A kind of method and device based on acceleration transducer predicted motion process heart rate
US20170164847A1 (en) * 2015-12-15 2017-06-15 Texas Instruments Incorporated Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals
JP2017136164A (en) * 2016-02-02 2017-08-10 富士通株式会社 Sensor information processing device, sensor unit, and sensor information processing program
JP6642055B2 (en) * 2016-02-02 2020-02-05 富士通株式会社 Sensor information processing device, sensor unit, and sensor information processing program
CN105943015A (en) * 2016-06-04 2016-09-21 浙江大学 Wearable heart rate variability (HRV) monitoring device with active noise reduction function
CN106264505A (en) * 2016-07-21 2017-01-04 浙江师范大学 A kind of heart rate spectral peak system of selection based on support vector machine

Also Published As

Publication number Publication date
CN108903929A (en) 2018-11-30

Similar Documents

Publication Publication Date Title
CN108903929B (en) Heart rate detection correction method, device, storage medium and system
CN109875529B (en) Vital sign detection method and system based on ultra-wideband radar
CN105249925B (en) A kind of traditional Chinese medical pulse manifestation collecting device and noise reduction system and noise-reduction method
CN111065323A (en) Device and method for determining the blood pressure of a subject
CN105748051A (en) Blood pressure measuring method and device
CN103800003A (en) ECG detection method and detector
KR20130092849A (en) Method and apparatus for eliminating motion artifact of biosignal using personalized biosignal pattern
WO2016187835A1 (en) Continuous blood pressure measurement method, apparatus and device
CN107273827B (en) Electrocardiosignal R wave detection method and device
CN108992053A (en) A method of real-time chainless detection heart rate and eartbeat interval
CN107249445A (en) Estimate heart rate by tracking optical signalling frequency component
CN109330582A (en) Heart rate and its characteristic index detection method based on ECG Signal Analysis
US20130190635A1 (en) Calibration method of electrocardiogram signals and the application program for the same
CN105105739B (en) Short-distance wireless heart rate and heart rate variability detection method
GB2552455A (en) Blood monitoring
CN111685760B (en) Human body respiratory frequency calculation method based on radar measurement
JP4122003B2 (en) Extraction method and device for biological signals such as heartbeat and respiration
US11457823B2 (en) Wearable blood pressure detecting device and detecting method thereof
KR101696791B1 (en) Pulmonary function test apparatus using chest impedance and thereof method
KR102059685B1 (en) Apparatus and method for estimating heart beat using photo-plethysmogram
CN110801214A (en) Heart rate real-time detection method and system
CN106504235B (en) Rhythm of the heart method based on image procossing
KR20230077607A (en) Apparatus and method for determining a distance for measuring heartbeat based on temporal phase coherency
CN112617786A (en) Heart rate detection device and method based on tof camera
KR101992763B1 (en) Device and method of abnormal electrocardiographic signals information output

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant