CN108903929B - Heart rate detection correction method, device, storage medium and system - Google Patents
Heart rate detection correction method, device, storage medium and system Download PDFInfo
- 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
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
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 secondWhere 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):
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;
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;
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.
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)
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 |
US20240252052A1 (en) * | 2021-11-01 | 2024-08-01 | Beijing Honor Device Co., Ltd. | 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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103298396B (en) * | 2010-10-12 | 2015-03-18 | 学校法人福冈大学 | 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 |
-
2018
- 2018-03-30 CN CN201810297974.9A patent/CN108903929B/en active Active
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 | |
CN103385702B (en) | A kind of non-invasive blood pressure continuous detection apparatus and method | |
CN109414203A (en) | Online heart rate estimation based on optical measurement | |
KR101910982B1 (en) | Method and apparatus for eliminating motion artifact of biosignal using personalized biosignal pattern | |
CN105249925B (en) | A kind of traditional Chinese medical pulse manifestation collecting device and noise reduction system and noise-reduction method | |
CN105748051A (en) | Blood pressure measuring method and device | |
CN103800003A (en) | ECG detection method and detector | |
KR101536361B1 (en) | Method for analyzing ecg and ecg apparatus therefor | |
WO2016187835A1 (en) | Continuous blood pressure measurement method, apparatus and device | |
CN107249445A (en) | Estimate heart rate by tracking optical signalling frequency component | |
CN107273827A (en) | A kind of R wave of electrocardiosignal detection method and device | |
KR101696791B1 (en) | Pulmonary function test apparatus using chest impedance and thereof method | |
US20190192082A1 (en) | Apparatus for eliminating motion artifacts by using ppg signal and method thereof | |
GB2552455A (en) | Blood monitoring | |
CN101897578B (en) | Method for segmenting arterial pressure signal by beats | |
CN104873184A (en) | Method and system for measuring heart rate | |
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 | |
CN116369888B (en) | Non-contact heart rate variability data acquisition method and device | |
KR102059685B1 (en) | Apparatus and method for estimating heart beat using photo-plethysmogram | |
CN110801214A (en) | Heart rate real-time detection method and system | |
KR101992763B1 (en) | Device and method of abnormal electrocardiographic signals information output | |
CN106504235B (en) | Rhythm of the heart method based on image procossing |
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 |