CN105286845A - Movement noise elimination method suitable for wearable heart rate measurement device - Google Patents

Movement noise elimination method suitable for wearable heart rate measurement device Download PDF

Info

Publication number
CN105286845A
CN105286845A CN201510873977.9A CN201510873977A CN105286845A CN 105286845 A CN105286845 A CN 105286845A CN 201510873977 A CN201510873977 A CN 201510873977A CN 105286845 A CN105286845 A CN 105286845A
Authority
CN
China
Prior art keywords
spectrum
heart rate
signal
motion
motion artifacts
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.)
Pending
Application number
CN201510873977.9A
Other languages
Chinese (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.)
Zhejiang Normal University CJNU
Original Assignee
Zhejiang Normal University CJNU
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 Zhejiang Normal University CJNU filed Critical Zhejiang Normal University CJNU
Priority to CN201510873977.9A priority Critical patent/CN105286845A/en
Publication of CN105286845A publication Critical patent/CN105286845A/en
Pending legal-status Critical Current

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

The invention discloses a movement noise elimination method suitable for a wearable heart rate measurement device. The method aims at improving the heart rate measurement accuracy of the wearable heart rate measurement device. In the method, the wearable heart rate measurement device collects multiple pulse oximeter signals and movement acceleration signals, during the same time period, of a user; the spectrum peak positions of movement noise in multiple pulse oximeter signal frequency spectrums are flush with the spectrum peak positions of movement acceleration signal frequency spectrums, and the pulse oximeter signal frequency spectrums can be removed through the spectral subtraction method; finally, the heart rate frequency point positions are accurately positioned according to the spectrum peak tracking mechanism. By means of the method, movement noise in heart rates is effectively eliminated, the non-peak situation, the multi-peak situation and the target spectrum peak tracking loss situation which happen to the pulse oximeter signal frequency spectrums after the spectral subtraction method is conducted are eliminated, and accurate measurement, based on the wearable device, of real-time heart rates is achieved.

Description

A kind of motion artifacts removing method being applicable to wearable heart rate measuring device
Technical field
The present invention relates to field of information processing, particularly relate to a kind of motion artifacts removing method being applicable to wearable heart rate measuring device.
Background technology
Heart rate refers to human heart number of times of beating per minute, is as the criterion with first sound.In human parameters detects, heart rate is an important physical signs, for medical diagnosis provides reference.Meanwhile, heart rate also can be used as the objective evaluation index of human motion physiological stress.Along with the rise of the wearable smart machines such as intelligent watch, Intelligent spire lamella, Intelligent bracelet, and people are for the attention of health status, based on photoplethysmographic signal, the method that heart rate is monitored received to the extensive concern of industrial quarters and academia.
Due to normal containing motion artifacts in photoplethysmographic signal, be difficult to Measurement accuracy heart rate, need corresponding noise-removed technology.The technology of existing many removal motion artifacts is suggested at present, and such as Independent Component Analysis (ICA), Wavelet noise-eliminating method, adaptive-filtering denoising method (ANC), classical mode decomposition method (EMD) etc. are all widely used.But above-mentioned algorithm is mainly for relaxing or inviolent motion, and such as hands moves, walk, jog (speed is lower than 8km/h).When motion artifacts is very strong, the effect of above-mentioned algorithm is then unsatisfactory.
The present invention proposes a kind of removing method for strong movements noise.In the method, in multiple photoplethysmographic signal frequency spectrum, the frequency location of motion artifacts aligns with the frequency location of acceleration of motion signal spectrum, utilize spectrum-subtraction easily can deduct motion artifacts spectrum peak from multiple original photoplethysmographic signal frequency spectrum, obtain multiple clean photoplethysmographic signal frequency spectrum.Meanwhile, this method propose spectrum peak follow-up mechanism, can process occur in multiple photoplethysmographic signal frequency spectrum after spectrum-subtraction compose peak by with the situation of losing without peak, multimodal and target.Invention effectively eliminates the motion artifacts in heart rate, achieve the Measurement accuracy based on the real-time heart rate of wearable device and calculating.
Summary of the invention
Technical problem to be solved by this invention be how when motion artifacts strongly provide a kind of method of effective removal motion artifacts, to obtain real-time heart rate value accurately.
In order to solve the problems of the technologies described above, the invention provides a kind of motion artifacts removing method being applicable to wearable heart rate measuring device, comprise spectrum-subtraction and spectrum peak follow-up mechanism two parts, it is characterized in that:
Described wearable heart rate measuring device gathers with the multiple photoplethysmographic signal in the time period and acceleration of motion signal at user's wrist place; Then, described spectrum-subtraction is utilized to remove motion artifacts in multiple photoplethysmographic signal; Finally, heart rate frequency point position is accurately located according to described spectrum peak follow-up mechanism.
The method comprises the steps:
Described wearable heart rate measuring device gathers the multiple photoplethysmographic signal of user at one time in section and acceleration of motion signal; Down-sampling process is carried out to above-mentioned multiple photoplethysmographic signal and acceleration of motion signal; Then the above-mentioned signal after down-sampling is carried out bandpass filtering operation.
Described spectrum-subtraction removes described motion artifacts signal effectively according to the strong correlation performance of motion artifacts signal in described acceleration of motion signal and described multiple photoplethysmographic signal, obtains multiple pure photoplethysmographic signal frequency spectrum; Each sub stage of described spectrum peak follow-up mechanism processes above-mentioned multiple pure photoplethysmographic signal frequency spectrum, the heart rate frequency point position of consumer positioning.
Preferably, the embedded multiple photoplethysmographic sensor of described wearable heart rate measuring device and three axis accelerometer; Multiple photoplethysmographic signal of described multiple photoplethysmographic sensor acquisition user; Described three axis accelerometer gather user at the same time between acceleration of motion signal in section.
Preferably, in described multiple photoplethysmographic signal, motion artifacts has more identical frequency with synchronous described acceleration of motion signal, the frequency location of motion artifacts in described multiple photoplethysmographic signal frequency spectrum is alignd with the frequency location of described acceleration of motion signal spectrum, utilize described spectrum-subtraction easily can deduct the spectrum peak of motion artifacts from described multiple photoplethysmographic signal frequency spectrum, obtain multiple clean photoplethysmographic signal frequency spectrum.
Preferably, for guaranteeing that described spectrum-subtraction is effective, described multiple photoplethysmographic signal frequency spectrum and described acceleration of motion signal spectrum needed to be operated by energy normalized before carrying out described spectrum-subtraction.
Preferably, described spectrum peak follow-up mechanism is mainly based on following two principles: the first, and in described multiple photoplethysmographic signal frequency spectrum, maximum spectrum peak position is corresponding with heart rate spectrum peak position in most cases; The second, the continuous window center rate overlapping in major part is very close; Described spectrum peak follow-up mechanism is selected and composes peak to find that three sub stages form by initializing, composing peak.
Compared with prior art, technical scheme provided by the invention effectively eliminates the motion artifacts in heart rate, and solve occur in multiple photoplethysmographic signal frequency spectrum after spectrum-subtraction without peak, multimodal and target spectrum peak by with the situation of losing.Thus improve the heart rate measurements precision of wearable heart rate measuring device, achieve the measurement and calculation of the real-time heart rate based on wearable device.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the motion artifacts removing method of the embodiment of the present invention;
Fig. 2 is the schematic flow sheet of the spectrum-subtraction of the embodiment of the present invention.
Detailed description of the invention
Describe embodiments of the present invention in detail below in conjunction with drawings and Examples, to the present invention, how application technology means solve technical problem whereby, and the implementation procedure reaching relevant art effect can fully understand and implement according to this.
In technical scheme of the present invention, in multiple photoplethysmographic signal frequency spectrum, the frequency location of motion artifacts aligns with the frequency location of acceleration of motion signal spectrum, utilize acceleration of motion signal spectrum easily can deduct the spectrum peak of motion artifacts from multiple original photoplethysmographic signal frequency spectrum, obtain multiple clean photoplethysmographic signal frequency spectrum.Meanwhile, this method propose spectrum peak follow-up mechanism, can process occur in multiple photoplethysmographic signal frequency spectrum after spectrum-subtraction compose peak by with the situation of losing without peak, multimodal and target.This technical scheme effectively eliminates the motion artifacts in heart rate, achieves the Measurement accuracy based on the real-time heart rate of wearable device and calculating.
The motion artifacts removing method of embodiment one, wearable heart rate measuring device
Fig. 1 is the schematic flow sheet of the motion artifacts removing method of the present embodiment, and Fig. 2 is the schematic flow sheet of the spectrum-subtraction of the present embodiment.
The present embodiment shown in Fig. 1, is the overall flow of the motion artifacts removing method of wearable heart rate measuring device, mainly comprises the steps:
Step S210, wearable heart rate measuring device utilizes two photoplethysmographic signal being distributed in photoplethysmographic sensor acquisition two passages of diverse location (hereinafter referred to as PPG 1and PPG 2), utilize three axis accelerometer collection with the acceleration of motion signal of the passage of three in the time period.
Step S220, the original sampling frequency of above-mentioned primary signal is 125Hz, and for reducing amount of calculation, needing to carry out down-sampling to sample frequency to above-mentioned primary signal is the operation of 25Hz.
Step S230, the second order Butterworth filter that the above-mentioned signal demand after down-sampling is 0.4Hz-4Hz by passband carries out filtering, to eliminate the interference of motion artifacts beyond certain frequency scope and other noise.
Step S240, in two photoplethysmographic signal frequency spectrums, the frequency location of motion artifacts aligns with the frequency location of acceleration of motion signal spectrum, utilizes spectrum-subtraction can obtain removing two photoplethysmographic signal frequency spectrums of motion artifacts, namely clean PPG 1, PPG 2signal spectrum.
In this step, typically, the concrete steps of spectrum-subtraction are as shown in Figure 2:
Step S310, for each Frequency point f i(i=1 ..., N), from the acceleration of motion signal spectrum of three passages, select maximum spectral coefficient, be defined as C i.
Step S320, PPG 1, PPG 2signal spectrum is at each Frequency point f i(i=1 ..., N) on spectral coefficient all deduct C i, PPG after above-mentioned process 1signal spectrum is at 0≤f iin≤199 scopes, spectral coefficient maximum is defined as p max1, PPG 2signal spectrum is at 0≤f iin≤199 scopes, spectral coefficient maximum is defined as p max2.
Step S330, PPG 1signal spectrum is at 0≤f iin≤199 scopes, spectral coefficient is less than p max1/ 4 be all set to 0, PPG 2signal spectrum is at 0≤f iin≤199 scopes, spectral coefficient is less than p max2/ 4 be all set to 0.
Step S340, obtains the clean photoplethysmographic signal frequency spectrum of two different passages after aforesaid operations, now starts spectrum peak follow-up mechanism, to reach the object of locating heart rate frequency point position exactly.
In order to better set forth described spectrum-subtraction, to following some be described:
The first, the Frequency point f of digital signal spectrum i(i=1 ..., N) from 0, the relation between itself and location index i as shown in formula (1),
f i=i-1(1)
The second, the Frequency point f of digital signal spectrum iwith the frequency f relation of analogue signal as shown in formula (2),
f = f i N f s = i - 1 N f s - - - ( 2 )
Wherein, f sfor sample frequency, N is sampled point;
3rd, the mankind have the maximum heart rate of record to be 230 beats/min, and the heart rate that (comprises strenuous exercise) is in most cases lower than 180 beats/min; F is set in the present embodiment s=25Hz, N=1024, therefore described spectrum-subtraction only analyzes 0≤f itwo photoplethysmographic signal frequency spectrums in≤199 scopes;
4th, effective in order to ensure spectrum-subtraction, should notice that two photoplethysmographic signal frequency spectrums and acceleration of motion signal spectrum needed to be operated by energy normalized before carrying out spectrum-subtraction process.
Step S250, obtains two clean photoplethysmographic signal frequency spectrums after aforesaid operations, the heart rate frequency point position of recycling spectrum peak follow-up mechanism consumer positioning.
In this step, typically, compose peak follow-up mechanism to be made up of three parts:
1) initialize: need user to reduce hand exercise in initial several seconds as far as possible, to ensure the accuracy of initial heart rate frequency point position, in the present embodiment, select PPG 1the heart rate spectrum peak position of maximum position, peak as correspondence is composed in signal spectrum;
2) compose peak to select: utilize the spectrum peak position in last time window corresponding to heart rate to go to find the spectrum peak corresponding to two photoplethysmographic signal spectral centroid rates of actual time window;
In actual applications, there will be some extreme cases.In latter two photoplethysmographic signal frequency spectrum of spectrum-subtraction, such as there is the situation without peak or multimodal.Wherein, two photoplethysmographic signal frequency spectrums are referred to all without the spectrum peak corresponding to heart rate or only have a photoplethysmographic signal frequency spectrum to contain spectrum peak corresponding to heart rate without peak; Multimodal refers to there is multiple spectrum peak near the spectrum peak position corresponding to two photoplethysmographic signal spectral centroid rates or only have near the spectrum peak position corresponding to a photoplethysmographic signal spectral centroid rate to there is multiple spectrum peak.
3) compose peak to find: this stage can effectively prevent multiple continuous time window from occurring unusual condition and causing spectrum peak corresponding to heart rate by with the situation of losing.
Compose peak follow-up mechanism in the present embodiment, can effectively process occur in latter two photoplethysmographic signal frequency spectrum of spectrum-subtraction without peak, multimodal and target spectrum peak by with the situation of losing.
Step S260, after above-mentioned steps process, can export the real-time heart rate value of user.
In the present embodiment, embedded two the photoplethysmographic sensors of wearable heart rate measuring device and three axis accelerometer, and gather with the photoplethysmographic signal of two in the time period and acceleration of motion signal at user's wrist place; In two photoplethysmographic signal frequency spectrums, the spectrum peak position of motion artifacts and the spectrum peak position of acceleration of motion signal spectrum align, and utilize spectrum-subtraction can obtain removing two photoplethysmographic signal frequency spectrums of motion artifacts; Finally, heart rate frequency point position can accurately be located according to spectrum peak follow-up mechanism.The method effectively eliminates the motion artifacts in heart rate, solve occur in latter two photoplethysmographic signal frequency spectrum of spectrum-subtraction without peak, multimodal and target spectrum peak by with the situation of losing.Thus improve the heart rate measurements precision of wearable heart rate measuring device, achieve the measurement and calculation of the real-time heart rate based on wearable device.
Although the embodiment disclosed by the present invention is as above, the embodiment that foregoing just adopts for the ease of understanding the present invention, and be not used to limit the present invention.Under the prerequisite of the spirit do not departed from disclosed by the present invention and scope, any modification and change can be done what implement in form and in details, but scope of patent protection of the present invention, the scope that still must define with appending claims is as the criterion.

Claims (5)

1. be applicable to a motion artifacts removing method for wearable heart rate measuring device, comprise spectrum-subtraction and spectrum peak follow-up mechanism two parts, it is characterized in that:
Described wearable heart rate measuring device gathers with the multiple photoplethysmographic signal in the time period and acceleration of motion signal at user's wrist place; Then, described spectrum-subtraction is utilized to remove motion artifacts in multiple photoplethysmographic signal; Finally, heart rate frequency point position is accurately located according to described spectrum peak follow-up mechanism;
The method comprises the steps:
Described wearable heart rate measuring device gathers the multiple photoplethysmographic signal of user at one time in section and acceleration of motion signal; Down-sampling process is carried out to above-mentioned multiple photoplethysmographic signal and acceleration of motion signal; Then the above-mentioned signal after down-sampling is carried out bandpass filtering operation;
Described spectrum-subtraction removes described motion artifacts signal effectively according to the strong correlation performance of motion artifacts signal in described acceleration of motion signal and described multiple photoplethysmographic signal, obtains multiple pure photoplethysmographic signal frequency spectrum; Each sub stage of described spectrum peak follow-up mechanism processes above-mentioned multiple pure photoplethysmographic signal frequency spectrum, the heart rate frequency point position of consumer positioning.
2. the motion artifacts removing method being applicable to wearable heart rate measuring device according to claim 1, is characterized in that:
The embedded multiple photoplethysmographic sensor of described wearable heart rate measuring device and three axis accelerometer; Multiple photoplethysmographic signal of described multiple photoplethysmographic sensor acquisition user; Described three axis accelerometer gather user at the same time between acceleration of motion signal in section.
3. the motion artifacts removing method being applicable to wearable heart rate measuring device according to claim 1, is characterized in that:
In described multiple photoplethysmographic signal, motion artifacts has more identical frequency with synchronous described acceleration of motion signal, the frequency location of motion artifacts in described multiple photoplethysmographic signal frequency spectrum is alignd with the frequency location of described acceleration of motion signal spectrum, utilize described spectrum-subtraction easily can deduct the spectrum peak of motion artifacts from described multiple photoplethysmographic signal frequency spectrum, obtain multiple clean photoplethysmographic signal frequency spectrum.
4. the motion artifacts removing method being applicable to wearable heart rate measuring device according to claim 3, is characterized in that:
For guaranteeing that described spectrum-subtraction is effective, described multiple photoplethysmographic signal frequency spectrum and described acceleration of motion signal spectrum needed to be operated by energy normalized before carrying out described spectrum-subtraction.
5. the motion artifacts removing method being applicable to wearable heart rate measuring device according to claim 1, is characterized in that:
Described spectrum peak follow-up mechanism is mainly based on following two principles: the first, and in described multiple photoplethysmographic signal frequency spectrum, maximum spectrum peak position is corresponding with heart rate spectrum peak position in most cases; The second, the continuous window center rate overlapping in major part is very close; Described spectrum peak follow-up mechanism is selected and composes peak to find that three sub stages form by initializing, composing peak.
CN201510873977.9A 2015-11-29 2015-11-29 Movement noise elimination method suitable for wearable heart rate measurement device Pending CN105286845A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510873977.9A CN105286845A (en) 2015-11-29 2015-11-29 Movement noise elimination method suitable for wearable heart rate measurement device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510873977.9A CN105286845A (en) 2015-11-29 2015-11-29 Movement noise elimination method suitable for wearable heart rate measurement device

Publications (1)

Publication Number Publication Date
CN105286845A true CN105286845A (en) 2016-02-03

Family

ID=55185177

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510873977.9A Pending CN105286845A (en) 2015-11-29 2015-11-29 Movement noise elimination method suitable for wearable heart rate measurement device

Country Status (1)

Country Link
CN (1) CN105286845A (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105816165A (en) * 2016-05-04 2016-08-03 上海轶卓信息科技有限公司 Real-time dynamic heart rate monitoring device and monitoring method
CN105877706A (en) * 2016-03-31 2016-08-24 济南大学 Heart-sound enhancement method based on improved spectral subtraction
CN105919584A (en) * 2016-06-23 2016-09-07 电子科技大学 Heart rate estimation method and device for wearable heart rate monitoring equipment
CN105943012A (en) * 2016-04-30 2016-09-21 浙江师范大学 Heart rate measurement method capable of removing motion noise in photoelectric plethysmography signals
CN106137167A (en) * 2016-07-21 2016-11-23 浙江师范大学 A kind of motion artifacts detection method based on photoplethysmographic signal
CN106264505A (en) * 2016-07-21 2017-01-04 浙江师范大学 A kind of heart rate spectral peak system of selection based on support vector machine
CN107223036A (en) * 2017-04-10 2017-09-29 深圳市汇顶科技股份有限公司 Object wearing device, the adaptive method and device for filtering out motion artifacts
WO2017202120A1 (en) * 2016-05-26 2017-11-30 华为技术有限公司 Method and device for collecting ppg signal
CN108294737A (en) * 2018-01-26 2018-07-20 深圳还是威健康科技有限公司 Method for measuring heart rate, device and intelligent wearable device
CN108478206A (en) * 2018-02-02 2018-09-04 北京邮电大学 Rhythm of the heart method based on pulse wave under motion state
CN108926333A (en) * 2017-05-25 2018-12-04 塔塔咨询服务有限公司 System and method for heart rate estimation
CN109222948A (en) * 2018-09-17 2019-01-18 歌尔科技有限公司 Eliminate method, apparatus, electronic equipment and the storage medium of motion artifacts noise
CN109222949A (en) * 2018-10-12 2019-01-18 杭州士兰微电子股份有限公司 Heart rate detection method and heartbeat detection device
WO2019014931A1 (en) * 2017-07-21 2019-01-24 深圳市汇顶科技股份有限公司 Interference analysis method and apparatus for biological signal, and wearable device
CN109528163A (en) * 2018-11-15 2019-03-29 深圳市国通世纪科技开发有限公司 A kind of method and apparatus of sleep monitor
WO2020042219A1 (en) * 2018-08-31 2020-03-05 深圳大学 Novel identity authentication method for wearable apparatus
CN110916639A (en) * 2019-12-23 2020-03-27 深圳市圆周率智能信息科技有限公司 Method, system, wearable device and computer-readable storage medium for acquiring exercise heart rate recovery rate
CN111643052A (en) * 2019-04-03 2020-09-11 上海铼锶信息技术有限公司 Method and system for reducing motion artifacts in pulse wave signals
CN111714110A (en) * 2020-05-19 2020-09-29 成都云卫康医疗科技有限公司 Real-time heart rate calculation method based on PPG waveform
CN112494001A (en) * 2020-11-25 2021-03-16 青岛歌尔智能传感器有限公司 PPG signal quality evaluation method and device, wearable device and storage medium
CN113271845A (en) * 2019-01-10 2021-08-17 株式会社电装 Biological information detection device

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5400409A (en) * 1992-12-23 1995-03-21 Daimler-Benz Ag Noise-reduction method for noise-affected voice channels
CN101972143A (en) * 2010-11-16 2011-02-16 电子科技大学 Blind source extraction-based atrial fibrillation monitoring method
CN104161505A (en) * 2014-08-13 2014-11-26 北京邮电大学 Motion noise interference eliminating method suitable for wearable heart rate monitoring device
CN104367310A (en) * 2013-08-14 2015-02-25 深圳市西微数字技术有限公司 Wearable heart rate detection device
TW201511735A (en) * 2013-08-06 2015-04-01 Laurence Richard Olivier A PPG-based physiological sensing system with a spatio-temporal sampling approach towards identifying and removing motion artifacts from optical signals
WO2015084376A1 (en) * 2013-12-05 2015-06-11 Apple Inc. Wearable multi-modal physiological sensing sysem
CN104739399A (en) * 2014-02-10 2015-07-01 北京金日吉通科贸有限公司 Method for detecting movement offset in pulse wave heart rate calculation
US20150196257A1 (en) * 2014-01-13 2015-07-16 The Board Of Regents, The University Of Texas System Systems and methods for physiological signal enhancement and biometric extraction using non-invasive optical sensors

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5400409A (en) * 1992-12-23 1995-03-21 Daimler-Benz Ag Noise-reduction method for noise-affected voice channels
CN101972143A (en) * 2010-11-16 2011-02-16 电子科技大学 Blind source extraction-based atrial fibrillation monitoring method
TW201511735A (en) * 2013-08-06 2015-04-01 Laurence Richard Olivier A PPG-based physiological sensing system with a spatio-temporal sampling approach towards identifying and removing motion artifacts from optical signals
CN104367310A (en) * 2013-08-14 2015-02-25 深圳市西微数字技术有限公司 Wearable heart rate detection device
WO2015084376A1 (en) * 2013-12-05 2015-06-11 Apple Inc. Wearable multi-modal physiological sensing sysem
US20150196257A1 (en) * 2014-01-13 2015-07-16 The Board Of Regents, The University Of Texas System Systems and methods for physiological signal enhancement and biometric extraction using non-invasive optical sensors
CN104739399A (en) * 2014-02-10 2015-07-01 北京金日吉通科贸有限公司 Method for detecting movement offset in pulse wave heart rate calculation
CN104161505A (en) * 2014-08-13 2014-11-26 北京邮电大学 Motion noise interference eliminating method suitable for wearable heart rate monitoring device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HAYATO FUKUSHIMA ET AL.: "Estimating Heart Rate using Wrist-type Photoplethysmography and Acceleration sensor while running", 《34TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE EMBS》 *
ZHILIN ZHANG ET AL.: "Photoplethysmography-BasedHeartRateMonitoringinPhysical ActivitiesviaJointSparseSpectrumReconstruction", 《IEEETRANSACTIONSONBIOMEDICALENGINEERING》 *
ZHILIN ZHANG ET AL.: "TROIKA: A General Framework for Heart Rate Monitoring Using Wrist-Type Photoplethysmographic Signals During Intensive Physical Exercise", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 *

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105877706A (en) * 2016-03-31 2016-08-24 济南大学 Heart-sound enhancement method based on improved spectral subtraction
CN105943012A (en) * 2016-04-30 2016-09-21 浙江师范大学 Heart rate measurement method capable of removing motion noise in photoelectric plethysmography signals
CN105816165B (en) * 2016-05-04 2019-12-13 把脉(上海)信息科技有限公司 real-time dynamic heart rate monitoring device and monitoring method
CN105816165A (en) * 2016-05-04 2016-08-03 上海轶卓信息科技有限公司 Real-time dynamic heart rate monitoring device and monitoring method
CN107432741B (en) * 2016-05-26 2019-07-09 华为终端有限公司 A kind of acquisition method and device of PPG signal
WO2017202120A1 (en) * 2016-05-26 2017-11-30 华为技术有限公司 Method and device for collecting ppg signal
CN107432741A (en) * 2016-05-26 2017-12-05 华为终端(东莞)有限公司 A kind of acquisition method and device of PPG signals
US11197620B2 (en) 2016-05-26 2021-12-14 Huawei Technologies Co., Ltd. PPG signal collection method and apparatus
CN105919584A (en) * 2016-06-23 2016-09-07 电子科技大学 Heart rate estimation method and device for wearable heart rate monitoring equipment
CN105919584B (en) * 2016-06-23 2018-10-16 电子科技大学 Heart rate method of estimation and device for wearable heart rate monitor apparatus
CN106137167A (en) * 2016-07-21 2016-11-23 浙江师范大学 A kind of motion artifacts detection method based on photoplethysmographic signal
CN106264505A (en) * 2016-07-21 2017-01-04 浙江师范大学 A kind of heart rate spectral peak system of selection based on support vector machine
CN107223036A (en) * 2017-04-10 2017-09-29 深圳市汇顶科技股份有限公司 Object wearing device, the adaptive method and device for filtering out motion artifacts
WO2018187895A1 (en) * 2017-04-10 2018-10-18 深圳市汇顶科技股份有限公司 Wearable device, method and device for self-adaptive filtering of motion interference
CN107223036B (en) * 2017-04-10 2019-03-29 深圳市汇顶科技股份有限公司 Object wearing device, the method and device for adaptively filtering out motion artifacts
CN108926333A (en) * 2017-05-25 2018-12-04 塔塔咨询服务有限公司 System and method for heart rate estimation
WO2019014931A1 (en) * 2017-07-21 2019-01-24 深圳市汇顶科技股份有限公司 Interference analysis method and apparatus for biological signal, and wearable device
CN108294737A (en) * 2018-01-26 2018-07-20 深圳还是威健康科技有限公司 Method for measuring heart rate, device and intelligent wearable device
CN108294737B (en) * 2018-01-26 2020-11-13 深圳市元征科技股份有限公司 Heart rate measuring method and device and intelligent wearable equipment
CN108478206B (en) * 2018-02-02 2021-08-13 北京邮电大学 Heart rate monitoring method based on pulse wave in motion state
CN108478206A (en) * 2018-02-02 2018-09-04 北京邮电大学 Rhythm of the heart method based on pulse wave under motion state
WO2020042219A1 (en) * 2018-08-31 2020-03-05 深圳大学 Novel identity authentication method for wearable apparatus
CN109222948A (en) * 2018-09-17 2019-01-18 歌尔科技有限公司 Eliminate method, apparatus, electronic equipment and the storage medium of motion artifacts noise
CN109222948B (en) * 2018-09-17 2021-07-13 歌尔科技有限公司 Method, device, electronic equipment and storage medium for eliminating motion interference noise
CN109222949A (en) * 2018-10-12 2019-01-18 杭州士兰微电子股份有限公司 Heart rate detection method and heartbeat detection device
CN109528163A (en) * 2018-11-15 2019-03-29 深圳市国通世纪科技开发有限公司 A kind of method and apparatus of sleep monitor
CN109528163B (en) * 2018-11-15 2021-11-05 深圳市苏仁智能科技有限公司 Sleep monitoring method and equipment
CN113271845A (en) * 2019-01-10 2021-08-17 株式会社电装 Biological information detection device
CN111643052A (en) * 2019-04-03 2020-09-11 上海铼锶信息技术有限公司 Method and system for reducing motion artifacts in pulse wave signals
CN110916639A (en) * 2019-12-23 2020-03-27 深圳市圆周率智能信息科技有限公司 Method, system, wearable device and computer-readable storage medium for acquiring exercise heart rate recovery rate
CN111714110A (en) * 2020-05-19 2020-09-29 成都云卫康医疗科技有限公司 Real-time heart rate calculation method based on PPG waveform
CN112494001A (en) * 2020-11-25 2021-03-16 青岛歌尔智能传感器有限公司 PPG signal quality evaluation method and device, wearable device and storage medium

Similar Documents

Publication Publication Date Title
CN105286845A (en) Movement noise elimination method suitable for wearable heart rate measurement device
CN110664390B (en) Heart rate monitoring system and method based on wrist strap type PPG and deep learning
Zhang et al. Combining ensemble empirical mode decomposition with spectrum subtraction technique for heart rate monitoring using wrist-type photoplethysmography
JP5369726B2 (en) Pulsation detection device and pulsation detection method
CN105105737B (en) Motion state rhythm of the heart method based on photoplethaysmography and spectrum analysis
CN105832289B (en) Apparatus and method for estimating biophysiological rate using Hilbert transform
CN106691425A (en) Wrist heart rate monitoring method of sports bracelet
CN104161505A (en) Motion noise interference eliminating method suitable for wearable heart rate monitoring device
WO2013179018A1 (en) Respiration rate extraction from cardiac signals
CN106264505A (en) A kind of heart rate spectral peak system of selection based on support vector machine
CN105943012B (en) Heart rate measuring method for removing motion noise in photoplethysmography signals
WO2011048729A1 (en) Pulse wave detection device and pulse wave detection method
CN106037743A (en) Method and device for extracting respiratory rate
CN205144547U (en) Motion noise eliminates system suitable for wearable heart rate measurement equipment
CN108042107A (en) A kind of PPG signals puppet difference correcting method
CN106236041B (en) A kind of algorithm and system measuring heart rate and respiratory rate in real time and accurately
CN105816165A (en) Real-time dynamic heart rate monitoring device and monitoring method
CN205144548U (en) Motion noise detecting system suitable for heart rate signal
TW201811261A (en) Signal detection method for accurately detecting the peak-to-peak interval of the photoplethysmography signal
US11452483B2 (en) Apparatus for eliminating motion artifacts by using PPG signal and method thereof
CN110801214A (en) Heart rate real-time detection method and system
JP2013106837A (en) Heart rate detection method, heart rate detector, and mental stress measuring apparatus
CN110151156B (en) Fetal heart rate extraction method and system based on in-window automatic weighted average algorithm
Shi et al. A Novel Method of sEMG Signal Segmentation
TWI655930B (en) Measuring device, measuring method and recording medium

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20160203

RJ01 Rejection of invention patent application after publication