CN102988036B - Method for measuring pulse rate - Google Patents

Method for measuring pulse rate Download PDF

Info

Publication number
CN102988036B
CN102988036B CN201210575419.0A CN201210575419A CN102988036B CN 102988036 B CN102988036 B CN 102988036B CN 201210575419 A CN201210575419 A CN 201210575419A CN 102988036 B CN102988036 B CN 102988036B
Authority
CN
China
Prior art keywords
pulse
frequency
time domain
pulse frequency
pulse signal
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
CN201210575419.0A
Other languages
Chinese (zh)
Other versions
CN102988036A (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.)
Zhongluan Polytron Technologies Inc
Original Assignee
Institute of Automation of Chinese Academy of Science
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 Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201210575419.0A priority Critical patent/CN102988036B/en
Publication of CN102988036A publication Critical patent/CN102988036A/en
Application granted granted Critical
Publication of CN102988036B publication Critical patent/CN102988036B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention provides a method for measuring pulse rate. The method comprises the steps of using a pulse sensor to acquire a time domain pulse signal; obtaining first pulse rate n1 according to the time domain pulse signal acquired by the pulse sensor; performing Fourier transform to the time domain pulse signal to obtain a frequent domain pulse signal; confirming a wave crest corresponding value of the frequency domain pulse signal in a frequency range from 0.65Hz-3.20Hz, namely a pulse frequency value, and obtaining a second pulse rate n2 through the pulse rate; and calculating final pulse rate of a tested object in unit time through the first pulse rate n1 and the second pulse rate n2. The method for measuring the pulse rate integrates a time domain analyzing method and a frequency domain analyzing method, uses acceleration information to perform feedback adjustment and control to a pulse value, and accordingly obtains accurate and reliable pulse information.

Description

Measure the method for pulse frequency
Technical field
The present invention relates to medical device field, relate in particular to a kind of method of measuring pulse frequency.
Background technology
Pulse signal is a kind of faint physiological signal of non-linear, non-stationary, and the integrated information of the aspects such as form, intensity, speed and the rhythm and pace of moving things that the ripple of fighting presents can reflect the flow characteristic of many physiological and pathologicals in cardiovascular system of human body.Often it is detected to analysis, can reach morning prevention to cardiovascular disease and the object of early treatment.At present, no matter be Traditional Chinese Medicine or modern doctor trained in Western medicine, in to the inspection of cardiovascular disease, all using its wave character and change the important evidence as evaluation cardiovascular system of human body physiology, pathological state.
Pulse frequency, mostly by static measurement, is carried out feature extraction in time domain to pulse signal at present.Temporal analysis carries out feature extraction to pulse wave time-domain signal, calculates Pulse Rate per minute, i.e. pulse frequency, and the advantage of his maximum is exactly that speed is fast, and real-time is high, easily meets the demands.But Time-Domain algorithm is very high to the prescription of signal, lower in pulse signal quality, especially, under kinestate, adopt method of the prior art can not meet Measurement accuracy pulse frequency.In addition, Static Detection correspondence has significant limitation, if continuous record human body pulse frequency dynamically, will inevitably disclose some rule of human vas pathological changes, is more convenient for monitoring healthy in people's daily life simultaneously.
Summary of the invention
(1) technical problem that will solve
For solving above-mentioned one or more problems, the invention provides a kind of method of measuring pulse frequency, to obtain exactly pulse frequency information.
(2) technical scheme
According to an aspect of the present invention, provide a kind of method of measuring pulse frequency.The method of this measurement pulse frequency comprises: gather time domain pulse signal by pulse transducer; Obtain the first pulse frequency n according to the time domain pulse signal of pulse transducer collection 1; Time domain pulse signal is carried out to Fourier transformation, obtain frequency domain pulse signal; For frequency domain pulse signal, determine its crest respective value in 0.65Hz-3.20Hz frequency range, i.e. pulse frequency value, obtains the second pulse frequency n by this pulse frequency 2; And by the first pulse frequency n 1with the second pulse frequency n 2calculate the final pulse frequency of the unit interval of tested object: N=n 1× w 1+ n 2× w 2, wherein, w 1+ w 2=1, w 1and w 2be respectively the weight of the first pulse frequency and the second pulse frequency, the acceleration statistical value class of living in of tested object is higher, the weight w of the second pulse frequency 2larger.
(3) beneficial effect
Can find out from technique scheme, the method that the present invention measures pulse frequency has following beneficial effect:
(1) temporal analysis is combined with frequency domain analysis, and pulse value is carried out to feedback regulation with acceleration information, thus the pulse information obtaining more accurately, reliable;
(2) can be under motion conditions and acquired signal situation of low quality under measurement pulse frequency, stable performance, result is accurate.
Brief description of the drawings
Fig. 1 is the flow chart of measuring pulse frequency method according to the embodiment of the present invention;
Fig. 2 measures in pulse frequency method described in Fig. 1, is identified the schematic diagram of pulse wave by the main crest of pulse wave (A), dicrotic notch crest (B), dicrotic wave crest (C) and dicrotic wave trough (D).
Detailed description of the invention
For making the object, technical solutions and advantages of the present invention clearer, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail.
It should be noted that, in accompanying drawing or description description, similar or identical part is all used identical figure number.The implementation that does not illustrate in accompanying drawing or describe is form known to a person of ordinary skill in the art in affiliated technical field.In addition, although the demonstration of the parameter that comprises particular value can be provided herein, should be appreciated that, parameter is without definitely equaling corresponding value, but can in acceptable error margin or design constraint, be similar to corresponding value.
In one exemplary embodiment of the present invention, provide a kind of method of measuring pulse frequency.As shown in Figure 1, this measurement pulse frequency method comprises:
Steps A, by the pulse transducer time domain pulse signal of collection 2s time;
This pulse transducer is installed with in the position of tested object pulse, as chest, wrist, finger etc., the pulse of the tested object of Real-time Collection.Preferably, this pulse transducer is attached at tested object chest, to facilitate tested object to carry out pulse frequency monitoring under kinestate.
In this step, the sampling duration of time domain pulse signal is 1-2s, and 2s is more quite a lot of, and because for frequency domain algorithm, sample frequency is fixed, and sampling number is more, and frequency resolution is higher.Comprehensive requirement of real-time again, 2s is more quite a lot of.Certainly, in the situation that requirement of real-time is not high, also can proper extension data acquisition time, for example: 4s, 8s, 10s, 20s, 60s etc.
Step B, the acceleration signal of synchronous acquisition 2s, measures the acceleration of tested object in X, Y and tri-directions of Z: a by 3-axis acceleration sensor x, a yand a z,, the counting as n of collection; Then add up the value of three axial acceleration in [3,3], [6,6], [10,10], [15,15], counting between [20,20], is designated as respectively a x1, a x2, a x3, a x4, a x5; a y1, a y2, a y3, a y4, a y5; a z1, a z2, a z3, a z4, a z5.Wherein a x1, a x2, a x3, a x4, a x5for tested object at the acceleration of directions X respectively at [a1, a1], [a2, a2], [a3, a3], [a4, a4], counting between [a5, a5]; a y1, a y2, a y3, a y4, a y5for tested object at the acceleration of Y-direction respectively at [a1, a1], [a2, a2], [a3, a3], [a4, a4], counting between [a5, a5]; a z1, a z2, a z3, a z4, a z5for tested object at the acceleration of Z direction respectively at [a1, a1], [a2, a2], [a3, a3], [a4, a4], counting between [a5, a5].
In the present embodiment, the 3-axis acceleration sensor of employing, can obtain the acceleration of tested object in three directions.It will be apparent to those skilled in the art that and also can adopt the sensor that can obtain one direction acceleration, and needn't strictly distinguish three directions, thereby reduce costs.
In this enforcement, pulse transducer and 3-axis acceleration sensor are integrated on an equipment.
Step C, judge tested object at the acceleration of X, Y and tri-directions of Z the statistical value a at the 5th class x5, a y5, a z5whether be all greater than predetermined threshold value T min, if so, execution step D, otherwise, re-execute steps A;
Under normal circumstances, this predetermined threshold value T minget 0.8n, n is the number that once gathers acceleration information on each axle.In this step, a x5, a y5, a z5whether be all greater than predetermined threshold value T, limit exactly acceleration and exceed [20,20] interval ratio data and be less than 0.2n, be also just equivalent to accekeration and be greater than 0.8n in [20,20] interval ratio.
Experiment showed, to exceed [20,20] interval ratio by acceleration more large-signal is poorer, in the time exceeding 0.2n, signal just can not have been differentiated, and the follow-up pulse value drawing will be very inaccurate.Therefore, by very necessary within the specific limits the acceleration restriction of tested object, otherwise will need resampling.
Step D, adopts the band filter that free transmission range is 0.65Hz-3.20Hz to carry out filtering to time domain pulse signal;
In the prior art, conventionally adopt the method for wavelet filtering to carry out filtering to time domain pulse signal, but Wavelets Filtering Algorithm is comparatively complicated, holding time is long.And method processing procedure of the present invention is simple, the time is short, is conducive to Real-time Collection pulse data.
Step e, to time domain pulse signal, identifies each pulse wave by crest and the trough of the main ripple of pulse wave, dicrotic notch and dicrotic wave, obtains the first pulse frequency n according to the time domain pulse signal gathering 1;
As shown in Figure 2, adopt and identify pulse wave by the crest (C) of the main crest of pulse wave (A), dicrotic notch crest (B), dicrotic wave and the trough (D) of dicrotic wave, its advantage is to realize accurate identification.It will be apparent to those skilled in the art that and also can take the mode of the main ripple of Direct Recognition pulse wave to identify pulse wave, doing like this accuracy rate can be slightly poor, but can't affect enforcement of the present invention.
Step F, carries out Fourier (FFT) conversion by time domain pulse signal, obtains frequency domain pulse signal;
Step G, for frequency domain pulse signal, determines its crest respective value in 0.65Hz-3.20Hz frequency range, i.e. pulse frequency value, obtains the second pulse frequency n by this pulse frequency 2;
Experiment showed, that the maximum feature of utilizing frequency domain method to carry out pulses measure is exactly good stability, capacity of resisting disturbance is strong, and still, if simple proportion method can make again real-time variation, therefore the present embodiment adopts the pulse algorithm that time-frequency mixes.
Step H, utilizes following formula, by the first pulse frequency n 1with the second pulse frequency n 2calculate the final pulse frequency of the unit interval of tested object:
N=n 1×w 1+n 2×w 2
Wherein, w 1and w 2be respectively the weight of the first pulse frequency and the second pulse frequency, both are relevant at each interval statistical value at the acceleration of X, Y and tri-directions of Z to tested object, and both sums are 1.
In the present embodiment, the present invention's degree of will speed up is divided into five class statistics, is respectively [3,3], and [6,6], [10,10], [15,15], [20,20], then according to three axial acceleration statistical value a on each class x1, a x2, a x3, a x4, a x5; a y1, a y2, a y3, a y4, a y5; a z1, a z2, a z3, a z4, a z5, and the sampling n that always counts determines the weights of corresponding the first pulse frequency and the second pulse frequency:
As min (a x1, a y1, a z1when)>=0.8n, w 1=0.9, w 2=0.1;
When being false above but min (a x2, a y2, a z2when)>=0.8n, w 1=0.7, w 2=0.3;
When being false above but min (a x3, a y3, a z3when)>=0.8n, w 1=0.5, w 2=0.5;
When being false above but min (a x4, a y4, a z4when)>=0.8n, w 1=0.3, w 2=0.7;
When being false above but min (a x5, a y5, a z5when)>=0.8n, w 1=0.1, w 2=0.9.
So far, the present embodiment measurement pulse frequency process finishes.
The invention provides one applicable in people's daily life, especially, under motion conditions, measure the method for pulse frequency.The data that the method utilizes pulse transducer collection to come up, then in conjunction with the acceleration information of a 3-axis acceleration sensor, time domain pulse information and frequency domain pulse information are combined to calculating Pulse Rate.The present invention has substantial degradation to pulse prescription, and stable performance is calculated simply, and the result of calculation of pulse is also more accurate.
It should be noted that, the above-mentioned definition to each element is not limited in various concrete structures or the shape in embodiment, mentioned, and those of ordinary skill in the art can know simply and replace it, for example:
The weighted value of (1) first pulse frequency and the second pulse frequency can be adjusted as required;
(2) number of acceleration class, upper lower limit value wherein can be adjusted as required;
(3) 3-axis acceleration sensor can replace with other multiaxis acceleration transducers.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the foregoing is only specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any amendment of making, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (7)

1. a method of measuring pulse frequency, is characterized in that, comprising:
Gather time domain pulse signal by pulse transducer;
Obtain the first pulse frequency n according to the described time domain pulse signal of pulse transducer collection 1;
Described time domain pulse signal is carried out to Fourier transformation, obtain frequency domain pulse signal;
For described frequency domain pulse signal, determine its crest respective value in 0.65Hz-3.20Hz frequency range, i.e. pulse frequency value, obtains the second pulse frequency n by this pulse frequency 2; And
By described the first pulse frequency n 1with described the second pulse frequency n 2calculate the final pulse frequency of the unit interval of tested object: N=n 1× w 1+ n 2× w 2, wherein, w 1+ w 2=1, w 1and w 2the weight that is respectively the first pulse frequency and the second pulse frequency, it is higher that the acceleration of tested object is added up class of living in, the weight w of the second pulse frequency 2larger;
Wherein, gather the acceleration of tested object in X, Y and tri-directions of Z: a by 3-axis acceleration sensor x, a yand a z, the counting as n of collection; Add up tested object at the acceleration of X, Y and tri-directions of Z at 5 class [a1, a1], [a2, a2], [a3, a3], [a4, a4], counting between [a5, a5], is designated as respectively a x1, a x2, a x3, a x4, a x5; a y1, a y2, a y3, a y4, a y5; a z1, a z2, a z3, a z4, a z5, wherein a1<a2<a3<a4LEssT. LTssT.LTa5, described by the first pulse frequency n 1with the second pulse frequency n 2calculate in the step of final pulse frequency of the unit interval of tested object:
As min (a x1, a y1, a z1) when≤0.8n, w 1=0.9, w 2=0.1;
When being false above but min (a x2, a y2, a z2) when≤0.8n, w 1=0.7, w 2=0.3;
When being false above but min (a x3, a y3, a z3) when≤0.8n, w 1=0.5, w 2=0.5;
When being false above but min (a x4, a y4, a z4) when≤0.8n, w 1=0.3, w 2=0.7; And
When being false above but min (a x5, a y5, a z5) when≤0.8n, w 1=0.1, w 2=0.9.
2. method according to claim 1, is characterized in that, the described described time domain pulse signal according to pulse transducer collection obtains the first pulse frequency n 1step before also comprise:
Judge tested object at the acceleration of X, Y and tri-directions of Z the statistical value a at [a5, a5] x5, a y5, a z5whether be all greater than predetermined threshold value T min, if not, again gather time domain pulse signal by pulse transducer.
3. method according to claim 2, is characterized in that, described T min=0.8n.
4. method according to claim 3, is characterized in that, described a1=3; A2=6; A3=10; A4=15; A5=20.
5. according to the method described in any one in claim 1 to 4, it is characterized in that, the described described time domain pulse signal according to pulse transducer collection obtains the first pulse frequency n 1step comprise:
To time domain pulse signal, identify each pulse wave by crest and the trough of the main ripple of pulse wave, dicrotic notch and dicrotic wave, obtain the first pulse frequency n wave period by pulse 1.
6. according to the method described in any one in claim 1 to 4, it is characterized in that, the described time domain pulse signal according to pulse transducer collection obtains the first pulse frequency n 1step before also comprise:
Adopt band filter to carry out filtering to described time domain pulse signal, the free transmission range of this band filter is 0.65Hz-3.20Hz.
7. according to the method described in any one in claim 1 to 4, it is characterized in that, in the described step that is gathered time domain pulse signal by pulse transducer, the duration of time domain pulse signal is between between 1s to 2s.
CN201210575419.0A 2012-12-26 2012-12-26 Method for measuring pulse rate Active CN102988036B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210575419.0A CN102988036B (en) 2012-12-26 2012-12-26 Method for measuring pulse rate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210575419.0A CN102988036B (en) 2012-12-26 2012-12-26 Method for measuring pulse rate

Publications (2)

Publication Number Publication Date
CN102988036A CN102988036A (en) 2013-03-27
CN102988036B true CN102988036B (en) 2014-08-06

Family

ID=47917495

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210575419.0A Active CN102988036B (en) 2012-12-26 2012-12-26 Method for measuring pulse rate

Country Status (1)

Country Link
CN (1) CN102988036B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109414203A (en) * 2016-06-30 2019-03-01 美国亚德诺半导体公司 Online heart rate estimation based on optical measurement

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10216893B2 (en) * 2010-09-30 2019-02-26 Fitbit, Inc. Multimode sensor devices
CN103230267B (en) * 2013-05-14 2015-06-03 北京理工大学 Anti-movement-interference extraction method for pulse rates
CN103549945B (en) * 2013-10-31 2015-07-15 广州视源电子科技股份有限公司 Method for recognizing pulse rate and blood oxygen saturation degree through cardiac contraction process characteristic
CN104739399A (en) * 2014-02-10 2015-07-01 北京金日吉通科贸有限公司 Method for detecting movement offset in pulse wave heart rate calculation
CN106606358A (en) * 2015-10-23 2017-05-03 富港电子(昆山)有限公司 Heart rate measuring method and apparatus
CN106419863B (en) * 2016-11-23 2019-05-17 深圳市维亿魄科技有限公司 It is drunk the monitoring method and wrist wearing devices of state based on pulse wave Energy distribution monitoring human body
CN108056769B (en) * 2017-11-14 2020-10-16 深圳市大耳马科技有限公司 Vital sign signal analysis processing method and device and vital sign monitoring equipment
CN109077711B (en) * 2018-08-20 2021-08-10 深圳市元征科技股份有限公司 Dynamic heart rate data acquisition method and device, wearable device and readable storage medium
CN109124610B (en) * 2018-09-26 2021-07-30 深圳星脉医疗仪器有限公司 Anti-interference method and device for non-invasive blood pressure measurement
WO2020133562A1 (en) * 2018-12-29 2020-07-02 深圳迈瑞生物医疗电子股份有限公司 Mobile monitoring device, and physiological signal adjustment and processing method
CN109528216A (en) * 2019-01-18 2019-03-29 京东方科技集团股份有限公司 The detection method and device of fetal hemoglobin saturation
CN110840413A (en) * 2019-10-14 2020-02-28 广州市倍尔康医疗器械有限公司 Temperature measurement method and system based on comprehensive physiological characteristics and storage medium
TWI777103B (en) * 2019-11-06 2022-09-11 達爾生技股份有限公司 Electronic device and blood oxygen correction method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6099478A (en) * 1997-03-18 2000-08-08 Seiko Epson Corporation Pulse counter and pulse display method
US6198951B1 (en) * 1997-09-05 2001-03-06 Seiko Epson Corporation Reflection photodetector and biological information measuring instrument
US20040086060A1 (en) * 2002-10-04 2004-05-06 Keisuke Tsubata Pulse wave detecting apparatus and fourier transform process apparatus
CN1531901A (en) * 2003-03-19 2004-09-29 ������������ʽ���� Information collector and pulse meter
CN1954771A (en) * 2003-03-19 2007-05-02 精工爱普生株式会社 Pulse meter, method for controlling pulse meter and wristwatch-type information device
CN101039617A (en) * 2004-10-15 2007-09-19 普尔塞特拉瑟技术有限公司 Motion cancellation of optical input signals for physiological pulse measurement
CN101991418A (en) * 2009-08-14 2011-03-30 深圳市理邦精密仪器股份有限公司 Method for improving respiratory rate detection accuracy

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5562805B2 (en) * 2010-11-05 2014-07-30 オータックス株式会社 Pulse rate measuring method and blood oxygen saturation measuring method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6099478A (en) * 1997-03-18 2000-08-08 Seiko Epson Corporation Pulse counter and pulse display method
US6198951B1 (en) * 1997-09-05 2001-03-06 Seiko Epson Corporation Reflection photodetector and biological information measuring instrument
US20040086060A1 (en) * 2002-10-04 2004-05-06 Keisuke Tsubata Pulse wave detecting apparatus and fourier transform process apparatus
CN1531901A (en) * 2003-03-19 2004-09-29 ������������ʽ���� Information collector and pulse meter
CN1954771A (en) * 2003-03-19 2007-05-02 精工爱普生株式会社 Pulse meter, method for controlling pulse meter and wristwatch-type information device
CN101039617A (en) * 2004-10-15 2007-09-19 普尔塞特拉瑟技术有限公司 Motion cancellation of optical input signals for physiological pulse measurement
CN101991418A (en) * 2009-08-14 2011-03-30 深圳市理邦精密仪器股份有限公司 Method for improving respiratory rate detection accuracy

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于卡尔曼滤波的动态脉搏波处理和脉率提取;蒋伟平等;《计算机与现代化》;20110228(第2期);全文 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109414203A (en) * 2016-06-30 2019-03-01 美国亚德诺半导体公司 Online heart rate estimation based on optical measurement

Also Published As

Publication number Publication date
CN102988036A (en) 2013-03-27

Similar Documents

Publication Publication Date Title
CN102988036B (en) Method for measuring pulse rate
CN109875529B (en) Vital sign detection method and system based on ultra-wideband radar
DK178263B1 (en) Method and system of detecting seizures
CN104382571A (en) Method and device for measuring blood pressure upon radial artery pulse wave conduction time
CN107970590B (en) A kind of running body-building data system and method based on Android platform
CN101732050A (en) Photoelectric volume wave-based breathing rate monitoring method
CN107788969A (en) The automatic testing method of QRS complex in a kind of electrocardiosignal
CN103315753A (en) Portable sentiment analyzer and sentiment analysis method
CN103932693A (en) Method for measuring human body heart rate on basis of mobile phone image
CN109009124A (en) Respiratory rate measurement and object localization method based on ULTRA-WIDEBAND RADAR
CN112263242A (en) Respiration detection and mode classification method based on FMCW radar
CN106344034A (en) Sleep quality evaluation system and method
CN106308801A (en) Method for detecting human breathing rate by utilizing smart phone
CN105243285A (en) Big data health forecast system
CN110115583A (en) The method and apparatus of monitoring of respiration
CN107495939A (en) Live biometric monitoring method, device and system
CN113017559A (en) Vital sign extraction algorithm and system based on piezoelectric film sensor
CN105138823A (en) Method for detecting quality of physiological signal based on self-correlation function
CN106419884A (en) Heart rate calculating method and system based on wavelet analysis
US20080281216A1 (en) Cardiac Arrhythmias Analysis of Electrophysiological Signals Based on Symbolic Dynamics
CN106913335A (en) A kind of detection method of apnea detection system
JPWO2004098409A1 (en) Extraction method and device for biological signals such as heartbeat and respiration
CN110801214A (en) Heart rate real-time detection method and system
CN104305958B (en) The photoelectricity volume ripple Multivariate analysis method of a kind of pole autonomic nerve state in short-term
CN106214143B (en) Transient motion interference identification method and heart rate confidence calculations method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: ZHONGLUAN TECHNOLOGY CO., LTD.

Free format text: FORMER OWNER: INSTITUTE OF AUTOMATION, CHINESE ACADEMY OF SCIENCES

Effective date: 20140728

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100190 HAIDIAN, BEIJING TO: 063021 TANGSHAN, HEBEI PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20140728

Address after: The more he Zhen Dong Liu Tun Cun East Hebei province Tangshan City Kaiping District 063021

Patentee after: Zhongluan Technology Co., Ltd.

Address before: 100190 Zhongguancun East Road, Beijing, No. 95, No.

Patentee before: Institute of Automation, Chinese Academy of Sciences

C56 Change in the name or address of the patentee
CP01 Change in the name or title of a patent holder

Address after: The more he Zhen Dong Liu Tun Cun East Hebei province Tangshan City Kaiping District 063021

Patentee after: Zhongluan Polytron Technologies Inc

Address before: The more he Zhen Dong Liu Tun Cun East Hebei province Tangshan City Kaiping District 063021

Patentee before: Zhongluan Technology Co., Ltd.