CN1251027A - 脉波诊断装置、血压监视装置、脉波形状监视装置和药理作用监视装置 - Google Patents

脉波诊断装置、血压监视装置、脉波形状监视装置和药理作用监视装置 Download PDF

Info

Publication number
CN1251027A
CN1251027A CN98803516A CN98803516A CN1251027A CN 1251027 A CN1251027 A CN 1251027A CN 98803516 A CN98803516 A CN 98803516A CN 98803516 A CN98803516 A CN 98803516A CN 1251027 A CN1251027 A CN 1251027A
Authority
CN
China
Prior art keywords
pulse
mentioned
pressure
waveform
blood pressure
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
Application number
CN98803516A
Other languages
English (en)
Other versions
CN1326489C (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.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
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 Seiko Epson Corp filed Critical Seiko Epson Corp
Publication of CN1251027A publication Critical patent/CN1251027A/zh
Application granted granted Critical
Publication of CN1326489C publication Critical patent/CN1326489C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4854Diagnosis based on concepts of traditional oriental medicine
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/021Measuring pressure in heart or blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02116Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave amplitude
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802Sensor mounted on worn items
    • A61B5/681Wristwatch-type devices
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6814Head
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6814Head
    • A61B5/6815Ear
    • A61B5/6816Ear lobe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6822Neck
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6825Hand
    • A61B5/6826Finger
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6828Leg
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/683Means for maintaining contact with the body
    • A61B5/6838Clamps or clips
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7207Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • A61B5/721Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/726Details of waveform analysis characterised by using transforms using Wavelet transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7239Details of waveform analysis using differentiation including higher order derivatives
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Vascular Medicine (AREA)
  • Otolaryngology (AREA)
  • Alternative & Traditional Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

当FFT处理部(40)对除去了体动的脉波波形MHj进行频谱分析并生成脉波分析数据MKD时,潮波特征抽出部(50)和重博波特征抽出部(60)分别生成表示潮波和重博波特征的潮波特征数据(TWD)和重博波特征数据(DWD)。然后,脉象判定部(70)根据这些数据(TWD、DWD)生成脉象数据(ZD),告状部(80)告知被测者的脉象。

Description

脉波诊断装置、血压监视装置、 脉波形状监视装置和药理作用监视装置
【技术领域】
本发明涉及很适合用来特定脉搏的类别的脉波诊断装置、将平均血压和脉压作为参数使用的血压监视装置、使用了与动脉血压波形的重博部有关的参数的脉波波形监视装置和药理作用监视装置。
【背景技术】
一般来说,脉波是从心脏发出的在血管内传送的血液的波。因此可知,通过检测脉波并进行分析可以得到各种医学信息。而且,随着脉波研究的深入人们发现,通过使用各种各样的方法对从人体采取的脉波进行分析,可以得到单从血压和心博数不能得到的各种各样的信息,能够根据这些信息进行诊断。
本发明者在PCT/JP96/01254(发明的名称:身体状态的诊断装置和控制装置)中,着眼于脉波波形的形状和其畸变率的关系,检测被测者的脉波波形并进行处理,由此算出该脉波波形的畸变率,根据该畸变率去特定脉波波形的形状,从而可以诊断该被测者的身体状态。
在此,就上述发明申请中叙述的脉波波形的形状和其畸变率的关系简单地进行说明。
首先,脉波波形有各种各样的分类,其形状也是各式各样,这里,说明按照东方传统医学之一的中国医学的分类方法所得到的具有代表性的脉波波形的形状。图45A~图45C是表示按照该分类方法所得到的具有代表性的脉波波形的形状的图。
图45A示出的脉波波形称为‘平脉’,是正常健康人的脉象。如图所示,该‘平脉’的特征是舒畅而缓和,节奏稳定而不乱。
其次,图45B示出的脉波波形称为‘滑脉’,是血流状态异常者的脉象,其特征是急剧上升急剧下降,大动脉重博深,其后的波峰也比正常值高出很多。该‘滑脉’认为是由于浮肿、肝肾疾病、呼吸道疾病、肠胃疾病、炎症等疾患引起的,使脉的往来非常流利、圆滑。
此外,图45B示出的脉波波形称为‘弦脉’,是血管壁紧张度上升者的脉象,其特征是急剧上升后不立刻下降而持续一定时间的高压状态。该‘弦脉’认为是肝胆疾病、皮肤病、高血压、疼痛性疾病等疾患的反映,原因是因自律神经系统紧张而使血管壁紧张和弹性减小,压出来的血液的有节奏的运动的影响难以表现出来。
再有,在图45A~45C中,纵轴和横轴分别代表血压(mmHg)和时间(秒)。
而且,这样的脉波波形的脉象和其畸变率d之间具有图46所示那样的关系。这里,脉波波形的畸变率d由下式(1)决定。 d = A 2 2 + A 3 2 + … … + A n 2 A 1 - - - - ( 1 )
再有,在该式(1)中,A1是脉波的基波成分的振幅,A2、A3、……An别是脉波的第2、第3、第n次谐波成分的振幅。
因此,若将被测者的脉波波形检测出来,进行例如FFT(高速付里叶变换)处理,分别求出振幅A1~An,再算出畸变率d,就可以利用图46所示的相关关系定量地特定出脉波波形的脉象。
如图46所示,当判定是滑脉时,其畸变率d在0.98~1.22的范围内,当判定是平脉时,其畸变率d在0.92~1.10的范围内,当判定是弦脉时,其畸变率d在0.73~0.94的范围内。
这时,当脉波波形的畸变率d在0.98~1.10范围内时,既可以判定为滑脉,也可以判定为平脉。此外,当脉波波形的畸变率d在0.92~0.94范围内时,既可以判定为平脉,也可以判定为弦脉。因此,先有的脉波诊断装置难以确切地判定脉象。
此外,在非侵袭血压监视中,使用测定最高血压和最低血压并进行显示的血压计。
但是,即使最高血压和最低血压相同,也存在各种各样的血压波形。因此,存在只用最高血压和最低血压不能充分表示各人的血压特征的问题。
此外,为了了解各人的血压状态,平均血压是很重要的参数,但若只测定最高血压和最低血压,还不能知道平均血压。
另一方面,中国医学和印度传统医学的诊脉是根据诊察者用手指适当按压患者手腕远位部的桡骨动脉所感知的脉波形状、即诊察者手指所感知的与桡骨动脉的血压变化对应的按压力的变化情况来进行诊察的。
例如,在中国医学中,对桡骨动脉适当施加按压力所感知的脉波形状大致分为3类,如前所述,分别称为平脉、滑脉和弦脉。平脉的感觉是舒畅而缓和节奏稳定而不乱,是健康人的脉象。滑脉的感觉是圆滑,来回流动,表示血流状态异常。弦脉的感觉是抻紧拉长的脉,起因于血管壁的紧张老化。
但是,这样的根据脉波形状的诊断法,因脉波形状的分类依赖诊察者的感觉故存在客观性和再现性问题。
本发明是鉴于上述问题点提出的,其目的在于提供一种脉波诊断装置,能够客观而正确地判定脉象。
本发明的另一个目的在于,提供一种血压监视装置,能够根据最高血压和最低血压的信息更详细地显示血压的状态,能够非侵袭地监视与血压有关的参数。
本发明的再一个目的在于,提供一种脉波波形监视装置,能够根据脉波形状进行客观的诊察,而且具有再现性。
【发明的公开】
(1)本发明的脉波诊断装置包括:
检测身体的脉波波形的脉波检测装置;
从上述脉波波形中抽出潮波的特征并生成潮波特征信息的潮波特征抽出装置;
从上述脉波波形中抽出重博波的特征并生成重博波特征信息的重博波特征抽出装置;
根据上述潮波特征信息和上述重博波特征信息判定上述身体的脉象的脉象判定装置。
(2)在第(1)条记载的脉波诊断装置中,上述潮波特征抽出装置最好根据上述潮波在时间域中的振幅变化生成潮波特征信息,上述重博波特征抽出装置最好根据上述重博波在时间域中的振幅变化生成重博波特征信息。
(3)在第(2)条记载的脉波诊断装置中,上述潮波和上述重博波在时间域中的振幅变化最好通过将上述脉波波形对时间进行一次微分或2次微分算出来。
(4)本发明的脉波诊断装置包括:
检测身体的脉波波形的脉波检测装置;
对上述脉波波形进行频谱分析的频谱分析装置;
根据上述频谱分析装置的分析结果从上述脉波波形中抽出潮波的特征并生成潮波特征信息的潮波特征抽出装置;
根据上述频谱分析装置的分析结果从上述脉波波形中抽出重博波的特征并生成重博波特征信息的重博波特征抽出装置;
根据上述潮波特征信息和上述重博波特征信息判定上述身体的脉象的脉象判定装置。
(5)在第(4)条记载的脉波诊断装置中,
上述潮波特征抽出装置最好先特定上述脉波波形中的潮波的持续时间,再根据该期间内的上述频谱分析装置的分析结果从上述脉波波形中抽出潮波的特征,并生成潮波特征信息,
上述重博波特征抽出装置最好先特定上述脉波波形中的重博波的持续时间,再根据该期间内的上述频谱分析装置的分析结果从上述脉波波形中抽出重博波的特征,并生成重博波特征信息。
(6)在第(4)或第(5)条记载的脉波诊断装置中,
上述频谱分析装置最好对上述脉波波形进行FFT处理。
(7)在第(4)或第(5)条记载的脉波诊断装置中,
上述频谱分析装置最好对上述脉波波形进行小波变换处理。
(8)在第(1)至第(7)条记载的脉波诊断装置中,
最好具有告知装置,告诉利用上述脉象判定装置判定的脉象。
(9)本发明的脉波诊断装置包括:
从身体的检测部位检测脉波波形的脉波检测装置;
自相关运算装置,计算表示利用上述脉波检测装置检测出的上述脉波波形的自相关的自相关数据。
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
(10)在第(9)条记载的脉波诊断装置中,
上述脉象数据生成装置最好通过将上述自相关数据与预先设定的阈值进行比较来生成上述脉象数据。
(11)在第(10)条记载的脉波诊断装置中,
上述脉象数据生成装置最好具有最小值检测部和比较部,最小值检测部检测一个心博期间上述自相关数据的最小值,
比较部将上述最小值检测部检测出的最小值与上述阈值进行比较并生成上述脉象数据。
(12)在第(10)条记载的脉波诊断装置中,
上述脉象数据生成装置最好具有最小值检测部和比较部,最小值检测部对在几个心博期间检测出上述自相关数据的各最小值进行平均后检测其平均最小值,
比较部将上述最小值检测部检测出的平均最小值与上述阈值进行比较并生成上述脉象数据。
(13)在第(9)条记载的脉波诊断装置中,
上述脉象数据生成装置最好具有时间测量部、运算部和比较部,时间测量部将上述自相关数据与预定的阈值进行比较,测量上述自相关数据高于上述阈值的时间间隔和低于上述阈值的时间间隔,
运算部计算由上述时间测量部测量的上述时间间隔在一个心博期间所占的比例,
比较部通过将上述运算部计算的结果与预定的阈值进行比较来生成上述脉象数据。
(14)在第(13)条记载的脉波诊断装置中,
上述运算部最好算出由上述时间测量部测量的上述时间间隔在一个心博期间所占的比例,并计算该算出结果的平均值。
(15)在第(9)条记载的脉波诊断装置中,
上述脉象数据生成装置最好具有变化率运算部和变化率比较部,变化率运算部根据上述自相关数据检测上述自相关数据的变化率,
变化率比较部通过将上述变化率运算部检测出的变化率与预定的阈值进行比较来生成上述脉象数据。
(16)在第(15)条记载的脉波诊断装置中,
上述变化率比较部最好检测上述变化率的最大值,通过将变化率最大值与上述阈值进行比较来生成上述脉象数据。
(17)在第(9)条记载的脉波诊断装置中,
上述脉象数据生成装置最好具有最小值检测部、第1比较部、时间测量部、运算部和第2比较部,最小值检测部检测1个心博期间上述自相关数据的最小值,
第1比较部将上述最小值检测部检测出的最小值与预定的第1阈值进行比较,当上述最小值低于上述阈值时,生成表示弦脉的脉象数据,
时间测量部将上述自相关数据与预定的第2阈值进行比较,测量上述自相关数据高于上述第2阈值的时间间隔或低于上述第2阈值的时间间隔,
运算部计算由上述时间测量部测量的上述时间间隔在一个心博期间所占的比例,
第2比较部通过将上述运算部计算的结果与预定的第3阈值进行比较来生成表示平脉或滑脉的上述脉象数据。
(18)在第(9)条记载的脉波诊断装置中,
上述脉象数据生成装置最好具有最小值检测部、第1比较部、变化率运算部和第2比较部,最小值检测部检测1个心博期间上述自相关数据的最小值,
第1比较部将上述最小值检测部检测出的最小值与上述阈值进行比较,当上述最小值低于上述阈值时,生成表示弦脉的脉象数据,
变化率运算部根据自相关数据检测自相关数据的变化率,
第2比较部通过将上述变化率运算部检测出的变化率与预定的阈值进行比较来生成表示平脉或滑脉的上述脉象数据。
(19)在第(17)或第(18)条记载的脉波诊断装置中,
上述自相关数据表示自相关系数,上述第1比较部的比较动作所用的上述第1阈值最好大约是0.25。
(20)在第(17)条记载的脉波诊断装置中,
上述自相关数据表示自相关系数,上述时间测量部的比较动作所用的上述第2阈值最好设定在0.4到0.8的范围内。
(21)在第(9)至第(20)任何一条记载的脉波诊断装置中,
最好具有检测表示上述身体活动的体动波形的体动检测装置和体动成分除去装置,该体动成分除去装置根据上述体动波形生成上述脉波波形中的体动成分,从上述脉波波形中除去上述体动成分并生成体动除去脉波波形,
上述自相关运算装置不是根据上述脉波波形而是根据上述体动除去脉波波形计算表示自相关的自相关数据。
(22)在第(21)条记载的脉波诊断装置中,
最好具有判定装置,根据上述体动检测装置检测出的上述体动波形判定上述身体有无活动,
上述体动成分除去装置在上述判定装置的判定结果是没有体动时停止体动除去动作,并输出上述脉波波形,而不输出上述体动除去脉波波形。
(23)在第(9)至第(20)任何一条记载的脉波诊断装置中,
最好具有第1小波变换装置、体动检测装置、第2小波变换装置、体动成分除去装置和逆小波变换装置,
第1小波变换装置对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动检测装置检测上述身体的活动并输出体动波形,
第2小波变换装置对上述体动检测装置检测出的上述体动波形进行小波变换并生成各频率区的体动分析数据,
体动成分除去装置从上述脉波分析数据中减去上述体动分析数据,并生成已除去了体动的体动除去脉波分析数据,
逆小波变换装置对上述体动除去脉波分析数据进行逆小波变换并生成体动除去脉波波形。
上述自相关运算装置不是根据上述脉波波形而是根据上述体动除去脉波波形计算表示自相关的自相关数据。
(24)在第(9)至第(20)任何一条记载的脉波诊断装置中,
最好具有小波变换装置、体动成分除去装置和逆小波变换装置,
小波变换装置对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动成分除去装置从上述脉波分析数据中除去与预定的体动对应的频率成分,并生成体动除去脉波分析数据,
逆小波变换装置对上述体动除去脉波分析数据进行逆小波变换并生成体动除去脉波波形。
上述自相关运算装置不是根据上述脉波波形而是根据上述体动除去脉波波形计算表示自相关的自相关数据。
(25)本发明的脉波诊断装置包括:
从身体的检测部位检测脉波波形的脉波检测装置;
小波变换装置,对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
自相关运算装置,对某频率区的上述脉波分析数据计算表示自相关的自相关数据,
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
(26)本发明的脉波诊断装置包括:
从身体的检测部位检测脉波波形的脉波检测装置;
第1小波变换装置,对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动检测装置,检测上述身体的活动并输出体动波形,
第2小波变换装置,对上述体动检测装置检测出的上述体动波形进行小波变换并生成各频率区的体动分析数据,
体动成分除去装置,从上述脉波分析数据中减去上述体动分析数据,并生成已除去了体动的体动除去脉波分析数据,
自相关运算装置,对某频率区的上述体动除去脉波分析数据计算表示自相关的自相关数据,
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
(27)本发明的脉波诊断装置包括:
从身体的检测部位检测脉波波形的脉波检测装置;
小波变换装置,对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动成分除去装置,从上述脉波分析数据中除去与预定的体动对应的频率成分,并生成体动除去脉波分析数据,
自相关运算装置,对某频率区的上述体动除去脉波分析数据计算表示自相关的自相关数据,
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
(28)在第(9)至第(27)任意一条记载的脉波诊断装置中,
最好具有告知装置,告诉利用上述脉象数据生成装置生成的上述脉象数据。
(29)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压。
本发明的血压监视装置具有根据动脉压波形检测部检测出的动脉压波形计算平均血压的平均血压算出部。因此,使用动脉压波形检测部检测出的动脉压波形,能够监视平均血压。
(30)第(29)条记载的血压监视装置,
最好进而具有脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压。
若按照本发明,使用动脉压波形检测部检测出的动脉压波形,能够监视脉压。
(31)第(30)条记载的血压监视装置,
最好进而具有血压变换部,将上述动脉压波形检测部检测出的动脉压波形变换成与心脏的高度对应的位置上的动脉压波形、即心脏位动脉压波形,上述平均血压算出部根据上述心脏位动脉压波形算出平均血压,上述脉压算出部根据上述心脏位动脉压波形算出脉压。
若按照本发明,血压变换部将动脉压波形检测部检测出的动脉压波形变换成与心脏的高度对应的位置上的动脉压波形、即心脏位动脉压波形。而且,根据心脏位动脉压波形,平均血压算出部算出平均血压,脉压算出部算出脉压。因此,使用动脉压波形检测部检测出的动脉压波形,至少能够监视与心脏的高度对应的动脉的平均血压和脉压当中的一个参数。
(32)第(30)或第(31)条记载的血压监视装置,
最好进而具有预先存储血压判定信息的血压判定信息存储部和血压判定部,血压判定部根据上述平均血压和上述脉压当中的至少一个参数以及上述血压判定信息进行血压判定。
若按照本发明,血压监视装置根据平均血压和脉压当中的至少一个参数以及预先存储的血压判定信息能够进行例如高血压、低血压和正常等的判定。
(33)第(30)至第(32)任何一条记载的血压监视装置,
最好进而具有输出部,输出与上述平均血压对应的信息、与上述脉压对应的信息和与上述血压判定对应的信息当中的至少一种信息。
若按照本发明,血压监视装置利用输出部能够输出与平均血压对应的信息、与脉压对应的信息和与上述血压判定对应的信息当中的至少一种信息,例如以数值、图形或电压等形式输出。
(34)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形计算最高血压和最低血压的压差、即脉压。
若按照本发明,使用动脉压波形检测部检测出的动脉压波形,能够监视脉压。
(35)(34)记载的血压监视装置,
最好具有根据上述脉压算出最高血压的最高血压算出部和根据上述脉压及上述最高血压算出最低血压的最低血压算出部。
若按照本发明,利用最高血压可以表示成脉压的线性函数的特点,最高血压算出部根据脉压决定最高血压。而且,脉压是最高血压和最低血压的压差,因此,若知道最高血压和脉压,最低血压算出部就能够求出最低血压。
(36)本发明的脉波形状监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
若按照本发明,脉波形状监视装置根据动脉压波形检测部连续测定动脉血压得到的动脉压波形,能够算出重博波高。
(37)第(36)条记载的脉波形状监视装置,
最好进而具有重博压差比算出部,计算从上述动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
若按照本发明,脉波形状监视装置根据动脉压波形检测部连续测定动脉血压得到的动脉压波形,能够算出重博波高和重博压差比。
(38)第(37)条记载的脉波形状监视装置,
最好进而具有平均血压脉压比算出部,计算从上述动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
若按照本发明,脉波形状监视装置根据动脉压波形检测部连续测定动脉血压得到的动脉压波形,能够算出重博波高、重博压差比和平均血压脉压比。
(39)第(38)条记载的脉波形状监视装置,
最好进而具有血压变换部,将上述动脉压波形检测部检测出的动脉压波形变换成与心脏的高度对应的位置上的动脉压波形、即心脏位动脉压波形,
上述重博波高算出部计算从上述心脏位动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
上述重博压差比算出部计算从上述心脏位动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
上述平均血压脉压比算出部计算从上述心脏位动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
若按照本发明,血压变换部将动脉压波形检测部检测出的动脉压波形变换成与心脏的高度对应的位置上的动脉压波形、即心脏位动脉压波形。而且根据心脏位动脉压波形,重博波高算出部算出重博波高,重博压差比算出部算出重博压差比,平均血压脉压比算出部算出平均血压脉压比。因此,使用动脉压波形检测部检测出的动脉压波形,至少能够监视与心脏的高度对应的动脉的重博波高、重博压差比和平均血压脉压比当中的一个参数。
(40)第(36)条记载的脉波形状监视装置,
最好进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和脉波形状判定部,脉波形状判定部根据上述重博波高和上述脉波形状判定信息判定脉波形状。
若按照本发明,根据重博波高和脉波形状判定信息,脉波形状判定部能够判定脉波形状。
(41)第(37)条记载的脉波形状监视装置,
最好进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和脉波形状判定部,
脉波形状判定部根据上述重博波高、上述重博压差比和上述脉波形状判定信息判定脉波形状。
若按照本发明,根据重博波高、重博压差比和脉波形状判定信息,脉波形状判定部能够判定脉波形状。
(42)第(38)条记载的脉波形状监视装置,
最好进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和脉波形状判定部,
脉波形状判定部根据上述重博波高、上述重博压差比、上述平均血压脉压比和上述脉波形状判定信息判定脉波形状。
若按照本发明,根据重博波高、重博压差比、平均血压脉压比和脉波形状判定信息,脉波形状判定部能够判定脉波形状。
(43)第(36)条记载的脉波形状监视装置,
最好进而具有输出部,输出与上述重博波高对应的信息和上述脉波形状对应的信息当中的至少一种信息。
若按照本发明,脉波形状监视装置从输出部能够输出与重博波高对应的信息、与脉波形状对应的信息当中的至少一种信息,例如以数值、图形或电压等形式输出。
(44)第(37)条记载的脉波形状监视装置,
最好进而具有输出部,输出与上述重博波高对应的信息、与上述重博压差比对应的信息和与上述脉波形状对应的信息当中的至少一种信息。
若按照本发明,脉波形状监视装置从输出部能够输出与重博波高对应的信息、与重博压差比对应的信息和与脉波形状对应的信息当中的至少一种信息,例如以数值、图形或电压等形式输出。
(45)第(38)条记载的脉波形状监视装置,
最好进而具有输出部,输出与上述重博波高对应的信息、与上述重博压差比对应的信息、与上述平均血压脉压比对应的信息和与上述脉波形状对应的信息当中的至少一种信息。
若按照本发明,脉波形状监视装置从输出部能够输出与重博波高对应的信息、与重博压差比对应的信息、与平均血压脉压比对应的信息和与脉波形状对应的信息当中的至少一种信息,例如以数值、图形或电压等形式输出。
(46)本发明的脉波形状监视装置具有:
脉波检测部,检测身体的脉波波形;
重博压差比算出部,计算从上述脉波波形得到的重博压和最低压的压差、即重博压差与最高压和最低压的压差、即脉压的比值、即重博压差比。
本发明与(36)记载的发明不同,重博压差比算出部因只需要求出比值故不需要绝对血压值。因此,可以不使用(36)记载的发明中的动脉压波形检测部,而使用脉波检测部,只需要检测出与其对应的波形(脉波)即可。
(47)本发明的脉波形状监视装置具有:
脉波检测部,检测身体的脉波波形;
平均压脉压比算出部,计算从上述脉波波形得到的平均压与最高压和最低压的压差、即脉压的比值、即平均压脉压比。
本发明与(16)记载的发明也不同,平均压脉压比算出部因只需要求出比值故不需要绝对血压值。因此,可以不使用(16)记载的发明中的动脉压波形检测部,而使用脉波检测部,只需要检测出与其对应的波形(脉波)即可。
(48)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压;
重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
(49)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压;
重博压差比算出部,计算从上述动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
(50)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压;
平均血压脉压比算出部,计算从上述动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
(51)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压;
重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
(52)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压;
重博压差比算出部,计算从上述动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
(53)本发明的血压监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压;
平均血压脉压比算出部,计算从上述动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
(54)本发明的脉波形状监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
驱出扩张压算出部,计算从上述动脉压波形得到的收缩期血压和重博之间的血压差、即驱出扩张压。
(55)第(54)条记载的脉波形状监视装置,
最好进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和根据上述驱出扩张压判定脉波形状的脉波形状判定部。
(56)本发明的脉波形状监视装置具有:
脉波检测部,检测身体的脉波波形;
驱出扩张压比算出部,计算从上述动脉压波形得到的重博压和收缩压之间的压差、即驱出扩张压与最高压和最低压之间的压差、即脉压的比。
(57)第(56)条记载的脉波形状监视装置,
最好进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和根据上述驱出扩张压比判定脉波形状的脉波形状判定部。
(58)本发明的药理作用监视装置具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
驱出扩张压算出部,计算从上述动脉压波形得到的收缩期血压和重博之间的血压差、即驱出扩张压。
(59)第(58)条记载的药理作用监视装置,
最好进而具有重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
(60)本发明的药理作用监视装置具有:
脉波检测部,从身体检测脉波波形;
驱出扩张压比算出部,计算从上述脉波波形得到的重博压和收缩压之间的压差、即驱出扩张压与收缩压和扩张压之间的压差、即脉压的比值、即驱出扩张压比。
(61)第(60)条记载的药理作用监视装置,
最好进而具有重博波高比算出部,计算从上述脉波波形得到的重博和重博波峰之间的压差、即重博波高与收缩压和扩张压之间的压差、即脉压的比值、即重博波高比。
【附图的简单说明】
图1是表示第1实施例的脉波诊断装置的构成的方框图。
图2是用于说明第1实施例的脉波诊断装置的动作的流程图。
图3是表示第2实施例的脉波诊断装置的构成的方框图。
图4是表示第2实施例的小波变换部的构成例的方框图。
图5是表示第2实施例的波形整形部的构成的方框图。
图6是表示第2实施例的波形整形部的动作的流程图。
图7是用于说明第2实施例的潮波、重博波检测部的动作的图。
图8是与心电波形、大动脉血压波形和末梢部血压波形对应示出的图。
图9是用于说明脉波波形和波形参数的对应关系的图。
图10是表示血压值差(y5-y4)和畸变率d的相关关系的图。
图11是表示滑脉的频谱分析结果的图。
图12是表示平脉的频谱分析结果的图。
图13是表示弦脉的频谱分析结果的图。
图14是表示各脉象的潮波和重博波的振幅的图表。
图15是表示各脉象的高次谐波成分的图表。
图16是表示弦脉的自相关系数的图。
图17是表示平脉的自相关系数的图。
图18是表示滑脉的自相关系数的图。
图19是表示第3实施例的脉波诊断装置的电路构成的方框图。
图20是表示第3实施例的脉波诊断装置的动作的流程图。
图21是表示第4实施例的脉象数据生成部的方框图。
图22是表示典型的脉波波形的自相关数据的变化率的图。
图23是表示第4实施例的脉波诊断装置的动作的流程图。
图24是表示第5实施例的脉波诊断装置的电路构成的方框图。
图25是表示第6实施例的脉波诊断装置的电路构成的方框图。
图26是表示脉波波形一段时间的脉波分析数据的图。
图27是用于说明第6实施例的体动除去部的动作的时序图。
图28是第6实施例中表示期间Tc内的脉波校正数据MKDa的图。
图29是第6实施例中表示期间Tc内的体动校正数据TKDa的图。
图30是第6实施例中表示体动成分除去后的体动除去脉波数据MKDaj的图。
图31是表示第7实施例的脉波诊断装置的电路构成的方框图。
图32是第7实施例的体动除去部的方框图。
图33是表示第7实施例的体动除去脉波数据的一例的图。
图34是表示第8实施例的脉波诊断装置的电路构成的方框图。
图35是第8实施例的第1小波变换部的方框图。
图36是表示第9实施例的脉波诊断装置的电路构成的方框图。
图37A是表示手表式脉波诊断装置的安装状态的图。
图37B是表示手表式脉波诊断装置的脉波检测部的图。
图37C是表示设在手表式脉波诊断装置的本体部中的连接部的图。
图38是表示脉波检测部的构成例的图。
图39A是表示另一种手表式脉波诊断装置的外观的图。
图39B是表示图39A所示的脉波诊断装置的安装状态的图。
图40是表示项链式脉波诊断装置的外观结构的图。
图41是表示将图40所示的脉波诊断装置的脉波检测部安放在颈动脉处的样子的图。
图42是表示眼镜式脉波诊断装置的外观结构的图。
图43是表示卡片式脉波诊断装置的外观结构的图。
图44A是表示计步器式脉波诊断装置的外观结构的图。
图44B是表示图44A所示的脉波诊断装置的安装状态的图。
图45A是表示典型的平脉脉波波形的图。
图45B是表示典型的滑脉脉波波形的图。
图45C是表示典型的弦脉脉波波形的图。
图46是表示畸变率和脉波波形的形状的相关关系的图。
图47是表示作为第10实施例的理论根据用来记录桡骨动脉的动脉压波形的结构的说明图。
图48是表示典型的动脉压波形的图。
图49是表示使用图47所示的结构得到的实验结果的曲线图,示出平均血压和扩张期血压的关系。
图50是表示使用图47所示的结构得到的实验结果的曲线图,示出平均血压和收缩期血压的关系。
图51是表示使用图47所示的结构得到的实验结果的曲线图,示出脉压和收缩期血压的关系。
图52是表示第10实施例的血压监视装置的方框图。
图53是表示作为第11实施例的理论根据的各脉波形状中的重博波高ΔBPD的分布的图。
图54是表示作为第11实施例的理论根据的重博压差比BPDd/ΔBP和重博波高ΔBPD之间的关系的图。
图55是表示作为第11实施例的理论根据的驱出扩张压ΔBPp和重博波高ΔBPD之间的关系的图。
图56是表示第11实施例的脉波形状监视装置的构成的方框图。
图57是表示第12实施例的药理作用监视装置的构成的方框图。
图58是表示用滤波器组构成小波变换部的变形例的方框图。
图59是表示用滤波器组构成逆小波变换部的变形例的方框图。
图60是表示作为告知装置的变形例中的一种方式的面孔图。
图61是表示变形例中的透过型光电式脉波传感器的例子的图。
图62是表示将光电式脉波传感器用于眼镜式脉波监视装置时的变形例的图。
图63是表示服用降压剂后动脉压波形中的各部分的压力变化的图。
图64是表示服用降压剂前和服用降压剂后的动脉压波形的图。
【实施发明的最佳形态】《1第1实施例》
首先,说明本发明的第1实施例的脉波诊断装置。《1.1第1实施例的理论根据》
不用说,心脏是利用反复收缩扩张来驱出血液的。这里,把因1个周期的收缩扩张而使血液从心脏流出的时间称为驱出时间。当由于运动而使单位时间心脏的收缩次数、即心博数变高时,因肾上腺素等儿茶酚胺被释放,结果,该驱出时间有变短的倾向。这意味着心肌的收缩力增大。
此外,随着该收缩时间的变长,心脏1个周期的收缩扩张所流出的血液量有增大的倾向。
人一旦进行运动等活动,因必须向心肌和骨骼肌供给大量的氧气,故心博数和博出量的乘积、即单位时间心脏送出的血液流量增加。这里,因心博数增加的结果使驱出时间变短,故博出量反而变小。只是,因心博数的增加率比博出量的减少率高,故心博数和博出量的乘积从整体上看增加了。
其次,说明心脏的运动和血液波形的关系。在图9所示的心电波形中,一般,将从R点到T波的终点U的期间称为心室收缩期,这与上述驱出时间相当。此外,将从U点到下一个R点的期间称为心室扩张期。这里,在心室收缩期中,心室的收缩不是一致发生的,而是从外侧向内侧慢慢地进行收缩。因此,心脏后面的大动脉起始部的血压波形便如图9所示,在从大动脉阀打开到关闭的心室收缩期中,呈上凸的形状。
与这样的大动脉起始部的血压波形对应,末梢部(桡骨动脉)中的血压波形、即末梢部的脉波波形也如图9所示那样。之所以变成这样的形状,认为是,首先,因心脏博出血液而产生称之为驱出波的第1波(ejection wave),接着,因心脏附近血管分叉部分的反射而产生称之为潮波的第2波(tidal wave),然后,伴随大动脉阀的关闭而产生重博,从而出现称之为重博波的第3波(dicrotic wave)。
因此,在脉波波形中,从血压值的最低点到重博这一段与心室收缩期相当,从重博到下一个周期的血压值的最低点这一段与心室扩张期相当。
这里,脉波波形中与大动脉阀打开相当的点是血压值的最小极小值的点。脉波波形中与大动脉阀关闭相当的点、即重博,若从时间轴来看,是从该最小极小值的点开始数的第3个极小值的点,若从血压值的大小来看,是从该最小极小值的点开始数的第2个极小值的点。
再有,图9所示的末梢部血压波形、即脉波波形实际上相对大动脉血压波形在时间上延迟了,这里,为了方便说明,忽略该时间延迟,使其相位一致。
下面,讨论末梢部血压波形、即脉波波形。从被测者的末梢部检测出的脉波波形,因为它就是血液经过由作为博动泵的心脏和作为导管的血管系统构成的封闭系统的压力波,所以,第1,除了取决于心脏的泵功能、即心功能状态之外,第2,还受到血管直径、血管的收缩伸展和血液粘性阻抗的影响。因此认为,若检测出脉波波形并进行分析,除了被测者动脉系统的状态外,还能够对心功能状态作出评价,东方医学的专业医生能够根据脉动的形态诊断身体的状态。
这里,讨论一下应该对脉波波形的哪个部分进行分析。
首先,本发明者定出了如图10所示那样的决定脉波波形的形状特征的波形参数。即,定出了以下波形参数。
(1)从脉波波形的血压值最小且为一个心博的上升谷点的P0(最小极小值点)开始到下一拍的上升起始点P6为止的时间t6
(2)脉波波形中依次出现的峰谷点(极大值点和极小值点)P1~P5的血压值(差)y1~y5
(3)从脉波开始时刻的谷点P0(最小极小值点)到出现上述各峰谷点P1~P5的经过时间t1~t5
再有,这时,y1~y5分别表示以谷点P0的血压值为基准的相对血压值。
而且,本发明者实际检测了74名从22岁到46岁的健康成人的脉波,分别求出了这些波形参数,另一方面,与前述的PCT/JP96/01254一样,对该脉波波形进行FFT处理并使用前述式(1)求出了该脉波波形的畸变率d。
而且,本发明者还对已求出的畸变率d和各波形参数及这些参数的差值等的相关关系分别进行了研究,结果表明,畸变率d和重博波(dicrodic wave)的从重博开始的振幅、即血压值差(y5-y4)有很强的相关关系,其相关系数(R2)为0.77。图11示出了该相关关系。
由此,本发明者在提出东方医学的专业医生通过用手指感觉重博波(dicrodic wave)和潮波(tidal wave)的特征来诊断脉象的假说的基础上进行了如下分析。
在分析中,对东方医学的专业医生判定脉象所用的各脉波波形进行了FFT处理,并算出相对基波成分的各高次谐波成分。图11表示滑脉的分析结果,图12表示平脉的分析结果,图13表示弦脉的分析结果。
在图11~图13中,f1是基波,f2是2次谐波,f3是3次谐波,…f10是10次谐波,分别示出了它们的振幅和相位,此外,波形Wf1是基波f1和2次谐波f2相加的波形,波形Wf2是基波f1到3次谐波f3相加的波形,…波形Wf9是基波f1到10次谐波f10相加的波形。
这里,若将图11所示的专业医生判定是滑脉脉象的原始波形与图12所示的判定为平脉脉象的原始波形进行比较,可知,两者虽然相似,但滑脉的重博(dicrotic notch)高度低,平脉的重博波(dicrodic wave)的振幅大。此外,若将它们的合成波形进行比较,可知,滑脉和平脉的重博波(dicrodic wave)利用从基波f1到4次谐波f4相加的波形Wf3基本上可以再现原波形。
另一方面,在图13所示的弦脉中,其潮波(tidal wave)利用从基波f1到7次谐波f7相加的波形Wf6基本上可以再现原波形。
此外,图14是表示图11~图13说明过的脉象的脉波波形中的重博波(dicrodic wave)的振幅和潮波(tidal wave)的振幅的图。重博波(dicrodic wave)的振幅,滑脉是7.3mmHg,平脉是10.6mmHg,比较大,而弦脉则比较小,是2.9mmHg。此外,滑脉和平脉的潮波(tidal wave)的振幅是0,而弦脉潮波(tidal wave)的振幅则是3.8mmHg。
由此,可以说,滑脉和平脉的重博波(dicrodic wave)有些特征,而且重博波(dicrodic wave)的特征可以从基波f1到4次谐波f4反映出来。而弦脉的潮波(tidal wave)有些特征,潮波(tidalwave)的特征可以从高频成分、例如5次谐波f5~7次谐波f7反映出来。
图15是按脉象以百分数表示各高次谐波的振幅对基波振幅之比的图。这里,若考察2次谐波f2~4次谐波f4的振幅之和对基波振幅的比(f2+f3+f4)/f1,则滑脉是1.74,平脉是1.5。因此,根据该值可以判别滑脉和平脉。若考察5次谐波f5~7次谐波f7的振幅之和对基波振幅的比(f5+f6+f7)/f1,则滑脉是0.36,平脉是0.26,弦脉是0.42,因此,根据该值可以判别是不是弦脉。《1.2脉波诊断装置的构成》
本实施例的脉波诊断装置是在上述理论根据的基础上构成的,对从被测者的身体检测出的脉波波形进行频谱分析,抽出潮波(tidalwave)成分和重博波(dicrodic wave),并根据抽出的结果去判定脉象。再有,关于脉波诊断装置的外观结构将在后面的《10上述各实施例的外观结构》一栏中叙述。
图1是表示本实施例脉波诊断装置的功能结构的方框图。在该图中,脉波检测部10例如检测被测者的末梢部(例如桡骨动脉)的脉波波形,将该检测信号作为MH输出给体动除去部30。
另一方面,体动检测部20例如由加速度传感器等构成,检测被测者身体的活动,将该检测信号作为TH输出给波形处理部21。波形处理部21由低通滤波器等构成,对体动检测部20输出的信号TH进行波形整形处理,作为表示体动成分的信号MHt输出。体动除去部30从脉波检测部10输出的信号MH中减去表示体动成份的信号MHt,作为表示脉波成分的信号MHj输出。
本实施例的脉波诊断装置虽然是对从被测者的身体检测出的脉波波形进行处理的装置,但当被测者伴随有任何动作时,在从脉波检测部10检测出的信号MH中,除了表示脉波成分的信号MHj之外,还叠加有表示被测者的体动成份的信号MHt。因此,MH=MHt+MHj,从脉波检测部10输出的信号MH不能正确地反映被测者的脉波波形。
另一方面,因血流受血管和组织的影响,故包含在信号MH中的体动成分MHt并不是表示被测者体动的信号TH本身,而是被钝化了的信号。
因此,体动成份除去部30利用波形处理部21对体动检测部20输出的直接表示被测者体动的信号TH进行波形整形,并作为体动成份的信号MHt使用,从脉波检测部10输出的信号MH中将其减去,据此来除去体动的影响,并作为表示脉波成份的信号MHj输出。再有,波形处理部21中的低通滤波器的形式、级数、常数等由实际测定数据决定。
但是,如果即使没有体动也使体动成分除去部30工作来除去体动,则因体动检测部20的噪声而使体动成分除去部30的输出信号的信噪比降低,还会因体动除去动作而消耗电力。因此,本实施例设有判定部22。判定部22根据体动波形TH生成判定有无体动的控制信号C。具体地说,通过体动波形TH与阈值进行比较来判定。考虑到体动检测部20的噪声电平,预先将该阈值设定成能够判断有无体动。而且,当控制信号C表示无体动时,停止波形处理部21和体动成分除去30的工作。这时,直接从体动成分除去部30输出脉波波形MH。因此,可以改善体动成分除去部30的输出信号的信噪比,还可以降低装置的功耗。
其次,FFT处理部40对表示脉波成分的信号MHj进行FFT处理,并进行脉波成分的频谱分析。作为FFT处理的分析结果,可以得到多根谱线,对每一根谱线算出其频率和能级。FFT处理部40比较这些频谱数据,将能级最高的谱线特定为脉波成分的基波f1。通过取基波f1的频率的整数倍去特定各高次谐波。而且,生成表示基波f1和2次谐波f2~10次谐波f10的能级的脉波分析数据MKD并输出。
其次,潮波特征抽出部50根据脉波分析数据MKD生成表示潮波(tidal wave)特征的潮波特征数据TWD。如上所述,因潮波(tidalwave)特征可以用脉波波形的5次谐波f5~7次谐波f7的和相对基波f1的比值来表示,故潮波特征抽出部50按照下式生成潮波特征数据TWD。
        TWD=(f5+f6+f7)/f1
其次,重博波特征抽出部60根据脉波分析数据MKD生成表示重博波(dicrotic wave)特征的重博波特征数据DWD。如上所述,因重博波(dicrotic wave)特征可以用脉波波形的2次谐波f2~4次谐波f4的和相对基波f1的比值来表示,故重博波特征抽出部60按照下式生成潮波特征数据DWD。
        TWD=(f2+f3+f4)/f1
其次,脉波判定部70根据潮波特征数据TWD和重博波特征数据DWD进行脉象的判定,生成表示被测者脉象的种类的脉象数据ZD。具体地说,首先,将潮波特征数据TWD与第1阈值进行比较,当潮波特征数据TWD高于第1阈值时,生成表示是弦脉的脉象数据ZD1。预先将第1阈值设定成能够判别是不是弦脉,在该例中将其设定为0.41。
另一方面,当潮波特征数据TWD低于第1阈值时,将重博波特征数据DWD与第2阈值进行比较,当重博波特征数据DWD低于第2阈值时,生成表示是平脉的脉象数据ZD2。另一方面,当重博波特征数据DWD高于第2阈值时,生成表示是滑脉的脉象数据ZD3。这里,预先将第2阈值设定成能够判定是平脉还是滑脉,在该例中将其设定为1.62。
告知部80是将脉象数据ZD显示出来或利用声音等向外部输出的装置,例如,显示与各脉象对应的文字,如‘滑脉’、‘平脉’、‘弦脉’,或显示图标等符号。因此,被测者和医生等第3者能够知道脉象。〖1.3脉波诊断装置的动作〗
下面,参照图2说明第1实施例的脉波诊断装置的动作。
首先,在脉波检测部10输出的信号MH上叠加了伴随被测者的体动的体动成分,通过体动成分除去部30除去该体动成分,变成只表示脉波成分的信号MHj,并供给FFT处理部40(步骤S1,S2)。
其次,在FFT处理部40中,对信号MHj进行FFT处理,生成脉波成分的基波f1和各高次谐波f2~f10,作为脉波分析数据MKD,并将其供给潮波特征抽出部50和重博波特征抽出部60(S3)。
然后,潮波特征抽出部50计算表示潮波的特征的脉波波形的5次谐波f5~7次谐波f7的和相对基波f1的比值,并生成潮波特征数据TWD。重博波特征抽出部60计算表示特征波的特征的脉波波形的2次谐波f2~4次谐波f4的和相对基波f1的比值,并生成重博波特征数据DWD(S4)。
这样,当抽出脉波波形的潮波(tidal wave)和重博波(dicroticwave)的特征时,脉象判定部70首先将潮波特征数据TWD与第1阈值比较(S5)。若潮波特征数据高于第1阈值(0.41),则进到S6,生成表示弦脉的脉象数据ZD1。
另一方面,若潮波特征数据低于第1阈值(0.41),S5的判定结果为否,则进到S7,根据重博波特征数据DWD判定脉象。这时,脉象判定部70判定重博波特征数据DWD是否低于第2阈值(1.62),若是低于,则进到S8,生成表示平脉的脉象数据ZD2。另一方面,若重博波特征数据DWD高于第2阈值(1.62),则S7的判定结果是否,进到S9,生成表示滑脉的脉象数据ZD3。
这样,在第1实施例中,诊脉熟练的专业医生从根据脉波波形的潮波和重博波判定脉象的观点出发,同时根据潮波和重博波的特征可以从规定的高次谐波上反映出来这一点,对脉波波形进行频谱分析,并抽出潮波和重博波的特征,从而能够客观而准确地判断脉象。《2第2实施例》
下面,说明本发明第2实施例的脉波诊断装置。
在上述第1实施例中,根据滑脉和平脉的特征部分在重博波(dicrotic wave),而弦脉的特征部分在潮波(tidal wave)的事实,对脉波波形进行FFT处理,从该分析结果中分别抽出与重博波和潮波对应的频率成分,并根据它去判定脉象。
在FFT处理的频谱分析中,当有相当长时间的数据没有时,就得不到分析结果。因此,一般的FFT处理,往往是对几个周期的脉波波形进行处理。
但是,重博波和潮波出现在脉波波形的一段时间内。因此,若只对脉波波形中的重博波和潮波的这一段时间进行频谱分析,就可以提高脉象的判定精度。
因此,第2实施例利用对脉波波形同时进行频谱分析和时间域分析的小波变换来特定脉象。《2.1脉波诊断装置的总体构成》
图3表示第2实施例的脉波诊断装置的构成。图3所示的脉波诊断装置除了使用小波变换部41代替FFT处理部40、设置了特定潮波和重博波的时间位置的潮波、重博波检测部42、潮波特征抽出部50和重博波特征抽出部60以及脉象判定部70的内部结构不同这些点之外,与图1所示的第1实施例的脉波诊断装置相同。下面,说明其不同点。《2.2小波变换部》
首先,与附图一起详细说明小波变换部41的构成。
一般,在从时间和频率两方面同时去捕捉信号的时间频率分析中,小波是将信号分割成各个部分的单位。小波变换表示以这样的单位分割的信号的各部分的大小。为了定义小波变换,作为基本函数将在时间上和频率上都定域化了的函数ψ(x)作为母小波导入。这里,函数f(x)的母小波ψ(x)的小波变换定义如下。 ( W ψ f ) ( b , a ) = ∫ - ∞ ∞ 1 a ψ ( x - b a ) f ( x ) dx - - - - ( 2 )
在式(2)中,b是移动(平行移动)母小波ψ(x)时使用的参数,a是伸缩时的参数。因此,在式(2)中,小波ψ((x-b)/a)是将母小波ψ(x)平行移动b并伸缩a的结果。这时,母小波ψ(x)的幅度与尺度参数a对应而伸长,所以,1/a是与频率对应的量。小波变换部41构成为能够执行式(2)的运算,其详细结构如图4所示。
在图4中,从体动成分除去部30输出的信号MHj供给波形整形部400和A/D变换部410。波形整形部400生成与脉波波形MHj同步的控制信号CS和时钟信号CK。
这里,图5示出波形整形部400的方框图。在图5中,振铃滤波器401是中心频率为2.2Hz、通频带为0.8Hz~3.5Hz的高Q值滤波器。脉波波形的基波成分通常在0.8Hz~3.5Hz的范围内,所以,当脉波波形MHj通过振铃滤波器401后,其基波成分被抽出。例如,如图6所示,脉波波形MHj输入振铃滤波器401后,其输出变成图6中的401out所示的波形。
其次,过零检测电路402由比较器等构成,将振铃滤波器401输出的信号与接地电平进行比较,生成矩形波。该矩形波是与脉搏同步的信号。例如,若振铃滤波器401的输出信号是图6中的401out所示的信号,则过零检测电路102的输出信号便成为该图中的402out所示的信号。
其次,比较部403、环路滤波器404、压控震荡电路405和分频电路406构成锁相环路。当比较部403的一个输入端供给过零检测电路402的输出信号,而另一端供给分频电路406的输出信号时,比较部403便输出与两者的相位差对应的误差信号。当误差信号经环路滤波器404供给压控震荡电路405时,压控震荡电路405便输出时钟信号CK。而且,时钟信号CK由分频电路406进行1/N分频,再反馈到比较部的另一个输入端。例如,若分频比是1/8,则时钟信号CK的频率便如图6的405out所示那样,是过零检测电路402的输出信号的频率的8倍。然后,时钟信号CK由分频电路407进行1/2分频,作为图6所示的控制信号CS输出。
此外,脉波波形MHj由图4所示的A/D变换器410变换成数字信号后,存储在第1存储器420和第2存储器430中。这里,第1存储器420的写使能端直接供给控制信号CS,第2存储器430的写使能端供给经反相器440反相后的控制信号CS。因此,第1、第2存储器420、430以时钟周期为单位交替存储脉波波形。
此外,450是多路切换器,有选择地将交替从第1、第2存储器420、430读出的脉波数据MD输出到基本函数展开部W。这样,在第1存储器420的写入期间,从第2存储器430读出脉波数据MD,在第1存储器420的读出期间,向第2存储器430写入脉波数据MD。
其次,基本函数展开部W是进行上述式(2)的运算的结构,当供给上述时钟信号CK时,便以时钟周期进行运算。基本函数展开部W由存储母小波ψ(x)的基本函数存储部W1、变换尺度参数a的尺度变换部W2、缓冲存储器W3、进行移动的平行移动部W4和乘法运算部W5构成。再有,作为存储在基本函数存储部W1中的母小波ψ(x),除了Gabal小波之外,还可以使用墨西哥草帽小波、Haar小波、Meyer小波、Shannon小波等,在该例中,使用称之为symlets5的母小波,能够很好地抽出脉波波形的特征。
首先,当从基本函数存储部W1读出母小波ψ(x)时,尺度变换部W2进行尺度参数a的变换。这里,尺度参数a是与周期对应的参数,故当a增大时,母小波ψ(x)便时间轴上伸长。这时,基本函数存储部W1存储的母小波ψ(x)的数据量是一定的,所以,a增大时,单位时间的数据量便减少。尺度变换部W2为了弥补这一点,进行内插处理,同时,当a变小时,进行间拔处理,生成函数ψ(x/a)。将该数据暂时存储在缓冲存储器W3中。
其次,平行移动部W4通过以与移动参数B对应的时序从缓冲存储器W3中读出函数ψ(x/a),生成将函数ψ(x/a)平行移动后的函数ψ(x-b/a)。
其次,乘法运算部W4将变量1/a1/2、函数ψ(x-b/a)和脉波数据MD相乘后,以心博单位进行小波变换,并生成脉波分析数据MKD。在该例中,脉波分析数据MKD在被分割成例如0Hz~0.5Hz、0.5Hz~1.0Hz、1.0Hz~1.5Hz、1.5Hz~2.0Hz、2.0Hz~2.5Hz、2.5Hz~3.0Hz、3.0Hz~3.5Hz、3.5Hz~4.0Hz频率区之后再输出。《2.3潮波、重博波检测部》
潮波、重博波检测部42通过将存在脉波分析数据MKD的某频率区与阈值进行比较来分别生成特定潮波(tidal wave)和重博波(dicrotic wave)的时间位置的控制信号ct、cd。
例如,在图7所示的脉波波形MHj中,tw是潮波,dw是重博波。当对该脉波波形MHj进行小波变换时,可以得到图7所示的小波分析结果。该分析结果将脉波分析数据MKD的示出值以浓淡的形式表示,越浓的部分脉波分析数据MKD的值越大,越淡的部分脉波分析数据MKD的值越小。
由该图可知,在与潮波(tidal wave)对应的时间区间内,存在被脉波分析数据MKD的值小的部分(白色)所包围的区域A和区域B。在表示频率的Y轴方向上脉波分析数据MKD的值小的部分表示在该时间内脉波波形MHj的能量低,该部分脉波波形MHj变得平坦。例如,区域A和区域B的交界处与潮波tw的峰值点一致。
因此,若着眼于某一个频率区,根据其变化可以检测潮波tw的时间位置。在该例中,通过将频率区X中的脉波分析数据MKD与阈值进行比较,生成图7所示那样的特定潮波tw的时间位置t1~t2的控制信号ct。另一方面,通过将频率区X中的脉波分析数据MKD与阈值进行比较,生成图7所示那样的特定重博波dw的时间位置t3~t4的控制信号cd。《2.4潮波特征抽出部、重博波特征抽出部》
首先,潮波特征抽出部50根据控制信号ct特定潮波的时间位置,从该期间的脉波分析数据MKD中,将特定的频率区域内的数据相加。
例如,就象图58所示的在后面《14其它实施例》(6)中叙述的那样,由滤波器组构成小波变换部41,利用与脉波波形MHj同步的时钟信号CK,使高通滤波器1A和低通滤波器1B的特性可变。这时,若与脉波波形MHj的基波f1对应的脉波分析数据是M*1,因潮波的特征部分由5次谐波f5~7次谐波f7表现,故将M*5、M*6、M*7相加后算出潮波特征抽出数据TWD。
其次,重博波特征抽出部60根据控制信号cd特定重博波的时间位置,从该期间的脉波分析数据MKD中,将特定的频率区域内的数据相加。
例如,就象在潮波特征抽出部50中说明过的那样,若与脉波波形MHj的基波f1对应的脉波分析数据是M*1,重博波的特征部分由2次谐波f2~4次谐波f4表现。因此,将M*2、M*3、M*4相加后算出重博波特征抽出数据DWD。
这样,第2实施例的潮波特征抽出部50和重博波特征抽出部60,因为巧妙地利用小波变换的分割频率域和时间域来进行波形分析的性质,不仅从频率域、而且从时间域抽出潮波(tidal wave)和重博波(dicrotic wave)的特征,所以,能够以高的精度抽出特征。《2.5脉象判定部、告知部》
脉象判定部70通过将潮波特征抽出数据TWD和重博波特征抽出数据DWD与阈值进行比较来判定滑脉、平脉和弦脉的脉象,并生成脉象数据ZD。此外,告知部80与第1实施例一样,通过显示脉象数据ZD和声音等告诉被测者或医生等的第3者。
如上所述,若按照第2实施例的脉波诊断装置,通过对从被测者的身上检测出的脉波波形进行小波变换,能够将脉波分析的时间限定在潮波(tidal wave)和重博波(dicrotic wave)存在的期间进行频谱分析。结果,能够高精度地抽出潮波(tidal wave)和重博波(dicrotic wave)的特征,进行正确的脉象判定。《3第3实施例》《3.1第3实施例的原理》
因为脉波波形表示因心脏的收缩而向大动脉送出的血流在动脉中传送时的脉动,所以具有与心脏的博动同步的一定的周期。作为分析周期性的波形并抽出其特征的分析方法,虽然有以FFT(快速付里叶变换)为代表的频谱分析等方法,但本发明者还注意到可以通过简单的运算进行处理的自相关函数。
这里,用x(t)表示不规则变动,若在x(t)中有周期T的周期变动,x(t)可由下式给出。
x(t)=x(t±nT)         这里,n=0,1,2,…
即,只要错开周期的整数倍,就与原来的波形重合了。若不规则变动x(t)的周期性强,只要在时间轴上错动周期的整数倍,就成为与原来波形相似的波形。因此,研究只错动了某一时间τ的波形与原来的波形到底有多大差别,为了判别变动中周期成分,只需要求x(t)和x(t+τ)的相关即可。
当设x(t)是与时间有关的不规则变量时,自相关函数就可以定义为时间相隔τ的两个变量的乘积的平均值,可有下式给出。
C(τ)=E[x(t)x(t+τ)]
这里,虽然E是总体平均,但在正常确立过程中可用时间平均来替换。因此,自相关函数C(τ)可以用下式来表示。 C ( τ ) = lim T → ∞ 1 T ∫ - T / 2 T / 2 x ( t ) x ( t + τ ) dt
此外,把用τ=0对自相关函数C(τ)归一化了的自相关函数称为自相关系数R(τ),自相关系数可由下式给出。
R(τ)=C(τ)/C(0)
     =E[x(t)x(t+τ)]/E[x2(t)]
这里,图16是图45C所示的弦脉的自相关系数,图17是图45A所示的平脉的自相关系数,图18是图45B所示的滑脉的自相关系数。若将这些图进行比较,可以看出随着弦脉→平脉→滑脉的变化自相关系数R(τ)存在起伏。具体地说,弦脉是平滑的波形,平脉是存在三个峰,峰的大小比较小,滑脉由两个峰构成。即,自相关系数R(τ)反映各脉波波形的特征。而且,因自相关系数R(τ)的波形形状是根据一定长时间的脉波波形得到的,故能够高精度地抽出脉波波形的特征。
此外,因自相关系数R(τ)的瞬时值与脉波波形的瞬时值不同,由时间平均给出,所以,即使从身体上检测出的脉波波形的信噪比较差,通过将噪声平均可以改善信噪比。进而,因自相关系数R(τ)是将自相关函数归一化后得出的,故对比较分析振幅值不同的脉波波形很有效。
因上述理由,本发明者以巧妙地利用自相关系数R(τ)的性质去特定脉象为根据,发明了由此判定身体状态的脉波诊断装置。《3.2脉波诊断装置的电路构成》
下面,参照图19说明脉波诊断装置的电路构成。图19是表示脉波诊断装置的电路构成的方框图。
脉波诊断装置1大致由检测脉波波形MH的脉波检测部10、根据脉波波形MH计算自相关数据RD的自相关运算部210、根据自相关数据RD生成脉象数据ZD的脉象数据生成部220和显示部230构成。
首先,自相关函数运算部210由波形处理部211、存储器212、乘法运算部213、平均值算出部214和归一化运算部215构成。存储器212至少存储1心博期间的脉波波形MH。但是,若从上述数1开始计算自相关函数C(τ),必须从-T/2到T/2对x(t)x(t+τ)进行积分(T为无限大),因该例中x(t)是与心博周期同步的脉波波形MH,所以,只要根据至少1个周期的脉波波形MH进行运算即可。因此,在该例中,利用波形处理部211对脉波波形MH进行整形,使其变成与脉波周期同步的方波,根据该方波生成控制存储器212的写入动作的写入控制信号WE。例如,利用该写入控制信号WE将4个周期的脉波波形MH写入存储器212中。其次,乘法运算部213读出与x(t)相当的脉波波形MH1和与x(t+τ)相当的脉波波形MH2,将两者相乘并输出相乘的结果。而且,平均值运算部214计算与x(t)x(t+τ)的总体平均相当的MH1·MH2的总体平均,并输出该计算结果。
其次,归一化运算部215用自相关函数C(0)将平均值算出部的计算结果归一化并生成自相关数据RD。具体地说,由x2(t)、即MH12的总体平均进行归一化。因此,自相关数据RD表示脉波波形MH的自相关系数R(τ)。
这里,对图45A~图45所示的各脉象的脉波波形MH研讨其自相关系数R(τ)。首先,图16所示的弦脉其自相关系数R(τ)的最小值是0.2,图17所示的平脉其自相关系数R(τ)的最小值是0.32,图1所示的滑脉其自相关系数R(τ)的最小值是0.3。因此,当自相关系数的最小值低于0.25时,可以判定它是弦脉。
其次,自相关系数R(τ)在0.5以上的期间相对1心博期间所占的比例,平脉是60%(参照图17),滑脉是34%(参照图18)。因此,可以判定,若超过47%是平脉,低于47%是滑脉。
脉象数据生成部220根据上述判定基准判定脉象,由判断是弦脉的最小值检测部221和第1比较部222、判断是平脉的时间测量部223、运算部224和第2比较部225以及生成脉象数据ZD的数据生成部226构成。
首先,最小值检测部221在与1心博周期对应的期间内检测自相关数据RD的最小值。第1比较部222判定最小值是否低于0.25,并将该结果输出给数据生成部226。
其次,时间测量部223将自相关数据RD与预定的阈值(0.5)比较,测量自相关数据RD高于阈值(0.5)的时间间隔。此外,运算部224计算测量出的时间间隔对1心博期间所占的比例。再有,1心博期间由上述波形处理部211供给。此外,第2比较部225判定计算结果是否超过了预定的阈值(47%)。
其次,数据生成部226根据第1比较部222和第2比较部225的判定结果生成脉象数据ZD。首先,当第1比较部222判断超过0.25时,生成表示弦脉的数据ZD1。此外,当第2比较部255判断自相关数据RD超过0.5的时间所占的比例超过47%时,生成表示平脉的脉象数据ZD2。当判断该时间所占的比例低于47%时,生成表示滑脉的脉象数据ZD3。
其次,显示部230由ROM、控制电路和液晶显示装置构成。当向显示部230供给脉象数据ZD时,控制电路对此进行检测,读出存储在ROM中的字符,并将其在液晶显示器上显示出来。作为字符除了使用‘平脉’、‘弦脉’、‘滑脉’等文字外,也可以使用特定的记号或图标。由此,能够向使用者或医生告知健康状况。
这样,就构成了脉波诊断装置的功能,但在实际装置中,自相关运算部210和脉象数据生成部220由CPU和存储器等构成。这时,CPU根据存储在一部分ROM中的控制程序执行各种运算处理和比较处理并生成脉象数据ZD。《3.3脉波诊断装置的动作》
下面,参照附图说明第3实施例的脉波诊断装置的动作。图20是表示脉波诊断装置的动作的流程图。
首先,利用脉波检测部10检测脉波波形MH(S1)。
当将该脉波波形MH供给自相关运算部210时,脉波波形MH被写入存储器212,延迟规定时间后再读出。根据读出的脉波波形MH进行自相关函数的运算,将该运算结果归一化并生成表示自相关系数R(τ)的自相关数据RD(S2)。
然后,如下所述,脉象数据生成部220根据自相关数据RD生成脉象数据ZD。首先,当最小值检测部221检测出与1心博期间对应的自相关数据RD的最小值时(S3),第1比较部222判定最小值是否低于0.25(S4)。当最小值低于0.25时,判定结果是是,进入S5,判定为弦脉,数据生成部226生成脉象数据ZD1。
若最小值高于0.25,判定结果是否,进行是平脉还是滑脉的判定。这时,时间测量部223在与1心博期间对应的期间内测量自相关数据RD超过0.5的时间间隔(S6),接着,运算部24计算测量的时间间隔在1心博期间内所占的比例。
然后,第2比较部225判定运算部224的计算结果是否超过47%(S8),当超过47%时,进入S9,判定为平脉,数据生成部226生成脉象数据ZD2。当没有超过47%时,进入S10,判定为滑脉,数据生成部226生成脉象数据ZD3。
这样,若按照第3实施例的脉波诊断装置1,因为根据脉波波形MH的自相关数据RD判定脉象,所以,即使没有诊脉的知识,被测者和医生也能够从脉波的形态中了解身体的状态。《3.4第3实施例的变形例》
(1)在上述第3实施例中,最小值检测部221是算出1心博期间的自相关数据RD,但也可以在多个各心博期间内检测出的自相关数据RD的各最小值对其进行平均,检测出平均最小值。这时,因最小值是平均后的值,故可以抑制因噪声引起的最小值的变动,能够提高是不是弦脉的判定精度。
(2)在上述第3实施例中,运算部224是算出测量的时间间隔在1心博期间所占的比例,但也可以在多个各心博期间检测出的各个比例进行平均并将其输出。这时,因比例是平均后的值故能够抑制由噪声引起的变动,所以,能够提高是平脉还是滑脉的判定精度。
(3)若将典型的平脉和滑脉的自相关系数R(τ)进行比较(参照图17和图18),发现自相关系数R(τ)的值在0.4到0.8的范围内自相关系数R(τ)在时间轴上的宽度都不相同。因此,时间测量部223也可以测量高于在0.4到0.8的范围内指定的任何阈值的时间间隔,第2比较部使用与该阈值对应的值来判别是平脉还是滑脉。《4第4实施例》
下面,说明第4实施例的脉波诊断装置。《4.1脉波诊断装置的电路构成》
第4实施例的脉波诊断装置的电路构成除了脉象数据生成部220的构成外与第3实施例相同。
下面,说明第4实施例的脉象数据生成部220的电路构成。图21是第4实施例的脉象数据生成部220的方框图。
本例的脉象数据生成部220由判定是不是弦脉的最小值检测部221和第1比较部222、判定是平脉还是滑脉的变化率运算部227、最大值检测部228和第2比较部225、以及根据该判定结果生成脉象数据ZD的数据生成部226构成。
首先,第1比较部222与第3实施例一样,通过判定由最小值检测部221检测出的自相关数据RD的最小值是否超过0.25来判定脉象是不是弦脉。
其次,变化率运算部227算出自相关数据RD的变化率RDd。例如,若自相关数据RD是如图16~图18所示那样,则弦脉、平脉和滑脉的变化率RDd如图22所示。从该图可知,滑脉变化率RDd的最大值大约是0.1,平脉的最大值是0.072。这与平脉的脉波波形MH和滑脉的脉波波形MH相比起伏较大是对应的,可以根据变化率RDd的最大值判定是平脉还是滑脉。因此,当变化率RDd的最大值超过0.085时,可以判定它是滑脉,当没有超过时,可以判定它是平脉。
其次,最大值检测部228检测至少比1心博期间长的规定时间内的变化率RDd的最大值。此外,第2比较部225将变化率的最大值与预定的阈值(0.85)进行比较,判定变化率RDd的最大值是否超过了阈值。
其次,数据生成部226根据第1比较部222和第2比较部225的判定结果生成脉象数据ZD。首先,当第1比较部222判定最小值超过0.25时,生成表示弦脉的脉象数据ZD1。此外,当第2比较部225判定变化率RDd的最大值低于0.085时,生成表示平脉的脉象数据ZD2,当判定高于0.085时,生成表示滑脉的脉象数据ZD3。《4.2脉波诊断装置的动作》
下面,参照附图说明第4实施例的脉波诊断装置的动作。图23是表示脉波诊断装置的动作的流程图。再有,因从S1到S5的动作与使用图20说明过的第3实施例的脉波诊断装置的动作相同,故在此省略其说明。
首先,在S11中,当变化率运算部227检测出自相关数据RD的变化率RDd时,最大值检测部228检测出与1心博期间对应的期间内变化率RDd的最大值(S12)。例如,若计算出的变化率RDd是图22所示的滑脉,则检测出来的最大值大约是0.1。
其次,第2比较部225判定检测出的最大值是否超过0.85,当低于0.85时,进入S14,判定为平脉,数据生成部226生成脉象数据ZD2。当高于0.85时,进入S15,判定为滑脉,数据生成部226生成脉象数据ZD3。
这样,在第4实施例中,根据平脉和滑脉其自相关数据RD的变化率RDd的最大值不同,将变化率RDd的最大值与阈值进行比较,从而判定是平脉还是滑脉,因此,即使是没有诊脉知识的人也能够知道正确的脉象。《5第5实施例》
在上述第3和第4实施例的脉波诊断装置1中,根据从脉波检测部10输出的脉波波形MH生成自相关数据RD并判定脉象。但是,在被测者步行或进行日常活动时,因身体的运动而使血流产生变动。因此,在从脉波检测部10输出的脉波波形MH上叠加了体动成分。因体动成分的大小依赖被测者运动的程度,当被测者的运动量大时重叠在脉波波形MH上的体动成分也增加,故难以正确地判定脉象。因此,在第5实施例中,从脉波波形MH中除去体动成分并据此进行脉象的判定。《5.1脉波诊断装置的电路构成》
图24是表示第5实施例的脉波诊断装置的电路构成的方框图。再有,因本例中的自相关运算部210、脉象数据生成部220和显示部230的构成与第3或第4实施例说明过的构成一样,故在此省略其说明。此外,因体动成分除去部30及其前级的各构成部、即脉波检测部10、体动检测部20、判定部22和波形处理部21与第1实施例中说明过的一样,故省略其说明。
本实施例的脉波诊断装置因为通过以上的构成来生成体动除去脉波波形MHj,所以,自相关运算部10能够根据体动除去脉波波形MHj生成自相关数据RD。因此,若按照该脉波诊断装置1,特定脉象时可以不受体动的影响。《6第6实施例》
下面,说明第6实施例的脉波诊断装置。
图25是第6实施例的脉波诊断装置1的方框图。第6实施例与第5实施例一样使用体动检测部20和波形处理部21检测体动成分TH,但与第5实施例说明过的使用小波变换进行体动除去这一点是不同的。《6.1第1、第2小波变换部和第1、第2频率修正部》
在图25中,第1小波变换部243对从脉波检测部10输出的脉波波形MH进行小波变换,并生成脉波分析数据MKD。第2小波变换部245对从脉波检测部10输出的脉波波形MH进行小波变换,并生成脉波分析数据TKD。再有,第1和第2小波变换部243、245与第2实施例说明过的一样。
图26是表示脉波波形MH部分时间的脉波分析数据MKD的图。在该图中,期间T在峰值P4的附近,脉波分析数据MKD是以将期间T分成8份的时间间隔得到的。但是,在小波变换中,频率分辨率和时间分辨率具有折衷选择的关系,所以,若牺牲频率分辨率,就能够以比较短的时间间隔得到脉波分析数据。
其次,第1频率修正部244对脉波分析数据MKD进行频率修正。在上述式(2)中,存在1/a1/2的项,当比较不同频率区间的数据时,有必要修正该项的影响。第1频率修正部244是为此而设的,将小波数据WD乘以系数a1/2,生成脉波修正数据MKDa。由此,可以根据对应的各频率进行修正,使单位频率的功率密度一定。此外,第2频率修正部246与第1频率修正部244一样,进行频率修正,从体动分析数据TKD生成体动修正数据TKDa。《6.2体动成分除去部》
体动成分除去部240从脉波修正数据MKDa中减去体动修正数据TKDa生成体动除去脉波数据MKDaj。下面,就这一点具体地进行说明。再有,在下面的说明中,设想被测者用手举起杯子再把它放回原来的位置。这时,由脉波检测部10检测出图27所示的脉波波形MH,此外,由波形处理部21检测出体动波形TH。
这里,体动波形TH从时间T1开始增加,在时间T2到达峰值,然后逐渐减小,在时间T3通过0电平,在时间T4到达负的峰值,在时间T5又回到0电平。但是,因体动波形TH是由体动检测部20使用加速度传感器检测的,故时间T3是与被测者将杯子举到最高的时刻对应的,时间T1与开始上举的时刻对应,时间T4与放下杯子的时刻对应。因此,从时间T1到时间T4的期间是体动存在的期间。再有,图27的脉波波形MHj假定是没有体动时脉波波形。此外,在本例中,脉波波形的基波频率是1.3Hz。
这里,将期间Tc(参照图27)的脉波修正数据MKDa示于图28,将期间Tc的体动修正数据TKDa示于图29。从该图看出,在体动波形TH上,在0.0Hz~1.0Hz的频带内存在电平较大的频率成分。当向体动成分除去部240供给脉波修正数据MKDa和体动修正数据TKDa时,体动成分除去部240从脉波修正数据MKDa中减去体动修正数据TKDa,生成图30所示的已除去体动成分的体动除去脉波数据MKDaj。由此,即使在体动存在的情况下也能消除其影响。《6.3判定部》
判定部22将体动波形TH与预定的阈值比较,生成表示有无体动的控制信号C,并将其供给波形处理部21、第2小波变换部245和第2频率修正部246。因此,当没有体动时,波形处理部21、第2小波变换部245和第2频率修正部246停止工作,可以实现降低运算处理时间、降低电力消耗和提高信噪比。《6.4逆小波变换部》
逆小波变换部247对体动除去脉波数据MKDaj进行逆小波变换,生成体动除去脉波波形MHj。这时,逆小波变换部247进行式(3)所示的逆小波变换后进行波形的再合成。
Figure A9880351600491
其次,当体动成分除去部240从检测出的脉波波形MH中减去体动成分TH生成体动除去脉波波形MHj时,与第5实施例一样,自相关运算部210根据体动除去脉波波形MHj生成自相关数据RD。
此后,当脉象数据生成部220根据自相关数据RD生成脉象数据ZD时,显示部230显示表示脉象数据ZD的文字(弦脉、平脉、滑脉)或与各脉象对应的符号。因此,被测者和医生等第3者能够认识脉象。
如上所述,若按照第6实施例,根据已除去体动的体动除去脉波波形MHj生成自相关数据RD,所以,被测者在日常生活中也能连续地检测脉象。《7第7实施例》
在第5、第6实施例中,利用体动检测部20检测体动波形TH,将脉波波形MH与体动波形TH比较,消除包含在脉波波形MH的频率成分中的体动成分,算出自相关函数RD,据此来特定脉象。但是,因必需要体动检测部20和波形处理部21等,故结构复杂。第7实施例是鉴于这一点提出的,它提供一种脉波诊断装置,即使在有体动的情况下,也能够以简单的结构正确地诊断脉象。
图31是第7实施例的脉波诊断装置的方框图,除了省略了体动检测部20、波形处理部21、第2小波变换部245和第2频率修正部246以及除了体动成分除去部240的内部结构之外,与图25所示的第7实施例的脉波诊断装置1相同。下面,就不同点进行说明。
体动成分除去部240从脉波修正数据MKDa中分离除去体动成分并生成体动除去脉波数据MKDaj。在此,体动除去部240利用了下述体动的性质。
体动是因手腕的上下运动和行走时手腕的摆动而产生的,在日常生活中,人体几乎没有瞬间的运动。因此,在日常生活中,体动波形TH的频率不太高,通常在0Hz~1Hz的范围内。这时,脉波波形MH的基波频率多数在1Hz~2Hz的范围内。因此,在日常生活中,体动波形TH的频率在比脉波波形MH的基波频率低的频率区内。
另一方面,在轻微的运动中,因主要受手腕摆动的影响,体动波形TH的频率有所提高,但因心率随运动量的增加而增加,故脉波波形MH的基波频率也同时提高了。因此,即使在运动中,通常,体动波形TH的频率也在比脉波波形MH的基波频率还低的频率区内。
体动成分除去部240是根据这一点去分离体动成分的,将比脉波波形MH的基波成分还低的频率区滤除。这时,若比脉波波形MH的基波成分高的频率区存在体动成分,则心功能的检测精度降低。但是,因上述体动成分在比脉波波形MH的基波成分还低的频率区的概率高,故能够以高精度诊断脉象的状态。
图32是体动成分除去部240的方框图。波形整形部301对脉波波形MH进行波形整形,生成与脉波波形MH同步的的复位脉冲。计数器302对未图示的时钟脉冲进行计数,利用上述复位脉冲使计数值复位。此外,平均值算出电路303算出计数器302的计数值的平均值。这时,利用平均值算出电路303算出的平均值与脉波波形MH的平均周期对应。因此,若参照平均值,就可以知道脉波波形MH的基波频率。
其次,置换电路304根据上述平均值特定包含脉波波形MH的基波频率的频率区。例如,当上述平均值是0.71秒时,基波频率是1.4Hz,所以,特定的频率区是1Hz~1.5Hz。然后,置换电路304在小于特定频率区的频率区内,将脉波修正数据MKDa置换成‘0’并生成体动除去脉波数据MKDaj。由此,能够滤除了比脉波波形MH的基波频率低的频率区的成分。这时,脉波成分与体动成分一起也被置换成‘0’,但因脉波波形MH的特征部分在比基波频率高的频率区内,故即使被置换成‘0’,对脉波波形几乎没有影响。
例如,若图27所示的脉波波形MH(基波频率是1.3Hz)是利用脉波检测部10检测出来的,则Tc期间的脉波修正数据MKDa如图28所示。
这时,因利用置换电路304特定的频率区是1Hz~1.5Hz,故作为置换对象的频率区成为0.5Hz~1Hz和0Hz~0.5Hz。因此,脉波修正数据MKDa中的0~1Hz被置换成‘0’,并生成图33所示的体动除去脉波数据MKDaj。
当利用逆小波变换247将这样得到脉波修正数据MKDa变换成体动除去脉波波形MHj时,则与第5实施例一样,自相关运算部210可以根据体动除去脉波波形MHj生成自相关数据RD。然后,当脉象生成部220根据自相关数据RD生成脉象数据ZD时,显示部230显示表示脉象数据ZD的文字(弦脉、平脉、滑脉)或与各脉象对应的符号。因此,被测者和医生等第3者能够认识脉象。
如上面说明过的那样,若按照第7实施例,根据除去了体动的体动除去脉波波形MHj生成自相关数据RD,所以,在被测者的日常生活中也能连续地检测脉象。
此外,若按照第7实施例,巧妙地利用体动成分在比脉波波形MH的基波频率成分还低的频率区内的概率高的体动的性质来除去体动成分。因此,能够省略第3、第4实施例必需要的体动检测部20和波形处理部21,而且即使存在体动时也能够正确地诊断脉象。《8第8实施例》
在上述第3和第4实施例中,根据脉波波形MH生成自相关数据RD,在上述第5~第7实施例中,根据体动除去脉波波形MHj生成自相关数据RD。而在上述小波变换中,能够按每一个分割后的时间频率得到分析结果。这里,脉波波形MH和体动除去脉波波形MHj的特征部分是正或负的峰值。此外,当出现这些峰值时,在比较高的频率域内,小波变换分析数据的值变大。因此,若着眼于某个频率区并对所注意的频率区的分析数据生成自相关数据RD,则可以据此特定脉象。第8实施例是根据这一点提出的。
图34示出第8实施例的脉波诊断装置的方框图。图中,当利用脉波检测部10检测出脉波MH时,第1小波变换部243对脉波波形MH进行小波变换,生成与特征的频率区对应的脉波分析数据MKDf。例如,象图26所示的小波变换那样,将频率区分割成8份,若所注意的频率区是3.0Hz~2.5Hz,将M16~M86作为脉波分析数据MKDf输出。
这时,第1小波变换努243的构成如图35所示。将图35与图4比较,可知省略了尺度变换部W2。这是因为,上述尺度变换部W2进行与周期对应的尺度参数a的变换,只要将与所注意的频率区对应的母小波存储在基本函数存储部W1中即可。
这样生成的脉波分析数据MKDf表示脉波波形MH的特征部分的时间变化,所以,通过算出脉波分析数据MKDf就能够有效地进行脉波波形MH的分析。因此,本例的自相关运算部210生成表示脉波分析数据MKDf的自相关系数的自相关数据RD。
其次,脉象数据生成部220根据自相关数据RD生成脉象数据ZD。这时,脉象数据生成部220对自相关数据RD进行运算处理,将该运算结果与设定的能够特定脉象的阈值比较并生成脉象数据ZD。当将这样得到的脉象数据ZD供给显示部230时,便显示所谓弦脉、平脉、滑脉的文字,由此,被测者和医生等第3者便能认识脉象。《9第9实施例》
第9实施例是将使用了小波的体动除去用于第8实施例的脉波诊断装置1的情况。
图36第9实施例的脉波诊断装置1的方框图。在该例中,第1小波变换部243和第2小波变换部245与上述第8实施例的情况一样,如图35那样构成。因此,第1、第2小波变换部243、245生成只与所注意的频率区对应的脉波分析数据MKDf和体动分析数据TKDf。此外,体动成分除去部240从脉波分析数据MKDf中减去体动分析数据TKDf并生成体动除去分析数据MKDaj。
这样,当生成体动除去分析数据MKDfaj时,与第8实施例一样,自相关运算部210根据体动除去分析数据MKDfaj生成表示自相关系数的自相关数据RD。其次,脉象数据生成部220根据自相关数据RD生成脉象数据ZD。这时,脉象数据生成部220对自相关数据RD进行运算处理,将该运算结果与设定的能够特定脉象的阈值比较并生成脉象数据ZD。当将这样得到的脉象数据ZD供给显示部230时,便显示所谓弦脉、平脉、滑脉的文字,由此,被测者和医生等第3者便能认识脉象。
这样,若按照第9实施例,着眼于某个频率区,进行小波变换来除去体动,所以,没有必要对不同频率区的小波变换的结果进行比较。因此,可以省去第1、第2频率修正部244、246。此外,因为从体动除去分析数据MKDfaj直接生成自相关数据RD,所以能够省略逆小波变换部247。《10上述各实施例的外观结构》
下面,说明上述第1~第9实施例的脉波诊断装置的几个外观结构例子。《10.1手表式A》
首先,参照图37A~图37C说明上述各实施例的手表式脉波诊断装置1的结构例。
如参照图37A和图37B所示,脉波诊断装置1主要由具有手表结构的装置本体100、与该装置本体100连接的电缆101和设在该电缆101的前端的脉波检测部10构成。
其中,装置本体100装有表带102。详细地说,表带102缠绕在被测者的左腕,其方向对着装置本体100的相当于手表的12点钟的方向,另一端在6点钟的方向上固定。
在装置本体100的6点钟的方向上还设有连接部103。该连接部103安装有可以自由装卸的连接件104,该连接件又是电缆101的端部。
再有,若把该连接件104拿掉,则如图37C所示,在连接部103上,除了与电缆101连接的接点111、112之外,还设有用于数据传送的LED113、光电管114。
另一方面,脉波检测部10如图37B所示,套在被测者的食指的根部,利用传感器固定带11来遮光。这样,当把脉波检测部10套在手指根部时,因电缆可以很短,故不妨碍活动。此外,若测量从手掌到手指的体温分布,与寒冷时手指尖的温度明显低的情况相比,指根的温度因血流量稳定而并不低。因此,若将脉波检测部10装在指根部,即使在冷天外出,也能正确地检测脉波波形。
此外,在装置本体100的表面设有由液晶面板形成的显示部110。该显示部110具有段显示区和点阵显示区等,显示当前时间和诊断内容。即,显示部110与各实施例中的告知部80或显示部230对应。
另一方面,在装置本体100的内部装有未图示的加速度传感器,检测由被测者手腕的摆动和身体上下运动引起的体动。即,该加速度传感器与各实施例的体动检测部20对应。
此外,在装置本体100的内部设有控制各种运算和变换的CPU(省略图示),在装置本体100的外沿分别设有用于指示各种操作的按钮开关SW1、SW2。《10.1.1脉波检测部的详细结构》
下面,参照图38说明脉波检测部10的构成。
如图所示,脉波检测部10由LED12、光电管13等构成。当开关SW变成on状态加上电源电压时,从LED12发出照射光。该照射光经被测者的血管和组织反射后,由光电管13接收该反射光。从而,将光电管13的光电流变换成电压后,作为脉波检测部10的信号输出。
这里,LED的发光波长选择在血液中的血红蛋白的吸收波长峰值的附近。因此,受光电平与血流量对应变化。因此,通过检测受光电平,可以检测脉波波形。
此外,作为LED12,InGaN族(铟—镓—氮族)的蓝色LED很适合。蓝色LED的发光光谱例如在450nm处有发光峰值,其发光波长范围在从350nm到600nm的范围内。作为与具有这样的发光特性的LED对应的光电管13,本实施例使用GaAsP族(镓—砷—磷族)。该光电管13的受光波长区,例如,其主要灵敏度区是在从300nm到600nm的范围内,300nm以下也属于灵敏度区。
若将这样的蓝色LED和光电管组合,在其重合区、即从300nm到600nm的波长范围内检测脉波,具有以下优点。
首先,在包含在外界光的光中,波长区为700nm以下的光难以透过手指的组织,所以,即使外界光照射传感器固定带没有覆盖的手指部分,也不能经手指组织到达光电管13,只有对检测没有影响的波长区的光方可到达光电管13。另一方面,比300nm长的波长区的光因几乎被皮肤表面全部吸收,故即使使受光波长区在700nm以下,实质的受光波长区是从300nm~700nm。因此,即使手指没有被大部覆盖,也可以抑制外界光的影响。此外,血液中血红蛋白对300nm到700nm波长的光吸收系数大,与对880nm波长的光的吸收系数相比,大几倍到几百倍以上。因此,若象该例那样,将与血红蛋白的吸收特性对应、吸收特性强的波长区(从300nm到700nm)的光作为检测光使用,其检测峰值与血流量变化对应,灵敏度常常变化,所以,可以提高随血流量变化的脉波波形MH的信噪比。《10.2手表式B》
下面,参照图39A和图39B说明手表式脉波诊断装置1的另一个结构例。在该结构中,不是利用LED和光电管等以光电的方式检测被测者的脉波波形,而是使用压力传感器进行检测。
如图39A所示,脉波诊断装置1设有一付表带102、102,在其一边的紧固件120的紧固侧设有压力传感器130的弹性橡胶131,向外突起。具有紧固件120的表带102是将供给压力传感器的检测信号的FPC(Falexible Printed Circuit)底板用软塑料包起来的结构(省略详细图示)。
此外,使用时,如图39B所示,设在紧固件120的弹性橡胶131应放在桡骨动脉140的附近,手表结构的脉波诊断装置1缠绕在被测者的左手腕150上。因此,可以经常地检测脉波。再有,至于缠绕的方式与通常手表的使用状态没有两样。
当把这样的弹性橡胶131按压在被测者的桡骨动脉140的附近时,该动脉的血流变动(即脉波)便通过弹性橡胶131传递到压力传感器130,压力传感器130将它作为血压检测出来。
此外,作为脉波检测部10的另一个例子,也可以象特愿平5-192620所记载的那样,将罗口带套在手指上,给罗口带持续加压,检测脉波波形的高次谐波成分,由此,抽出潮波和重博波的特征。《10.3项圈式》
还考虑将各实施例的脉波诊断装置1做成图40所示那样的项圈形式。
在该图中,压力传感器130设在电缆101的前端,例如,如图41所示那样,使用粘性带170贴在被测者的颈动脉处。此外,在图40中,在具有中空部、形状象胸针的装置本体100内装有该装置的主要部件,同时,在其前面还设有显示部110、开关SW1、SW2。再有,电缆101的一部分埋在锁链160内,将压力传感器130输出的信号MH送给装置本体100。《10.4眼镜式》
作为上述各实施例的脉波诊断装置1的形态例,进而考虑图42所示那样的眼镜形式。
如该图所示,装置本体分成盒子100a和盒子100b,分别装在眼镜腿181上,经埋在眼镜腿181内的引线将它们互相电连接。在盒子100a的镜片182一侧的侧面装有液晶面板183,同时,在该侧面的一端,以给定的角度固定一面镜子184。此外,在盒子100a内还装有包含光源(省略图示)的液晶面板183的驱动电路和用来作成显示数据的电路,它们构成显示部110。从该光源发出的光经液晶面板183被镜子184反射,投射在镜片182上。此外,在盒子100b内装有装置的主要部分,在其上面,设有开关SW1、SW2。
另一方面,压力传感器130经电缆101与盒子100b电连接,和项圈式结构的情况一样,贴在颈动脉处。再有,连接盒子100a和盒子100b的引线也可以沿眼镜腿181布设。此外,在该例中,虽然是将装置本体分成盒子100a和盒子100b的结构,但也可以使它们一体化,做成一个盒子。进而,镜子183也可以做成角度可调的可动式。《10.5卡片式》
作为另一个形态例子,还考虑图43所示那样的卡片形式。该卡片式的装置本体100例如可以放在被测者的左胸口袋里。压力传感器130经电缆101与装置本体100电连接,与项圈式和眼镜式的情况一样,贴在被测者的颈动脉处。《10.6计步器式》
作为另一个形态例子,进而考虑图44A所示那样的计步器形式。该计步器的装置本体100如图44B所示,装在被测者的腰带191上。压力传感器130经电缆101与装置本体100电连接,利用粘性带将其固定在被测者的股关节的大腿动脉处,进而,利用护身带192将其保护起来。这时,最好采取将电缆101缝在衣服内等措施,以免妨碍被测者的日常生活。《11第10实施例》《11.1第10实施例的理论根据》
图47示出本发明者为了记录桡骨动脉的动脉压波形而使用的结构。该图所示的连续血压监视装置580(コ-リン社制CBM-2000)的构成包括臂带血压测定部582、桡骨动脉传感器部584、它们的控制部586和与控制部586连接的个人计算机588。
图48是表示用这样的装置测定的典型的动脉压波形、即桡骨动脉的血压波形的图。该图所示的前面所述那样的动脉血压波形通常包括具有最高峰的驱出波(ejection wave)、具有次高峰的潮波(tidalwave)和第3峰的重博波(dicrodic wave)。而且,驱出波的峰值与收缩血压BPsys对应。扩张血压BPdir与血压波形中的最低血压对应。而且,把收缩血压BPsys和扩张血压BPdir的压差称为脉压ΔBP。进而,平均血压BPmean是将血压波形积分求时间的平均得到的。
在图所示的连续血压监视装置580中,臂带血压测定部582用来测量收缩血压BPsys和扩张血压BPdir。桡骨动脉传感器部584检测与桡骨动脉的血压波形对应的脉波波形,利用臂带血压测定部582测定的收缩血压BPsys和扩张血压BPdir进行校正,得到动脉压波形。
本发明者利用上述连续血压监视装置580对74名22岁到46岁的健康成人测定了桡骨动脉的动脉压波形,测定是在空腹时、经过15分钟的安静后在座位上进行的。
图49~图51是表示将通过这样的测定得到的数据划成图形的结果的图。即,图49表示平均血压BPmean和扩张血压BPdir的关系有很强的线性关系,它们之间的相关系数r=0.95,具有很强的相关关系。图50表示平均血压BPmean和收缩血压BPsys的关系也有很强的线性关系(相关系数r=0.87)。图51表示收缩血压BPsys和扩张血压BPdir的压差、即脉压ΔBP与收缩血压BPsys之间存在很强的线性关系(相关系数r=0.86)。因此可知,通过平均血压BPmean和脉压ΔBP可以反映血压的状态。《11.2血压监视装置的构成》
图52是表示本实施例的血压监视装置500的构成的方框图。如该图所示,血压监视装置500的构成包括动脉压波形检测部504、血压变换部516、平均血压算出部508、脉压算出部512、血压判定信息存储部520、血压判定部524和输出部528。
动脉压检测部504连续测定动脉的血压并检测出动脉压波形。动脉压波形检测部504例如可以由图47所示的桡骨动脉传感器部584、图37b所示的脉波检测部10、图39A、图40、图41、图42、图43和图44A所示的压力传感器130形成。
血压变换部516将动脉压波形检测部504检测出的动脉压波形变换成与心脏的高度对应位置上的动脉压波形、即心脏位动脉压波形。例如,血压变换部516可以构成为将动脉压波形检测部504检测出的动脉压波形作为通过由臂带血压计测定的最高血压和最低血压校正过的动脉压波形。
平均血压算出部508根据动脉压波形或心脏位动脉压波形算出平均血压BPmean
脉压算出部512根据动脉压波形算出最高血压(收缩血压BPsys)和最低血压(扩张血压BPdir)的压差、即脉压ΔBP。
血压判定信息存储部520预先存储血压判定信息。血压判定信息例如是对平均血压和脉压来说、作为高血压和正常血压的临界点的阈值血压值和作为低血压和正常血压的临界点的阈值血压值的信息。再有,平均血压和脉压与收缩血压BPsys(最高血压)和扩张血压BPdir(最低血压)有很高的相关,能够作为表示血压状态的指标,这一点与前面所述的图49~图51所示一样。
血压判定部524根据平均血压BPmean和脉压ΔBP中的至少一个和血压判定信息存储部520存储的血压判定信息进行血压判定。例如,利用所得到的平均血压BPmean和脉压ΔBP中的至少一个和预先存储的血压判定信息进行例如高血压、低血压、正常等的判定。
输出部528输出与平均血压BPmean对应的信息、与脉压ΔBP对应的信息和与血压判定对应的信息当中的至少一个信息。输出部528可以将与平均血压对应的信息、与脉压对应的信息和与血压判定对应的信息当中的至少一个信息例如作为数值、图形等向液晶显示装置、CRT、打印机等输出,也可以作为与这些信息对应的电压或数字信息等输出。《11.3血压监视装置的动作》
根据图52说明本实施例的血压监视装置500的动作。
首先,动脉压波形检测部504利用前述的压力传感器130等连续测定动脉的血压并检测出动脉压波形。
接着,血压变换部516将动脉压波形检测部504检测出的动脉压波形变换成与心脏的高度对应位置的动脉压波形、即心脏位动脉压波形。
其次,平均血压算出部508根据从动脉压波形检测部504或血压变换部516输出的动脉压波形或心脏位动脉压波形算出平均血压BPmean
与此平行,脉压算出部512根据从动脉压波形检测部504或血压变换部516输出的动脉压波形或心脏位动脉压波形算出最高血压(收缩血压BPsys)和最低血压(扩张血压BPdir)的压差、即脉压ΔBP。
接下来,将平均血压BPmean和脉压ΔBP中的至少一个输入血压判定部524,血压判定部524根据这些数据和血压判定信息存储部520存储的血压判定信息进行血压判定,例如判定是高血压、低血压还是正常等。
而且,利用输出部528输出与平均血压BPmean对应的信息、与脉压ΔBP对应的信息和与血压判定对应的信息当中的至少一个信息。输出部例如包括液晶显示装置、CRT、打印机等,可以将这些信息作为数值、图形等显示,也可以作为与这些信息对应的电压或数字信息等输出。
这样一来,若按照本发明的血压监视装置500,根据动脉压波形检测部504检测出的动脉压波形算出平均血压或脉压,可以对它们进行监视。此外,血压监视装置500利用所得到的平均血压BPmean和脉压ΔBP中的至少一个和预先存储的血压判定信息进行例如高血压、低血压、正常等的判定。《12第11实施例》《12.1第11实施例的理论根据》
本发明者使用根据在《第10实施例的理论根据》拦中说明过的实验所得到的数据进一步进行了下面的实验。
即,参考中国医学的脉波形状的形状模型(神户中医学研究会:中医临床的舌诊和脉诊;费兆馥(主编):中国脉诊研究,上海中医学院出版社(1991)),将得到的桡骨动脉压波形区别为图45A、图45B、图45C所示的典型的动脉压波形、即平脉、滑脉和弦脉。
结果知道,利用重博和重博波峰之间的血压差、即重博波高ΔBPD(参照图48)能够象图53所示那样人为地分类成滑脉、平脉和弦脉。即,滑脉的ΔBPD是11±4mmHg,平脉的ΔBPD是7±2mmHg,弦脉的ΔBPD是3±1mmHg,并确认在3个群类之间有1%的人为判别误差。
此外,图54示出重博压差比BPDd/ΔBP与重博波高ΔBPD之间具有线性关系,该重博压差比是重博血压和最低血压的压差、即重博压差BPDd(参照图48)与最高血压和最低血压的压差、即脉压ΔBP的比值。从该图可知,它们之间的相关系数是-0.86,具有很强的相关关系。因此可知,重博压差比BPDd/ΔBP也可以作为对中国医学脉诊中的典型的脉波形状、即滑脉、平脉、弦脉进行人为分类的指标。
进而,图55示出重博血压和最高血压的压差、即驱出扩张压ΔBPp与重博波高ΔBPD之间具有线性关系。从该图可知,它们之间的相关系数是0.77,具有很强的相关关系。因此可知,也可以将驱出扩张压ΔBPp作为对中国医学脉诊中的典型的脉波形状、即滑脉、平脉、弦脉进行人为分类的指标。
而且,图45A、图45B和图45C表示相对分类成中国医学脉诊中的滑脉、平脉、弦脉的脉波形状的典型的动脉压波形,从这些图中可知,平均血压BPmean对脉压ΔBP的比、即平均血压脉压比BPmean/ΔBP按滑脉、平脉、弦脉的顺序增大。因此,利用平均血压脉压比BPmean/ΔBP能够对脉波形状、即滑脉、平脉、弦脉进行分类。《12.2脉波形状监视装置的构成》
图56是表示本实施例的脉波形状监视装置540的构成的方框图。如该图所示,脉波形状监视装置540的构成包括动脉压波形检测部504、重博波高算出部544、重博压差比算出部548、平均血压脉压比算出部552、驱出扩张压算出部554、脉波形状判定信息存储部556、脉波形状判定部560和输出部564。再有,图56虽未示出,也可以在动脉压波形检测部504的后级设置重博波高算出部544、重博压差比算出部548和在平均血压脉压比算出部552的前级设置血压变换部。
动脉压波形检测部504连续测定动脉血压并检测出动脉压波形。动脉压波形检测部504例如可以由图47所示的桡骨动脉传感器部584、图37B所示的脉波检测部10、图39A、图40、图41、图42、图43和图44A所示的压力传感器130形成。
血压变换部将动脉压波形检测部检测出的动脉压波形变换成与心脏的高度对应位置上的动脉压波形、即心脏位动脉压波形。
重博波高算出部544算出从动脉压波形或心脏位动脉压波形得到的重博(dicrotic notch)和重博波(dicrotic wave)的峰值之间的血压差、即重博波高ΔBPD(参照图48)。
重博压差比算出部548算出从动脉压波形或心脏位动脉压波形得到的重博(dicrotic notch)血压和最低血压(扩展血压BPdir)的压差、即重博压差BPDd与最高血压和最低血压的压差、即脉压ΔBP的比值、即重博压差比BPDd/ΔBP。
平均血压脉压算出部552算出从动脉压波形或心脏位动脉压波形得到的平均血压BPmean与最高血压和最低血压的压差、即脉压ΔBP的比值、即平均血压脉压比BPmean/ΔBP。
驱出扩张压ΔBPp算出部554算出从动脉压波形或心脏位动脉压波形得到的重博(dicrotic notch)和收缩血压BPsys之间的血压差、即驱出扩张压ΔBPp(参照图48)。
脉波形状判定信息存储部556预先存储脉波形状判定信息。脉波形状判定信息例如是对重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、或驱出扩张压ΔBPp来说、作为各脉波形状、例如中国医学诊脉中的滑脉、平脉、弦脉的临界值的阈值信息。再有,这些脉波形状可以利用重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、或驱出扩张压ΔBPp进行分类,这一点与前面提到的图53、图54、图45A、图45B、图45C、图55所示的一样。
脉波形状判定部560根据重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、或驱出扩张压ΔBPp中的至少一个和脉波形状判定信息进行脉波形状判定。因此,利用重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、或驱出扩张压ΔBPp中的至少一个和预先存储的脉波形状判定信息可以进行例如中国医学诊脉中的滑脉、平脉、弦等的判定。
输出部564输出与重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、驱出扩张压ΔBPp对应的信息和与脉波形状判定的信息当中的至少一个信息。输出部528可以例如作为数值、图形等向液晶显示装置、CRT、打印机等输出,也可以作为与这些信息对应的电压或数字信息等输出。《12.3脉波形状监视装置的动作》
根据图55说明本实施例的脉波形状监视装置540的动作。
首先,动脉压波形检测部504利用前述的压力传感器130等连续测定动脉的血压并检测出动脉压波形。
接着,血压变换部将动脉压波形检测部504检测出的动脉压波形变换成与心脏的高度对应位置的动脉压波形、即心脏位动脉压波形。再有,在下面的步骤中,当使用动脉压波形时,可以省略该步骤。
其次,重博波高算出部544算出从动脉压波形或心脏位动脉压波形得到的重博(dicrotic notch)和重博波(dicrotic wave)的峰值之间的血压差、即重博波高ΔBPD(参照图48)。
重博压差比算出部548与重博波高ΔBPD的算出平行,算出从动脉压波形或心脏位动脉压波形得到的重博(dicrotic notch)血压和最低血压(扩展血压BPdir)的压差、即重博压差BPDd与最高血压和最低血压的压差、即脉压ΔBP的比值、即重博压差比BPDd/ΔBP。
平均血压脉压算出部552与上述动作平行,算出从动脉压波形或心脏位动脉压波形得到的平均血压BPmean与最高血压和最低血压的压差、即脉压ΔBP的比值、即平均血压脉压比BPmean/ΔBP。
驱出扩张压算出部554同样与上述动作平行,算出从动脉压波形或心脏位动脉压波形得到的重博(dicrotic notch)和收缩血压BPsys之间的血压差、即驱出扩张压ΔBPp(参照图48)。
其次,脉波形状判定部560根据重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、或驱出扩张压ΔBPp中的至少一个和脉波形状判定信息进行脉波形状、例如中国医学诊脉中的滑脉、平脉、弦等的判定。
输出部564输出与重博波高ΔBPD对应的信息、与重博压差比BPDd/ΔBP对应的信息、与平均血压脉压比BPmean/ΔBP对应的信息、与驱出扩张压ΔBPp对应的信息和与脉波形状判定信息当中的至少一个信息。该输出部可以例如作为数值、图形等向液晶显示装置、CRT、打印机等输出,也可以作为与这些信息对应的电压或数字信息等输出。
这样一来,若按照本发明的脉波形状监视装置540,根据动脉压波形检测部504检测出的动脉压波形可以算出重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、或驱出扩张压ΔBPp中的至少一个。而且,脉波形状监视装置540利用重博波高ΔBPD、重博压差比BPDd/ΔBP、平均血压脉压比BPmean/ΔBP、或驱出扩张压ΔBPp中的至少一个和脉波形状判定信息可以判定脉波形状。进而,脉波形状监视装置540可以将与重博波高ΔBPD对应的信息、与重博压差比BPDd/ΔBP对应的信息、与平均血压脉压比BPmean/ΔBP对应的信息、与驱出扩张压ΔBPp对应的信息和与脉波形状对应的信息当中的至少一个信息例如作为数值、图形或电压输出。《13第12实施例》《13.1原理》
本发明者利用与第11实施例使用的脉波形状监视装置大致相同的构成监视了前面所述的驱出扩张压ΔBPp(从动脉压波形得到的重博与收缩血压之间的压差)和重博波高ΔBPD(重博和重博波峰之间的压差)。结果,确认了通过投入某种药品、例如降压剂可以象下述那样引起驱出扩张压的变化。
图63是表示把降压剂利心平在0分时投入、在前60分和后60分投入时分别测定平均血压BPmean、驱出扩张压ΔBPp、重博压差BPDd(重博血压与最低血压之间的压差)和重博波高ΔBPD(重博和重博波峰之间的压差)的结果的图。此外,在图64中,带有control的波形是投入降压剂前60分钟平均的1次博动的动脉压波形,带有NF(15min)的波形是投入降压剂后15分钟每10秒钟平均的1次博动的动脉压波形。从这些图可知,通过投入降压剂,平均血压BPmean上升了,驱出扩张压ΔBPp下降了,而且,重博波高ΔBPD增加了。本实施例涉及基于这样的的药理作用的观察的药理作用监视装置。《13.2药理作用监视装置的构成和作用》
药理作用监视装置570如图57所示,没有第11实施例的重博压差比算出部548和平均血压脉压比算出部552。而且,在药理作用监视装置570中,脉波形状判定部被换成药理作用判定部574,脉波形状判定信息存储部被换成药理作用判定信息存储部572。除此之外,药理作用监视装置570与第11实施例的脉波形状监视装置540的构成相同。
利用这样的的构成,药理作用监视装置570在动脉压波形检测部504检测出动脉压波形时,将该波形信息输入给驱出扩张压算出部554和重博波高算出部544。而且,驱出扩张压算出部554算出收缩血压(最高血压)和重博的差、即驱出扩张压,该驱出扩张压数据对药理作用判定部574和输出部564输出。此外,重博波高算出部544算出重博和重博波峰之间的血压差、即重博波高,并将该数据对药理作用判定部574和输出部564输出。
药理作用判定部574根据药理作用判定信息存储部572存储的药理作用判定信息,从输入的驱出扩张压的数据和重博波高的数据判定药理作用,将该结果输出给输出部564。输出部564通过例如LCD将驱出扩张压和重博波高的值和药理作用判定部的判定结果以图像信息、从扬声器来的声音信息和电压输出等形式告知给使用者。
再有,在药理作用监视装置570中,也可以不使用动脉压波形检测部504,而使用脉波检测部,检测从心脏博出的在血管中流动的血液的波、即脉波。再有,这时,因驱出扩张压算出部554和重博波高算出部544不能算出绝对血压,故利用算出用脉压、即收缩血压和扩张血压的压差归一化了的驱出扩张压比和重博波高比的驱出扩张压算出部和重博波高比算出部可以正确地进行数据的比较。
再有,上面示出了具有驱出扩张压算出部554和重博波高比算出部544的药理作用监视装置570,但作为药理作用监视装置的构成,也可以只具有驱出扩张压算出部554和重博波高比算出部544当中的任何一个。《14其它的变形例》
以上,说明了本发明的实施例,但本发明不限于前面所述的各实施例,在本发明的宗旨范围内或在与权利要求的范围均等的范围内,可以有各种变形实施例。
(1)在上述第1实施例中,潮波特征抽出部50利用(f2+f3+f4)/f1生成潮波特征数据TWD,重博波特征抽出部60利用(f5+f6+f7)/f1生成重博波特征数据DWD,但本发明不限于此,只要能抽出潮波、重博波的特征,不管什么样的谐波成分,都可以用来抽出特征。
(2)在上述第2实施例中,根据脉波分析数据特定潮波和重博波的时间位置,但本发明不限于此,只要能够特定各时间位置,用什么样的方法都可以。例如,将信号MHj微分,算出其过零点,示出信号MHj的峰点,据此来特定潮波、重博波的时间位置也可以。
(3)在上述第2实施例中,使用小波变换对信号MHj进行频谱分析,但也可以利用窗函数从信号MHj中抽出潮波和重博波,并对它门进行FFT处理。
(4)在第1和第2实施例中,以抽出潮波和重博波的特征为目的,利用FFT处理进行频谱分析(第1实施例),利用小波变换进行时间频谱分析(第2实施例),但本发明不限于此,只要能够根据潮波和重博波的振幅抽出它们的特征,用什么样的方法都可以。例如,也可以通过对信号MHj进行2次微分,强调信号MHj的波峰,并据此求出潮波和重博波的振幅,将各振幅作为特征抽出。
(5)在上述第6、第7实施例中,使用了第1频率校正部或第2频率校正部,但也可以省略这一部分。
(6)上述实施例所进行的小波变换和逆小波变换也可以使用滤波器组。小波变换使用的滤波器组的构成例如图58所示。图中,滤波器组由3级构成,其基本单元是高通滤波器1A和抽取滤波器1C、低通滤波器1B和抽取滤波器1C。高通滤波器1A和低通滤波器1B将信号分割成规定的频带,分别输出高频成分和低频成分。在该例中,因设想脉波数据MD的频带是0Hz~4Hz,故设定一级高通滤波器的通带为2Hz~4Hz,另一方面,设定一级低通滤波器的通带为0Hz~2Hz。此外,抽取滤波器1C每隔1个抽样拔掉1个数据。
当把这样生成的数据供给次级时,重复进行频带分割和间拔,最终得到将0Hz~4Hz的频带分成8份的数据M1~M8。
此外,高通滤波器1A和低通滤波器1B可以由其内部包含延迟元件(D触发器)的横向滤波器构成。但是,人的心率是在40~200的范围之内,脉波波形MH的基波频率随身体状态的不同时时刻刻都在变动。这时,若使分割的频带与基波频率同步改变,便可以得到动态地跟踪身体状态的信息。因此,也可以通过使脉波波形MH作为供给横向滤波器的时钟,自适应地改变分割的频带。
此外,在脉波分析数据MKD中,反映脉波波形MH的特征的代表性的频率成分是基波、2次谐波和3次谐波的各频率成分。因此,也可以使用滤波器组输出的数据M*1~M*8中的一部分来判定脉象。这时,若象上述那样使滤波器组与脉波波形MH同步,则可以省去高通滤波器1A、低通滤波器1B和抽取滤波器1C,使其构成简单。
其次,图58表示逆小波变换所用的逆滤波器组的构成例。图中,滤波器组由3级构成,其基本单元是高通滤波器2A和内插滤波器2C、低通滤波器2B和内插滤波器2C、加法器2D。高通滤波器2A和低通滤波器2B将信号分割成规定的频带,分别输出高频成分和低频成分。此外,内插滤波器2C每2次抽样内插1个样品。
这里,为了再现波形,有必要对图58所示的滤波器组和图59所示的滤波器组使用完全重构的滤波器组。这时,高通滤波器1A、2A和低通滤波器1B、2B的特性具有如下关系是有必要的。
H0(-Z)F0(Z)+H1(-Z)F1(Z)=0
H0(Z)F0(Z)+H1(-Z)F1(Z)=2Z-L
此外,高通滤波器2A和低通滤波器2B可以由其内部包含延迟元件(D触发器)的横向滤波器构成。再有,为了使小波变换部10使用的滤波器组与脉波波形MH的基波频率同步,使分割的频带可变,也可以在供给的时钟信号与脉波波形MH同步时将该时钟信号供给高通滤波器2A和低通滤波器2B。
(7)此外,在上述实施例中,将显示部30作为告知装置的一个例子进行了说明,但作为装置告知人的手段,可以列举以下说明的手段。将这些手段按五官进行分类是合适的。再有,这些手段当然不仅可以单独使用也可以将多个手段组合起来使用。而且,如以下说明的那样,例如,若使用非视觉的手段,则能够使视觉障碍者理解告知的内容,同样,若使用非听觉的手段,则能够使听觉障碍者理解告知的内容,能够构成对残疾人使用者也很方便的装置。
首先,作为听觉告知手段,有以心功能的分析诊断结果等的告知为目的或以警告为目的的装置等。例如,除了蜂鸣器之外,还有压电元件和扬声器。此外,作为特殊的例子,可以考虑让作为告知对象的人携带一部便携式无线呼叫接收机,当进行告知时从装置一侧对该便携式无线呼叫接收机进行呼叫。此外,在使用这些仪器进行告知时,多数情况是不单单为了告知,还想一起传达某些信息。这时,可以与要传达的信息的内容对应改变下面所示的音量等的电平。例如,音调、音量、音色和音乐的种类(曲子等)。
其次,使用视觉告知手段是以装置通知各种信息和测量结果为目的,或是为了进行警告。以这些为目的的手段可以考虑以下仪器。例如,显示装置、CRT(阴极射线管显示装置)、LCD(液晶显示装置)、打印机、X-Y绘图仪、和灯等。再有,作为特殊的显示装置有眼镜式幻灯机。此外,告知的形式考虑有以下几种。例如,在数值告知中,除数字显示和模拟显示之外,还有图形显示、显示色的浓淡、直接数值显示或在显示数值附加等级时的棒图显示、圆图和面孔表情图等。面孔表情图例如有如图60所示那样的图。
其次,触觉告知手段可以考虑以警告为目的来使用。这样的手段有如下几种。首先,有电刺激的手段,在手表等便携式仪器中设置有从里面向外突出的形状记忆合金,对该形状记忆合金通电进行刺激。此外,还有机械刺激的手段,其中又分为利用装在手表等便携式仪器中的可以从里面进出的突起物(例如不太尖的针等)进行刺激的结构和利用手表的装置本体100的振动进行刺激的结构等。
其次,嗅觉告知手段可以构成为,在装置内设置香料等的吐出机构,使香气与告知的内容对应,吐出与告知内容对应的香料。顺便提一下,作为香料等的吐出机构最好使用微型泵等。
再有,这些手段当然不仅可以单独使用也可以将多个手段组合使用。
(8)在上述各实施例中,作为脉波检测装置的一个例子举出脉波检测部10进行了说明,但本发明并不限于此,只要能检测出脉动,什么样的装置都可以。
例如,脉波检测部10利用了反射光,但也可以利用透过光。但波长区在700nm以下的光难以透过手指的组织。因此,当利用透过光时,从发光部照射600nm~1000nm波长的光,依组织→血管→组织的顺序使照射光透过,并检测出该透过光的光量变化。因为透过光受血液中血红蛋白的吸收,所以,通过检测透过光光量的变化可以检测出脉波波形。
这时,发光部最好使用InGaAs族(铟—镓—砷)和GaAs族(镓—砷)的激光发光二极管。但因600nm~1000nm波长的外界光容易透过组织,当外界光入射到受光部时会使脉波信号的信噪比降低。因此,可以从发光部照射已偏振了的激光,受光部经偏振光滤光器接受透过光。由此,可以不受外界光的影响,以很高的信噪比检测脉波信号。
这时,如图61所示,将发光部146设在固定件145的固定侧,而受光部147设在手表本体一侧。这时,从发光部146照射的光透过血管143后,通过桡骨148和尺骨142之间到达受光部147。再有,当使用透过光时,因必需使照射光透过组织,故若考虑组织的吸收,其波长最好是600nm~1000nm。
此外,图62是检测部位是耳朵的例子。夹持件190和夹持件191由弹簧192供给能量,能够以轴193为中心转动。此外,在夹持件190和夹持件191上设有发光部194和受光部195。当使用该脉波检测部时,由夹持件190和夹持件191将耳朵夹住后再检测脉波。
(9)在第9实施例中,利用体动检测部20检测体动波形TH,将脉波分析数据MKDf与体动分析数据TKDf比较,消除体动成分,算出自相关数据RD,并据此来特定脉象。但是,因为体动成分多数发生在比脉波波形MH的基波频率低的低频区内,所以,若将要监视的频率区选择在比脉波波形MH的基波频率高的高频区内,就可以省略体动检测部20、波形处理部21、判定部22、第2小波变换部45和体动成分除去部240。即,在图34所示的脉波诊断装置1中,若将要监视的频率区选择在比脉波波形MH的基波频率高的高频区内,那么即使有体动也能正确地特定脉象。
(10)在第10和第11实施例中,示出了将血压监视装置和脉波形状监视装置作为单独的装置形成的例子。但是,也可以使第11实施例示出的脉波形状监视装置包含在第10实施例示出的血压监视装置中。
(11)此外,在第11实施例中,示出了设有重博波高算出部544、重博压差比算出部548、平均血压脉压比算出部和驱出扩张压算出部的脉波形状监视装置的例子,但是,也可以具有这些算出部当中的至少任何一个。

Claims (61)

1.一种脉波诊断装置,其特征在于,包括:
检测身体的脉波波形的脉波检测装置;
从上述脉波波形中抽出潮波的特征并生成潮波特征信息的潮波特征抽出装置;
从上述脉波波形中抽出重博波的特征并生成重博波特征信息的重博波特征抽出装置;
根据上述潮波特征信息和上述重博波特征信息判定上述身体的脉象的脉象判定装置。
2.权利要求1所述的脉波诊断装置,其特征在于,上述潮波特征抽出装置根据上述潮波在时间域中的振幅变化生成潮波特征信息,上述重博波特征抽出装置根据上述重博波在时间域中的振幅变化生成重博波特征信息。
3.权利要求2所述的脉波诊断装置,其特征在于,上述潮波和上述重博波在时间域中的振幅变化通过将上述脉波波形对时间进行一次微分或2次微分算出来。
4.一种脉波诊断装置,其特征在于,包括:
检测身体的脉波波形的脉波检测装置;
对上述脉波波形进行频谱分析的频谱分析装置;
根据上述频谱分析装置的分析结果从上述脉波波形中抽出潮波的特征并生成潮波特征信息的潮波特征抽出装置;
根据上述频谱分析装置的分析结果从上述脉波波形中抽出重博波的特征并生成重博波特征信息的重博波特征抽出装置;
根据上述潮波特征信息和上述重博波特征信息判定上述身体的脉象的脉象判定装置。
5.权利要求4所述的脉波诊断装置,其特征在于,
上述潮波特征抽出装置先特定上述脉波波形中的潮波的持续时间,再根据该期间内的上述频谱分析装置的分析结果从上述脉波波形中抽出潮波的特征,并生成潮波特征信息,
上述重博波特征抽出装置先特定上述脉波波形中的重博波的持续时间,再根据该期间内的上述频谱分析装置的分析结果从上述脉波波形中抽出重博波的特征,并生成重博波特征信息。
6.权利要求4或5所述的脉波诊断装置,其特征在于,
上述频谱分析装置对上述脉波波形进行FFT处理。
7.权利要求4或5所述的脉波诊断装置,其特征在于,
上述频谱分析装置对上述脉波波形进行小波变换处理。
8.权利要求1至7任一项所述的脉波诊断装置,其特征在于,
具有告知装置,告诉利用上述脉象判定装置判定的脉象。
9.一种脉波诊断装置,其特征在于,包括:
从身体的检测部位检测脉波波形的脉波检测装置;
自相关运算装置,计算表示利用上述脉波检测装置检测出的上述脉波波形的自相关的自相关数据;
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
10.权利要求9所述的脉波诊断装置,其特征在于,
上述脉象数据生成装置通过将上述自相关数据与预先设定的阈值进行比较来生成上述脉象数据。
11.权利要求10所述的脉波诊断装置,其特征在于,
上述脉象数据生成装置具有最小值检测部和比较部,最小值检测部检测一个心博期间上述自相关数据的最小值,
比较部将上述最小值检测部检测出的最小值与上述阈值进行比较并生成上述脉象数据。
12.权利要求10所述的脉波诊断装置,其特征在于,
上述脉象数据生成装置具有最小值检测部和比较部,最小值检测部对在几个心博期间检测出上述自相关数据的各最小值进行平均后检测其平均最小值,
比较部将上述最小值检测部检测出的平均最小值与上述阈值进行比较并生成上述脉象数据。
13.权利要求(9)所述的脉波诊断装置,其特征在于,
上述脉象数据生成装置最好具有时间测量部、运算部和比较部,时间测量部将上述自相关数据与预定的阈值进行比较,测量上述自相关数据高于上述阈值的时间间隔和低于上述阈值的时间间隔,
运算部计算由上述时间测量部测量的上述时间间隔在一个心博期间所占的比例,
比较部通过将上述运算部计算的结果与预定的阈值进行比较来生成上述脉象数据。
14.权利要求13所述的脉波诊断装置,其特征在于,
上述运算部算出由上述时间测量部测量的上述时间间隔在一个心博期间所占的比例,并计算该算出结果的平均值。
15.权利要求9所述的脉波诊断装置,其特征在于,
上述脉象数据生成装置具有变化率运算部和变化率比较部,变化率运算部根据上述自相关数据检测上述自相关数据的变化率,
变化率比较部通过将上述变化率运算部检测出的变化率与预定的阈值进行比较来生成上述脉象数据。
16.权利要求5所述的脉波诊断装置,其特征在于,
上述变化率比较部最好检测上述变化率的最大值,通过将变化率最大值与上述阈值进行比较来生成上述脉象数据。
17.权利要求9所述的脉波诊断装置,其特征在于,
上述脉象数据生成装置具有最小值检测部、第1比较部、时间测量部、运算部和第2比较部,最小值检测部检测1个心博期间上述自相关数据的最小值,
第1比较部将上述最小值检测部检测出的最小值与预定的第1阈值进行比较,当上述最小值低于上述阈值时,生成表示弦脉的脉象数据,
时间测量部将上述自相关数据与预定的第2阈值进行比较,测量上述自相关数据高于上述第2阈值的时间间隔或低于上述第2阈值的时间间隔,
运算部计算由上述时间测量部测量的上述时间间隔在一个心博期间所占的比例,
第2比较部通过将上述运算部计算的结果与预定的第3阈值进行比较来生成表示平脉或滑脉的上述脉象数据。
18.权利要求9所述的脉波诊断装置,其特征在于,
上述脉象数据生成装置具有最小值检测部、第1比较部、变化率运算部和第2比较部,最小值检测部检测1个心博期间上述自相关数据的最小值,
第1比较部将上述最小值检测部检测出的最小值与上述阈值进行比较,当上述最小值低于上述阈值时,生成表示弦脉的脉象数据,
变化率运算部根据自相关数据检测上述自相关数据的变化率,
第2比较部通过将上述变化率运算部检测出的变化率与预定的阈值进行比较来生成表示平脉或滑脉的上述脉象数据。
19.权利要求17或18所述的脉波诊断装置,其特征在于,
上述自相关数据表示自相关系数,上述第1比较部的比较动作所用的上述第1阈值大约是0.25。
20.权利要求17所述的脉波诊断装置,其特征在于,
上述自相关数据表示自相关系数,上述所述时间测量部的比较动作所用的上述第2阈值设定在0.4到0.8的范围内。
21.权利要求9至20中任何一项所述的脉波诊断装置,其特征在于,
具有检测表示上述身体活动的体动波形的体动检测装置和体动成分除去装置,该体动成分除去装置根据上述体动波形生成上述脉波波形中的体动成分,从上述脉波波形中除去上述体动成分并生成体动除去脉波波形,
上述自相关运算装置不是根据上述脉波波形而是根据上述体动除去脉波波形计算表示自相关的自相关数据。
22.权利要求21所述的脉波诊断装置,其特征在于,
具有判定装置,根据上述体动检测装置检测出的上述体动波形判定上述身体有无活动,
上述体动成分除去装置在上述判定装置的判定结果是没有体动时停止体动除去动作,并输出上述脉波波形,而不输出上述体动除去脉波波形。
23.权利要求9至20中任何一项所述的脉波诊断装置,其特征在于,
具有第1小波变换装置、体动检测装置、第2小波变换装置、体动成分除去装置和逆小波变换装置,
第1小波变换装置对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动检测装置检测上述身体的活动并输出体动波形,
第2小波变换装置对上述体动检测装置检测出的上述体动波形进行小波变换并生成各频率区的体动分析数据,
体动成分除去装置从上述脉波分析数据中减去上述体动分析数据,并生成已除去了体动的体动除去脉波分析数据,
逆小波变换装置对上述体动除去脉波分析数据进行逆小波变换并生成体动除去脉波波形。
上述自相关运算装置不是根据上述脉波波形而是根据上述体动除去脉波波形计算表示自相关的自相关数据。
24.权利要求9至20中任何一项所述的脉波诊断装置,其特征在于,
具有小波变换装置、体动成分除去装置和逆小波变换装置,
小波变换装置对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动成分除去装置从上述脉波分析数据中除去与预定的体动对应的频率成分,并生成体动除去脉波分析数据,
逆小波变换装置对上述体动除去脉波分析数据进行逆小波变换并生成体动除去脉波波形。
上述自相关运算装置不是根据上述脉波波形而是根据上述体动除去脉波波形计算表示自相关的自相关数据。
25.一种脉波诊断装置,其特征在于,包括:
从身体的检测部位检测脉波波形的脉波检测装置;
小波变换装置,对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
自相关运算装置,对某频率区的上述脉波分析数据计算表示自相关的自相关数据,
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
26.一种脉波诊断装置,其特征在于,包括:
从身体的检测部位检测脉波波形的脉波检测装置;
第1小波变换装置,对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动检测装置,检测上述身体的活动并输出体动波形,
第2小波变换装置,对上述体动检测装置检测出的上述体动波形进行小波变换并生成各频率区的体动分析数据,
体动成分除去装置,从上述脉波分析数据中减去上述体动分析数据,并生成已除去了体动的体动除去脉波分析数据,
自相关运算装置,对某频率区的上述体动除去脉波分析数据计算表示自相关的自相关数据,
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
27.一种脉波诊断装置,其特征在于,包括:
从身体的检测部位检测脉波波形的脉波检测装置;
小波变换装置,对上述脉波检测装置检测出的上述脉波波形进行小波变换并生成各频率区的脉波分析数据,
体动成分除去装置,从上述脉波分析数据中除去与预定的体动对应的频率成分,并生成体动除去脉波分析数据,
自相关运算装置,对某频率区的上述体动除去脉波分析数据计算表示自相关的自相关数据,
脉象数据生成装置,根据上述自相关数据生成表示上述脉波波形的种类的脉象数据。
28.权利要求9至27中任何一项所述的脉波诊断装置,其特征在于,
具有告知装置,告诉利用上述脉象数据生成装置生成的上述脉象数据。
29.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压。
30.权利要求29所述的血压监视装置,其特征在于,
进而具有脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压。
31.权利要求30所述的血压监视装置,其特征在于,
进而具有血压变换部,将上述动脉压波形检测部检测出的动脉压波形变换成与心脏的高度对应的位置上的动脉压波形、即心脏位动脉压波形,
上述平均血压算出部根据上述心脏位动脉压波形算出平均血压,上述脉压算出部根据上述心脏位动脉压波形算出脉压。
32.权利要求30或31所述的血压监视装置,其特征在于,
进而具有预先存储血压判定信息的血压判定信息存储部和血压判定部,血压判定部根据上述平均血压和上述脉压当中的至少一个参数以及上述血压判定信息进行血压判定。
33.权利要求30至32中任何一项所述的血压监视装置,其特征在于,
进而具有输出部,输出与上述平均血压对应的信息、与上述脉压对应的信息和与上述血压判定对应的信息当中的至少一种信息。
34.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形计算最高血压和最低血压的压差、即脉压。
35.权利要求34所述的血压监视装置,其特征在于,
具有根据上述脉压算出最高血压的最高血压算出部和根据上述脉压及上述最高血压算出最低血压的最低血压算出部。
36.一种脉波形状监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
37.权利要求36所述的脉波形状监视装置,其特征在于,
进而具有重博压差比算出部,计算从上述动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
38.权利要求37所述的脉波形状监视装置,其特征在于,
进而具有平均血压脉压比算出部,计算从上述动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
39.权利要求38所述的脉波形状监视装置,其特征在于,
进而具有血压变换部,将上述动脉压波形检测部检测出的动脉压波形变换成与心脏的高度对应的位置上的动脉压波形、即心脏位动脉压波形,
上述重博波高算出部计算从上述心脏位动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
上述重博压差比算出部计算从上述心脏位动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
上述平均血压脉压比算出部计算从上述心脏位动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
40.权利要求36所述的脉波形状监视装置,其特征在于,
进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和脉波形状判定部,脉波形状判定部根据上述重博波高和上述脉波形状判定信息判定脉波形状。
41.权利要求37所述的脉波形状监视装置,其特征在于,
进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和脉波形状判定部,
脉波形状判定部根据上述重博波高、上述重博压差比和上述脉波形状判定信息判定脉波形状。
42.权利要求38所述的脉波形状监视装置,其特征在于,
进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和脉波形状判定部,
脉波形状判定部根据上述重博波高、上述重博压差比、上述平均血压脉压比和上述脉波形状判定信息判定脉波形状。
43.权利要求36所述的脉波形状监视装置,其特征在于,
进而具有输出部,输出与上述重博波高对应的信息和上述脉波形状对应的信息当中的至少一种信息。
44.权利要求37所述的脉波形状监视装置,其特征在于,
进而具有输出部,输出与上述重博波高对应的信息、与上述重博压差比对应的信息和与上述脉波形状对应的信息当中的至少一种信息。
45.权利要求38所述的脉波形状监视装置,其特征在于,
进而具有输出部,输出与上述重博波高对应的信息、与上述重博压差比对应的信息、与上述平均血压脉压比对应的信息和与上述脉波形状对应的信息当中的至少一种信息。
46.一种脉波形状监视装置,其特征在于,具有:
脉波检测部,检测身体的脉波波形;
重博压差比算出部,计算从上述脉波波形得到的重博压和最低压的压差、即重博压差与最高压和最低压的压差、即脉压的比值、即重博压差比。
47.一种脉波形状监视装置,其特征在于,具有:
脉波检测部,检测身体的脉波波形;
平均压脉压比算出部,计算从上述脉波波形得到的平均压与最高压和最低压的压差、即脉压的比值、即平均压脉压比。
48.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压;
重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
49.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压;
重博压差比算出部,计算从上述动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
50.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
平均血压算出部,根据上述动脉压波形计算平均血压;
平均血压脉压比算出部,计算从上述动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
51.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压;
重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
52.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压;
重博压差比算出部,计算从上述动脉压波形得到的重博的血压和最低血压的压差、即重博压差与最高血压和最低血压的压差、即脉压的比值、即重博压差比。
53.一种血压监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
脉压算出部,根据上述动脉压波形算出最高血压和最低血压的压差、即脉压;
平均血压脉压比算出部,计算从上述动脉压波形得到的平均血压与最高血压和最低血压的压差、即脉压的比值、即平均血压脉压比。
54.一种脉波形状监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
驱出扩张压算出部,计算从上述动脉压波形得到的收缩期血压和重博之间的血压差、即驱出扩张压。
55.权利要求54所述的脉波形状监视装置,其特征在于,
进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和根据上述驱出扩张压判定脉波形状的脉波形状判定部。
56.一种脉波形状监视装置,其特征在于,具有:
脉波检测部,检测身体的脉波波形;
驱出扩张压比算出部,计算从上述动脉压波形得到的重博压和收缩压之间的压差、即驱出扩张压与收缩压和扩张压之间的压差、即脉压的比,即驱出扩张压比。
57.权利要求56所述的脉波形状监视装置,其特征在于,
进而具有预先存储脉波形状判定信息的脉波形状判定信息存储部和根据上述驱出扩张压比判定脉波形状的脉波形状判定部。
58.一种药理作用监视装置,其特征在于,具有:
动脉压波形检测部,连续地测定动脉的血压并检测动脉压波形;
驱出扩张压算出部,计算从上述动脉压波形得到的收缩期血压和重博之间的血压差、即驱出扩张压。
59.权利要求(58)所述的药理作用监视装置,其特征在于,
进而具有重博波高算出部,计算从上述动脉压波形得到的重博和重博波峰之间的血压差、即重博波高。
60.一种药理作用监视装置,其特征在于,具有:
脉波检测部,从身体检测脉波波形;
驱出扩张压比算出部,计算从上述脉波波形得到的重博压和收缩压之间的压差、即驱出扩张压与收缩压和扩张压之间的压差、即脉压的比值、即驱出扩张压比。
61.权利要求60所述的药理作用监视装置,其特征在于,
进而具有重博波高比算出部,计算从上述脉波波形得到的重博和重博波峰之间的压差、即重博波高与收缩压和扩张压之间的压差、即脉压的比值、即重博波高比。
CNB988035162A 1997-11-20 1998-11-20 脉波诊断装置 Expired - Fee Related CN1326489C (zh)

Applications Claiming Priority (10)

Application Number Priority Date Filing Date Title
JP320149/1997 1997-11-20
JP32014997 1997-11-20
JP320149/97 1997-11-20
JP321768/97 1997-11-21
JP32176897 1997-11-21
JP321768/1997 1997-11-21
JP21349498 1998-07-13
JP213494/1998 1998-07-13
JP213494/98 1998-07-13
PCT/JP1998/005259 WO1999026529A1 (fr) 1997-11-20 1998-11-20 Appareil de diagnostic lie aux formes d'impulsion, dispositif de controle de pression arterielle, dispositif de controle de la configuration des formes d'impulsion et dispositif de controle d'effet pharmacologique

Publications (2)

Publication Number Publication Date
CN1251027A true CN1251027A (zh) 2000-04-19
CN1326489C CN1326489C (zh) 2007-07-18

Family

ID=27329505

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB988035162A Expired - Fee Related CN1326489C (zh) 1997-11-20 1998-11-20 脉波诊断装置

Country Status (8)

Country Link
US (1) US6293915B1 (zh)
EP (2) EP0968681B1 (zh)
JP (1) JP3843462B2 (zh)
CN (1) CN1326489C (zh)
DE (2) DE69835843T2 (zh)
HK (1) HK1027496A1 (zh)
TW (1) TW425278B (zh)
WO (1) WO1999026529A1 (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100396233C (zh) * 2002-07-09 2008-06-25 健资国际私人有限公司 监测血压的方法和设备
CN102727191A (zh) * 2011-03-29 2012-10-17 精工爱普生株式会社 脉波计以及信号处理方法
CN102955891A (zh) * 2011-08-30 2013-03-06 华晶科技股份有限公司 具有资料检测功能的电子装置及其资料检测方法
CN104254275A (zh) * 2012-02-22 2014-12-31 阿克拉里斯医疗有限责任公司 生理信号检测装置和系统
CN104490373A (zh) * 2014-12-17 2015-04-08 辛勤 脉搏信号的判断方法、判断装置以及生理参数测量设备
CN106539570A (zh) * 2016-07-04 2017-03-29 悦享趋势科技(北京)有限责任公司 定位潮波的方法及装置
CN107708533A (zh) * 2015-07-01 2018-02-16 浜松光子学株式会社 粘弹特性获取装置、粘弹特性获取方法、粘弹特性获取程序及存储该程序的存储介质
CN107735018A (zh) * 2015-07-01 2018-02-23 浜松光子学株式会社 血压比计算装置、血压比计算方法、血压比计算程序及存储该程序的存储介质
CN108289624A (zh) * 2015-11-27 2018-07-17 浜松光子学株式会社 血压信息计算装置、血压信息计算方法、血压信息计算程序及存储该程序的存储介质
CN108403094A (zh) * 2018-03-23 2018-08-17 中国科学院微电子研究所 一种识别脉搏波波峰的方法

Families Citing this family (148)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6760416B1 (en) 1999-04-30 2004-07-06 Hewlett-Packard Development Company, L.P. Metering data transfer between end systems
US6859527B1 (en) 1999-04-30 2005-02-22 Hewlett Packard/Limited Communications arrangement and method using service system to facilitate the establishment of end-to-end communication over a network
US6526131B1 (en) 1999-04-30 2003-02-25 Hewlett-Packard Company Initiation of communication between network service system and customer-premises equipment
US6684336B1 (en) 1999-04-30 2004-01-27 Hewlett-Packard Development Company, L.P. Verification by target end system of intended data transfer operation
US9248306B2 (en) 1999-09-30 2016-02-02 Physio-Control, Inc. Pulse detection apparatus, software, and methods using patient physiological signals
US20030109790A1 (en) * 2001-12-06 2003-06-12 Medtronic Physio-Control Manufacturing Corp. Pulse detection method and apparatus using patient impedance
US20040039419A1 (en) * 1999-09-30 2004-02-26 Stickney Ronald E. Apparatus, software, and methods for cardiac pulse detection using a piezoelectric sensor
US6440082B1 (en) 1999-09-30 2002-08-27 Medtronic Physio-Control Manufacturing Corp. Method and apparatus for using heart sounds to determine the presence of a pulse
US6592528B2 (en) 2000-03-23 2003-07-15 Seiko Epson Corporation Biological information evaluation apparatus
JP2002172095A (ja) * 2000-12-06 2002-06-18 K & S:Kk 脈波測定装置
US6697658B2 (en) 2001-07-02 2004-02-24 Masimo Corporation Low power pulse oximeter
FR2827496B1 (fr) 2001-11-09 2004-06-04 Univ Clermont Auvergne Procedes et equipements destines a l'analyse de signaux biologiques representatifs des variations de la pression intracranienne et de la pression sanguine
SE0102920D0 (sv) * 2001-08-31 2001-08-31 St Jude Medical Medical device
EP1297784B8 (en) * 2001-09-28 2011-01-12 CSEM Centre Suisse d'Electronique et de Microtechnique SA - Recherche et Développement Method and device for pulse rate detection
WO2003040731A1 (en) * 2001-11-06 2003-05-15 Wireless Republic Group Apparatus and method for capturing and working acceleration, and application thereof, and computer readable recording medium storing programs for realizing the acceleration capturing and working methods
TW524670B (en) * 2002-04-01 2003-03-21 Ind Tech Res Inst Non-invasive apparatus system for monitoring autonomic nervous system and uses thereof
US20060165596A1 (en) * 2002-08-23 2006-07-27 Nachiket Kharalkar Method and apparatus for noninvasively evaluating endothelial function
US20070225614A1 (en) * 2004-05-26 2007-09-27 Endothelix, Inc. Method and apparatus for determining vascular health conditions
US20040116969A1 (en) * 2002-08-26 2004-06-17 Owen James M. Pulse detection using patient physiological signals
US20040039420A1 (en) * 2002-08-26 2004-02-26 Medtronic Physio-Control Manufacturing Corp. Apparatus, software, and methods for cardiac pulse detection using accelerometer data
US20050096557A1 (en) * 2003-01-08 2005-05-05 Frederick Vosburgh Noninvasive cardiovascular monitoring methods and devices
US8255029B2 (en) * 2003-02-27 2012-08-28 Nellcor Puritan Bennett Llc Method of analyzing and processing signals
JP5104842B2 (ja) * 2003-03-19 2012-12-19 セイコーエプソン株式会社 脈拍計、制御プログラムおよび記録媒体
US7455643B1 (en) 2003-07-07 2008-11-25 Nellcor Puritan Bennett Ireland Continuous non-invasive blood pressure measurement apparatus and methods providing automatic recalibration
US7292245B2 (en) * 2004-01-20 2007-11-06 Sensitron, Inc. Method and apparatus for time series graph display
US7442169B2 (en) * 2004-03-05 2008-10-28 Atcor Medical Pty Limited Methods of distinguishing between vasoconstriction and vasodilation as a cause of hypotension
US7039538B2 (en) * 2004-03-08 2006-05-02 Nellcor Puritant Bennett Incorporated Pulse oximeter with separate ensemble averaging for oxygen saturation and heart rate
EP2286721B1 (en) 2005-03-01 2018-10-24 Masimo Laboratories, Inc. Physiological Parameter Confidence Measure
US20070112275A1 (en) * 2005-08-15 2007-05-17 Cooke William H Medical Intervention Indicator Methods and Systems
US20070100220A1 (en) 2005-10-28 2007-05-03 Baker Clark R Jr Adjusting parameters used in pulse oximetry analysis
US20070225606A1 (en) * 2006-03-22 2007-09-27 Endothelix, Inc. Method and apparatus for comprehensive assessment of vascular health
US20080027330A1 (en) * 2006-05-15 2008-01-31 Endothelix, Inc. Risk assessment method for acute cardiovascular events
DE102006024459A1 (de) * 2006-05-24 2007-11-29 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Sensor, Verarbeitungseinrichtung, Verfahren und Computerprogramm zum Liefern einer Information über einen Vitalparameter eines Lebewesens
CN100413464C (zh) * 2006-05-26 2008-08-27 中国人民解放军空军航空医学研究所 在脉搏波法动脉血压连续测量中的脉搏波传导时间的获取方法和装置
US20080081963A1 (en) * 2006-09-29 2008-04-03 Endothelix, Inc. Methods and Apparatus for Profiling Cardiovascular Vulnerability to Mental Stress
JP2010535574A (ja) 2007-08-07 2010-11-25 カウンシル オブ サイエンティフィック アンド インダストリアル リサーチ 動脈のナンディパルス波形の定量的検出に有用な非侵襲装置「ナンディ・トラジーニ」
JP5555164B2 (ja) * 2007-09-19 2014-07-23 コーニンクレッカ フィリップス エヌ ヴェ 異常状態検出方法及び装置
CA2710994A1 (en) * 2007-12-26 2009-07-09 Nellcor Puritan Bennett Llc Historical trend icons for physiological parameters
US8827917B2 (en) * 2008-06-30 2014-09-09 Nelleor Puritan Bennett Ireland Systems and methods for artifact detection in signals
US7944551B2 (en) 2008-06-30 2011-05-17 Nellcor Puritan Bennett Ireland Systems and methods for a wavelet transform viewer
US8077297B2 (en) * 2008-06-30 2011-12-13 Nellcor Puritan Bennett Ireland Methods and systems for discriminating bands in scalograms
US8295567B2 (en) 2008-06-30 2012-10-23 Nellcor Puritan Bennett Ireland Systems and methods for ridge selection in scalograms of signals
US20090324033A1 (en) * 2008-06-30 2009-12-31 Nellcor Puritan Bennett Ireland Signal Processing Systems and Methods for Determining Slope Using an Origin Point
US8660799B2 (en) 2008-06-30 2014-02-25 Nellcor Puritan Bennett Ireland Processing and detecting baseline changes in signals
US20090326402A1 (en) * 2008-06-30 2009-12-31 Nellcor Puritan Bennett Ireland Systems and methods for determining effort
US8398556B2 (en) 2008-06-30 2013-03-19 Covidien Lp Systems and methods for non-invasive continuous blood pressure determination
US8437825B2 (en) 2008-07-03 2013-05-07 Cercacor Laboratories, Inc. Contoured protrusion for improving spectroscopic measurement of blood constituents
US8285352B2 (en) 2008-07-15 2012-10-09 Nellcor Puritan Bennett Llc Systems and methods for identifying pulse rates
US8358213B2 (en) 2008-07-15 2013-01-22 Covidien Lp Systems and methods for evaluating a physiological condition using a wavelet transform and identifying a band within a generated scalogram
US20100016676A1 (en) * 2008-07-15 2010-01-21 Nellcor Puritan Bennett Ireland Systems And Methods For Adaptively Filtering Signals
US8370080B2 (en) * 2008-07-15 2013-02-05 Nellcor Puritan Bennett Ireland Methods and systems for determining whether to trigger an alarm
US20100016692A1 (en) * 2008-07-15 2010-01-21 Nellcor Puritan Bennett Ireland Systems and methods for computing a physiological parameter using continuous wavelet transforms
US8226568B2 (en) 2008-07-15 2012-07-24 Nellcor Puritan Bennett Llc Signal processing systems and methods using basis functions and wavelet transforms
US8679027B2 (en) 2008-07-15 2014-03-25 Nellcor Puritan Bennett Ireland Systems and methods for pulse processing
US8761855B2 (en) 2008-07-15 2014-06-24 Nellcor Puritan Bennett Ireland Systems and methods for determining oxygen saturation
US8082110B2 (en) * 2008-07-15 2011-12-20 Nellcor Puritan Bennett Ireland Low perfusion signal processing systems and methods
US8506498B2 (en) 2008-07-15 2013-08-13 Nellcor Puritan Bennett Ireland Systems and methods using induced perturbation to determine physiological parameters
US8660625B2 (en) * 2008-07-15 2014-02-25 Covidien Lp Signal processing systems and methods for analyzing multiparameter spaces to determine physiological states
US8385675B2 (en) * 2008-07-15 2013-02-26 Nellcor Puritan Bennett Ireland Systems and methods for filtering a signal using a continuous wavelet transform
US8515509B2 (en) 2008-08-04 2013-08-20 Cercacor Laboratories, Inc. Multi-stream emitter for noninvasive measurement of blood constituents
JP2010048120A (ja) * 2008-08-20 2010-03-04 Olympus Corp 形状記憶合金アクチュエータシステム
US9314168B2 (en) 2008-09-30 2016-04-19 Nellcor Puritan Bennett Ireland Detecting sleep events using localized blood pressure changes
US8532751B2 (en) 2008-09-30 2013-09-10 Covidien Lp Laser self-mixing sensors for biological sensing
US8410951B2 (en) 2008-09-30 2013-04-02 Covidien Lp Detecting a signal quality decrease in a measurement system
US8696585B2 (en) * 2008-09-30 2014-04-15 Nellcor Puritan Bennett Ireland Detecting a probe-off event in a measurement system
US9301697B2 (en) 2008-09-30 2016-04-05 Nellcor Puritan Bennett Ireland Systems and methods for recalibrating a non-invasive blood pressure monitor
US9687161B2 (en) 2008-09-30 2017-06-27 Nellcor Puritan Bennett Ireland Systems and methods for maintaining blood pressure monitor calibration
US9011347B2 (en) 2008-10-03 2015-04-21 Nellcor Puritan Bennett Ireland Methods and apparatus for determining breathing effort characteristics measures
US9155493B2 (en) 2008-10-03 2015-10-13 Nellcor Puritan Bennett Ireland Methods and apparatus for calibrating respiratory effort from photoplethysmograph signals
US20100087714A1 (en) * 2008-10-03 2010-04-08 Nellcor Puritan Bennett Ireland Reducing cross-talk in a measurement system
US8216136B2 (en) 2009-03-05 2012-07-10 Nellcor Puritan Bennett Llc Systems and methods for monitoring heart rate and blood pressure correlation
US8364225B2 (en) * 2009-05-20 2013-01-29 Nellcor Puritan Bennett Ireland Estimating transform values using signal estimates
US20100298728A1 (en) * 2009-05-20 2010-11-25 Nellcor Puritan Bennett Ireland Signal Processing Techniques For Determining Signal Quality Using A Wavelet Transform Ratio Surface
US8444570B2 (en) * 2009-06-09 2013-05-21 Nellcor Puritan Bennett Ireland Signal processing techniques for aiding the interpretation of respiration signals
US20100324827A1 (en) * 2009-06-18 2010-12-23 Nellcor Puritan Bennett Ireland Fluid Responsiveness Measure
US20100331716A1 (en) * 2009-06-26 2010-12-30 Nellcor Puritan Bennett Ireland Methods and apparatus for measuring respiratory function using an effort signal
US8290730B2 (en) 2009-06-30 2012-10-16 Nellcor Puritan Bennett Ireland Systems and methods for assessing measurements in physiological monitoring devices
US9198582B2 (en) 2009-06-30 2015-12-01 Nellcor Puritan Bennett Ireland Determining a characteristic physiological parameter
US20100331715A1 (en) * 2009-06-30 2010-12-30 Nellcor Puritan Bennett Ireland Systems and methods for detecting effort events
US8636667B2 (en) * 2009-07-06 2014-01-28 Nellcor Puritan Bennett Ireland Systems and methods for processing physiological signals in wavelet space
US20110021892A1 (en) * 2009-07-23 2011-01-27 Nellcor Puritan Bennett Ireland Systems and methods for respiration monitoring
US20110021941A1 (en) * 2009-07-23 2011-01-27 Nellcor Puritan Bennett Ireland Systems and methods for respiration monitoring
US8346333B2 (en) 2009-07-30 2013-01-01 Nellcor Puritan Bennett Ireland Systems and methods for estimating values of a continuous wavelet transform
US8594759B2 (en) * 2009-07-30 2013-11-26 Nellcor Puritan Bennett Ireland Systems and methods for resolving the continuous wavelet transform of a signal
US8478376B2 (en) * 2009-07-30 2013-07-02 Nellcor Puritan Bennett Ireland Systems and methods for determining physiological information using selective transform data
US8755854B2 (en) 2009-07-31 2014-06-17 Nellcor Puritan Bennett Ireland Methods and apparatus for producing and using lightly filtered photoplethysmograph signals
US8628477B2 (en) * 2009-07-31 2014-01-14 Nellcor Puritan Bennett Ireland Systems and methods for non-invasive determination of blood pressure
US9220440B2 (en) 2009-09-21 2015-12-29 Nellcor Puritan Bennett Ireland Determining a characteristic respiration rate
WO2011037699A2 (en) 2009-09-24 2011-03-31 Nellcor Puritan Bennett Llc Determination of a physiological parameter
US8923945B2 (en) 2009-09-24 2014-12-30 Covidien Lp Determination of a physiological parameter
US8400149B2 (en) * 2009-09-25 2013-03-19 Nellcor Puritan Bennett Ireland Systems and methods for gating an imaging device
US9066660B2 (en) 2009-09-29 2015-06-30 Nellcor Puritan Bennett Ireland Systems and methods for high-pass filtering a photoplethysmograph signal
US20110077484A1 (en) * 2009-09-30 2011-03-31 Nellcor Puritan Bennett Ireland Systems And Methods For Identifying Non-Corrupted Signal Segments For Use In Determining Physiological Parameters
US8463347B2 (en) 2009-09-30 2013-06-11 Nellcor Puritan Bennett Ireland Systems and methods for normalizing a plethysmograph signal for improved feature analysis
US20110098933A1 (en) * 2009-10-26 2011-04-28 Nellcor Puritan Bennett Ireland Systems And Methods For Processing Oximetry Signals Using Least Median Squares Techniques
KR101081659B1 (ko) * 2010-01-29 2011-11-09 이병훈 병명이 표시되는 혈압기
US9451887B2 (en) 2010-03-31 2016-09-27 Nellcor Puritan Bennett Ireland Systems and methods for measuring electromechanical delay of the heart
US8898037B2 (en) 2010-04-28 2014-11-25 Nellcor Puritan Bennett Ireland Systems and methods for signal monitoring using Lissajous figures
US9050043B2 (en) 2010-05-04 2015-06-09 Nellcor Puritan Bennett Ireland Systems and methods for wavelet transform scale-dependent multiple-archetyping
WO2011141765A1 (en) * 2010-05-14 2011-11-17 Centre For Development Of Advanced Computing Diagnostic classifications of pulse signal waveform data
US8834378B2 (en) 2010-07-30 2014-09-16 Nellcor Puritan Bennett Ireland Systems and methods for determining respiratory effort
US9532721B2 (en) 2010-08-06 2017-01-03 The United States Of America As Represented By The Secretary Of The Army Patient care recommendation system
JP5733499B2 (ja) * 2010-10-29 2015-06-10 株式会社デルタツーリング 生体状態推定装置及びコンピュータプログラム
US8825428B2 (en) 2010-11-30 2014-09-02 Neilcor Puritan Bennett Ireland Methods and systems for recalibrating a blood pressure monitor with memory
US9357934B2 (en) 2010-12-01 2016-06-07 Nellcor Puritan Bennett Ireland Systems and methods for physiological event marking
US9259160B2 (en) 2010-12-01 2016-02-16 Nellcor Puritan Bennett Ireland Systems and methods for determining when to measure a physiological parameter
US9113830B2 (en) 2011-05-31 2015-08-25 Nellcor Puritan Bennett Ireland Systems and methods for detecting and monitoring arrhythmias using the PPG
US20120316458A1 (en) * 2011-06-11 2012-12-13 Aliphcom, Inc. Data-capable band for medical diagnosis, monitoring, and treatment
US9597022B2 (en) 2011-09-09 2017-03-21 Nellcor Puritan Bennett Ireland Venous oxygen saturation systems and methods
WO2013052896A1 (en) * 2011-10-07 2013-04-11 Edwards Lifesciences Corporation Detecting a vasoactive agent in the bloodstream
US9060695B2 (en) 2011-11-30 2015-06-23 Covidien Lp Systems and methods for determining differential pulse transit time from the phase difference of two analog plethysmographs
US20130158417A1 (en) * 2011-12-16 2013-06-20 General Electric Company Method, apparatus and computer program for automatic non-invasive blood pressure measurement
US20150065896A1 (en) * 2012-03-30 2015-03-05 Seiko Epson Corporation Pulsation detecting device, electronic apparatus, and program
JP6010979B2 (ja) * 2012-03-30 2016-10-19 セイコーエプソン株式会社 拍動検出装置、電子機器及びプログラム
US20140135632A1 (en) * 2012-11-15 2014-05-15 Pulsecor Limited Method and apparatus for determining cardiac medical parameters from supra-systolic signals obtained from an oscillometric blood pressure system
US20140288885A1 (en) * 2013-03-19 2014-09-25 Seiko Epson Corporation Signal processing device, pulse wave measuring apparatus, and signal processing method
US10226188B2 (en) 2013-08-23 2019-03-12 Covidien Lp Systems and methods for monitoring blood pressure
US10022068B2 (en) 2013-10-28 2018-07-17 Covidien Lp Systems and methods for detecting held breath events
US9955894B2 (en) 2014-01-28 2018-05-01 Covidien Lp Non-stationary feature relationship parameters for awareness monitoring
JP2017520286A (ja) 2014-05-09 2017-07-27 ボストン サイエンティフィック サイムド,インコーポレイテッドBoston Scientific Scimed,Inc. 心臓組織をマッピングする医療用デバイス
ES2747822T3 (es) 2014-09-25 2020-03-11 Aseptika Ltd Dispositivo médico
JP2017535318A (ja) 2014-10-03 2017-11-30 ボストン サイエンティフィック サイムド,インコーポレイテッドBoston Scientific Scimed,Inc. 心臓組織をマッピングするための医療システム
JP2017534344A (ja) * 2014-10-03 2017-11-24 ボストン サイエンティフィック サイムド,インコーポレイテッドBoston Scientific Scimed,Inc. 心臓組織をマッピングするための医療装置
EP3206586B1 (en) 2014-10-16 2021-03-24 Viewcare Technologies 1 ApS A method of detecting dicrotic notch
DE102014225483B3 (de) * 2014-12-10 2016-05-04 Gert Küchler Verfahren und Gerät zur Bestimmung mindestens eines physiologischen Parameters
EP3287069B1 (en) 2015-04-21 2019-11-13 Shinano Kenshi Co., Ltd. Biological information reading device
JP6524858B2 (ja) * 2015-08-24 2019-06-05 オムロンヘルスケア株式会社 脈波測定装置
JP6528602B2 (ja) * 2015-08-24 2019-06-12 オムロンヘルスケア株式会社 圧脈波センサ及び生体情報測定装置
US10226215B2 (en) 2015-08-28 2019-03-12 Covidien Lp Cable management feature for wearable medical monitor
WO2017179694A1 (ja) * 2016-04-15 2017-10-19 オムロン株式会社 生体情報分析装置、システム、プログラム、及び、生体情報分析方法
US10751004B2 (en) * 2016-07-08 2020-08-25 Edwards Lifesciences Corporation Predictive weighting of hypotension profiling parameters
JP6202510B1 (ja) * 2017-01-31 2017-09-27 株式会社Arblet 血圧情報測定システム、血圧情報測定方法、血圧情報測定プログラム、血圧情報測定装置、サーバ装置、演算方法及び演算プログラム
JP6747344B2 (ja) * 2017-03-14 2020-08-26 オムロンヘルスケア株式会社 血圧データ処理装置、血圧データ処理方法および血圧データ処理プログラム
JP7187824B2 (ja) * 2018-05-31 2022-12-13 セイコーエプソン株式会社 生体解析装置、生体解析方法およびプログラム
JP7069598B2 (ja) * 2017-08-16 2022-05-18 セイコーエプソン株式会社 生体解析装置、生体解析方法およびプログラム
CN107898443B (zh) * 2017-11-21 2020-11-24 深圳先进技术研究院 重搏波检测方法、装置及计算机存储介质
CN109833035B (zh) * 2017-11-28 2021-12-07 深圳市岩尚科技有限公司 脉搏波血压测量装置的分类预测数据处理方法
JP7285843B2 (ja) * 2017-12-22 2023-06-02 アシスタンス パブリック-ホピトー デ パリ 平均動脈圧を測定するためのシステム
US11660053B2 (en) * 2018-04-16 2023-05-30 Samsung Electronics Co., Ltd. Apparatus and method for monitoring bio-signal measuring condition, and apparatus and method for measuring bio-information
KR102655739B1 (ko) * 2018-04-16 2024-04-05 삼성전자주식회사 생체신호 측정 상태 모니터링 장치 및 방법과, 생체정보 측정 장치 및 방법
KR102147568B1 (ko) * 2018-06-01 2020-08-24 한국 한의학 연구원 맥파 신호의 평활도를 결정하는 장치 및 방법
CN110575198B (zh) * 2018-06-08 2022-07-01 佳能医疗系统株式会社 解析装置及解析方法
TWI731339B (zh) * 2018-06-11 2021-06-21 當代漢雲企業有限公司 脈診量測裝置及方法
CN115066203A (zh) 2020-01-13 2022-09-16 梅西莫股份有限公司 具有生理参数监测的可穿戴设备
CN111528821A (zh) * 2020-05-27 2020-08-14 福州数据技术研究院有限公司 一种脉搏波中重搏波特征点识别方法
CN113288075B (zh) * 2021-06-01 2022-07-19 山西和脉互联网智慧中医产业技术研究院有限公司 浮中沉脉位识别方法及系统
WO2023223087A1 (en) * 2022-05-20 2023-11-23 L&T Technology Services Limited Non-invasive blood pressure measurement
CN115381412B (zh) * 2022-09-23 2023-07-14 广东省新黄埔中医药联合创新研究院 脉搏波阵列信号包络图结合时域信号识别弦脉特征的方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59189830A (ja) 1983-04-11 1984-10-27 コーリン電子株式会社 脈波信号処理装置
US4958638A (en) * 1988-06-30 1990-09-25 Georgia Tech Research Corporation Non-contact vital signs monitor
US5027641A (en) * 1989-02-23 1991-07-02 Costello Jr Leo F Oscillometric non-invasive blood pressure simulator
JPH04250132A (ja) 1991-01-28 1992-09-07 Matsushita Electric Works Ltd 電子血圧計
JP2877270B2 (ja) 1992-01-20 1999-03-31 三菱重工業株式会社 自動塗装装置
JP3047712B2 (ja) * 1993-11-19 2000-06-05 セイコーエプソン株式会社 脈波診断装置
WO1994015526A1 (en) * 1993-01-07 1994-07-21 Seiko Epson Corporation Pulse wave analyzer, and diagnosis apparatus using the same
JP2979933B2 (ja) * 1993-08-03 1999-11-22 セイコーエプソン株式会社 脈波解析装置
JP3605216B2 (ja) 1995-02-20 2004-12-22 セイコーエプソン株式会社 脈拍計
JP3649464B2 (ja) 1995-03-27 2005-05-18 コーリンメディカルテクノロジー株式会社 心機能評価装置
EP0778001B1 (en) * 1995-05-12 2004-04-07 Seiko Epson Corporation Apparatus for diagnosing condition of living organism and control unit
JP3301294B2 (ja) 1995-09-13 2002-07-15 セイコーエプソン株式会社 健康状態管理装置
GB9600209D0 (en) * 1996-01-05 1996-03-06 Monitoring Tech Ltd Improved method and apparatus for the measurement of cardiac output
DE69723946T2 (de) * 1996-04-17 2004-07-15 Seiko Epson Corp. Arrhythmiedetektor

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100396233C (zh) * 2002-07-09 2008-06-25 健资国际私人有限公司 监测血压的方法和设备
CN102727191A (zh) * 2011-03-29 2012-10-17 精工爱普生株式会社 脉波计以及信号处理方法
CN102727191B (zh) * 2011-03-29 2015-08-05 精工爱普生株式会社 脉波计以及信号处理方法
CN102955891A (zh) * 2011-08-30 2013-03-06 华晶科技股份有限公司 具有资料检测功能的电子装置及其资料检测方法
CN104254275A (zh) * 2012-02-22 2014-12-31 阿克拉里斯医疗有限责任公司 生理信号检测装置和系统
CN104490373A (zh) * 2014-12-17 2015-04-08 辛勤 脉搏信号的判断方法、判断装置以及生理参数测量设备
CN107735018A (zh) * 2015-07-01 2018-02-23 浜松光子学株式会社 血压比计算装置、血压比计算方法、血压比计算程序及存储该程序的存储介质
CN107708533A (zh) * 2015-07-01 2018-02-16 浜松光子学株式会社 粘弹特性获取装置、粘弹特性获取方法、粘弹特性获取程序及存储该程序的存储介质
CN107708533B (zh) * 2015-07-01 2020-07-31 浜松光子学株式会社 粘弹特性获取装置、粘弹特性获取方法、粘弹特性获取程序及存储该程序的存储介质
US10772513B2 (en) 2015-07-01 2020-09-15 Hamamatsu Photonics K.K. Blood pressure ratio calculation device, blood pressure ratio calculation method, blood pressure ratio calculation program, and recording medium recording said program
CN107735018B (zh) * 2015-07-01 2020-12-08 浜松光子学株式会社 血压比计算装置、血压比计算方法、血压比计算程序及存储该程序的存储介质
US11154206B2 (en) 2015-07-01 2021-10-26 Hamamatsu Photonics K.K. Viscoelasticity characteristics acquisition device, viscoelasticity characteristics acquisition method, viscoelasticity characteristics acquisition program, and recording medium recording said program
CN108289624A (zh) * 2015-11-27 2018-07-17 浜松光子学株式会社 血压信息计算装置、血压信息计算方法、血压信息计算程序及存储该程序的存储介质
US11089965B2 (en) 2015-11-27 2021-08-17 Hamamatsu Photonics K.K. Blood pressure information calculating device, blood pressure information calculating method, blood pressure information calculating program, and recording medium for recording said program
CN106539570A (zh) * 2016-07-04 2017-03-29 悦享趋势科技(北京)有限责任公司 定位潮波的方法及装置
CN108403094A (zh) * 2018-03-23 2018-08-17 中国科学院微电子研究所 一种识别脉搏波波峰的方法

Also Published As

Publication number Publication date
JP3843462B2 (ja) 2006-11-08
DE69832782D1 (de) 2006-01-19
CN1326489C (zh) 2007-07-18
US6293915B1 (en) 2001-09-25
EP1477112A1 (en) 2004-11-17
HK1027496A1 (en) 2001-01-19
DE69835843D1 (de) 2006-10-19
DE69832782T2 (de) 2006-09-21
EP0968681A1 (en) 2000-01-05
EP0968681B1 (en) 2005-12-14
EP0968681A4 (en) 2000-08-16
WO1999026529A1 (fr) 1999-06-03
DE69835843T2 (de) 2007-09-20
TW425278B (en) 2001-03-11
EP1477112B1 (en) 2006-09-06

Similar Documents

Publication Publication Date Title
CN1251027A (zh) 脉波诊断装置、血压监视装置、脉波形状监视装置和药理作用监视装置
CN1154432C (zh) 用于测量生理状态的装置
CN1103580C (zh) 脉波分析装置
CN1155331C (zh) 锻炼测验用辅助装置
CN1255078C (zh) 体温测量装置
CN1182890C (zh) 放松指导装置和生物反馈指导装置
CN1268282C (zh) 脉搏检测装置
CN1227843C (zh) 生物体用偏光通信系统
CN1226959C (zh) 活体状态的诊断装置和控制装置
CN1242693A (zh) 脉波诊断装置、运动指标检测装置、运动强度检测装置、心输出量检测装置、每搏输出量检测装置、心功能诊断装置及其检测方法
CN1102036A (zh) 脉波分析装置和使用该装置的诊断装置
CN1155332C (zh) 心律失常检测设备
CN1292706C (zh) 脉搏计及其控制方法、以及手表型信息装置
CN1123319C (zh) 诊断心血管病状的装置
CN1193708C (zh) 带步数测定功能的内脏脂肪计
CN1185985C (zh) 锻炼量测量设备
CN1688247A (zh) 控制无创参数测量的方法和装置
CN1911163A (zh) 计算出血压平均值的电子血压计
CN1461200A (zh) 中枢血压波形推断装置以及末梢血压波形检测装置
CN1162138C (zh) 按摩机
CN1723842A (zh) 人体信息检测装置、使用该装置的健康管理系统
CN1852676A (zh) 生物体信息检测装置和血压计
CN1723848A (zh) 使用加速度传感器测量体育运动量的设备和方法
CN1226811A (zh) 触觉检测器,触觉报告器,信息输入器,触觉复现器,触觉传输系统,脉博诊断器,脉博诊断训练器以及脉博诊断信息传输器
CN1657005A (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

Granted publication date: 20070718

Termination date: 20151120

EXPY Termination of patent right or utility model