CN108652604B - 一种基于心电信号的无气囊血压检测方法及系统 - Google Patents

一种基于心电信号的无气囊血压检测方法及系统 Download PDF

Info

Publication number
CN108652604B
CN108652604B CN201810119321.1A CN201810119321A CN108652604B CN 108652604 B CN108652604 B CN 108652604B CN 201810119321 A CN201810119321 A CN 201810119321A CN 108652604 B CN108652604 B CN 108652604B
Authority
CN
China
Prior art keywords
blood pressure
electrocardiosignals
characteristic
electrocardiosignal
prrx
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
CN201810119321.1A
Other languages
English (en)
Other versions
CN108652604A (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.)
Peking University Shenzhen Graduate School
Original Assignee
Peking University Shenzhen Graduate School
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 Peking University Shenzhen Graduate School filed Critical Peking University Shenzhen Graduate School
Priority to CN201810119321.1A priority Critical patent/CN108652604B/zh
Priority to PCT/CN2018/088430 priority patent/WO2019153579A1/zh
Publication of CN108652604A publication Critical patent/CN108652604A/zh
Application granted granted Critical
Publication of CN108652604B publication Critical patent/CN108652604B/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/021—Measuring pressure in heart or blood vessels
    • 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/02405—Determining heart rate variability
    • 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/7235—Details of waveform analysis
    • A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Vascular Medicine (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种基于心电信号的无气囊血压检测方法及系统。本发明通过获取心电信号，及获取对应心电信号的pRRx序列，通过对pRRx序列进行线性分析和/或非线性分析，得到相应特征指标，将计算得到的特征指标和对应的血压值作为输入和标签，进行机器学习，训练得到心电信号的特征指标与血压值对应关系的模型函数；当要检测某个时间点的血压值，通过获取该时间点之前的心电信号，计算并根据心电信号的特征指标，通过所述模型函数，得到该时间点的血压值。与现有技术相比，通过无创、无气囊采集心电信号作为源信号，成本低、安全有效、易于连续测量操作、用户体验好，而且本方法的血压检测过程计算量较小，算法复杂程度低，效率高。

Description

一种基于心电信号的无气囊血压检测方法及系统
技术领域
本发明涉及无气囊血压检测技术领域，具体涉及一种基于心电信号的无气囊血压检测方法及系统。
背景技术
血压(Blood Pressure，BP)作为人体一个重要生理参数，是人体心血管疾病的重要判断依据。血压测量方法分为直接测量法和间接测量法两种。
直接测量为测量系统直接与血液接触的测量，因其损伤皮肤和血管又被称为损伤性测量。
间接测量法，也称非损伤性测量，又可细分为间歇式测量法和连续式测量法。间歇式气囊测量法已经历了100多年的历史，所测血压基本接近主动脉内压力(动脉张力法)，是临床诊断过程中最常用、最普遍的检查方法，同时，因气囊对肌肉和血管有较大的挤压力，该方法无法扩展应用到连续测量。
在无创连续血压测量方法中，主要包括：动脉张力法与容积补偿法。其中，动脉张力法长时间测量血压时要保持传感器测量位置相对固定较为困难，容积补偿法长时间测量会导致静脉充血，给被检测者带来不舒适感甚至压痛感，且测量装置较为复杂。
发明内容
本发明主要解决的技术问题是常用的气囊间接测量血压方法，因气囊对肌肉和血管有较大的挤压力，无法扩展应用到连续测量，现有间接无创连续血压检测方法也存在操作困难、用户体验差等问题。
为解决上述技术问题，本发明提出一种新的间接血压测量方法，即：一种基于心电信号的无气囊血压检测方法，包括：获取心电信号；根据所述心电信号，计算对应血压值。
另一方面，本发明还提出一种基于心电信号的无气囊血压检测系统，包括：心电信号采集装置，用于采集待检测者的心电信号；处理器，用于执行如上所述的方法。
另一方面，本发明还提出一种基于心电信号的无气囊血压检测产品，包括：存储器，用于存储程序；处理器，用于通过执行所述存储器存储的程序以实现如上所述的方法。
另一方面，本发明还提出一种计算机可读存储介质，包括程序，所述程序能够被处理器执行以实现如上所述的方法。
本发明采用的基于心电信号的无气囊血压检测方法与现有技术相比，通过无创、无气囊采集心电信号作为源信号，成本低、安全有效、易于连续测量操作、用户体验好，而且本方法的血压检测过程计算量较小，算法复杂程度低，效率高。
附图说明
图1为一种基于心电信号的无气囊血压检测方法流程图；
图2为一种心电信号的特征指标与血压值对应关系的模型函数建立方法流程图；
图3为一种基于心电信号的无气囊血压检测系统示意图；
图4为一种基于心电信号的无气囊血压检测产品示意图。
具体实施方式
下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中，很多细节描述是为了使得本申请能被更好的理解。然而，本领域技术人员可以毫不费力的认识到，其中部分特征在不同情况下是可以省略的，或者可以由其他元件、材料、方法所替代。在某些情况下，本申请相关的一些操作并没有在说明书中显示或者描述，这是为了避免本申请的核心部分被过多的描述所淹没，而对于本领域技术人员而言，详细描述这些相关操作并不是必要的，他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。
另外，说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时，方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此，说明书和附图中的各种顺序只是为了清楚描述某一个实施例，并不意味着是必须的顺序，除非另有说明其中某个顺序是必须遵循的。
本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。而本申请所说“连接”、“联接”，如无特别说明，均包括直接和间接连接(联接)。
本发明提出的基于心电信号的无气囊血压检测方法主要基于心电信号的RR间隔序列，所述RR间隔是指心电信号波形中相邻的R峰和R峰之间的时间间隔，RR间隔序列包括一段心电信号中的所有RR间隔。
本发明实施例一：请参照图1，一种基于心电信号的无气囊血压检测方法，其包括A000步骤～A100步骤，下面具体说明：
A000：获取待检测者的心电信号。
A100：根据所述心电信号，计算对应血压值。
在一实施例中，A100步骤包括：根据心电信号，计算心电信号的一个或多个特征指标，根据心电信号的特征指标，计算对应血压值。
在一实施例中，心电信号的特征指标，包括：对心电信号的pRRx序列进行线性分析以得到一个或多个线性的特征指标，和/或进行非线性分析，以得到一个或多个非线性的特征指标。其中任意一段心电信号的pRRx序列通过以下方式计算得到：计算该段心电信号中相邻RR间期之差大于阈值x毫秒的数量与全部RR间期的数量的比值，通过设置值不同的阈值x，得到每一个阈值x对应的比值，这些比值构成了所述pRRx序列。在本实施例中，该比值用百分比表示，如式(1)所示：
Figure BDA0001571596950000031
根据所述心电信号的pRRx序列进行线性分析和/或非线性分析，可以得到一个或多个特征指标。
例如，线性分析获得的特征指标可以包括：pRRx序列的均值AVRR、pRRx序列的标准差SDRR、pRRx序列中相邻pRRx差值的均方根rMSSD、pRRx序列中相邻pRRx差值的标准差SDSD。
对每段心电信号的pRRx序列进行非线性分析，采用熵值分析法，即：根据现有技术，对于概率分布函数p(x)的随机变量集A，熵的定义如式(2)所示：
H(A)＝-∑pA(x)logpA(x) (2)
可以获得的特征指标包括：
(1)pRRx序列直方分布信息熵Sdh是对pRRx序列的数值分布信息熵；
(2)pRRx序列功率谱直方分布信息熵Sph是对pRRx序列进行离散傅里叶变换得到功率谱，然后根据功率谱序列的数值分布计算其信息熵；
(3)pRRx序列功率谱全频段分布信息熵Spf是对pRRx序列进行离散傅里叶变换得到功率谱，在全频段[fs/N,fs/2](信号的采样频率为fs，采样点数为N)内插入i-1个分点f1，f2，...，fm-1，将全频段分割成i个子频段。把每个频段内的功率密度之和作为该频段的功率密度，则得到m个功率密度。将这i个功率密度归一化得到每个频段出现的概率pi，则∑ipi＝1，相应的功率谱全频段熵如式(3)所示：
Figure BDA0001571596950000041
对每段心电信号的pRRx序列进行非线性分析，也可以采用下面四种分形维数计算分析方法可以得到如下的特征指标：
(1)结构函数法计算所得的分形维数Dsf，其中，结构函数法是指对于给定的序列z(x)，定义增量方差为结构函数，其关系为：
Figure BDA0001571596950000042
对于若干个标度τ，对序列z(x)的离散值计算出相应的S(τ)，然后画出logS(τ)-logτ的函数曲线，在无标度区进行线性拟合，得到斜率α，则对应分形维数Dsf与斜率α的转化关系如式(5)所示：
Figure BDA0001571596950000043
(2)相关函数法计算所得的分形维数Dcf，其中，相关函数法是指对于给定的序列z(x)，相关函数C(τ)定义为式(6)所示：
C(τ)＝AVE(z(x+τ)*z(x))，τ＝1,2,3,...,N-1 (6)
其中，AVE(·)表示平均，τ表示两点距离。此时相关函数为幂型，由于不存在特征长度，则分布为分形，有C(τ)ατ-α。这时，画出logC(τ)-logτ的函数曲线，在无标度区进行线性拟合，得到斜率α，则对应分形维数Dcf与斜率α的转化关系如式(7)所示：
Dcf＝2-α (7)
(3)变差法计算所得的分形维数Dvm，其中，变差法用宽为τ的矩形框首尾相接的将分形曲线覆盖起来，令第i个框内曲线的最大值和最小值之差为H(i)，即为矩形的高度。将所有矩形的高和宽相乘得到总面积S(τ)。改变τ的大小，得到一系列的S(τ)。如式(8)所示：
Figure BDA0001571596950000044
画出logN(τ)-logτ的函数曲线，在无标度区进行线性拟合得到斜率α，则对应分形维数Dvm与斜率α的转化关系如式(7)所示。
(4)均方根法计算所得的分形维数Drms，其中，均方根法用宽为τ的矩形框首尾相接的将分形曲线覆盖起来，令第i个框内曲线的最大值和最小值之差为H(i)，即为矩形的高度。计算这些矩形高度的均方根值S(τ)。改变τ的大小，得到一系列的S(τ)。画出logS(τ)-logτ的函数曲线，在无标度区进行线性拟合得到斜率α，则对应分形维数Drms与斜率α的转化关系如式(7)所示。
用于进行血压值计算的心电信号特征指标是上述线性和/或非线性分析得到的特征指标中的一个、多个，或者是其中几个的集合，也可以是除本实施例所罗列之外的现有分析方法所得到的相应特征指标。
在一实施例中，A100步骤在根据心电信号的特征指标来计算对应血压值时，可以预先建立心电信号的特征指标与血压值对应关系的模型函数，将心电信号的特征指标输入模型函数，得到对应血压值。例如，A100步骤可以通过机器学习和训练，来建立心电信号的特征指标与血压值对应关系的模型函数，请参照图2所示。
如图2所示，A100步骤建立上述模型函数，可以包括A110～A112步骤，下面具体说明。
A110：预先获取若干个血压值，及每个血压值的时间点之前的一段心电信号。其中，所述获取若干个血压值，包括运动与静坐、不同情绪状态、食用降压药前后、早晨与下午、不同睡眠状态等多个时间点的血压值，也可以根据需要增加获取血压值的时间点；这个步骤中所述获取血压值的方法可以采用现有技术中常用的、精准度高的方法，例如有创或气囊血压仪的检测结果，同时，对应每个血压值需要获取心电信号，由于个体新陈代谢情况存在差异，每个采样者所需的心电信号时间长度并不相同，以实际建模效果为准，本实施例选取1～30分钟不同时间长度的心电信号。
A111：获取这些心电信号的特征指标。
A112：将这些心电信号的特征指标作为输入，这些心电信号对应的血压值作为标签，进行机器学习，训练得到心电信号的特征指标与血压值对应关系的模型函数。其中，所述预先获取的若干血压值和心电信号都取自同一待测者，得到的血压检测模型也用于同一待测者的无气囊血压检测。此外，针对每个人的舒张压和收缩压需要单独建立模型函数。而且，在本实施例中，进行机器学习时，输出的血压值按照步长为5mmHg进行区间划分，舒张压区间为40～150mmHg，收缩压区间为50～300mmHg，从而将收缩压和舒张压划分成了若干个区间。
根据上述步骤得到心电信号的特征指标与血压值对应关系的模型函数后，再将A000步骤所获取待检测者的心电信号输入该模型函数，即可得到血压值，完成无气囊血压检测。连续获取几段心电信号输入该模型函数，即可得到对应连续检测的血压值。在本实施例中，收缩压和舒张压具体输出取值为此区间的中值并向上取整，例如：通过模型测得收缩压[70，75]，区间中值为72.5，向上取整为73，则输出的收缩压为数值73mmHg。
实施例二：一种基于心电信号的无气囊血压检测系统，如图3所示，包括心电信号采集装置B00和处理器B10，下面具体说明：
心电信号采集装置B00，用于采集待检测者的心电信号；
处理器B10，用于执行上述任一实施例所述的基于心电信号的无气囊血压检测方法。例如，处理器B10可以根据心电信号，计算心电信号的一个或多个特征指标，根据心电信号的特征指标，计算对应血压值。另一方面，处理器B10可以预先建立心电信号的特征指标与血压值对应关系的模型函数，将心电信号的特征指标输入模型函数，得到对应血压值。处理器B10通过预先获取若干个血压值，及每个血压值的时间点之前的心电信号；获取这些心电信号的特征指标；将这些心电信号的特征指标作为输入，这些心电信号对应的血压值作为标签，进行机器学习，训练得到心电信号的特征指标与血压值对应关系的模型函数。
实施例三：一种基于心电信号的无气囊血压检测产品C00，如图4所示，包括存储器C01和处理器C02，下面具体说明：
存储器C01，用于存储程序；
处理器C02，用于通过执行所述存储器存储的程序以实现上述任一实施例所述的基于心电信号的无气囊血压检测方法。例如，处理器C02执行存储器C01中存储的程序，可以根据心电信号，计算心电信号的一个或多个特征指标，根据心电信号的特征指标，计算对应血压值。另一方面，存储器C01中存储的程序还可以用于预先建立心电信号的特征指标与血压值对应关系的模型函数，将心电信号的特征指标输入模型函数，得到对应血压值。另一方面，处理器C02执行存储器C01中存储的程序，通过预先获取若干个血压值，及每个血压值的时间点之前的心电信号；获取这些心电信号的特征指标；将这些心电信号的特征指标作为输入，这些心电信号对应的血压值作为标签，进行机器学习，训练得到心电信号的特征指标与血压值对应关系的模型函数。
通过结合实施例一所述方法，使用实施例二系统中的装置，可以基于心电信号无创、无气囊检测得到血压值。这样的装置成本低、安全有效、易于连续测量操作、用户体验好，而且血压检测过程计算量较小，算法复杂程度低，效率高。
本领域技术人员可以理解，上述实施方式中各种方法的全部或部分功能可以通过硬件的方式实现，也可以通过计算机程序的方式实现。当上述实施方式中全部或部分功能通过计算机程序的方式实现时，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器、随机存储器、磁盘、光盘、硬盘等，通过计算机执行该程序以实现上述功能。例如，将程序存储在设备的存储器中，当通过处理器执行存储器中程序，即可实现上述全部或部分功能。另外，当上述实施方式中全部或部分功能通过计算机程序的方式实现时，该程序也可以存储在服务器、另一计算机、磁盘、光盘、闪存盘或移动硬盘等存储介质中，通过下载或复制保存到本地设备的存储器中，或对本地设备的系统进行版本更新，当通过处理器执行存储器中的程序时，即可实现上述实施方式中全部或部分功能。
以上应用了具体个例对本发明进行阐述，只是用于帮助理解本发明，并不用以限制本发明。对于本发明所属技术领域的技术人员，依据本发明的思想，还可以做出若干简单推演、变形或替换。

Claims (15)

1.一种基于心电信号的无气囊血压检测系统，其特征在于，包括：
心电信号采集装置，用于采集待检测者的心电信号；
处理器，用于执行如下方法：
获取心电信号；
根据所述心电信号，计算对应血压值；其中，所述根据心电信号，计算对应血压值包括：根据心电信号，计算心电信号的一个或多个特征指标，根据心电信号的特征指标，计算对应血压值；其中，心电信号的特征指标，包括：对心电信号的pRRx序列进行线性分析以得到一个或多个线性的特征指标，和/或进行非线性分析，以得到一个或多个非线性的特征指标；其中任意一段心电信号的pRRx序列通过以下方式计算得到：计算该段心电信号中相邻RR间期之差大于阈值x毫秒的数量与全部RR间期的数量的比值，通过设置值不同的阈值x，得到每一个阈值x对应的比值，这些比值构成了所述pRRx序列。
2.如权利要求1所述系统，其特征在于，所述处理器还用于执行如下方法：
预先建立心电信号的特征指标与血压值对应关系的模型函数，将心电信号的特征指标输入模型函数，得到对应血压值。
3.如权利要求1或2所述系统，其特征在于，心电信号的特征指标，还包括：
所述线性分析获得的特征指标：pRRx序列的均值AVRR、pRRx序列的标准差SDRR、pRRx序列中相邻pRRx差值的均方根rMSSD、pRRx序列中相邻pRRx差值的标准差SDSD中的至少一者；和/或，
所述非线性的特征指标包括对所述pRRx序列进行熵值分析法所得到的特征指标，包括：pRRx序列直方分布信息熵Sdh、pRRx序列功率谱直方分布信息熵Sph、pRRx序列功率谱全频段分布信息熵Spf中的至少一者；和/或，所述非线性的特征指标包括所述pRRx序列进行分形维数计算分析所得到的特征指标，包括：结构函数法计算所得的分形维数Dsf、相关函数法计算所得的分形维数Dcf、变差法计算所得的分形维数Dvm、均方根法计算所得的分形维数Drms中的至少一者。
4.如权利要求2所述系统，其特征在于，所述预先建立心电信号的特征指标与血压值对应关系的模型函数，包括：
预先获取若干个血压值，及每个血压值的时间点之前的心电信号；
获取这些心电信号的特征指标；
将这些心电信号的特征指标作为输入，这些心电信号对应的血压值作为标签，进行机器学习，训练得到心电信号的特征指标与血压值对应关系的模型函数。
5.如权利要求4所述系统，其特征在于，通过机器学习建立模型函数时，输出的血压值按照步长为5mmHg进行区间划分，具体输出的血压值为此区间的中值并向上取整。
6.一种基于心电信号的无气囊血压检测产品，其特征在于，包括：
存储器，用于存储程序；
处理器，用于通过执行所述存储器存储的程序以实现如下方法：
获取心电信号；
根据所述心电信号，计算对应血压值；其中，所述根据心电信号，计算对应血压值包括：根据心电信号，计算心电信号的一个或多个特征指标，根据心电信号的特征指标，计算对应血压值；其中，心电信号的特征指标，包括：对心电信号的pRRx序列进行线性分析以得到一个或多个线性的特征指标，和/或进行非线性分析，以得到一个或多个非线性的特征指标；其中任意一段心电信号的pRRx序列通过以下方式计算得到：计算该段心电信号中相邻RR间期之差大于阈值x毫秒的数量与全部RR间期的数量的比值，通过设置值不同的阈值x，得到每一个阈值x对应的比值，这些比值构成了所述pRRx序列。
7.如权利要求6所述产品，其特征在于，所述处理器还用于执行如下方法：
预先建立心电信号的特征指标与血压值对应关系的模型函数，将心电信号的特征指标输入模型函数，得到对应血压值。
8.如权利要求6或7所述产品，其特征在于，心电信号的特征指标，还包括：
所述线性分析获得的特征指标：pRRx序列的均值AVRR、pRRx序列的标准差SDRR、pRRx序列中相邻pRRx差值的均方根rMSSD、pRRx序列中相邻pRRx差值的标准差SDSD中的至少一者；和/或，
所述非线性的特征指标包括对所述pRRx序列进行熵值分析法所得到的特征指标，包括：pRRx序列直方分布信息熵Sdh、pRRx序列功率谱直方分布信息熵Sph、pRRx序列功率谱全频段分布信息熵Spf中的至少一者；和/或，所述非线性的特征指标包括所述pRRx序列进行分形维数计算分析所得到的特征指标，包括：结构函数法计算所得的分形维数Dsf、相关函数法计算所得的分形维数Dcf、变差法计算所得的分形维数Dvm、均方根法计算所得的分形维数Drms中的至少一者。
9.如权利要求7所述产品，其特征在于，所述预先建立心电信号的特征指标与血压值对应关系的模型函数，包括：
预先获取若干个血压值，及每个血压值的时间点之前的心电信号；
获取这些心电信号的特征指标；
将这些心电信号的特征指标作为输入，这些心电信号对应的血压值作为标签，进行机器学习，训练得到心电信号的特征指标与血压值对应关系的模型函数。
10.如权利要求9所述产品，其特征在于，通过机器学习建立模型函数时，输出的血压值按照步长为5mmHg进行区间划分，具体输出的血压值为此区间的中值并向上取整。
11.一种计算机可读存储介质，其特征在于，包括程序，所述程序能够被处理器执行以实现如下方法：
获取心电信号；
根据所述心电信号，计算对应血压值；其中，所述根据心电信号，计算对应血压值包括：根据心电信号，计算心电信号的一个或多个特征指标，根据心电信号的特征指标，计算对应血压值；其中，心电信号的特征指标，包括：对心电信号的pRRx序列进行线性分析以得到一个或多个线性的特征指标，和/或进行非线性分析，以得到一个或多个非线性的特征指标；其中任意一段心电信号的pRRx序列通过以下方式计算得到：计算该段心电信号中相邻RR间期之差大于阈值x毫秒的数量与全部RR间期的数量的比值，通过设置值不同的阈值x，得到每一个阈值x对应的比值，这些比值构成了所述pRRx序列。
12.如权利要求11所述介质，其特征在于，所述处理器还用于执行如下方法：
预先建立心电信号的特征指标与血压值对应关系的模型函数，将心电信号的特征指标输入模型函数，得到对应血压值。
13.如权利要求11或12所述介质，其特征在于，心电信号的特征指标，还包括：
所述线性分析获得的特征指标：pRRx序列的均值AVRR、pRRx序列的标准差SDRR、pRRx序列中相邻pRRx差值的均方根rMSSD、pRRx序列中相邻pRRx差值的标准差SDSD中的至少一者；和/或，
所述非线性的特征指标包括对所述pRRx序列进行熵值分析法所得到的特征指标，包括：pRRx序列直方分布信息熵Sdh、pRRx序列功率谱直方分布信息熵Sph、pRRx序列功率谱全频段分布信息熵Spf中的至少一者；和/或，所述非线性的特征指标包括所述pRRx序列进行分形维数计算分析所得到的特征指标，包括：结构函数法计算所得的分形维数Dsf、相关函数法计算所得的分形维数Dcf、变差法计算所得的分形维数Dvm、均方根法计算所得的分形维数Drms中的至少一者。
14.如权利要求12所述介质，其特征在于，所述预先建立心电信号的特征指标与血压值对应关系的模型函数，包括：
预先获取若干个血压值，及每个血压值的时间点之前的心电信号；
获取这些心电信号的特征指标；
将这些心电信号的特征指标作为输入，这些心电信号对应的血压值作为标签，进行机器学习，训练得到心电信号的特征指标与血压值对应关系的模型函数。
15.如权利要求14所述介质，其特征在于，通过机器学习建立模型函数时，输出的血压值按照步长为5mmHg进行区间划分，具体输出的血压值为此区间的中值并向上取整。
CN201810119321.1A 2018-02-06 2018-02-06 一种基于心电信号的无气囊血压检测方法及系统 Active CN108652604B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201810119321.1A CN108652604B (zh) 2018-02-06 2018-02-06 一种基于心电信号的无气囊血压检测方法及系统
PCT/CN2018/088430 WO2019153579A1 (zh) 2018-02-06 2018-05-25 一种基于心电信号的无气囊血压检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810119321.1A CN108652604B (zh) 2018-02-06 2018-02-06 一种基于心电信号的无气囊血压检测方法及系统

Publications (2)

Publication Number Publication Date
CN108652604A CN108652604A (zh) 2018-10-16
CN108652604B true CN108652604B (zh) 2020-03-10

Family

ID=63785078

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810119321.1A Active CN108652604B (zh) 2018-02-06 2018-02-06 一种基于心电信号的无气囊血压检测方法及系统

Country Status (2)

Country Link
CN (1) CN108652604B (zh)
WO (1) WO2019153579A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113520350B (zh) * 2021-07-27 2024-05-17 香港心脑血管健康工程研究中心有限公司 获取血压图信号的相关特征参数和指标信息的处理方法及装置
CN113854985A (zh) * 2021-08-27 2021-12-31 联卫医疗科技(上海)有限公司 一种用于血压预测的机器学习模型样本获取的方法及装置
CN117503085B (zh) * 2023-11-07 2024-06-14 西康软件有限责任公司 血压数据评估方法、装置、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105725999A (zh) * 2016-01-26 2016-07-06 四川长虹电器股份有限公司 一种基于脉搏波传导时间的血压测量方法
WO2016122054A1 (ko) * 2015-01-27 2016-08-04 주식회사 휴이노 혈압을 실시간으로 모니터링하기 위한 방법, 시스템 및 비일시성의 컴퓨터 판독 가능한 기록 매체
CN106539572A (zh) * 2016-11-02 2017-03-29 中国科学院电子学研究所 一种多参数融合的连续性血压测量方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100467056B1 (ko) * 2002-08-31 2005-01-24 (주)유인바이오테크 자동혈압측정장치 및 방법
CN100361625C (zh) * 2004-03-26 2008-01-16 香港中文大学 非侵入式血压测量装置及方法
CN1698536A (zh) * 2004-05-20 2005-11-23 香港中文大学 采用自动补偿的无袖带式连续血压测量方法
CN201320166Y (zh) * 2008-12-18 2009-10-07 北京京精医疗设备有限公司 指触式无袖带数字血压计
JP6388447B2 (ja) * 2015-03-23 2018-09-12 国立大学法人九州工業大学 生体信号処理装置及び血圧測定システム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016122054A1 (ko) * 2015-01-27 2016-08-04 주식회사 휴이노 혈압을 실시간으로 모니터링하기 위한 방법, 시스템 및 비일시성의 컴퓨터 판독 가능한 기록 매체
CN105725999A (zh) * 2016-01-26 2016-07-06 四川长虹电器股份有限公司 一种基于脉搏波传导时间的血压测量方法
CN106539572A (zh) * 2016-11-02 2017-03-29 中国科学院电子学研究所 一种多参数融合的连续性血压测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Measuring of Systolic Blood Pressure Based On Heart Rate;M.K. Ali Hassan 等;《4th Kuala Lumpur International Conference on Biomedical Engineering 2008》;20080628;第595-598页 *

Also Published As

Publication number Publication date
CN108652604A (zh) 2018-10-16
WO2019153579A1 (zh) 2019-08-15

Similar Documents

Publication Publication Date Title
CN110251105B (zh) 一种无创血压测量方法、装置、设备及系统
EP3117766B1 (en) Processing biological data
CN108652640B (zh) 一种基于心电信号的无创血糖检测方法及系统
EP0818175B1 (en) Living body condition measuring apparatus
CN110236573B (zh) 心理压力状态的检测方法及相关装置
CN102834047A (zh) 在心血管参数确定中不规则心动周期影响的排除
CN101765398B (zh) 前负荷依赖性和流体反应性的评估
CN106691406A (zh) 一种基于单探头光电容积脉搏波的血管弹性和血压的检测方法
CN108652604B (zh) 一种基于心电信号的无气囊血压检测方法及系统
EP2854620A1 (en) Narrow band feature extraction from cardiac signals
EP3785610B1 (en) Method of processing a signal representing a physiological rhythm
RU2011122787A (ru) Способы и системы для неинвазивного измерения уровней глюкозы
Wang et al. An advanced signal processing technique for impedance cardiography
Li et al. A new signal decomposition to estimate breathing rate and heart rate from photoplethysmography signal
CN111714088B (zh) 基于中医原理的人体特征指标检测方法和系统
CN112294272A (zh) 监护仪及其不规则脉率识别方法
CN105105734A (zh) 一种基于心音信号的无创连续血压测量方法、装置及系统
Park et al. Prediction of daily mental stress levels using a wearable photoplethysmography sensor
US20110237962A1 (en) Use of the frequency spectrum of artifact in oscillometry
CN104644151B (zh) 一种基于光电容积脉搏信号的压力脉搏波波形传播预测方法
C. Dela Cruz et al. Deriving heart rate and respiratory rate from ecg using wavelet transform
CN102271574B (zh) 推导中心主动脉收缩压力值的方法和用于分析动脉数据集以推导中心主动脉收缩压力值的方法
Peters et al. Beat-to-beat heart rate detection in multi-lead abdominal fetal ECG recordings
Soni et al. Effect of physical activities on heart rate variability and skin conductance
CN114145725B (zh) 一种基于无创连续血压测量的ppg采样率估算方法

Legal Events

Date Code Title Description
PB01 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