CN102711597B - 脉搏波传播速度测定装置以及脉搏波传播速度测定程序 - Google Patents
脉搏波传播速度测定装置以及脉搏波传播速度测定程序 Download PDFInfo
- Publication number
- CN102711597B CN102711597B CN201080050888.XA CN201080050888A CN102711597B CN 102711597 B CN102711597 B CN 102711597B CN 201080050888 A CN201080050888 A CN 201080050888A CN 102711597 B CN102711597 B CN 102711597B
- Authority
- CN
- China
- Prior art keywords
- mentioned
- pulse wave
- blood
- fiducial time
- 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.)
- Expired - Fee Related
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/026—Measuring blood flow
- A61B5/0285—Measuring or recording phase velocity of blood waves
-
- 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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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/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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7239—Details of waveform analysis using differentiation including higher order derivatives
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Vascular Medicine (AREA)
- Hematology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明的脉搏波传播速度测定装置(10)用基准时间检测部(2)来检测脉搏波(S)的射血波分量(S1)、反射分量(S2)的基准时间(T1、T2),并且用脉搏波振幅检测部(3)来检测与射血波分量(S1)、反射分量(S2)的基准时间(T1、T2)相对应的脉搏波(S)的振幅(W1、W2)。并且,传播速度检测部(5)基于射血波分量(S1)、反射波分量(S2)的基准时间(T1、T2)、和与射血波分量(S1)、反射波分量(S2)的基准时间(T1、T2)相对应的脉搏波(S)的振幅(W1、W2),来求出脉搏波(S)的传播速度(PWV1)。因此,能够考虑到射血波分量(S1)的振幅与反射分量(S2)的振幅的差异所导致的射血波与反射波的脉搏波传播速度的差异来高精度地对脉搏波传播速度进行测定。
Description
技术领域
本发明涉及通过对生物体的脉搏波进行测定,来算出脉搏波的传播速度的脉搏波传播速度测定装置以及脉搏波传播速度测定程序。
背景技术
在现有技术中,已知脉搏波在把握生物体的循环系统的状态方面具有各种各样的重要的信息。尤其是,脉搏波在生物体的2个部位进行传播的速度以及时间,被指出具有能够把握动脉硬化状态等的可能性,是在医疗现场也受到关注的生物体指标,分别被称作脉搏波传播速度(PWV:PulseWave Velocity)、脉搏波传播时间(PTT:Pulse Transit Time)等。
在该脉搏波传播速度中,根据测定部位提出了多种测定手法,例如,颈动脉和大腿动脉间的脉搏波传播速度被称作cfPWV(carotid-femoralPWV,颈动脉-股动脉脉搏波传播速度),作为脉搏波传播速度(PWV)中的黄金标准(gold standard)而被利用。
一般来说,为了求出脉搏波传播速度(PWV),需要2个部位的脉搏波测定点。
但是,已知无论是在生物体上的哪个部位测定出的脉搏波,都是来自心脏的射血波和从生物体内的各个部位反射的反射波的合成波。并且,通过将该射血波和反射波分离,即使脉搏波的测定部位为1处,也存在能够求出脉搏波传播速度或者脉搏波传播时间的可能性。
因此,在非专利文献1中,公开了一种用于将射血波和反射波分离的技术。此外,在专利文献1和专利文献2中,公开了一种通过将射血波和反射波分离,即使脉搏波的测定部位为1处,也能求出脉搏波传播速度或者脉搏波传播时间的技术。
一般来说,脉搏波在生物体内的各处产生反射。该脉搏波的反射,主要起因于血管的阻抗不整合,例如,在存在血管的分支或血管的弹性力的变化等的部位产生脉搏波的反射。
在此,如同在专利文献1中也记述了的那样,在将射血波和反射波分离时,若假定主要的反射点位于肠骨动脉或者腹部大动脉周边,则针对在生物体各部测定出的脉搏波,射血波和反射波的分离顺利地进行。
此外,在专利文献2中,公开了如下技术:脉搏波传播速度或者脉搏波传播时间与血压存在相关关系,根据从1个测定部位得到的脉搏波来算出脉搏波传播速度或者脉搏波传播时间,并识别血压。此外,在非专利文献2中公开了如下技术:根据从2个测定部位得到的脉搏波来算出脉搏波传播速度或者脉搏波传播时间,并识别血压。
另外,在上述专利文献1以及2的任意一者中,根据1个测定部位的脉搏波都能够求出脉搏波传播速度或者脉搏波传播时间,但其前提是射血波和反射波各自的脉搏波传播速度相同。按照该前提,(1)将1个脉搏波分离为射血波和反射波并求出两者的时间差,(2)求出射血波和反射波各自的脉搏波所传播的距离的差。由此,求出了脉搏波传播速度或者脉搏波传播时间。
但是,实际上,射血波和反射波的脉搏波传播速度不一致。原因在于,射血波和反射波的振幅不同。如前述的专利文献2中记载的那样,脉搏波传播速度与血压存在相关关系,而血压与脉搏波的振幅有关。即,脉搏波传播速度与脉搏波的振幅有关。
如上所述,反射波通过将反射点假定为腹部大动脉周边能够很好地说明脉搏波的波形形状。发生脉搏波的反射的理由在于,存在大动脉和腹部大动脉的阻抗不整合,但该不整合的程度因人而异,而且即使是同一人物,也根据健康状态/血管状态等而不同。其结果,反射波的振幅与行波不同,而且其不同程度因人而异。因此,射血波和反射波的脉搏波传播速度不同。
在先技术文献
专利文献
专利文献1:JP特开2003-010139号公报
专利文献2:JP特开2007-007075号公报
非专利文献1:Takazawa K et al.“Underestimation of vasodilatoreffects of nitroglycerin by upper limb blood pressure”,Hypertension 1995;26:520-3
非专利文献2:McCombie,Devin“Development of a wearable bloodpressure monitor using adaptive calibration of peripheral pulse transit timemeasurements”,Ph.D.Thesis,Massachusetts Institute of Technology,Dept.of Mechanical Engineering,2008.
发明内容
发明要解决的课题
因此,本发明的课题在于,提供一种能够根据1个测定部位的脉搏波来更准确地求出脉搏波传播速度的脉搏波传播速度测定装置以及脉搏波传播速度测定程序。
用于解决课题的手段
为了解决上述课题,本发明的脉搏波传播速度测定装置的特征在于,具备:
脉搏波检测部,其检测生物体的某一个部位的脉搏波;
基准时间检测部,其检测用于确定由上述脉搏波检测部检测出的上述一个部位的脉搏波中包含的射血波分量的基准时间和用于确定上述脉搏波中包含的反射波分量的基准时间;
脉搏波振幅检测部,其检测与由上述基准时间检测部检测出的上述射血波分量的基准时间相对应的上述脉搏波的振幅,并且检测与由上述基准时间检测部检测出的上述反射波分量的基准时间相对应的上述脉搏波的振幅;和
传播速度检测部,其基于由上述基准时间检测部检测出的上述射血波分量的基准时间以及上述反射波分量的基准时间、和由上述脉搏波振幅检测部检测出的与上述射血波分量的基准时间相对应的脉搏波的振幅以及与上述反射波分量的基准时间相对应的脉搏波的振幅,来求出上述脉搏波的传播速度。
本发明的脉搏波传播速度测定装置,用上述基准时间检测部来检测脉搏波的射血波分量以及反射分量的基准时间,并且用上述脉搏波振幅检测部来检测与上述射血波分量以及反射波分量的基准时间相对应的脉搏波的振幅。并且,上述传播速度检测部基于上述射血波分量、反射波分量的基准时间、和与该射血波分量、反射波分量的基准时间相对应的脉搏波的振幅,来求出上述脉搏波的传播速度。因此,能够考虑到射血波分量的振幅与反射波分量的振幅的差异所导致的射血波与反射波的脉搏波传播速度的差异来高精度地对脉搏波传播速度进行测定。
此外,在一个实施方式的脉搏波传播速度测定装置中,上述脉搏波振幅检测部具有射血波分量除去部,
该射血波分量除去部求出与上述反射波分量的基准时间相对应的上述脉搏波的振幅中包含的射血波分量,从与上述反射波分量的基准时间相对应的上述脉搏波的振幅中除去上述脉搏波的振幅中包含的射血波分量,来求出与上述反射波分量的基准时间相对应的上述反射波分量的振幅,
上述传播速度检测部,
基于由上述基准时间检测部检测出的上述射血波分量的基准时间以及上述反射波分量的基准时间、和由上述脉搏波振幅检测部检测出的与上述射血波分量的基准时间相对应的脉搏波的振幅、以及由上述射血波分量除去部求出的上述反射波分量的振幅,来求出上述脉搏波的传播速度。
根据本实施方式,由上述射血波分量除去部从与上述反射波分量的基准时间相对应的上述脉搏波的振幅中除去上述脉搏波的振幅中包含的射血波分量,来求出与上述反射波分量的基准时间相对应的上述反射波分量的振幅,因此能够更准确地求出与上述反射波分量的基准时间相对应的上述反射波分量的振幅。因此,能够通过传播速度检测部,更准确地求出上述脉搏波的传播速度。
此外,在一个实施方式的脉搏波传播速度测定装置中,上述传播速度检测部,
基于由上述基准时间检测部检测出的上述射血波分量的基准时间与上述反射波分量的基准时间的时间差、和由上述脉搏波振幅检测部检测出的与上述射血波分量的基准时间相对应的脉搏波的振幅与由上述射血波分量除去部求出的上述反射波分量的振幅的振幅差,来求出上述脉搏波的传播速度。
根据本实施方式,能够基于由上述基准时间检测部检测出的上述射血波分量的基准时间与上述反射波分量的基准时间的时间差、和由上述脉搏波振幅检测部检测出的与上述射血波分量的基准时间相对应的脉搏波的振幅与由上述射血波分量除去部求出的上述反射波分量的振幅的振幅差,来更准确地求出上述脉搏波的传播速度。
此外,在一个实施方式的脉搏波传播速度测定程序中,使计算机执行如下功能:
基准时间导出功能,其求出用于确定生物体的某部位的脉搏波中包含的射血波分量的基准时间和用于确定上述脉搏波中包含的反射波分量的基准时间;
脉搏波振幅导出功能,其求出与上述射血波分量的基准时间相对应的上述脉搏波的振幅,并且求出与上述反射波分量的基准时间相对应的上述脉搏波的振幅;和
传播速度导出功能,其基于上述射血波分量的基准时间以及上述反射波分量的基准时间、和与上述射血波分量的基准时间相对应的脉搏波的振幅以及与上述反射波分量的基准时间相对应的脉搏波的振幅,来算出上述脉搏波的传播速度。
根据本实施方式的脉搏波传播速度测定程序,使计算机执行上述基准时间导出功能,来求出用于确定脉搏波的射血波分量和反射分量的基准时间,并且通过上述脉搏波振幅导出功能来求出与用于确定上述射血波分量以及反射分量的基准时间相对应的脉搏波的振幅。并且,通过上述传播速度导出功能,基于用于确定上述射血波分量、反射分量的基准时间、和与用于确定上述射血波分量、反射分量的基准时间相对应的脉搏波的振幅,来算出上述脉搏波的传播速度。因此,能够考虑到射血波分量的振幅与反射分量的振幅的差异所导致的射血波与反射波的脉搏波传播速度的差异来更高精度地检测脉搏波传播速度。
此外,在一个实施方式的脉搏波传播速度测定程序中,上述脉搏波振幅导出功能包含射血波分量除去功能,
该射血波分量除去功能求出与上述反射波分量的基准时间相对应的上述脉搏波的振幅中包含的射血波分量,从与上述反射波分量的基准时间相对应的上述脉搏波的振幅中除去上述脉搏波的振幅中包含的射血波分量,来求出与上述反射波分量的基准时间相对应的上述反射波分量的振幅,
上述传播速度导出功能,
基于上述射血波分量的基准时间以及上述反射波分量的基准时间、和与上述射血波分量的基准时间相对应的脉搏波的振幅、以及通过上述射血波分量除去功能求出的上述反射波分量的振幅,来算出上述脉搏波的传播速度。
根据本实施方式的脉搏波传播速度测定程序,通过上述射血波分量除去功能来从与上述反射波分量的基准时间相对应的上述脉搏波的振幅中除去上述脉搏波的振幅中包含的射血波分量,来求出与上述反射波分量的基准时间相对应的上述反射波分量的振幅。由此,能够更准确地求出与上述反射波分量的基准时间相对应的上述反射波分量的振幅。因此,能够通过传播速度导出功能,更准确地求出上述脉搏波的传播速度。
此外,在一个实施方式的脉搏波传播速度测定程序中,上述传播速度导出功能,
基于上述射血波分量的基准时间与上述反射波分量的基准时间的时间差、和与上述射血波分量的基准时间相对应的脉搏波的振幅与通过上述射血波分量除去功能求出的上述反射波分量的振幅的振幅差,来求出上述脉搏波的传播速度。
根据本实施方式,能够基于上述射血波分量的基准时间与上述反射波分量的基准时间的时间差、和与上述射血波分量的基准时间相对应的脉搏波的振幅与通过上述射血波分量除去功能求出的上述反射波分量的振幅的振幅差,来更准确地求出上述脉搏波的传播速度。
发明的效果
根据本发明的脉搏波传播速度测定装置,基于用于确定射血波分量、反射分量的基准时间、和与该用于确定射血波分量、反射分量的基准时间相对应的脉搏波的振幅,来求出脉搏波的传播速度,因此能够考虑到射血波分量的振幅与反射分量的振幅的差异所导致的射血波与反射波的脉搏波传播速度的差异,来高精度地检测脉搏波传播速度。
附图说明
图1是本发明的实施方式的脉搏波传播速度测定装置的框图。
图2A是表示用上述实施方式的脉搏波检测部检测出的脉搏波的波形的波形图。
图2B是表示上述脉搏波的加速度波形的波形图。
图2C是表示用上述实施方式的脉搏波检测部检测出的脉搏波的3次微分波形的波形图。
图2D是表示上述脉搏波的4次微分波形的波形图。
图3是表示上述脉搏波的波形以及上述脉搏波的基准点Q1、Q2上的脉搏波的振幅W1、W2的波形图。
图4是表示上述脉搏波的波形以及上述脉搏波的射血波分量S1和反射波分量S2的波形图。
图5是示意性地表示上述脉搏波S的射血波分量S1在人的生物体内传播的路径d和反射波分量S2在人的生物体内传播的路径h的示意图。
具体实施方式
以下,通过图示的实施方式对本发明更详细地进行说明。
在图1的框图中,示出作为本发明的实施方式的脉搏波传播速度测定装置10。该脉搏波传播速度测定装置10具备脉搏波检测部1,该脉搏波检测部1检测人类的生物体的某一个部位的脉搏波。作为该脉搏波检测部1,例如,存在用受光元件来测定从发光元件输出的红外光根据血管内的血液量而被反射或者吸收的程度的部件(光电容积脉搏波法),但除此之外,也可以是利用将血管内的血液压迫血管的压力的变化作为电气信号而取出的压力脉搏波法来检测脉搏波的部件等。此外,用上述脉搏波检测部1来测定脉搏波的生物体部位,虽然没有特别大的限制事项,但优选尽量为非侵入·非束缚的部位,例如,优选指尖·手腕·耳朵等。
此外,该脉搏波传播速度测定装置10,具备:射血波·反射波信息提取部6,其检测用于分别确定由上述脉搏波检测部1检测出的上述脉搏波中包含的射血波分量、反射波分量的基准时间、和与上述基准时间相对应的上述脉搏波的振幅;和脉搏波传播速度检测部5,其基于来自该信息提取部6的表示上述基准时间以及上述振幅的信息,来求出上述脉搏波的传播速度。
上述射血波·反射波信息提取部6具有基准时间检测部2,该基准时间检测部2检测用于确定由上述脉搏波检测部1检测出的上述脉搏波中包含的射血波分量的基准时间、和用于确定上述脉搏波中包含的反射波分量的基准时间。
以下对该基准时间检测部2求出上述基准时间的过程的一例进行说明。首先,在图2A中示出由上述脉搏波检测部1检测出的脉搏波的波形的一例。图2A中的纵轴是与脉搏波的振幅(mmHg)相对应的测定电压值(V)。在本实施方式的脉搏波传播速度测定装置10中,作为一例,利用由一般的袖带式血压计测定出的血压(mmHg),来进行由脉搏波检测部1测定出的电压值(V)与脉搏波的振幅(mmHg)如何对应的校正(校准)。另外,该校正(校准)只要在最初的使用开始时进行即可,在之后的测定中只要使用上述校准的结果即可。
上述基准时间检测部2,例如,在上述脉搏波为Murgo等人进行的血压波形分类的TypeC的情况下,检测图2A所示的脉搏波的波形上的收缩期的极大点Q1的时间T1作为射血波分量的基准时间T1。在图2B中,示出上述脉搏波的加速度波,在图2C中,示出上述脉搏波的3次微分波。并且,上述基准时间检测部2,例如,在上述脉搏波为Murgo等人所进行的血压波形分类的TypeC的情况下,检测图2D所示的脉搏波的4次微分波的第3零交叉点(zero cross point)Q2作为反射波分量的基准时间T2。另外,该第3零交叉点Q2,意味着在图2A所示的脉搏波成为极小值以后,图2D所示的4次微分波形第3次向下进行零交叉的点。
此外,在上述的说明中,对上述检测出的脉搏波为上述血压波形分类的TypeC的情况进行了说明,但在上述检测出的脉搏波为上述血压波形分类的TypeC以外的波形形状的情况下,如果存在能够更合适地确定脉搏波的射血波分量的基准时间和反射波分量的基准时间的手法则也可以采用。例如,脉搏波若没有较大的血压变动,则基本上不会示出很大的波形变化,因此通过对测定出的多个脉搏波进行叠加(加法平均),能够实现脉搏波检测精度的改善。
此外,脉搏波根据噪声电平、年龄·性别·疾病的有无·健康状况等,而显示出各种各样的波形形状,因此如果存在能够更合适地确定脉搏波的射血波分量的基准时间和反射波分量的基准时间的手法则也可以采用。例如,通过将检测出的脉搏波与该时间点的被测定者的状态一起作为历史记录进行保留,能够实现对脉搏波的射血波分量的基准时间和反射波分量的基准时间进行确定的精度的改善。
此外,上述射血波·反射波信息提取部6具有脉搏波振幅检测部3。该脉搏波振幅检测部3,如图3的波形图所例示的那样,检测与由上述基准时间检测部2检测出的上述射血波分量的基准时间T1相对应的上述脉搏波S的振幅W1,并且检测与由上述基准时间检测部2检测出的上述反射波分量的基准时间T2相对应的上述脉搏波S的振幅W2。
并且,上述射血波·反射波信息提取部6具有射血波分量除去部4。该射血波分量除去部4,如图4所示,从与上述反射波分量S2的基准时间T2相对应的上述脉搏波S的振幅W2中除去上述脉搏波S的振幅W2中包含的射血波分量S1,来求出与上述反射波分量S2的基准时间T2相对应的上述反射波分量S2的振幅W3。作为求出该与反射波分量S2的基准时间T2相对应的上述反射波分量S2的振幅W3的手法,例如,可以考虑利用Windkessel模型等,对脉搏波中的射血波的减少程度进行模型化,并从在反射波分量S2的基准时间T2的周边测定出的脉搏波振幅中减去射血波S1的残存分量的方法等。当然,如果存在能够更准确地确定射血波S1的残存分量的手法则也可以采用。
并且,上述射血波·反射波信息提取部6将表示由上述基准时间检测部2检测出的上述射血波分量的基准时间T1以及上述反射波分量的基准时间T2的信息输入到上述脉搏波传播速度检测部5,并且将表示由上述脉搏波振幅检测部3检测出的与上述射血波分量的基准时间T1相对应的脉搏波S的振幅W1以及与上述反射波分量的基准时间T2相对应的上述脉搏波S的反射波分量S2的振幅W3的信息输入到上述脉搏波传播速度检测部5。于是,上述脉搏波传播速度检测部5基于上述射血波分量的基准时间T1、反射波分量的基准时间T2、与上述射血波分量的基准时间T1相对应的脉搏波S的振幅W1、和与上述反射波分量S2的基准时间T2相对应的上述脉搏波S的反射波分量S2的振幅W3,来求出上述脉搏波S的传播速度。
接着,对上述脉搏波传播速度检测部5求出上述脉搏波S的传播速度PWV的过程进行说明。
一般来说,所谓脉搏波传播速度,是对生物体部位的2个部位各自的脉搏波进行测定,并求出各自的脉搏波的射血波分量进行传播的速度,其根本原理是基于Moens-Korteweg的原理。此外,在非专利文献2中,根据Moens-Korteweg的关系导出了脉搏波传播速度与血压的关系,并用下式(1)来表示。
(PWV)2=α·exp(β×P) ...(1)
在上式(1)中,PWV是脉搏波传播速度(m/秒),P是血压(mmHg),α、β是按照每个个人、在个人内也按照每个测定时间而稍有变化的常数。
在此,对射血波S1的脉搏波传播速度PWV1(m/秒)和反射波S2的脉搏波传播速度PWV2(m/秒)分别应用上式(1),则能够得到下式(2)、(3)。
(PWV1)2=α·exp(β×W1) ...(2)
(PWV2)2=α·exp(β×W3) ...(3)
在上式(2)中,W1是射血波S1的脉搏波振幅、即图4的与射血波S1的基准时间T1相对应的脉搏波S的振幅W1,该振幅W1与射血波S1的脉搏波压力相对应。此外,W3是反射波S2的脉搏波振幅、即图4的与反射波S2的基准时间T2相对应的反射波S2的振幅W3,该振幅W3与反射波S2的脉搏波压力相对应。
另一方面,如图5所示,设从人体的心脏51到脉搏波测定部位52的距离为d(m),从心脏51到反射点53为止的距离为h(m),从心脏51的跳动时到射血波S1的基准时间T1为止的时间为ΔT1(秒),射血波S1的基准时间T1与反射波S2的基准时间T2的时间差(T2-T1)为ΔT2(秒),则能够得到下式(4)、(5)。
(PWV1)2=(d/ΔT1)2 ...(4)
(PWV2)2=((2h+d)/(ΔT1+ΔT2))2 ...(5)
根据上式(2)和上式(4)能够得到下式(6),根据上式(3)和上式(5)能够得到下式(7)。
α·exp(β×W1)=(d/ΔT1)2 ...(6)
α·exp(β×W3)=((2h+d)/(ΔT1+ΔT2))2 ...(7)
根据上式(6)和上式(7),能够消去时间ΔT1。即,根据上式(7)能够得到下式(8)。
(ΔT1+ΔT2)2=(2h+d)2/α·exp(β×W3)
ΔT1+ΔT2=(2h+d)/{α·exp(β×W3)}1/2
ΔT1=(2h+d)/{α·exp(β×W3)}1/2-ΔT2 ...(8)
将该式(8)代入到上式(4)中,能够得到下式(9)。
(PWV1)2=d2×[(2h+d)/{α·exp(β×W3)}1/2-ΔT2]-2
PWV1=d×[(2h+d)/α·exp(β×W3)}1/2-ΔT2]-1...(9)
即,能够根据已知的常数α、β、已知的测定值即距离d、h、由基准时间检测部2求出的基准时间T2与T1的差(T2-T1)=ΔT2、和由射血波分量除去部4求出的与反射波S2的基准时间T2相对应的反射波S2的振幅W3,来算出射血波S1的脉搏波传播速度PWV1。
即,上述脉搏波传播速度检测部5,通过基于上述那样的式(4)~(7)而导出的上式(9),基于上述射血波分量S1、反射分量S2的基准时间T1、T2的差ΔT2、和与上述反射波分量S2的基准时间T2相对应的脉搏波S的射血波分量S1、反射分量S2的振幅W1、W3,来求出上述脉搏波S的射血波S1的传播速度PWV1。因此,能够考虑到射血波分量S1的振幅W1与反射分量S2的振幅W3的差异所导致的射血波与反射波的脉搏波传播速度的差异来高精度地对脉搏波传播速度进行测定。
此外,在上述实施方式中,说明了脉搏波振幅检测部3具有射血波分量除去部4的情况,但在上述脉搏波振幅检测部3不具有射血波分量除去部4的情况下,使用在上式(3)中采用了与反射波S2的基准时间T2相对应的脉搏波S的振幅W2来取代与反射波S2的基准时间T2相对应的反射波S2的振幅W3的下式(10)。
(PWV2)2=α·exp(β×W2) ...(10)
在此情况下,根据上式(10)和上式(5)能够得到下式(11)。
α·exp(β×W2)=((2h+d)/(ΔT1+ΔT2))2 ...(11)
因此,根据前述的式(6)和上式(11),能够消去时间ΔT1,从而得到下式(12)。
ΔT1=(2h+d)/{α·exp(β×W2)}1/2-ΔT2 ...(12)
将该式(12)代入到前述的式(4)中,能够得到下式(13)。
(PWV1)2=d2×[(2h+d)/{α·exp(β×W2)}1/2-ΔT2]-2
PWV1=d×[(2h+d)/α·exp(β×W2)}1/2-ΔT2]-1...(13)
即,能够根据已知的常数α、β、已知的测定值即距离d、h、由基准时间检测部2求出的基准时间T2与T1的差(T2-T1)=ΔT2、和由脉搏波振幅检测部3求出的与反射波S2的基准时间T2相对应的脉搏波S的振幅W2,来算出射血波S1的脉搏波传播速度PWV1。
即,上述脉搏波传播速度检测部5,通过如上述那样导出的上式(13),基于上述射血波分量S1、反射分量S2的基准时间T1、T2的差ΔT2、和与上述反射波分量S2的基准时间T2相对应的脉搏波S的基准时间T2的振幅W2,来求出上述脉搏波S的射血波S1的传播速度PWV1。因此,能够考虑到脉搏波S的射血波S1与反射波S2的脉搏波传播速度的差异来高精度地对脉搏波传播速度进行测定。
另外,在上述实施方式中,对测定脉搏波传播速度的装置进行了说明,但从上式(1)可知,也能够根据在上述实施方式中测定出的脉搏波传播速度来测定血压。即,在利用者使用上述脉搏波传播速度测定装置的情况下,也可以显示为血压值来取代直接显示由该装置得到的脉搏波传播速度。这是因为,上述脉搏波传播速度对于医疗工作者来说是非常熟悉的生物体指标,而对于非医疗工作者来说,比起脉搏波传播速度,血压是更容易熟悉的生物体指标。
使用了袖带的一般的血压计,花费数十秒来将一组最高血压·最低血压·平均血压·脉搏数等提示给利用者,但根据基于本发明的实施方式的血压计,能够按照每1拍来检测血压,因此有可能更详细地把握生物体的状态。此外,通过取得数拍单位的加法平均·移动平均等,能够更高精度地检测脉搏波传播速度、血压。
此外,如上所述,也可以通过脉搏波传播速度测定程序来使计算机执行如下功能:由上述基准时间检测部2求出上述基准时间T1和T2的基准时间导出功能;由上述脉搏波振幅检测部3求出与上述基准时间T1、T2相对应的脉搏波S的振幅W1、W2的脉搏波振幅导出功能;由上述射血波分量除去部4求出与上述反射波分量S2的基准时间T2相对应的上述反射波分量S2的振幅W3的射血波分量除去功能;和由上述脉搏波传播速度检测部5,基于上述射血波分量S1、反射分量S2的基准时间T1、T2的差ΔT2、和与上述反射波分量S2的基准时间T2相对应的脉搏波S的反射分量S2的振幅W3,来求出上述脉搏波S的射血波S1的传播速度PWV1的功能。
此外,如上所述,也可以通过脉搏波传播速度测定程序来使计算机执行如下功能:由上述基准时间检测部2求出上述基准时间T1和T2的基准时间导出功能;由上述脉搏波振幅检测部3求出与上述基准时间T1、T2相对应的脉搏波S的振幅W1、W2的脉搏波振幅导出功能;和由上述脉搏波传播速度检测部5,基于上述射血波分量S1、反射分量S2的基准时间T1、T2的差ΔT2、和与上述反射波分量S2的基准时间T2相对应的脉搏波S的振幅W2,来求出上述脉搏波S的射血波S1的传播速度PWV1的功能。
符号说明
1脉搏波检测部
2基准时间检测部
3脉搏波振幅检测部
4射血波分量除去部
5脉搏波传播速度检测部
6射血波·反射波信息提取部
10脉搏波传播速度测定装置
Q1极大点
Q24次微分波的第3零交叉点
S脉搏波
S1射血波
S2反射波
T1射血波分量的基准时间
T2反射波分量的基准时间
W1与基准时间T1相对应的脉搏波S的振幅
W2与基准时间T2相对应的脉搏波S的振幅
W3与基准时间T2相对应的反射波S2的振幅
PWV1射血波S1的传播速度
PWV2反射波S2的传播速度
Claims (3)
1.一种脉搏波传播速度测定装置,具备:
脉搏波检测部(1),其检测生物体的某一个部位的脉搏波;
基准时间检测部(2),其检测用于确定由上述脉搏波检测部(1)检测出的上述一个部位的脉搏波中包含的射血波分量的基准时间和用于确定上述脉搏波中包含的反射波分量的基准时间;
脉搏波振幅检测部(3),其检测与由上述基准时间检测部(2)检测出的上述射血波分量的基准时间相对应的上述脉搏波的振幅,并且检测与由上述基准时间检测部检测出的上述反射波分量的基准时间相对应的上述脉搏波的振幅;和
传播速度检测部(5),其求出上述脉搏波的传播速度,
该脉搏波传播速度测定装置的特征在于,
上述传播速度检测部(5)基于由上述基准时间检测部(2)检测出的上述射血波分量的基准时间以及上述反射波分量的基准时间、和由上述脉搏波振幅检测部(3)检测出的与上述射血波分量的基准时间相对应的脉搏波的振幅以及与上述反射波分量的基准时间相对应的脉搏波的振幅,根据下式(1)来求出上述脉搏波的传播速度,
PWV1=d×[(2h+d)/α·exp(β×W2)}1/2-ΔT2]-1......(1)
其中,在上式中,PWV1表示射血波的脉搏波传播速度,α、β表示已知的常数,d表示从人体的心脏到脉搏波测定部位的距离,h表示从心脏到反射点为止的距离,ΔT2表示由基准时间检测部(2)求出的、反射波分量的基准时间与射血波分量的基准时间的差,W2表示由脉搏波振幅检测部(3)求出的与反射波分量的基准时间相对应的脉搏波的振幅。
2.根据权利要求1所述的脉搏波传播速度测定装置,其特征在于,
上述脉搏波振幅检测部(3)具有射血波分量除去部(4),
该射血波分量除去部(4)求出与上述反射波分量的基准时间相对应的上述脉搏波的振幅中包含的射血波分量,从与上述反射波分量的基准时间相对应的上述脉搏波的振幅中除去上述脉搏波的振幅中包含的射血波分量,来求出与上述反射波分量的基准时间相对应的上述反射波分量的振幅,
上述传播速度检测部(5)基于由上述基准时间检测部(2)检测出的上述射血波分量的基准时间以及上述反射波分量的基准时间、和由上述脉搏波振幅检测部(3)检测出的与上述射血波分量的基准时间相对应的脉搏波的振幅、以及由上述射血波分量除去部(4)求出的上述反射波分量的振幅,代替上述式(1)而由下式(2)来求出上述脉搏波的传播速度,
PWV1=d×[(2h+d)/α·exp(β×W3)}1/2-ΔT2]-1...(2)
其中,上式中,W3表示由射血波分量除去部(4)求出的与反射波的基准时间相对应的反射波(S2)的振幅。
3.根据权利要求2所述的脉搏波传播速度测定装置,其特征在于,
上述传播速度检测部(5)基于由上述基准时间检测部(2)检测出的上述射血波分量的基准时间与上述反射波分量的基准时间的时间差、和由上述脉搏波振幅检测部(3)检测出的与上述射血波分量的基准时间相对应的脉搏波的振幅与由上述射血波分量除去部(4)求出的上述反射波分量的振幅的振幅差,来求出上述脉搏波的传播速度。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009257217A JP5328613B2 (ja) | 2009-11-10 | 2009-11-10 | 脈波伝播速度測定装置および脈波伝播速度測定プログラム |
JP2009-257217 | 2009-11-10 | ||
PCT/JP2010/069408 WO2011058900A1 (ja) | 2009-11-10 | 2010-11-01 | 脈波伝播速度測定装置および脈波伝播速度測定プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102711597A CN102711597A (zh) | 2012-10-03 |
CN102711597B true CN102711597B (zh) | 2015-09-30 |
Family
ID=43991559
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080050888.XA Expired - Fee Related CN102711597B (zh) | 2009-11-10 | 2010-11-01 | 脉搏波传播速度测定装置以及脉搏波传播速度测定程序 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20120226174A1 (zh) |
EP (1) | EP2499967B1 (zh) |
JP (1) | JP5328613B2 (zh) |
CN (1) | CN102711597B (zh) |
WO (1) | WO2011058900A1 (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635110A (zh) * | 2015-01-23 | 2015-05-20 | 昆明理工大学 | 一种基于实测数据的波速校核方法 |
CN104720777B (zh) * | 2015-03-10 | 2017-11-03 | 中国科学院电子学研究所 | 一种无创连续血压生理监测系统 |
KR102501837B1 (ko) * | 2015-11-09 | 2023-02-21 | 삼성전자주식회사 | 신호 특징 추출 방법 및 장치 |
GB2552455B8 (en) * | 2016-06-16 | 2022-06-15 | Digital & Future Tech Limited | Blood monitoring |
TWI644651B (zh) * | 2017-11-17 | 2018-12-21 | 國立成功大學 | 訊號同步處理裝置及具有訊號同步處理功能的聽診器、聽診資訊輸出系統與病徵診斷系統 |
DE102018209198B3 (de) * | 2018-06-08 | 2019-08-29 | Gert Küchler | Verfahren und Gerät zur Bestimmung mindestens eines physiologischen Parameters |
KR102576126B1 (ko) * | 2019-09-19 | 2023-09-06 | 삼성전자주식회사 | 생체정보 추정 장치 및 방법 |
CN110960199B (zh) * | 2019-12-24 | 2022-05-27 | 中国人民解放军陆军军医大学第一附属医院 | 一种双变量测量动脉硬化程度的系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1348373A1 (en) * | 2002-03-01 | 2003-10-01 | Colin Corporation | Arteriosclerosis inspecting apparatus |
EP1402814A1 (en) * | 2002-09-27 | 2004-03-31 | Colin Corporation | Arteriosclerosis-degree evaluating apparatus |
CN1504166A (zh) * | 2002-12-05 | 2004-06-16 | ŷķ������ҽ����ҵ��ʽ���� | 脉搏波测定装置 |
CN1909829A (zh) * | 2004-02-18 | 2007-02-07 | 伊里斯·米科诺斯 | 测量血液动力学参数的装置及方法 |
WO2007062456A1 (en) * | 2005-12-01 | 2007-06-07 | Atcor Medical Pty Ltd | A method of estimating pulse wave velocity |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3495348B2 (ja) * | 2001-07-02 | 2004-02-09 | 日本コーリン株式会社 | 脈波伝播速度情報測定装置 |
JP3533406B2 (ja) * | 2001-07-02 | 2004-05-31 | コーリンメディカルテクノロジー株式会社 | 動脈硬化評価装置 |
JP3621379B2 (ja) * | 2002-01-09 | 2005-02-16 | コーリンメディカルテクノロジー株式会社 | 動脈硬化評価装置 |
JP3643567B2 (ja) * | 2002-04-17 | 2005-04-27 | コーリンメディカルテクノロジー株式会社 | 振幅増加指数測定装置 |
JP2003305011A (ja) * | 2002-04-17 | 2003-10-28 | Nippon Colin Co Ltd | 振幅増加指数測定装置 |
DE602005012204D1 (de) * | 2004-02-18 | 2009-02-26 | Miklos Illyes | Gerät und verfahren zur messung von hämodynamischen parametern |
JP4581496B2 (ja) * | 2004-06-14 | 2010-11-17 | オムロンヘルスケア株式会社 | 脈波解析装置および脈波解析プログラム |
JP2006043146A (ja) * | 2004-08-04 | 2006-02-16 | Sony Corp | 脈波解析方法及び装置 |
JP4855721B2 (ja) * | 2005-06-29 | 2012-01-18 | フクダ電子株式会社 | 血圧測定装置 |
-
2009
- 2009-11-10 JP JP2009257217A patent/JP5328613B2/ja not_active Expired - Fee Related
-
2010
- 2010-11-01 EP EP10829860.5A patent/EP2499967B1/en not_active Not-in-force
- 2010-11-01 WO PCT/JP2010/069408 patent/WO2011058900A1/ja active Application Filing
- 2010-11-01 US US13/508,933 patent/US20120226174A1/en not_active Abandoned
- 2010-11-01 CN CN201080050888.XA patent/CN102711597B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1348373A1 (en) * | 2002-03-01 | 2003-10-01 | Colin Corporation | Arteriosclerosis inspecting apparatus |
EP1402814A1 (en) * | 2002-09-27 | 2004-03-31 | Colin Corporation | Arteriosclerosis-degree evaluating apparatus |
CN1504166A (zh) * | 2002-12-05 | 2004-06-16 | ŷķ������ҽ����ҵ��ʽ���� | 脉搏波测定装置 |
CN1909829A (zh) * | 2004-02-18 | 2007-02-07 | 伊里斯·米科诺斯 | 测量血液动力学参数的装置及方法 |
WO2007062456A1 (en) * | 2005-12-01 | 2007-06-07 | Atcor Medical Pty Ltd | A method of estimating pulse wave velocity |
Also Published As
Publication number | Publication date |
---|---|
JP5328613B2 (ja) | 2013-10-30 |
WO2011058900A1 (ja) | 2011-05-19 |
EP2499967B1 (en) | 2014-06-25 |
JP2011101696A (ja) | 2011-05-26 |
US20120226174A1 (en) | 2012-09-06 |
CN102711597A (zh) | 2012-10-03 |
EP2499967A4 (en) | 2013-08-21 |
EP2499967A1 (en) | 2012-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102711597B (zh) | 脉搏波传播速度测定装置以及脉搏波传播速度测定程序 | |
CN102843962B (zh) | 脉搏波传播速度测定装置、脉搏波传播速度的测定方法以及脉搏波传播速度的测定程序 | |
US7137955B2 (en) | Methods and systems for distal recording of phonocardiographic signals | |
EP2260886B1 (en) | Blood pressure estimation apparatus and blood pressure estimation method | |
KR102346873B1 (ko) | 피검자의 맥파를 측정하기 위한 시스템 및 방법 | |
US20140148702A1 (en) | Method, device and system for determining the open/closed switch moment of an artery of interest under a changing pressure | |
US20210369235A1 (en) | Method for determining pulse transmission time, arteriosclerosis detection apparatus and system | |
US10925516B2 (en) | Method and apparatus for estimating the aortic pulse transit time from time intervals measured between fiducial points of the ballistocardiogram | |
JP2012070877A (ja) | 血圧情報測定装置および血圧情報測定方法 | |
JP2011212364A (ja) | 心音測定装置 | |
CN107049270A (zh) | 用于测量脉搏波传导时间的系统及方法 | |
Rasool et al. | Continuous and noninvasive blood pressure estimation by two-sensor measurement of pulse transit time | |
JP6158324B2 (ja) | 血管異常検知装置 | |
JP5681434B2 (ja) | 脈波伝播速度測定装置および脈波伝播速度測定プログラム | |
Stone et al. | Acute changes in carotid-femoral pulse-wave velocity are tracked by heart-femoral pulse-wave velocity | |
CN104042200A (zh) | 一种逐拍动脉血压的无创监测装置及其方法 | |
Li et al. | Initial study on pulse wave velocity acquired from one hand using two synchronized wireless reflectance pulse oximeters | |
CN102413760A (zh) | 监控外周退耦 | |
JP5690550B2 (ja) | 脈波伝播速度測定装置 | |
CN101332084B (zh) | 脉搏波的面积重心位移分析方法 | |
Manoj et al. | Evaluation of Nonlinear Wave Separation Method to Assess Reflection Transit Time: A Virtual Patient Study | |
CN104602597B (zh) | 循环器官功能运算装置 | |
Spigulis et al. | Micro-circulation of skin blood: optical monitoring by advanced photoplethysmography techniques | |
JP2017158610A (ja) | 血圧測定装置および血圧測定方法 | |
Loza et al. | Measuring Regional Wave Speed using a Low Cost Portable Device |
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: 20150930 Termination date: 20201101 |