CN105559766A - 一种基于ppg的手腕式动态心率实时测量方法 - Google Patents
一种基于ppg的手腕式动态心率实时测量方法 Download PDFInfo
- Publication number
- CN105559766A CN105559766A CN201510972536.4A CN201510972536A CN105559766A CN 105559766 A CN105559766 A CN 105559766A CN 201510972536 A CN201510972536 A CN 201510972536A CN 105559766 A CN105559766 A CN 105559766A
- Authority
- CN
- China
- Prior art keywords
- heart rate
- signal
- ppg
- real
- module
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 210000000707 wrist Anatomy 0.000 title claims abstract description 15
- 230000003595 spectral effect Effects 0.000 claims abstract description 24
- 238000001228 spectrum Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims description 14
- 238000010183 spectrum analysis Methods 0.000 claims description 6
- 230000002452 interceptive effect Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000009532 heart rate measurement Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 230000010247 heart contraction Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004118 muscle contraction Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
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
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- 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)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及心率测量领域,具体涉及一种基于PPG的手腕式动态心率实时测量方法。该方法由发射器控制模块、接收器模块、带通滤波模块,信号分析模块和频谱选择模块组成,心率信号依次通过上述模块的处理后得到最终精确的心率谱线。本方法无需额外噪声参考源、运算量小、运算速度快,可以在人体保持运动,甚至是强烈运动的情形下更准确、更快速的提取出心率信号。
Description
技术领域
本发明涉及心率测量领域,具体涉及一种基于PPG的手腕式动态心率实时测量方法。
背景技术
PPG(光电容积描记)技术进行人体运动心率的检测是光学无损检测技术在生物医学中的一个应用.它利用光电传感器,检测经过人体血液和组织吸收后的反射光强度的不同,描记出血管容积在心动周期内的变化,从得到的脉搏波形中计算出心率。
PPG信号是一种容易遭受运动干扰的信号,运动越大,干扰就越强烈。由于干扰噪声大多数发生在和有用信号相同的频率范围内,普通的滤波方法很难将心率信号提取出来。
目前已经有了大量使用PPG测量人体运动时心率方法的研究,但总体来说还没有一种可行的方法可以有效的,在运算能力低下和运算速度有限的单片机系统中实时运行。
CN201410397295提出使用了6轴运动传感器(加速度和陀螺仪)的方法提供噪声参考源,来提取有用的心率信号,这种方法主要有两个方面的不足。第一,噪声参考源的轴数量越多所需要的计算量就越大;其二,自适应滤波要求被分析的信号噪声与参考源有较高的相关性,而PPG信号的干扰源除了运动干扰之外,还有其他比如肌肉收缩、电源纹波、佩戴松紧程度等干扰,要有效的去除PPG信号的干扰信号,不能仅仅只考虑运动干扰,还应考虑其他干扰。
发明内容
为了解决上述技术问题,本发明的发明人提供了一种基于PPG的手腕式动态心率实时测量方法。
本方法无需额外噪声参考源、运算量小、运算速度快,可以在人体保持运动,甚至是强烈运动的情形下更准确、更快速的提取出心率信号。
为了实现上述目的,本发明提供了一种基于PPG的手腕式动态心率实时测量方法,具体实现方案如下。
本方法由发射器控制模块、接收器模块、带通滤波模块,信号分析模块和频谱选择模块组成,心率信号依次通过上述模块。
所述发射器控制模块用于调节驱动电流,按分时复用的方式打开驱动1路和驱动2路,驱动1路的驱动电流I1为20毫安,驱动2路的驱动电流I2为8毫安。
所述接收器模块用于接收信号并进行AD转换得到驱动1路和驱动2路的AD值,其中驱动1路的采样周期为2毫秒,驱动2路的采样周期为2微秒。
所述带通滤波模块用于屏蔽干扰信号,其通带频率为0.5Hz-4Hz。
所述信号分析模块用于频谱分析和相关性分析,频谱分析时将信号频率做傅里叶变换,使其分辨率提高4倍。
所述频谱选择模块用于选择心率谱线。
通过上述方案,本方法可以在人体保持运动,甚至是强烈运动的情形下更准确、更快速的提取出心率信号,并且无需额外噪声参考源、运算量小、运算速度快。
附图说明
图1是本发明实施例的心率实时测量方法的框架及信号流序图。
图2是本发明实施例的心率实时测量方法的分时复用演示图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
本发明实施例提出一种基于PPG的手腕式动态心率实时测量方法。
如图1所示,本发明的方法由发射器控制模块、接收器模块、带通滤波模块,信号分析模块和频谱选择模块组成,心率信号依次通过上述模块。
其中,所述发射器控制模块用于调节驱动1路和驱动2路的电流,具体的,是按照分时复用的方式打开驱动1路和驱动2路。驱动1路的驱动电流I1为20毫安,驱动2路的驱动电流I2为8毫安。驱动1路的电流设置目标是最大化具体硬件的信噪比,驱动2路的电流设置目标是降低信噪比。
所述接收器模块用于接收信号并进行AD(模数)转换得到驱动1路和驱动2路的AD值。如图2所示,所述发射器控制模块按照分时复用的方式分时打开驱动电路,如t1时间内,驱动1路电路打开,所述接收器接收信号并进行AD(模数)转换得到驱动1路AD值(模数转换结果);t3时间内驱动2路电路打开,所述接收器接收信号并进行AD转换得到驱动2路AD值。
其中采样周期为Ts,Ts满足Ts=t1+t2+t3+t4。t1=t3=2毫秒,t2=t4=2微妙。
本实施例要求t3时间内采集到的信号需要与t1时间内采集到的信号具备较高的相关性,最好的二者的相关系数大于0.9。
得到两路采样信号后,送入带通滤波器进行滤波,所述带通滤波模块用于屏蔽干扰信号,带通滤波器的通带可以根据人体正常心跳范围在每分钟30-240次,其对应的通带频率为0.5Hz-4Hz来设定。
经过带通滤波后,将两路信号送入信号分析模块,所述信号分析模块用于频谱分析和相关性分析。
频谱分析时为了提高频率分辨率,本实施例采用了ZOOMFFT(游标傅里叶变换),即对信号在0.5Hz-4Hz内进行放大,放大倍数D定义为4。
进行傅里叶变换点数选择为N=256,采样频率Fs=40Hz,因此原始分辨率为Fs/N=40Hz/256=0.15625Hz,对应的每分钟心跳次数的分辨率0.15625*60=9.375BPM(BPM即每分钟心跳次数)。
如果不对频率加以细化,无法直接作为心率结果输出的,采用ZOOMFFT之后,其频率分辨率提高了4倍,因此频率分辨率Fs/N/D=0.15625Hz/4=0.0390625Hz,对应的心跳分辨率提高为0.0390625*60=2.34375BPM。
所述频谱选择模块:用于选择心率谱线。
经过傅里叶变换得到两路信号的频谱后,需要配合两路信号的相关程度进行分析。
如果两路信号频谱的相似度大于0.8,则说明此时手腕使用者正在休息或只有少量轻微运动时的测量,或者手腕使用者运动时的干扰信号的频率和心跳频率谱线重合,此时选择驱动1路频率轴上对应频谱的能量值最高的谱线即为心率谱线。
如果两路信号的频谱相似度小于0.8,那么可以用前一次的心率谱线选择结果作为一个参考中心,以此为参考,选中能量谱线最高的2根谱线,选择的标准是:第二高谱线要在第一高谱线能量的60%以上,并且其他未选中的谱线能量要在最高谱线能量的30%以下。
根据范围内最高谱线的选择情况做如下分析:
如果范围内,只有一根最高能量谱线,第二高的能量谱线选择不成功,那么此谱线就是心率谱线。
如果范围内,两根最高的能量谱线选择都未成功,那么以前一次结果作为本次心率谱线。
如果两根谱线S1、S2都选择成功,要做的就是决策这两根谱线S1、S2中,哪一根是心率谱线,哪一根是干扰谱线。这时需要运算运动1路上S1、S2两根谱线的比值R1,再将两根谱线的索引在运动2路频率轴上找到对应的两根谱线,以相同方式求取比值R2,如果R1>R2,则以S1作为本次心率谱线,否则以S2作为心率谱线。
本发明无需额外噪声参考源、运算量小、运算速度快,可以在人体保持运动,甚至是强烈运动的情形下更准确、更快速的提取出心率信号。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。
Claims (6)
1.一种基于PPG的手腕式动态心率实时测量方法,其特征在于,本方法由发射器控制模块、接收器模块、带通滤波模块,信号分析模块和频谱选择模块组成,心率信号依次通过上述模块。
2.根据权利要求1所述的基于PPG的手腕式动态心率实时测量方法,其特征在于,所述发射器控制模块用于调节驱动电流,按分时复用的方式打开驱动1路和驱动2路,驱动1路的驱动电流I1为20毫安,驱动2路的驱动电流I2为8毫安。
3.根据权利要求1所述的基于PPG的手腕式动态心率实时测量方法,其特征在于,所述接收器模块用于接收信号并进行AD转换得到驱动1路和驱动2路的AD值,其中驱动1路的采样周期为2毫秒,驱动2路的采样周期为2微秒。
4.根据权利要求1所述的基于PPG的手腕式动态心率实时测量方法,其特征在于,所述带通滤波模块用于屏蔽干扰信号,其通带频率为0.5Hz-4Hz。
5.根据权利要求1所述的基于PPG的手腕式动态心率实时测量方法,其特征在于,所述信号分析模块用于频谱分析和相关性分析,频谱分析时将信号频率做傅里叶变换,使其分辨率提高4倍。
6.根据权利要求1所述的基于PPG的手腕式动态心率实时测量方法,其特征在于,所述频谱选择模块用于选择心率谱线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510972536.4A CN105559766A (zh) | 2015-12-23 | 2015-12-23 | 一种基于ppg的手腕式动态心率实时测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510972536.4A CN105559766A (zh) | 2015-12-23 | 2015-12-23 | 一种基于ppg的手腕式动态心率实时测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105559766A true CN105559766A (zh) | 2016-05-11 |
Family
ID=55870819
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510972536.4A Pending CN105559766A (zh) | 2015-12-23 | 2015-12-23 | 一种基于ppg的手腕式动态心率实时测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105559766A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107468227A (zh) * | 2017-09-12 | 2017-12-15 | 王沛 | 一种信号发射器和ppg信息采集设备 |
CN108294737A (zh) * | 2018-01-26 | 2018-07-20 | 深圳还是威健康科技有限公司 | 心率测量方法、装置及智能穿戴设备 |
CN108814935A (zh) * | 2018-04-20 | 2018-11-16 | 史琳 | 一种音律疾病诊疗系统 |
CN109044316A (zh) * | 2018-07-11 | 2018-12-21 | 歌尔科技有限公司 | 纯净动态心率信号计算方法、装置及智能穿戴设备 |
CN110108929A (zh) * | 2019-05-22 | 2019-08-09 | 电子科技大学 | 一种抗干扰型雷电流采集装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007037100A1 (ja) * | 2005-09-27 | 2007-04-05 | Citizen Holdings Co., Ltd. | 心拍計および心拍検出方法 |
CN101730503A (zh) * | 2007-03-15 | 2010-06-09 | 帝国创新有限公司 | 心率测量 |
CN105105737A (zh) * | 2015-08-03 | 2015-12-02 | 南京盟联信息科技有限公司 | 基于光电容积描记和谱分析的运动状态心率监测方法 |
US20150351646A1 (en) * | 2014-06-09 | 2015-12-10 | Stmicroelectronics S.R.L. | Method for the estimation of the heart-rate and corresponding system |
-
2015
- 2015-12-23 CN CN201510972536.4A patent/CN105559766A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007037100A1 (ja) * | 2005-09-27 | 2007-04-05 | Citizen Holdings Co., Ltd. | 心拍計および心拍検出方法 |
CN101272731A (zh) * | 2005-09-27 | 2008-09-24 | 西铁城控股株式会社 | 心率计以及心率检测方法 |
CN101730503A (zh) * | 2007-03-15 | 2010-06-09 | 帝国创新有限公司 | 心率测量 |
US20150351646A1 (en) * | 2014-06-09 | 2015-12-10 | Stmicroelectronics S.R.L. | Method for the estimation of the heart-rate and corresponding system |
CN105105737A (zh) * | 2015-08-03 | 2015-12-02 | 南京盟联信息科技有限公司 | 基于光电容积描记和谱分析的运动状态心率监测方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107468227A (zh) * | 2017-09-12 | 2017-12-15 | 王沛 | 一种信号发射器和ppg信息采集设备 |
CN107468227B (zh) * | 2017-09-12 | 2020-10-16 | 南京邮电大学 | 一种信号发射器和ppg信息采集设备 |
CN108294737A (zh) * | 2018-01-26 | 2018-07-20 | 深圳还是威健康科技有限公司 | 心率测量方法、装置及智能穿戴设备 |
CN108294737B (zh) * | 2018-01-26 | 2020-11-13 | 深圳市元征科技股份有限公司 | 心率测量方法、装置及智能穿戴设备 |
CN108814935A (zh) * | 2018-04-20 | 2018-11-16 | 史琳 | 一种音律疾病诊疗系统 |
CN109044316A (zh) * | 2018-07-11 | 2018-12-21 | 歌尔科技有限公司 | 纯净动态心率信号计算方法、装置及智能穿戴设备 |
CN109044316B (zh) * | 2018-07-11 | 2021-02-02 | 歌尔科技有限公司 | 纯净动态心率信号计算方法、装置及智能穿戴设备 |
CN110108929A (zh) * | 2019-05-22 | 2019-08-09 | 电子科技大学 | 一种抗干扰型雷电流采集装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105559766A (zh) | 一种基于ppg的手腕式动态心率实时测量方法 | |
KR101406130B1 (ko) | 박동 주기 산출장치 및 이것을 포함한 생체 센서 | |
CN102144916B (zh) | 一种可自动调压的多通道脉搏信号检测方法及装置 | |
CN205083467U (zh) | 估计心率的设备和系统 | |
FI114282B (fi) | Menetelmä, järjestely ja sykemittari sydämen lyönnin tunnistamiseksi | |
CN103099615B (zh) | 一种消除运动心电信号干扰的方法和装置 | |
CN101991410B (zh) | 一种脉率搜索和计算方法 | |
KR20100008239A (ko) | 피피지 신호의 동잡음 제거방법 | |
CN104873186A (zh) | 一种可穿戴的动脉检测装置及其数据处理方法 | |
CN102421354A (zh) | 具有脉冲和通道切换adc噪声滤波器和衍生导联纠错器的ecg设备 | |
US20120283535A1 (en) | Method and system for pulse measurement | |
CN108937957A (zh) | 检测方法、装置及检测设备 | |
JP2018505737A (ja) | 心拍数及び心拍変動を検出する装置 | |
CN115702782A (zh) | 基于深度学习的心率检测方法及可穿戴设备 | |
CN104068841B (zh) | 一种测量心脏收缩时间参数的测量方法及装置 | |
CN107468232A (zh) | 胎心监护装置和方法 | |
TWI505816B (zh) | 血氧飽和度檢測方法及裝置 | |
CN109330582A (zh) | 基于心电信号分析的心率及其特征指标检测方法 | |
CN102670182A (zh) | 脉搏波采集质量实时分析装置 | |
CN109620195B (zh) | 基于可穿戴设备的多波长信号融合心率检测方法及系统 | |
US20160235371A1 (en) | Pulsimeter, frequency analysis device, and pulse measurement method | |
CN106214143B (zh) | 瞬时运动干扰识别方法及心率置信度计算方法 | |
CN111329463A (zh) | 一种基于ppg心率测量的运动伪差消除系统及其实现方法 | |
CN109620198A (zh) | 心血管指数检测、模型训练方法及装置 | |
CN110730630B (zh) | 心率检测方法、装置、芯片、电子装置及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160511 |
|
WD01 | Invention patent application deemed withdrawn after publication |