[go: up one dir, main page]

CN105816165A - Real-time dynamic heart rate monitoring device and monitoring method - Google Patents

Real-time dynamic heart rate monitoring device and monitoring method Download PDF

Info

Publication number
CN105816165A
CN105816165A CN201610288485.8A CN201610288485A CN105816165A CN 105816165 A CN105816165 A CN 105816165A CN 201610288485 A CN201610288485 A CN 201610288485A CN 105816165 A CN105816165 A CN 105816165A
Authority
CN
China
Prior art keywords
heart rate
signal
frequency
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.)
Granted
Application number
CN201610288485.8A
Other languages
Chinese (zh)
Other versions
CN105816165B (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.)
Baomai Shanghai Information Technology Co ltd
Original Assignee
Shanghai Tiezhuo Information Technology Co ltd
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 Shanghai Tiezhuo Information Technology Co ltd filed Critical Shanghai Tiezhuo Information Technology Co ltd
Priority to CN201610288485.8A priority Critical patent/CN105816165B/en
Priority to CN201910951368.9A priority patent/CN110876615B/en
Publication of CN105816165A publication Critical patent/CN105816165A/en
Application granted granted Critical
Publication of CN105816165B publication Critical patent/CN105816165B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring pulse rate or heart rate
    • A61B5/0245Measuring pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analogue processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details 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)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Physiology (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

The invention discloses a real-time dynamic heart rate monitoring device and a monitoring method, wherein the device comprises a baseline drift elimination module, a band-pass filtering module, a frequency domain analysis module and a heart rate frequency point selection module; the baseline wander elimination module is used for eliminating interference signals causing wander of a baseline line of the pulse wave signals; the band-pass filtering module is used for acquiring frequency signal components belonging to a heart rate frequency band and eliminating noise signals outside the heart rate frequency band; the frequency domain analysis module is used for acquiring a signal frequency spectrum of the pulse wave signal in a preset specific frequency band interval; the heart rate frequency point selection module is used for acquiring a frequency point corresponding to the heart rate and outputting the frequency point. The invention eliminates the influence of human body movement and muscle activity on heart rate analysis, reduces hardware cost, can flexibly select the concerned frequency band interval, improves the frequency domain precision of signals only in the concerned frequency band interval, and avoids causing excessive calculation and storage expenses.

Description

一种实时动态心率监测装置及监测方法A real-time dynamic heart rate monitoring device and monitoring method

技术领域 technical field

本发明涉及电子领域,尤其涉及一种实时动态心率监测装置及监测方法。 The invention relates to the field of electronics, in particular to a real-time dynamic heart rate monitoring device and a monitoring method.

背景技术 Background technique

随着移动互联网技术在医疗健康领域的应用与发展,在市场上出现了智能手表、智能手环、智能腕带等大量形态各异,具备测量心率、血压、血氧浓度、呼吸频率等医用生理参数的可穿戴式移动健康产品。心率定义为人体每分钟的心跳次数,是评估人健康状态的一个重要医用常规生理参数。心率监测对于疾病风险预警、病症诊断、年度例行体检具有十分重要的意义。特别地,健身活动、户外跑步等运动方式对于实时动态心率监测有着广泛的应用需求。 With the application and development of mobile Internet technology in the medical and health field, a large number of different forms of smart watches, smart bracelets, and smart wristbands have appeared on the market, capable of measuring heart rate, blood pressure, blood oxygen concentration, respiratory rate, etc. Parameters of wearable mobile health products. Heart rate is defined as the number of heartbeats per minute of the human body, and is an important medical routine physiological parameter for evaluating the health status of a person. Heart rate monitoring is of great significance for disease risk early warning, disease diagnosis, and annual routine physical examination. In particular, sports such as fitness activities and outdoor running have extensive application requirements for real-time dynamic heart rate monitoring.

目前,绝大多数移动健康产品采用的心率监测技术原理是光电透射测量法。在产品的硬件设计上,与人体皮肤接触的传感器会发出一束光打在皮肤上,同时测量经皮肤反射或透射的光。因为血液对特定波长的光有吸收作用,心脏泵血的过程直接影响传感器测得的光信号强度改变,硬件按照设定的采样率记录信号强度变化采集原始数据,即脉搏波信号。数据分析软件单元运行心率监测算法处理脉搏波信号,输出心率值。心率监测算法是心率测量产品的关键核心技术,决定了心率测量值的准确性和可靠性。实际应用中,心率监测包括静态心率监测和实时动态心率监测,后者具有更广的应用空间,同时也对现有技术提出了很大的挑战。 At present, the principle of heart rate monitoring technology adopted by most mobile health products is photoelectric transmission measurement. In the hardware design of the product, the sensor in contact with the human skin will emit a beam of light to hit the skin, and measure the light reflected or transmitted through the skin at the same time. Because the blood absorbs light of a specific wavelength, the process of heart pumping directly affects the change of the light signal intensity measured by the sensor. The hardware records the signal intensity change according to the set sampling rate and collects the original data, that is, the pulse wave signal. The data analysis software unit runs the heart rate monitoring algorithm to process the pulse wave signal, and outputs the heart rate value. Heart rate monitoring algorithm is the key core technology of heart rate measurement products, which determines the accuracy and reliability of heart rate measurement values. In practical applications, heart rate monitoring includes static heart rate monitoring and real-time dynamic heart rate monitoring. The latter has a wider application space, but also poses a great challenge to the existing technology.

通过实际测试发现,目前绝大多数移动健康产品的动态心率监测算法的技术现状是传统的信号时域波形分析或信号频域分析结合加速度计读数辅助判断,在运动状态下的人体心率实时监测这一应用场景中暴露了以下缺点。第一,由于噪声干扰导致信号时域波形上的特征点不太明显,会造成算法无法获取完整输入信息;第二,信号时域波形匹配的准则设置过多,算法参数具体数值设定存在困难;第三,对于嵌入式模块而言,波形匹配算法的计算复杂度较大;第四,传统的信号频域分析方法在提高频谱精度时会增大计算和数据存储开销;第五,加速度计增加了硬件成本,同时也增加计算、存储和能耗方面的资源开销。 Through actual tests, it is found that the current technical status of the dynamic heart rate monitoring algorithm of most mobile health products is traditional signal time-domain waveform analysis or signal frequency domain analysis combined with accelerometer readings to assist judgment, real-time monitoring of human heart rate in the state of exercise. The following disadvantages are exposed in an application scenario. First, due to noise interference, the feature points on the signal time-domain waveform are not obvious, which will cause the algorithm to fail to obtain complete input information; second, there are too many criteria for signal time-domain waveform matching, and it is difficult to set the specific value of the algorithm parameters ; Third, for embedded modules, the computational complexity of the waveform matching algorithm is large; Fourth, the traditional signal frequency domain analysis method will increase the calculation and data storage overhead when improving the frequency spectrum accuracy; Fifth, the accelerometer This increases hardware costs and also increases resource overhead in computing, storage, and energy consumption.

因此,现有技术还有待于改进和发展。 Therefore, the prior art still needs to be improved and developed.

发明内容 Contents of the invention

鉴于现有技术的不足,本发明目的在于提供一种实时动态心率监测装置及监测方法,旨在解决现有技术中动态心率监测装置在监测心率时,算法计算复杂度高,硬件成本高,计算、存储和能源方面的资源开销大的缺陷。 In view of the deficiencies in the prior art, the purpose of the present invention is to provide a real-time dynamic heart rate monitoring device and a monitoring method, aiming at solving the problem of high computational complexity of the algorithm, high hardware cost, and computational complexity of the dynamic heart rate monitoring device in the prior art when monitoring heart rate , storage, and energy resource overhead are large drawbacks.

本发明的技术方案如下: Technical scheme of the present invention is as follows:

一种实时动态心率监测装置,其中,装置包括基线漂移消除模块,带通滤波模块,频域分析模块,心率频点选择模块; A real-time dynamic heart rate monitoring device, wherein the device includes a baseline drift elimination module, a band-pass filter module, a frequency domain analysis module, and a heart rate frequency point selection module;

所述基线漂移消除模块用于消除导致脉搏波信号的基准线出现漂移的干扰信号; The baseline drift elimination module is used to eliminate interference signals that cause the baseline of the pulse wave signal to drift;

所述带通滤波模块用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号; The band-pass filter module is used to obtain frequency signal components belonging to the heart rate frequency band while eliminating noise signals outside the heart rate frequency band;

所述频域分析模块用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱; The frequency domain analysis module is used to obtain the signal spectrum of the pulse wave signal in a preset specific frequency band interval;

所述心率频点选择模块用于获取心率对应的频点并输出频点; The heart rate frequency point selection module is used to obtain the frequency point corresponding to the heart rate and output the frequency point;

所述基线漂移消除模块与所述带通滤波模块连接,所述带通滤波模块与所述频域分析模块连接,所述频域分析模块还与所述心率频点选择模块连接。 The baseline drift elimination module is connected to the band-pass filter module, the band-pass filter module is connected to the frequency domain analysis module, and the frequency domain analysis module is also connected to the heart rate frequency point selection module.

所述的实时动态心率监测装置,其中,所述基线漂移消除模块具体包括基线漂移趋势项信号提取单元和信号线性叠加单元; The real-time dynamic heart rate monitoring device, wherein the baseline drift elimination module specifically includes a baseline drift trend item signal extraction unit and a signal linear superposition unit;

所述基线漂移趋势项信号提取单元用于获取与原脉搏波信号等长的基线漂移趋势项信号; The baseline drift trend item signal extraction unit is used to obtain a baseline drift trend item signal of the same length as the original pulse wave signal;

所述信号线性叠加单元用于将原脉搏波信号减掉基线漂移趋势项信号,得到去掉基线漂移的脉搏波信号; The signal linear superposition unit is used to subtract the baseline drift trend item signal from the original pulse wave signal to obtain the pulse wave signal with the baseline drift removed;

所述基线漂移趋势项信号提取单元与所述信号线性叠加单元连接。 The baseline drift trend item signal extraction unit is connected to the signal linear superposition unit.

所述的实时动态心率监测装置,其中,所述带通滤波模块具体包括滤波参数设置单元和滤波单元, The real-time dynamic heart rate monitoring device, wherein the bandpass filter module specifically includes a filter parameter setting unit and a filter unit,

所述滤波参数设置单元用于设置通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数; The filter parameter setting unit is used to set the lower limit and upper limit frequency of the passband, the lower limit and upper limit frequency of the stopband, the attenuation coefficient in the passband, and the attenuation coefficient in the stopband;

所述滤波单元用于获取滤波模块阶数和系数后,将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波,获取去除噪声后的脉搏波信号; The filter unit is used to filter the pulse wave signal after the baseline drift is eliminated by the baseline drift elimination module after obtaining the order and coefficient of the filter module, and obtain the pulse wave signal after the noise is removed;

所述滤波参数设置单元和滤波单元连接。 The filtering parameter setting unit is connected to the filtering unit.

所述的实时动态心率监测装置,其中,所述频域分析模块具体包括特定频段区间设置单元、频点功率计算单元和信号频谱获取单元; The real-time dynamic heart rate monitoring device, wherein the frequency domain analysis module specifically includes a specific frequency band interval setting unit, a frequency point power calculation unit and a signal spectrum acquisition unit;

所述特定频段区间设置单元用于设置特定频段区间的起始频点、结束频点、频点细分数目; The specific frequency band interval setting unit is used to set the start frequency point, the end frequency point, and the number of frequency point subdivisions of the specific frequency band interval;

所述频点功率计算单元用于对于特定频段区间的每个频点,分别用第一变换多项式计算频域信号的实部,用第二变换多项式计算频域信号的虚部,根据频域信号的实部和虚部计算得到脉搏波信号在频点的功率; The frequency point power calculation unit is used to calculate the real part of the frequency domain signal with the first transformation polynomial and the imaginary part of the frequency domain signal with the second transformation polynomial for each frequency point of the specific frequency band interval, according to the frequency domain signal Calculate the real and imaginary parts of the pulse wave signal at the frequency;

所述信号频谱获取单元用于获取特定频段区间内的脉搏波信号在每个频点的功率,叠加后生成脉搏波信号在整个关注频段区间的信号频谱; The signal spectrum acquisition unit is used to acquire the power of the pulse wave signal at each frequency point in a specific frequency band interval, and generate the signal spectrum of the pulse wave signal in the entire frequency band interval of interest after superposition;

所述频点功率计算单元分别与所述特定频段区间设置单元、所述信号频谱获取单元连接。 The frequency point power calculation unit is respectively connected with the specific frequency band interval setting unit and the signal spectrum acquisition unit.

所述的实时动态心率监测装置,其中,所述心率频点选择模块具体包括心率频点选择单元和心率频点动态跟踪单元, The real-time dynamic heart rate monitoring device, wherein the heart rate frequency point selection module specifically includes a heart rate frequency point selection unit and a heart rate frequency point dynamic tracking unit,

所述心率频点选择单元用于在初始状态下获取脉搏波信号频谱中谱峰所在的频点位置作为心率频点; The heart rate frequency point selection unit is used to obtain the frequency point where the peak in the pulse wave signal spectrum is located as the heart rate frequency point in the initial state;

所述心率频点动态跟踪单元用于在预先设置的固定跟踪周期内,以上一跟踪周期输出的心率频点为中心,在中心左右一个特定宽度的频谱范围内进行高精度的频域分析,选择谱峰所在的频点位置作为心率频点输出; The heart rate frequency point dynamic tracking unit is used to perform high-precision frequency domain analysis in a spectrum range of a specific width around the center with the heart rate frequency point output from the previous tracking cycle as the center within the preset fixed tracking cycle, and select The frequency point where the spectrum peak is located is output as the heart rate frequency point;

所述心率频点选择单元与所述心率频点动态跟踪单元连接。 The heart rate frequency point selection unit is connected with the heart rate frequency point dynamic tracking unit.

所述的实时动态心率监测装置,其中,所述第一变换多项式和所述第二变换多项式为正交多项式。 The real-time ambulatory heart rate monitoring device, wherein, the first transformation polynomial and the second transformation polynomial are orthogonal polynomials.

所述的实时动态心率监测装置,其中,所述第一变换多项式为勒让德多项式,雅可比多项式,拉盖尔多项式,切比雪夫多项式,埃尔米特多项式中的一种; The real-time dynamic heart rate monitoring device, wherein the first transformation polynomial is one of Legendre polynomials, Jacobi polynomials, Laguerre polynomials, Chebyshev polynomials, and Hermite polynomials;

所述第二变换多项式为正交多项式为勒让德多项式,雅可比多项式,拉盖尔多项式,切比雪夫多项式,埃尔米特多项式中的一种。 The second transformation polynomial is an orthogonal polynomial and is one of Legendre polynomials, Jacobi polynomials, Laguerre polynomials, Chebyshev polynomials, and Hermitian polynomials.

一种基于所述的实时动态心率监测装置的监测方法,其中,方法包括步骤: A monitoring method based on the real-time dynamic heart rate monitoring device, wherein the method comprises the steps of:

A、获取脉搏波信号,通过基线漂移消除模块消除基线漂移干扰信号; A. Obtain the pulse wave signal, and eliminate the baseline drift interference signal through the baseline drift elimination module;

B、带通滤波模块对基线漂移消除模块输出的信号进行滤波,保留属于心率频段的信号分量; B. The band-pass filter module filters the signal output by the baseline drift elimination module, and retains the signal components belonging to the heart rate frequency band;

C、频域分析模块获取脉搏波信号在预先设置的特定频段区间内的信号频谱; C. The frequency domain analysis module acquires the signal spectrum of the pulse wave signal in the preset specific frequency band interval;

D、心率频点选择模块根据信号频率中获取心率对应的频点并输出。 D. The heart rate frequency point selection module obtains and outputs the frequency point corresponding to the heart rate according to the signal frequency.

所述的实时动态心率监测方法,其中,所述步骤A具体包括: The real-time dynamic heart rate monitoring method, wherein, the step A specifically includes:

A1、通过信号滤波法或曲线拟合法获取与原脉搏波信号等长的基线漂移趋势项信号; A1. Obtain the baseline drift trend item signal with the same length as the original pulse wave signal by signal filtering method or curve fitting method;

A2、将原脉搏波信号,基线漂移趋势项信号使用行向量或者列向量存储,然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信,得到去掉基线漂移的脉搏波信号。 A2. Store the original pulse wave signal and the baseline drift trend item signal in a row vector or column vector, and then subtract the baseline drift trend item information from the original pulse wave signal according to the matrix addition rule to obtain the pulse wave signal with the baseline drift removed.

所述的实时动态心率监测方法,其中,所述步骤B具体包括: The real-time dynamic heart rate monitoring method, wherein, the step B specifically includes:

B1、获取预先设置的通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数; B1. Obtain the preset lower limit and upper limit frequency of the passband, the lower limit and upper limit frequency of the stopband, the attenuation coefficient in the passband, and the attenuation coefficient in the stopband;

B2、获取滤波模块阶数和系数后,将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波,获取去除噪声后的脉搏波信号。 B2. After obtaining the order and coefficient of the filtering module, filter the pulse wave signal after the baseline drift is eliminated by the baseline drift elimination module, and obtain the pulse wave signal after the noise is removed.

本发明提供了一种实时动态心率监测装置及监测方法,本发明可排除人体运动、肌肉活动对心率分析的影响,减少了硬件成本,可以灵活地选择关注频段区间,仅在关注频段区间内提高信号的频域精度,避免带来过大的计算和存储开销,选择信号频谱分析的技术路线,避免了时域信号波形特征点难找,算法参数具体值难确定的问题。 The present invention provides a real-time dynamic heart rate monitoring device and a monitoring method. The present invention can eliminate the influence of human body movement and muscle activity on heart rate analysis, reduce hardware costs, and can flexibly select the frequency range of interest, and only improve the heart rate within the frequency range of interest. The frequency domain accuracy of the signal avoids excessive calculation and storage overhead, and the technical route of signal spectrum analysis is selected to avoid the problems of finding the characteristic points of the time domain signal waveform and determining the specific value of the algorithm parameters.

附图说明 Description of drawings

图1为本发明的一种实时动态心率监测装置的较佳实施例的功能原理框图。 Fig. 1 is a functional principle block diagram of a preferred embodiment of a real-time dynamic heart rate monitoring device of the present invention.

图2为本发明的一种实时动态心率监测装置的监测方法的较佳实施例的流程图。 Fig. 2 is a flowchart of a preferred embodiment of a monitoring method of a real-time dynamic heart rate monitoring device of the present invention.

具体实施方式 detailed description

为使本发明的目的、技术方案及效果更加清楚、明确,以下对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。 In order to make the object, technical solution and effect of the present invention more clear and definite, the present invention will be further described in detail below. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

本发明还提供了一种实时动态心率监测装置的较佳实施例的功能原理图,如图1所示,其中,装置包括基线漂移消除模块100,带通滤波模块200,频域分析模块300,心率频点选择模块400; The present invention also provides a functional schematic diagram of a preferred embodiment of a real-time dynamic heart rate monitoring device, as shown in Figure 1, wherein the device includes a baseline drift elimination module 100, a bandpass filter module 200, a frequency domain analysis module 300, Heart rate frequency selection module 400;

所述基线漂移消除模块100用于消除导致脉搏波信号的基准线出现漂移的干扰信号;所述带通滤波模块200用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号;所述频域分析模块300用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱;所述心率频点选择模块400用于获取心率对应的频点并输出频点;所述基线漂移消除模块100与所述带通滤波模块200连接,所述带通滤波模块200与所述频域分析模块300连接,所述频域分析模块300还与所述心率频点选择模块400连接。 The baseline drift elimination module 100 is used to eliminate the interference signal that causes the baseline of the pulse wave signal to drift; the bandpass filter module 200 is used to obtain frequency signal components belonging to the heart rate frequency band while eliminating noise signals outside the heart rate frequency band; The frequency domain analysis module 300 is used to obtain the signal spectrum of the pulse wave signal in a preset specific frequency range; the heart rate frequency point selection module 400 is used to obtain the frequency point corresponding to the heart rate and output the frequency point; the baseline drift elimination The module 100 is connected to the bandpass filter module 200 , the bandpass filter module 200 is connected to the frequency domain analysis module 300 , and the frequency domain analysis module 300 is also connected to the heart rate frequency point selection module 400 .

具体实施时,所述基线漂移消除模块具体包括基线漂移趋势项信号提取单元和信号线性叠加单元; During specific implementation, the baseline drift elimination module specifically includes a baseline drift trend item signal extraction unit and a signal linear superposition unit;

所述基线漂移趋势项信号提取单元用于获取与原脉搏波信号等长的基线漂移趋势项信号; The baseline drift trend item signal extraction unit is used to obtain a baseline drift trend item signal of the same length as the original pulse wave signal;

所述信号线性叠加单元用于将原脉搏波信号减掉基线漂移趋势项信号,得到去掉基线漂移的脉搏波信号; The signal linear superposition unit is used to subtract the baseline drift trend item signal from the original pulse wave signal to obtain the pulse wave signal with the baseline drift removed;

所述基线漂移趋势项信号提取单元与所述信号线性叠加单元连接。 The baseline drift trend item signal extraction unit is connected to the signal linear superposition unit.

具体实施时,因为呼吸、肢体活动或运动会导致脉搏波信号的基准线呈现上下漂移的情况,本质上这种基线漂移是一种低频信号,会对关键信息的判断带来干扰。消除基线漂移模块的主要作用就是降低呼吸、运动对脉搏波信号带来的基线漂移干扰,避免对加速度计辅助信息的依赖。 During the specific implementation, the baseline of the pulse wave signal will drift up and down because of breathing, body movement or movement. Essentially, this baseline drift is a low-frequency signal that will interfere with the judgment of key information. The main function of the elimination of baseline drift module is to reduce the baseline drift interference caused by respiration and exercise on the pulse wave signal, and to avoid dependence on the auxiliary information of the accelerometer.

消除基线漂移模块由两个功能单元组成:基线漂移趋势项信号提取单元,信号线性叠加单元。 The module for eliminating baseline drift consists of two functional units: a baseline drift trend item signal extraction unit, and a signal linear superposition unit.

基线漂移趋势项信号提取单元,具体的实施方式有两类,第一类是信号滤波方法,第二类是曲线拟合方法。信号滤波方法的具体实施方式为中值滤波或均值滤波,使用滑动窗口遍历脉搏波信号,计算窗口内的全部信号值的中值或均值,最终得到与原脉搏波信号等长的基线漂移趋势项信号。曲线拟合方法的具体实施方式,将基线漂移信号视作可以表示为一个高次多项式的时间函数,用脉搏波信号作为输入数据样本,采用非线性拟合方法得到高次多项式的各项系数,最终通过多项式函数计算得到与原脉搏波信号等长的基线漂移趋势项信号。 The baseline drift trend item signal extraction unit has two specific implementation methods, the first type is a signal filtering method, and the second type is a curve fitting method. The specific implementation of the signal filtering method is median filtering or mean filtering, using a sliding window to traverse the pulse wave signal, calculating the median or mean value of all signal values in the window, and finally obtaining a baseline drift trend item with the same length as the original pulse wave signal Signal. The specific implementation of the curve fitting method considers the baseline drift signal as a time function that can be expressed as a high-order polynomial, uses the pulse wave signal as an input data sample, and adopts a nonlinear fitting method to obtain the coefficients of the high-order polynomial, Finally, the baseline drift trend item signal with the same length as the original pulse wave signal is obtained through polynomial function calculation.

信号线性叠加单元,具体实施方式就是将原脉搏波信号,基线漂移趋势项信号使用行向量或者列向量存储,然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信号,最终得到去掉基线漂移的脉搏波信号。 The signal linear superposition unit, the specific implementation method is to store the original pulse wave signal and the baseline drift trend item signal in a row vector or column vector, and then subtract the baseline drift trend item signal from the original pulse wave signal according to the matrix addition rule, and finally get the baseline drift signal Drifting pulse wave signal.

进一步地,所述带通滤波模块具体包括滤波参数设置单元和滤波单元, Further, the bandpass filtering module specifically includes a filtering parameter setting unit and a filtering unit,

所述滤波参数设置单元用于设置通带下限和上限频率、阻带下限和上限频率和上限频率、通带内衰减系数、阻带内衰减系数; The filter parameter setting unit is used to set the lower limit and upper limit frequency of the passband, the lower limit and upper limit frequency and the upper limit frequency of the stopband, the attenuation coefficient in the passband, and the attenuation coefficient in the stopband;

所述滤波单元用于获取滤波模块阶数和系数后,将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波,获取去除噪声后的脉搏波信号; The filter unit is used to filter the pulse wave signal after the baseline drift is eliminated by the baseline drift elimination module after obtaining the order and coefficient of the filter module, and obtain the pulse wave signal after the noise is removed;

所述滤波参数设置单元和滤波单元连接。 The filtering parameter setting unit is connected to the filtering unit.

具体实施时,大量医学实践表明,普遍情况下人体心率的范围是40次/分钟至220次/分钟,这表明心率信号的频率范围是0.6Hz至3.7Hz。因此,可以将0Hz至0.6Hz,以及3.7Hz以上的信号频段视作噪声频段,这些噪声包括传感器电路导致的毛刺噪声、人体活动引入的噪声等等。带通滤波器模块的主要作用就是仅保留属于心率频段的信号分量,同时消除带外噪声对于发现心率信号的干扰。 During specific implementation, a large number of medical practices have shown that the range of the human heart rate is generally 40 beats/minute to 220 beats/minute, which indicates that the frequency range of the heart rate signal is 0.6 Hz to 3.7 Hz. Therefore, the signal frequency bands from 0 Hz to 0.6 Hz and above 3.7 Hz can be regarded as noise frequency bands. These noises include glitch noise caused by sensor circuits, noise introduced by human activities, and so on. The main function of the band-pass filter module is to retain only the signal components belonging to the heart rate frequency band, and at the same time eliminate the interference of out-of-band noise to the detection of the heart rate signal.

带通滤波器模块,具体的实施方式可以采用巴特沃斯滤波器,通过设置通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数,参照巴特沃斯滤波器设计流程计算得到滤波器阶数和系数。将去掉基线漂移的脉搏波信号输入带通滤波器模块,最终得到去除噪声的脉搏波信号。 Band-pass filter module, the specific embodiment can adopt Butterworth filter, by setting passband lower limit and upper limit frequency, stopband lower limit and upper limit frequency, passband attenuation coefficient, stopband attenuation coefficient, refer to Butterworth The filter design process calculates the filter order and coefficients. The pulse wave signal from which the baseline drift has been removed is input into the band-pass filter module, and finally a pulse wave signal from which noise has been removed is obtained.

在具体实施时,不排除存在极端情况下(疾病、剧烈运动)的心率信号,为保证全面地监测到心率变化,可以适当增大心率频段的动态范围,例如,设置成0.5Hz至4Hz;带通滤波器的实施方式除巴特沃斯滤波器外,还可以选择切比雪夫滤波器等。 In the specific implementation, it is not ruled out that there are heart rate signals in extreme cases (disease, strenuous exercise). In order to ensure that heart rate changes are fully monitored, the dynamic range of the heart rate frequency band can be appropriately increased, for example, set to 0.5Hz to 4Hz; In addition to the Butterworth filter, Chebyshev filter and the like can also be selected for the implementation of the pass filter.

进一步的实施例中,所述频域分析模块具体包括特定频段区间设置单元、频点功率计算单元和信号频谱获取单元; In a further embodiment, the frequency domain analysis module specifically includes a specific frequency band interval setting unit, a frequency point power calculation unit, and a signal spectrum acquisition unit;

所述特定频段区间设置单元用于设置特定频段区间的起始频点、结束频点、频点细分数目; The specific frequency band interval setting unit is used to set the start frequency point, the end frequency point, and the number of frequency point subdivisions of the specific frequency band interval;

所述频点功率计算单元用于对于特定频段区间的每个频点,分别用第一变换多项式计算频域信号的实部,用第二变换多项式计算频域信号的虚部,根据频域信号的实部和虚部计算得到脉搏波信号在频点的功率; The frequency point power calculation unit is used to calculate the real part of the frequency domain signal with the first transformation polynomial and the imaginary part of the frequency domain signal with the second transformation polynomial for each frequency point of the specific frequency band interval, according to the frequency domain signal Calculate the real and imaginary parts of the pulse wave signal at the frequency;

所述信号频谱获取单元用于获取特定频段区间内的脉搏波信号在每个频点的功率,叠加后生成脉搏波信号在整个关注频段区间的信号频谱; The signal spectrum acquisition unit is used to acquire the power of the pulse wave signal at each frequency point in a specific frequency band interval, and generate the signal spectrum of the pulse wave signal in the entire frequency band interval of interest after superposition;

所述频点功率计算单元分别与所述特定频段区间设置单元、所述信号频谱获取单元连接。 The frequency point power calculation unit is respectively connected with the specific frequency band interval setting unit and the signal spectrum acquisition unit.

具体实施时,去除噪声的脉搏波信号是时域信号,为了更精准的提取心率信息,需要采用信号变换方法获得脉搏波信号对应的频域信号,并且保证信号频谱有足够高的频域精度。经典的时频信号变换方法是快速傅里叶变换即FFT,它的缺点是为提高信号频域精度需要付出嵌入式设备难以承担的计算和存储开销。高精度频域分析模块的主要作用是,可以灵活地选择关注频段区间,仅在关注频段区间内提高信号频域精度,避免带来过大的计算和存储开销。 In specific implementation, the pulse wave signal from which the noise is removed is a time domain signal. In order to extract heart rate information more accurately, it is necessary to use a signal transformation method to obtain the frequency domain signal corresponding to the pulse wave signal, and to ensure that the signal spectrum has a sufficiently high frequency domain accuracy. The classic time-frequency signal transformation method is the fast Fourier transform or FFT. Its disadvantage is that in order to improve the frequency domain accuracy of the signal, it needs to pay the calculation and storage overhead that the embedded device cannot bear. The main function of the high-precision frequency domain analysis module is to flexibly select the frequency band interval of interest, and improve the frequency domain accuracy of the signal only in the frequency band interval of interest, avoiding excessive calculation and storage overhead.

高精度频域分析模块,具体的实施方式为,设置关注频段区间的起始频点f1、结束频点f2和频点细分数目Nf。 The specific implementation of the high-precision frequency domain analysis module is to set the starting frequency point f1, the ending frequency point f2, and the number of subdivided frequency points Nf of the frequency band interval of interest.

用q(n),n=0,1,…,N-1表示数据长度为N的脉搏波信号序列。对于关注频段区间内的每个频点f,使用下面的变换多项式计算频域信号实部, Use q(n), n=0, 1, . . . , N-1 to represent a pulse wave signal sequence with a data length of N. For each frequency point f in the frequency band interval of interest, the real part of the frequency domain signal is calculated using the following transformation polynomial,

Qr=q(0)+q(1)*cos(f)+q(2)*T(2)+…+q(N-1)*T(N-1) Qr=q(0)+q(1)*cos(f)+q(2)*T(2)+…+q(N-1)*T(N-1)

和使用下面的变换多项式计算频域信号虚部, and compute the imaginary part of the frequency-domain signal using the following transformation polynomial,

Qi=-q(1)*sin(f)–q(2)*sin(f)*U(2)-…-q(N-1)*sin(f)*U(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)为第二变换多项式,第一变换多项式和所述第二变换多项式为正交多项式。所述第一变换多项式为勒让德多项式,雅可比多项式,拉盖尔多项式,切比雪夫多项式,埃尔米特多项式中的一种;所述第二变换多项式为正交多项式为勒让德多项式,雅可比多项式,拉盖尔多项式,切比雪夫多项式,埃尔米特多项式中的一种。 T(n) is a first transformation polynomial, U(n) is a second transformation polynomial, and the first transformation polynomial and the second transformation polynomial are orthogonal polynomials. The first transformation polynomial is one of Legendre polynomials, Jacobi polynomials, Laguerre polynomials, Chebyshev polynomials, and Hermitian polynomials; the second transformation polynomial is an orthogonal polynomial that is Legendre One of polynomials, Jacobi polynomials, Laguerre polynomials, Chebyshev polynomials, and Hermite polynomials.

根据Qr和Qi,可以计算得到脉搏波信号在频点f的功率。 According to Qr and Qi, the power of the pulse wave signal at the frequency point f can be calculated.

当关注频段区间内的每个频点都经过以上计算后,则可以得到脉搏波信号在整个关注频段区间的信号频谱。高频域精度的量化指标与频点细分数目Nf有关,频点细分数目越大频域精度越高。 After each frequency point in the frequency band of interest is calculated above, the signal spectrum of the pulse wave signal in the entire frequency band of interest can be obtained. The quantitative index of high-frequency domain accuracy is related to the number of frequency subdivisions Nf, and the greater the number of frequency subdivisions, the higher the frequency domain accuracy.

进一步地,所述心率频点选择模块具体包括心率频点选择单元和心率频点动态跟踪单元, Further, the heart rate frequency point selection module specifically includes a heart rate frequency point selection unit and a heart rate frequency point dynamic tracking unit,

所述心率频点选择单元用于在初始状态下获取脉搏波信号频谱中谱峰所在的频点位置作为心率频点; The heart rate frequency point selection unit is used to obtain the frequency point where the peak in the pulse wave signal spectrum is located as the heart rate frequency point in the initial state;

所述心率频点动态跟踪单元用于在预先设置的固定跟踪周期内,以上一跟踪周期输出的心率频点为中心,在中心左右一个特定宽度的频谱范围内进行高精度的频域分析,选择谱峰所在的频点位置作为心率频点输出; The heart rate frequency point dynamic tracking unit is used to perform high-precision frequency domain analysis in a spectrum range of a specific width around the center with the heart rate frequency point output from the previous tracking cycle as the center within the preset fixed tracking cycle, and select The frequency point where the spectrum peak is located is output as the heart rate frequency point;

所述心率频点选择单元与所述心率频点动态跟踪单元连接。 The heart rate frequency point selection unit is connected with the heart rate frequency point dynamic tracking unit.

具体实施时,得到脉搏波信号在关注频段区间的信号频谱后,并不能很直接地计算心率值,而是需要通过分析判断发现代表心率的频点,才能计算心率值,这是心率频点选择模块的作用。 In the specific implementation, after obtaining the signal spectrum of the pulse wave signal in the attention frequency range, the heart rate value cannot be calculated directly, but the heart rate value needs to be found through analysis and judgment to find the frequency point representing the heart rate, which is the selection of the heart rate frequency point The role of the module.

心率频点选择模块,具体实施方式,分为两个阶段。第一阶段是,初始状态精确心率频点选择;第二阶段是心率频点动态跟踪。 The heart rate frequency point selection module is specifically implemented in two stages. The first stage is the precise selection of heart rate frequency points in the initial state; the second stage is dynamic tracking of heart rate frequency points.

初始状态精确心率频点选择,具体实施方式是,保证人体处于安静状态下,通过消除基线漂移模块、带通滤波器模块和高精度频域分析模块得到脉搏波信号频谱,其中谱峰所在频点位置就是代表心率的频点,可以计算出心率值。 Precise heart rate frequency point selection in the initial state. The specific implementation method is to ensure that the human body is in a quiet state, and obtain the pulse wave signal spectrum by eliminating the baseline drift module, band-pass filter module and high-precision frequency domain analysis module. The position is the frequency point representing the heart rate, and the heart rate value can be calculated.

心率频点动态跟踪,具体实施方式是,以固定的跟踪周期进行循环,在每个跟踪周期内,通过消除基线漂移模块和带通滤波器模块得到去除噪声的脉搏波信号,以上一跟踪周期输出的心率频点为中心,在中心左右一个狭小宽度的频谱范围内进行高精度的频域分析,选择谱峰所在的频点位置作为心率频点输出。 Heart rate frequency point dynamic tracking, the specific implementation method is to cycle with a fixed tracking period, and in each tracking period, the pulse wave signal with noise removed is obtained by eliminating the baseline drift module and the band-pass filter module, and output in the previous tracking period The heart rate frequency point is the center, and high-precision frequency domain analysis is performed in a narrow width spectrum range around the center, and the frequency point where the spectrum peak is located is selected as the heart rate frequency point output.

本发明还提供了一种基于所述的实时动态心率监测装置的监测方法的较佳实施例,如图2所示,其中,方法包括: The present invention also provides a preferred embodiment of a monitoring method based on the real-time dynamic heart rate monitoring device, as shown in Figure 2, wherein the method includes:

步骤S100、获取脉搏波信号,通过基线漂移消除模块消除基线漂移干扰信号; Step S100, acquire the pulse wave signal, and eliminate the baseline drift interference signal through the baseline drift elimination module;

步骤S200、带通滤波模块对基线漂移消除模块输出的信号进行滤波,保留属于心率频段的信号分量; Step S200, the bandpass filter module filters the signal output by the baseline drift elimination module, and retains the signal components belonging to the heart rate frequency band;

步骤S300、频域分析模块获取脉搏波信号在预先设置的特定频段区间内的信号频谱; Step S300, the frequency domain analysis module acquires the signal spectrum of the pulse wave signal within a preset specific frequency band interval;

步骤S400、心率频点选择模块根据信号频率中获取心率对应的频点并输出。 Step S400, the heart rate frequency point selection module obtains and outputs the frequency point corresponding to the heart rate according to the signal frequency.

具体实施时,通过多种滤波器技术去除信号采集过程中混入的高频噪声,去除由肌肉抖动和呼吸引起的信号基线漂移,采用一种比FFT计算复杂度小的高精度信号频域分析方法,从信号频谱中筛选定位代表心率的频点,对心率进行实时动态地监测。减少了加速度计导致的硬件成本和资源上的额外开销。第二,避免了时域信号波形特征点难找,算法参数具体值难确定的问题。第三,算法的计算复杂度是普通嵌入式模块可以承担的。具体的监测方法如上监测装置的具体实施例所述。 In the specific implementation, the high-frequency noise mixed in the signal acquisition process is removed through a variety of filter techniques, the signal baseline drift caused by muscle shaking and breathing is removed, and a high-precision signal frequency domain analysis method with less computational complexity than FFT is adopted. , select and locate the frequency point representing the heart rate from the signal spectrum, and monitor the heart rate dynamically in real time. The hardware cost and resource overhead caused by the accelerometer is reduced. Second, it avoids the problems of finding the characteristic points of the time-domain signal waveform and determining the specific values of the algorithm parameters. Third, the computational complexity of the algorithm can be borne by ordinary embedded modules. The specific monitoring method is as described in the specific embodiment of the monitoring device above.

进一步地,所述步骤S100具体包括: Further, the step S100 specifically includes:

步骤S101、通过信号滤波法或曲线拟合法获取与原脉搏波信号等长的基线漂移趋势项信号; Step S101. Obtain a baseline drift trend item signal with the same length as the original pulse wave signal by signal filtering method or curve fitting method;

步骤S102、将原脉搏波信号,基线漂移趋势项信号使用行向量或者列向量存储,然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信,得到去掉基线漂移的脉搏波信号。 Step S102: Store the original pulse wave signal and the baseline drift trend item signal in a row vector or column vector, and then subtract the baseline drift trend item information from the original pulse wave signal according to the matrix addition rule to obtain the pulse wave signal with the baseline drift removed.

具体实施时,其中所述步骤S101中信号滤波方法的具体实施方式为中值滤波或均值滤波,使用滑动窗口遍历脉搏波信号,计算窗口内的全部信号值的中值或均值,最终得到与原脉搏波信号等长的基线漂移趋势项信号。曲线拟合方法的具体实施方式,将基线漂移信号视作可以表示为一个高次多项式的时间函数,用脉搏波信号作为输入数据样本,采用非线性拟合方法得到高次多项式的各项系数,最终通过多项式函数计算得到与原脉搏波信号等长的基线漂移趋势项信号。 During specific implementation, the specific implementation of the signal filtering method in the step S101 is median filtering or mean filtering, using a sliding window to traverse the pulse wave signal, calculating the median or mean value of all signal values in the window, and finally obtaining the same value as the original A baseline drift trend term signal equal in length to the pulse wave signal. The specific implementation of the curve fitting method considers the baseline drift signal as a time function that can be expressed as a high-order polynomial, uses the pulse wave signal as an input data sample, and adopts a nonlinear fitting method to obtain the coefficients of the high-order polynomial, Finally, the baseline drift trend item signal with the same length as the original pulse wave signal is obtained through polynomial function calculation.

进一步的实施例中,所述步骤S200具体包括: In a further embodiment, the step S200 specifically includes:

步骤S201、获取预先设置的通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数; Step S201, obtaining the preset passband lower limit and upper limit frequency, stopband lower limit and upper limit frequency, passband attenuation coefficient, and stopband attenuation coefficient;

步骤S202、获取滤波模块阶数和系数后,将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波,获取去除噪声后的脉搏波信号。 Step S202, after obtaining the order and coefficient of the filtering module, filter the pulse wave signal after the baseline drift is eliminated by the baseline drift elimination module, and obtain the pulse wave signal after the noise is removed.

具体实施时,采用巴特沃斯滤波器,通过设置通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数,参照巴特沃斯滤波器设计流程计算得到滤波器阶数和系数。将去掉基线漂移的脉搏波信号输入带通滤波器模块,最终得到去除噪声的脉搏波信号。 In the specific implementation, the Butterworth filter is used, by setting the lower limit and upper limit frequency of the pass band, the lower limit and upper limit frequency of the stop band, the attenuation coefficient in the pass band, and the attenuation coefficient in the stop band, and referring to the design process of the Butterworth filter to calculate the filter order and coefficients. The pulse wave signal from which the baseline drift has been removed is input into the band-pass filter module, and finally a pulse wave signal from which noise has been removed is obtained.

在具体实施时,不排除存在极端情况下(疾病、剧烈运动)的心率信号,为保证全面地监测到心率变化,可以适当增大心率频段的动态范围,例如,设置成0.5Hz至4Hz;带通滤波器的实施方式除巴特沃斯滤波器外,还可以选择切比雪夫滤波器等。 In the specific implementation, it is not ruled out that there are heart rate signals in extreme cases (disease, strenuous exercise). In order to ensure that heart rate changes are fully monitored, the dynamic range of the heart rate frequency band can be appropriately increased, for example, set to 0.5Hz to 4Hz; In addition to the Butterworth filter, Chebyshev filter and the like can also be selected for the implementation of the pass filter.

综上所述,本发明提供了一种实时动态心率监测装置及监测方法,装置包括基线漂移消除模块,带通滤波模块,频域分析模块,心率频点选择模块;所述基线漂移消除模块用于消除导致脉搏波信号的基准线出现漂移的干扰信号;所述带通滤波模块用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号;所述频域分析模块用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱;所述心率频点选择模块用于获取心率对应的频点并输出频点。本发明排除人体运动、肌肉活动对心率分析的影响,减少了硬件成本,可以灵活地选择关注频段区间,仅在关注频段区间内提高信号的频域精度,避免带来过大的计算和存储开销。 In summary, the present invention provides a real-time dynamic heart rate monitoring device and monitoring method. The device includes a baseline drift elimination module, a bandpass filter module, a frequency domain analysis module, and a heart rate frequency point selection module; the baseline drift elimination module uses It is used to eliminate the interference signal that causes the baseline of the pulse wave signal to drift; the bandpass filter module is used to obtain frequency signal components belonging to the heart rate frequency band while eliminating noise signals outside the heart rate frequency band; the frequency domain analysis module is used to obtain the pulse wave signal The signal spectrum of the wave signal in a preset specific frequency band interval; the heart rate frequency point selection module is used to obtain the frequency point corresponding to the heart rate and output the frequency point. The present invention eliminates the impact of human body movement and muscle activity on heart rate analysis, reduces hardware costs, can flexibly select the frequency range of interest, and only improves the frequency domain accuracy of the signal within the frequency range of interest, avoiding excessive calculation and storage overhead .

应当理解的是,本发明的应用不限于上述的举例,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,所有这些改进和变换都应属于本发明所附权利要求的保护范围。 It should be understood that the application of the present invention is not limited to the above examples, and those skilled in the art can make improvements or transformations according to the above descriptions, and all these improvements and transformations should belong to the protection scope of the appended claims of the present invention.

Claims (10)

1.一种实时动态心率监测装置,其特征在于,装置包括基线漂移消除模块,带通滤波模块,频域分析模块,心率频点选择模块; 1. A real-time dynamic heart rate monitoring device, characterized in that the device includes a baseline drift elimination module, a bandpass filter module, a frequency domain analysis module, and a heart rate frequency point selection module; 所述基线漂移消除模块用于消除导致脉搏波信号的基准线出现漂移的干扰信号; The baseline drift elimination module is used to eliminate interference signals that cause the baseline of the pulse wave signal to drift; 所述带通滤波模块用于获取属于心率频段的频率信号分量同时消除心率频段外的噪声信号; The band-pass filter module is used to obtain frequency signal components belonging to the heart rate frequency band while eliminating noise signals outside the heart rate frequency band; 所述频域分析模块用于获取脉搏波信号在预先设置的特定频段区间内的信号频谱; The frequency domain analysis module is used to obtain the signal spectrum of the pulse wave signal in a preset specific frequency band interval; 所述心率频点选择模块用于获取心率对应的频点并输出频点; The heart rate frequency point selection module is used to obtain the frequency point corresponding to the heart rate and output the frequency point; 所述基线漂移消除模块与所述带通滤波模块连接,所述带通滤波模块与所述频域分析模块连接,所述频域分析模块还与所述心率频点选择模块连接。 The baseline drift elimination module is connected to the band-pass filter module, the band-pass filter module is connected to the frequency domain analysis module, and the frequency domain analysis module is also connected to the heart rate frequency point selection module. 2.根据权利要求1所述的实时动态心率监测装置,其特征在于,所述基线漂移消除模块具体包括基线漂移趋势项信号提取单元和信号线性叠加单元; 2. The real-time dynamic heart rate monitoring device according to claim 1, wherein the baseline drift elimination module specifically includes a baseline drift trend item signal extraction unit and a signal linear superposition unit; 所述基线漂移趋势项信号提取单元用于获取与原脉搏波信号等长的基线漂移趋势项信号; The baseline drift trend item signal extraction unit is used to obtain a baseline drift trend item signal of the same length as the original pulse wave signal; 所述信号线性叠加单元用于将原脉搏波信号减掉基线漂移趋势项信号,得到去掉基线漂移的脉搏波信号; The signal linear superposition unit is used to subtract the baseline drift trend item signal from the original pulse wave signal to obtain the pulse wave signal with the baseline drift removed; 所述基线漂移趋势项信号提取单元与所述信号线性叠加单元连接。 The baseline drift trend item signal extraction unit is connected to the signal linear superposition unit. 3.根据权利要求2所述的实时动态心率监测装置,其特征在于,所述带通滤波模块具体包括滤波参数设置单元和滤波单元, 3. The real-time dynamic heart rate monitoring device according to claim 2, wherein the bandpass filtering module specifically includes a filtering parameter setting unit and a filtering unit, 所述滤波参数设置单元用于设置通带下限和上限频率、阻带下限和上限频率和上限频率、通带内衰减系数、阻带内衰减系数; The filter parameter setting unit is used to set the lower limit and upper limit frequency of the passband, the lower limit and upper limit frequency and the upper limit frequency of the stopband, the attenuation coefficient in the passband, and the attenuation coefficient in the stopband; 所述滤波单元用于获取滤波模块阶数和系数后,将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波,获取去除噪声后的脉搏波信号; The filter unit is used to filter the pulse wave signal after the baseline drift is eliminated by the baseline drift elimination module after obtaining the order and coefficient of the filter module, and obtain the pulse wave signal after the noise is removed; 所述滤波参数设置单元和滤波单元连接。 The filtering parameter setting unit is connected to the filtering unit. 4.根据权利要求3所述的实时动态心率监测装置,其特征在于,所述频域分析模块具体包括特定频段区间设置单元、频点功率计算单元和信号频谱获取单元; 4. The real-time dynamic heart rate monitoring device according to claim 3, wherein the frequency domain analysis module specifically includes a specific frequency band interval setting unit, a frequency point power calculation unit and a signal spectrum acquisition unit; 所述特定频段区间设置单元用于设置特定频段区间的起始频点、结束频点、频点细分数目; The specific frequency band interval setting unit is used to set the start frequency point, the end frequency point, and the number of frequency point subdivisions of the specific frequency band interval; 所述频点功率计算单元用于对于特定频段区间的每个频点,分别用第一变换多项式计算频域信号的实部,用第二变换多项式计算频域信号的虚部,根据频域信号的实部和虚部计算得到脉搏波信号在频点的功率; The frequency point power calculation unit is used to calculate the real part of the frequency domain signal with the first transformation polynomial and the imaginary part of the frequency domain signal with the second transformation polynomial for each frequency point of the specific frequency band interval, according to the frequency domain signal Calculate the real and imaginary parts of the pulse wave signal at the frequency; 所述信号频谱获取单元用于获取特定频段区间内的脉搏波信号在每个频点的功率,叠加后生成脉搏波信号在整个关注频段区间的信号频谱; The signal spectrum acquisition unit is used to acquire the power of the pulse wave signal at each frequency point in a specific frequency band interval, and generate the signal spectrum of the pulse wave signal in the entire frequency band interval of interest after superposition; 所述频点功率计算单元分别与所述特定频段区间设置单元、所述信号频谱获取单元连接。 The frequency point power calculation unit is respectively connected with the specific frequency band interval setting unit and the signal spectrum acquisition unit. 5.根据权利要求4所述的实时动态心率监测装置,其特征在于,所述心率频点选择模块具体包括心率频点选择单元和心率频点动态跟踪单元, 5. The real-time dynamic heart rate monitoring device according to claim 4, wherein the heart rate frequency point selection module specifically includes a heart rate frequency point selection unit and a heart rate frequency point dynamic tracking unit, 所述心率频点选择单元用于在初始状态下获取脉搏波信号频谱中谱峰所在的频点位置作为心率频点; The heart rate frequency point selection unit is used to obtain the frequency point where the peak in the pulse wave signal spectrum is located as the heart rate frequency point in the initial state; 所述心率频点动态跟踪单元用于在预先设置的固定跟踪周期内,以上一跟踪周期输出的心率频点为中心,在中心左右一个特定宽度的频谱范围内进行高精度的频域分析,选择谱峰所在的频点位置作为心率频点输出; The heart rate frequency point dynamic tracking unit is used to perform high-precision frequency domain analysis in a spectrum range of a specific width around the center, centered on the heart rate frequency point output in the previous tracking cycle, within a preset fixed tracking period, and select The frequency point where the spectrum peak is located is output as the heart rate frequency point; 所述心率频点选择单元与所述心率频点动态跟踪单元连接。 The heart rate frequency point selection unit is connected with the heart rate frequency point dynamic tracking unit. 6.根据权利要求5所述的实时动态心率监测装置,其特征在于,所述第一变换多项式和所述第二变换多项式为正交多项式。 6. The real-time ambulatory heart rate monitoring device according to claim 5, characterized in that, the first transformation polynomial and the second transformation polynomial are orthogonal polynomials. 7.根据权利要求6所述的实时动态心率监测装置,其特征在于,所述第一变换多项式为勒让德多项式,雅可比多项式,拉盖尔多项式,切比雪夫多项式,埃尔米特多项式中的一种; 7. The real-time dynamic heart rate monitoring device according to claim 6, wherein the first transformation polynomial is Legendre polynomial, Jacobi polynomial, Laguerre polynomial, Chebyshev polynomial, Hermitian polynomial one of 所述第二变换多项式为正交多项式为勒让德多项式,雅可比多项式,拉盖尔多项式,切比雪夫多项式,埃尔米特多项式中的一种。 The second transformation polynomial is an orthogonal polynomial and is one of Legendre polynomials, Jacobi polynomials, Laguerre polynomials, Chebyshev polynomials, and Hermitian polynomials. 8.一种基于权利要求1所述的实时动态心率监测装置的监测方法,其特征在于,方法包括步骤: 8. A monitoring method based on the real-time dynamic heart rate monitoring device according to claim 1, wherein the method comprises the steps of: A、获取脉搏波信号,通过基线漂移消除模块消除基线漂移干扰信号; A. Obtain the pulse wave signal, and eliminate the baseline drift interference signal through the baseline drift elimination module; B、带通滤波模块对基线漂移消除模块输出的信号进行滤波,保留属于心率频段的信号分量; B. The band-pass filter module filters the signal output by the baseline drift elimination module, and retains the signal components belonging to the heart rate frequency band; C、频域分析模块获取脉搏波信号在预先设置的特定频段区间内的信号频谱; C. The frequency domain analysis module acquires the signal spectrum of the pulse wave signal in the preset specific frequency band interval; D、心率频点选择模块根据信号频率中获取心率对应的频点并输出。 D. The heart rate frequency point selection module obtains and outputs the frequency point corresponding to the heart rate according to the signal frequency. 9.根据权利要求8所述的实时动态心率监测方法,其特征在于,所述步骤A具体包括: 9. The real-time dynamic heart rate monitoring method according to claim 8, wherein said step A specifically comprises: A1、通过信号滤波法或曲线拟合法获取与原脉搏波信号等长的基线漂移趋势项信号; A1. Obtain the baseline drift trend item signal with the same length as the original pulse wave signal by signal filtering method or curve fitting method; A2、将原脉搏波信号,基线漂移趋势项信号使用行向量或者列向量存储,然后按照矩阵加法规则用原脉搏波信号减掉基线漂移趋势项信,得到去掉基线漂移的脉搏波信号。 A2. Store the original pulse wave signal and the baseline drift trend item signal in a row vector or column vector, and then subtract the baseline drift trend item information from the original pulse wave signal according to the matrix addition rule to obtain the pulse wave signal with the baseline drift removed. 10.根据权利要求9所述的实时动态心率监测方法,其特征在于,所述步骤B具体包括: 10. The real-time dynamic heart rate monitoring method according to claim 9, wherein said step B specifically comprises: B1、获取预先设置的通带下限和上限频率、阻带下限和上限频率、通带内衰减系数、阻带内衰减系数; B1. Obtain the preset lower limit and upper limit frequency of the passband, the lower limit and upper limit frequency of the stopband, the attenuation coefficient in the passband, and the attenuation coefficient in the stopband; B2、获取滤波模块阶数和系数后,将经所述基线漂移消除模块消除基线漂移后的脉搏波信号进行滤波,获取去除噪声后的脉搏波信号。 B2. After obtaining the order and coefficient of the filtering module, filter the pulse wave signal after the baseline drift is eliminated by the baseline drift elimination module, and obtain the pulse wave signal after the noise is removed.
CN201610288485.8A 2016-05-04 2016-05-04 A real-time dynamic heart rate monitoring device and monitoring method Expired - Fee Related CN105816165B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610288485.8A CN105816165B (en) 2016-05-04 2016-05-04 A real-time dynamic heart rate monitoring device and monitoring method
CN201910951368.9A CN110876615B (en) 2016-05-04 2016-05-04 Real-time dynamic heart rate monitoring device and monitoring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610288485.8A CN105816165B (en) 2016-05-04 2016-05-04 A real-time dynamic heart rate monitoring device and monitoring method

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN201910951368.9A Division CN110876615B (en) 2016-05-04 2016-05-04 Real-time dynamic heart rate monitoring device and monitoring method

Publications (2)

Publication Number Publication Date
CN105816165A true CN105816165A (en) 2016-08-03
CN105816165B CN105816165B (en) 2019-12-13

Family

ID=56527978

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201610288485.8A Expired - Fee Related CN105816165B (en) 2016-05-04 2016-05-04 A real-time dynamic heart rate monitoring device and monitoring method
CN201910951368.9A Expired - Fee Related CN110876615B (en) 2016-05-04 2016-05-04 Real-time dynamic heart rate monitoring device and monitoring method

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201910951368.9A Expired - Fee Related CN110876615B (en) 2016-05-04 2016-05-04 Real-time dynamic heart rate monitoring device and monitoring method

Country Status (1)

Country Link
CN (2) CN105816165B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107088062A (en) * 2017-05-08 2017-08-25 歌尔科技有限公司 A kind of heart rate module, the electronic installation and heart rate acquisition method for gathering heart rate
CN109924966A (en) * 2017-12-19 2019-06-25 郝振龙 It is a kind of for moving when measure heart rate Intelligent bracelet device
CN110801210A (en) * 2019-11-06 2020-02-18 心核心科技(北京)有限公司 Pulse wave signal filtering method and device, readable medium and electronic equipment
CN113616217A (en) * 2021-10-12 2021-11-09 深圳市倍轻松科技股份有限公司 Method and device for generating baseline drift curve

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115605128B (en) * 2022-05-19 2024-09-27 道本妙用科技(北京)有限公司 A pulse wave intelligent analysis method and system based on human order parameter model
CN114916934B (en) * 2022-05-19 2025-08-08 道本妙用科技(北京)有限公司 A method and system for emotion analysis based on pulse wave sequence parameter phase space
CN117951505B (en) * 2024-03-27 2024-06-04 剑博微电子(南京)有限公司 Noise reduction method and system for medical chip of Internet of things

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040232965A1 (en) * 2003-05-21 2004-11-25 Terry Kuo Pulse interval to voltage converter and conversion method thereof
CN102028457A (en) * 2010-11-24 2011-04-27 北京麦邦光电仪器有限公司 Pulse rate measuring method and ring type pulse rate measuring meter
US20150080670A1 (en) * 2011-03-04 2015-03-19 Flint Hills Scientific, L.L.C. Detecting, assessing and managing a risk of death in epilepsy
CN105286845A (en) * 2015-11-29 2016-02-03 浙江师范大学 Movement noise elimination method suitable for wearable heart rate measurement device

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101194834A (en) * 2006-12-05 2008-06-11 重庆博恩富克医疗设备有限公司 Bio-electrical impedance measuring method and apparatus
CN101385644B (en) * 2008-07-25 2011-06-15 沈阳中国医科大学医疗器械研制中心(有限公司) 12 lead wireless remote electrocardiograph monitoring system
WO2010140241A1 (en) * 2009-06-04 2010-12-09 富士通株式会社 Awaking degree judgment device, method for judging degree of awaking and awaking degree judgment program
EP2713867A4 (en) * 2011-05-24 2015-01-21 Univ California IMAGING SOURCE OF MAGNEETHONEPHALOGRAPHY
CN102217931A (en) * 2011-06-09 2011-10-19 李红锦 Method and device for acquiring heart rate variation characteristic parameter
CN102984103B (en) * 2011-09-07 2015-08-05 华为技术有限公司 Signal processing method in spread spectrum system and device
CN103417206B (en) * 2012-05-22 2015-08-26 中国科学院深圳先进技术研究院 ECG removes the method and system of Hz noise
JP2014122875A (en) * 2012-11-26 2014-07-03 Canon Inc Device and method for measuring layered object
CN103654758A (en) * 2013-12-23 2014-03-26 韩山师范学院 Anti-jamming heart rate measurement method
US9936886B2 (en) * 2014-06-09 2018-04-10 Stmicroelectronics S.R.L. Method for the estimation of the heart-rate and corresponding system
CN104095630B (en) * 2014-07-29 2016-02-10 杭州电子科技大学 A kind of based on the phase locked fatigue detection method of brain electricity
CN104125579B (en) * 2014-08-07 2017-07-11 桂林电子科技大学 A kind of frequency spectrum sensing method and device based on time domain energy Yu frequency domain spectra entropy
CN104569902B (en) * 2014-11-21 2017-04-19 国家电网公司 Digital type electric energy meter power consumption measuring device and method
CN104490373B (en) * 2014-12-17 2016-12-07 辛勤 The determination methods of pulse signal, judgment means and physiological parameter measuring device
CN104655929B (en) * 2015-01-04 2017-11-21 中国科学院物理研究所 A kind of digital time-frequency measuring method of time-domain signal and corresponding target identification method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040232965A1 (en) * 2003-05-21 2004-11-25 Terry Kuo Pulse interval to voltage converter and conversion method thereof
CN102028457A (en) * 2010-11-24 2011-04-27 北京麦邦光电仪器有限公司 Pulse rate measuring method and ring type pulse rate measuring meter
US20150080670A1 (en) * 2011-03-04 2015-03-19 Flint Hills Scientific, L.L.C. Detecting, assessing and managing a risk of death in epilepsy
CN105286845A (en) * 2015-11-29 2016-02-03 浙江师范大学 Movement noise elimination method suitable for wearable heart rate measurement device

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107088062A (en) * 2017-05-08 2017-08-25 歌尔科技有限公司 A kind of heart rate module, the electronic installation and heart rate acquisition method for gathering heart rate
CN107088062B (en) * 2017-05-08 2023-09-01 歌尔科技有限公司 Heart rate module, electronic device for collecting heart rate and heart rate collection method
CN109924966A (en) * 2017-12-19 2019-06-25 郝振龙 It is a kind of for moving when measure heart rate Intelligent bracelet device
CN110801210A (en) * 2019-11-06 2020-02-18 心核心科技(北京)有限公司 Pulse wave signal filtering method and device, readable medium and electronic equipment
CN113616217A (en) * 2021-10-12 2021-11-09 深圳市倍轻松科技股份有限公司 Method and device for generating baseline drift curve

Also Published As

Publication number Publication date
CN105816165B (en) 2019-12-13
CN110876615B (en) 2022-08-02
CN110876615A (en) 2020-03-13

Similar Documents

Publication Publication Date Title
CN105816165B (en) A real-time dynamic heart rate monitoring device and monitoring method
EP3478166B1 (en) On-demand heart rate estimation based on optical measurements
Yang et al. Estimation and validation of arterial blood pressure using photoplethysmogram morphology features in conjunction with pulse arrival time in large open databases
CN106037694B (en) A kind of continuous blood pressure measurer based on pulse wave
CN103027690B (en) Hypoperfusion oxyhemoglobin saturation measuring method based on self-correlation modeling method
Schäck et al. Computationally efficient heart rate estimation during physical exercise using photoplethysmographic signals
US20190298209A1 (en) Heartbeat detection
KR102532764B1 (en) Apparatus and method for estimating biophysiological rates
Li et al. Design of a continuous blood pressure measurement system based on pulse wave and ECG signals
Liu et al. Filtering-induced time shifts in photoplethysmography pulse features measured at different body sites: The importance of filter definition and standardization
EP3292813B1 (en) Method and device for processing bio-signals
CN107949321A (en) Temporal Interference Removal and Improved Heart Rate Measurement Tracking Mechanism
CN112089405B (en) Pulse wave characteristic parameter measuring and displaying device
CN106994010A (en) A kind of heart rate detection method and system based on PPG signals
CN105286845A (en) Movement noise elimination method suitable for wearable heart rate measurement device
JP2023510943A (en) System and method for pulse transit time measurement with optical data
Tang et al. PPG signal reconstruction using a combination of discrete wavelet transform and empirical mode decomposition
Hong et al. Aging index using photoplethysmography for a healthcare device: comparison with brachial-ankle pulse wave velocity
Li et al. A new signal decomposition to estimate breathing rate and heart rate from photoplethysmography signal
CN106264505A (en) A kind of heart rate spectral peak system of selection based on support vector machine
CN105852863B (en) Respiration rate measuring method and device
CN110292372B (en) Detection device
CN106539580B (en) A Continuous Monitoring Method for Dynamic Changes of Autonomic Nervous System
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
CN103027692B (en) Dynamic spectrum data processing method based on uncertainty

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