CN108272446B - 无创连续血压测量系统及其校准方法 - Google Patents
无创连续血压测量系统及其校准方法 Download PDFInfo
- Publication number
- CN108272446B CN108272446B CN201810092502.XA CN201810092502A CN108272446B CN 108272446 B CN108272446 B CN 108272446B CN 201810092502 A CN201810092502 A CN 201810092502A CN 108272446 B CN108272446 B CN 108272446B
- Authority
- CN
- China
- Prior art keywords
- calibration
- current
- blood pressure
- time
- formula
- 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
Links
- 230000036772 blood pressure Effects 0.000 title claims abstract description 110
- 238000000034 method Methods 0.000 title claims abstract description 68
- 230000035487 diastolic blood pressure Effects 0.000 claims description 63
- 230000035488 systolic blood pressure Effects 0.000 claims description 63
- 238000012937 correction Methods 0.000 claims description 14
- 239000000284 extract Substances 0.000 claims description 4
- 238000009530 blood pressure measurement Methods 0.000 abstract description 12
- 230000004962 physiological condition Effects 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 5
- 230000007547 defect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 230000000903 blocking effect Effects 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- 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/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
- A61B5/02125—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
-
- 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/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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Vascular Medicine (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Pulmonology (AREA)
- Ophthalmology & Optometry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明公开了一种无创连续血压测量系统及其校准方法。本发明在系统启动时先用普通的最小二乘法进行公式系数和常数项的初始化,使系统立即具有计算正确的无创连续血压值的能力,在之后的校准中再采用递归最小二乘法。本发明一方面克服了直接应用递归最小二乘法将会导致的系统启动过程中计算公式系数和常数项的收敛消耗较长时间的技术障碍,使系统自启动开始立即具有计算正确的无创连续血压值的能力;另一方面克服了单纯使用最小二乘法所得的血压计算公式不能反映使用者的一般生理情况,而只能反映其较短时间内的生理情况,且受每次间歇血压测量的偶然误差影响大,连续无创血压结果的准确度难以达到临床要求的缺陷。
Description
技术领域
本发明涉及一种无创连续血压测量系统,以及使用测量系统监测血压时的系统校准方法,属于血压监测技术领域。
背景技术
血压是人体的重要生理参数之一,在高血压等各种心血管疾病的诊疗中起着重要作用,也是日常健康监护的重要项目。临床上最常用的血压值是收缩压和舒张压。目前常用的无创血压监测是基于示波法或柯氏音法的间歇血压监测。间歇血压监测除了不能提供连续的血压数据外,还存在着袖带频繁充放气给使用者带来的不适的问题。
现有的无创连续血压测量方法,主要有扁平张力法、恒定容积法和脉搏波传导时间法。扁平张力法对测量位点选取和传感器放置的要求严格,若测量位点选取或传感器放置不当则会造成大的测量误差;恒定容积法需要一套精密的伺服装置,且伺服装置的佩戴给使用者的日常生活行为带来较大不便。脉搏波传导时间法对测量位点选取和传感器放置的要求较宽松,故使用方便,且给使用者日常生活造成的不便少,因而越来越受到重视。
脉搏波传导时间法需要建立血压和包括脉搏波传导时间在内的变量组之间的函数关系。医学研究表明这个函数关系是因人而异且时变的,故需要进行系统校准,即利用间歇血压测量对计算公式的系数和常数项进行修正。
现有的校准方法,只利用最近的两次或若干次测量的间歇血压值,通过最小二乘法进行系数和常数项的校准,而没有利用较早测量的间歇血压值,这就导致所得血压计算公式不能反映使用者的一般生理情况,而只能反映其较短时间内的生理情况,且受每次间歇血压测量的偶然误差影响大,故求得的连续无创血压结果的准确度难以达到临床要求。
递归最小二乘法(RLS法),作为一种在最优滤波器中应用的自适应方法,在数字信号处理、系统辨识等方面有广泛的应用。它带有遗忘因子,能令较新的数据在自适应过程中以更大的权重,而较早的数据的作用则逐步除去,从而使滤波器可以工作在平稳状态下。其对非平稳信号的适应性好,收敛速度快,精度高。而脉搏波传导时间法的校准与最优滤波器的自适应分属于两个不同的领域,目前未见将递归最小二乘法用于脉搏波传导时间法的计算公式校准中。
综上所述,脉搏波传导时间法相比其他的无创连续血压测量方法存在较大优势,但基于此方法的测量系统必须妥善解决系统校准问题,否则测量准确度难以达到临床要求。这一问题亟待解决。
发明内容
本发明的目的是提供一种无创连续血压测量系统及其校准方法,以获得准确的无创连续收缩压和舒张压。
本发明的发明构思如下:递归最小二乘法(RLS法)作为一种在最优滤波器中应用的自适应方法,具有对非平稳信号的适应性好、收敛速度快、精度高等优点。而脉搏波传导时间法的校准与最优滤波器的自适应分属于两个不同的领域,若要在脉搏波传导时间法的计算公式校准中直接采用这种方法,尤其是从系统启动时就采用递归最小二乘法,则计算公式系数和常数项的收敛将会消耗较长的时间,系统无法自始即具备计算正确的无创连续血压值的能力。为解决这一问题,本发明在系统启动时先用普通的最小二乘法进行公式系数和常数项的初始化,使系统立即具有计算正确的无创连续血压值的能力,在之后的校准中再采用递归最小二乘法。
具体地说,为实现上述发明目的,本发明所采用的技术方案如下:
本发明无创连续血压测量系统包括连续生理信号采集单元(1)、间歇校准血压信号采集单元(2)、连续血压计算与系统校准单元(3);
所述连续信号采集单元(1)包括连续心电信号采集单元(11)和连续脉搏波信号采集单元(12);所述连续心电信号采集单元(11)连续采集心电信号,并提取每搏心电信号的特征点的时刻发送至连续血压计算与系统校准单元(3);所述连续脉搏波信号采集单元(12)连续采集脉搏波信号,并提取每搏脉搏波信号的特征点的时刻发送至连续血压计算与系统校准单元(3);
所述间歇校准血压信号采集单元(2)在系统启动时和连续血压计算与系统校准单元(3)向其发出校准开始命令时,测量第一舒张压和第一收缩压并发送至连续血压计算与系统校准单元(3);
所述连续血压计算与系统校准单元(3)在系统启动时,根据所接收的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到脉搏波传导时间和心率,再根据脉搏波传导时间和心率以及对应的第一舒张压、第一收缩压,利用最小二乘法对公式(1)和公式(2)的系数和常数项进行初始化,使系统达到能连续输出第二收缩压和第二舒张压的状态,并设置第1个校准时间点;所述连续血压计算与系统校准单元(3)在不执行校准时,根据当前获得的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到当前脉搏波传导时间和当前心率,并根据当前脉搏波传导时间和当前心率分别按照当前公式(1)、公式(2)进行计算,得到当前第二收缩压和当前第二舒张压;当所述连续血压计算与系统校准单元(3)执行校准时,根据当前第一收缩压、当前第一舒张压、当前第二收缩压、当前第二舒张压、当前脉搏波传导时间、当前心率以及当前公式(1)和当前公式(2)的各系数和常数项,利用递归最小二乘法,计算当前公式(1)和当前公式(2)的各系数和常数项的修正值,再将各系数和常数项分别与其修正值相加后得到的值作为更新后的公式(1)和公式(2)的系数和常数项,随后利用更新后的公式(1)和公式(2)计算当前第二收缩压和当前第二舒张压,并设置下一校准时间点;
式中,k≥0,h≥1,m≥1,n≥1;当k=0时,SBPk,h表示第1次校准前的第h个第二收缩压,DBPk,h表示第1次校准前的第h个第二舒张压,ai,k和bj,k分别表示第1次校准前的系数,a0,k和b0,k分别表示第1次校准前的常数项,PTTk,h表示第1次校准前的第h个脉搏波传导时间,HRk,h表示第1次校准前的第h个心率;当k≥1时,SBPk,h表示自第k次校准到第k+1次校准前得到的第h个第二收缩压,DBPk,h表示自第k次校准到第k+1次校准前得到的第h个第二舒张压,ai,k和bj,k分别表示自第k次校准到第k+1次校准前的系数,a0,k和b0,k分别表示自第k次校准到第k+1次校准前的常数项,PTTk,h表示自第k次校准到第k+1次校准前的第h个脉搏波传导时间,HRk,h表示自第k次校准到第k+1次校准前的第h个心率;fi(PTTk,h,HRk,h)是关于PTTk,h和HRk,h两个变量中至少一个的函数,该函数在系统启动后不再发生变化;gj(PTTk,h,HRk,h)是关于PTTk,h和HRk,h两个变量中至少一个的函数,该函数在系统启动后不再发生变化。
本发明无创连续血压测量系统的校准方法包括如下步骤:
1)系统启动时,先由间歇校准血压信号采集单元(2)采集第一收缩压和第一舒张压发送给连续血压计算与系统校准单元(3);连续血压计算与系统校准单元(3)根据来自所述连续信号采集单元(1)的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到脉搏波传导时间和心率,并根据所接收的第一收缩压和第一舒张压以及对应的脉搏波传导时间和心率,利用最小二乘法,对公式(1)和公式(2)的系数和常数项进行初始化,使系统达到能连续输出第二收缩压和第二舒张压的状态,并设置第1个校准时间点;
2)当校准时间点到来时,连续血压计算与系统校准单元(3)向间歇校准血压信号采集单元(2)发送校准开始命令,间歇校准血压信号采集单元(2)采集当前第一收缩压和当前第一舒张压并发送到连续血压计算与系统校准单元(3);连续血压计算与系统校准单元(3)根据当前第一收缩压、当前第一舒张压、当前第二收缩压、当前第二舒张压、当前脉搏波传导时间、当前心率以及当前公式(1)和当前公式(2)的各系数和常数项,利用递归最小二乘法,计算当前公式(1)和当前公式(2)的各系数和常数项的修正值,再将各系数和常数项分别与其修正值相加后得到的值作为更新后的公式(1)和公式(2)的系数和常数项;
3)连续血压计算与系统校准单元(3)根据当前公式(1)和当前公式(2)计算当前第二收缩压和当前第二舒张压,且设置下一个校准时间点,并在下一个校准时间点到来时返回执行步骤2)。
与现有技术相比,本发明的有益效果是:
1)本发明无创连续血压测量系统,在由连续生理信号采集单元、间歇校准血压信号采集单元和连续血压计算与系统校准单元组成的硬件平台的基础上,实现了基于脉搏波传导时间PTT和心率HR两种参数进行的无创连续血压计算。
2)本发明无创连续血压测量系统在采用最小二乘法进行系统初始化的基础上,采用递归最小二乘(RLS)法进行系统校准,不断在原有公式系数和常数项的基础上加修正值,这就在利用新的间歇血压值的同时,通过原有公式系数和常数项以较低的权重间接地利用了较早测量的间歇血压值,进而在公式中同时反映使用者的一般的和当前的生理情况,降低每次间歇血压测量的偶然误差对公式系数和常数项造成的影响,提高了连续无创血压测量的准确度。由此,本发明不仅将本应用于最优滤波器自适应的递归最小二乘法应用到脉搏波传导时间法的计算公式校准中,而且通过在系统启动时先用普通的最小二乘法进行公式系数和常数项的初始化,使系统立即具有计算正确的无创连续血压值的能力,在之后的校准中再采用递归最小二乘法,从而一方面克服了直接应用递归最小二乘法将会导致的系统启动过程中计算公式系数和常数项的收敛消耗较长时间的技术障碍,使系统自启动开始立即具有计算正确的无创连续血压值的能力;另一方面克服了单纯使用最小二乘法所得的血压计算公式不能反映使用者的一般生理情况,而只能反映其较短时间内的生理情况,且受每次间歇血压测量的偶然误差影响大,连续无创血压结果的准确度难以达到临床要求的缺陷,获得了预料不到的技术效果。
附图说明
图1是本发明的一种无创连续血压测量系统的结构示意图;
图2是本发明的连续心电信号采集单元的一种实施方式的结构示意图;
图3是本发明的连续脉搏波信号采集单元的一种实施方式的结构示意图;
图4是本发明的间歇校准血压信号采集单元的一种实施方式的结构示意图;
图5是本发明的连续血压计算与系统校准单元的一种实施方式的结构示意图;
图6是本发明无创连续血压测量系统佩戴于人体时的一种实施方式的示意图。
具体实施方式
下面结合附图与具体实施方式对本发明做进一步详细描述。
如图1所示,本发明无创连续血压测量系统包括连续生理信号采集单元1、间歇校准血压信号采集单元2、以及连续血压计算与系统校准单元3。连续信号采集单元1包括连续心电信号采集单元11和连续脉搏波信号采集单元12。连续心电信号采集单元11连续采集心电信号,进一步求出其R波波峰点所在的时刻并发送至连续血压计算与系统校准单元3;连续脉搏波信号采集单元12连续采集脉搏波信号,进一步求出其主波波峰点所在的时刻并发送至连续血压计算与系统校准单元3。间歇校准血压信号采集单元2在系统启动时以及在连续血压计算与系统校准单元3给出校准开始命令时,测量第一舒张压和第一收缩压并发送至连续血压计算与系统校准单元3。连续血压计算与系统校准单元3在系统启动时,接收间歇校准血压信号采集单元2发来的第一舒张压和第一收缩压,并根据连续信号采集单元1发来的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到脉搏波传导时间和心率,再根据接收到的第一舒张压、第一收缩压、以及与第一舒张压、第一收缩压对应的脉搏波传导时间和心率,利用最小二乘法对公式(1)和公式(2)的系数和常数项进行初始化,使系统达到能连续输出第二收缩压和第二舒张压的状态,并设置第1个校准时间点。此后,当连续血压计算与系统校准单元3不执行校准时,根据当前获得的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到当前脉搏波传导时间和当前心率,并根据当前脉搏波传导时间和当前心率分别按照当前公式(1)、公式(2)进行计算,得到当前第二收缩压和当前第二舒张压。当连续血压计算与系统校准单元3执行校准时,根据当前第一收缩压、当前第一舒张压、当前第二收缩压、当前第二舒张压、当前脉搏波传导时间、当前心率以及当前公式(1)和当前公式(2)的各系数和常数项,利用递归最小二乘法,计算当前公式(1)和当前公式(2)的各系数和常数项的修正值,再将各系数和常数项分别与其修正值相加后得到的值作为更新后的公式(1)和公式(2)的系数和常数项,随后利用更新后的公式(1)和公式(2)计算当前第二收缩压和第二舒张压,并设置下一校准时间点。
以下以具体的实施例进一步详细说明本发明。
如图2所示,连续心电信号采集单元11可采用双电极心电采集电路。心电电极1和心电电极2采集到的心电信号通过前置放大器进行差分放大,再通过滤波器进行带通滤波和工频陷波,最后通过主放大器放大后经第一模拟/数字转换器(ADC)变换为数字信号送入第一微处理器,由第一微处理器求出心电信号的R波波峰点所在的时刻,并通过第一无线发射电路将该时刻发送给连续血压计算与系统校准单元3。前置放大器、滤波器和主放大器可采用由通用运放和阻容元件实现的普通放大电路和滤波电路。第一模拟/数字转换器(ADC)和第一微处理器可采用将两者集成的单片机,如TI公司的MSP430系列16位低功耗单片机。第一无线发射电路可采用基于现有的Zigbee、蓝牙等协议的无线发射电路。
如图3所示,连续脉搏波信号采集单元12可采用指端式光电脉搏波采集电路。该电路中的指夹式脉搏波传感器采集到的脉搏波信号通过隔直电路除去直流成分,再通过放大器放大,进而通过低通滤波器滤去高频成分后经第二模拟/数字转换器(ADC)变换为数字信号送入第二微处理器,由第二微处理器求出脉搏波的主波波峰点所在的时刻,并通过第二无线发射电路将该时刻发送给连续血压计算与系统校准单元3。隔直电路、放大器和低通滤波器可采用由通用运放和阻容元件实现的普通隔直电路、放大电路和滤波电路。第二模拟/数字转换器(ADC)和第二微处理器可采用将这两者集成的单片机,如TI公司的MSP430系列16位低功耗单片机。第二无线发射电路可采用基于现有的Zigbee、蓝牙等协议的无线发射电路。
如图4所示,作为本发明的优选实施方式,间歇校准血压信号采集单元2可采用通用的袖带式血压计结构,包括袖带、充/放气泵、压力传感器、第三模拟/数字转换器(ADC)和第三微处理器。压力传感器测量袖带内的压力,发送给第三模拟/数字转换器(ADC)进行数字化。第三微处理器具备以下功能:接收连续血压计算与系统校准单元3的校准开始命令;控制充/放气泵对袖带进行充/放气;在袖带放气过程中接收由第三模拟/数字转换器(ADC)数字化的压力信号,据此计算出第一收缩压和第一舒张压并发送给连续血压计算与系统校准单元3。袖带、充/放气泵、压力传感器可采用市售常规产品。第三模拟/数字转换器(ADC)和第三微处理器可采用将两者集成的单片机,如TI公司的MSP430系列16位低功耗单片机。
如图5所示,连续血压计算与系统校准单元3包括第四微处理器、无线收发电路和输出模块。第四微处理器在系统启动时,通过无线收发电路接收间歇校准血压信号采集单元2发来的第一舒张压和第一收缩压;并根据连续信号采集单元1发来的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到脉搏波传导时间和心率,再根据第一舒张压、第一收缩压、以及与第一舒张压、第一收缩压对应的脉搏波传导时间和心率,利用最小二乘法对公式(1)和公式(2)的系数和常数项进行初始化,使系统达到能连续输出第二收缩压和第二舒张压的状态,并设置第1个校准时间点;此后,当不执行校准时,根据通过无线收发电路获得的当前心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到当前脉搏波传导时间和心率,并根据当前脉搏波传导时间和心率分别按照当前公式(1)、公式(2)进行计算,得到当前第二收缩压和第二舒张压;当执行校准时,根据当前第一收缩压、当前第一舒张压、当前第二收缩压、当前第二舒张压、当前脉搏波传导时间、当前心率以及当前公式(1)和当前公式(2)的各系数和常数项,利用递归最小二乘法,计算当前公式(1)和当前公式(2)的各系数和常数项的修正值,再将各系数和常数项分别与其修正值相加后得到的值作为更新后的公式(1)和公式(2)的系数和常数项,随后利用更新后的公式(1)和公式(2)计算当前第二收缩压和第二舒张压,并设置下一校准时间点。得到的第二收缩压和第二舒张压由输出模块输出。第四微处理器可采用通用单片机,如TI公司的MSP430系列16位低功耗单片机。无线收发电路可采用基于现有的Zigbee、蓝牙等协议的无线收发电路。输出模块可采用串口、无线发射电路、液晶显示屏等。
根据以下公式(1)和公式(2):
上式中,k≥0,h≥1,m≥1,n≥1;当k=0时,SBPk,h表示第1次校准前的第h个第二收缩压,DBPk,h表示第1次校准前的第h个第二舒张压,a1,k、a2,k和b1,k分别表示第1次校准前的系数,a0,k和b0,k分别表示第1次校准前的常数项,PTTk,h表示第1次校准前的第h个脉搏波传导时间,HRk,h表示第1次校准前的第h个心率;当k≥1时,SBPk,h表示自第k次校准时(含)到第k+1次校准前得到的第h个第二收缩压,DBPk,h表示自第k次校准时(含)到第k+1次校准前得到的第h个第二舒张压,a1,k、a2,k和b1,k分别表示自第k次校准时(含)到第k+1次校准前的系数,a0,k和b0,k分别表示自第k次校准时(含)到第k+1次校准前的常数项,PTTk,h表示自第k次校准时(含)到第k+1次校准前的第h个脉搏波传导时间,HRk,h表示自第k次校准时(含)到第k+1次校准前的第h个心率。
fi(PTTk,h,HRk,h)是关于PTTk,h和HRk,h两个变量中至少一个的函数,该函数在系统启动后不再发生变化;gj(PTTk,h,HRk,h)是关于PTTk,h和HRk,h两个变量中至少一个的函数,该函数在系统启动后不再发生变化。fi(PTTk,h,HRk,h)可以取PTTk,h、ln(PTTk,h)、HRk,h、(PTTk,h×HRk,h)等具体形式;gj(PTTk,h,HRk,h)可以取的具体形式与fi(PTTk,h,HRk,h)相同。
优选地,本实施例中所用的血压计算公式取m=2,n=1,f1(PTTk,h,HRk,h)取PTTk,h,f2(PTTk,h,HRk,h)取HRk,h;g1(PTTk,h,HRk,h)取PTTk,h,即可得到以下式(3)和式(4):
SBPk,h=a1,k·PTTk,h+a2,k·HRk,h+a0,k (3)
DBPk,h=b1,k·PTTk,h+b0,k (4)
公式(3)和(4)的系数和常数项a1,k、a2,k、a0,k、b1,k和b0,k在系统初始化时采用最小二乘法获得初值a1,0、a2,0、a0,0、b1,0和b0,0,并在每次校准时进行调整。
由上可见,由于公式(1)和公式(2)所采用的公式可取多种具体形式,使系统灵活性更强,可以更好地适应不同个体的需求。
本发明无创连续血压测量系统的工作流程如下:
步骤1):系统开始工作时,首先进行初始化。间歇校准血压信号采集单元2进行间歇血压测量,测量次数取m+1和n+1中较大者。每测量一次,就将测得的第一收缩压和第一舒张压发送给连续血压计算与系统校准单元3。在间歇血压测量全部完成后,连续血压计算与系统校准单元3根据所获得的所有第一收缩压和所有第一舒张压以及与之对应的的脉搏波传导时间和心率,利用最小二乘法,求出初始公式系数和常数项a1,0、a2,0、a0,0、b1,0、b0,0,并设定2min后为第1个校准时间点。
此后,连续血压计算与系统校准单元3接收连续信号采集单元1发送的心电信号的R波波峰点所在的时刻以及脉搏波信号主波波峰点所在的时刻。连续血压计算与系统校准单元3即可依据式(3)和式(4)所示的血压计算模型,根据由心电信号的R波波峰点所在的时刻以及脉搏波信号主波波峰点所在的时刻得到的脉搏波传导时间和心率,计算第二收缩压和第二舒张压并输出。
步骤2):当第k个校准时间点到来时(此处k≥1),连续血压计算与系统校准单元3发送校准开始命令给间歇校准血压信号采集单元2,间歇校准血压信号采集单元2进行1次间歇血压测量;连续血压计算与系统校准单元3接收间歇校准血压信号采集单元2测得的当前第一收缩压和当前第一舒张压,并根据当前第一收缩压、当前第一舒张压、当前第二收缩压、当前第二舒张压、当前脉搏波传导时间、当前心率,以及当前公式系数和常数项,利用递归最小二乘(RLS)法,计算公式(3)和公式(4)系数和常数项的修正值,并分别加在当前两公式的公式系数和常数项上,对两公式系数和常数项进行修正,这样就进行了系统校准。
步骤3):连续血压计算与系统校准单元3根据当前公式(3)和当前公式(4)计算当前第二收缩压SBPk,1和当前第二舒张压DBPk,1,且设置第k+1个校准时间点(此处k≥1),并在第k+1个校准时间点到来时返回执行步骤2)。优选地,第k+1个校准时间点与第k个校准时间点的时间间隔可由以下公式(5)算出:
公式(5)中,ΔTk+1为第k+1个校准时间点与第k准时间点的时间间隔;|ΔSBPk|为第二收缩压SBPk,1与当前第一收缩压的误差的绝对值;|ΔDBPk|为第二舒张压DBPk,1与当前第一舒张压的误差的绝对值。
如图6所示,本发明无创连续血压测量系统使用时,可将连续心电信号采集单元11通过胸带固定于胸前;将连续脉搏波信号采集单元12夹于左手任一手指的指端;连续血压计算与系统校准单元3与间歇校准血压信号采集单元2安装在一起,佩戴于右上臂。这样的左右分工可以避免对使用者的右手造成干扰,同时也可避免利用袖带进行校准血压的测量时,对右臂血流的暂时性阻滞导致测得的脉搏波传导时间不正确。
Claims (1)
1.一种无创连续血压测量系统,其特征在于,包括连续生理信号采集单元(1)、间歇校准血压信号采集单元(2)、连续血压计算与系统校准单元(3);
所述连续信号采集单元(1)包括连续心电信号采集单元(11)和连续脉搏波信号采集单元(12);所述连续心电信号采集单元(11)连续采集心电信号,并提取每搏心电信号的特征点的时刻发送至连续血压计算与系统校准单元(3);所述连续脉搏波信号采集单元(12)连续采集脉搏波信号,并提取每搏脉搏波信号的特征点的时刻发送至连续血压计算与系统校准单元(3);
所述间歇校准血压信号采集单元(2)在系统启动时和连续血压计算与系统校准单元(3)向其发出校准开始命令时,测量第一舒张压和第一收缩压并发送至连续血压计算与系统校准单元(3);
所述连续血压计算与系统校准单元(3)在系统启动时,根据所接收的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到脉搏波传导时间和心率,再根据脉搏波传导时间和心率以及对应的第一舒张压、第一收缩压,利用最小二乘法对公式(1)和公式(2)的系数和常数项进行初始化,使系统达到能连续输出第二收缩压和第二舒张压的状态,并设置第1个校准时间点;所述连续血压计算与系统校准单元(3)在不执行校准时,根据当前获得的心电信号的特征点的时刻以及脉搏波信号的特征点的时刻,得到当前脉搏波传导时间和当前心率,并根据当前脉搏波传导时间和当前心率分别按照当前公式(1)、公式(2)进行计算,得到当前第二收缩压和当前第二舒张压;当所述连续血压计算与系统校准单元(3)执行校准时,根据当前第一收缩压、当前第一舒张压、当前第二收缩压、当前第二舒张压、当前脉搏波传导时间、当前心率以及当前公式(1)和当前公式(2)的各系数和常数项,利用递归最小二乘法,计算当前公式(1)和当前公式(2)的各系数和常数项的修正值,再将各系数和常数项分别与其修正值相加后得到的值作为更新后的公式(1)和公式(2)的系数和常数项,随后利用更新后的公式(1)和公式(2)计算当前第二收缩压和当前第二舒张压,并设置下一校准时间点;
式中,k≥0,h≥1,m≥1,n≥1;当k=0时,SBPk,h表示第1次校准前的第h个第二收缩压,DBPk,h表示第1次校准前的第h个第二舒张压,ai,k和bj,k分别表示第1次校准前的系数,a0,k和b0,k分别表示第1次校准前的常数项,PTTk,h表示第1次校准前的第h个脉搏波传导时间,HRk,h表示第1次校准前的第h个心率;当k≥1时,SBPk,h表示自第k次校准到第k+1次校准前得到的第h个第二收缩压,DBPk,h表示自第k次校准到第k+1次校准前得到的第h个第二舒张压,ai,k和bj,k分别表示自第k次校准到第k+1次校准前的系数,a0,k和b0,k分别表示自第k次校准到第k+1次校准前的常数项,PTTk,h表示自第k次校准到第k+1次校准前的第h个脉搏波传导时间,HRk,h表示自第k次校准到第k+1次校准前的第h个心率;fi(PTTk,h,HRk,h)是关于PTTk,h和HRk,h两个变量中至少一个的函数,该函数在系统启动后不再发生变化;gj(PTTk,h,HRk,h)是关于PTTk,h和HRk,h两个变量中至少一个的函数,该函数在系统启动后不再发生变化;
下一校准时间点与当前校准时间点的时间间隔由以下公式(5)计算:
公式(5)中,k≥1,ΔTk+1为下一校准时间点与当前校准时间点的时间间隔;|ΔSBPk|为当前第二收缩压与当前第一收缩压的误差的绝对值;|ΔDBPk|为当前第二舒张压与当前第一舒张压的误差的绝对值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810092502.XA CN108272446B (zh) | 2018-01-30 | 2018-01-30 | 无创连续血压测量系统及其校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810092502.XA CN108272446B (zh) | 2018-01-30 | 2018-01-30 | 无创连续血压测量系统及其校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108272446A CN108272446A (zh) | 2018-07-13 |
CN108272446B true CN108272446B (zh) | 2021-03-26 |
Family
ID=62807023
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810092502.XA Active CN108272446B (zh) | 2018-01-30 | 2018-01-30 | 无创连续血压测量系统及其校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108272446B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108968937A (zh) * | 2018-07-23 | 2018-12-11 | 山东美泽赫生理技术有限公司 | 一种血压监测装置及监测方法 |
CN109222942A (zh) * | 2018-11-09 | 2019-01-18 | 中科数字健康科学研究院(南京)有限公司 | 一种无创连续血压测量方法和设备 |
CN110292370B (zh) * | 2019-07-03 | 2020-12-15 | 浙江大学 | 一种基于脉搏波传导时间的胸口无创血压检测方法 |
CN113171069B (zh) * | 2021-03-05 | 2024-01-26 | 上海立阖泰医疗科技有限公司 | 一种嵌入式图形血压测量系统 |
CN113180621B (zh) * | 2021-03-05 | 2024-01-26 | 上海立阖泰医疗科技有限公司 | 基于freeRTOS的连续无创血压测量系统 |
CN113925477B (zh) * | 2021-07-08 | 2023-10-27 | 哈尔滨铂云医疗器械有限公司 | 电子血压计自动校准系统 |
CN113940645B (zh) * | 2021-10-21 | 2024-01-23 | 上海海尔医疗科技有限公司 | 基于最小二乘法的血压测量系统及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010037056A1 (en) * | 2008-09-26 | 2010-04-01 | Qualcomm Incorporated | Method and apparatus for non-invasive cuff-less blood pressure estimation using pulse arrival time and heart rate with adaptive calibration |
TW201021763A (en) * | 2008-09-29 | 2010-06-16 | Qualcomm Inc | Method and apparatus for non-invasive cuff-less blood pressure estimation using pulse arrival time and heart rate with adaptive calibration |
CN104739395A (zh) * | 2015-03-25 | 2015-07-01 | 华中科技大学 | 一种基于脉搏波的人体血压预测方法 |
CN105455798A (zh) * | 2015-10-19 | 2016-04-06 | 东南大学 | 基于Android手机端的连续血压测量系统和校准测量方法 |
CN106343976A (zh) * | 2016-09-14 | 2017-01-25 | 京东方科技集团股份有限公司 | 建立血压模型的方法和装置以及确定血压的方法和装置 |
WO2017028011A1 (zh) * | 2015-08-14 | 2017-02-23 | 华为技术有限公司 | 一种血压测量数据的处理方法及装置 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4669485A (en) * | 1984-02-17 | 1987-06-02 | Cortronic Corporation | Apparatus and method for continuous non-invasive cardiovascular monitoring |
JP3581975B2 (ja) * | 1997-01-06 | 2004-10-27 | 日本光電工業株式会社 | 血圧監視装置 |
CN1698536A (zh) * | 2004-05-20 | 2005-11-23 | 香港中文大学 | 采用自动补偿的无袖带式连续血压测量方法 |
MX2007014997A (es) * | 2005-06-01 | 2008-04-22 | Innosurance Inc | Recoleccion y analisis de datos de operacion de vehiculo de motor. |
US20090093222A1 (en) * | 2007-10-03 | 2009-04-09 | Qualcomm Incorporated | Calibration and beamforming in a wireless communication system |
CN101286059B (zh) * | 2008-05-29 | 2010-06-02 | 上海交通大学 | 在线测量数据获取、校正及过失误差的序列补偿方法 |
JP5884256B2 (ja) * | 2010-05-19 | 2016-03-15 | セイコーエプソン株式会社 | 血圧測定装置及び血圧測定方法 |
EP2551650B1 (en) * | 2011-07-27 | 2019-09-25 | Endress+Hauser Consult AG | Calibration method |
US8951204B2 (en) * | 2012-05-04 | 2015-02-10 | The Guy P. Curtis And Frances L. Curtis Trust | Method for using a pulse oximetry signal to monitor blood pressure |
EP3383256B1 (en) * | 2015-12-02 | 2024-04-10 | Spry Health, Inc. | Systems and methods for detecting photoplethysmographic device usage |
WO2018095083A1 (zh) * | 2016-11-22 | 2018-05-31 | 浙江脉联医疗设备有限公司 | 脉搏波传播时间的校正方法 |
CN106618504A (zh) * | 2017-01-19 | 2017-05-10 | 陈妍 | 一种具备自适应校准的高准确度连续无创血压测量装置 |
CN107137067A (zh) * | 2017-05-04 | 2017-09-08 | 深圳市松恩电子科技有限公司 | 一种脉搏波连续血压监测装置 |
-
2018
- 2018-01-30 CN CN201810092502.XA patent/CN108272446B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010037056A1 (en) * | 2008-09-26 | 2010-04-01 | Qualcomm Incorporated | Method and apparatus for non-invasive cuff-less blood pressure estimation using pulse arrival time and heart rate with adaptive calibration |
TW201021763A (en) * | 2008-09-29 | 2010-06-16 | Qualcomm Inc | Method and apparatus for non-invasive cuff-less blood pressure estimation using pulse arrival time and heart rate with adaptive calibration |
CN104739395A (zh) * | 2015-03-25 | 2015-07-01 | 华中科技大学 | 一种基于脉搏波的人体血压预测方法 |
WO2017028011A1 (zh) * | 2015-08-14 | 2017-02-23 | 华为技术有限公司 | 一种血压测量数据的处理方法及装置 |
CN105455798A (zh) * | 2015-10-19 | 2016-04-06 | 东南大学 | 基于Android手机端的连续血压测量系统和校准测量方法 |
CN106343976A (zh) * | 2016-09-14 | 2017-01-25 | 京东方科技集团股份有限公司 | 建立血压模型的方法和装置以及确定血压的方法和装置 |
Non-Patent Citations (1)
Title |
---|
Noninvasive cuffless blood pressure estimation using pulse transit time and Hilbert-Huang transform;Younhee Choi,Qiao Zhang,Seokbum Ko;《Computers and Electrical Engineering》;20121108(第39期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108272446A (zh) | 2018-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108272446B (zh) | 无创连续血压测量系统及其校准方法 | |
CN107865647B (zh) | 血压检测装置以及血压检测装置的校正方法 | |
AU2016300616B2 (en) | Blood pressure calculation method based on pulse reflected wave transit time, and blood pressure meter | |
US20200245877A1 (en) | Continuous Blood Pressure Measurement | |
US20020038090A1 (en) | Central-artery-pressure-waveform estimating apparatus | |
JP2018508273A (ja) | 被験体の平均動脈圧を導き出すための方法および装置 | |
US20120157791A1 (en) | Adaptive time domain filtering for improved blood pressure estimation | |
EP2501278B1 (en) | Apparatus for sensing blood flow and hemodynamic parameters | |
WO2003007816A1 (en) | Arterial line emulator | |
Dastjerdi et al. | Non-invasive blood pressure estimation using phonocardiogram | |
US20130060152A1 (en) | Apparatus and method for continuous oscillometric blood pressure measurement | |
US20140081159A1 (en) | Non-invasive continuous blood pressure monitoring system and method | |
WO1996039074A1 (en) | Improved calibration for blood pressure pulses | |
EP3648664A1 (en) | Self-calibrating systems and methods for blood pressure wave form analysis and diagnostic support | |
CN112890790B (zh) | 一种穿戴式无创血压动态跟踪监测方法 | |
US20200253564A1 (en) | Continuous Blood Pressure Measurement | |
TWI832029B (zh) | 非施壓連續血壓量測裝置及方法 | |
CN108926334A (zh) | 基于脉搏波的血压获取方法及其系统和装置 | |
Zhang et al. | A LabVIEW based measure system for pulse wave transit time | |
US20220218217A1 (en) | Blood pressure measurement system and blood pressure measurement method using the same | |
CN210095711U (zh) | 一种无创连续血压测量设备 | |
US20150032012A1 (en) | Non-invasive Blood Pressure Measurement System and Methods of Use | |
Al-Masri et al. | Designing a biosignal extraction device for measuring a blood pressure | |
Victor et al. | Blood presure monitoring using programmable system on chip | |
AU5806296A (en) | Improved calibration for blood pressure pulses |
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 |