CN102469942A - 脉搏波解析装置以及记录介质 - Google Patents

脉搏波解析装置以及记录介质 Download PDF

Info

Publication number
CN102469942A
CN102469942A CN2010800342085A CN201080034208A CN102469942A CN 102469942 A CN102469942 A CN 102469942A CN 2010800342085 A CN2010800342085 A CN 2010800342085A CN 201080034208 A CN201080034208 A CN 201080034208A CN 102469942 A CN102469942 A CN 102469942A
Authority
CN
China
Prior art keywords
pulse wave
arteries
veins
degree
index
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
CN2010800342085A
Other languages
English (en)
Other versions
CN102469942B (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.)
Fukuda Denshi Co Ltd
Original Assignee
Omron Healthcare Co Ltd
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 Omron Healthcare Co Ltd filed Critical Omron Healthcare Co Ltd
Publication of CN102469942A publication Critical patent/CN102469942A/zh
Application granted granted Critical
Publication of CN102469942B publication Critical patent/CN102469942B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021—Measuring pressure in heart or blood vessels
    • A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02125—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7242—Details of waveform analysis using integration
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/30—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to physical therapies or activities, e.g. physiotherapy, acupressure or exercising
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/053—Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0535—Impedance plethysmography

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Cardiology (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Hematology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Data Mining & Analysis (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Databases & Information Systems (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

脉搏波解析装置(100)存储多个脉拍的脉搏波波形，并通过解析该多个脉拍的脉搏波波形来计算脉搏波解析指标。在脉搏波解析指标的计算中，对用于构成多个脉拍的脉搏波波形的每个脉拍的脉搏波形状进行累积，并从脉搏波解析指标的计算对象中，排除所累积的脉搏波形状与每个脉拍的脉搏波形状的近似度低的脉拍。

Description

脉搏波解析装置以及记录介质
技术领域
本发明涉及一种脉搏波解析装置，特别涉及一种能够通过解析多个脉拍的脉搏波波形来计算出预定的脉搏波解析指标的脉搏波解析装置以及脉搏波解析用程序。
背景技术
脉搏波解析用于测定脉搏波传播速度等脉搏波解析指标。脉搏波传播速度作为评价动脉硬度的非侵入性指标，应用于治疗现场。
作为用于精确测量脉搏波解析指标的技术，有以下一些技术。
在日本特开2006-247221号公报(专利文献1)中记载了，通过自相关函数波形，来判断脉搏波中是否含有干扰波。
在日本特开2001-128946号公报(专利文献2)中记载了，为了测定正确的脉搏波传播速度信息，检测触点并根据其时间差来计算脉搏波速度。
在日本特开平10-328150号(专利文献3)中记载了，为了高精度的测定脉搏波传播速度，用心跳同步波的最大倾斜线与基线来计算脉搏波传播速度。
在日本特开2008-168073号(专利文献4)中记载了，为了提高动脉硬化检查的可靠性及效率性，对取得的脉搏波的特征点进行检测，明确显示该特征点并实时地在画面上显示脉搏波波形。
现有技术文献
专利文献1：日本特开2006-247221号
专利文献2：日本特开2001-128946号
专利文献3：日本特开平10-328150号
专利文献4：日本特开2008-168073号
发明内容
技术问题
脉搏波传播速度是脉搏波解析指标的一种，例如，用如下方法求出。臂踝脉搏波传播速度(baPWV：brachial-ankle Pulse Wave Velocity)作为脉搏波传播速度的一种形态，记录几个脉拍至十几个脉拍的使缠在上臂和脚踝的袖带保持恒定压而选取的PVR(Pulse Volume Recording：脉搏波量记录)波形。然后，针对上臂以及脚踝各自的PVR波形，通过检测每一脉拍脉搏波的上升位置，来计算脉搏波传播速度。
在这种方法中，因为全部脉拍的脉搏波被用于计算脉搏波解析指标，所以在选取PVR波形的过程中引起心律不齐或身体移动等、脉搏波紊乱而使指标的测量精度降低。另外，其结果，还存在将错误的测定值(低精度的脉搏波解析指标)用于诊断的忧虑。
另外，上述各专利文献中的提案不足以高精度地计算出脉搏波解析指标。
本发明是为了解决如上所述的问题而作出的，本发明的目的是提供一种能够高精度地计算脉搏波解析指标的脉搏波解析装置以及记录脉搏波解析用程序的记录介质。
解决问题的手段
本发明的一个方面的脉搏波解析装置具有：存储部，其用于存储多个脉拍的脉搏波波形；解析处理部，其进行特定处理，该特定处理是指，通过对多个脉拍的脉搏波波形进行解析，来计算脉搏波解析指标的处理；输出部，其用于将计算出的脉搏波解析指标作为解析结果来输出。解析处理部对用于构成多个脉拍的脉搏波波形的每个脉拍的脉搏波形状进行累积，并从计算对象中排除所累积的脉搏波形状与每个脉拍的脉搏波形状之间的近似度低的脉拍，来计算脉搏波解析指标。
优选的是，解析处理部还对用于计算脉搏波解析指标的脉搏波形状的近似度进行累积，由此计算搏动的稳定度；输出部还将稳定度作为表示脉搏波解析指标的可靠性的指标来输出。
优选的是，存储部对应每个肢体部位来存储多个脉拍的脉搏波波形；解析处理部对应每个肢体部位来累积每个脉拍的脉搏波形状，并计算近似度、脉搏波解析指标以及稳定度；输出部将稳定度高的脉搏波解析指标作为解析结果来输出。
优选的是，存储部针对左右的肢体部位存储多个脉拍的脉搏波波形；解析处理部对应肢体部位计算近似度，并利用近似度高的肢体部位的脉搏波形，来计算脉搏波解析指标。
优选的是，解析处理部仅限于在每个脉拍的脉搏波形状中对脉搏波解析指标的计算产生影响的范围内，计算近似度。
优选的是，脉搏波解析指标表示动脉硬化的程度和/或血管的狭窄程度。
优选的是，脉搏波解析指标包含作为表示动脉硬化的程度的指标的脉搏波传播速度。
进而优选的是，脉搏波解析装置还具有脉搏波检测部，该脉搏波检测部用于对肢体部位的脉搏波进行检测；解析处理部基于来自脉搏波检测部的检测信号，来测定多个脉拍的脉搏波波形。
本发明的另一个方面的记录介质，记录有脉搏波解析用程序，该脉搏波解析用程序使计算机发挥用于解析脉搏波的装置的功能，脉搏波解析用程序使计算机执行以下步骤：对用于构成存储于存储部中的多个脉拍的脉搏波波形的每个脉拍的脉搏波形状进行累积的步骤；从计算对象中排除所累积的脉搏波形状与每个脉拍的脉搏波形状之间的近似度低的脉拍，来计算脉搏波解析指标的步骤；将计算出的脉搏波解析指标作为解析结果来输出的步骤。
发明的效果
根据本发明，只利用很少受到身体移动等影响的稳定的脉拍便能够计算脉搏波解析指标。其结果，能够将高精度的脉搏波解析指标作为解析结果来输出。
附图说明
图1是本发明实施方式的脉搏波解析装置的简要结构图。
图2是表示本发明实施方式的脉搏波解析装置的功能结构的功能模块图。
图3是表示每个肢体部位的脉搏波测定结果的一例图。
图4是用于说明本发明实施方式中的、累积的脉搏波形状和每个脉拍的近似度的计算方法的图。
图5是用于说明脉搏波传播距离的计算方法的图。
图6是表示本发明实施方式中的脉搏波解析处理的流程图。
图7是表示图6中的步骤S108中的排除处理结果的一例的图。
图8是表示本发明实施方式中的解析结果信息的输出例的图。
图9是表示本发明实施方式中的解析结果信息的再一输出例的图。
图10是表示本发明实施方式中的解析结果信息的另一输出例的图。
图11A是表示测定ID1的多个脉拍的脉搏波波形的图。
图11B是以上升位置为起点重叠图11A的每个脉拍的脉搏波形状而示出的图。
图12A是表示测定ID2的多个脉拍的脉搏波波形的图。
图12B是以上升位置为起点重叠图12A的每个脉拍的脉搏波形状而示出的图。
国13A是表示测定ID3的多个脉拍的脉搏波波形的图。
图13B是以上升位置为起点重叠图13A的每个脉拍的脉搏波形状而示出的图。
图14A是表示测定ID4的多个脉拍的脉搏波波形的图。
图14B是以上升位置为起点重叠图14A的每个脉拍的脉搏波形状而示出的图。
图15A是表示测定ID5的多个脉拍的脉搏波波形的图。
图15B是以上升位置为起点重叠图15A的每个脉拍的脉搏波形状而示出的图。
图16A是表示测定ID6的多个脉拍的脉搏波波形的图。
图16B是以上升位置为起点重叠图16A的每个脉拍的脉搏波形状而示出的图。
图17A是表示测定ID7的多个脉拍的脉搏波波形的图。
園17B是以上升位置为起点重叠图17A的每个脉拍的脉搏波形状而示出的图。
图18A是表示测定ID8的多个脉拍的脉搏波波形的图。
图18B是以上升位置为起点重叠图18A的每个脉拍的脉搏波形状而示出的图。
图19A是表示测定ID9的多个脉拍的脉搏波波形的图。
图19B是以上升位置为起点重叠图19A的每个脉拍的脉搏波形状而示出的图。
图20A是表示测定ID10的多个脉拍的脉搏波波形的图。
图20B是以上升位置为起点重叠图20A的每个脉拍的脉搏波形状而示出的图。
图21A是表示测定ID11的多个脉拍的脉搏波波形的图。
图21B是以上升位置为起点重叠图21A的每个脉拍的脉搏波形状而示出的图。
图22A是表示测定ID12的多个脉拍的脉搏波波形的图。
图22B是以上升位置为起点重叠图22A的每个脉拍的脉搏波形状而示出的图。
图23A是表示测定ID13的多个脉拍的脉搏波波形的图。
图23B是以上升位置为起点重叠图23A的每个脉拍的脉搏波形状而示出的图。
图24A是表示测定ID14的多个脉拍的脉搏波波形的图。
图24B是以上升位置为起点重叠图24A的每个脉拍的脉搏波形状而示出的图。
图25A是表示测定ID15的多个脉拍的脉搏波波形的图。
图25B是以上升位置为起点重叠图25A的每个脉拍的脉搏波形状而示出的图。
图26A是表示测定ID16的多个脉拍的脉搏波波形的图。
图26B是以上升位置为起点重叠图26A的每个脉拍的脉搏波形状而示出的图。
图27A是表示测定ID17的多个脉拍的脉搏波波形的图。
图27B是以上升位置为起点重叠图27A的每个脉拍的脉搏波形状而示出的图。
图28是表示将图11A至图27B的脉搏波波形作为对象的情况下的，由装置计算的近似度的位次和由观测员观测的近似程度的位次之间的关系的图。
图29是表示图28中所示的、基于装置的近似的位次和基于观测员的近似的位次之间的相关关系的图。
图30是用于说明本发明实施方式的变形例2的脉搏波解析指标的计算方法的图。
具体实施方式
以下，参照附图对本发明的实施方式进行详细说明。其中，图中的相同的附图标记表示相同或者相应的部分，不再重复进行说明。
(简要结构)
图1是本发明实施方式的脉搏波解析装置的简要结构图。
参照图1，脉搏波解析装置100包含：信息处理单元1；四个检测单元20ar、20al、20br、20bl；四个袖带24ar、24al、24br、24bl。
袖带24ar、24al、24br、24bl，分别在被测定者200的肢体部位佩戴。具体地说，分别在右上臂(右上肢)、左上臂(左上肢)、右脚踝(右下肢)以及左脚踝(左上肢)佩戴。其中，“肢体部位”表示包含四肢的部位，也可以是手腕及手指部等。袖带24ar、24al、24br、24bl只要没有特别区分，就统称为“袖带24”。
检测单元20ar、20al、20br、20bl分别包含必要的硬件，以检测被测定者200的肢体部位的脉搏波。因为检测单元20ar、20al、20br、20bl的结构全部一样，所以只要没有必特别区分，就统称为“检测单元20”。
信息处理单元1包含：控制部2、输出部4、操作部6、存储装置8。
控制部2是进行控制脉搏波解析装置100整体的装置，代表性的是由包括CPU(CentraI Processing Unit：中央处理单元)10、ROM(Read OnlyMemory：只读存储器)12、RAM(Random Access Memory：随机存取存储器)14的计算机构成。
CPU10相当于演算处理部，读取预先存储在ROM12的程序，将RAM14作为工作存储器进行使用的同时执行该程序。
另外，在控制部21连接有输出部4、操作部6以及存储装置8。输出部4输出所测定的脉搏波及脉搏波解析结果等。输出部4可以是由LED(LightEmitting Diode：发光二极管)或者LCD(Liquid Crystal Display：液晶显示器)等构成的显示设备，也可以是打印机(驱动程序)。
操作部6接收来自用户的指示。存储装置8保存各种数据及程序。控制部2的CPU10，对记录在存储装置8的数据及程序进行读出及写入。存储装置8例如可以由一个硬盘、非易失性存储器(例如，闪存)或者可拆卸的外部记录介质等构成。
在这里，具体说明各个检测单元20的结构。
检测单元20br通过对在被测定者200的右上臂佩戴的袖带24br的内部压力(以下，称其为“袖带压”)进行调整以及检测，来检测出右上臂的脉搏波。袖带24br内置有有未图示的流体袋(例如空气袋)。
检测单元20br包含压力传感器28br、调压阀26br、压力泵25br、A/D(anaIog to digital，模拟/数字)转换部29br、配管27br。袖带24br、压力传感器28br和调压阀26br通过配管22br相连接。
压力传感器28br是对经由配管22br传递的压力变化进行检测的检测部位。作为一例，包含在由单晶硅等构成的半导体芯片以预定间隔排列的多个传感器元件。由压力传感器28br检测出的压力变化信号通过A/D转换部29br转换为数字信号，作为脉搏波信号pbr(t)输入至控制部2。
调压阀26br插入在压力泵25br和袖带24br之间，在测定时，使用于袖带24br的加压的压力保护在预定的范围内。压力泵25br按照来自控制部2的检测指令进行动作，并为了对袖带24br进行加压而向袖带24br内的流体袋(未图示)供给空气。
通过该加压，袖带24br对测定部位施加压力，与右上臂的脉搏波对应的压力变化分别经由各个配管22br传递至检测单元20br。检测单元20br，通过检测该传递的压力变化，检测右上臂的脉搏波。
检测单元20bl也同样包含压力传感器28bl、调压阀26bl、压力泵25bl、A/D转换部29bl、配管27bI。袖带24bl、压力传感器28bl和调压阀26bI通过配管22bI相连接。
另外，检测单元20ar包含压力传感器28ar、调压阀26ar、压力泵25ar、A/D转换部29ar、配管27ar。袖带24ar、压力传感器28ar和调压阀26ar通过配管22ar相连接。
检测单元20al也同样包含压力传感器28aI、调压阀26al、压力泵25al、A/D转换部29aI、配管27aI。袖带24al、压力传感器28al和调压阀26al通过配管22aI相连接。
因为检测单元20bl、20ar、20al内的各个部分的功能与检测单元20br是一样的，所以不反复进行详细说明。另外，关于检测单元20内的各个部分，只要没有特别区分，就省略“ar”、“br”等标记而说明。
此外，在本实施方式中，虽然对用压力传感器28检测脉搏波的结构进行了说明，但也可以是用动脉容量传感器(未图示)检测脉搏波的结构。此时，动脉容量传感器例如可包含：发光元件，对动脉照射光；光接收元件，接收由发光元件照射的光的透射光或者反射光。或者，可包含多个电极，使极少的恒定电流在被测定者200的测定部位流动，并且检测电压变化，该电压变化是随着与脉搏波的传播对应地发生的阻抗(生物阻抗)变化而发生的。
(功能结构)
图2是表示本发明实施方式的脉搏波解析装置的功能结构的功能模块图。
参照图2，本实施方式的脉搏波解析装置100包含调整部30、脉搏波测定部102、解析处理部104、血压测定部106、血压指标计算部108以及输出部4作为其功能结构。此外，血压测定部106以及血压指标计算部108可以不包含在脉搏波解析装置100的功能结构内。
调整部30是调整袖带24内的压力的功能部。例如，如图1所示，调整部30的功能通过压力泵25以及调压阀26达成。
脉搏波测定部102与调整部30以及A/D转换部29相连接，进行用于测定各个肢体部位的脉搏波(pVR)的处理。脉搏波测定部102通过向调整部30发送指令信号来调整袖带24的内部压力，并接收响应该指令信号而检测出的袖带压信号Par(t)、PaI(t)、Pbr(t)、PbI(t)。然后，按时间序列记录所接收的袖带压信号Par(t)、Pal(t)、Pbr(t)、Pbl(t)，对应每个肢体部位，得到多个脉拍的脉搏波波形。例如，在预定时间内(例如10秒左右)进行脉搏波的测定。
由脉搏波测定部102测定的脉搏波测定结果输出至输出部4。
图3是表示每个肢体部位的脉搏波测定结果的一例的图。在图3中，沿相同的时间轴表示各个肢体部位的脉搏波波形。如图3所示，可以用虚线等表示每个脉拍的脉搏波的上升位置。
解析处理部104通过解析由脉搏波测定部102测定的每个肢体部位的脉搏波，来计算预定的脉搏波解析指标(以下简称“解析指标”)作为被测定者200(图1)的脉搏波的特征量。本实施方式中的“解析指标”表示，与动脉硬化和/或血管的狭窄具有相关性的指标。就是说，“解析指标”表示动脉硬化程度和/或血管狭窄程度。
作为表示动脉硬化程度的解析指标，例如可例举脉搏波传播速度、脉搏波传播时间(PTT：Pulse Transit Time)、A1(Augmentation Index：增强指数)、以及TR(Traveling time to Reflectedwave：反射波传播时间)等。此外，脉搏波传播速度不限定于通过上臂脉搏波以及脚踝脉搏波计算出的脉搏波传播速度(baPWV)，也可通过其他两个测定部位的脉搏波计算出，或者仅通过一个测定部位(肢体部位)的脉搏波计算。
作为表示血管狭窄程度的解析指标，例如可例举脚踝脉搏波的上升特征值以及脉搏波的锐度等。例如，踝脉搏波的上升特征值作为UT(UT：upstrokeTime：上升时间)而被计算出。UT作为从上升点至顶点的脚踝脉搏波上升的期间而被计算出。脉搏波的锐度作为％MAP(正常化脉搏波面积)而被计算出。％MAP例如作为使脉搏波面积均等时的最低血压的高度M与脉搏波的顶点高度H即脉搏压的比例(＝M/H×100)而被计算出。
在本实施方式中，说明的虽然是计算出baPWV来作为解析指标，但也可计算出如上所述的其他特征量来作为解析指标。
解析处理部104进行从多个脉拍的脉搏波波形识别每个脉拍的脉搏波形状(脉搏波波形的形状)的处理。具体地说，进行脉搏波的分离处理，从每个脉拍提取脉搏波波形。由此，识别每个脉拍的脉搏波形状。脉搏波的分离处理可以通过基于特定频率的滤波处理或者微分处理等公知的手法实现。
解析处理部104，累积所识别的每个脉拍的脉搏波形状，对应每个脉拍计算出与所累积的脉搏波形状(以下称为“累积形状”)的近似度。在本实施方式中，“累积每个脉拍的脉搏波形状”是指，虽然所表示的是对每个脉拍的脉搏波形状进行平均化，但也可进行与平均化相等的处理。
在本实施方式中，“近似度”是指表示两个波形的近似程度的值，更具体的是，用数值表示两个波形相匹配的程度。例如，近似度通过下面公式(1)求出。
公式1
Figure BPA00001499128100101
图4是用于说明本发明实施方式中的、累积形状和每个脉拍的近似度的计算方法的图。
参照图4，近似度作为因累积形状Wa和实测的第i脉拍的脉搏波形状Wi的偏差而产生的面积的倒数被计算出。就是说，近似度能够作为以脉搏的上升作为起点时的、对应每个抽样时间的振幅值Pa、Pi的差的总和的倒数而被求出。
另外，也可以作为对应每个抽样时间的振幅值Pa、Pi的差的定积分值的倒数而被求出。
另外，对振幅值Pa、Pi的差赋予加权，例如，如下面公式(2)所示，也可由振幅值Pa、Pi的差的2次方和的倒数求出近似度。
公式2
近似度的计算公式是基于预先进行的试验结果确定的。有关近似度的计算公式的设定方法(原理)将在后面说明。
此外，在本实施方式中，虽然以脉搏的上升为起点而求出累积脉搏波和1脉拍的脉搏波整体的近似度，但也可以将用于计算近似度的区间限定在对解析指标造成强烈影响的范围内。例如，可以限定在从脉搏波波形的上升点至顶点的范围，也可以限定在脉搏波波形的前1/2部分。就是说，可以在1脉拍的脉搏波形状之中，将不影响计算解析指标的后半部分不用于计算近似度。
另外，在计算近似度时，虽然采用了以脉搏的上升为起点的振幅值的差，但起点的位置(使两个波形匹配的位置)可不限定于脉搏的上升处。例如，以脉搏波波形的顶点为起点等，可将规定的基准位置确定为起点。或者，不用规定的基准位置，而是，例如，对应每脉拍地将与累积波形的近似度最大的位置确定为起点。
另外，在本实施方式中，虽然是计算近似度，但也可计算与累积形状的“偏差度”。偏差度可作为公式(1)、公式(2)的近似度的倒数而被计算出。
解析处理部104指定被判定为与累积形状的近似度低(＝偏差度高)的脉搏波形状的脉拍，从解析指标的计算对象中排除所指定的脉拍。这样，通过排除与累积形状的近似度低的脉搏波形状，适当的排除心律不齐及身体移动导致突发性紊乱的可能性高的、不稳定的脉搏波。
以往，若测定如图3所示的多个脉拍的脉搏波波形，由医生肉眼判断波形是否紊乱。根据本实施方式，计算与累积形状的近似度，从而能够通过演算来辨别波形是否紊乱。
解析处理部104通过仅解析未排除的脉拍的脉搏波波形，即稳定的脉搏波波形，来计算baPWV。在本实施方式中，以左右两边的脚踝为测定部位，计算两种baPWV，例如右上臂-左脚踝的脉搏波传播速度(以下称为baPWV-RL)以及右上臂-右脚踝的脉搏波传播速度(以下称为baPWV-RR)。如此，计算出两种baPWV，是因为它们的差异也用于診断左下腿部以及右下腿部的动脉狭窄。
此外，在本实施方式中，上臂的测定部位是右边，因为这些已被确定为默认，但也可将左边定为基准。另外，例如，在右上臂的血压相比左上臂的血压低预定值(例如16-20mmHg)以下的情况下，可用左上臂代替右上臂作为测定部位。或者，在经由操作部6输入表示将左边作为测定部位指示的情况下，也可用左上臂代替右上臂作为测定部位。
在计算baPWV-RL以及baPWV-RR时，仅解析未排除的脉拍的脉搏波波形。解析处理部104针对作为对象的各个脉拍，计算上臂脉搏波以及脚踝脉搏波的上升位置的时间差(图5的时间Tr、TI)，通过将计算出的时间的平均值除以血管长度(＝脉搏波传播距离)，来计算各个baPW。
图5中，时间Tr、TI是表示与右上臂的脉搏波的上升点P0r的时间差。在将左上臂作为测定部位的情况下，用与左上臂的脉搏波的上升点POI的时间差，来计算左上臂-左脚踝的脉搏波传播速度(以下称为rbaPWV-LL)，以及左上臂-右脚踝的脉搏波传播速度(以下称为rbaPWV-LR)。图5中的时间Tr、TI表示脉搏波传播时间(PTT)。
此外，计算各个baPWV所需的血管长度通过将被测定者的身高代入预定的换算公式来计算出。例如，被测定者的身高通过操作部6输入。
解析处理部104进一步地通过对与用于计算解析指标的脉搏波形状相左关的近似度进行累积，来计算用于计算解析指标的脉搏波整体的搏动的稳定度。在本实施方式中，针对baPWV-RL以及baPWVRR，分别计算成为计算对象的脉搏波整体的搏动的稳定度。搏动的稳定度并不是所测定的全部脉拍的脉搏波波形，而是从用于计算baPWV的脉搏波波形本身导出的，因此与baPWV的可靠性有直接关联。因此，可以说由解析处理部104计算出的各个搏动的稳定度表示对应的baPWV的可靠性(可靠程度)。
此外，在不计算近似度而计算偏差度的情况下，用于计算解析指标的脉搏波形状的偏差度的累积值作为用于计算解析指标的脉搏波全体搏动的紊乱度而被计算出。
由解析处理部104计算出的各个baPWV以及各个稳定度输出至输出部4。输出部4输出baPWV-RL以及baPWVRR，并与各个baPWV关联，输出(显示或者打印)表示该值的可靠性的指标。作为表示各个baPWV的可靠性的指标，可输出作为搏动的稳定度而被计算出的值本身，也可以将计算出的值转换为电平值、标记或者符号等进行输出。
与脉搏波测定部102一样，血压测定部106与调整部30以及A/D转换部29相连接，进行测定各个肢体部位血压的处理。血压测定部106通过向调整部30发送指令信号来调整袖带24的内部压力，并且接收响应该指令信号而检测的袖带压信号par(t)、Pal(t)、Pbr(t)、Pbl(t)。血压测定部106，例如按照已知的示波法，测定被测定者200的每个肢体部位的最高血压以及最低血压。更具体地说，对应每个肢体部位，快速施加袖带压至预定的压力值，在缓慢地减压时，通过对要检测的时间序列的袖带压信号采用预定的算法，来计算最高血压以及最低血压。血压测定部106还可以测定脉搏脉拍数、平均血压以及脉搏压。
血压指标计算部108基于由血压测定部106测定的血压值，来计算预定的血压指标。本实施方式中的“血压指标”表示与血管的堵塞程度(动脉狭窄程度)相关的指标，具体的说，例如可举出ABI(Ankle Brachial Index：踝臂指数)。在本实施方式中，将左右两边的脚踝和上臂作为测定部位来计算ABI。例如，将右上臂的最高血压和右脚踝的最高血压之比以及右上臂的最高血压和左脚踝的最高血压之比分别作为每个ABI-RR以及ABI-RL来进行计算。此外，在计算各个ABI的过程中，上臂的测定部位为右边的情况是右上臂的最高血压高于左上臂的情况。在左上臂的最高血压高于右上臂的情况下，可用左上臂作为测定部位。另外，可以将右上臂和左上臂的最高血压的平均作为计算ABI时的上臂血压。另外，在本实施方式中，虽然将ABI作为血压指标来进行计算，但也可用其他的血压特征量。
血压测定部106的测定结果以及由血压指标计算部108计算出的ABI-RR以及ABI-RL输出至输出部4。输出部4除了输出表示baPWV-RL、baPWV-RR以及这些可靠性的指标之外，还输出每个肢体部位的血压值、ABI-RR以及ABI-RL。由此，医生等医疗从业人员能够更准确地的诊断是否存在动脉硬化疑虑。
以上说明的脉搏波测定部102、解析处理部104、血压测定部106以及血压指标计算部108的动作可通过执行存储在ROM12中的软件实现，或就其中至少一个部分而言，还可以通过硬件实现。另外，解析处理部104执行的处理中的一部分也可通过硬件实现。
(动作)
接着，有关本实施方式的脉搏波解析装置100的动作进行说明。在说明动作时，对本实施方式中最有特点的部分即由解析处理部104执行的处理进行说明。
图6是表示本发明实施方式中的脉搏波解析处理的流程图。图6的流程图表示的处理，作为程序预先存储在ROM12内，通过CPU10读出此程序并执行，来实现脉搏波解析处理的功能。
此外，在开始脉搏波解析处理的时间点，由脉搏波测定部102测定的每个肢体部位的脉搏波波形，存储在RAM14或者存储装置8上。就是说，本实施方式的脉搏波解析处理，不限于脉搏波测定之后的处理，也可对存储在存储装置8内的、过去测定的脉搏波波形进行处理。
参照图6，解析处理部104对作为测定部位的每个肢体部位，对存储的多个脉拍的脉搏波波形进行分离处理(步骤S102)。由此，多个脉拍的脉搏波波形以一脉拍单位分离，识别每个脉拍的脉搏波形状。
接着，解析处理部104对在每个肢体部位识别的全部脉拍的脉搏波形状进行平均化(步骤S104)，计算所平均化的脉搏波形状(累积形状)和各个脉搏波形状的近似度(步骤S106)。在计算近似度时，例如，使用上述公式(1)。对应每个肢体部位计算出的各个脉搏波的近似度暂时记录在RAM14中。
解析处理部104进行排除近似度低的脉拍的处理(步骤S108)。具体地说，首先，对应每个肢体部位，指定近似度不满足预定基准的(就是说，近似度为预先确定的阈值以下的)脉搏波形状，并从baPWV的计算对象中排除所指定的脉搏波形状。此外，用于指定脉搏波波形的阈值，不限定于预先确定的值。例如，可以根据全部脉拍的平均近似度确定阈值，来指定不满足该基准(确定的阈值)的脉搏波形状作为排除对象。排除处理的结果，即是否用于计算baPWV的信息被暂时记录在RAM14中。
图7是表示图6中的步骤S108中的排除处理结果的一例的图。
参照图7，与各个脉拍(i＝1、2、3、……、n)对应地，记录表示右上臂、左上臂、右脚踝以及左脚踝各自的脉搏波形状是否已从计算对象中排除的信息。在本实施方式中，例如，在与判定为应该排除的脉搏波波形的脉拍编号对应的栏中记录“1”，其他为“0”。此外，各个肢体部位的排除处理结果的保存方法不限定于图7的示例。
若排除处理结束，则反映该处理结果而计算baPWV(脉搏波传播速度)(步骤S110)。在本实施方式中，步骤S108的排除处理的结果，排除被判定为近似度低的脉拍，而计算baPWVRL以及baPWVRR。
这里，假设记录有如图7所示的结果，来说明各个baPWV的计算方法。在计算baPWV-RL的情况下，利用右上臂的脉搏波和左脚踝的脉搏波。参照图7时，右上臂的第3脉拍的脉搏波和左脚踝的第6脉拍的脉搏波作为排除对象进行记录。因此，排除第3脉拍以及第6脉拍的脉搏波并计算baPWV-RL。更具体地说，针对第3脉拍以及第6脉拍之外的各个脉拍，计算右上臂以及左脚踝的脉搏波的脉搏波传播时间，并根据计算出的脉搏波传播时间的平均值和血管长度的估计值，来计算出baPWV-RL。
在计算baPWV-RR的情况下，利用右上臂的脉搏波与右脚踝的脉搏波。参照图7时，右上臂的第3脉拍的脉搏波与右脚踝的第5脉拍、第6脉拍的脉搏波作为排除对象进行记录。因此，通过排除第3脉拍、第5脉拍以及第6脉拍的脉搏波来计算baPWV-RR。更具体的说，针对3脉拍、第5脉拍以及第6脉拍之外的各个脉拍，计算右上臂以及右脚踝的脉搏波的脉搏波传播时间，并根据计算出的脉搏波传播时间的平均值和血管长度的估计值，来计算baPWV-RR。
这样，根据本实施方式，由于大大紊乱的脉搏波波形不用于计算各个baPWV，因而能够高精度地计算解析指标。另外，通过累积实测的波形，求出成为计算近似度的基准的脉搏波形状。因此，根据测定时的被测定者的病情及症状或者测定条件(吃药之后等)，能够适当地指定应排除的脉搏波形状。
接着，解析处理部104计算baPWV的可靠性(步骤S112)。在本实施方式中，分别计算用于计算baPWV-RL以及baPWV-RR的脉搏波整体的搏动的稳定度。具体地说，用于计算baPWV-RL的脉搏波全体搏动的稳定度，即baPWV-RL的可靠性通过对用于计算baPWV-RL的脉拍的近似度进行累积(例如平均)来计算。在计算baPWV-RL时，在图7中的例子中，由于排除了第3脉拍以及第6脉拍的脉搏波，因而利用除第3脉拍以及第6脉拍之外的、有关右上臂以及左脚踝的脉拍的近似度的平均值来表示baPWV-RL的可靠性。
更具体地说，求出右上臂的脉搏波形状中的、除第3脉拍以及第6脉拍之外的脉拍近似度的平均值以及左脚踝的脉搏波形状中的、除第3脉拍以及第6脉拍之外的脉拍近似度的平均值。将这些平均值进行一步平均化的值作为用于计算的脉搏波整体的稳定度而被计算出。
或者，计算右上臂的脉搏波形状中的、除第3脉拍以及第6脉拍之外的脉拍以及左脚踝的脉搏波形状中的、除第3脉拍以及第6脉拍之外的所有脉拍的近似度的平均值来作为用于计算的脉搏波整体的稳定度。
baPWV-RR的可靠性也用同样的方法计算。在计算baPWV-RR过程中，在图7的例子中，因为排除了第3脉拍、第5脉拍以及第6脉拍的脉搏波，通过累积(例如平均)除第3脉拍、第5脉拍以及第6脉拍之外的各个脉拍累积形状的近似度，计算baPWV-RR的可靠性。
这样，以对累积形状的近似度进行个别评价后，通过累积用于计算指标的全部脉拍的近似度，求出各个baPWV的可靠性。因此，能够使1脉拍对全部脉拍造成的影响度与以往开始进行过来的baPWV的计算方法(对应每个脉拍计算脉搏波传播时间，并将血管长度的估计值除以平均估计值)相同。
若上述的解析处理结束，则向输出部4输出各个解析结果(步骤S114)。在本实施方式中，具有打印机(驱动程序)功能的输出部4将得到的解析结果打印在纸质介质上。在纸质介质上，例如，打印有如图8所示的解析结果信息。
图8是表示本发明实施方式中的解析结果信息的输出例的图。
参照图8，作为脉搏波解析的结果，在纸质介质上打印出baPWV-RR值91、表示baPWV-RR的可靠性的指标92、baPWV-RL值93以及表示baPWV-RL的可靠性的指标94。例如，表示可靠性的指标是，以可靠性高(稳定度高)低的顺序，用“+++”、“++”、“+”、“±”、“-”的5个阶段的符号表示。假设与稳定度的数值范围对应地，应显示的这些符号预先存储在ROM12中。
这样，通过在各个baPWV值91、93的正下方配置它们的值的可靠性的指标92、94，来将这些指标关联起来后进行打印。其结果，医生等医疗从事人员不仅考虑作为解析指标而输出的baPWV值，还可以考虑该指标的可靠性高度，来准确地进行诊断。
如上所述，本实施方式的脉搏波解析装置100也进行各个肢体部位的血压测定，并计算ABI-RR、ABI-RL。因此，在图8中，进一步输出右上臂的血压81、右脚踝的血压82、左上臂的血压83、左脚踝的血压84、ABI-RR值85以及ABI-RL值86。
图8所示的baPWV以及血压值的单位分别是“cm/s”以及“mmHg″。
作为解析结果信息，也可以进一步输出上述的UT及％MAP。另外，输出如图9所示的图标。
图9是表示本发明实施方式中的解析结果信息的再一输出例的图。
在图9中，是以baPWV为纵轴、ABI为横轴的图表。在该图表中，(例如，用颜色区分)示有预先能够识别的、根据baPWV和ABI来确定的动脉硬化度的程度。在该图表中，根据预先确定的标记、文字、符号等表示右下腿部以及左下腿部各自的动脉硬化程度。
在图9中，利用在由解析处理部104计算出的baPWV-RR和由血压指标计算部108计算出的ABI-RR的交叉点标绘的涂黑的三角标记位置，来表示右下腿部的动脉硬化度的程度。利用在由解析处理部104计算出的baPWV-RL和由血压指标计算部108计算出的ABI-RL的交点上标绘的白色三角标记位置，来表示左下腿部的动脉硬化程度。
此外，在本实施方式中，作为解析结果信息，输出所测定的baPWV-RR和baPWV-RL。但是，在所测定的baPWV-RR以及baPWV-RL中，作为脉搏波解析结果，仅输出在图6的步骤S112中计算出的可靠性(稳定度)高的一方的baPWV。例如，由于提示患者用的报告通常只打印简单的信息，因而在患者用的报告中，可以只输出一方的baPWV。在该情况下，代替图9的图表，可打印如图10所示的图表。在图10中，在以baPWV为纵轴、ABI为横轴的图表中，仅标绘出了稳定度高的一方的baPWV。这样，由于仅输出可靠性(稳定度)高的一方的baPWV，比只标绘值高的一方的baPWV及左右的baPWV的平均值的情况，可进行更适当的判断和诊断。
(近似度的计算公式)
如上所述，在本实施方式中，为了对高精度解析指标进行计算，采用与每个脉拍的累积脉搏波的近似度。因此，需要通过实验适当地确定近似度的计算公式。
图11A至图27B是表示测定ID1至ID17各自的每个脉拍的脉搏波形状以上升位置为起点重叠的例子的图。在各个图中，图11A至图27A表示沿着时间轴测定的多个脉拍的脉搏波波形。图11B至图27B表示图11A至图27A所示的各个脉拍的脉搏波形状以上升位置作为起点重叠的情况。图11A至图27A的纵轴值表示对压力进行数字转换的输出值。图11B至图27B的纵轴值表示振幅。另外，图11B至图27B的横轴的分割符表示抽样的点。
图28是表示将图11A至图27B的脉搏波波形作为对象的情况下的，由装置计算的近似度的位次和由观测员观测的近似程度的位次之间的关系的图。
图28中的“指数”是，对应每个测定ID，通过预定的换算公式用百分比表示的、通过上述计算公式(1)计算出的每个脉拍的(与未图示的累积形状的)近似度的平均值。“装置的位次”是，用百分比表示的、在显示的指数的升序上的位次。“观测员位次”表示，由医疗从业人员及脉搏波解析装置100的开发人员等具有丰富的脉搏波知识的人员，对应每个测定ID，用手册判断每个脉拍的(与未图示的累积形状的)平均近似度的位次。
图29是表示图28中所示的、基于装置的近似的位次和基于观测员的近似的位次之间的相关关系的图。如图29所示，在二维的坐标平面上，取Y轴上的“人员的判定位次(观测员位次)”与X轴上的“装置的判定位次”的情况下的与两者相关的测定系数R2用0.6844表示。
这样，在通过上述公式(1)计算近似度的情况下，虽可得到与观测员位次相近的结果，但近似度的计算公式是使观测员位次与装置的判定位次相匹配的公式，即，希望通过实验来确定使测定係数R2接近1.0的公式。这是因为以往是由医疗从业等人员通过肉眼来判定波形的紊乱(心律不齐及身体移动)。
(变形例1)
在上述实施方式中，计算baPWV时，上臂的测定部位是将定为默认的、一边的手臂(例如右上臂)定为测定部位。或者，根据右上臂和左上臂的血压差判定将右上臂作为用于计算baPWV的上臂还是将左上臂作为用于计算baPWV的上臂。
但是，还可以基于解析处理部104计算出的近似度，来确定用于计算baPWV的上臂的测定部位。
具体地说，在图6的步骤S110(计算脉搏波传播速度)中，也可进行如下处理。在下面的说明中，假设用右上臂的脉搏波计算出的baPWV为默认的。
解析处理部104，在前一个步骤S108(排除处理)中，判断右上臂的脉搏波之中的与未从计算对象中排除的脉拍“BTr”的近似度有关的平均值“AVr”是否为预先确定的阈值以上。在判断为平均值AVr小于阈值的情况下，用左上臂的脉搏波计算baPWV。或者，若判断为左上臂的脉搏波之中的与未从计算对象中排除的脉拍“BTI”的近似度有关的平均值“AVI”也小于预先确定的阈值，则可告知用户重新测定。
或者，对平均值AVr与平均值AVI进行比较，可将其值高的一方的部位的脉搏波用于baPWV的计算。
在图7的例子的情况下，上述的脉拍BTr表示除与右上臂的累积脉搏波的近似度低的第3脉拍之外的脉拍。或者，脉拍BTr也可考虑脚踝脉搏波中的排除结果，而表示除第3脉拍、第5脉拍以及第6脉拍之外的脉拍。同样，上述的脉拍BTI表示除与左上臂的累积脉搏波的近似度低的第4脉拍之外的脉拍。或者，脉拍BTI也可考虑脚踝脉搏波除外结果，而表示除第4脉拍至第6脉拍之外的脉拍。
脚踝脉搏波与上臂脉搏波一样，仅用左右任意一边的脉搏波计算baPWV，根据脉搏波处于稳定的上臂脉搏波以及脚踝脉搏波仅计算出1个baPWV。
(变形例2)
在上述实施方式中，虽测定baPWV-RL以及baPWV-RR，但也可根据一个测定部位的脉搏波来计算PWV。
在根据一个测定部位的脉搏波来计算PWV的情况下，通过将脉搏波传播时间(PTT)除以脉搏波的传播距离(Lpt)来求出PWV。传播距离为，所有称之为躯体高度的、心脏与作为反射部位的髂内动脉分叉部之间的距离的2倍。躯体高度是与身高成正比的高度。虽不能直接测定脉搏波传播距离，但可用预定的换算公式求出。将图30所示的Tpp和TR代入预定的换算公式来计算PTT。Tpp表示作为行波的射血波的顶点(最大点)的出现时间和反射波顶点(最大点)的出现时间的时间间隔。TR表示射血波的出现时间和行波从髂内动脉的分叉部反射回来的反射波的出现时间之间的时间间隔。另外，这些也可作为用于进行动脉硬化度的判定的指标。在图30中，用作为射血波顶点的A点和作为反射波顶点的B点之间的时间间隔表示Tpp。在图30中，用从射血波的上升点至反射波的上升点的时间间隔表示TR。
或者，可计算出AI作为解析指标。在该情况下，参照图30，解析处理部104，提取相对于A点的振幅P1的B点的振幅P2，将振幅P1除以振幅P2得到AI值。
在计算这种PWV及AI作为解析指标的情况下，作为影响计算指标的范围，限定在包含从脉搏的上升点至反射波顶点的范围内，可计算近似度。
(变形例3)
虽然说明的是上述实施方式的脉搏波解析装置100包含检测单元20、袖带24以及信息处理单元1，但脉搏波解析装置也可以在不包含检测单元20以及袖带24的通用计算机中实现。就是说，在本实施方式中，脉搏波解析装置若包含通常由CPU10实现的解析处理部104的功能，则能实现如图6所示的脉搏波解析处理。通用计算机具有与信息处理单元1相同的硬盘。
可以将本实施方式以及各个变形例的脉搏波解析装置所执行的脉搏波解析方法作为程序提供。这种程序被记录在计算机可读取的、非暂时性(non-transitory)的记录介质上。例如，这种“计算机可读取的记录介质”包含CD-ROM(Compact Disc-ROM：只读光盘)等光学介质、存储卡之类的磁性记录介质等。另外，可以使这些程序记录在计算机可读取的记录介质中，作为程序产品提供。另外，也可通过网络下载的形式提供这些程序。
此外，可以将与本实施方式相关的程序设为，在规定的排列、规定的时间序列读出作为计算机操作系统(OS)的一部分提供的程序模块之中的所需的模块并执行处理。在这种情况下，程序本身不包含上述模块，而是与OS合作执行处理。不包含这种模块的程序也包括在与本实施方式相关的程序内。
另外，提供与本实施方式相关的程序可以并入其他程序的一部分。在这种情况下，程序本身不包含包括上述其他程序的模块，而与其他程序合作执行处理。这种并入其他程序的程序也包括在与本实施方式相关的程序内。
所提供的程序产品被安装在硬盘等程序存储部中而被执行。此外，程序产品包含程序本身和存储有程序的记录介质。
本次公开的上述实施方式以及它们的变形例在所有方面都是例示性的，并不是限制性的。本发明的技术范围由权利要求书限定，而不是由上述说明限定，并且包括与权利要求书等同的意思或等同范围内的所有变更。
产业上的可利用性
根据本发明，脉搏波解析装置仅采用身体移动等影响少的、稳定的脉拍便能够计算脉搏波解析指标。
附图标记的说明
1 信息处理单元，
2 控制部，
4 输出部，
6 操作部，
8 存储装置，
20ar，20al，20br，20bl 检测单元，
22al 22ar，22bl，22br 配管
24ar，24al，24br，24bl 袖带，
25al，25ar，25bl，25br 压力泵，
26al，26ar，26bl，26br 调压阀，
27al，27ar，27bl，27br 配管，
28al，28ar，28bl，28br 压力传感器，
29al，29ar，29bl，29br A/D转换部，
30 调整部，
100 脉搏波解析装置，
102 脉搏波浪1定部，
104 解析处理部，
106 血压测定部，
108 血压指标计算部，
200 被测定者。

Claims (9)

1.一种脉搏波解析装置(1)，其特征在于，
该脉搏波解析装置(1)具有：
存储部(8)，其用于存储多个脉拍的脉搏波波形，
解析处理部(2)，其进行特定处理，该特定处理是指，通过对所述多个脉拍的脉搏波波形进行解析，来计算脉搏波解析指标的处理；
所述解析处理部(2)对用于构成所述多个脉拍的脉搏波波形的每个脉拍的脉搏波形状进行累积，并从计算对象中排除所累积的脉搏波形状与所述每个脉拍的脉搏波形状之间的近似度低的脉拍，来计算所述脉搏波解析指标；
所述脉搏波解析装置(1)还具有输出部(4)，该输出部(4)用于将计算出的所述脉搏波解析指标作为解析结果来输出。
2.根据权利要求1所述的脉搏波解析装置(1)，其特征在于，
所述解析处理部(2)还对用于计算所述脉搏波解析指标的脉搏波形状的所述近似度进行累积，由此计算搏动的稳定度；
所述输出部(4)还将所述稳定度作为表示所述脉搏波解析指标的可靠性的指标来输出。
3.根据权利要求2所述的脉搏波解析装置(1)，其特征在于，
所述存储部(8)对应每个肢体部位来存储所述多个脉拍的脉搏波波形；
所述解析处理部(2)对应每个肢体部位来累积所述每个脉拍的脉搏波形状，并计算所述近似度、所述脉搏波解析指标以及所述稳定度；
所述输出部(4)将所述稳定度高的所述脉搏波解析指标作为所述解析结果来输出。
4.根据权利要求1所述的脉搏波解析装置(1)，其特征在于，
所述存储部(8)针对左右的肢体部位存储多个脉拍的脉搏波波形；
所述解析处理部(2)对应所述每个肢体部位来计算所述近似度，并利用所述近似度高的肢体部位的脉搏波形状，来计算所述脉搏波解析指标。
5.根据权利要求1所述的脉搏波解析装置(1)，其特征在于，所述解析处理部(2)仅限于在所述每个脉拍的脉搏波形状中对所述脉搏波解析指标的计算产生影响的范围内，计算所述近似度。
6.根据权利要求1所述的脉搏波解析装置(1)，其特征在于，所述脉搏波解析指标表示动脉硬化的程度和/或血管的狭窄程度。
7.根据权利要求6所述的脉搏波解析装置(1)，其特征在于，所述脉搏波解析指标包含作为表示动脉硬化的程度的指标的脉搏波传播速度。
8.根据权利要求1所述的脉搏波解析装置(1)，其特征在于，
还具有脉搏波检测部(20ar、20al、20br、20bl)，所述脉搏波检测部(20ar、20al、20br、20bl)用于对肢体部位的脉搏波进行检测；
所述脉搏波解析部(2)基于来自所述脉搏波检测部(20ar、20al、20br、20bl)的检测信号，来测定所述多个脉拍的脉搏波波形。
9.一种记录介质，记录有脉搏波解析用程序，该脉搏波解析用程序使计算机发挥用于解析脉搏波的装置的功能，其特征在于，
所述脉搏波解析用程序使所述计算机执行以下步骤：
对用于构成存储于存储部中的多个脉拍的脉搏波波形的每个脉拍的脉搏波形状进行累积的步骤；
从计算对象中排除所累积的脉搏波形状与所述每个脉拍的脉搏波形状之间的近似度低的脉拍，来计算所述脉搏波解析指标的步骤；
将计算出的所述脉搏波解析指标作为解析结果来输出的步骤。
CN201080034208.5A 2009-10-30 2010-10-27 脉搏波解析装置以及脉搏波解析方法 Active CN102469942B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2009-250928 2009-10-30
JP2009250928A JP5287668B2 (ja) 2009-10-30 2009-10-30 脈波解析装置および脈波解析プログラム
PCT/JP2010/069093 WO2011052651A1 (ja) 2009-10-30 2010-10-27 脈波解析装置および記録媒体

Publications (2)

Publication Number Publication Date
CN102469942A true CN102469942A (zh) 2012-05-23
CN102469942B CN102469942B (zh) 2015-11-25

Family

ID=43922076

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201080034208.5A Active CN102469942B (zh) 2009-10-30 2010-10-27 脉搏波解析装置以及脉搏波解析方法

Country Status (6)

Country Link
US (1) US9044145B2 (zh)
JP (1) JP5287668B2 (zh)
CN (1) CN102469942B (zh)
DE (1) DE112010004170T5 (zh)
RU (1) RU2570282C2 (zh)
WO (1) WO2011052651A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104545802A (zh) * 2013-10-15 2015-04-29 日本光电工业株式会社 生物信息显示装置和生物信息显示方法
CN105246399A (zh) * 2013-06-26 2016-01-13 英特尔公司 前导中风风险指标的检测
CN108778104A (zh) * 2016-03-16 2018-11-09 福田电子株式会社 血压脉搏波测定装置及程序
CN112754445A (zh) * 2015-03-31 2021-05-07 深圳市长桑技术有限公司 一种生理参数监测的系统和方法
CN113349751A (zh) * 2020-03-05 2021-09-07 深圳市理邦精密仪器股份有限公司 一种充放气控制方法、装置及医疗设备

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5887836B2 (ja) * 2011-10-28 2016-03-16 オムロンヘルスケア株式会社 測定装置、指標算出方法、および指標算出プログラム
JP2013094261A (ja) * 2011-10-28 2013-05-20 Omron Healthcare Co Ltd 測定装置、評価方法、および評価プログラム
JP6078753B2 (ja) * 2012-05-14 2017-02-15 株式会社 ライフインターフェイス 四肢装着型生体情報計測装置
CN103941873B (zh) * 2014-04-30 2017-05-10 北京智谷睿拓技术服务有限公司 识别方法和设备
JP6226822B2 (ja) * 2014-06-11 2017-11-08 日本光電工業株式会社 生体情報測定装置、作動方法、及びプログラム
JP2016214303A (ja) * 2015-05-14 2016-12-22 日本光電工業株式会社 指標出力方法、指標出力装置および指標出力プログラム
JP6862093B2 (ja) * 2016-03-16 2021-04-21 フクダ電子株式会社 血圧脈波測定装置
WO2017171802A1 (en) * 2016-03-31 2017-10-05 Edwards Lifesciences Corporation Aortic stenosis screening
JP6683034B2 (ja) * 2016-06-24 2020-04-15 オムロンヘルスケア株式会社 血圧脈波測定装置
WO2020061887A1 (zh) * 2018-09-27 2020-04-02 深圳市伊欧乐科技有限公司 心率测量方法、装置及计算机可读存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1504166A (zh) * 2002-12-05 2004-06-16 ŷķ������ҽ����ҵ��ʽ���� 脉搏波测定装置
US6755792B2 (en) * 2001-12-17 2004-06-29 Colin Medical Technology Corporation Arteriosclerosis examining apparatus
CN1537508A (zh) * 2003-04-16 2004-10-20 ŷķ������ҽ����ҵ��ʽ���� 脉搏波测定装置和活体波分析程序产品
JP2006247221A (ja) * 2005-03-11 2006-09-21 Omron Healthcare Co Ltd 脈波検出装置
CN101156771A (zh) * 2007-09-28 2008-04-09 天津市先石光学技术有限公司 基于脉搏波频谱分析提高血管硬度测量精度的方法及装置
US20090018408A1 (en) * 2007-07-05 2009-01-15 Kabushiki Kaisha Toshiba Apparatus and method for processing pulse waves

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10328150A (ja) 1997-05-30 1998-12-15 Nippon Colin Co Ltd 脈波伝播速度情報測定装置
JP3213296B2 (ja) 1999-11-01 2001-10-02 日本コーリン株式会社 脈波伝播速度情報測定装置
RU2207801C2 (ru) * 2001-02-21 2003-07-10 Самарский государственный медицинский университет Способ интегральной оценки регионарного кровообращения конечностей
JP3621379B2 (ja) * 2002-01-09 2005-02-16 コーリンメディカルテクノロジー株式会社 動脈硬化評価装置
JP3643565B2 (ja) * 2002-02-21 2005-04-27 コーリンメディカルテクノロジー株式会社 動脈波形検査装置
JP3918811B2 (ja) * 2003-12-22 2007-05-23 オムロンヘルスケア株式会社 動脈硬化検査装置
RU2270609C1 (ru) * 2005-02-02 2006-02-27 Дмитрий Николаевич Майстренко Способ определения показаний к проведению коррекции нарушений регионарной гемодинамики после реконструктивных операций на артериях нижних конечностей
JP4469746B2 (ja) * 2005-03-29 2010-05-26 株式会社東芝 心拍計測装置及び心拍計測装置の作動方法
JP5255771B2 (ja) 2007-01-15 2013-08-07 フクダ電子株式会社 生体情報処理装置および生体情報処理方法
RU2361527C1 (ru) * 2008-03-24 2009-07-20 ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ УЧРЕЖДЕНИЕ РОССИЙСКИЙ НАУЧНЫЙ ЦЕНТР РАДИОЛОГИИ И ХИРУРГИЧЕСКИХ ТЕХНОЛОГИЙ ФЕДЕРАЛЬНОГО АГЕНТСТВА ПО ВЫСОКОТЕХНОЛОГИЧНОЙ МЕДИЦИНСКОЙ ПОМОЩИ (ФГУ "РНЦРХТ Росмедтехнологий") Способ лечения облитерирующего атеросклероза сосудов нижних конечностей

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6755792B2 (en) * 2001-12-17 2004-06-29 Colin Medical Technology Corporation Arteriosclerosis examining apparatus
CN1504166A (zh) * 2002-12-05 2004-06-16 ŷķ������ҽ����ҵ��ʽ���� 脉搏波测定装置
CN1537508A (zh) * 2003-04-16 2004-10-20 ŷķ������ҽ����ҵ��ʽ���� 脉搏波测定装置和活体波分析程序产品
JP2006247221A (ja) * 2005-03-11 2006-09-21 Omron Healthcare Co Ltd 脈波検出装置
US20090018408A1 (en) * 2007-07-05 2009-01-15 Kabushiki Kaisha Toshiba Apparatus and method for processing pulse waves
CN101156771A (zh) * 2007-09-28 2008-04-09 天津市先石光学技术有限公司 基于脉搏波频谱分析提高血管硬度测量精度的方法及装置

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105246399A (zh) * 2013-06-26 2016-01-13 英特尔公司 前导中风风险指标的检测
CN105246399B (zh) * 2013-06-26 2017-10-03 英特尔公司 前导中风风险指标的检测
CN104545802A (zh) * 2013-10-15 2015-04-29 日本光电工业株式会社 生物信息显示装置和生物信息显示方法
CN104545802B (zh) * 2013-10-15 2019-06-28 日本光电工业株式会社 生物信息显示装置和生物信息显示方法
CN112754445A (zh) * 2015-03-31 2021-05-07 深圳市长桑技术有限公司 一种生理参数监测的系统和方法
CN108778104A (zh) * 2016-03-16 2018-11-09 福田电子株式会社 血压脉搏波测定装置及程序
CN113349751A (zh) * 2020-03-05 2021-09-07 深圳市理邦精密仪器股份有限公司 一种充放气控制方法、装置及医疗设备

Also Published As

Publication number Publication date
US20120095353A1 (en) 2012-04-19
JP2011092556A (ja) 2011-05-12
DE112010004170T5 (de) 2012-10-04
WO2011052651A1 (ja) 2011-05-05
CN102469942B (zh) 2015-11-25
US9044145B2 (en) 2015-06-02
RU2012103188A (ru) 2013-08-10
JP5287668B2 (ja) 2013-09-11
RU2570282C2 (ru) 2015-12-10

Similar Documents

Publication Publication Date Title
CN102469942A (zh) 脉搏波解析装置以及记录介质
EP3427654A1 (en) Biological information analysis device and system, and program
JP6308742B2 (ja) 血圧情報出力装置、血圧情報出力プログラム、媒体、血圧情報出力方法
CN104684459B (zh) 光体积描记术测量中识别运动伪影并且改善测量结果和警报的可靠性的方法和系统
JP4855721B2 (ja) 血圧測定装置
JP2006326293A (ja) 健康状態に応じて指標を与えるために心臓血管の働きを評価するための装置
CN106560155A (zh) 用于测量生物信息的设备和方法
US20130310691A1 (en) System, method and device for automatic and autonomous determination of hemodynamic and cardiac parameters using ultrasound
WO2017089921A1 (en) Method to quantify photoplethysmogram (ppg) signal quality
KR101656740B1 (ko) 동맥 혈압 파형의 특징점 검출 장치 및 방법
US10905381B2 (en) Blood pressure correction information generating device, blood pressure measurement device and blood pressure correction information generating method
US11147461B2 (en) Blood pressure analyzing apparatus, blood pressure measuring apparatus, and blood pressure analyzing method
KR20180108841A (ko) 혈압 맥파 측정 장치 및 프로그램
CN110520048A (zh) 情绪估计装置
CN107106046B (zh) 血管弹性率评价装置
JP6414981B2 (ja) 生体情報推定装置、生体情報推定方法、及びコンピュータプログラム
US20190298190A1 (en) Pulse detection, measurement and analysis based health management system, method and apparatus
JP2022160608A (ja) 生体信号処理装置およびその制御方法
US20200275845A1 (en) Blood-pressure measurement apparatus and blood-pressure measurement method
JP2021069615A (ja) 生体情報測定装置、及び生体情報測定プログラム
JP4453431B2 (ja) 生体情報測定装置
US20080082011A1 (en) Computer Input Apparatus with Health Care Function and Health Care Method
CN112601490B (zh) 用于产生指示心脏异常的信息的装置
EP3939500B1 (en) Apparatus and method for estimating bio-information
CN114246568A (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
TR01 Transfer of patent right

Effective date of registration: 20170630

Address after: Tokyo, Japan, Japan

Patentee after: Fukuda Denshi KK

Address before: Kyoto Japan

Patentee before: Omron Healthcare Co., Ltd.