CN105816165B - 一种实时动态心率监测装置及监测方法 - Google Patents
一种实时动态心率监测装置及监测方法 Download PDFInfo
- Publication number
- CN105816165B CN105816165B CN201610288485.8A CN201610288485A CN105816165B CN 105816165 B CN105816165 B CN 105816165B CN 201610288485 A CN201610288485 A CN 201610288485A CN 105816165 B CN105816165 B CN 105816165B
- Authority
- CN
- China
- Prior art keywords
- heart rate
- frequency
- signal
- frequency point
- pulse wave
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012544 monitoring process Methods 0.000 title claims abstract description 27
- 238000012806 monitoring device Methods 0.000 title claims abstract description 26
- 238000001914 filtration Methods 0.000 claims abstract description 56
- 238000001228 spectrum Methods 0.000 claims abstract description 47
- 238000004458 analytical method Methods 0.000 claims abstract description 40
- 230000008030 elimination Effects 0.000 claims abstract description 31
- 238000003379 elimination reaction Methods 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 230000009466 transformation Effects 0.000 claims description 18
- 238000000605 extraction Methods 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 5
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 7
- 230000033001 locomotion Effects 0.000 abstract description 7
- 210000003205 muscle Anatomy 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000036541 health Effects 0.000 description 4
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009532 heart rate measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011426 transformation method Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000037081 physical activity Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000036391 respiratory frequency Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000000707 wrist 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/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
-
- 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
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- 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
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- 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
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Power Engineering (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明公开了一种实时动态心率监测装置及监测方法，装置包括基线漂移消除模块，带通滤波模块，频域分析模块，心率频点选择模块；所述基线漂移消除模块用于消除导致脉搏波信号的基准线出现漂移的干扰信号；所述带通滤波模块用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号；所述频域分析模块用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱；所述心率频点选择模块用于获取心率对应的频点并输出频点。本发明排除人体运动、肌肉活动对心率分析的影响，减少了硬件成本，可以灵活地选择关注频段区间，仅在关注频段区间内提高信号的频域精度，避免带来过大的计算和存储开销。
Description
技术领域
本发明涉及电子领域，尤其涉及一种实时动态心率监测装置及监测方法。
背景技术
随着移动互联网技术在医疗健康领域的应用与发展，在市场上出现了智能手表、智能手环、智能腕带等大量形态各异，具备测量心率、血压、血氧浓度、呼吸频率等医用生理参数的可穿戴式移动健康产品。心率定义为人体每分钟的心跳次数，是评估人健康状态的一个重要医用常规生理参数。心率监测对于疾病风险预警、病症诊断、年度例行体检具有十分重要的意义。特别地，健身活动、户外跑步等运动方式对于实时动态心率监测有着广泛的应用需求。
目前，绝大多数移动健康产品采用的心率监测技术原理是光电透射测量法。在产品的硬件设计上，与人体皮肤接触的传感器会发出一束光打在皮肤上，同时测量经皮肤反射或透射的光。因为血液对特定波长的光有吸收作用，心脏泵血的过程直接影响传感器测得的光信号强度改变，硬件按照设定的采样率记录信号强度变化采集原始数据，即脉搏波信号。数据分析软件单元运行心率监测算法处理脉搏波信号，输出心率值。心率监测算法是心率测量产品的关键核心技术，决定了心率测量值的准确性和可靠性。实际应用中，心率监测包括静态心率监测和实时动态心率监测，后者具有更广的应用空间，同时也对现有技术提出了很大的挑战。
通过实际测试发现，目前绝大多数移动健康产品的动态心率监测算法的技术现状是传统的信号时域波形分析或信号频域分析结合加速度计读数辅助判断，在运动状态下的人体心率实时监测这一应用场景中暴露了以下缺点。第一，由于噪声干扰导致信号时域波形上的特征点不太明显，会造成算法无法获取完整输入信息；第二，信号时域波形匹配的准则设置过多，算法参数具体数值设定存在困难；第三，对于嵌入式模块而言，波形匹配算法的计算复杂度较大；第四，传统的信号频域分析方法在提高频谱精度时会增大计算和数据存储开销；第五，加速度计增加了硬件成本，同时也增加计算、存储和能耗方面的资源开销。
因此，现有技术还有待于改进和发展。
发明内容
鉴于现有技术的不足，本发明目的在于提供一种实时动态心率监测装置及监测方法，旨在解决现有技术中动态心率监测装置在监测心率时，算法计算复杂度高，硬件成本高，计算、存储和能源方面的资源开销大的缺陷。
本发明的技术方案如下：
一种实时动态心率监测装置，其中，装置包括基线漂移消除模块，带通滤波模块，频域分析模块，心率频点选择模块；
所述基线漂移消除模块用于消除导致脉搏波信号的基准线出现漂移的干扰信号；
所述带通滤波模块用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号；
所述频域分析模块用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱；
所述心率频点选择模块用于获取心率对应的频点并输出频点；
所述基线漂移消除模块与所述带通滤波模块连接，所述带通滤波模块与所述频域分析模块连接，所述频域分析模块还与所述心率频点选择模块连接。
所述的实时动态心率监测装置，其中，所述基线漂移消除模块具体包括基线漂移趋势项信号提取单元和信号线性叠加单元；
所述基线漂移趋势项信号提取单元用于获取与原脉搏波信号等长的基线漂移趋势项信号；
所述信号线性叠加单元用于将原脉搏波信号减掉基线漂移趋势项信号，得到去掉基线漂移的脉搏波信号；
所述基线漂移趋势项信号提取单元与所述信号线性叠加单元连接。
所述的实时动态心率监测装置，其中，所述带通滤波模块具体包括滤波参数设置单元和滤波单元，
所述滤波参数设置单元用于设置通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数；
所述滤波单元用于获取滤波模块阶数和系数后，将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波，获取去除噪声后的脉搏波信号；
所述滤波参数设置单元和滤波单元连接。
所述的实时动态心率监测装置，其中，所述频域分析模块具体包括特定频段区间设置单元、频点功率计算单元和信号频谱获取单元；
所述特定频段区间设置单元用于设置特定频段区间的起始频点、结束频点、频点细分数目；
所述频点功率计算单元用于对于特定频段区间的每个频点，分别用第一变换多项式计算频域信号的实部，用第二变换多项式计算频域信号的虚部，根据频域信号的实部和虚部计算得到脉搏波信号在频点的功率；
所述信号频谱获取单元用于获取特定频段区间内的脉搏波信号在每个频点的功率，叠加后生成脉搏波信号在整个关注频段区间的信号频谱；
所述频点功率计算单元分别与所述特定频段区间设置单元、所述信号频谱获取单元连接。
所述的实时动态心率监测装置，其中，所述心率频点选择模块具体包括心率频点选择单元和心率频点动态跟踪单元，
所述心率频点选择单元用于在初始状态下获取脉搏波信号频谱中谱峰所在的频点位置作为心率频点；
所述心率频点动态跟踪单元用于在预先设置的固定跟踪周期内，以上一跟踪周期输出的心率频点为中心，在中心左右一个特定宽度的频谱范围内进行高精度的频域分析，选择谱峰所在的频点位置作为心率频点输出；
所述心率频点选择单元与所述心率频点动态跟踪单元连接。
所述的实时动态心率监测装置，其中，所述第一变换多项式和所述第二变换多项式为正交多项式。
所述的实时动态心率监测装置，其中，所述第一变换多项式为勒让德多项式，雅可比多项式，拉盖尔多项式，切比雪夫多项式，埃尔米特多项式中的一种；
所述第二变换多项式为正交多项式为勒让德多项式，雅可比多项式，拉盖尔多项式，切比雪夫多项式，埃尔米特多项式中的一种。
一种基于所述的实时动态心率监测装置的监测方法，其中，方法包括步骤：
A、获取脉搏波信号，通过基线漂移消除模块消除基线漂移干扰信号；
B、带通滤波模块对基线漂移消除模块输出的信号进行滤波，保留属于心率频段的信号分量；
C、频域分析模块获取脉搏波信号在预先设置的特定频段区间内的信号频谱；
D、心率频点选择模块根据信号频率中获取心率对应的频点并输出。
所述的实时动态心率监测方法，其中，所述步骤A具体包括：
A1、通过信号滤波法或曲线拟合法获取与原脉搏波信号等长的基线漂移趋势项信号；
A2、将原脉搏波信号，基线漂移趋势项信号使用行向量或者列向量存储，然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信，得到去掉基线漂移的脉搏波信号。
所述的实时动态心率监测方法，其中，所述步骤B具体包括：
B1、获取预先设置的通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数；
B2、获取滤波模块阶数和系数后，将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波，获取去除噪声后的脉搏波信号。
本发明提供了一种实时动态心率监测装置及监测方法，本发明可排除人体运动、肌肉活动对心率分析的影响，减少了硬件成本，可以灵活地选择关注频段区间，仅在关注频段区间内提高信号的频域精度，避免带来过大的计算和存储开销，选择信号频谱分析的技术路线，避免了时域信号波形特征点难找，算法参数具体值难确定的问题。
附图说明
图1为本发明的一种实时动态心率监测装置的较佳实施例的功能原理框图。
图2为本发明的一种实时动态心率监测装置的监测方法的较佳实施例的流程图。
具体实施方式
为使本发明的目的、技术方案及效果更加清楚、明确，以下对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明还提供了一种实时动态心率监测装置的较佳实施例的功能原理图，如图1所示，其中，装置包括基线漂移消除模块100，带通滤波模块200，频域分析模块300，心率频点选择模块400；
所述基线漂移消除模块100用于消除导致脉搏波信号的基准线出现漂移的干扰信号；所述带通滤波模块200用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号；所述频域分析模块300用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱；所述心率频点选择模块400用于获取心率对应的频点并输出频点；所述基线漂移消除模块100与所述带通滤波模块200连接，所述带通滤波模块200与所述频域分析模块300连接，所述频域分析模块300还与所述心率频点选择模块400连接。
具体实施时，所述基线漂移消除模块具体包括基线漂移趋势项信号提取单元和信号线性叠加单元；
所述基线漂移趋势项信号提取单元用于获取与原脉搏波信号等长的基线漂移趋势项信号；
所述信号线性叠加单元用于将原脉搏波信号减掉基线漂移趋势项信号，得到去掉基线漂移的脉搏波信号；
所述基线漂移趋势项信号提取单元与所述信号线性叠加单元连接。
具体实施时，因为呼吸、肢体活动或运动会导致脉搏波信号的基准线呈现上下漂移的情况，本质上这种基线漂移是一种低频信号，会对关键信息的判断带来干扰。消除基线漂移模块的主要作用就是降低呼吸、运动对脉搏波信号带来的基线漂移干扰，避免对加速度计辅助信息的依赖。
消除基线漂移模块由两个功能单元组成：基线漂移趋势项信号提取单元，信号线性叠加单元。
基线漂移趋势项信号提取单元，具体的实施方式有两类，第一类是信号滤波方法，第二类是曲线拟合方法。信号滤波方法的具体实施方式为中值滤波或均值滤波，使用滑动窗口遍历脉搏波信号，计算窗口内的全部信号值的中值或均值，最终得到与原脉搏波信号等长的基线漂移趋势项信号。曲线拟合方法的具体实施方式，将基线漂移信号视作可以表示为一个高次多项式的时间函数，用脉搏波信号作为输入数据样本，采用非线性拟合方法得到高次多项式的各项系数，最终通过多项式函数计算得到与原脉搏波信号等长的基线漂移趋势项信号。
信号线性叠加单元，具体实施方式就是将原脉搏波信号，基线漂移趋势项信号使用行向量或者列向量存储，然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信号，最终得到去掉基线漂移的脉搏波信号。
进一步地，所述带通滤波模块具体包括滤波参数设置单元和滤波单元，
所述滤波参数设置单元用于设置通带下限和上限频率、阻带下限和上限频率和上限频率、通带内衰减系数、阻带内衰减系数；
所述滤波单元用于获取滤波模块阶数和系数后，将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波，获取去除噪声后的脉搏波信号；
所述滤波参数设置单元和滤波单元连接。
具体实施时，大量医学实践表明，普遍情况下人体心率的范围是40次/分钟至220次/分钟，这表明心率信号的频率范围是0.6 Hz至3.7 Hz。因此，可以将0 Hz至0.6 Hz，以及3.7 Hz以上的信号频段视作噪声频段，这些噪声包括传感器电路导致的毛刺噪声、人体活动引入的噪声等等。带通滤波器模块的主要作用就是仅保留属于心率频段的信号分量，同时消除带外噪声对于发现心率信号的干扰。
带通滤波器模块，具体的实施方式可以采用巴特沃斯滤波器，通过设置通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数，参照巴特沃斯滤波器设计流程计算得到滤波器阶数和系数。将去掉基线漂移的脉搏波信号输入带通滤波器模块，最终得到去除噪声的脉搏波信号。
在具体实施时，不排除存在极端情况下（疾病、剧烈运动）的心率信号，为保证全面地监测到心率变化，可以适当增大心率频段的动态范围，例如，设置成0.5 Hz至4 Hz；带通滤波器的实施方式除巴特沃斯滤波器外，还可以选择切比雪夫滤波器等。
进一步的实施例中，所述频域分析模块具体包括特定频段区间设置单元、频点功率计算单元和信号频谱获取单元；
所述特定频段区间设置单元用于设置特定频段区间的起始频点、结束频点、频点细分数目；
所述频点功率计算单元用于对于特定频段区间的每个频点，分别用第一变换多项式计算频域信号的实部，用第二变换多项式计算频域信号的虚部，根据频域信号的实部和虚部计算得到脉搏波信号在频点的功率；
所述信号频谱获取单元用于获取特定频段区间内的脉搏波信号在每个频点的功率，叠加后生成脉搏波信号在整个关注频段区间的信号频谱；
所述频点功率计算单元分别与所述特定频段区间设置单元、所述信号频谱获取单元连接。
具体实施时，去除噪声的脉搏波信号是时域信号，为了更精准的提取心率信息，需要采用信号变换方法获得脉搏波信号对应的频域信号，并且保证信号频谱有足够高的频域精度。经典的时频信号变换方法是快速傅里叶变换即FFT，它的缺点是为提高信号频域精度需要付出嵌入式设备难以承担的计算和存储开销。高精度频域分析模块的主要作用是，可以灵活地选择关注频段区间，仅在关注频段区间内提高信号频域精度，避免带来过大的计算和存储开销。
高精度频域分析模块，具体的实施方式为，设置关注频段区间的起始频点f1、结束频点f2和频点细分数目Nf。
用q(n)，n=0，1，…，N-1表示数据长度为N的脉搏波信号序列。对于关注频段区间内的每个频点f，使用下面的变换多项式计算频域信号实部，
Qr = q(0) + q(1) * cos(f) + q(2) * T(2) + … + q(N-1) * T(N-1)
和使用下面的变换多项式计算频域信号虚部，
Qi= -q(1) * sin(f) – q(2) * sin(f) * U(2) - … - q(N-1) * sin(f) * U(N-1)
T(n)为第一变换多项式，U(n)为第二变换多项式，第一变换多项式和所述第二变换多项式为正交多项式。所述第一变换多项式为勒让德多项式，雅可比多项式，拉盖尔多项式，切比雪夫多项式，埃尔米特多项式中的一种；所述第二变换多项式为正交多项式为勒让德多项式，雅可比多项式，拉盖尔多项式，切比雪夫多项式，埃尔米特多项式中的一种。
根据Qr和Qi，可以计算得到脉搏波信号在频点f的功率。
当关注频段区间内的每个频点都经过以上计算后，则可以得到脉搏波信号在整个关注频段区间的信号频谱。高频域精度的量化指标与频点细分数目Nf有关，频点细分数目越大频域精度越高。
进一步地，所述心率频点选择模块具体包括心率频点选择单元和心率频点动态跟踪单元，
所述心率频点选择单元用于在初始状态下获取脉搏波信号频谱中谱峰所在的频点位置作为心率频点；
所述心率频点动态跟踪单元用于在预先设置的固定跟踪周期内，以上一跟踪周期输出的心率频点为中心，在中心左右一个特定宽度的频谱范围内进行高精度的频域分析，选择谱峰所在的频点位置作为心率频点输出；
所述心率频点选择单元与所述心率频点动态跟踪单元连接。
具体实施时，得到脉搏波信号在关注频段区间的信号频谱后，并不能很直接地计算心率值，而是需要通过分析判断发现代表心率的频点，才能计算心率值，这是心率频点选择模块的作用。
心率频点选择模块，具体实施方式，分为两个阶段。第一阶段是，初始状态精确心率频点选择；第二阶段是心率频点动态跟踪。
初始状态精确心率频点选择，具体实施方式是，保证人体处于安静状态下，通过消除基线漂移模块、带通滤波器模块和高精度频域分析模块得到脉搏波信号频谱，其中谱峰所在频点位置就是代表心率的频点，可以计算出心率值。
心率频点动态跟踪，具体实施方式是，以固定的跟踪周期进行循环，在每个跟踪周期内，通过消除基线漂移模块和带通滤波器模块得到去除噪声的脉搏波信号，以上一跟踪周期输出的心率频点为中心，在中心左右一个狭小宽度的频谱范围内进行高精度的频域分析，选择谱峰所在的频点位置作为心率频点输出。
本发明还提供了一种基于所述的实时动态心率监测装置的监测方法的较佳实施例，如图2所示，其中，方法包括：
步骤S100、获取脉搏波信号，通过基线漂移消除模块消除基线漂移干扰信号；
步骤S200、带通滤波模块对基线漂移消除模块输出的信号进行滤波，保留属于心率频段的信号分量；
步骤S300、频域分析模块获取脉搏波信号在预先设置的特定频段区间内的信号频谱；
步骤S400、心率频点选择模块根据信号频率中获取心率对应的频点并输出。
具体实施时，通过多种滤波器技术去除信号采集过程中混入的高频噪声，去除由肌肉抖动和呼吸引起的信号基线漂移，采用一种比FFT计算复杂度小的高精度信号频域分析方法，从信号频谱中筛选定位代表心率的频点，对心率进行实时动态地监测。减少了加速度计导致的硬件成本和资源上的额外开销。第二，避免了时域信号波形特征点难找，算法参数具体值难确定的问题。第三，算法的计算复杂度是普通嵌入式模块可以承担的。具体的监测方法如上监测装置的具体实施例所述。
进一步地，所述步骤S100具体包括：
步骤S101、通过信号滤波法或曲线拟合法获取与原脉搏波信号等长的基线漂移趋势项信号；
步骤S102、将原脉搏波信号，基线漂移趋势项信号使用行向量或者列向量存储，然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信，得到去掉基线漂移的脉搏波信号。
具体实施时，其中所述步骤S101中信号滤波方法的具体实施方式为中值滤波或均值滤波，使用滑动窗口遍历脉搏波信号，计算窗口内的全部信号值的中值或均值，最终得到与原脉搏波信号等长的基线漂移趋势项信号。曲线拟合方法的具体实施方式，将基线漂移信号视作可以表示为一个高次多项式的时间函数，用脉搏波信号作为输入数据样本，采用非线性拟合方法得到高次多项式的各项系数，最终通过多项式函数计算得到与原脉搏波信号等长的基线漂移趋势项信号。
进一步的实施例中，所述步骤S200具体包括：
步骤S201、获取预先设置的通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数；
步骤S202、获取滤波模块阶数和系数后，将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波，获取去除噪声后的脉搏波信号。
具体实施时，采用巴特沃斯滤波器，通过设置通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数，参照巴特沃斯滤波器设计流程计算得到滤波器阶数和系数。将去掉基线漂移的脉搏波信号输入带通滤波器模块，最终得到去除噪声的脉搏波信号。
在具体实施时，不排除存在极端情况下（疾病、剧烈运动）的心率信号，为保证全面地监测到心率变化，可以适当增大心率频段的动态范围，例如，设置成0.5 Hz至4 Hz；带通滤波器的实施方式除巴特沃斯滤波器外，还可以选择切比雪夫滤波器等。
综上所述，本发明提供了一种实时动态心率监测装置及监测方法，装置包括基线漂移消除模块，带通滤波模块，频域分析模块，心率频点选择模块；所述基线漂移消除模块用于消除导致脉搏波信号的基准线出现漂移的干扰信号；所述带通滤波模块用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号；所述频域分析模块用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱；所述心率频点选择模块用于获取心率对应的频点并输出频点。本发明排除人体运动、肌肉活动对心率分析的影响，减少了硬件成本，可以灵活地选择关注频段区间，仅在关注频段区间内提高信号的频域精度，避免带来过大的计算和存储开销。
应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (8)
1.一种实时动态心率监测装置，其特征在于，装置包括基线漂移消除模块，带通滤波模块，频域分析模块，心率频点选择模块；
所述基线漂移消除模块用于消除导致脉搏波信号的基准线出现漂移的干扰信号；
所述带通滤波模块用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号；
所述频域分析模块用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱；
所述心率频点选择模块用于获取心率对应的频点并输出频点；
所述基线漂移消除模块与所述带通滤波模块连接，所述带通滤波模块与所述频域分析模块连接，所述频域分析模块还与所述心率频点选择模块连接；
所述基线漂移消除模块具体包括基线漂移趋势项信号提取单元和信号线性叠加单元；
所述基线漂移趋势项信号提取单元用于获取与原脉搏波信号等长的基线漂移趋势项信号；
所述信号线性叠加单元用于将原脉搏波信号减掉基线漂移趋势项信号，得到去掉基线漂移的脉搏波信号；
所述基线漂移趋势项信号提取单元与所述信号线性叠加单元连接；
所述信号线性叠加单元用于将原脉搏波信号，基线漂移趋势项信号使用行向量或者列向量存储，然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信号，得到去掉基线漂移的脉搏波信号；
所述频域分析模块具体包括特定频段区间设置单元、频点功率计算单元和信号频谱获取单元；
所述特定频段区间设置单元用于设置特定频段区间的起始频点、结束频点、频点细分数目；
所述频点功率计算单元用于对于特定频段区间的每个频点，分别用第一变换多项式计算频域信号的实部，用第二变换多项式计算频域信号的虚部，根据频域信号的实部和虚部计算得到脉搏波信号在频点f的功率；
用q(n)，n＝0，1，……，N-1表示数据长度为N的脉搏波信号序列，
所述用第一变换多项式计算频域信号的实部具体为：
Qr＝q(0)+q(1)*cos(f)+q(2)*T(2)+…+q(N-1)*T(N-1)；
所述用第二变换多项式计算频域信号的虚部具体为：
Qi＝-q(1)*sin(f)–q(2)*sin(f)*U(2)-…-q(N-1)*sin(f)*U(N-1)；
T(n)为第一变换多项式，U(n)为第二变换多项式，f为频点；
所述信号频谱获取单元用于获取特定频段区间内的脉搏波信号在每个频点的功率，叠加后生成脉搏波信号在整个关注频段区间的信号频谱；
所述心率频点选择模块具体用于两个阶段，第一阶段是初始状态精确心率频点选择，第二阶段是心率频点动态跟踪；
所述初始状态精确心率频点选择具体为：保证人体处于安静状态下，通过消除基线漂移模块、带通滤波器模块和频域分析模块得到脉搏波信号频谱，其中谱峰所在频点位置就是代表心率的频点。
2.根据权利要求1所述的实时动态心率监测装置，其特征在于，所述带通滤波模块具体包括滤波参数设置单元和滤波单元，
所述滤波参数设置单元用于设置通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数；
所述滤波单元用于获取滤波模块阶数和系数后，将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波，获取去除噪声后的脉搏波信号；
所述滤波参数设置单元和滤波单元连接。
3.根据权利要求2所述的实时动态心率监测装置，其特征在于，
所述频点功率计算单元分别与所述特定频段区间设置单元、所述信号频谱获取单元连接。
4.根据权利要求3所述的实时动态心率监测装置，其特征在于，所述心率频点选择模块具体包括心率频点选择单元和心率频点动态跟踪单元，
所述心率频点选择单元用于在初始状态下获取脉搏波信号频谱中谱峰所在的频点位置作为心率频点；
所述心率频点动态跟踪单元用于在预先设置的固定跟踪周期内，以上一跟踪周期输出的心率频点为中心，在中心左右一个特定宽度的频谱范围内进行高精度的频域分析，选择谱峰所在的频点位置作为心率频点输出；
所述心率频点选择单元与所述心率频点动态跟踪单元连接。
5.根据权利要求4所述的实时动态心率监测装置，其特征在于，所述第一变换多项式和所述第二变换多项式为正交多项式。
6.根据权利要求5所述的实时动态心率监测装置，其特征在于，所述第一变换多项式为勒让德多项式，雅可比多项式，拉盖尔多项式，切比雪夫多项式，埃尔米特多项式中的一种；
所述第二变换多项式为勒让德多项式，雅可比多项式，拉盖尔多项式，切比雪夫多项式，埃尔米特多项式中的一种。
7.一种基于权利要求1所述的实时动态心率监测装置的监测方法，其特征在于，方法包括步骤：
A、获取脉搏波信号，通过基线漂移消除模块消除基线漂移干扰信号；
B、带通滤波模块对基线漂移消除模块输出的信号进行滤波，保留属于心率频段的信号分量；
C、频域分析模块获取脉搏波信号在预先设置的特定频段区间内的信号频谱；
D、心率频点选择模块根据信号频率中获取心率对应的频点并输出；
所述步骤A具体包括：
A1、通过信号滤波法或曲线拟合法获取与原脉搏波信号等长的基线漂移趋势项信号；
A2、将原脉搏波信号，基线漂移趋势项信号使用行向量或者列向量存储，然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信，得到去掉基线漂移的脉搏波信号；
所述频域分析模块用于设置特定频段区间的起始频点、结束频点、频点细分数目；对于特定频段区间的每个频点，分别用第一变换多项式计算频域信号的实部，用第二变换多项式是计算频域信号的虚部，根据频域信号的实部和虚部计算得到脉搏波信号在频点f的功率；获取特定频段区间内的脉搏波信号在每个频点的功率，叠加后生成脉搏波信号在整个关注频段区间的信号频谱；
用q(n)，n＝0，1，……，N-1表示数据长度为N的脉搏波信号序列，
所述用第一变换多项式计算频域信号的实部具体为：
Qr＝q(0)+q(1)*cos(f)+q(2)*T(2)+…+q(N-1)*T(N-1)；
所述用第二变换多项式计算频域信号的虚部具体为：
Qi＝-q(1)*sin(f)–q(2)*sin(f)*U(2)-…-q(N-1)*sin(f)*U(N-1)
T(n)为第一变换多项式，U(n)为第二变换多项式，f为频点；
所述心率频点选择模块具体用于两个阶段，第一阶段是初始状态精确心率频点选择，第二阶段是心率频点动态跟踪；
所述初始状态精确心率频点选择具体为：保证人体处于安静状态下，通过消除基线漂移模块、带通滤波器模块和频域分析模块得到脉搏波信号频谱，其中谱峰所在频点位置就是代表心率的频点。
8.根据权利要求7所述的监测方法，其特征在于，所述步骤B具体包括：
B1、获取预先设置的通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数；
B2、获取滤波模块阶数和系数后，将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波，获取去除噪声后的脉搏波信号。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910951368.9A CN110876615B (zh) | 2016-05-04 | 2016-05-04 | 一种实时动态心率监测装置及监测方法 |
| CN201610288485.8A CN105816165B (zh) | 2016-05-04 | 2016-05-04 | 一种实时动态心率监测装置及监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610288485.8A CN105816165B (zh) | 2016-05-04 | 2016-05-04 | 一种实时动态心率监测装置及监测方法 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910951368.9A Division CN110876615B (zh) | 2016-05-04 | 2016-05-04 | 一种实时动态心率监测装置及监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105816165A CN105816165A (zh) | 2016-08-03 |
| CN105816165B true CN105816165B (zh) | 2019-12-13 |
Family
ID=56527978
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610288485.8A Expired - Fee Related CN105816165B (zh) | 2016-05-04 | 2016-05-04 | 一种实时动态心率监测装置及监测方法 |
| CN201910951368.9A Active CN110876615B (zh) | 2016-05-04 | 2016-05-04 | 一种实时动态心率监测装置及监测方法 |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910951368.9A Active CN110876615B (zh) | 2016-05-04 | 2016-05-04 | 一种实时动态心率监测装置及监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN105816165B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107088062B (zh) * | 2017-05-08 | 2023-09-01 | 歌尔科技有限公司 | 一种心率模组、采集心率的电子装置和心率采集方法 |
| CN109924966A (zh) * | 2017-12-19 | 2019-06-25 | 郝振龙 | 一种用于运动时测量心率的智能手环装置 |
| CN110801210A (zh) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | 一种脉搏波信号的滤波方法、装置、可读介质及电子设备 |
| CN113616217B (zh) * | 2021-10-12 | 2021-12-21 | 深圳市倍轻松科技股份有限公司 | 基线漂移曲线的生成方法和装置 |
| CN117951505B (zh) * | 2024-03-27 | 2024-06-04 | 剑博微电子(南京)有限公司 | 一种物联网医疗芯片的降噪方法及*** |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102028457A (zh) * | 2010-11-24 | 2011-04-27 | 北京麦邦光电仪器有限公司 | 脉率测量方法及指环式脉率测量仪 |
| CN105286845A (zh) * | 2015-11-29 | 2016-02-03 | 浙江师范大学 | 一种适用于可穿戴式心率测量设备的运动噪声消除方法 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI261415B (en) * | 2003-05-21 | 2006-09-01 | Leadtek Research Inc | Pulse interval to voltage converter and converting method thereof |
| CN101194834A (zh) * | 2006-12-05 | 2008-06-11 | 重庆博恩富克医疗设备有限公司 | 一种生物电阻抗测量方法及测量装置 |
| CN101385644B (zh) * | 2008-07-25 | 2011-06-15 | 沈阳中国医科大学医疗器械研制中心（有限公司） | 12导联无线远程心电监护*** |
| WO2010140241A1 (ja) * | 2009-06-04 | 2010-12-09 | 富士通株式会社 | 覚醒度判定装置、覚醒度判定方法および覚醒度判定プログラム |
| US9504390B2 (en) * | 2011-03-04 | 2016-11-29 | Globalfoundries Inc. | Detecting, assessing and managing a risk of death in epilepsy |
| CN103717129A (zh) * | 2011-05-24 | 2014-04-09 | 加利福尼亚大学董事会 | 脑磁图源成像 |
| CN102217931A (zh) * | 2011-06-09 | 2011-10-19 | 李红锦 | 一种心率变异性特征参数的获取方法和装置 |
| CN102984103B (zh) * | 2011-09-07 | 2015-08-05 | 华为技术有限公司 | 扩频***中的信号处理方法及装置 |
| CN103417206B (zh) * | 2012-05-22 | 2015-08-26 | 中国科学院深圳先进技术研究院 | Ecg去除工频干扰的方法及*** |
| JP2014122875A (ja) * | 2012-11-26 | 2014-07-03 | Canon Inc | 層状物体の測定装置および方法 |
| CN103654758A (zh) * | 2013-12-23 | 2014-03-26 | 韩山师范学院 | 一种抗干扰的心率检测方法 |
| US9936886B2 (en) * | 2014-06-09 | 2018-04-10 | Stmicroelectronics S.R.L. | Method for the estimation of the heart-rate and corresponding system |
| CN104095630B (zh) * | 2014-07-29 | 2016-02-10 | 杭州电子科技大学 | 一种基于脑电相位同步的疲劳检测方法 |
| CN104125579B (zh) * | 2014-08-07 | 2017-07-11 | 桂林电子科技大学 | 一种基于时域能量与频域谱熵的频谱感知方法及装置 |
| CN104569902B (zh) * | 2014-11-21 | 2017-04-19 | 国家电网公司 | 一种数字式电能表功耗测量装置及方法 |
| CN104490373B (zh) * | 2014-12-17 | 2016-12-07 | 辛勤 | 脉搏信号的判断方法、判断装置以及生理参数测量设备 |
| CN104655929B (zh) * | 2015-01-04 | 2017-11-21 | 中国科学院物理研究所 | 一种时域信号的数字时频测量方法及相应的目标识别方法 |
-
2016
- 2016-05-04 CN CN201610288485.8A patent/CN105816165B/zh not_active Expired - Fee Related
- 2016-05-04 CN CN201910951368.9A patent/CN110876615B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102028457A (zh) * | 2010-11-24 | 2011-04-27 | 北京麦邦光电仪器有限公司 | 脉率测量方法及指环式脉率测量仪 |
| CN105286845A (zh) * | 2015-11-29 | 2016-02-03 | 浙江师范大学 | 一种适用于可穿戴式心率测量设备的运动噪声消除方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110876615B (zh) | 2022-08-02 |
| CN105816165A (zh) | 2016-08-03 |
| CN110876615A (zh) | 2020-03-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105816165B (zh) | 一种实时动态心率监测装置及监测方法 | |
| CN113349752B (zh) | 一种基于传感融合的可穿戴设备实时心率监测方法 | |
| US9247911B2 (en) | Devices and methods for real-time denoising of electrocardiograms | |
| US9675274B2 (en) | Systems and methods for determining respiration information from a photoplethysmograph | |
| Chung et al. | Deep learning for heart rate estimation from reflectance photoplethysmography with acceleration power spectrum and acceleration intensity | |
| CN105832289A (zh) | 用于使用希尔伯特变换估计生物生理速率的设备和方法 | |
| Yang et al. | Estimation and validation of arterial blood pressure using photoplethysmogram morphology features in conjunction with pulse arrival time in large open databases | |
| WO2000022408A9 (en) | Method, apparatus and system for removing motion artifacts from measurements of bodily parameters | |
| US20190298209A1 (en) | Heartbeat detection | |
| EP3473171A1 (en) | Method for calculating physiological parameters and corresponding medical equipment | |
| CN104757955A (zh) | 一种基于脉搏波的人体血压预测方法 | |
| US20140081088A1 (en) | Computer-implemented method for determining physical movements of a body organ | |
| CN112089405B (zh) | 一种脉搏波特征参数测量及显示装置 | |
| Hong et al. | Aging index using photoplethysmography for a healthcare device: comparison with brachial-ankle pulse wave velocity | |
| EP3292813A1 (en) | Method and device for processing bio-signals | |
| CN106539580B (zh) | 一种自主神经***动态变化的连续监测方法 | |
| Kumar et al. | STSR: spectro-temporal super-resolution analysis of a reference signal less photoplethysmogram for heart rate estimation during physical activity | |
| Zheng et al. | Adaptive notch-filtration to effectively recover photoplethysmographic signals during physical activity | |
| De Giovanni et al. | Ultra-low power estimation of heart rate under physical activity using a wearable photoplethysmographic system | |
| CN105852863A (zh) | 一种呼吸率测量方法及装置 | |
| Guo et al. | An effective photoplethysmography heart rate estimation framework integrating two-level denoising method and heart rate tracking algorithm guided by finite state machine | |
| Amhia et al. | Stability and Phase Response Analysis of Optimum Reduced‐Order IIR Filter Designs for ECG R‐Peak Detection | |
| Gupta et al. | Exploiting tunable Q-factor wavelet transform domain sparsity to denoise wrist PPG signals | |
| CN116502124A (zh) | 脉搏波信号的质量评估方法及装置、计算机可读存储介质 | |
| CN116028809A (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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20191108 Address after: Room 202, 2f, building 5, 1180 Xingxian Road, Jiading District, Shanghai 200000 Applicant after: Baomai (Shanghai) Information Technology Co.,Ltd. Address before: 200810 Shanghai City, north of the city of Jiading District Road No. 235, building 2, floor 2 Applicant before: Shanghai Tiezhuo Information Technology Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xu Weiwei Inventor after: Deng Hanlin Inventor after: Zhang Lei Inventor after: Chen Huanting Inventor after: Sun Rongxue Inventor before: Xu Weiwei Inventor before: Deng Hanlin |
|
| CB03 | Change of inventor or designer information | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191213 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |