CN105943012B - 一种去除光电容积脉搏波信号中运动噪声的心率测量方法 - Google Patents

一种去除光电容积脉搏波信号中运动噪声的心率测量方法 Download PDF

Info

Publication number
CN105943012B
CN105943012B CN201610293793.XA CN201610293793A CN105943012B CN 105943012 B CN105943012 B CN 105943012B CN 201610293793 A CN201610293793 A CN 201610293793A CN 105943012 B CN105943012 B CN 105943012B
Authority
CN
China
Prior art keywords
heart rate
spectrum matrix
signals
photoplethysmographic
motion
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
CN201610293793.XA
Other languages
English (en)
Other versions
CN105943012A (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 CN201610293793.XA priority Critical patent/CN105943012B/zh
Publication of CN105943012A publication Critical patent/CN105943012A/zh
Application granted granted Critical
Publication of CN105943012B publication Critical patent/CN105943012B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • 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/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
    • 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/02438—Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1118—Determining activity level
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802—Sensor mounted on worn items
    • A61B5/681—Wristwatch-type devices
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813—Specially adapted to be attached to a specific body part
    • A61B5/6824—Arm or wrist
    • 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
    • A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • 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
    • 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)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明公开了一种去除光电容积脉搏波信号中运动噪声的心率测量方法，能有效的降低运动噪声对心率测量的影响。该方法中，选用脉搏血氧仪和三轴加速度计采集用户在同时间段内的多个光电容积脉搏波信号及运动加速度信号；利用多个光电容积脉搏波信号和运动加速度信号构成频谱矩阵；然后依据频谱矩阵的整体稀疏和行稀疏的结构特征构建频谱矩阵分解模型，并通过近端梯度加速优化算法求解频谱矩阵分解模型的最优解；最后，利用谱峰跟踪方法准确定位心率频率点位置。本发明能有效地消除光电容积脉搏波信号中的运动噪声，实现基于可穿戴设备的心率准确测量。

Description

一种去除光电容积脉搏波信号中运动噪声的心率测量方法
技术领域
本发明涉及运动状态的心率测量，尤其涉及一种去除光电容积脉搏波信号中运动噪声的心率测量方法。
背景技术
随着电子技术微型化的发展以及人们对于健康的重视，使用光电容积脉搏波信号进行心率测量的方法引起了学术界和工业界的高度关注。光电容积脉搏波信号是借助一种无创检测方法从皮肤表面提取的生物信号，它的信号强度弱、易受噪声干扰。在运动状态下，存在由于组织干扰、静脉血容量以及光程变化产生的运动噪声干扰，并且运动噪声的频率在很多情况下会与心率的频率十分接近，因此使得心率测量变得较为困难。
针对光电容积脉搏波信号中运动噪声干扰的消除，研究人员已进行了相关研究并取得了一些研究成果。例如，发明专利“一种适用于可穿戴式心率监测设备的运动和噪声干扰消除方法”中将归一化最小均方差(NLMS)自适应滤波器和Mallat算法融合，用于消除运动部分噪声干扰，其中自适应滤波器的运动参考信号选用由三轴加速度计及陀螺仪组成的六轴加速度信号；发明专利“基于光电容积描记和谱分析的运动状态心率监测方法”中首先对源信号进行零相位基线滤波，再用滑动窗方法分割处理后的信号并使用短时傅立叶变换得到频谱，然后利用谱峰追踪方法找到代表心率数据的谱峰，最后用循环移动平均滤波器对心率数据进行处理，得出最终心率结果；张坤等人在“运用小波模极大值滤波算法消除光电容积脉搏波中的运动干扰”中先利用朗伯-比尔定理和指尖半球形模型分析运动干扰和正常信号之间的关系，然后引入小波模极大值重构滤波算法去除信号中的运动干扰。然而上述算法主要针对缓和或者不剧烈的运动，比如手移动、走路、慢跑(速度低于8km/h)。
针对剧烈运动的情况，发明专利“一种适用于心率信号的运动噪声检测方法”中提出了联合稀疏谱重构模型，该模型中对整个频谱矩阵进行了行稀疏和全局稀疏的限制，并通过不精确增广拉格朗日乘子法求解该模型的最优解；本发明提出的频谱矩阵分解模型是将整体的频谱矩阵分解成行稀疏限制的运动噪声信号频谱矩阵和全局稀疏限制的真实光电容积脉搏波信号频谱矩阵，并通过近端梯度加速算法求解该模型的最优解。发明专利“一种适用于可穿戴式心率测量设备的运动噪声消除方法”中先利用谱减法去除光电容积脉搏波信号频谱中的运动噪声，然后根据谱峰跟踪机制定位心率频率点位置；而本发明中含有运动噪声的光电容积脉搏波信号经过频谱矩阵分解模型处理后，获得真实光电容积脉搏波信号频谱，可直接利用谱峰跟踪方法定位心率频率点位置，减少了计算复杂度，缩短了运行时间。
目前部分去噪方法中仅使用单个光电容积脉搏波信号，当该光电容积脉搏波信号频谱中不存在心率所对应的谱峰时，只能预测心率，易影响心率测量的准确性。本发明引入多个光电容积脉搏波信号，使各个光电容积脉搏波信号之间“相互辅助”，弥补了基于单个光电容积脉搏波信号去噪方法的不足。因此，本发明有效地消除了光电容积脉搏波信号中的运动噪声，实现了基于可穿戴设备的实时心率的准确测量。
发明内容
本发明所要解决的技术问题是如何在运动噪声非常强烈的情况下提供一种有效去除运动噪声的方法，以达到准确获得实时心率的目的。
为了解决上述技术问题，本发明提供了一种去除光电容积脉搏波信号中运动噪声的心率测量方法，包括信号采集、频谱矩阵分解模型和谱峰跟踪方法三个部分，其特征在于：
所述脉搏血氧仪和所述三轴加速计采集多个光电容积脉搏波信号及运动加速度信号；利用所述多个光电容积脉搏波信号和所述运动加速度信号构成频谱矩阵；然后，依据上述频谱矩阵的结构特征构建频谱矩阵分解模型，并求解所述频谱矩阵分解模型的最优解；最后，利用所述谱峰跟踪方法准确定位心率频率点位置；
该方法包括如下步骤：
所述脉搏血氧仪和所述三轴加速计在用户手腕处采集同一时间段内的多个光电容积脉搏波信号及运动加速度信号；对上述多个光电容积脉搏波信号和运动加速度信号进行预处理；
同时，利用所述多个光电容积脉搏波信号和所述运动加速度信号构成频谱矩阵；依据上述频谱矩阵的整体稀疏和行稀疏的结构特征建立所述频谱矩阵分解模型，并通过近端梯度加速优化算法求解所述频谱矩阵分解模型的最优解；
所述谱峰跟踪方法的各个子阶段对上述去噪后的多个光电容积脉搏波信号频谱进行处理，定位用户的心率频率点位置；
优选地，所述预处理的过程包括对所述多个光电容积脉搏波信号和所述运动加速度信号进行下采样处理以及对下采样后的上述信号进行带通滤波操作；
优选地，所述频谱矩阵由所述多个光电容积脉搏波信号和所述运动加速度信号构成；所述频谱矩阵分解模型的构造是依据所述频谱矩阵的全局稀疏和行稀疏的结构特征，并将所述频谱矩阵分解为运动噪声信号频谱矩阵和真实光电容积脉搏波信号频谱矩阵，所述频谱矩阵分解模型的目标函数如下：
Figure BSA0000129702270000031
其中，
Figure BSA0000129702270000032
分别用来约束运动噪声信号频谱矩阵行稀疏和真实光电容积脉搏波信号频谱矩阵全局稀疏，pi，j是运动噪声信号频谱矩阵P的第i行第j列元素，qi，j是真实光电容积脉搏波信号频谱矩阵Q的第i行第j列元素，λ1、λ2是用来权衡各项重要性的权值；约束条件中Y∈RM×H是一个观测矩阵，
Figure BSA0000129702270000033
是一个冗余离散傅里叶变换基；
优选地，所述近端梯度加速优化算法与传统的次梯度方法相比，具有更好的收敛速度，能以二次收敛速度获得全局最优解，即所述近端梯度加速优化算法经过m次迭代后可得到残差为O(1/m2)的最优解；
优选地，所述谱峰跟踪方法主要包括初始化、谱峰选择、谱峰预测和谱峰发现四个子阶段；所述初始化阶段是为了准确定位首个时间窗口的心率频率点位置；所述谱峰发现阶段根据上一时间窗口的心率历史信息，寻找当前时间窗口的心率频率点位置；所述谱峰预测阶段依赖前几个时间窗口的心率历史信息，预测当前时间窗口的心率频率点位置；所述谱峰发现阶段能有效地处理目标谱峰被跟丢的情况。
与现有技术相比，本发明提供的技术方案有效地消除了光电容积脉搏波信号中的运动噪声，并对心率所对应的谱峰进行了准确地跟踪，提高了心率测量的精度。
附图说明
图1为本发明实施例的去除光电容积脉搏波信号中运动噪声的心率测量方法的流程示意图；
图2为本发明实施例的谱峰跟踪方法的流程示意图。
具体实施方式
以下结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。
本发明的技术方案中，由于真实光电容积脉搏波信号的周期性和心率信号的周期性一致，故在频域表现出稀疏性；且运动加速度信号和运动噪声信号的强相关性在频率表现为光电容积脉搏波信号频谱的谱峰位置和运动加速度信号频谱的谱峰位置大部分相同，所以由多个光电容积脉搏波信号和运动加速度信号构成的频谱矩阵可分解为行稀疏限制的运动噪声信号频谱矩阵和全局稀疏限制的真实光电容积脉搏波信号频谱矩阵。同时，该方法提出谱峰跟踪方法，用于准确定位心率频率点位置。此技术方案有效地消除了光电容积脉搏波信号中的运动噪声，实现了基于可穿戴设备的实时心率的准确测量及计算。
实施例一、去除光电容积脉搏波信号中运动噪声的心率测量方法
图1为本实施例的去除光电容积脉搏波信号中运动噪声的心率测量方法的流程示意图，图2为本实施例的谱峰跟踪方法的流程示意图。
图1所示的本实施例，是去除光电容积脉搏波信号中运动噪声的心率测量方法的整体流程，主要包括如下步骤：
步骤S210，利用分布在不同位置的两个脉搏血氧仪采集两个通道的光电容积脉搏波信号(以下简称PPG1和PPG2)，再利用三轴加速度计采集同时间段内的三个通道的运动加速度信号；
步骤S220，上述原始信号的初始采样频率为125Hz，为减少计算量，需要对上述原始信号进行下采样至采样频率为25Hz的操作；
步骤S230，经下采样后的上述信号需要通过通带为0.4Hz-4Hz的二阶巴特沃斯滤波器进行滤波，以消除一定频率范围以外的运动噪声及其它噪声的干扰；
步骤S240，利用上述两个光电容积脉搏波信号和三个运动加速度信号构成频谱矩阵；
步骤S250，依据上述频谱矩阵的全局稀疏和行稀疏的结构特征，构造频谱矩阵分解模型；
本步骤中，典型地，公式(1)为频谱矩阵分解模型的目标函数：
Figure BSA0000129702270000051
其中，
Figure BSA0000129702270000052
分别用来约束运动噪声信号频谱矩阵行稀疏和真实光电容积脉搏波信号频谱矩阵全局稀疏，pi，j是运动噪声信号频谱矩阵P的第i行第j列元素，qi，j是真实光电容积脉搏波信号频谱矩阵Q的第i行第j列元素，λ1、λ2是用来权衡各项重要性的权值；约束条件中Y∈RM×H是一个观测矩阵，本实施例中H＝5，
Figure BSA0000129702270000053
是一个冗余离散傅里叶变换基；
步骤S260，通过近端梯度加速优化算法可求解上述频谱矩阵分解模型的目标函数的最优解；
本步骤中，典型地，上述频谱矩阵分解模型的目标函数由可微分经验代价函数部分L(P，Q)和凸非光滑正则化部分G(P，Q)组成，即：
Figure BSA0000129702270000054
G(P，Q)＝λ1||P||1，2+λ2||Q||1，1 (3)
利用复合梯度映射将上述目标函数构造成如下表达式：
Figure BSA0000129702270000055
其中，F(P，Q；Θ，Ψ)由正则化项G(P，Q)和L(P，Q)在点(Θ，Ψ)处的一阶泰勒展开式两部分组成，泰勒展开式的余项为(P，Q)和(Θ，Ψ)之间欧几里得距离的平方，
Figure BSA0000129702270000056
和
Figure BSA0000129702270000057
分别为L(P，Q)关于P和Q的偏导数，η是控制步长的参数；
为了得到上述目标函数的最优解，需要对公式(4)中各变量进行独立地并行更新，即各变量的更新公式分别为公式(5)、(6)、(7)、(8)：
Figure BSA0000129702270000061
Figure BSA0000129702270000062
Figure BSA0000129702270000063
Figure BSA0000129702270000064
步骤S270，经上述操作后得到两个去噪后的光电容积脉搏波信号频谱，再利用谱峰跟踪方法定位用户的心率频率点位置；
本步骤中，典型地，谱峰跟踪方法的具体步骤如图2所示：
步骤S310，初始化是指需要用户在最初的几秒内尽量减少手部运动，以保证初始心率频率点位置的准确性，本实施例中选择PPG1信号频谱中谱峰最大的位置作为对应心率的谱峰位置；
步骤S320，通过初始化确定首个时间窗口中心率所对应的谱峰位置后，再确定后续时间窗口中心率所对应谱峰的搜索范围；
步骤S330，判断确定的搜索范围内PPG1存在的谱峰数，当谱峰数在正常范围(本实施例中正常范围设定为0＜谱峰数≤4)内，跳转至步骤S340，启动谱峰选择；当出现无峰或者多峰情况，跳转至步骤S350，缩小心率所对应谱峰的搜索范围；
步骤S360，在缩小的搜索范围内判断PPG2存在的谱峰数，若谱峰数在正常范围内(本实施例中正常范围设定为0＜谱峰数≤4)或者出现多峰情况，则返回至步骤S340，启动谱峰选择；若仍出现无峰情况，则跳转至步骤S370，启动谱峰预测；
步骤S380，经过对多个时间窗口中心率所对应谱峰的连续跟踪后，有可能出现目标谱峰被跟丢的情况，此时跳转至步骤S390，启动谱峰发现。
本实施例中，无峰情况是指两个光电容积脉搏波信号频谱中都不存在心率所对应的谱峰或者只有一个光电容积脉搏波信号频谱中含有心率所对应的谱峰；多峰情况是指两个光电容积脉搏波信号频谱中心率所对应的谱峰位置附近存在多个其它谱峰或者只有一个光电容积脉搏波信号频谱中心率所对应的谱峰位置附近存在多个其它谱峰。
步骤S280，经上述步骤处理后，输出用户的实时心率。
本实施例中，采用了两个脉搏血氧仪和一个三轴加速度计，从用户手腕处采集同时间段内的两个光电容积脉搏波信号和运动加速度信号；利用采集的信号构成频谱矩阵，并依据频谱矩阵所具有的结构特征将运动噪声的去除过程建模为频谱矩阵分解模型，通过近端梯度加速优化算法求解该模型的最优解；最后，根据谱峰跟踪方法精确地定位心率频率点位置。此方法有效地消除了光电容积脉搏波信号中的运动噪声，并对心率所对应的谱峰进行了准确地跟踪，实现了基于可穿戴设备的实时心率测量及计算。
虽然本发明所揭露的实施方式如上，但上述内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。在不脱离本发明所揭露的精神及范围的前提下，可在实施的形式上及细节上作任何的修饰与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (5)

1.一种去除光电容积脉搏波信号中运动噪声的心率测量方法，包括信号采集、频谱矩阵分解模型和谱峰跟踪方法三个部分，其特征在于：
所述脉搏血氧仪和三轴加速计采集多个光电容积脉搏波信号及运动加速度信号；利用所述多个光电容积脉搏波信号和所述运动加速度信号构成频谱矩阵；然后，依据上述频谱矩阵的结构特征构建频谱矩阵分解模型，并求解所述频谱矩阵分解模型的最优解；最后，利用所述谱峰跟踪方法准确定位心率频率点位置；
该方法包括如下步骤，
所述脉搏血氧仪和三轴加速计在用户手腕处采集同一时间段内的多个光电容积脉搏波信号及运动加速度信号；对上述多个光电容积脉搏波信号和运动加速度信号进行预处理；
同时，利用所述多个光电容积脉搏波信号和所述运动加速度信号构成频谱矩阵；依据上述频谱矩阵的整体稀疏和行稀疏的结构特征建立所述频谱矩阵分解模型，并通过近端梯度加速优化算法求解所述频谱矩阵分解模型的最优解，经上述操作后得到两个去噪后的光电容积脉搏波信号频谱；
所述谱峰跟踪方法主要包括初始化、谱峰选择、谱峰预测和谱峰发现四个子阶段，各个子阶段对上述去噪后的多个光电容积脉搏波信号频谱进行处理，定位用户的心率频率点位置。
2.根据权利要求1所述的去除光电容积脉搏波信号中运动噪声的心率测量方法，其特征在于：
所述预处理的过程包括对所述多个光电容积脉搏波信号和所述运动加速度信号进行下采样处理以及对下采样后的上述信号进行带通滤波操作。
3.根据权利要求1所述的去除光电容积脉搏波信号中运动噪声的心率测量方法，其特征在于：
所述频谱矩阵由所述多个光电容积脉搏波信号和所述运动加速度信号构成；所述频谱矩阵分解模型的构造是依据所述频谱矩阵的全局稀疏和行稀疏的结构特征，并将所述频谱矩阵分解为运动噪声信号频谱矩阵和真实光电容积脉搏波信号频谱矩阵，所述频谱矩阵分解模型的目标函数如下，
Figure FSB0000185129510000011
s.t：Y＝Φ(P+Q)
其中，
Figure FSB0000185129510000021
分别用来约束运动噪声信号频谱矩阵行稀疏和真实光电容积脉搏波信号频谱矩阵全局稀疏，pi，j是运动噪声信号频谱矩阵P的第i行第j列元素，qi，j是真实光电容积脉搏波信号频谱矩阵Q的第i行第j列元素，λ1、λ2是用来权衡各项重要性的权值；约束条件中Y∈RM×H是一个观测矩阵，
Figure FSB0000185129510000022
是一个冗余离散傅里叶变换基。
4.根据权利要求1所述的去除光电容积脉搏波信号中运动噪声的心率测量方法，其特征在于：
所述近端梯度加速优化算法与传统的次梯度方法相比，具有更好的收敛速度，能以二次收敛速度获得全局最优解，即所述近端梯度加速优化算法经过m次迭代后可得到残差为O(1/m2)的最优解。
5.根据权利要求1所述的去除光电容积脉搏波信号中运动噪声的心率测量方法，其特征在于：
所述谱峰跟踪方法主要包括初始化、谱峰选择、谱峰预测和谱峰发现四个子阶段；所述初始化阶段是为了准确定位首个时间窗口的心率频率点位置；所述谱峰发现阶段根据上一时间窗口的心率历史信息，寻找当前时间窗口的心率频率点位置；所述谱峰预测阶段依赖前几个时间窗口的心率历史信息，预测当前时间窗口的心率频率点位置；所述谱峰发现阶段能有效地处理目标谱峰被跟丢的情况。
CN201610293793.XA 2016-04-30 2016-04-30 一种去除光电容积脉搏波信号中运动噪声的心率测量方法 Active CN105943012B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610293793.XA CN105943012B (zh) 2016-04-30 2016-04-30 一种去除光电容积脉搏波信号中运动噪声的心率测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610293793.XA CN105943012B (zh) 2016-04-30 2016-04-30 一种去除光电容积脉搏波信号中运动噪声的心率测量方法

Publications (2)

Publication Number Publication Date
CN105943012A CN105943012A (zh) 2016-09-21
CN105943012B true CN105943012B (zh) 2020-06-16

Family

ID=56914028

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610293793.XA Active CN105943012B (zh) 2016-04-30 2016-04-30 一种去除光电容积脉搏波信号中运动噪声的心率测量方法

Country Status (1)

Country Link
CN (1) CN105943012B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019006631A1 (zh) * 2017-07-03 2019-01-10 深圳市汇顶科技股份有限公司 质量评估方法及装置、模型建立方法及模块、穿戴设备
WO2019036908A1 (zh) * 2017-08-22 2019-02-28 深圳市汇顶科技股份有限公司 心率检测方法及装置、电子终端
CN108937878B (zh) * 2018-06-06 2020-10-02 北京邮电大学 一种脉搏波信号运动噪声消除的方法
CN109222948B (zh) * 2018-09-17 2021-07-13 歌尔科技有限公司 消除运动干扰噪声的方法、装置、电子设备和存储介质
CN109222949B (zh) * 2018-10-12 2021-07-09 杭州士兰微电子股份有限公司 心率检测方法和心率检测装置
CN109875526B (zh) * 2018-12-28 2022-10-11 北京津发科技股份有限公司 一种基于压力与反射式综合测量脉搏装置
CN114305355B (zh) * 2022-01-05 2023-08-22 北京科技大学 基于毫米波雷达的呼吸心跳检测方法、系统及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105286846B (zh) * 2015-11-29 2018-03-02 浙江师范大学 一种适用于心率信号的运动噪声检测方法
CN205144548U (zh) * 2015-11-29 2016-04-13 浙江师范大学 适用于心率信号的运动噪声检测系统
CN105286845A (zh) * 2015-11-29 2016-02-03 浙江师范大学 一种适用于可穿戴式心率测量设备的运动噪声消除方法

Also Published As

Publication number Publication date
CN105943012A (zh) 2016-09-21

Similar Documents

Publication Publication Date Title
CN105943012B (zh) 一种去除光电容积脉搏波信号中运动噪声的心率测量方法
Zhang Heart rate monitoring from wrist-type photoplethysmographic (PPG) signals during intensive physical exercise
CN113349752B (zh) 一种基于传感融合的可穿戴设备实时心率监测方法
CN103034837B (zh) 特征参数与脉象要素的关联
CN111310570B (zh) 一种基于vmd和wpd的脑电信号情感识别方法及系统
CN109907752A (zh) 一种去除运动伪影干扰与心电特征检测的心电诊断与监护方法及系统
CN109875543B (zh) 用于可穿戴心率监测设备的多种健身运动状态下的心率估计方法及装置
Bashar et al. A machine learning approach for heart rate estimation from PPG signal using random forest regression algorithm
CN105286846B (zh) 一种适用于心率信号的运动噪声检测方法
CN111481192B (zh) 一种基于改进U-Net的心电信号R波检测方法
CN105997043B (zh) 一种基于腕式可穿戴设备的脉率提取方法
CN106137167A (zh) 一种基于光电容积脉搏波信号的运动噪声检测方法
Sun et al. PPG signal motion artifacts correction algorithm based on feature estimation
CN106264505A (zh) 一种基于支持向量机的心率谱峰选择方法
CN115836846B (zh) 一种基于自监督迁移学习的无创血压估计方法
Chang et al. Deepheart: accurate heart rate estimation from PPG signals based on deep learning
CN114052693A (zh) 心率分析方法、装置及设备
Brophy et al. An interpretable machine vision approach to human activity recognition using photoplethysmograph sensor data
Zheng et al. Adaptive notch-filtration to effectively recover photoplethysmographic signals during physical activity
CN114027804A (zh) 一种脉象诊断的方法、装置及可读存储介质
Bieri et al. Beliefppg: Uncertainty-aware heart rate estimation from ppg signals via belief propagation
CN111278353A (zh) 一种生命体征信号噪声的检测方法与系统
Pankaj et al. Blood pressure estimation and classification using a reference signal-less photoplethysmography signal: a deep learning framework
CN113616194B (zh) 一种监测手部震颤频率和强度的装置及方法
Hao et al. A novel sleep staging algorithm based on hybrid neural network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant