CN101732040B - 一种无创多路脉搏波检测方法 - Google Patents

一种无创多路脉搏波检测方法 Download PDF

Info

Publication number
CN101732040B
CN101732040B CN 200910243060 CN200910243060A CN101732040B CN 101732040 B CN101732040 B CN 101732040B CN 200910243060 CN200910243060 CN 200910243060 CN 200910243060 A CN200910243060 A CN 200910243060A CN 101732040 B CN101732040 B CN 101732040B
Authority
CN
China
Prior art keywords
pulse wave
pulse
point
pressure
signal
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.)
Expired - Fee Related
Application number
CN 200910243060
Other languages
English (en)
Other versions
CN101732040A (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.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
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 Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN 200910243060 priority Critical patent/CN101732040B/zh
Publication of CN101732040A publication Critical patent/CN101732040A/zh
Application granted granted Critical
Publication of CN101732040B publication Critical patent/CN101732040B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种无创多路脉搏波检测装置,包括至少两个传感器,每个传感器上均设置有固定装置,并且各个传感器均与信号处理单元连接,传感器测得的信号通过信号处理单元进行处理、储存。本发明对人体不同部位的脉搏波采集点采用多路同步采集脉搏波和波速、波形处理技术,可针对不同人体部位分别进行波速和波形分析,为临床实时监测提供更加具体的反应心血管系统状态的信息。

Description

一种无创多路脉搏波检测方法
技术领域
本发明涉及一种无创多路脉搏波检测装置、方法及分析系统。 
背景技术
人体的心脏、动脉、毛细血管及静脉组成心血管系统。在心血管系统中,随着心脏周期性射血,不仅有周期性脉动的血液流速流量在动脉血管中传播,还伴随着血液压力时高时低的脉动以及血管壁时张时缩的振动在动脉血管系中向前传播。这种血液流速、流量以及压力脉动或血管壁振动在动脉血管系中的传播过程就称为脉搏波在动脉中的传播。从力学的角度看,脉搏波是由血液流动、血管壁振动复合而成的,一旦心脏、血液或血管壁组织结构发生变异,必然导致脉搏波出现异常,通过对人体脉搏波的无创检测和分析,可以得到表征心血管系统功能的参数,评价心血管系统的状态,从而实现心血管疾病的早发现、早治疗。 
临床上对脉搏波的关注普遍集中在两点血压的测量上,即舒张压(DP)和收缩压(SP)的测量。研究表明:不同的生理病理现象使脉搏形状发生不同的变化;脉搏波速度(PWV)的大小反映了血管的弹性程度;PWV与舒张压、收缩压、导致动脉硬化的参数、血清总胆固醇含量、高浓度脂蛋白胆固醇含量、血液粘度等直接相关;与舒张压和收缩压相比,PWV与心血管病发病率和死亡率的相关性更大。可见,仅仅停留在两点血压的测量上是不够的,应该在不增加被检测者痛苦的前提下,更深层次的检测和挖掘脉搏波中蕴涵的信息,及早的发现病情。 
目前国内外对无创脉搏检测分析,已经取得众多研究成果,如公开号为CN101176663A的发明专利《基于脉搏波的无创心输出量检测方法和装置》中所公开的内容,公开号为CN201244022Y的实用新型专利《脉搏波测量与生理特征参数分析仪》中所公开的内容等,类似的脉搏波检测方法和装置功能单一,一般只能单一的检测和分析某一个或某几个心血管参数,只能侧面反映心血管系统的状态,而不能够给出明显反应血管 弹性程度的脉搏波传播速度(PWV)。又例如公开号为WO2008/135892A2的国际专利《APPARATUS FOR PERFORMING PULSE WAVE VELOCITYMEASUREMENT》中给出的测量脉搏波传播速度的装置,又只能测量PWV,而没有分析其他心血管参数的功能,且没有给出相应的计算方法和相对成熟的软硬件系统。 
其次,目前的脉搏波检测分析研究成果多集中于一路或两路脉搏波检测,例如本实验室以前的研究成果,如专利公开号为FPI06152的《一种测量脉搏波波速的装置和方法》,以及发表于《生物医学工程学杂志》第25卷第5期的《一种无创脉搏波检测分析系统的研制》中所公开的内容等。心血管病变多发生在某些部位,所以只靠检测一两个检测点的脉搏波信号,难以得到能够全面描述心血管状态的信息,所以有必要在以前的研究内容基础上对分散于人体不同部位的监测点进行多点同步检测,以对比不同部位的脉搏波传播速度。对多路脉搏波信号取平均波形较之一路或两路脉搏波去平均波形进行心血管参数的提取更能够确切的体现人体的心血管系统健康状况。 
发明内容
针对现有技术存在的问题,本发明的目的在于提供一种无创多路脉搏波检测装置、方法及分析系统,通过对脉搏信号的全面检测和综合分析,得到表征心血管系统功能的参数,评价心血管系统的状态,从而实现心血管疾病的早发现、早治疗。 
为实现上述目的,本发明一种无创多路脉搏波检测装置,包括至少两个传感器,每个传感器上均设置有固定装置,并且各个传感器均与信号处理单元连接,传感器测得的信号通过信号处理单元进行处理、储存。 
进一步,所述传感器为生理压力传感器。 
进一步,所述信号处理单元包括顺序连接的放大电路,滤波电路,A/D转换电路和主控制电路,所述传感器分别与相应的放大电路连接。 
进一步,所述信号处理单元还包括用于将传感器信号进行分离的信号分离电路,该信号分离电路将信号分离为静压信号和脉动信号,所述A/D转换电路与计算机连接。 
进一步,所述放大电路由放大芯片构成,其分别对所述传感器采集 的脉搏波信号进行放大。 
进一步,所述信号分离电路由电解电容构成。 
进一步,所述静压信号直接送入所述A/D转换电路进行A/D转换和处理,所述静压信号将直接实时的显示在所述计算机的显示屏幕上。 
进一步,所述脉动信号分别由两级放大电路组成,实现对相应脉动信号的放大。 
进一步,所述滤波电路为RC低通滤波电路,用于对所述放大电路的输出信号进行滤波,滤除频率为50Hz的工频干扰信号。 
进一步,所述A/D转换电路完成对所述静压信号和所述滤波电路的输出信号的模数转换及多次采集求平均值,最终通过数据传输将信号送往所述计算机;所述计算机利用配套的脉搏波信号分析软件完成对脉搏信号的接收、存储、中值滤波、数据分析、波形显示和数据回放工作。 
进一步,所述固定装置包括桡动脉传感器固定装置、颈动脉传感器固定装置、踝动脉传感器固定装置;所述固定装置为单触头固定装置和/或多触头固定装置。 
一种应用上述装置进行脉搏波波速和波形分析的方法,具体为: 
1)调节传感器对动脉血管施加压力的大小,得到最大脉动信号; 
2)同时记录相隔一定距离Δs的多路脉搏波信号; 
3)提取特征参考点; 
4)根据要计算脉搏波速度的位置选择其中两路脉搏波,计算两路脉搏波的时间差Δt; 
5)计算脉搏波速度 v = Δs Δt , 其中Δs为体表测得量检测点间的血管长度,Δt为同时测得的相隔距离为Δs的两路脉搏波对应特征参考点间的时间差; 
6)对脉搏波进行波长平均化,采用线性插值法确定平均化后各点的压力值,并利用水银血压计测量舒张压和收缩压对脉压值进行标定,以计算脉压PP、平均压MP、心率HR、波形特征参数K、心搏出量SV、心输出量CO、外周阻抗R、心指数CI。 
进一步,所述步骤2)中的测量采集点为颈动脉、两侧桡动脉、股动脉和两侧踝动脉中任意两个或多个。
进一步,所述步骤3)中采用上升沿法来确定最终波速参考点,首先 识别出动脉波形中最显著的共同点以确定上升沿部分,然后沿着上升沿部分由中间向两侧来寻找波峰点和波根点位置,并以此作为最终特征点的参考点。 
一种应用上述无创多路脉搏波检测装置和方法的分析系统,包括信息管理模块、数据采集模块、检查分析模块和报告打印模块;信息管理模块用于患者信息登记和查询、波形信息存储、以往检查结果再现;数据采集模块用于对多路脉搏波信号进行同步采集;检查分析模块用于对采集的信号进行波速及波形分析;报告打印模块用于检查结果打印。 
本发明的优点在于: 
1)本发明对人体不同部位的脉搏波采集点采用多路同步采集脉搏波和波速、波形处理技术,可针对不同人体部位分别进行波速和波形分析,为临床实时监测提供更加具体的反应心血管系统状态的信息。 
2)本发明在传统的提取脉搏波特征参考点方法的基础上,提出了改进提取效率和精度的方法。该方法与现有方法相比具有以下优点:适应性好,能够适用于各种脉搏波波形;识别性强;算法简单,对硬件要求低,在采样率一定的情况下保证了一定的精确度;精度高,在波足附近,P值随t变化缓慢,t对P比较敏感,因此P的一个小误差ΔP就有可能带来较大的t误差。在上升沿中提取特征值,t值对P值灵敏度较低,从而能够减小P值误差所带来的t值误差,则精确度会提高很多。 
3)本发明针对不同人体部位脉搏波采集点设计几种不同的传感器固定装置,便于传感器准确和以适当大小的压力接触于检测点上,以准确快速的得到脉搏波形。具体有颈动脉传感器固定装置,桡动脉传感器固定装置,踝动脉传感器固定装置。其中桡动脉固定装置又分为单触头和多触头传感器固定装置。 
4)本发明具有实现所述各种功能的计算机程序产品,可以实现多路脉搏波信号同步采集,患者信息登记和查询、波形信息存储、以往检查结果再现、波速及波形分析,检查结果打印等步骤。 
附图说明
图1为本发明一种无创多路脉搏波检测装置中传感器固定装置与信号处理单元连接结构示意图; 
图2A为桡动脉传感器固定装置示意图; 
图2B为桡动脉两触头传感器固定装置示意图; 
图2C为颈动脉传感器固定装置示意图; 
图2D为踝动脉传感器固定装置示意图; 
图3为本发明一种无创多路脉搏波检测方法中测量及分析脉搏波波速和波形的流程图; 
图4为本发明一种无创多路脉搏波检测装置结构示意图; 
图5为本发明中对获取的脉搏波形判断上升沿方法示意图; 
图6为本发明实施例1患者查询界面图; 
图7为本发明实施例1两路脉搏波波形分析界面示意图; 
图8为本发明一实施例1两路脉搏波波形分析结果及结果打印界面示意图; 
图9为本发明实施例1两路脉搏波分析报告打印结果示意图; 
图10为本发明实施例2六路脉搏波波形显示及选择界面示意图。 
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细描述: 
实施例一: 
如图1、图2A、图2B、图2C、图2D和图4所示的一实施例中一种无创多路脉搏波检测装置,取n=2,即选择采集两路脉搏波。包括: 
第一传感器、第二传感器采用PT14M3生理压力传感器,第一传感器与放大电路相连接,第二传感器与放大电路相连接。
第一传感器放置于颈动脉脉搏波检测传感器固定装置,第二传感器放置于桡动脉脉搏波检测单触头传感器固定装置。分别用于检测颈动脉脉搏波和桡动脉脉搏波。 
放大电路和放大电路主要由放大芯片AD620构成,分别对第一传感器和第二传感器采集的脉搏波信号进行放大。 
第一信号分离电路和第二信号分离电路主要由电解电容构成,第一信号分离电路将放大电路的输出信号分离为第一路静压信号和第一路脉动信号,第二信号分离电路将放大电路的输出信号分离为第二路静压信号和第二路脉动信号。 
第一路静压信号和第二路静压信号直接送入A/D转换器进行模数转换和处理。这两路静压信号值将直接实时的显示在计算机的显示器屏幕上,以便操作者根据当前静压值的大小实时调节静压值,最终得到幅值最大的脉动信号。 
第一脉动放大电路和第二脉动放大电路分别由两级以放大器HA17324A为核心的放大电路组成,分别实现对第一路脉动信号和第二路脉动信号的放大。 
第一滤波电路和第二滤波电路为RC低通滤波电路,分别对第一脉动放大电路和第二脉动放大电路的输出信号进行滤波。滤除频率为50Hz的工频干扰信号。 
A/D转换器完成对4路数据(包括第一路静压信号,第二路静压信号,第一滤波电路3的输出信号,第二滤波电路的输出信号)的模数转换及多次采集求平均,最终通过数据传输设备将信号送往计算机16,此处单片机15可采用单片机16F877A。 
计算机为普通的个人计算机,主要工作是利用配套的脉搏信号分析软件完成脉搏信号的接收、存储、中值滤波、数据分析、波形显示和数据回放等工作。 
一种应用上述装置进行脉搏波波速和波形分析的方法,如图3所示,具体如下: 
1)调节传感器对动脉血管施加压力的大小,得到最大脉动信号;测量时,在计算机屏幕上实时显示当前各脉搏传感器对血管施加压力的大小及各路脉搏波形曲线。观察传感器对动脉血管施加的压力值,调节压力大小,使脉动信号的幅值达到最大。 
2)同时记录相隔一定距离Δs的两路脉搏波信号; 
3)提取特征参考点,首先识别出每一个波形的上升沿,其次通过识别的上升沿来确定特征点。具体方法如下: 
(1)分析判断整个信号数据,找出最大值最小值,利用数字信号分析中的穿级计数法,求出一个能够穿过所有上升沿的数值Pmean。 
(2)在波形数据中寻找每个上升沿中最靠近值Pmean的点,如图5中的A、B、C、D、E等点。 
(3)从A点开始向两侧逐点判断寻找波峰点P波根点P。在此可以根据具体系统和采样率采用合适参数来判断波峰及波根点。 
(4)判断每一对波根波峰代表的上升沿是否属于真实的脉搏波主波,从而剔除错误判断出来的非主波,如B、G点。 
(5)为了减小用此波根作为特征点带来的误差,可以取上升沿的1/5高度点M为特征点。在离散数据中,M点的横坐标位置可以采用线性插值的方法。如图4所示,通过信号数据中最接近M点的两个点L、N的位置来插值求得tM,图中L、N是采样点,M是对应于PM的点。 
4)计算两路脉搏间的时间差Δt: 
对于同时测得的两路脉搏波信号,若第一路脉搏波中完整波形的个数为n1,第二路脉搏中完整波形的个数为n2,两路波形中每一个完整波形对都对应着一个时间差,因此对应于完整波形对的时间差个数为n=min(n1,n2),而实际测得的这n个时间差并不相同,而真实的脉搏波传播时间只有一个,这些测量值只是真实传播时间的一个反映,它们受到被测者情绪或者肌肉和血管当前状态的影响而产生波动,考虑到它们中的大部分只是在一个小范围内波动,只是偶尔有一些波动较大的时间差,因此对这n个时间差进行处理,去掉误差较大的值,最终得到一个最接近真实脉搏传播时间的值。 
若第一对脉搏波的时间差为Δt1,第二对脉搏波的时间差为Δt2,…,第n对脉搏波的时间差为Δtn,最终Δt的确定方式如下: 
(1)求序列Δt1,Δt2,…,Δtn的平均值δt 
δt = 1 n Σ i = 1 n Δ t i
(2)求序列中每个时间差Δti对于均值δt的偏差的绝对值 
p[i]=|Δti-δt|
(3)求序列Δt1,Δt2,…,Δtn的标准偏差σΔt
σ Δt = Σ i = 1 n ( δt - Δt i ) 2 n ( n - 1 )
(4)剔除误差较大的时间差 
若|p[i]|>2σΔt,则认为偏差过大,剔除Δti,回到第一步重新计算,若|p[i]|≤2σΔt,则进入下一步 
(5)确定最终的Δt: 
若最终剩余时间差序列为Δt1,Δt2,…,Δtm,则 
Δt = 1 m Σ i = 1 m Δt i
5)计算脉搏波速度; 
脉搏波速度v定义为两脉搏测量点间的距离Δs与脉搏传播时间Δt的比值,即 
v = Δs Δt
其中Δs为体表测得两检测点间的血管长度。计算脉搏速度关键在于脉搏波传播时间Δt的计算,同时测得的相隔一定距离的两路脉搏波对应参考点间的时间差即脉搏波的传播时间Δt。 
6)分析脉搏波波形 
在确定脉搏波特征点后便可以对波形作进一步的处理和分析。这部分的工作主要包括波形转换、波形平均以及脉压PP、平均压MP、心率HR、波形特征参数K、心搏出量SV、心输出量CO、外周阻抗R、心指数CI等心血管系统参数的计算。 
波形转换 
心脏节律性的收缩和舒张使脉搏波呈周期性的变化,但这里所指的周期性变化仅是一个近似说法,脉搏信号并不是绝对的周期信号,每个心动周期的长度并不完全一样,而是在一个小范围内波动。为方便对脉搏波的分析和比较,需要把不同长度的脉搏波转换为相同长度。若一段脉搏信号中每个波的长度依次为T1,T2,T3,…,Tn,则需把它们转化为长度均为平均脉搏长度T的一段波形。 
系统对脉搏信号的采集频率较高(667Hz),因此数据量较大,若被测者心率为80,则其每个脉搏波的数据量约为60÷80×667=500个点, 完全可以采用线性插值的方法对波形进行转化。具体实现步骤如下: 
1)计算平均脉搏长度T 
T = 1 n Σ i = 1 n T i = 1 n ( t end - t start )
2)线性插值法计算压力值 
若t1时刻压力值为P1,t2时刻压力值为P2,则位于t1和t2之间的t时刻压力值P为: 
P = P 1 + ( P 2 - P 1 ) ( t 1 - t 2 ) ( t - t 1 )
其中(t1,P1)和(t2,P2)为原始波形上的点,(t,P)为转化后波形上的点。 
3)平均波形计算 
平均波是所有脉搏波的平均,平均波上每个点均为所有脉搏波对应点的平均值,因此计算平均波要求每个脉搏波长度相同,此处利用转化为相同长度后的波形进行计算。若转化后每个脉搏波的数据长度为nT,完整脉搏波个数为n,则平均波形上每个点的压力值MP(k)为: 
MP ( k ) = 1 n Σ i = 0 n - 1 P ( k + n T × i ) , ( 1 ≤ k ≤ n T )
波形参数的计算 
波形参数主要包括脉压(PP)、平均压(MP)、心率(HR)、波形特征参数(K)、心搏出量(SV)、心输出量(CO)、外周阻抗(R)和心指数(CI)。下面分别介绍这几项参数的提取。 
(1)脉压 
本系统采用水银柱血压计测得舒张压和收缩压,将测得的血压值输入系统,用于校正数据采集器得到的压力值并参与波形分析。若测得的收缩压为SP,舒张压为DP,则脉压PP为: 
PP=SP-DP 
(2)平均压 
平均压为一段脉博波数据中所有完整波形数据点压力值的平均。若完整波形的起始点为jstart,结束点为jend,则平均压为: 
MP = 1 j end - j start + 1 Σ j = j start j end P [ j ]
此处P[j]是经提取并放大的压力输出值,并不是实际血压值,因此计算所得的MP也不是实际平均压,需要利用输入系统的收缩压SP和舒张压DP对MP进行校正,最终得到实际平均压。数据采集电路对脉动信号的处理如放大、抬高电平等都是线性的,因此在校正压力值时采用线性方法即可。若一段脉搏中有n个完整波形,波峰点为P[jpeak(1)],P[jpeak(2)],…,P[jpeak(n)],波足点为P[jfoot(1)],P[jfoot(2)],…,P[jfoot(n)],具体校正方法如下: 
计算实际测得波峰平均压力值(对应于SP) 
SP ′ = 1 n Σ i = 1 n P [ j peak ( i ) ]
计算实际测得波足平均压力值(对应于DP) 
DP ′ = 1 n Σ i = 1 n P [ j foot ( i ) ]
线性校正,平均压的实际值为: 
MP = MP ′ - DP ′ SP ′ - DP ′ · ( SP - DP ) + DP
(3)心率 
心率为每分钟的心跳次数,是衡量人体心血管系统状态的一个重要参数。若脉搏平均长度为T(单位为s),则心率HR为: 
HR = 60 T
(4)外周阻力 
对于脉搏波数字信号序列P[j],其中j∈[0,N],其外周阻力的计算方法如下: 
分别计算每个独立脉搏波对应的外周阻力: 
R [ i ] = Σ j = j foot ( i ) j foot ( i + 1 ) { P [ j ] - DP ′ SP ′ - DP ′ · ( SP - DP ) } SV
计算序列R[1],R[2],R[3],…,R[n]均值 
R ‾ = 1 n Σ i = 1 n R [ i ]
当n≤8时,以数据序列的平均值作为外周阻力R=R 
当n>8时,采用拉依达判定准则确定R: 
计算序列偏差绝对值:p[i]=|R-R[i]| 
计算序列标准偏差: σ R = Σ i = 1 n ( R ‾ - R [ i ] ) 2 n ( n - 1 )
若p[i]≥2σR,去除R[i],重新计算均值、偏差及偏差绝对值,直到剩下的所有数据点满足p[i]<2σR为止。 
将最终序列均值作为外周阻力R=R 
(5)波形特征参数K 
舒张压、收缩压和平均压已知的情况下,带入公式可直接计算出波形特征参数: 
K = MP - DP SP - DP
(6)心搏出量SV 
测得脉压、体重、年龄、心率,带入经验公式得到心搏出量: 
SV=PP×(0.013×体重-0.007×年龄-0.004×心率+1.307) 
(7)心输出量CO 
根据心搏出量和心率,计算心输出量: 
CO=SV×HR
(8)心指数CI 
根据心输出量和体表面积,计算心指数: 
CI = CO BSA
其中,BSA根据经验公式计算得到: 
中国成年男性:BSA=0.00607×身高+0.0127×体重-0.0698 
中国成年女性:BSA=0.00586×身高+0.0126×体重-0.0461 
本发明中应用无创多路脉搏波检测装置和方法的分析系统的操作方法如下: 
1)连接并启动系统 
连接好计算机的主机和显示器,将计算机与数据获取装置用RS232串行接口连接(不接电源),将数据获取装置和脉搏波探头连接,为计算机和数据获取装置接上电源。打开计算机启动脉搏检测与分析系统程序,进入主界面。 
2)检测脉搏波 
点击“进入系统”按钮,进入“引导检查界面”; 
对于未登记患者,点击“新患者登记”按钮->进入“新患者登记界面”; 
对于已登记患者,点击“患者查询”按钮->进入“患者查询界面”,患者查询界面如图6所示; 
“病史及药物说明”栏:同时按住“Ctrl+Enter”换行; 
点击“保存”按钮->保存填写的数据; 
点击“检查”按钮->进入“患者检查界面”; 
选择“检查类别”; 
点击“开始|停止”按钮,开始检查; 
点击“开始|停止”按钮,停止检查; 
点击“保存”按钮,保存数据; 
3)分析波形,分析波形界面如图7所示 
点击“选择分析范围”按钮; 
移动鼠标到要选择的起始点,点击鼠标右键; 
移动鼠标到要选择的终止点,点击鼠标右键(如果略去第五步,默认分析全部数据); 
点击“分析”按钮,进入结果分析界面; 
4)打印分析结果,分析结果显示界面如图8所示; 
5)点击“退出”按钮->退出到主界面。 
实施例2 
本实施例不同于实施例1的特点在于,本实施例中选择n=6,即采用六路脉搏波进行脉搏波波速及波形的分析。本实施例中一种无创多路脉搏波检测装置,包括: 
第一传感器放置于颈动脉脉搏波检测传感器固定装置,第二、三传感器放置于桡动脉脉搏波检测单触头传感器固定装置,第四、五传感器放置于踝动脉脉搏波检测传感器固定装置,第六传感器待其他传感器固定好后用手持方式固定在检测位置。分别用于检测颈动脉脉搏波和两侧桡动脉脉搏波,两侧踝动脉脉搏波,和股动脉脉搏波。 
六路放大电路主要由放大芯片AD620构成,分别对六路采集的脉搏波信号进行放大。 
各信号分离电路主要由电解电容构成,将对应的前置放大电路的输出信号分离为相应的静压信号和脉动信号。 
第一至第六路各路静压信号直接送入单片机进行A/D转换和处理。静压信号值将直接实时的显示在计算机的显示器屏幕上,以便操作者根据当前静压值的大小实时调节静压值,最终得到幅值最大的脉动信号。 
第一至第六路各路脉动信号分别由两级以放大器HA17324A为核心的放大电路组成,分别实现对个路脉动信号的放大。 
第一至第六路滤波电路为RC低通滤波电路,分别对各路分理处的脉动信号的输出信号进行滤波。滤除频率为50Hz的工频干扰信号。 
A/D转换器完成对12路数据(包括六路静压信号,六路滤波电路的输出信号)的模数转换及多次采集求平均,最终通过USB接口将信号送往计算机,此处A/D转换器可采用单片机16F877A。 
计算机为普通的个人计算机,主要工作是利用配套的脉搏信号分析软件完成脉搏信号的接收、存储、中值滤波、数据分析、波形显示和数据回放等工作。 
一种应用上述装置进行脉搏波波速和波形分析的方法,如图3所示,具体如下: 
1)调节各传感器对动脉血管施加压力的大小,得到最大脉动信号;测量时,在计算机屏幕上实时显示当前各脉搏传感器对血管施加压力的大小及各路脉搏波形曲线。观察传感器对动脉血管施加的压力值, 调节压力大小,使脉动信号的幅值达到最大。 
2)同时记录位于不同检测点的六路脉搏波信号; 
3)对每路放大滤波后的脉动信号提取特征参考点,首先识别出每一个波形的上升沿,其次通过识别的上升沿来确定特征点。具体方法参考实施例1。 
4)选择要分析波速的部位,例如颈-股间脉搏波波速,颈-左桡动脉间脉搏波波速等,针对已选择的两路脉搏波信号计算两路脉搏间的时间差Δt:Δt的计算方法参考实施例1 
5)计算脉搏波速度; 
脉搏波速度v定义为两脉搏测量点间的距离Δs与脉搏传播时间Δt的比值,即 
v = Δs Δt
其中Δs为体表测得两检测点间的血管长度。 
6)分析脉搏波波形的方法与实施例1相同。 
本发明中应用无创多路脉搏波检测装置和方法的分析系统的操作方法如下: 
1)连接并启动系统 
连接好计算机的主机和显示器,将计算机与数据获取装置用USB数据线连接(不接电源),将数据获取装置和脉搏波探头连接,为计算机和数据获取装置接上电源。打开计算机启动脉搏检测与分析系统程序,进入主界面。 
2)检测脉搏波 
点击“进入系统”按钮,进入“引导检查界面”; 
对于未登记患者,点击“新患者登记”按钮->进入“新患者登记界面”; 
对于已登记患者,点击“患者查询”按钮->进入“患者查询界面”,患者查询界面如图6所示; 
“病史及药物说明”栏:同时按住“Ctrl+Enter”换行; 
点击“保存”按钮->保存填写的数据; 
点击“检查”按钮->进入“患者检查界面”; 
选择“检查类别”; 
点击“开始|停止”按钮,开始检查; 
点击“开始|停止”按钮,停止检查; 
点击“保存”按钮,保存数据.; 
3)分析波形,如图9、图10所示 
点击选择要进行分析的两路脉搏波波形; 
点击“选择分析范围”按钮; 
移动鼠标到要选择的起始点,点击鼠标右键; 
移动鼠标到要选择的终止点,点击鼠标右键(如果略去第五步,默认分析全部数据); 
点击“分析”按钮,进入结果分析界面; 
重复以上步骤重新选择要分析的两路脉搏波波形; 
4)打印分析结果; 
5)点击“退出”按钮->退出到主界面。 

Claims (2)

1.一种脉搏波波速和波形分析的非诊断性方法,具体为:
1)调节传感器对动脉血管施加压力的大小,得到最大脉动信号;
2)同时记录相隔一定距离Δs的多路脉搏波信号;
3)提取特征参考点;
4)根据要计算脉搏波速度的位置选择其中两路脉搏波,计算两路脉搏波的时间差Δt;
5)计算脉搏波速度
Figure FSB00000925699400011
其中Δs为体表测得量检测点间的血管长度,Δt为同时测得的相隔距离为Δs的两路脉搏波对应特征参考点间的时间差;
6)对脉搏波进行波长平均化,采用线性插值法确定平均化后各点的压力值,并利用水银血压计测量舒张压和收缩压对脉压值进行标定,以计算脉压PP、平均压MP、心率HR、波形特征参数K、心搏出量SV、心输出量CO、外周阻抗R、心指数CI;
其中,步骤3)中采用上升沿法来确定最终波速参考点,首先识别出动脉波形中最显著的共同点以确定上升沿部分,然后沿着上升沿部分由中间向两侧来寻找波峰点和波根点位置,并以此作为最终特征点的参考点。
2.如权利要求1所述的非诊断性方法,其特征在于,所述步骤2)中的测量采集点为颈动脉、两侧桡动脉、股动脉和两侧踝动脉中任意两个或多个。
CN 200910243060 2009-12-24 2009-12-24 一种无创多路脉搏波检测方法 Expired - Fee Related CN101732040B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200910243060 CN101732040B (zh) 2009-12-24 2009-12-24 一种无创多路脉搏波检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200910243060 CN101732040B (zh) 2009-12-24 2009-12-24 一种无创多路脉搏波检测方法

Publications (2)

Publication Number Publication Date
CN101732040A CN101732040A (zh) 2010-06-16
CN101732040B true CN101732040B (zh) 2013-03-27

Family

ID=42456437

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200910243060 Expired - Fee Related CN101732040B (zh) 2009-12-24 2009-12-24 一种无创多路脉搏波检测方法

Country Status (1)

Country Link
CN (1) CN101732040B (zh)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102525431B (zh) * 2010-12-10 2014-01-08 李祥 一种心血管生理信号的检测装置及方法
CN102334981B (zh) * 2011-07-23 2013-11-27 浙江大学 一种人体多路脉搏波信号采集系统
CN103099610B (zh) * 2011-11-11 2015-05-13 杭州电子科技大学 基于左右肱动脉脉搏波传输时间差的动态血压测量装置及方法
CN102551698B (zh) * 2011-12-16 2015-03-18 广州医学院第二附属医院 颈、脑动脉脉搏波速度测量系统
CN105188522B (zh) * 2013-03-08 2019-06-11 富士胶片株式会社 脉搏波传播速度的测定方法、测定系统及摄像装置
CN103400069B (zh) * 2013-07-29 2016-01-13 东北大学 一种基于功率谱分析的脉搏波身份识别装置及方法
CN104735297B (zh) * 2013-12-20 2017-10-10 南京视威电子科技股份有限公司 一种广播级监视器及监视图像生成方法
CN106231995B (zh) * 2014-04-14 2020-07-03 株式会社村田制作所 脉搏波传播时间测量装置及生命体状态推定装置
CN204515353U (zh) 2015-03-31 2015-07-29 深圳市长桑技术有限公司 一种智能手表
WO2016107607A1 (en) 2015-01-04 2016-07-07 Vita-Course Technologies Co.,Ltd System and method for health monitoring
CN107847153B (zh) * 2015-07-03 2020-12-04 深圳市长桑技术有限公司 一种生理参数监测的系统和方法
CN104958064A (zh) * 2015-07-15 2015-10-07 四川宇峰科技发展有限公司 一种可穿戴式动脉硬化检测仪及脉搏波传导速度检测方法
TWI669096B (zh) * 2017-07-13 2019-08-21 國立臺灣大學 具有確定頸動脈血壓的多功能量測裝置
CN107961000A (zh) * 2017-12-29 2018-04-27 西安城派信息科技有限公司 缺血预适应仪
CN108742588B (zh) * 2018-04-12 2020-11-20 中国医学科学院生物医学工程研究所 一种脉搏波波速检测方法及装置
CN108836281A (zh) * 2018-06-28 2018-11-20 中国科学院微电子研究所 用于脉搏特征参数采集的系统
WO2020132808A1 (zh) * 2018-12-24 2020-07-02 深圳迈瑞生物医疗电子股份有限公司 针对血流动力学的生理体征监测方法和医疗监护设备
CN109497991B (zh) * 2019-01-04 2021-07-02 济南汇医融工科技有限公司 基于机器学习的心衰检测装置
CN109875544B (zh) * 2019-01-31 2021-09-17 湖北省体育科学研究所 用于运动员精准选材的心脏机能评定方法
CN110811599A (zh) * 2019-11-19 2020-02-21 北京东方逸腾数码医疗设备技术有限公司 一种信号同步采集系统
CN110974193A (zh) * 2019-12-30 2020-04-10 北京中科芯健医疗科技有限公司 一种人体功能状态评估装置
CN111150385A (zh) * 2020-01-15 2020-05-15 上海科技大学 一种基于压电薄膜的可穿戴式血液流速传感器
CN111558148B (zh) * 2020-03-23 2023-06-20 未来穿戴技术有限公司 颈部按摩仪的健康检测方法及颈部按摩仪
CN112274126B (zh) * 2020-10-28 2022-11-29 河北工业大学 一种基于多路脉搏波的无创连续血压检测方法、装置
CN112842291B (zh) * 2021-01-29 2022-11-18 清华大学深圳国际研究生院 一种脉搏波速测量系统及无创式血流状况评估系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2766772Y (zh) * 2005-01-14 2006-03-29 孙汉钧 一种人体健康信息监测传感器及其装置
CN101084828A (zh) * 2006-06-05 2007-12-12 中国科学院力学研究所 一种测量脉搏波速度的装置和方法
CN101496718A (zh) * 2009-02-27 2009-08-05 山东大学 颈动脉脉搏波检测传感器调整装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2766772Y (zh) * 2005-01-14 2006-03-29 孙汉钧 一种人体健康信息监测传感器及其装置
CN101084828A (zh) * 2006-06-05 2007-12-12 中国科学院力学研究所 一种测量脉搏波速度的装置和方法
CN101496718A (zh) * 2009-02-27 2009-08-05 山东大学 颈动脉脉搏波检测传感器调整装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李婷,虞钢.一种无创脉搏波检测分析系统的研制.《生物医学工程学杂志》.2008,第25卷(第5期), *

Also Published As

Publication number Publication date
CN101732040A (zh) 2010-06-16

Similar Documents

Publication Publication Date Title
CN101732040B (zh) 一种无创多路脉搏波检测方法
CN105943005B (zh) 基于光电绿光脉搏与心电图混合的无创血压检测装置
CN101357062B (zh) 基于容积脉搏波检测的生命体征参数测评装置
CN103892818B (zh) 一种无创中心动脉血压测量方法和设备
CN103479343B (zh) 基于振荡式血压计信号的中心动脉压检测系统
CN102302361B (zh) 一种心率及身体状态监测装置及方法
WO2017024457A1 (zh) 连续血压测量装置、测量模型建立方法和系统
CN103054562B (zh) 基于多路脉搏波波形分析的心血管功能检测方法及其装置
TWI250867B (en) Pulse wave analysis device
CN101176663B (zh) 基于脉搏波的无创心输出量检测方法和装置
US20190175042A1 (en) Method and apparatus for determining blood pressure
CN109512412B (zh) 一种中心动脉血压测量装置
CN201213789Y (zh) 人体上、下肢及主动脉pwv值同时检测显示系统
CN201088579Y (zh) 一种动脉硬化检测和评估装置
CN104757959A (zh) 一种基于影像叠影的脉搏波传输速率检测方法和系统
Kumar et al. Traditional practices and recent advances in Nadi Pariksha: A comprehensive review
CN106539562A (zh) 一种由脉搏波评估心功能的方法
CN112155514A (zh) 一种基于检测脉搏生物电信息的非接触眼压测量方法
Yan et al. A non-invasive blood pressure prediction method based on pulse wave feature fusion
CN113069091A (zh) 一种ppg信号的脉象分类装置及方法
KR20080030189A (ko) 혈관의 건강 상태를 감시하는 방법 및 장치
CN104027097A (zh) 血管功能无创检测方法及装置
CN114515147B (zh) 一种基于bcg信号与ppg信号融合的生理监测系统
CN105310678A (zh) 一种基于脉波分析法计算心脏每搏血量的检测方法
CN115316966A (zh) 一种基于片上光学微腔的血压探测装置及方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130327

Termination date: 20181224