JP2007105131A - 脈波診断装置及び脈波診断装置制御方法 - Google Patents

脈波診断装置及び脈波診断装置制御方法 Download PDF

Info

Publication number
JP2007105131A
JP2007105131A JP2005297097A JP2005297097A JP2007105131A JP 2007105131 A JP2007105131 A JP 2007105131A JP 2005297097 A JP2005297097 A JP 2005297097A JP 2005297097 A JP2005297097 A JP 2005297097A JP 2007105131 A JP2007105131 A JP 2007105131A
Authority
JP
Japan
Prior art keywords
pulse wave
pulse
poincare
wave amplitude
interval
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.)
Pending
Application number
JP2005297097A
Other languages
English (en)
Inventor
Hiroshi Koizumi
弘 小泉
Shinji Mino
真司 美野
Junichi Shimada
純一 嶋田
Naoyoshi Tatara
尚愛 多々良
Shoichi Hayashida
尚一 林田
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.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone 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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2005297097A priority Critical patent/JP2007105131A/ja
Publication of JP2007105131A publication Critical patent/JP2007105131A/ja
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

【課題】本発明は、長時間にわたる脈波の変動を検出することのできる脈波診断装置及び脈波診断装置制御方法の提供を目的とする。
【解決手段】本発明に係る脈波診断装置及び脈波診断装置制御方法は、光電センサを用いて脈波を検出し、検出した脈波から脈波の変動を算出することを特徴とする。具体的には、脈拍間隔のゆらぎの周波数成分を抽出する周波数成分抽出部か、連続して算出された２つの脈拍間隔同士で形成される直交座標平面上での脈拍間隔の点をポアンカレ座標として算出する脈拍間隔ポアンカレ算出部か、脈波振幅を算出する脈波振幅算出部か、連続して算出された２つの脈波振幅同士で形成される直交座標平面上での脈波振幅の点をポアンカレ座標として算出する脈波振幅ポアンカレ算出部か、を備えることを特徴とする。
【選択図】図１

Description

本発明は、脈波を検出して脈波間隔の変動や脈波振幅の変動を計測する脈波診断装置及び脈波診断装置制御方法に関する。
心拍変動が精神的ストレスの指標となることが知られており、この心拍変動を用いた精神状態判定装置が提案されている（例えば、非特許文献１又は特許文献１参照。）。従来の心拍変動を用いた判定装置は、電極を皮膚表面に装着して心電を検出する。電極を皮膚表面に装着する際、皮膚との接触抵抗を下げるために電極と皮膚との間隙をジェルなどの整合剤で充填する必要がある。このため、整合剤と皮膚表面とが長時間接触するので、長時間心電を検出するとかぶれなどの皮膚疾患を起こしやすかった。また、時間経過により整合剤が乾燥するので、長時間心電を検出すると電極が剥離する問題もあった。また、整合剤の乾燥とは逆に、被検者の発汗によっても電極が剥離しやすくなるという問題があった。さらに、心電を用いて心拍変動を検出するためには、このような電極を少なくとも３個装着する必要があるので、長時間にわたる心拍変動の検出が困難であった。
特開２００４−３４４２６９号公報 林博史編集：「心拍変動の臨床応用」医学書院、１５頁〜１６頁
上記のように、従来は、皮膚表面に電極を装着して心電を検出していたので、長時間にわたる心拍変動の検出が困難であった。心拍変動は脈波の変動によっても検出することができる。よって、本発明は、長時間にわたる脈波の変動を検出することのできる脈波診断装置及び脈波診断装置制御方法の提供を目的とする。
上記課題を解決するために、本発明に係る脈波診断装置及び脈波診断装置制御方法は、光電センサを用いて脈波を検出して脈波の変動を算出することを特徴とする。
光電センサを用いて脈波を検出するので、ジェルなどの整合剤を用いることなく脈波を検出することができる。また、光電センサを用いて脈波を検出するので、人体の外耳などの抹消血管から脈波を検出することができる。よって、長時間にわたる脈波の変動を検出することができる。
具体的には、本発明に係る脈波診断装置は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、前記光電脈波検出部の検出した脈波の１拍ごとの脈拍間隔を算出し、前記脈拍間隔のゆらぎの周波数成分を抽出する周波数成分抽出部と、を備えることを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。周波数成分抽出部が脈拍間隔のゆらぎの周波数成分を抽出するので、長時間にわたる脈波の変動を検出することができる。
前記周波数成分抽出部は、前記脈拍間隔のゆらぎの周波数成分を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに抽出し、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部をさらに備え、前記周波数成分抽出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈拍間隔のゆらぎの周波数成分を抽出することが好ましい。
ゆらぎを算出する脈拍数を可変できるので、所望のデータ数での脈拍間隔の変動を算出することができる。
前記周波数成分抽出部の抽出するゆらぎの周波数成分をスペクトル表示する表示部をさらに備えることが好ましい。
ゆらぎの周波数成分をスペクトル表示することができるので、長時間にわたる脈拍間隔の変動を表示することができる。
前記周波数成分抽出部の抽出するゆらぎの周波数成分のうち特定の２つの周波数の周波数成分の比を算出する周波数成分比算出部をさらに備えることが好ましい。
特定の２つの周波数の周波数成分の比を算出するので、脈拍間隔のゆらぎの周波数成分の変動を算出することができる。
前記周波数成分比算出部は、前記周波数成分抽出部の抽出するゆらぎの周波数成分のうち０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比を算出することが好ましい。
０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比の変化を算定することができる。０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比の変化は精神的ストレスの判定に役立つとされているので、精神的ストレスの推定に資することができる。
また、本発明に係る脈波診断装置は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、前記光電脈波検出部の検出する脈波の１拍ごとの脈拍間隔を算出し、連続して算出された２つの前記脈拍間隔同士で形成される直交座標平面上での前記脈拍間隔の点をポアンカレ座標として算出する脈拍間隔ポアンカレ算出部と、を備えることを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。脈拍間隔ポアンカレ算出部が脈拍間隔のポアンカレ座標を算出するので、時刻に依存しない脈拍間隔の分布を算出することができる。
前記脈拍間隔ポアンカレ算出部は、前記脈拍間隔のポアンカレ座標を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部をさらに備え、前記脈拍間隔ポアンカレ算出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈拍間隔のポアンカレ座標を算出することが好ましい。
脈拍間隔のポアンカレ座標を算出する脈拍数を可変できるので、所望のデータ数での脈拍間隔のポアンカレ座標を算出することができる。
前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標を表示する表示部をさらに備えることが好ましい。
脈拍間隔のポアンカレ座標を表示することができるので、時刻に依存しない脈拍間隔の分布を表示することができる。
前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈拍間隔統計処理部をさらに備えることが好ましい。
ポアンカレ座標の標準偏差、中央値又は平均値を算出することができるので、時刻に依存しない脈拍間隔の分布の統計値を算出することができる。
前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出する脈拍間隔分布算出部をさらに備えることが好ましい。
脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出すれば、脈拍間隔のポアンカレ座標の分布が狭まったことを検出できる。
前記脈拍間隔分布算出部は、脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することが好ましい。
予め定められた座標範囲内での最大及び最小の変化幅を算出するので、脈拍間隔の診断に有用な座標範囲内での変化幅を検出することができる。
また、本発明に係る脈波診断装置は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、前記光電脈波検出部の検出する脈波の脈波振幅を算出する脈波振幅算出部と、前記脈波振幅算出部の算出する脈波振幅が所定の大きさ以上か否かを判定する判定部と、を備えることを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。脈波振幅算出部が脈波振幅を算出するので、長時間にわたる脈波の変動を検出することができる。さらに、光電センサを用いて血管の容積変化を検出できるので、脈波振幅を算出することによって血管のスティフネス等の血管の状態の変化を検出することができる。
前記脈波振幅算出部は、前記光電脈波検出部の検出する脈波の脈波振幅の平均値を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部をさらに備え、前記脈波振幅算出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈波振幅の平均値を算出することが好ましい。
脈波振幅の平均値を算出する脈拍数を可変できるので、所望のデータ数での脈波振幅の変動を算出することができる。
前記脈波振幅算出部の算出する脈波振幅の時間推移か、又は、前記脈波振幅算出部の算出する脈波振幅の平均値の時間推移を表示する表示部をさらに備えることが好ましい。
脈波振幅の時間推移又は脈波振幅の平均値の時間推移を表示することができるので、長時間にわたる脈波振幅の変動を表示することができる。
前記脈波振幅算出部の算出する平均値のうち特定の２つの平均値の比を算出する脈波振幅比算出部をさらに備えることが好ましい。
脈波振幅算出部の算出する２つの脈波振幅の平均値の比を算出するので、２つの脈波振幅の平均値の増減を判定することができる。
また、本発明に係る脈波診断装置は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、前記光電脈波検出部の検出する脈波の１拍ごとの脈波振幅を算出し、連続して算出された２つの前記脈波振幅同士で形成される直交座標平面上での前記脈波振幅の点をポアンカレ座標として１拍ごとに算出する脈波振幅ポアンカレ算出部と、を備えることを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。脈波振幅ポアンカレ算出部が脈波振幅のポアンカレ座標を算出するので、時刻に依存しない脈波振幅の分布を算出することができる。
前記脈波振幅ポアンカレ算出部は、脈波振幅のポアンカレ座標を予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、外部からの入力により前記第１脈拍数を新たな第１脈拍数として設定する脈拍数可変部をさらに備え、前記脈波振幅ポアンカレ算出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈波振幅のポアンカレ座標を算出することが好ましい。
脈波振幅のポアンカレ座標を算出する脈拍数を可変できるので、所望のデータ数での脈波振幅のポアンカレ座標を算出することができる。
前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標を表示する表示部をさらに備えることが好ましい。
脈波振幅のポアンカレ座標を表示することができるので、時刻に依存しない脈波振幅の分布を表示することができる。
前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理部をさらに備えることが好ましい。
ポアンカレ座標の標準偏差、中央値又は平均値を算出するので、時刻に依存しない脈波振幅の分布の統計値を算出することができる。
前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出する脈波振幅分布算出部をさらに備えることが好ましい。
脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出すれば、脈波振幅のポアンカレ座標の分布が狭まったことを検出できる。
前記脈波振幅分布算出部は、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することが好ましい。
予め定められた座標範囲内での最大及び最小の変化幅を算出するので、脈波振幅の診断に有用な座標範囲内での変化幅を検出することができる。
前記光電脈波検出部は、人体の外耳への装着機構をさらに備えることが好ましい。
装着機構を備えるので、光電脈波検出部を外耳に装着し、外耳の脈波を検出することができる。
前記光電脈波検出部は、検出する脈波を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。
光電脈波検出部が５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で脈波を抽出するので、脈波の変化を精度よく識別することができる。
本発明に係る脈波診断装置制御方法は、光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、周波数成分抽出部が、前記光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を抽出し、前記脈拍間隔のゆらぎの周波数成分を抽出する周波数成分抽出過程と、を有することを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。周波数成分抽出部が脈拍間隔のゆらぎの周波数成分を抽出するので、長時間にわたる脈波の変動を検出することができる。
前記周波数成分抽出過程において、前記周波数成分抽出部が、前記脈拍間隔のゆらぎの周波数成分を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに抽出し、脈拍数可変部が、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変過程を前記周波数成分抽出過程の前にさらに有し、前記周波数成分抽出過程において、前記周波数成分抽出部が、前記脈拍数可変過程で設定した新たな第１脈拍数を用いて脈拍間隔のゆらぎの周波数成分を抽出することが好ましい。
ゆらぎを算出する脈拍数を可変できるので、所望のデータ数での脈拍間隔の変動を算出することができる。
表示部が、前記周波数成分抽出過程で抽出したゆらぎの周波数成分をスペクトル表示する表示過程をさらに有することが好ましい。
ゆらぎの周波数成分をスペクトル表示することができるので、長時間にわたる脈拍間隔の変動を表示することができる。
周波数成分比算出部が、前記周波数成分抽出過程で抽出したゆらぎの周波数成分のうち特定の２つの周波数の周波数成分の比を算出する周波数成分比算出過程をさらに有することが好ましい。
特定の２つの周波数の周波数成分の比を算出するので、脈拍間隔のゆらぎの周波数成分の変動を算出することができる。
前記周波数成分比算出過程において、前記周波数成分比算出部が、前記周波数成分抽出過程で抽出した周波数成分のうち０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比を算出することが好ましい。
０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比の変化を算定することができる。０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比の変化は精神的ストレスの判定に役立つとされているので、精神的ストレスの推定に資することができる。
また、本発明に係る脈波診断装置制御方法は、光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、脈拍間隔ポアンカレ算出部が、前記光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を算出し、連続して算出された２つの前記脈拍間隔同士で形成される直交座標平面上での前記脈拍間隔の点をポアンカレ座標として算出する脈拍間隔ポアンカレ算出過程と、を有することを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。脈拍間隔ポアンカレ算出部が脈拍間隔のポアンカレ座標を算出するので、時刻に依存しない脈拍間隔の分布を算出することができる。
前記脈拍間隔ポアンカレ算出過程において、前記脈拍間隔ポアンカレ算出部が、前記脈拍間隔のポアンカレ座標を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、脈拍数可変部が、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変過程を前記脈拍間隔ポアンカレ算出過程の前にさらに有し、前記脈拍間隔ポアンカレ算出過程において、前記脈拍間隔ポアンカレ算出部が、前記脈拍数可変過程で設定した新たな第１脈拍数を用いて脈拍間隔のポアンカレ座標を算出することが好ましい。
脈拍間隔のポアンカレ座標を算出する脈拍数を可変できるので、所望のデータ数での脈拍間隔のポアンカレ座標を算出することができる。
表示部が、前記脈拍間隔ポアンカレ算出過程で算出した脈拍間隔のポアンカレ座標を表示するポアンカレプロット表示過程をさらに有することが好ましい。
脈拍間隔のポアンカレ座標を表示することができるので、時刻に依存しない脈拍間隔の分布を表示することができる。
脈拍間隔統計処理部が、前記脈拍間隔ポアンカレ算出過程で算出した脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈拍間隔統計処理過程をさらに有することが好ましい。
ポアンカレ座標の標準偏差、中央値又は平均値を算出することができるので、時刻に依存しない脈拍間隔の分布の統計値を算出することができる。
前記脈拍間隔分布算出部が、前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出する脈拍間隔分布算出過程をさらに有することが好ましい。
脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出すれば、脈拍間隔のポアンカレ座標の分布が狭まったことを検出できる。
前記脈拍間隔分布算出過程において、前記脈拍間隔分布算出部が、前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することが好ましい。
予め定められた座標範囲内での最大及び最小の変化幅を算出するので、脈拍間隔の診断に有用な座標範囲内での変化幅を検出することができる。
また、本発明に係る脈波診断装置制御方法は、光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、脈波振幅算出部が、前記光電脈波検出過程で検出した脈波の脈波振幅を算出する脈波振幅算出過程と、判定部が、前記脈波振幅算出過程で算出した脈波振幅が所定の大きさ以上か否かを判定する判定過程と、を有することを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。脈波振幅算出部が脈波振幅を算出するので、長時間にわたる脈波の変動を検出することができる。さらに、光電センサを用いて血管の容積変化を検出できるので、脈波振幅を算出することによって血管のスティフネス等の血管の状態の変化を検出することができる。
前記脈波振幅算出過程において、前記脈波振幅算出部が、前記光電脈波検出過程で検出した脈波の脈波振幅の平均値を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、脈拍数可変部が、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変過程を前記脈波振幅算出過程の前にさらに有し、前記脈波振幅算出過程において、前記脈波振幅算出部が、前記脈拍数可変過程で設定した新たな第１脈拍数を用いて脈波振幅の平均値を算出することが好ましい。
脈波振幅の平均値を算出する脈拍数を可変できるので、所望のデータ数での脈波振幅の変動を算出することができる。
表示部が、前記脈波振幅算出過程で算出された脈波振幅の時間推移か、又は、前記脈波振幅算出過程で算出された脈波振幅の平均値の時間推移を表示する表示過程をさらに有することが好ましい。
脈波振幅の時間推移又は脈波振幅の平均値の時間推移を表示することができるので、長時間にわたる脈波振幅の変動を表示することができる。
脈波振幅比算出部が、前記脈波振幅算出過程で算出した平均値のうち特定の２つの平均値の比を算出する脈波振幅比算出過程をさらに有することが好ましい。
脈波振幅算出部の算出する２つの脈波振幅の平均値の比を算出するので、２つの脈波振幅の平均値の増減を判定することができる。
また、本発明に係る脈波診断装置制御方法は、光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、脈波振幅ポアンカレ算出部が、前記光電脈波検出過程で検出した脈波の１拍ごとの脈波振幅を算出し、連続して算出された２つの前記脈波振幅同士で形成される直交座標平面上での前記脈波振幅の点をポアンカレ座標として１拍ごとに算出する脈波振幅ポアンカレ算出過程と、を有することを特徴とする。
光電センサを用い、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。脈波振幅ポアンカレ算出部が脈波振幅のポアンカレ座標を算出するので、時刻に依存しない脈波振幅の分布を算出することができる。
前記脈波振幅ポアンカレ算出過程において、前記脈波振幅ポアンカレ算出部が、前記ポアンカレ座標を予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、脈拍数可変部が、外部からの入力により前記第１脈拍数を新たな第１脈拍数として設定する脈拍数可変過程を前記脈波振幅ポアンカレ算出過程の前にさらに有し、前記脈波振幅ポアンカレ算出過程において、前記脈波振幅ポアンカレ算出部が、前記脈拍数可変過程で設定する新たな第１脈拍数を用いて脈波振幅のポアンカレ座標を算出することが好ましい。
脈波振幅のポアンカレ座標を算出する脈拍数を可変できるので、所望のデータ数での脈波振幅のポアンカレ座標を算出することができる。
表示部が、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標を表示する表示過程をさらに有することが好ましい。
脈波振幅のポアンカレ座標を表示することができるので、時刻に依存しない脈波振幅の分布を表示することができる。
脈波振幅統計処理部が、前記脈波振幅ポアンカレ算出過程で算出した脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理過程をさらに有することが好ましい。
ポアンカレ座標の標準偏差、中央値又は平均値を算出するので、時刻に依存しない脈波振幅の分布の統計値を算出することができる。
前記脈波振幅分布算出部が、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出する脈波振幅分布算出過程をさらに有することが好ましい。
脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出すれば、脈波振幅のポアンカレ座標の分布が狭まったことを検出できる。
前記脈波振幅分布算出過程において、前記脈波振幅分布算出部が、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することが好ましい。
予め定められた座標範囲内での最大及び最小の変化幅を算出するので、脈波振幅の診断に有用な座標範囲内での変化幅を検出することができる。
前記光電脈波検出過程において、前記光電脈波検出部が、検出する脈波を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。
光電脈波検出部が５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で脈波を抽出するので、脈波の変化を精度よく識別することができる。
本発明によれば、光電センサを用いて脈波を検出し、検出した脈波の脈拍間隔の変動又は脈波振幅の変動を算出するので、長時間にわたる脈波の変動を検出することのできる脈波診断装置及び脈波診断装置制御方法を提供することができる。
添付の図面を参照して本発明の実施の形態を説明する。以下に説明する実施の形態は本発明の構成の例であり、本発明は、以下の実施の形態に制限されるものではない。
（実施形態１）
図１は、本実施形態に係る脈波診断装置の構成概略図である。図１に示す脈波診断装置９１は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部１１と、光電脈波検出部１１の検出した脈波１の１拍ごとの脈拍間隔を算出し、脈拍間隔のゆらぎの周波数成分を抽出する周波数成分抽出部１２と、を備えることを特徴とする。
脈波診断装置９１は、外部からの入力により、第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部１３と、周波数成分抽出部１２の抽出するゆらぎの周波数成分をスペクトル表示する表示部１４と、周波数成分抽出部１２の抽出するゆらぎの周波数成分のうち特定の２つの周波数の周波数成分の比を算出する周波数成分比算出部１５と、をさらに備えていてもよい。また、所定時刻に光電脈波検出部１１に脈波の検出を始めさせ、所定時刻に光電脈波検出部１１に脈波の検出を終了させるタイマー１６と、周波数成分比算出部１５の算出する周波数成分の比があらかじめ定められた閾値以上であるか否かを判定する判定部１７と、をさらに備えていてもよい。
光電脈波検出部１１は、動脈を透過した透過光を検出するものである。例えば、透過型の光電センサである。透過型の光電センサとしては、例えば、発光素子と、前記発光素子から出射された光が被検者の一部を透過した透過光を受光する受光素子とを用いることができる。また、光電脈波検出部１１は、動脈で散乱された散乱光を受光して脈波を検出するものである。例えば、反射型の光電センサである。反射型の光電センサとしては、例えば、発光素子と、前記発光素子から出射された光が被検者の一部で反射された反射光を受光する受光素子とを備えるものである。発光素子としては、例えば、ＬＥＤ（半導体発光素子）又は半導体レーザを用いることができる。発光素子の出射する光の波長は、例えば、ヘモグロビンで吸収される波長や、血管壁で散乱される波長とすることができる。受光素子としては、例えば、フォトダイオード又はフォトトランジスタを用いることができる。発光素子、受光素子とも複数有してもよい。光電センサを用いて脈波を検出することにより、被検体の振動の影響の少ない脈波検出ができる。ここで、被検体とは、例えば、耳介１０１である。また、末梢血管での脈波検出の精度を向上することができる。
光電脈波検出部１１は、人体の外耳１０３への装着機構３１をさらに備えることが好ましい。装着機構３１としては、例えば、挟持する機構を用いることができる。また、人体の外耳道１０２の内壁を押圧することによって脱離を防ぐ機構であってもよい。また、装着機構３１は、耳介１０１の上部に掛ける機構を用いてもよい。また、外耳道１０２に挿入する機構を用いてもよい。このように、光電脈波検出部１１が装着機構３１を備えれば、外耳１０３の脈波を検出することができる。
光電脈波検出部１１は、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。脈波には１００Ｈｚまでの脈波の周波数成分が含まれている。よって、ナイキスト条件から、サンプリング間隔を５ｍｓｅｃ以下とすることによって、脈波を精度よく検出することができる。５００μｓｅｃ未満とすると、脈波と無関係の微細なノイズを検出する可能性がある。さらに、光電脈波検出部１１は、脈波１を１ｍｓｅｃのサンプリング間隔で抽出することが好ましい。５００μｓｅｃ以上１ｍｓｅｃ以下であれば、脈波の微細な変動を検出することができる。このように、光電脈波検出部１１のサンプリング間隔が５００μｓｅｃ以上５ｍｓｅｃ以下であれば、光電脈波検出部１１の検出した脈波１を精度よく検出することができる。
周波数成分抽出部１２の構成について図２を用いて説明する。図２は、光電脈波検出部の検出する脈波の一例を示す模式図である。横軸は時間、縦軸は光電脈波検出部１１の出力する脈波１の信号出力を示す。図２に示す脈波信号出力は、山４１ａ、谷４２ａ、山４１ｂ、谷４２ｂ、山４１ｃ、谷４２ｃ、山４１ｄ及び谷４２ｄが連続している。
周波数成分抽出部１２は、光電脈波検出部１１の検出した脈波１の１拍ごとの脈拍間隔を算出するものである。脈拍間隔は、例えば、隣接する山４１同士の時間間隔４３（Ｒ−Ｒ間隔）である。また、脈拍間隔は、隣接する谷４２同士の時間間隔（Ｐ−Ｐ間隔）としてもよい。又、脈拍間隔は、脈波の立ち上がり点の間隔としてもよい。又、脈拍間隔は、脈波が予め定めた閾値を超える時点の間隔としてもよい。
さらに、図１に示す周波数成分抽出部１２は、算出した脈拍間隔から、脈拍間隔のゆらぎの周波数成分を抽出するものである。図３は、脈拍間隔のゆらぎの周波数成分の一例を示すグラフである。横軸が周波数、縦軸が周波数成分である。ここで、周波数成分は、例えばパワースペクトル密度とすることができる。
さらに、周波数成分抽出部１２は、脈拍間隔のゆらぎの周波数成分を、予め定められた第１脈拍数で抽出することが好ましい。第１脈拍数とは、予め定められた脈拍数であり、例えば脈拍間隔のゆらぎを算出するデータ数である。例えば、２５６〜５１２拍程度の連続したデータである。さらに、第１脈拍数は、脈拍数可変部１３によって増減可能であることが好ましい。第１脈拍数が増減可能であれば、脈拍間隔のゆらぎを算出する時間幅を可変できるので、所望のデータ数での脈拍間隔の変動を算出することができる。さらに、周波数成分抽出部１２は、脈拍間隔のゆらぎの周波数成分を、予め定められた第２脈拍数ごとに抽出することが好ましい。第２脈拍数とは、予め定められた脈拍数であり、例えば脈拍間隔のゆらぎを算出するタイミングである。第２脈拍数は、例えば、１拍とすることができる。また、１０分程度の間にわたり検出される脈拍数としてもよい。例えば、５００拍以上１０００拍以下の範囲内のいずれかの脈拍数とすることができる。また、６０分程度の間にわたり検出される脈拍数としてもよい。また、第２脈拍数は、第１脈拍数と同一であってもよい。また、第２脈拍数も可変であってもよい。例えば、１０００拍と６０００拍となどの所定の脈拍数が数段にわたり切り替えられることが好ましい。拍間隔のゆらぎの周波数成分を第２脈拍数ごとに算出することができるので、最新のゆらぎの周波数成分を第２脈拍数ごとに取得することができる。
脈波診断装置９１は、外部からの入力により、第１脈拍数を増加又は減少して新たな第１脈拍数を設定する脈拍数可変部１３をさらに備えることが好ましい。脈拍数可変部１３は、例えば、外部から入力された数値を新たな第１脈拍数として設定するものである。又、脈拍数可変部１３は、複数の脈拍数を記憶しているテーブルを格納しており、外部からの入力信号によりテーブルに記憶されている脈拍数のいずれかを選択し、新たな第１脈拍数として設定するものでもよい。
図１に示す脈波診断装置９１は、周波数成分抽出部１２の抽出するゆらぎの周波数成分のうち特定の２つの周波数の周波数成分の比を算出する周波数成分比算出部１５をさらに備えることが好ましい。特定の２つの周波数とは、あらかじめ定められた周波数である。例えば、０．１Ｈｚと０．３Ｈｚである。特定の２つの周波数は、０．１Ｈｚと０．３Ｈｚに限定するものではなく、０．１Ｈｚ付近と０．３Ｈｚ付近であってもよい。０．１Ｈｚ付近と０．３Ｈｚ付近とは、例えば、０．０５〜０．１５Ｈｚの範囲のいずれかの周波数と、０．２〜０．４Ｈｚの範囲のいずれかの周波数であってもよい。このように、０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比の変化を算定することができる。０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比の変化は精神的ストレスの判定に役立つとされているので、精神的ストレスの推定に資することができる。
図１に示す脈波診断装置９１は、所定時刻に光電脈波検出部１１に脈波の検出を始めさせ、所定時刻に光電脈波検出部１１に脈波の検出を終了させるタイマー１６をさらに備えていてもよい。
図１に示す脈波診断装置９１は、周波数成分比算出部１５の算出する脈拍間隔の周波数成分があらかじめ定められた閾値以上であるか否かを判定する判定部１７を備えていてもよい。この場合、あらかじめ定められた閾値とは、例えば、ある時間安静を保ったときに検出された脈波の脈拍間隔の周波数成分とすることができる。閾値は、被検者に応じて可変であることが好ましい。例えば、使用する被検者の脈拍間隔の周波数成分を予め算出し、メモリーに記憶させることができる。判定部１７は、メモリーから閾値を取得する。このように、判定部１７が、周波数成分比算出部１５の算出する脈拍間隔の周波数成分の変動を、記憶されている脈拍間隔の周波数成分を基準に判定すれば、医師又は被検者は、精神的ストレスの有無又はその強度を推定することができる。
脈波診断装置９１は、周波数成分抽出部１２の抽出する脈拍間隔のゆらぎの周波数成分をスペクトル表示する表示部１４をさらに備えることが好ましい。表示部１４としては、例えば、液晶表示装置又はＬＥＤディスプレイを用いることができる。表示部１４は、さらに、脈拍数可変部１３の設定した新たな第１脈拍数を表示してもよい。又、第２脈拍数を表示してもよい。又、周波数成分比算出部１５の算出する特定の２つの周波数の周波数成分の比をさらに表示してもよい。又、判定部１７の判定結果を表示してもよい。このように、表示部１４は、ゆらぎの周波数成分をスペクトル表示することができるので、長時間にわたる脈拍間隔の変動を表示することができる。
次に、脈波診断装置９１の動作について説明する。本実施形態に係る脈波診断装置９１を用いた脈波診断装置制御方法は、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波１を検出する光電脈波検出過程と、周波数成分抽出部１２が、光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を抽出し、脈拍間隔のゆらぎの周波数成分を抽出する周波数成分抽出過程と、を有する。
脈波診断装置９１を用いた脈波診断装置制御方法は、脈拍数可変部１３が、脈拍間隔のゆらぎの周波数成分を算出する予め定められた第１脈拍数を、外部からの入力により増減して新たな第１脈拍数を設定する脈拍数可変過程を周波数成分抽出過程の前にさらに有していてもよい。又、表示部１４が、周波数成分抽出過程で抽出したゆらぎの周波数成分をスペクトル表示する表示過程をさらに有していてもよい。周波数成分抽出過程の後に周波数成分比算出過程をさらに有していてもよいし、周波数成分比算出過程の後に判定過程をさらに有していてもよい。
光電脈波検出過程では、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波１を検出する。脈波１の検出は、例えば、光電脈波検出部１１へ信号が入力されたのを受けて脈波１の検出を開始し、再度信号が入力されたのを受けて脈波１の検出を終了する。また、脈波診断装置９１がタイマー１６を備え、タイマー１６からの入力信号によって脈波１の検出を開始又は終了してもよい。
光電脈波検出過程において、光電脈波検出部１１が、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。光電脈波検出部１１が５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で脈波１を抽出するので、ナイキスト条件から、サンプリング間隔を５ｍｓｅｃ以下とすることによって、脈波１を精度よく検出することができる。さらに、光電脈波検出部１１は、脈波１を１ｍｓｅｃのサンプリング間隔で抽出することが好ましい。
周波数成分抽出過程では、周波数成分抽出部１２が、光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を抽出し、脈拍間隔のゆらぎの周波数成分を抽出する。周波数成分抽出部１２が脈拍間隔のゆらぎの周波数成分を抽出するので、脈波１の変動を検出することができる。
脈拍間隔のゆらぎの周波数成分は、精神的ストレスが印加されないときではピークが２つであるのに対し、精神的ストレスが印加されるとピークが１つになるとされている。脈拍間隔のゆらぎの周波数成分のピークが２つから１つになったことを用いて、精神的ストレスの判定に資することができる。例えば、表示部１４を備え、表示部１４がゆらぎの周波数成分をスペクトル表示すれば、ピークが２つから１つになったか、ピークのいずれかが小さくなったことを視認することができる。また、脈拍間隔のゆらぎの周波数成分をコンピュータへ転送し、コンピュータで脈拍間隔のゆらぎの周波数成分のピークの１つが小さくなったことを判定してもよい。このように、周波数成分抽出部１２が脈拍間隔のゆらぎの周波数成分を算出することによって、精神的ストレスの推定に資することができる。
周波数成分抽出過程において、周波数成分抽出部１２が、脈拍間隔のゆらぎの周波数成分を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに抽出することが好ましい。さらに、周波数成分抽出過程において、周波数成分抽出部１２が、脈拍数可変過程で設定した新たな第１脈拍数を用いて脈拍間隔のゆらぎの周波数成分を抽出することが好ましい。第１脈拍数及び第２脈拍数は、前述の周波数成分抽出部１２で説明したものと同様とすることができる。このように、ゆらぎを算出する脈拍数を可変できるので、所望のデータ数での脈拍間隔の変動を算出することができる。
脈拍数可変過程では、脈拍数可変部１３が、外部からの入力により、第１脈拍数を増減して新たな第１脈拍数を設定する。外部から脈拍数可変部１３に信号が入力されたときに脈拍数可変過程が行われ、新たな第１脈拍数を設定する。その後、光電脈波検出部１１の検出した脈波１が周波数成分抽出部１２に入力されると、周波数成分抽出部１２が新たな第１脈拍数を用いて脈拍間隔のゆらぎの周波数成分を抽出する。
周波数成分比算出過程では、周波数成分比算出部１５が、周波数成分抽出過程で抽出したゆらぎの周波数成分のうち特定の２つの周波数の周波数成分の比を算出する。特定の２つの周波数の周波数成分の比を算出するので、いずれか一方のみの周波数の周波数成分が著しく減少したことを判定することができる。これにより、脈拍間隔のゆらぎの周波数成分の変動を算出することができる。
また、周波数成分比算出過程において、周波数成分比算出部１５が、周波数成分抽出過程で抽出した周波数成分のうち０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比を算出することが好ましい。０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比の変化は精神的ストレスの判定に役立つとされているので、精神的ストレスの推定に資することができる。
判定過程では、判定部１７が、周波数成分比算出過程において周波数成分比算出部１５の算出する脈拍間隔の周波数成分があらかじめ定められた閾値以上であるか否かを判定する。
表示部１４が、周波数成分抽出過程で抽出したゆらぎの周波数成分をスペクトル表示する表示過程をさらに有することが好ましい。スペクトル表示とは、本実施形態では、周波数ごとの周波数成分を表示するグラフ表示である。表示過程は、算出された脈拍間隔のゆらぎの周波数成分を入力されたときは常時表示することが好ましい。また、外部から信号が入力されたときに表示してもよい。このように、表示過程を有すれば、表示部１４がゆらぎの周波数成分をスペクトル表示することができるので、長時間にわたる脈拍間隔の変動を表示することができる。
以上説明したように、本実施形態に係る脈波診断装置９１及び脈波診断装置９１を用いた脈波診断装置制御方法は、光電脈波検出部１１が光電センサを用いて外耳１０３の脈波１を検出することができるので、長時間にわたり脈波を検出算出することができる。さらに、周波数成分抽出部１２が脈拍間隔のゆらぎの周波数成分を抽出することができるので、脈波の変動を検出することができる。よって、脈波診断装置９１及び脈波診断装置９１を用いた脈波診断装置制御方法は、長時間にわたる脈波の変動を検出することができる。
（実施形態２）
図４は、本実施形態に係る脈波診断装置の構成概略図である。図４に示す脈波診断装置９２は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部１１と、光電脈波検出部１１の検出する脈波の１拍ごとの脈拍間隔を算出し、連続して算出された２つの脈拍間隔同士で形成される直交座標平面上での脈拍間隔の点をポアンカレ座標として算出する脈拍間隔ポアンカレ算出部２１と、を備える。
脈波診断装置９２は、外部からの入力により、第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部１３と、脈拍間隔ポアンカレ算出部２１の算出する脈拍間隔のポアンカレ座標を表示する表示部１４と、脈拍間隔ポアンカレ算出部２１の算出する脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈拍間隔統計処理部２３と、をさらに備えていてもよい。また、所定時刻に光電脈波検出部１１に脈波の検出を始めさせ、所定時刻に光電脈波検出部１１に脈波の検出を終了させるタイマー１６と、脈拍間隔のポアンカレ座標の分布があらかじめ定められた座標範囲内であるか否かを判定する判定部２４と、をさらに備えていてもよい。
光電脈波検出部１１は、動脈を透過した透過光を検出するものである。また、光電脈波検出部１１は、動脈で散乱された散乱光を受光して脈波を検出するものである。例えば、前述の実施形態１と同様のものを用いることができる。光電脈波検出部１１は、人体の外耳１０３への装着機構３１をさらに備えることが好ましい。また、光電脈波検出部１１は、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。
脈拍間隔ポアンカレ算出部２１は、光電脈波検出部１１の検出する脈波の１拍ごとの脈拍間隔を算出するものである。光電脈波検出部１１の検出する脈波の１拍ごとの脈拍間隔を算出する構成は、前述の実施形態１で説明した周波数成分抽出部１２と同様のものを用いることができる。
脈拍間隔ポアンカレ算出部２１は、さらに、連続して算出された２つの脈拍間隔同士で形成される直交座標平面上での脈拍間隔の点をポアンカレ座標として算出するものである。図５は、脈拍間隔の第１例を示すポアンカレプロットである。図６は、脈拍間隔の第２例を示すポアンカレプロットである。図５は、精神的ストレスの印加されていない状態での脈拍間隔の一例を示す。図６は、精神的ストレスの印加された状態での脈拍間隔の一例を示す。図５及び図６において、横軸の脈拍間隔（ｉ）はｉ番目に算出された脈拍間隔、縦軸の脈拍間隔（ｉ＋１）は（ｉ＋１）番目に算出された脈拍間隔を示す。脈拍間隔のポアンカレ座標は、例えばこの座標平面上の座標とすることができる。なお、脈拍間隔のポアンカレ座標は図５、図６に示す座標系に限定されるものではない。例えば、縦軸及び横軸の目盛が異なっていてもよい。
脈拍間隔のポアンカレ座標を算出した場合、精神的ストレスが印加されていなければ、脈拍間隔のポアンカレ座標の分布は図５に示すような楕円形若しくは長円形又はこれらに近い形状となる。しかし、精神的ストレスが印加されると、脈拍間隔のポアンカレ座標の分布は図６に示すような円に近い形状となる。よって、脈拍間隔のポアンカレ座標の分布の形状に基づいて精神的ストレスの有無又はその強度を推定することができる。また、このように、脈拍間隔のポアンカレ座標を算出すれば、医師又は被検者は、精神的ストレスの有無又はその強度を推定することができる。
脈拍間隔ポアンカレ算出部２１は、脈拍間隔のポアンカレ座標を、予め定められた第１脈拍数で算出することが好ましい。さらに、脈拍間隔ポアンカレ算出部２１は、脈拍間隔のポアンカレ座標を、予め定められた第２脈拍数ごとに算出することが好ましい。第１脈拍数及び第２脈拍数は、前述の実施形態１の周波数成分抽出部１２で説明したものと同様とすることができる。脈拍間隔のポアンカレ座標を第２脈拍数ごとに算出することができるので、最新のポアンカレ座標を第２脈拍数ごとに取得することができる。
脈波診断装置９２は、外部からの入力により、第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部１３をさらに備えることが好ましい。脈拍数可変部１３は、前述の実施形態１で説明した脈拍数可変部１３と同様のものを用いることができる。
図４に示す脈波診断装置９２は、脈拍間隔ポアンカレ算出部２１の算出する脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈拍間隔統計処理部２３をさらに備えることが好ましい。脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出することができれば、時刻に依存しない脈拍間隔の分布の統計値を算出することができる。
脈拍間隔統計処理部２３に代えて、脈拍間隔ポアンカレ算出部２１の算出する脈拍間隔のポアンカレ座標の座標範囲内での最大及び最小の変化幅を算出する脈拍間隔分布算出部が配置されていてもよい。図１２は、脈拍間隔の第３例を示すポアンカレプロットである。横軸の脈拍間隔（ｉ）はｉ番目に算出された脈拍間隔、縦軸の脈拍間隔（ｉ＋１）は（ｉ＋１）番目に算出された脈拍間隔を示す。ポアンカレ座標上にＹ＝Ｘを示す直線４５と、Ｙ＝Ｘ＋ａを示す直線４６と、Ｙ＝Ｘ−ａを示す直線４７とが示されている。この場合、脈拍間隔分布算出部は、脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出する。例えば、ポアンカレ座標の最大と最小の差を算出する。さらに、脈拍間隔分布算出部は、予め定められた座標範囲内での最大及び最小の変化幅４８を算出することが好ましい。予め定められた座標範囲内とは、例えば、図中に示す直線４６及び直線４７で囲まれた座標範囲である。このように、最大及び最小の変化幅４８を算出すれば、脈拍間隔のポアンカレ座標の分布が狭まったことを検出できる。さらに、脈拍間隔分布算出部は、脈拍間隔のポアンカレ座標が直線４６及び直線４７で囲まれる座標範囲内にあるか否かを判定してもよい。また、脈拍間隔分布算出部は、直線４６及び直線４７で囲まれる座標範囲に入った脈拍間隔のポアンカレ座標の最大又は最小を算出するにとどまるものでもよい。
ここで、ａは、予め定められた値である。例えば、１０〜２０ｍｓｅｃとすることができる。また、ａは可変であることが好ましい。例えば、ある時間安静を保ったときに検出された脈波の脈拍間隔のポアンカレ座標が分布する座標範囲と、ストレス印加時に検出した脈波の脈拍間隔のポアンカレ座標が分布する座標範囲とを比較し、ストレス印加時では収まるが安静時では収まらないポアンカレ座標の座標範囲からａを定めることが好ましい。
また、Ｘ軸データとＹ軸データの積の最大値から最小値を差し引いてもよい。この場合、前述の図５及び図６に示す脈拍間隔（ｉ）及び脈拍間隔（ｉ＋１）を用いて表すと、（脈拍間隔（ｉ）×脈拍間隔（ｉ＋１））ｍａｘ−（脈拍間隔（ｉ）×脈拍間隔（ｉ＋１））ｍｉｎとなる。このような算出を行うことでも、脈拍間隔のポアンカレ座標の分布が狭まったことを検出できる。このように、ポアンカレ座標の分布範囲を算出すれば、医師又は被検者は、精神的ストレスの有無又はその強度を推定することができる。
図４に示す脈波診断装置９２は、脈拍間隔のポアンカレ座標の分布があらかじめ定められた座標範囲内であるか否かを判定する判定部２４をさらに備えることが好ましい。この場合、あらかじめ定められた座標範囲とは、例えば、ある時間安静を保ったときに検出された脈波の脈拍間隔のポアンカレ座標とすることができる。あらかじめ定められた座標範囲は、前述の図１２で説明した座標範囲と同様の座標範囲を用いることができる。あらかじめ定められた座標範囲は、被検者に応じて可変であることが好ましい。例えば、被検者の脈拍間隔のポアンカレ座標を予め算出し、メモリーに記憶させることができる。判定部２４は、メモリーからあらかじめ定められた座標範囲を取得する。このように、判定部２４が脈拍間隔のポアンカレ座標の分布があらかじめ定められた座標範囲内であるか否かを判定すれば、医師又は被検者は、精神的ストレスの有無又はその強度を推定することができる。
図４に示す脈波診断装置９２は、脈拍間隔ポアンカレ算出部２１の算出する脈拍間隔のポアンカレ座標を表示する表示部１４をさらに備えることが好ましい。表示部１４としては、例えば、前述の実施形態１で説明した表示部１４と同様のものを用いることができる。表示部１４は、さらに、脈拍数可変部１３の設定した新たな第１脈拍数を表示してもよい。又、第２脈拍数を表示してもよい。又、脈拍間隔統計処理部２３の算出する脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を表示してもよい。このように、表示部１４は、脈拍間隔のポアンカレ座標を表示することができるので、時刻に依存しない脈拍間隔の分布を表示することができる。
次に、脈波診断装置９２の動作について説明する。本実施形態に係る脈波診断装置９２を用いた脈波診断装置制御方法は、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、脈拍間隔ポアンカレ算出部２１が、光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を算出し、連続して算出された２つの脈拍間隔同士で形成される直交座標平面上での脈拍間隔の点をポアンカレ座標として算出する脈拍間隔ポアンカレ算出過程と、を有する。
脈波診断装置９２を用いた脈波診断装置制御方法は、表示部１４が、脈拍間隔ポアンカレ算出過程で算出した脈拍間隔のポアンカレ座標を表示するポアンカレプロット表示過程をさらに有していてもよい。又、脈拍間隔統計処理部２３が、脈拍間隔ポアンカレ算出過程で算出した脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈拍間隔統計処理過程をさらに有していてもよい。又、脈拍数可変過程を脈拍間隔ポアンカレ算出過程の前にさらに有していてもよい。又、脈拍間隔ポアンカレ算出過程の後に脈拍間隔統計処理過程をさらに有していてもよい。又、脈拍間隔ポアンカレ算出過程の後に判定過程をさらに有していてもよい。
光電脈波検出過程では、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。光電脈波検出過程は前述の実施形態１で説明した光電脈波検出過程と同様としてもよい。又、光電脈波検出過程において、光電脈波検出部１１が、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。
脈拍間隔ポアンカレ算出過程では、脈拍間隔ポアンカレ算出部２１が、光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を算出する。
さらに、脈拍間隔ポアンカレ算出過程では、連続して算出された２つの脈拍間隔同士で形成される直交座標平面上での脈拍間隔の点をポアンカレ座標として算出する。
脈拍間隔ポアンカレ算出過程において、脈拍間隔ポアンカレ算出部２１が、脈拍間隔のポアンカレ座標を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出することが好ましい。さらに、脈拍間隔ポアンカレ算出過程において、脈拍間隔ポアンカレ算出部２１が、脈拍数可変過程で設定した新たな第１脈拍数を用いて脈拍間隔のポアンカレ座標を算出することが好ましい。第１脈拍数及び第２脈拍数は、前述の実施形態１の周波数成分抽出部１２で説明したものと同様とすることができる。
脈拍数可変過程では、脈拍数可変部１３が、外部からの入力により、第１脈拍数を増減して新たな第１脈拍数を設定する。脈拍間隔のポアンカレ座標を算出する脈拍数を増減することができるので、所望のデータ数での脈拍間隔のポアンカレ座標を算出することができる。
脈拍間隔統計処理過程では、脈拍間隔統計処理部２３が、脈拍間隔ポアンカレ算出過程で算出したポアンカレ座標の標準偏差、中央値又は平均値を算出する。脈拍間隔統計処理過程において、脈拍間隔統計処理部２３は、脈拍間隔ポアンカレ算出過程で算出したポアンカレ座標の標準偏差、中央値及び平均値を表示してもよい。脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出することができるので、時刻に依存しない脈拍間隔の分布の統計値を算出することができる。
脈拍間隔統計処理部２３に脈拍間隔分布算出部が配置される場合、前記脈拍間隔統計処理過程は、脈拍間隔分布算出部が、脈拍間隔ポアンカレ算出部２１の算出する脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出する脈拍間隔分布算出過程となる。この場合、脈拍間隔分布算出過程では、予め定められた時間範囲内の脈拍間隔のポアンカレ座標の最大値と脈拍間隔のポアンカレ座標の最小値との差を算出する。さらに、脈拍間隔分布算出過程において、予め定められた座標範囲内での最大及び最小の変化幅を算出することが好ましい。予め定められた座標範囲内とは、例えば、前述の図１２で説明した直線４６及び直線４７で囲まれた座標範囲内である。このように、予め定められた座標範囲内での最大及び最小の変化幅を算出すれば、脈拍間隔の診断に有用な座標範囲内での変化幅を検出することができる。このように、脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出すれば、脈拍間隔のポアンカレ座標の分布が狭まったことを検出できる。
判定過程では、判定部２４が、脈拍間隔のポアンカレ座標の分布があらかじめ定められた座標範囲内であるか否かを判定する。判定過程を有することで、脈拍間隔のポアンカレ座標の分布があらかじめ定められた座標範囲内であるか否かの判定を自動化することができる。
表示部１４が、脈拍間隔ポアンカレ算出過程で算出した脈拍間隔のポアンカレ座標を表示するポアンカレプロット表示過程をさらに有することが好ましい。表示過程は、算出された脈拍間隔のポアンカレ座標を入力されたときは常時表示することが好ましい。また、外部から信号が入力されたときに脈拍間隔のポアンカレ座標を表示してもよい。表示過程において、表示部１４は、例えば、予め定められた時間範囲内に算出されたポアンカレ座標を表示する。表示する時間範囲は可変であることが好ましい。このように、表示過程において、表示部１４は脈拍間隔のポアンカレ座標を表示することができるので、時刻に依存しない脈拍間隔の分布を表示することができる。
以上説明したように、本実施形態に係る脈波診断装置９２及び脈波診断装置９２を用いた脈波診断装置制御方法は、光電脈波検出部１１が光電センサを用いて外耳１０３の脈波１を検出することができるので、長時間にわたり脈波を検出することができる。さらに、脈拍間隔ポアンカレ算出部２１が脈拍間隔のポアンカレ座標を算出するので、脈波の変動を検出することができる。よって、脈波診断装置９２及び脈波診断装置９２を用いた脈波診断装置制御方法は、長時間にわたる脈波の変動を検出することができる。
（実施形態３）
図７は、本実施形態に係る脈波診断装置の構成概略図である。図７に示す脈波診断装置９３は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部１１と、光電脈波検出部１１の検出する脈波の脈波振幅を算出する脈波振幅算出部２５と、脈波振幅算出部２５の算出する脈波振幅が所定の大きさ以上か否かを判定する判定部２６と、を備える。
脈波診断装置９３は、脈波振幅算出部２５の算出する脈波振幅の時間推移を表示する表示部１４と、脈波振幅算出部２５の算出する平均値のうち特定の２つの平均値の比を算出する脈波振幅比算出部２７と、外部からの入力により、第１脈拍数を増加又は減少して新たな第１脈拍数を設定する脈拍数可変部１３と、をさらに備えていてもよい。また、所定時刻に光電脈波検出部１１に脈波１の検出を始めさせ、所定時刻に光電脈波検出部１１に脈波１の検出を終了させるタイマー１６をさらに備えていてもよい。
光電脈波検出部１１は、動脈を透過した透過光を検出するものである。また、光電脈波検出部１１は、動脈で散乱された散乱光を受光して脈波を検出するものである。例えば、前述の実施形態１と同様のものを用いることができる。光電脈波検出部１１は、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。また、光電脈波検出部１１は、人体の外耳１０３への装着機構３１をさらに備えることが好ましい。さらに、装着機構３１は、光電センサを耳介１０１の付け根に装着することが好ましい。光電センサを耳介１０１の付け根に装着すれば、後耳介動脈及び迷走神経の耳介枝に近接しているため、後耳介動脈の血管状態を検出することができる。
脈波振幅算出部２５は、光電脈波検出部１１の検出する脈波の脈波振幅を算出するものである。脈波振幅は、例えば、前述の図２に示す隣接する山４１と谷４２との信号出力の差異とすることができる。
図８は、光電脈波検出部の検出する脈波振幅の一例を示す時系列グラフである。図８に示す脈波振幅は、点線で示す略１００ｓｅｃの時点以降に精神的ストレスを印加したときの脈波振幅の時間推移を示す。ストレスが印加されていない状態では略１．２Ｖ以下から略１．８Ｖ以上の脈波振幅が算出されているのに対し、精神的ストレスを印加後の略２００ｓｅｃ以降に脈波振幅が略１．２Ｖから略１．８Ｖまでの範囲内しか算出されていない。このように、精神的ストレスの印加を脈波振幅の算出によって推定することができる。これは、ストレスが加わると、「防御反応」により、脳、筋肉の血管は拡張するが、皮膚近傍や内臓の血管は収縮する。光電センサで捕らえる血管は皮膚表面に近いため、血管の収縮による脈波振幅の減少が現れる。心電は心臓を駆動するタイミング信号であるのに対し、脈波振幅は人体の血管の拡張又は収縮の影響を受けるので、脈波振幅から血管の収縮を検出すれば、精神的ストレスが印加されたことを検出することができる。
さらに、血管が体表面に多く出ている指や外耳１０３のうち、指では血管運動（Ｖａｓｏｍｏｔｉｏｎ）が多く、人体からの情報が隠されてしまう。一方の外耳１０３では、血管運動の影響が少ないので、血管からの情報を容易に取得することができる。さらに、脈波振幅の変化は、ストレスの印加に対する応答が速いので、リアルタイムでのストレスの評価に用いることができる。
脈波振幅算出部２５は、光電脈波検出部１１の検出する脈波の脈波振幅の平均値を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出することが好ましい。さらに、脈波振幅算出部２５は、脈拍数可変部１３の設定する新たな第１脈拍数を用いて脈波振幅の平均値を算出することが好ましい。第１脈拍数及び第２脈拍数は、前述の実施形態１の周波数成分抽出部１２で説明したものと同様とすることができる。脈波振幅の平均値を第２脈拍数ごとに算出することができるので、最新の脈波振幅の平均値を第２脈拍数ごとに取得することができる。
脈波診断装置９３は、外部からの入力により、第１脈拍数を増加又は減少して新たな第１脈拍数を設定する脈拍数可変部１３をさらに備えることが好ましい。脈拍数可変部１３は、前述の実施形態１で説明した脈拍数可変部１３と同様のものを用いることができる。
脈波診断装置９３は、脈波振幅算出部２５の算出する平均値のうち特定の２つの平均値の比を算出する脈波振幅比算出部２７をさらに備えることが好ましい。脈波振幅算出部２５の算出する２つの脈波振幅の平均値の比を算出するので、２つの脈波振幅の平均値の増減を判定することができる。
判定部２６は、脈波振幅算出部２５の算出する脈波振幅が所定の大きさ以上か否かを判定するものである。所定の大きさとは、あらかじめ定められた脈波振幅の大きさである。例えば、ストレスが印加されている時の脈波振幅と、ストレスが印加されていないときの脈波振幅との境界の値である。ストレスが印加されていないときの脈波振幅とは、例えば、ある時間安静を保ったときに検出された脈波の脈波振幅である。脈波振幅は個人によって大きさが異なるので、被検者ごとにあらかじめ測定し、閾値を定めることが好ましい。さらに、脈波診断装置９３が脈波振幅比算出部２７をさらに備える場合、判定部２６は脈波振幅算出部２５の算出する２つの脈波振幅の平均値の比が所定の大きさ以上か否かを判定することが好ましい。
脈波振幅算出部２５の算出する脈波振幅の時間推移を表示する表示部１４をさらに備えることが好ましい。又は、脈波振幅算出部２５の算出する脈波振幅の平均値の時間推移を表示する表示部１４をさらに備えることが好ましい。表示部１４としては、例えば、前述の実施形態１で説明したものと同様のものを用いることができる。表示部１４は、脈波振幅算出部２５が脈波振幅の平均値を算出する場合は、脈波振幅の平均値の時間推移を表示することが好ましい。表示部１４は、さらに、脈拍数可変部１３の設定した新たな第１脈拍数を表示してもよい。又、第２脈拍数を表示してもよい。又、判定部２６の判定結果を表示してもよい。又、脈波振幅比算出部２７の算出する特定の２つの平均値の比を表示してもよい。このように、表示部１４は、脈波振幅の時間推移又は脈波振幅の平均値の時間推移を表示することができるので、時刻に依存しない脈拍間隔の分布を表示することができる。
次に、脈波診断装置９３の動作について説明する。本実施形態に係る脈波診断装置９３を用いた脈波診断装置制御方法は、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、脈波振幅算出部２５が、光電脈波検出過程で検出した脈波の脈波振幅を算出する脈波振幅算出過程と、判定部２６が、脈波振幅算出過程で算出した脈波振幅が所定の大きさ以上か否かを判定する判定過程と、を有する。
脈波診断装置９３を用いた脈波診断装置制御方法は、脈拍数可変部１３が、外部からの入力により、第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変過程を脈波振幅算出過程の前にさらに有していてもよい。又、表示部１４が、脈波振幅算出過程で算出された脈波振幅の時間推移か、又は、脈波振幅算出過程で算出された脈波振幅の平均値の時間推移を表示する表示過程をさらに有していてもよい。又、脈波振幅算出過程の後に、判定過程をさらに有していてもよい。
光電脈波検出過程では、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する。光電脈波検出過程は前述の実施形態１で説明した光電脈波検出過程と同様としてもよい。光電脈波検出過程において、光電脈波検出部１１が、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。
脈波振幅算出過程では、脈波振幅算出部２５が、光電脈波検出過程で検出した脈波の脈波振幅を算出する。
判定過程では、判定部２６が、脈波振幅算出過程で算出した脈波振幅が所定の大きさ以上か否かを判定する。所定の大きさは、前述の判定部２６で説明したものと同様である。判定過程は、例えば、脈波振幅が算出されるたびに行うことができる。
脈波振幅算出過程において、脈波振幅算出部２５が、光電脈波検出過程で検出した脈波の脈波振幅の平均値を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出することが好ましい。さらに、脈波振幅算出過程において、脈波振幅算出部２５が、脈拍数可変過程で設定した新たな第１脈拍数を用いて脈波振幅の平均値を算出することが好ましい。第１脈拍数及び第２脈拍数は、前述の実施形態１の周波数成分抽出部１２で説明したものと同様とすることができる。このように、脈波振幅の平均値を第２脈拍数ごとに算出することができるので、最新の脈波振幅の平均値を第２脈拍数ごとに取得することができる。
脈拍数可変過程では、脈拍数可変部１３が、外部からの入力により、第１脈拍数を増減して新たな第１脈拍数を設定する。脈波振幅の平均値を算出する脈拍数を可変できるので、所望のデータ数での脈波振幅の変動を算出することができる。
脈波振幅比算出部２７が、脈波振幅算出過程で算出した平均値のうち特定の２つの平均値の比を算出する脈波振幅比算出過程をさらに有することが好ましい。脈波振幅算出部２５の算出する２つの脈波振幅の平均値の比を算出するので、２つの脈波振幅の平均値の増減を判定することができる。
判定過程では、判定部２６は、脈波振幅算出過程において脈波振幅算出部２５の算出する脈波振幅が所定の大きさ以上か否かを判定する。所定の大きさとは、前述の判定部２６で説明したとおりである。
表示部１４が、脈波振幅算出過程で算出された脈波振幅の時間推移を表示する表示過程をさらに有することが好ましい。又、表示部１４が、脈波振幅算出過程で算出された脈波振幅の平均値の時間推移を表示する表示過程をさらに有することが好ましい。表示過程において表示部１４が表示する脈波振幅の時間推移は、予め定められた時間範囲のものとすることができる。脈波振幅の平均値の時間推移についても同様である。このように、表示過程において、表示部１４は脈波振幅の時間推移又は脈波振幅の平均値の時間推移を表示することができる。
以上説明したように、本実施形態に係る脈波診断装置９３及び脈波診断装置９３を用いた脈波診断装置制御方法は、光電脈波検出部１１が光電センサを用いて外耳１０３の脈波１を検出することができるので、長時間にわたり脈波を検出算出することができる。さらに、脈波振幅算出部２５が脈波の脈波振幅を算出するので、脈波の変動を検出することができる。よって、脈波診断装置９３及び脈波診断装置９３を用いた脈波診断装置制御方法は、長時間にわたる脈波の変動を検出することができる。さらに、光電センサを用いて抹消血管の脈波を検出することができるので、血管運動が大きい場合でも、脈波振幅を算出することによって血管状態の変化を検出することができる。
（実施形態４）
図９は、本実施形態に係る脈波診断装置の構成概略図である。図９に示す脈波診断装置９４は、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部１１と、光電脈波検出部１１の検出する脈波の１拍ごとの脈波振幅を算出し、連続して算出された２つの脈波振幅同士で形成される直交座標平面上での脈波振幅の点をポアンカレ座標として１拍ごとに算出する脈波振幅ポアンカレ算出部２８と、を備えることを特徴とする。
脈波診断装置９４は、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標を表示する表示部１４と、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理部２９と、外部からの入力により、第１脈拍数を増加又は減少して新たな第１脈拍数を設定する脈拍数可変部１３と、をさらに備えていてもよい。また、所定時刻に光電脈波検出部１１に脈波１の検出を始めさせ、所定時刻に光電脈波検出部１１に脈波１の検出を終了させるタイマー１６と、脈波振幅ポアンカレ算出部２８の算出する周波数成分の比があらかじめ定められた閾値以上であるか否かを判定する判定部３０と、をさらに備えていてもよい。
光電脈波検出部１１は、動脈を透過した透過光を検出するものである。また、光電脈波検出部１１は、動脈で散乱された散乱光を受光して脈波を検出するものである。例えば、前述の実施形態１と同様のものを用いることができる。光電脈波検出部１１は、人体の外耳１０３への装着機構３１をさらに備えることが好ましい。また、光電脈波検出部１１は、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。
脈波振幅ポアンカレ算出部２８は、光電脈波検出部１１の検出する脈波１の１拍ごとの脈波振幅を算出するものである。脈波振幅は、例えば、前述の図２に示す隣接する山４１と谷４２との信号出力の差異とすることができる。
脈波振幅ポアンカレ算出部２８は、連続して算出された２つの脈波振幅同士で形成される直交座標平面上での脈波振幅の点をポアンカレ座標として１拍ごとに算出するものである。図１０は、脈波振幅の第１例を示すポアンカレプロットである。図１１は、脈波振幅の第２例を示すポアンカレプロットである。図１０は、精神的ストレスの印加されていない状態での脈波振幅の一例を示す。図１１は、精神的ストレスの印加された状態での脈波振幅の一例を示す。図１０及び図１１において、横軸はｉ番目に算出された脈波振幅、縦軸は（ｉ＋１）番目に算出された脈波振幅を示す。脈波振幅のポアンカレ座標は、例えばこの座標平面上の座標とすることができる。なお、脈波振幅のポアンカレ座標は図１０、図１１に示す座標系に限定されるものではない。例えば、縦軸及び横軸の目盛が異なっていてもよい。
脈波振幅のポアンカレ座標を算出した場合、精神的ストレスが印加されていなければ、脈波振幅のポアンカレ座標の分布は楕円形若しくは長円形又はこれらに近い形状となる。しかし、精神的ストレスが印加されると、脈波振幅のポアンカレ座標の分布は円に近くなる。よって、脈波振幅のポアンカレ座標の分布の形状に基づいて精神的ストレスの有無又はその強度を推定することができる。このように、脈波振幅のポアンカレ座標を算出すれば、医師又は被検者は、精神的ストレスの有無又はその強度を推定することができる。
脈波振幅ポアンカレ算出部２８は、脈波振幅のポアンカレ座標を予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出することが好ましい。さらに、脈波振幅ポアンカレ算出部２８は、脈波振幅のポアンカレ座標を、予め定められた第２脈拍数ごとに算出することが好ましい。第１脈拍数及び第２脈拍数は、前述の実施形態１の周波数成分抽出部１２で説明したものと同様とすることができる。脈波振幅のポアンカレ座標を第２脈拍数ごとに算出することができるので、脈波振幅のポアンカレ座標を第２脈拍数ごとに取得することができる。
脈波診断装置９４は、外部からの入力により、第１脈拍数を増加又は減少して新たな第１脈拍数を設定する脈拍数可変部１３をさらに備えることが好ましい。脈拍数可変部１３は、前述の実施形態１で説明した脈拍数可変部１３と同様のものを用いることができる。
脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理部２９をさらに備えることが好ましい。脈波振幅統計処理部２９が脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出するので、時刻に依存しない脈波振幅の分布の統計値を算出することができる。
脈波振幅統計処理部２９に代えて、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出する脈波振幅分布算出部が配置されていてもよい。図１３は、脈波振幅の第３例を示すポアンカレプロットである。横軸はｉ番目に算出された脈波振幅、縦軸は（ｉ＋１）番目に算出された脈波振幅を示す。ポアンカレ座標上にＹ＝Ｘを示す直線５５と、Ｙ＝Ｘ＋ａを示す直線５６と、Ｙ＝Ｘ−ａを示す直線５７とが示されている。この場合、脈波振幅分布算出部は、脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出する。例えば、ポアンカレ座標の最大と最小の差を算出する。さらに、脈波振幅分布算出部は、予め定められた座標範囲内での最大及び最小の変化幅５８を算出することが好ましい。予め定められた座標範囲内とは、例えば、図中に示す直線５６及び直線５７で囲まれた座標範囲である。このように、最大及び最小の変化幅５８を算出すれば、脈波振幅のポアンカレ座標の分布が狭まったことを検出できる。さらに、脈波振幅分布算出部は、脈波振幅のポアンカレ座標が直線５６及び直線５７で囲まれる座標範囲内にあるか否かを判定してもよい。また、脈波振幅分布算出部は、直線５６及び直線５７で囲まれる座標範囲に入った脈波振幅のポアンカレ座標の最大又は最小を算出するにとどまるものでもよい。
ここで、ａは、予め定められた値である。ａは可変であることが好ましい。例えば、ある時間安静を保ったときに検出された脈波の脈波振幅のポアンカレ座標が分布する座標範囲と、ストレス印加時に検出した脈波の脈波振幅のポアンカレ座標が分布する座標範囲とを比較し、ストレス印加時では収まるが安静時では収まらない座標範囲をａとすることが好ましい。
また、脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出してもよい。また、Ｘ軸データとＹ軸データの積の最大値から最小値を差し引いてもよい。この場合、前述の図１０及び図１１に示す脈波振幅（ｉ）及び脈波振幅（ｉ＋１）を用いてあらわすと、（脈波振幅（ｉ）×脈波振幅（ｉ＋１））ｍａｘ−（脈波振幅（ｉ）×脈波振幅（ｉ＋１））ｍｉｎとなる。このような算出を行うことで、脈波振幅のポアンカレ座標の分布が狭まったことを検出できる。これにより、医師又は被検者は、精神的ストレスの有無又はその強度を推定することができる。
脈波診断装置９４は、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標があらかじめ定められた座標範囲内であるか否かを判定する判定部３０をさらに備えることが好ましい。この場合、あらかじめ定められた座標範囲とは、例えば、ある時間安静を保ったときに検出された脈波の脈波振幅のポアンカレ座標とすることができる。あらかじめ定められた座標範囲は、前述の図１３で説明した座標範囲と同様の座標範囲を用いることができる。あらかじめ定められた座標範囲は、被検者に応じて可変であることが好ましい。例えば、被検者の脈波振幅のポアンカレ座標を予め算出し、メモリーに記憶させることができる。判定部３０は、メモリーからあらかじめ定められた座標範囲を取得する。このように、判定部３０が、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標の分布があらかじめ定められた座標範囲内であるか否かを判定すれば、医師又は被検者は、精神的ストレスの有無又はその強度を推定することができる。
脈波診断装置９４は、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標を表示する表示部１４をさらに備えることが好ましい。表示部１４は、さらに、脈波振幅統計処理部２９の算出する脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を表示してもよい。又、判定部３０の判定結果を表示してもよい。このように、表示部１４は、脈波振幅のポアンカレ座標を表示することができるので、時刻に依存しない脈波振幅の分布を表示することができる。
次に、脈波診断装置９４の動作について説明する。本実施形態に係る脈波診断装置９４を用いた脈波診断装置制御方法は、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、脈波振幅ポアンカレ算出部２８が、光電脈波検出過程で検出した脈波の１拍ごとの脈波振幅を算出し、連続して算出された２つの脈波振幅同士で形成される直交座標平面上での脈波振幅の点をポアンカレ座標として１拍ごとに算出する脈波振幅ポアンカレ算出過程と、を有する。
脈波診断装置９４を用いた脈波診断装置制御方法は、脈拍数可変部１３が、外部からの入力により第１脈拍数を新たな第１脈拍数として設定する脈拍数可変過程を脈波振幅ポアンカレ算出過程の前にさらに有していてもよい。又、表示部１４が、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標を表示する表示過程をさらに有していてもよい。又、脈波振幅統計処理部２９が、脈波振幅ポアンカレ算出過程で算出した脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理過程をさらに有していてもよい。又、脈波振幅ポアンカレ算出過程の後に判定過程をさらに有していてもよい。
光電脈波検出過程では、光電脈波検出部１１が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波１を検出する。光電脈波検出過程は前述の実施形態１で説明した光電脈波検出過程と同様としてもよい。光電脈波検出過程において、光電脈波検出部１１が、検出する脈波１を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することが好ましい。
脈波振幅ポアンカレ算出過程では、脈波振幅ポアンカレ算出部２８が、光電脈波検出過程で検出した脈波１の１拍ごとの脈波振幅を算出し、連続して算出された２つの脈波振幅同士で形成される直交座標平面上での脈波振幅の点をポアンカレ座標として１拍ごとに算出する。
脈波振幅ポアンカレ算出過程において、脈波振幅ポアンカレ算出部２８が、脈波振幅のポアンカレ座標を予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出することが好ましい。ここで、第１脈拍数は、脈拍数可変過程で設定する新たな第１脈拍数である。第１脈拍数及び第２脈拍数は、前述の実施形態１の周波数成分抽出部１２で説明したものと同様とすることができる。脈波振幅のポアンカレ座標を第２脈拍数ごとに算出することができるので、最新の脈波振幅のポアンカレ座標を第２脈拍数ごとに取得することができる。
脈拍数可変過程では、脈拍数可変部１３が、外部からの入力により第１脈拍数を新たな第１脈拍数として設定する。脈波振幅のポアンカレ座標を算出する脈拍数を増減することができるので、所望のデータ数での脈波振幅のポアンカレ座標を算出することができる。
脈波振幅統計処理部２９が、脈波振幅ポアンカレ算出過程で算出した脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理過程をさらに有することが好ましい。脈波振幅統計処理過程では、脈波振幅統計処理部２９が脈波振幅のポアンカレ座標の標準偏差、中央値及び平均値を算出してもよい。このように、脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出するので、時刻に依存しない脈波振幅の分布の統計値を算出することができる。
脈波振幅統計処理部２９に脈波振幅分布算出部が配置される場合、前記脈波振幅統計処理過程は、脈波振幅分布算出部が、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出する脈波振幅分布算出過程となる。この場合、脈波振幅分布算出過程では、予め定められた時間範囲内の脈波振幅のポアンカレ座標の最大値と脈波振幅のポアンカレ座標の最小値との差を算出する。さらに、脈波振幅分布算出過程において、予め定められた座標範囲内での最大及び最小の変化幅を算出することが好ましい。予め定められた座標範囲内とは、例えば、前述の図１３で説明した直線５６及び直線５７である。このように、予め定められた座標範囲内での最大及び最小の変化幅を算出すれば、脈波振幅の診断に有用な座標範囲内での変化幅を検出することができる。このように、脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出すれば、脈波振幅のポアンカレ座標の分布が狭まったことを検出できる。
判定過程では、判定部３０が、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標があらかじめ定められた座標範囲内であるか否かを判定する。あらかじめ定められた座標範囲とは、前述の図９の判定部３０で説明したとおりである。判定過程では、脈波振幅分布算出過程で算出した最大及び最小の変化幅があらかじめ定められた座標範囲内であるか否かを判定してもよい。
表示部１４が、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標を表示する表示過程をさらに有することが好ましい。表示過程は、算出された脈波振幅のポアンカレ座標を入力されたときは常時表示することが好ましい。また、外部から信号が入力されたときに脈波振幅のポアンカレ座標を表示してもよい。表示過程において、表示部１４は、例えば、予め定められた時間範囲内に算出されたポアンカレ座標を表示する。表示する時間範囲は可変であることが好ましい。なお、表示部１４が、脈波振幅ポアンカレ算出部２８の算出する脈波振幅のポアンカレ座標の平均値の時間推移を表示してもよい。このように、表示過程において、表示部１４は脈波振幅のポアンカレ座標の平均値の時間推移を表示することができる。
以上説明したように、本実施形態に係る脈波診断装置９４及び脈波診断装置９４を用いた脈波診断装置制御方法は、光電脈波検出部１１が光電センサを用いて外耳１０３の脈波１を検出することができるので、長時間にわたり脈波を検出算出することができる。さらに、脈波振幅ポアンカレ算出部２８が脈波振幅のポアンカレ座標を算出するので、脈波１の変動を検出することができる。よって、脈波診断装置９４及び脈波診断装置９４を用いた脈波診断装置制御方法は、長時間にわたる脈波の変動を検出することができる。
本発明は、脈波を検出できるので、美容や健康の用途にも用いることができる。
実施形態１に係る脈波診断装置の構成概略図である。 光電脈波検出部の検出する脈波の一例を示す模式図である。 脈拍間隔のゆらぎの周波数成分の一例を示すグラフである。 実施形態２に係る脈波診断装置の構成概略図である。 精神的ストレスの印加されていない状態での脈拍間隔の一例を示すポアンカレプロットである。 精神的ストレスの印加された状態での脈拍間隔の一例を示すポアンカレプロットである。 実施形態３に係る脈波診断装置の構成概略図である。 光電脈波検出部の検出する脈波振幅の一例を示す時系列グラフである。 実施形態４に係る脈波診断装置の構成概略図である。 精神的ストレスの印加されていない状態での脈波振幅の一例を示すポアンカレプロットである。 精神的ストレスの印加された状態での脈波振幅の一例を示すポアンカレプロットである。 脈拍間隔の第３例を示すポアンカレプロットである。 脈波振幅の第３例を示すポアンカレプロットである。
符号の説明
１ 脈波
１１ 光電脈波検出部
１２ 周波数成分抽出部
１３ 脈拍数可変部
１４ 表示部
１５ 周波数成分比算出部
１６ タイマー
１７、２４、２６、３０ 判定部
２１ 脈拍間隔ポアンカレ算出部
２３ 脈拍間隔統計処理部
２５ 脈波振幅算出部
２７ 脈波振幅比算出部
２８ 脈波振幅ポアンカレ算出部
２９ 脈波振幅統計処理部
３１ 装着機構
４１ 山
４２ 谷
４３ 脈拍間隔
４５ 脈拍間隔のポアンカレ座標平面状のＹ＝Ｘを示す直線
４６ Ｙ＝Ｘ＋ａを示す直線
４７ Ｙ＝Ｘ−ａを示す直線
４８ 脈拍間隔のポアンカレ座標の最大及び最小の変化幅
５５ 脈波振幅のポアンカレ座標平面上のＹ＝Ｘを示す直線
５６ 脈波振幅のポアンカレ座標平面上のＹ＝Ｘ＋ａを示す直線
５７ 脈波振幅のポアンカレ座標平面上のＹ＝Ｘ−ａを示す直線
５８ 脈波振幅のポアンカレ座標の最大及び最小の変化幅
９１、９２、９３、９４ 脈波診断装置
１０１ 耳介
１０２ 外耳道
１０３ 外耳
Ｆ１、Ｆ２ 周波数

Claims (45)

  1. 動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、
    前記光電脈波検出部の検出した脈波の１拍ごとの脈拍間隔を算出し、前記脈拍間隔のゆらぎの周波数成分を抽出する周波数成分抽出部と、を備えることを特徴とする脈波診断装置。
  2. 前記周波数成分抽出部は、前記脈拍間隔のゆらぎの周波数成分を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに抽出し、
    外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部をさらに備え、
    前記周波数成分抽出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈拍間隔のゆらぎの周波数成分を抽出することを特徴とする請求項１に記載の脈波診断装置。
  3. 前記周波数成分抽出部の抽出するゆらぎの周波数成分をスペクトル表示する表示部をさらに備えることを特徴とする請求項１又は２に記載の脈波診断装置。
  4. 前記周波数成分抽出部の抽出するゆらぎの周波数成分のうち特定の２つの周波数の周波数成分の比を算出する周波数成分比算出部をさらに備えることを特徴とする請求項１から３に記載のいずれかの脈波診断装置。
  5. 前記周波数成分比算出部は、前記周波数成分抽出部の抽出するゆらぎの周波数成分のうち０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比を算出することを特徴とする請求項４に記載の脈波診断装置。
  6. 動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、
    前記光電脈波検出部の検出する脈波の１拍ごとの脈拍間隔を算出し、連続して算出された２つの前記脈拍間隔同士で形成される直交座標平面上での前記脈拍間隔の点をポアンカレ座標として算出する脈拍間隔ポアンカレ算出部と、を備えることを特徴とする脈波診断装置。
  7. 前記脈拍間隔ポアンカレ算出部は、前記脈拍間隔のポアンカレ座標を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、
    外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部をさらに備え、
    前記脈拍間隔ポアンカレ算出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈拍間隔のポアンカレ座標を算出することを特徴とする請求項６に記載の脈波診断装置。
  8. 前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標を表示する表示部をさらに備えることを特徴とする請求項６又は７に記載の脈波診断装置。
  9. 前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈拍間隔統計処理部をさらに備えることを特徴とする請求項６から８に記載のいずれかの脈波診断装置。
  10. 前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出する脈拍間隔分布算出部をさらに備えることを特徴とする請求項６から８に記載のいずれかの脈波診断装置。
  11. 前記脈拍間隔分布算出部は、脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することを特徴とする請求項１０に記載の脈波診断装置。
  12. 動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、
    前記光電脈波検出部の検出する脈波の脈波振幅を算出する脈波振幅算出部と、
    前記脈波振幅算出部の算出する脈波振幅が所定の大きさ以上か否かを判定する判定部と、を備えることを特徴とする脈波診断装置。
  13. 前記脈波振幅算出部は、前記光電脈波検出部の検出する脈波の脈波振幅の平均値を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、
    外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変部をさらに備え、
    前記脈波振幅算出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈波振幅の平均値を算出することを特徴とする請求項１２に記載の脈波診断装置。
  14. 前記脈波振幅算出部の算出する脈波振幅の時間推移か、又は、前記脈波振幅算出部の算出する脈波振幅の平均値の時間推移を表示する表示部をさらに備えることを特徴とする請求項１２又は１３に記載の脈波診断装置。
  15. 前記脈波振幅算出部の算出する平均値のうち特定の２つの平均値の比を算出する脈波振幅比算出部をさらに備えることを特徴とする請求項１２から１４に記載のいずれかの脈波診断装置。
  16. 動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出部と、
    前記光電脈波検出部の検出する脈波の１拍ごとの脈波振幅を算出し、連続して算出された２つの前記脈波振幅同士で形成される直交座標平面上での前記脈波振幅の点をポアンカレ座標として１拍ごとに算出する脈波振幅ポアンカレ算出部と、を備えることを特徴とする脈波診断装置。
  17. 前記脈波振幅ポアンカレ算出部は、脈波振幅のポアンカレ座標を予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、
    外部からの入力により前記第１脈拍数を新たな第１脈拍数として設定する脈拍数可変部をさらに備え、
    前記脈波振幅ポアンカレ算出部は、前記脈拍数可変部の設定する新たな第１脈拍数を用いて脈波振幅のポアンカレ座標を算出することを特徴とする請求項１６に記載の脈波診断装置。
  18. 前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標を表示する表示部をさらに備えることを特徴とする請求項１６又は１７に記載の脈波診断装置。
  19. 前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理部をさらに備えることを特徴とする請求項１６から１８に記載のいずれかの脈波診断装置。
  20. 前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出する脈波振幅分布算出部をさらに備えることを特徴とする請求項１６から１８に記載のいずれかの脈波診断装置。
  21. 前記脈波振幅分布算出部は、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することを特徴とする請求項２０に記載の脈波診断装置。
  22. 前記光電脈波検出部は、人体の外耳への装着機構をさらに備えることを特徴とする請求項１から２１に記載のいずれかの脈波診断装置。
  23. 前記光電脈波検出部は、検出する脈波を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することを特徴とする請求項１から２２に記載のいずれかの脈波診断装置。
  24. 光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、
    周波数成分抽出部が、前記光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を抽出し、前記脈拍間隔のゆらぎの周波数成分を抽出する周波数成分抽出過程と、を有することを特徴とする脈波診断装置制御方法。
  25. 前記周波数成分抽出過程において、前記周波数成分抽出部が、前記脈拍間隔のゆらぎの周波数成分を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに抽出し、
    脈拍数可変部が、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変過程を前記周波数成分抽出過程の前にさらに有し、
    前記周波数成分抽出過程において、前記周波数成分抽出部が、前記脈拍数可変過程で設定した新たな第１脈拍数を用いて脈拍間隔のゆらぎの周波数成分を抽出することを特徴とする請求項２４に記載の脈波診断装置制御方法。
  26. 表示部が、前記周波数成分抽出過程で抽出したゆらぎの周波数成分をスペクトル表示する表示過程をさらに有することを特徴とする請求項２４又は２５に記載の脈波診断装置制御方法。
  27. 周波数成分比算出部が、前記周波数成分抽出過程で抽出したゆらぎの周波数成分のうち特定の２つの周波数の周波数成分の比を算出する周波数成分比算出過程をさらに有することを特徴とする請求項２４から２６に記載のいずれかの脈波診断装置制御方法。
  28. 前記周波数成分比算出過程において、前記周波数成分比算出部が、前記周波数成分抽出過程で抽出した周波数成分のうち０．１Ｈｚ付近と０．３Ｈｚ付近との周波数成分の比を算出することを特徴とする請求項２７に記載の脈波診断装置制御方法。
  29. 光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、
    脈拍間隔ポアンカレ算出部が、前記光電脈波検出過程で検出した脈波の１拍ごとの脈拍間隔を算出し、連続して算出された２つの前記脈拍間隔同士で形成される直交座標平面上での前記脈拍間隔の点をポアンカレ座標として算出する脈拍間隔ポアンカレ算出過程と、を有することを特徴とする脈波診断装置制御方法。
  30. 前記脈拍間隔ポアンカレ算出過程において、前記脈拍間隔ポアンカレ算出部が、前記脈拍間隔のポアンカレ座標を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、
    脈拍数可変部が、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変過程を前記脈拍間隔ポアンカレ算出過程の前にさらに有し、
    前記脈拍間隔ポアンカレ算出過程において、前記脈拍間隔ポアンカレ算出部が、前記脈拍数可変過程で設定した新たな第１脈拍数を用いて脈拍間隔のポアンカレ座標を算出することを特徴とする請求項２９に記載の脈波診断装置制御方法。
  31. 表示部が、前記脈拍間隔ポアンカレ算出過程で算出した脈拍間隔のポアンカレ座標を表示するポアンカレプロット表示過程をさらに有することを特徴とする請求項２９又は３０に記載の脈波診断装置制御方法。
  32. 脈拍間隔統計処理部が、前記脈拍間隔ポアンカレ算出過程で算出した脈拍間隔のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈拍間隔統計処理過程をさらに有することを特徴とする請求項２９から３１に記載のいずれかの脈波診断装置制御方法。
  33. 前記脈拍間隔分布算出部が、前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の最大及び最小の変化幅を算出する脈拍間隔分布算出過程をさらに有することを特徴とする請求項２９から３１に記載のいずれかの脈波診断装置制御方法。
  34. 前記脈拍間隔分布算出過程において、前記脈拍間隔分布算出部が、前記脈拍間隔ポアンカレ算出部の算出する脈拍間隔のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することを特徴とする請求項３３に記載の脈波診断装置制御方法。
  35. 光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、
    脈波振幅算出部が、前記光電脈波検出過程で検出した脈波の脈波振幅を算出する脈波振幅算出過程と、
    判定部が、前記脈波振幅算出過程で算出した脈波振幅が所定の大きさ以上か否かを判定する判定過程と、を有することを特徴とする脈波診断装置制御方法。
  36. 前記脈波振幅算出過程において、前記脈波振幅算出部が、前記光電脈波検出過程で検出した脈波の脈波振幅の平均値を、予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、
    脈拍数可変部が、外部からの入力により、前記第１脈拍数を増減して新たな第１脈拍数を設定する脈拍数可変過程を前記脈波振幅算出過程の前にさらに有し、
    前記脈波振幅算出過程において、前記脈波振幅算出部が、前記脈拍数可変過程で設定した新たな第１脈拍数を用いて脈波振幅の平均値を算出することを特徴とする請求項３５に記載の脈波診断装置制御方法。
  37. 表示部が、前記脈波振幅算出過程で算出された脈波振幅の時間推移か、又は、前記脈波振幅算出過程で算出された脈波振幅の平均値の時間推移を表示する表示過程をさらに有することを特徴とする請求項３５又は３６に記載の脈波診断装置制御方法。
  38. 脈波振幅比算出部が、前記脈波振幅算出過程で算出した平均値のうち特定の２つの平均値の比を算出する脈波振幅比算出過程をさらに有することを特徴とする請求項３５から３７に記載のいずれかの脈波診断装置制御方法。
  39. 光電脈波検出部が、動脈を透過した透過光又は動脈で散乱された散乱光を受光して脈波を検出する光電脈波検出過程と、
    脈波振幅ポアンカレ算出部が、前記光電脈波検出過程で検出した脈波の１拍ごとの脈波振幅を算出し、連続して算出された２つの前記脈波振幅同士で形成される直交座標平面上での前記脈波振幅の点をポアンカレ座標として１拍ごとに算出する脈波振幅ポアンカレ算出過程と、を有することを特徴とする脈波診断装置制御方法。
  40. 前記脈波振幅ポアンカレ算出過程において、前記脈波振幅ポアンカレ算出部が、前記ポアンカレ座標を予め定められた第１脈拍数で予め定められた第２脈拍数ごとに算出し、
    脈拍数可変部が、外部からの入力により前記第１脈拍数を新たな第１脈拍数として設定する脈拍数可変過程を前記脈波振幅ポアンカレ算出過程の前にさらに有し、
    前記脈波振幅ポアンカレ算出過程において、前記脈波振幅ポアンカレ算出部が、前記脈拍数可変過程で設定する新たな第１脈拍数を用いて脈波振幅のポアンカレ座標を算出することを特徴とする請求項３９に記載の脈波診断装置制御方法。
  41. 表示部が、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標を表示する表示過程をさらに有することを特徴とする請求項３９又は４０に記載の脈波診断装置制御方法。
  42. 脈波振幅統計処理部が、前記脈波振幅ポアンカレ算出過程で算出した脈波振幅のポアンカレ座標の標準偏差、中央値又は平均値を算出する脈波振幅統計処理過程をさらに有することを特徴とする請求項３９から４１に記載のいずれかの脈波診断装置制御方法。
  43. 前記脈波振幅分布算出部が、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の最大及び最小の変化幅を算出する脈波振幅分布算出過程をさらに有することを特徴とする請求項３９から４１に記載のいずれかの脈波診断装置制御方法。
  44. 前記脈波振幅分布算出過程において、前記脈波振幅分布算出部が、前記脈波振幅ポアンカレ算出部の算出する脈波振幅のポアンカレ座標の予め定められた座標範囲内での最大及び最小の変化幅を算出することを特徴とする請求項４３に記載の脈波診断装置制御方法。
  45. 前記光電脈波検出過程において、前記光電脈波検出部が、検出する脈波を５００μｓｅｃ以上５ｍｓｅｃ以下のサンプリング間隔で抽出することを特徴とする請求項２４から４４に記載のいずれかの脈波診断装置制御方法。

JP2005297097A 2005-10-12 2005-10-12 脈波診断装置及び脈波診断装置制御方法 Pending JP2007105131A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005297097A JP2007105131A (ja) 2005-10-12 2005-10-12 脈波診断装置及び脈波診断装置制御方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005297097A JP2007105131A (ja) 2005-10-12 2005-10-12 脈波診断装置及び脈波診断装置制御方法

Related Child Applications (2)

Application Number Title Priority Date Filing Date
JP2010094405A Division JP2010188152A (ja) 2010-04-15 2010-04-15 脈波診断装置及び脈波診断装置制御方法
JP2010094406A Division JP5160586B2 (ja) 2010-04-15 2010-04-15 脈波診断装置及び脈波診断装置制御方法

Publications (1)

Publication Number Publication Date
JP2007105131A true JP2007105131A (ja) 2007-04-26

Family

ID=38031469

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005297097A Pending JP2007105131A (ja) 2005-10-12 2005-10-12 脈波診断装置及び脈波診断装置制御方法

Country Status (1)

Country Link
JP (1) JP2007105131A (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011519653A (ja) * 2008-05-08 2011-07-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 生理学的状態を決定する方法及びシステム
WO2014184868A1 (ja) * 2013-05-14 2014-11-20 株式会社 東芝 電子機器および生体信号測定方法
JP2018511378A (ja) * 2015-03-11 2018-04-26 プレコルディール オサケユイチア 心不全を示す情報を生成するための方法及び装置

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05329116A (ja) * 1992-05-29 1993-12-14 Matsushita Electric Works Ltd 脈拍計
JPH05329118A (ja) * 1992-05-31 1993-12-14 Shimadzu Corp 生体データ収集装置
JPH06142073A (ja) * 1992-10-30 1994-05-24 Fukuda Denshi Co Ltd 加算平均時パワースペクトル表示方法
JPH07143972A (ja) * 1993-11-25 1995-06-06 Matsushita Electric Works Ltd 睡眠状態判定方法及びその装置
JPH09308770A (ja) * 1996-03-18 1997-12-02 Ikyo Kk 心身状態を反映させたゲーム機及び心身状態の判断のためのデータを表示する方法及び心身状態を判定して表明する方法及び心身状態に基づいて制御対象を制御する方法
JP2000333919A (ja) * 1999-05-25 2000-12-05 Nec Corp 生体情報計測装置
JP2001212095A (ja) * 2000-02-07 2001-08-07 Terumo Corp リアルタイム拍動モニター
JP2001212089A (ja) * 2000-02-07 2001-08-07 Terumo Corp 心拍間隔表示方法および心拍間隔表示装置
JP2002191569A (ja) * 2000-12-26 2002-07-09 Seiko Precision Inc 脈拍計及び脈拍数の測定方法
JP2003290164A (ja) * 2002-03-29 2003-10-14 Denso Corp 睡眠支援装置
JP2004016248A (ja) * 2002-06-12 2004-01-22 Junichiro Hayano 心電図に基づく発作性心房細動の検出・評価システム
JP2004528104A (ja) * 2001-04-30 2004-09-16 センソメトリクス エイエス 体腔内の圧力監視のための装置と方法とシステム
WO2005000119A1 (ja) * 2003-06-27 2005-01-06 Soiken Inc. 疲労度評価方法、疲労度評価装置、およびデータベース

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05329116A (ja) * 1992-05-29 1993-12-14 Matsushita Electric Works Ltd 脈拍計
JPH05329118A (ja) * 1992-05-31 1993-12-14 Shimadzu Corp 生体データ収集装置
JPH06142073A (ja) * 1992-10-30 1994-05-24 Fukuda Denshi Co Ltd 加算平均時パワースペクトル表示方法
JPH07143972A (ja) * 1993-11-25 1995-06-06 Matsushita Electric Works Ltd 睡眠状態判定方法及びその装置
JPH09308770A (ja) * 1996-03-18 1997-12-02 Ikyo Kk 心身状態を反映させたゲーム機及び心身状態の判断のためのデータを表示する方法及び心身状態を判定して表明する方法及び心身状態に基づいて制御対象を制御する方法
JP2000333919A (ja) * 1999-05-25 2000-12-05 Nec Corp 生体情報計測装置
JP2001212095A (ja) * 2000-02-07 2001-08-07 Terumo Corp リアルタイム拍動モニター
JP2001212089A (ja) * 2000-02-07 2001-08-07 Terumo Corp 心拍間隔表示方法および心拍間隔表示装置
JP2002191569A (ja) * 2000-12-26 2002-07-09 Seiko Precision Inc 脈拍計及び脈拍数の測定方法
JP2004528104A (ja) * 2001-04-30 2004-09-16 センソメトリクス エイエス 体腔内の圧力監視のための装置と方法とシステム
JP2003290164A (ja) * 2002-03-29 2003-10-14 Denso Corp 睡眠支援装置
JP2004016248A (ja) * 2002-06-12 2004-01-22 Junichiro Hayano 心電図に基づく発作性心房細動の検出・評価システム
WO2005000119A1 (ja) * 2003-06-27 2005-01-06 Soiken Inc. 疲労度評価方法、疲労度評価装置、およびデータベース

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011519653A (ja) * 2008-05-08 2011-07-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 生理学的状態を決定する方法及びシステム
JP2015027532A (ja) * 2008-05-08 2015-02-12 コーニンクレッカ フィリップス エヌ ヴェ 人の生理学的状態を決定するシステムの作動方法
WO2014184868A1 (ja) * 2013-05-14 2014-11-20 株式会社 東芝 電子機器および生体信号測定方法
JP2018511378A (ja) * 2015-03-11 2018-04-26 プレコルディール オサケユイチア 心不全を示す情報を生成するための方法及び装置

Similar Documents

Publication Publication Date Title
KR102622403B1 (ko) 생물학적 데이터 처리
JP7191159B2 (ja) コンピュータプログラム、及び、被験者の感情状態を提供する方法
US9943237B2 (en) Analysis of direct and indirect heartbeat data variations
JP4855721B2 (ja) 血圧測定装置
JP4229919B2 (ja) 脈波検出装置及びその方法
JP5979604B2 (ja) 生体情報検出装置及び生体情報検出方法、生体情報検出プログラム
US20150073239A1 (en) Continuous cuffless blood pressure measurement using a mobile device
US20160262690A1 (en) Method for managing sleep quality and apparatus utilizing the same
EP3878352A1 (en) Apparatus and method for estimating bio-information
KR20190029889A (ko) 생체정보 추정 장치 및 방법
US20210085259A1 (en) Apparatus and method for estimating bio-information
US11844631B2 (en) Apparatus and method for estimating bio-information
JP6522327B2 (ja) 脈波解析装置
JP2001198094A (ja) 脈拍数検出装置
WO2016006250A1 (ja) 生体情報測定装置
JP6547838B2 (ja) 診断支援装置、診断支援方法、診断支援プログラム
KR20200050530A (ko) 혈관 건강도 측정 시스템 및 그 방법
JP2007105131A (ja) 脈波診断装置及び脈波診断装置制御方法
JP5160586B2 (ja) 脈波診断装置及び脈波診断装置制御方法
JP2010188152A (ja) 脈波診断装置及び脈波診断装置制御方法
EP3854301B1 (en) Apparatus and method for estimating bio-information
EP2979632A1 (en) Neurogenic baroreflex sensitivity measurement device, neurogenic baroreflex sensitivity measurement program and neurogenic baroreflex sensitivity measurement method
JP6197926B2 (ja) 生体情報検出装置及び生体情報検出方法、生体情報検出プログラム
EP4173571B1 (en) Apparatus for estimating bio-information, and method of determining false detection of bio-signal peaks
JP2017225839A (ja) 生体情報検出装置及び生体情報検出方法、生体情報検出プログラム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070222

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100216

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100413

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100525

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20101005