CN109124610A - 一种无创血压测量的抗干扰方法和装置 - Google Patents
一种无创血压测量的抗干扰方法和装置 Download PDFInfo
- Publication number
- CN109124610A CN109124610A CN201811120261.1A CN201811120261A CN109124610A CN 109124610 A CN109124610 A CN 109124610A CN 201811120261 A CN201811120261 A CN 201811120261A CN 109124610 A CN109124610 A CN 109124610A
- Authority
- CN
- China
- Prior art keywords
- value
- pulse wave
- pulse
- wave signal
- frequency
- 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
Links
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
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
-
- 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
- A61B5/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
- A61B5/0225—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers the pressure being controlled by electric signals, e.g. derived from Korotkoff sounds
-
- 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
-
- 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/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Mathematical Physics (AREA)
- Ophthalmology & Optometry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明公开了一种无创血压测量的抗干扰方法和装置,包括:获取脉搏数据中的压力波信号和脉搏波信号;对脉搏波进行分析,并判断脉搏波信号是否受到干扰,若是,则延长脉搏波信号的采样时间并对脉搏波信号进行频域分析获得频域幅度的极大值,将极大值作为脉搏波的第一包络序列值,否则对脉搏波进行时域搜索分析,获得脉搏波的第二包络序列值;根据第一包络序列值和/或第二包络序列值求取血压值。本发明通过综合分析时域形态和频域幅度,剔除掉干扰信息,提取出有效的脉搏波信息,提高了抗干扰性能,得到更为准确的无创血压测量结果。
Description
技术领域
本发明涉及血压测量技术领域,更具体地说,涉及一种无创血压测量的抗干扰方法和装置。
背景技术
血压作为人体重要的生理参数,在临床具有极其重要的意义。目前血压测量方法主要分为有创血压测量和无创血压测量两大类。其中有创血压测量准确,但是由于需要侵入人体,有无菌要求,所以其应用场景局限于手术室,ICU等特殊场合。而无创血压具备操作简单,使用方便的特点,适合各种安静场合下的血压测量。
然而,由于无创血压测量不侵入人体内,在测量过程中若测量患者不安静,则无创血压非常容易受到测量患者走动、肢体动作等所带来的干扰,从而影响血压的正常测量,无法获得准确的血压值。
发明内容
针对上述无创血压测量容易受到干扰而导致测量结果不准确的缺陷,本发明提供了一种无创血压测量的抗干扰方法和装置。
一方面,本发明提供的一种无创血压测量的抗干扰方法,包括:
获取脉搏数据中的压力波信号和脉搏波信号;
对所述脉搏波进行分析,并判断所述脉搏波信号是否受到干扰,若是,则延长所述脉搏波信号的采样时间并对所述脉搏波信号进行频域分析获得频域幅度的极大值,将所述极大值作为所述脉搏波的第一包络序列值,否则对所述脉搏波进行时域搜索分析,获得所述脉搏波的第二包络序列值;
根据所述第一包络序列值和/或所述第二包络序列值求取血压值。
可选地,所述对所述脉搏波信号进行频域分析获得频域幅度的极大值具体包括:
通过基2时间抽取快速傅里叶变换对所述脉搏波信号进行频域分析获得频域幅度的极大值。
可选地,所述通过基2时间抽取快速傅里叶变换对所述脉搏波信号进行频域分析获得频域幅度的极大值包括:
对所述脉搏波信号通过第一公式进行傅里叶变换,所述第一公式为:
其中,脉搏波序列为x(n),N=2L,L为自然数且为蝶形迭代运算的阶数;
将所述第一公式按n的奇偶性分解为第二公式:
令x1(m)=x(2m),x2(m)=x(2m+1),且由于将所述第二公式转换为第三公式:
用F1(k)和F2(k)分别表示x1(m)和x2(m)的N/2点时间抽取快速傅里叶变换,获得第四公式和第五公式:
将所述第四公式和所述第五公式代入所述第三公式,并根据以及F1(k)、F2(k)的隐含周期性得到第六公式:
计算所述第六公式,并进行迭代分解为N/4点的抽取快速傅里叶变换,N/8点的抽取快速傅里叶变换,…,依次类推,总共迭代分解L次,完成整个基2抽取快速傅里叶变换,获得频域幅度的极大值。
可选地,在对同一个患者的脉搏波进行分析时,获取所述患者血压测量的历史记录中的第一脉率值作为参考脉率值,并在所述参考脉率值的±5bpm范围内,搜索脉搏波频域幅度的极大值,作为第三包络序列值;
根据所述第三包络序列值求取血压值。
可选地,在对无干扰的脉搏波进行时域搜索分析并获取到第二包络序列值的周期信息后,求取对应的第二脉率值;
将所述第二脉率值与历史记录中的参考脉率值进行比较,若所述第二脉率值和所述参考脉率值之间的偏差值大于预设差值,则将所述第二脉率值更正为新的参考脉率值。
可选地,所述判断所述脉搏波信号是否受到干扰包括:
判断同一压力平台期内的预设数量个脉搏波的幅度值、周期值、脉宽值和基线值是否均满足预设相似条件,若是,则判定为未受到干扰,否则判定为受到干扰;所述预设数量个为至少两个。
可选地,所述预设相似条件为脉搏波两两之间的幅度值差、周期值差、脉宽值差和基线值差均分别小于对应的预设差值。
可选地,所述延长所述脉搏波信号的采样时间具体为延长受到干扰的脉搏波信号的压力平台期的时间,直至频域分析中获得预置数量个采样点。
可选地,所述获取脉搏数据中的压力波信号和脉搏波信号包括:
获取袖套中的原始压力信号,并对所述原始压力信号进行低通滤波获得压力波信号,对所述原始压力信号进行带通滤波获得脉搏波信号;
或,
通过AD采样器采集获取袖套中的原始压力信号,并对所述原始压力信号进行数字低通滤波分离获得压力波信号,对所述原始压力信号进行数字带通滤波分离获得脉搏波。
另一方面,本发明提供的一种无创血压测量的抗干扰装置,包括:
获取模块,用于获取脉搏数据中的压力波信号和脉搏波信号;
分析模块,用于对所述脉搏波进行分析,并判断所述脉搏波信号是否受到干扰,若是,则延长所述脉搏波信号的采样时间并对所述脉搏波信号进行频域分析获得频域幅度的极大值,将所述极大值作为所述脉搏波的第一包络序列值,否则对所述脉搏波进行时域搜索分析,获得所述脉搏波的第二包络序列值;
血压值求取模块,用于根据所述第一包络序列值和/或所述第二包络序列值求取血压值。
本发明具有以下有益效果:
本发明中在获取到压力波信号和脉搏波信号之后,对脉搏波是否受到干扰进行判断,当脉搏波未受到干扰时对脉搏波进行时域的形态分析,获得脉搏波的时域波形特征;当脉搏波受到干扰时对脉搏波进行频域分析,获得脉搏波的幅频特征,通过综合分析时域形态和频域幅度,剔除掉干扰信息,提取出有效的脉搏波信息,提高了抗干扰性能,得到更为准确的无创血压测量结果。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1为本发明实施例提供的一种无创血压测量的抗干扰方法的流程示意图;
图2为本发明实施例提供的一种时域搜索脉搏波的示意图;
图3为本发明实施例提供的一种频域分析脉搏波的幅频特征的示意图;
图4为本发明实施例提供的一种无创血压测量的抗干扰装置的结构示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
经发明人研究发现,目前常用的无创血压测量方法主要是振荡法,振荡法是通过袖带充气升压超过收缩压以阻塞手臂的血流,再逐步放气降压,利用电子传感器采集放气过程中的压力数据,手臂动脉血流的脉动传递到袖带内产生压力波和脉搏波,从而得到一组幅度由小变大,再由大变小的压力脉搏波数据及对应的压力波数据,找到最大脉搏幅度Am,其所对应的压力为平均压Pm,收缩压脉搏幅度As=Am*Ks,其对应压力为收缩压,舒张压脉搏幅度Ad=Am*Kd,其对应压力为舒张压。振荡法可以自动测量,方便操作,但测量准确性依赖于经验系数Ks和Kd的选定。此外,在传统常用的振荡法测量中,脉搏波是幅度很小的交流信号,波形形态表征复杂,非常容易受到受试者走动,肢体动作等等各种干扰。并且,传统的振荡法都是对压力平台期脉搏波的形态进行识别和判断,即通过识别分析脉搏波的特征形态,比如幅度、周期、到峰值时间、斜率、收缩面积、舒张面积等等特征量,来测定用于脉搏的特征测量的各个品质值,最终利用品质值确定血压值和脉率值。诸如此类的无创血压的抗干扰方法,都是利用时域的脉搏波形态进行分析和甄别,抑制干扰的效果并不明显,而且误识别的概率较大,因为仅仅依靠时域的波形形态,很难分离信号与干扰并识别出有效的脉搏幅度。
有鉴于此,本发明实施例中提供一种基于阶梯式放气的无创血压测量的抗干扰方法,以实现在干扰环境下实现血压值的准确测量。无创血压测量的抗干扰方法可以应用于电子血压计,通过将电子血压计上的袖带绑于患者的手臂上,并且通过袖带上的传感器获取袖带在充气绑紧患者的手臂时的压力信号,对压力信号进行分析以获得血压值。其中,阶梯式放气为在充气阶段袖带压力连续增长,当达到一定值后,通过气阀进入放气阶段。在放气阶段中,气阀的短期开放,会导致袖带的压力突然下降,然后气阀再闭合一段时间,袖带的压力维持在一个稳定的数值上,通过连续开放气阀和闭合气阀,实现阶梯式放气。而气阀每次闭合的时候,袖带的压力均维持于一个对应的压力值上,即为对应的压力平台,而气阀在该压力值下闭合的时间段为压力平台期。
如图1所示,图1为本发明实施例提供的一种无创血压测量的抗干扰方法的流程示意图;本发明实施例提供的一种无创血压测量的抗干扰方法,包括:
S101、获取脉搏数据中的压力波信号和脉搏波信号;
在患者绑上袖带之后,可以通过袖带上的压力传感器获取袖带在充气绑紧患者的手臂时的压力信号,即获取袖套中的原始压力信号作为患者的脉搏数据,对所述原始压力信号进行低通滤波可以获得压力波信号,对所述原始压力信号进行带通滤波可以获得脉搏波信号。另外,还可以通过AD(Analog-Digita,模数)采样器采集获取袖套中的原始压力信号,对所述原始压力信号进行数字低通滤波分离获得压力波信号,对所述原始压力信号进行数字带通滤波分离获得脉搏波。
S102、对所述脉搏波进行分析,并判断所述脉搏波信号是否受到干扰;
S103、若是,则延长所述脉搏波信号的采样时间并对所述脉搏波信号进行频域分析获得频域幅度的极大值,将所述极大值作为所述脉搏波的第一包络序列值;
S104、否则对所述脉搏波进行时域搜索分析,获得所述脉搏波的第二包络序列值;
在脉搏波信号未受到明显的干扰的时候,可以通过时域的搜索分析脉搏波波形的形态,并且进行相关的匹配操作,从而获取在对应的压力平台期的包络序列值;而在脉搏波信号受到明显的干扰的时候,对脉搏波信号进行时域的搜索分析容易得到会影响最终的血压值计算的包络序列伪值,即很难依靠时域的搜索波形分析提取出有效的包络序列值,因此,此时可以通过延长脉搏波信号的采样时间以获得足够的采样点数,并对脉搏波信号进行频域分析,得到有效的包络序列值。具体地,可以参阅图2和图3,图2为本发明实施例提供的一种时域搜索脉搏波的示意图,图3为本发明实施例提供的一种频域分析脉搏波的幅频特征的示意图。
需要说明的是,血压测量是基于阶梯式放气过程进行的,即整个血压测量过程包括了阶梯式放气中的多个压力平台期。在多个压力平台期中,在脉搏波未受到干扰的情况下采用时域搜索分析,在脉搏波受到干扰的情况下采用频域分析,则在同一个压力平台期内可能存在同时进行的时域搜索分析和频域分析的情况。当在其中一个压力平台期进行频域分析时,由于进行频域分析所需的采样点较多,因此需要延长脉搏波信号的采样时间以获得足够的采样点数,即延长该压力平台期以使得采样时间足够长,具体为通过控制气阀闭合的时间,延长受到干扰的脉搏波信号的压力平台期的时间,直至频域分析中获得预置数量个采样点。其中,预置数量个采样点可以根据在实际测量情况下进行频域分析所需的采样点进行设置和调整,此处不做具体限定。
在本发明实施例提供的一种实施方式中,可以通过基2时间抽取快速傅里叶变换对所述脉搏波信号进行频域分析获得频域幅度的极大值。
具体地,所述通过基2时间抽取快速傅里叶变换(Decimation in time-FastFourier Transform,DFT)对所述脉搏波信号进行频域分析获得频域幅度的极大值包括:对所述脉搏波信号通过第一公式进行傅里叶变换,所述第一公式为:
其中,脉搏波序列为x(n),N为变换区间且N=2L,L为自然数且为蝶形迭代运算的阶数;
将所述第一公式按n的奇偶性分解为第二公式:
令x1(m)=x(2m),x2(m)=x(2m+1),且由于将所述第二公式转换为第三公式:
即,按n的奇偶性可以将x(n)分解为两个N/2长的序列x1(m)和x2(m),则N点DFT可分解为两个N/2点DFT来计算。用F1(k)和F2(k)分别表示x1(m)和x2(m)的N/2点DFT,获得第四公式和第五公式:
将所述第四公式和所述第五公式代入所述第三公式,并根据以及F1(k)、F2(k)的隐含周期性得到第六公式:
计算所述第六公式,并进行迭代分解为N/4点的抽取快速傅里叶变换,N/8点的抽取快速傅里叶变换,…,依次类推,总共迭代分解L次,完成整个基2抽取快速傅里叶变换,获得频域幅度的极大值。
通过迭代分解运算,可以使得DFT的运算量下降几个数量级,从而大大提高了数字信号处理的速度,使得在无创血压测量过程中,实现了压力平台期实时同步地分析受干扰的脉搏波的频域幅度。
需要说明的是,在本发明实施例中采用了基2时间抽取快速傅立叶变换对所述脉搏波信号进行频域分析,但不限于该具体方法,本发明实施例中还可以采用对等采用基2频分抽取快速傅立叶、多基时分或频分快速傅立叶,均可达到上述的效果。
S105、根据所述第一包络序列值和/或所述第二包络序列值求取血压值。
最后,结合压力波信号以及所述第一包络序列值和/或所述第二包络序列值求取血压值。根据在不同的压力波下第一包络序列值和/或所述第二包络序列值的幅度变化特征,利用经验准则来识别动脉收缩压、平均压和舒张压,获得血压值。
进一步地,在同一个患者的多次、长时的动态血压测量场景中,可以利用电子血压计中所储存的患者历史测量记录里的脉率值信息,提供更准确的参比信息,能够更有效的利用脉搏波频域分析进行干扰的剔除,提高测量准确性。因此,在对同一个患者的脉搏波进行分析时,可以获取所述患者血压测量的历史记录中的第一脉率值作为参考脉率值,并在所述参考脉率值的±5bpm范围内,搜索脉搏波频域幅度的极大值,作为第三包络序列值;最后再根据所述第三包络序列值求取血压值。需要说明的是,患者血压测量的历史记录中的第一脉率值为电子血压计在测量患者的血压值过程中,对无干扰的脉搏波进行时域搜索分析并获得的相应的脉率值,即该脉率值为未受干扰下测量获得的脉率值。
可以理解的是,在血压测量过程中,对无干扰平台期的脉搏波的形态特征进行分析后,可以与历史记录中的脉搏波形态特征相互比对,如果脉搏波形态有较大变异,则可以判定已更换患者或者血压有较大的变化,则需要对当前的参考脉率值进行更正,以便对有干扰平台期的脉搏波频域幅度进行更准确的提取。因此,在本发明实施例中,在每次测量过程中,对无干扰的脉搏波进行时域搜索分析并获取到第二包络序列值的周期信息后,求取对应的第二脉率值;将所述第二脉率值与历史记录中的参考脉率值进行比较,若所述第二脉率值和所述参考脉率值之间的偏差值大于预设差值,则将所述第二脉率值更正为新的参考脉率值。其中,预设差值可以为根据人体在正常情况下其脉率值的波动范围进行设置,其范围可以为15-25bpm,较优地可以采取20bpm,此处不做具体限定。
进一步地,在本发明实施例中,所述判断所述脉搏波信号是否受到干扰包括:判断同一压力平台期内的预设数量个脉搏波的幅度值、周期值、脉宽值和基线值是否均满足预设相似条件,若是,则判定为未受到干扰,否则判定为受到干扰;所述预设数量个为至少两个。具体地,所述预设相似条件为脉搏波两两之间的幅度值差、周期值差、脉宽值差和基线值差均分别小于对应的预设差值,其中预设差值包括有预设幅度值差、预设周期值差、预设脉宽值差和预设基线值差。可以理解的是,若在同一个平台期内多个脉搏波之间的幅度值、周期值、脉宽值和基线值的差值均不大,可以认为脉搏波未受到干扰,若多个脉搏波之间的幅度值、周期值、脉宽值和基线值的差值较大,则可以认为脉搏波受到了干扰从而导致脉搏波的波动较大。具体地,幅度值偏差不超出30%-50%,较优地,可以采用40%作为幅度值偏差值进行判断;周期值偏差不超出10-20bpm,较优地,可以采用15bpm作为周期值偏差值进行判断;脉宽值偏差不超出60-80ms,较优地,可以采用70ms作为脉宽值偏差值进行判断;基线值则不超出标准基线值的5%-15%,较优地,可以采用10%作为基线值偏差值进行判断。
另一方面,可以参阅图4,图4为本发明实施例提供的一种无创血压测量的抗干扰装置的结构示意图。本发明提供的一种无创血压测量的抗干扰装置,包括:
获取模块401,用于获取脉搏数据中的压力波信号和脉搏波信号;
分析模块402,用于对所述脉搏波进行分析,并判断所述脉搏波信号是否受到干扰,若是,则延长所述脉搏波信号的采样时间并对所述脉搏波信号进行频域分析获得频域幅度的极大值,将所述极大值作为所述脉搏波的第一包络序列值,否则对所述脉搏波进行时域搜索分析,获得所述脉搏波的第二包络序列值;
血压值求取模块403,用于根据所述第一包络序列值和/或所述第二包络序列值求取血压值。
可选地,分析模块402还包括:分析单元,用于通过基2时间抽取快速傅里叶变换对所述脉搏波信号进行频域分析获得频域幅度的极大值。
可选地,所述分析单元具体用于对所述脉搏波信号通过第一公式进行傅里叶变换,所述第一公式为:
其中,脉搏波序列为x(n),N=2L,L为自然数且为蝶形迭代运算的阶数;
将所述第一公式按n的奇偶性分解为第二公式:
令x1(m)=x(2m),x2(m)=x(2m+1),且由于将所述第二公式转换为第三公式:
用F1(k)和F2(k)分别表示x1(m)和x2(m)的N/2点时间抽取快速傅里叶变换,获得第四公式和第五公式:
将所述第四公式和所述第五公式代入所述第三公式,并根据以及F1(k)、F2(k)的隐含周期性得到第六公式:
计算所述第六公式,并进行迭代分解为N/4点的抽取快速傅里叶变换,N/8点的抽取快速傅里叶变换,…,依次类推,总共迭代分解L次,完成整个基2抽取快速傅里叶变换,获得频域幅度的极大值。
可选地,还包括:参考分析模块,用于在对同一个患者的脉搏波进行分析时,获取所述患者血压测量的历史记录中的第一脉率值作为参考脉率值,并在所述参考脉率值的±5bpm范围内,搜索脉搏波频域幅度的极大值,作为第三包络序列值;
根据所述第三包络序列值求取血压值。
可选地,还包括比较模块,用于在对无干扰的脉搏波进行时域搜索分析并获取到第二包络序列值的周期信息后,求取对应的第二脉率值;
将所述第二脉率值与历史记录中的参考脉率值进行比较,若所述第二脉率值和所述参考脉率值之间的偏差值大于预设差值,则将所述第二脉率值更正为新的参考脉率值。
可选地,分析模块402还包括:判断单元,用于判断同一压力平台期内的预设数量个脉搏波的幅度值、周期值、脉宽值和基线值是否均满足预设相似条件,若是,则判定为未受到干扰,否则判定为受到干扰;所述预设数量个为至少两个。
可选地,所述预设相似条件为脉搏波两两之间的幅度值差、周期值差、脉宽值差和基线值差均分别小于对应的预设差值。
可选地,所述延长所述脉搏波信号的采样时间具体为延长受到干扰的脉搏波信号的压力平台期的时间,直至频域分析中获得预置数量个采样点。
可选地,获取模块401包括:第一获取单元,用于获取袖套中的原始压力信号,并对所述原始压力信号进行低通滤波获得压力波信号,对所述原始压力信号进行带通滤波获得脉搏波信号;
第二获取单元,用于通过AD采样器采集获取袖套中的原始压力信号,并对所述原始压力信号进行数字低通滤波分离获得压力波信号,对所述原始压力信号进行数字带通滤波分离获得脉搏波。
本发明实施例提供的一种无创血压测量的抗干扰装置用于执行上述的无创血压测量的抗干扰方法,具体的执行步骤可以参阅上述的方法实施例,此处不再赘述。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。
Claims (10)
1.一种无创血压测量的抗干扰方法,其特征在于,包括:
获取脉搏数据中的压力波信号和脉搏波信号;
对所述脉搏波进行分析,并判断所述脉搏波信号是否受到干扰,若是,则延长所述脉搏波信号的采样时间并对所述脉搏波信号进行频域分析获得频域幅度的极大值,将所述极大值作为所述脉搏波的第一包络序列值,否则对所述脉搏波进行时域搜索分析,获得所述脉搏波的第二包络序列值;
根据所述第一包络序列值和/或所述第二包络序列值求取血压值。
2.根据权利要求1所述的无创血压测量的抗干扰方法,其特征在于,所述对所述脉搏波信号进行频域分析获得频域幅度的极大值具体包括:
通过基2时间抽取快速傅里叶变换对所述脉搏波信号进行频域分析获得频域幅度的极大值。
3.根据权利要求2所述的无创血压测量的抗干扰方法,其特征在于,所述通过基2时间抽取快速傅里叶变换对所述脉搏波信号进行频域分析获得频域幅度的极大值包括:
对所述脉搏波信号通过第一公式进行傅里叶变换,所述第一公式为:
其中,脉搏波序列为x(n),N=2L,L为自然数且为蝶形迭代运算的阶数;
将所述第一公式按n的奇偶性分解为第二公式:
令x1(m)=x(2m),x2(m)=x(2m+1),且由于将所述第二公式转换为第三公式:
用F1(k)和F2(k)分别表示x1(m)和x2(m)的N/2点时间抽取快速傅里叶变换,获得第四公式和第五公式:
将所述第四公式和所述第五公式代入所述第三公式,并根据以及F1(k)、F2(k)的隐含周期性得到第六公式:
计算所述第六公式,并进行迭代分解为N/4点的抽取快速傅里叶变换,N/8点的抽取快速傅里叶变换,…,依次类推,总共迭代分解L次,完成整个基2抽取快速傅里叶变换,获得频域幅度的极大值。
4.根据权利要求1所述的无创血压测量的抗干扰方法,其特征在于,在对同一个患者的脉搏波进行分析时,获取所述患者血压测量的历史记录中的第一脉率值作为参考脉率值,并在所述参考脉率值的±5bpm范围内,搜索脉搏波频域幅度的极大值,作为第三包络序列值;
根据所述第三包络序列值求取血压值。
5.根据权利要求4所述的无创血压测量的抗干扰方法,其特征在于,在对无干扰的脉搏波进行时域搜索分析并获取到第二包络序列值的周期信息后,求取对应的第二脉率值;
将所述第二脉率值与历史记录中的参考脉率值进行比较,若所述第二脉率值和所述参考脉率值之间的偏差值大于预设差值,则将所述第二脉率值更正为新的参考脉率值。
6.根据权利要求1所述的无创血压测量的抗干扰方法,其特征在于,所述判断所述脉搏波信号是否受到干扰包括:
判断同一压力平台期内的预设数量个脉搏波的幅度值、周期值、脉宽值和基线值是否均满足预设相似条件,若是,则判定为未受到干扰,否则判定为受到干扰;所述预设数量个为至少两个。
7.根据权利要求6所述的无创血压测量的抗干扰方法,其特征在于,所述预设相似条件为脉搏波两两之间的幅度值差、周期值差、脉宽值差和基线值差均分别小于对应的预设差值。
8.根据权利要求1所述的无创血压测量的抗干扰方法,其特征在于,所述延长所述脉搏波信号的采样时间具体为延长受到干扰的脉搏波信号的压力平台期的时间,直至频域分析中获得预置数量个采样点。
9.根据权利要求1所述的无创血压测量的抗干扰方法,其特征在于,所述获取脉搏数据中的压力波信号和脉搏波信号包括:
获取袖套中的原始压力信号,并对所述原始压力信号进行低通滤波获得压力波信号,对所述原始压力信号进行带通滤波获得脉搏波信号;
或,
通过AD采样器采集获取袖套中的原始压力信号,并对所述原始压力信号进行数字低通滤波分离获得压力波信号,对所述原始压力信号进行数字带通滤波分离获得脉搏波。
10.一种无创血压测量的抗干扰装置,其特征在于,包括:
获取模块,用于获取脉搏数据中的压力波信号和脉搏波信号;
分析模块,用于对所述脉搏波进行分析,并判断所述脉搏波信号是否受到干扰,若是,则延长所述脉搏波信号的采样时间并对所述脉搏波信号进行频域分析获得频域幅度的极大值,将所述极大值作为所述脉搏波的第一包络序列值,否则对所述脉搏波进行时域搜索分析,获得所述脉搏波的第二包络序列值;
血压值求取模块,用于根据所述第一包络序列值和/或所述第二包络序列值求取血压值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811120261.1A CN109124610B (zh) | 2018-09-26 | 2018-09-26 | 一种无创血压测量的抗干扰方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811120261.1A CN109124610B (zh) | 2018-09-26 | 2018-09-26 | 一种无创血压测量的抗干扰方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109124610A true CN109124610A (zh) | 2019-01-04 |
CN109124610B CN109124610B (zh) | 2021-07-30 |
Family
ID=64812298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811120261.1A Active CN109124610B (zh) | 2018-09-26 | 2018-09-26 | 一种无创血压测量的抗干扰方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109124610B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109770989A (zh) * | 2019-01-24 | 2019-05-21 | 中国人民解放军南部战区总医院 | 一种智能型气压式止血带及其控制方法 |
CN112914528A (zh) * | 2021-03-01 | 2021-06-08 | 出门问问信息科技有限公司 | 一种无袖带血压测量方法、装置及计算机可读介质 |
CN112914531A (zh) * | 2021-01-19 | 2021-06-08 | 深圳市杰纳瑞医疗仪器股份有限公司 | 确定血压包络波的方法、电子设备及存储介质 |
CN113349751A (zh) * | 2020-03-05 | 2021-09-07 | 深圳市理邦精密仪器股份有限公司 | 一种充放气控制方法、装置及医疗设备 |
CN113453615A (zh) * | 2019-02-18 | 2021-09-28 | 欧姆龙健康医疗事业株式会社 | 血压值分析辅助装置、血压值分析辅助系统、血压值分析辅助方法、以及程序 |
CN113811238A (zh) * | 2019-05-08 | 2021-12-17 | 深圳市长桑技术有限公司 | 血压计、血压计控制方法及检测有效的脉搏波的方法 |
CN114176546A (zh) * | 2021-08-07 | 2022-03-15 | 广东乐心医疗电子股份有限公司 | 血压测量方法、装置和电子设备 |
WO2023016378A1 (zh) * | 2021-08-07 | 2023-02-16 | 广东乐心医疗电子股份有限公司 | 血压测量方法、装置和电子设备 |
CN115721275A (zh) * | 2022-11-11 | 2023-03-03 | 深圳市捷美瑞科技有限公司 | 脉搏波数据采集和血压测量方法、装置、设备及介质 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101095612A (zh) * | 2006-06-28 | 2008-01-02 | 株式会社东芝 | 监测生物学信息的设备和方法 |
CN101334834A (zh) * | 2007-06-29 | 2008-12-31 | 北京交通大学 | 一种自底向上的注意信息提取方法 |
CN102499692A (zh) * | 2011-11-30 | 2012-06-20 | 沈阳工业大学 | 一种超声波步态检测装置与方法 |
CN102988036A (zh) * | 2012-12-26 | 2013-03-27 | 中国科学院自动化研究所 | 测量脉率的方法 |
CN103690152A (zh) * | 2014-01-06 | 2014-04-02 | 山东大学 | 一种基于脉搏解析的动脉弹性功能评估装置 |
CN103961082A (zh) * | 2013-01-29 | 2014-08-06 | 鸿富锦精密工业(深圳)有限公司 | 具有心率监测功能的便携式电子产品 |
CN106725415A (zh) * | 2016-11-15 | 2017-05-31 | 广州视源电子科技股份有限公司 | 电生理信号的处理方法和装置 |
CN107072561A (zh) * | 2014-10-31 | 2017-08-18 | 欧姆龙健康医疗事业株式会社 | 血压测量装置 |
CN107320089A (zh) * | 2017-06-27 | 2017-11-07 | 西南大学 | 自校准的人体血压测量方法 |
WO2018044059A1 (ko) * | 2016-08-31 | 2018-03-08 | (주)그린콤 | 피트니스 모니터링 시스템 |
CN107920763A (zh) * | 2015-07-16 | 2018-04-17 | 普瑞温迪克斯有限责任公司 | 处理生物数据 |
CN108272448A (zh) * | 2018-03-29 | 2018-07-13 | 合肥工业大学 | 一种非接触式的婴儿生理参数昼夜监护方法 |
CN108366744A (zh) * | 2015-12-28 | 2018-08-03 | 莱夫毕恩技术有限公司 | 用于检测生理参数的系统及方法 |
CN108464834A (zh) * | 2017-02-22 | 2018-08-31 | 清华大学深圳研究生院 | 一种基于心电信号的驾驶学员管理方法与系统 |
US20180249964A1 (en) * | 2017-03-02 | 2018-09-06 | Logos Care, Inc. | Deep learning algorithms for heartbeats detection |
-
2018
- 2018-09-26 CN CN201811120261.1A patent/CN109124610B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101095612A (zh) * | 2006-06-28 | 2008-01-02 | 株式会社东芝 | 监测生物学信息的设备和方法 |
CN101334834A (zh) * | 2007-06-29 | 2008-12-31 | 北京交通大学 | 一种自底向上的注意信息提取方法 |
CN102499692A (zh) * | 2011-11-30 | 2012-06-20 | 沈阳工业大学 | 一种超声波步态检测装置与方法 |
CN102988036A (zh) * | 2012-12-26 | 2013-03-27 | 中国科学院自动化研究所 | 测量脉率的方法 |
CN103961082A (zh) * | 2013-01-29 | 2014-08-06 | 鸿富锦精密工业(深圳)有限公司 | 具有心率监测功能的便携式电子产品 |
CN103690152A (zh) * | 2014-01-06 | 2014-04-02 | 山东大学 | 一种基于脉搏解析的动脉弹性功能评估装置 |
CN107072561A (zh) * | 2014-10-31 | 2017-08-18 | 欧姆龙健康医疗事业株式会社 | 血压测量装置 |
CN107920763A (zh) * | 2015-07-16 | 2018-04-17 | 普瑞温迪克斯有限责任公司 | 处理生物数据 |
CN108366744A (zh) * | 2015-12-28 | 2018-08-03 | 莱夫毕恩技术有限公司 | 用于检测生理参数的系统及方法 |
WO2018044059A1 (ko) * | 2016-08-31 | 2018-03-08 | (주)그린콤 | 피트니스 모니터링 시스템 |
CN106725415A (zh) * | 2016-11-15 | 2017-05-31 | 广州视源电子科技股份有限公司 | 电生理信号的处理方法和装置 |
CN108464834A (zh) * | 2017-02-22 | 2018-08-31 | 清华大学深圳研究生院 | 一种基于心电信号的驾驶学员管理方法与系统 |
US20180249964A1 (en) * | 2017-03-02 | 2018-09-06 | Logos Care, Inc. | Deep learning algorithms for heartbeats detection |
CN107320089A (zh) * | 2017-06-27 | 2017-11-07 | 西南大学 | 自校准的人体血压测量方法 |
CN108272448A (zh) * | 2018-03-29 | 2018-07-13 | 合肥工业大学 | 一种非接触式的婴儿生理参数昼夜监护方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109770989A (zh) * | 2019-01-24 | 2019-05-21 | 中国人民解放军南部战区总医院 | 一种智能型气压式止血带及其控制方法 |
CN113453615A (zh) * | 2019-02-18 | 2021-09-28 | 欧姆龙健康医疗事业株式会社 | 血压值分析辅助装置、血压值分析辅助系统、血压值分析辅助方法、以及程序 |
CN113811238A (zh) * | 2019-05-08 | 2021-12-17 | 深圳市长桑技术有限公司 | 血压计、血压计控制方法及检测有效的脉搏波的方法 |
CN113811238B (zh) * | 2019-05-08 | 2024-04-19 | 长桑数字技术(青岛)有限公司 | 血压计、血压计控制方法及检测有效的脉搏波的方法 |
CN113349751A (zh) * | 2020-03-05 | 2021-09-07 | 深圳市理邦精密仪器股份有限公司 | 一种充放气控制方法、装置及医疗设备 |
CN112914531A (zh) * | 2021-01-19 | 2021-06-08 | 深圳市杰纳瑞医疗仪器股份有限公司 | 确定血压包络波的方法、电子设备及存储介质 |
CN112914531B (zh) * | 2021-01-19 | 2022-09-30 | 深圳市杰纳瑞医疗仪器股份有限公司 | 确定血压包络波的方法、电子设备及存储介质 |
CN112914528B (zh) * | 2021-03-01 | 2023-06-27 | 出门问问信息科技有限公司 | 一种无袖带血压测量的模型生成方法、装置及计算机可读介质 |
CN112914528A (zh) * | 2021-03-01 | 2021-06-08 | 出门问问信息科技有限公司 | 一种无袖带血压测量方法、装置及计算机可读介质 |
CN114176546A (zh) * | 2021-08-07 | 2022-03-15 | 广东乐心医疗电子股份有限公司 | 血压测量方法、装置和电子设备 |
WO2023016378A1 (zh) * | 2021-08-07 | 2023-02-16 | 广东乐心医疗电子股份有限公司 | 血压测量方法、装置和电子设备 |
CN114176546B (zh) * | 2021-08-07 | 2022-12-09 | 广东乐心医疗电子股份有限公司 | 血压测量装置和电子设备 |
CN115721275A (zh) * | 2022-11-11 | 2023-03-03 | 深圳市捷美瑞科技有限公司 | 脉搏波数据采集和血压测量方法、装置、设备及介质 |
CN115721275B (zh) * | 2022-11-11 | 2023-11-24 | 深圳市捷美瑞科技有限公司 | 脉搏波数据采集和血压测量方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109124610B (zh) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109124610A (zh) | 一种无创血压测量的抗干扰方法和装置 | |
Kachuee et al. | Cuffless blood pressure estimation algorithms for continuous health-care monitoring | |
Esmaili et al. | Nonlinear cuffless blood pressure estimation of healthy subjects using pulse transit time and arrival time | |
Kurylyak et al. | A Neural Network-based method for continuous blood pressure estimation from a PPG signal | |
Chon et al. | Estimation of respiratory rate from photoplethysmogram data using time–frequency spectral estimation | |
CN108185996B (zh) | 动脉血管年龄估算模型构建方法和装置 | |
EP1740095B1 (en) | Non-invasive monitoring of respiratory rate | |
CN108498089A (zh) | 一种基于深度神经网络的无创连续血压测量方法 | |
CN103027690B (zh) | 一种基于自相关建模法的低灌注血氧饱和度测量方法 | |
CN112089405B (zh) | 一种脉搏波特征参数测量及显示装置 | |
CN106037694A (zh) | 一种基于脉搏波的连续血压测量装置 | |
CN112274127B (zh) | 一种基于单路脉搏波的无创连续血压检测方法、装置 | |
CN110327036A (zh) | 从穿戴式心电中提取呼吸信号和呼吸频率的方法 | |
Dastjerdi et al. | Non-invasive blood pressure estimation using phonocardiogram | |
CN107028603A (zh) | 使用脉搏触诊信号来检测人体内的糖尿病的装置和方法 | |
Chen et al. | Assessment of algorithms for oscillometric blood pressure measurement | |
CN101919704A (zh) | 一种心音信号定位、分段方法 | |
CN104068841B (zh) | 一种测量心脏收缩时间参数的测量方法及装置 | |
CN102197998B (zh) | 示波测量法中伪影的频谱的使用 | |
CN105852863B (zh) | 一种呼吸率测量方法及装置 | |
Qin et al. | Advances in cuffless continuous blood pressure monitoring technology based on PPG signals | |
Ghosh et al. | Introduction of boosting algorithms in continuous non-invasive cuff-less blood pressure estimation using pulse arrival time | |
Madhav et al. | Monitoring respiratory activity using PPG signals by order reduced-modified covariance AR technique | |
Fleming et al. | Non-invasive measurement of respiratory rate in children using the photoplethysmogram | |
Sidhu et al. | Performance analysis of oscillometric blood pressure estimation techniques in cardiac patients |
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 |