JP2001245859A - 生体信号検出装置及び非観血血圧計 - Google Patents
生体信号検出装置及び非観血血圧計Info
- Publication number
- JP2001245859A JP2001245859A JP2000270294A JP2000270294A JP2001245859A JP 2001245859 A JP2001245859 A JP 2001245859A JP 2000270294 A JP2000270294 A JP 2000270294A JP 2000270294 A JP2000270294 A JP 2000270294A JP 2001245859 A JP2001245859 A JP 2001245859A
- Authority
- JP
- Japan
- Prior art keywords
- disturbance
- signal
- biological information
- biological
- separation
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/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/021—Measuring pressure in heart or blood vessels
- A61B5/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7239—Details of waveform analysis using differentiation including higher order derivatives
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Vascular Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Physiology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ophthalmology & Optometry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
(57)【要約】
【課題】 センサ２で検出される生体情報に体動による
外乱が含まれている場合に、その外乱を生体情報から分
離すること。 【解決手段】 生体信号検出装置１は、Ａ／Ｄ変換器
３、分離マトリクス計算手段４、信号分離計算手段５、
及び外部出力手段６より構成される。分離マトリクス計
算手段４は、Ａ／Ｄ変換器３から送られる現在時点のデ
ジタル信号を用いてメモリにある現在時点の分離マトリ
クスＷ(n) をΔＷ(n) だけ修正してメモリの値を更新す
る。信号分離計算手段５は、分離マトリクス計算手段４
から送られた分離マトリクスとデジタル信号との積を計
算することで、体動による外乱を含まない生体信号と体
動による外乱を示す体動信号とから成るベクトルｙ(n)
を求め、外部出力手段６へ送る。外部出力手段６は、サ
ンプリング周波数１００Ｈｚに同期して、体動による外
乱が分離された生体信号と体動信号を外部へ出力する。
外乱が含まれている場合に、その外乱を生体情報から分
離すること。 【解決手段】 生体信号検出装置１は、Ａ／Ｄ変換器
３、分離マトリクス計算手段４、信号分離計算手段５、
及び外部出力手段６より構成される。分離マトリクス計
算手段４は、Ａ／Ｄ変換器３から送られる現在時点のデ
ジタル信号を用いてメモリにある現在時点の分離マトリ
クスＷ(n) をΔＷ(n) だけ修正してメモリの値を更新す
る。信号分離計算手段５は、分離マトリクス計算手段４
から送られた分離マトリクスとデジタル信号との積を計
算することで、体動による外乱を含まない生体信号と体
動による外乱を示す体動信号とから成るベクトルｙ(n)
を求め、外部出力手段６へ送る。外部出力手段６は、サ
ンプリング周波数１００Ｈｚに同期して、体動による外
乱が分離された生体信号と体動信号を外部へ出力する。
Description
【０００１】
【発明の属する技術分野】本発明は、生体情報から外乱
を含まない生体信号を検出できる生体信号検出装置、及
び血圧変動のみを計測する原理・方式の非観血血圧計に
関する。
を含まない生体信号を検出できる生体信号検出装置、及
び血圧変動のみを計測する原理・方式の非観血血圧計に
関する。
【０００２】
【従来の技術】従来より、カフを用いて連続的に血圧測
定を行える血圧計が知られているが、この血圧計では、
測定間隔を短くすると、カフによる締め付けが頻繁に繰
り返されることから、被験者に大きな負担を与え、時に
はカフを巻いた部分で内出血を起こすこともある。これ
に対し、血圧の絶対値を直接測定するのではなく、血圧
変動のみを計測する原理・方式の血圧計がある。例え
ば、脈波が動脈中を進行する速度（脈波伝播速度）や脈
波波形の特徴量等を捉えてパラメータ化し、その変化量
を血圧の変化量に換算して表示する血圧計がある。この
種の血圧計は、カフによる締め付けを繰り返す必要がな
いため、被験者に大きな負担を与えず、原理的に１拍の
心電や脈波信号によって血圧を算出できることから、連
続的な血圧測定を無侵襲で行うことが可能である。
定を行える血圧計が知られているが、この血圧計では、
測定間隔を短くすると、カフによる締め付けが頻繁に繰
り返されることから、被験者に大きな負担を与え、時に
はカフを巻いた部分で内出血を起こすこともある。これ
に対し、血圧の絶対値を直接測定するのではなく、血圧
変動のみを計測する原理・方式の血圧計がある。例え
ば、脈波が動脈中を進行する速度（脈波伝播速度）や脈
波波形の特徴量等を捉えてパラメータ化し、その変化量
を血圧の変化量に換算して表示する血圧計がある。この
種の血圧計は、カフによる締め付けを繰り返す必要がな
いため、被験者に大きな負担を与えず、原理的に１拍の
心電や脈波信号によって血圧を算出できることから、連
続的な血圧測定を無侵襲で行うことが可能である。
【０００３】しかし、上記の血圧計（血圧変動のみを計
測する血圧計）は、脈波の変化量を元に血圧を算出する
ため、体動（体の動き）等により脈波波形に外乱が載る
と、正確な血圧測定ができないという問題がある。この
対策として、以下の従来技術が知られている。 体動ノイズの発生を判断する判断手段を有し、この判
断手段により体動ノイズの発生が判断された場合は、セ
ンサ信号の読み込みを所定時間待機させることで、体動
ノイズの読み込みによる血圧測定の精度低下を防止する
（特許番号第２６６４９８３号参照）。
測する血圧計）は、脈波の変化量を元に血圧を算出する
ため、体動（体の動き）等により脈波波形に外乱が載る
と、正確な血圧測定ができないという問題がある。この
対策として、以下の従来技術が知られている。 体動ノイズの発生を判断する判断手段を有し、この判
断手段により体動ノイズの発生が判断された場合は、セ
ンサ信号の読み込みを所定時間待機させることで、体動
ノイズの読み込みによる血圧測定の精度低下を防止する
（特許番号第２６６４９８３号参照）。
【０００４】複数の脈波の周期を用いて現在の脈波の
基本周期を算出し、この基本周期から算出される脈波周
期のばらつきに応じて異常脈波の振幅を補正する方法
（特開平５−１８４５４７号参照）。 人体の生体情報から体動による外乱を除去する方法。
つまり、体動を含む一般のアーチファクト（artifact）
の影響を脈波や心電図から取り除くもので、脈波の観測
時点から６４拍分の過去の脈波を見て、脈波振幅の分散
が大きい部分を除去し、それ以外の信号部分を解析に使
用する（異なる測定部位の脈波伝播時間と血圧の関係の
比較、第３８回日本ＭＥ学会大会論文集 Sendai,21-23
April,1999 pp.373）。
基本周期を算出し、この基本周期から算出される脈波周
期のばらつきに応じて異常脈波の振幅を補正する方法
（特開平５−１８４５４７号参照）。 人体の生体情報から体動による外乱を除去する方法。
つまり、体動を含む一般のアーチファクト（artifact）
の影響を脈波や心電図から取り除くもので、脈波の観測
時点から６４拍分の過去の脈波を見て、脈波振幅の分散
が大きい部分を除去し、それ以外の信号部分を解析に使
用する（異なる測定部位の脈波伝播時間と血圧の関係の
比較、第３８回日本ＭＥ学会大会論文集 Sendai,21-23
April,1999 pp.373）。
【０００５】
【発明が解決しようとする課題】ところが、上記の方
法では、頻繁に体動ノイズが載る自由行動下での血圧測
定には適さない。の方法では、単発的な異常脈波には
対応できるが、頻繁に脈波が乱れる自由行動下では、脈
波の基本周波数自体を求めることが困難であり、の方
法と同様に、自由行動下での血圧測定には適さない。
の方法では、体動などの外乱が入った部分を使用しない
ため、解析に使用可能な生体信号が限られてしまう。ま
た、生体信号の時間的な連続性も失われてしまう。更
に、体動などの外乱が長時間にわたって加わると、生体
信号の解析が不可能となる。以上の結果、血圧変動のみ
を計測する原理・方式の非観血血圧計では、自由行動下
で頻繁に発生する体動により脈波波形が大きく乱れた場
合に、十分な血圧推定精度を確保できない。
法では、頻繁に体動ノイズが載る自由行動下での血圧測
定には適さない。の方法では、単発的な異常脈波には
対応できるが、頻繁に脈波が乱れる自由行動下では、脈
波の基本周波数自体を求めることが困難であり、の方
法と同様に、自由行動下での血圧測定には適さない。
の方法では、体動などの外乱が入った部分を使用しない
ため、解析に使用可能な生体信号が限られてしまう。ま
た、生体信号の時間的な連続性も失われてしまう。更
に、体動などの外乱が長時間にわたって加わると、生体
信号の解析が不可能となる。以上の結果、血圧変動のみ
を計測する原理・方式の非観血血圧計では、自由行動下
で頻繁に発生する体動により脈波波形が大きく乱れた場
合に、十分な血圧推定精度を確保できない。
【０００６】本発明は、上記事情に基づいて成されたも
ので、その第１の目的は、生体情報から外乱を含まない
生体信号を検出でき、且つその生体信号の時間的な連続
性を維持できる生体信号検出装置を提供すること。ま
た、第２の目的は、血圧変動のみを計測する原理・方式
の非観血血圧計において、心拍に同期する脈動成分から
体動によるノイズを分離して、自由行動下でも正確に且
つ連続的に血圧測定を可能にすることである。
ので、その第１の目的は、生体情報から外乱を含まない
生体信号を検出でき、且つその生体信号の時間的な連続
性を維持できる生体信号検出装置を提供すること。ま
た、第２の目的は、血圧変動のみを計測する原理・方式
の非観血血圧計において、心拍に同期する脈動成分から
体動によるノイズを分離して、自由行動下でも正確に且
つ連続的に血圧測定を可能にすることである。
【０００７】
【課題を解決するための手段】本発明の生体信号検出装
置は、上記の課題を解決する技術的手段として以下の構
成を有している。請求項１では、生体情報を検出する２
つ以上のセンサから生体情報をアナログ信号として入力
し、そのアナログ信号をデジタル信号に変換するＡ／Ｄ
変換手段と、センサで検出された生体情報に体動による
外乱が含まれており、その外乱を含んだ生体情報がセン
サで検出されるまでに振幅のみ変化して、体動に伴う外
乱と外乱を含んでいない生体情報とが混ざっていると見
做される場合に、外乱を含んでいない生体情報と外乱と
の統計的独立性を利用して、外乱を含んだ生体情報から
外乱を分離するための１つの分離マトリクスを計算する
分離マトリクス計算手段と、Ａ／Ｄ変換手段で変換され
たデジタル信号に分離マトリクスを作用させて、デジタ
ル信号を外乱を含まない生体信号と外乱を示す体動信号
とに分離する信号分離計算手段とを備える。
置は、上記の課題を解決する技術的手段として以下の構
成を有している。請求項１では、生体情報を検出する２
つ以上のセンサから生体情報をアナログ信号として入力
し、そのアナログ信号をデジタル信号に変換するＡ／Ｄ
変換手段と、センサで検出された生体情報に体動による
外乱が含まれており、その外乱を含んだ生体情報がセン
サで検出されるまでに振幅のみ変化して、体動に伴う外
乱と外乱を含んでいない生体情報とが混ざっていると見
做される場合に、外乱を含んでいない生体情報と外乱と
の統計的独立性を利用して、外乱を含んだ生体情報から
外乱を分離するための１つの分離マトリクスを計算する
分離マトリクス計算手段と、Ａ／Ｄ変換手段で変換され
たデジタル信号に分離マトリクスを作用させて、デジタ
ル信号を外乱を含まない生体信号と外乱を示す体動信号
とに分離する信号分離計算手段とを備える。
【０００８】請求項２では、生体情報を検出する２つ以
上のセンサから生体情報をアナログ信号として入力し、
そのアナログ信号をデジタル信号に変換するＡ／Ｄ変換
手段と、センサで検出された生体情報に体動による外乱
が含まれており、その外乱を含んだ生体情報がセンサで
検出されるまでに周波数と位相が線形的に変化して、体
動に伴う外乱と外乱を含んでいない生体情報とが混ざっ
ていると見做される場合に、外乱を含んでいない生体情
報と外乱との統計的独立性を利用して、外乱を含んだ生
体情報から外乱を分離するための２つ以上の分離マトリ
クスを計算する分離マトリクス計算手段と、Ａ／Ｄ変換
手段で変換されたデジタル信号に２つ以上の分離マトリ
クスを作用させて、デジタル信号を外乱を含まない生体
信号と外乱を示す体動信号とに分離する信号分離計算手
段とを備える。
上のセンサから生体情報をアナログ信号として入力し、
そのアナログ信号をデジタル信号に変換するＡ／Ｄ変換
手段と、センサで検出された生体情報に体動による外乱
が含まれており、その外乱を含んだ生体情報がセンサで
検出されるまでに周波数と位相が線形的に変化して、体
動に伴う外乱と外乱を含んでいない生体情報とが混ざっ
ていると見做される場合に、外乱を含んでいない生体情
報と外乱との統計的独立性を利用して、外乱を含んだ生
体情報から外乱を分離するための２つ以上の分離マトリ
クスを計算する分離マトリクス計算手段と、Ａ／Ｄ変換
手段で変換されたデジタル信号に２つ以上の分離マトリ
クスを作用させて、デジタル信号を外乱を含まない生体
信号と外乱を示す体動信号とに分離する信号分離計算手
段とを備える。
【０００９】請求項３では、生体情報を検出する２つ以
上のセンサから生体情報をアナログ信号として入力し、
そのアナログ信号をデジタル信号に変換するＡ／Ｄ変換
手段と、センサから体動による外乱を含んでいない生体
情報を入力し、その生体情報の確率密度関数を推定する
確率密度関数推定手段と、センサで検出された生体情報
に体動による外乱が含まれており、その外乱を含んだ生
体情報がセンサで検出されるまでに振幅のみ変化して、
体動に伴う外乱と外乱を含んでいない生体情報とが混ざ
っていると見做される場合に、確率密度関数を用い、外
乱を含んでいない生体情報と外乱との統計的独立性を利
用して、外乱を含んだ生体情報から外乱を分離するため
の１つの分離マトリクスを計算する分離マトリクス計算
手段と、Ａ／Ｄ変換手段で変換されたデジタル信号に分
離マトリクスを作用させて、デジタル信号を外乱を含ま
ない生体信号と外乱を示す体動信号とに分離する信号分
離計算手段とを備える。
上のセンサから生体情報をアナログ信号として入力し、
そのアナログ信号をデジタル信号に変換するＡ／Ｄ変換
手段と、センサから体動による外乱を含んでいない生体
情報を入力し、その生体情報の確率密度関数を推定する
確率密度関数推定手段と、センサで検出された生体情報
に体動による外乱が含まれており、その外乱を含んだ生
体情報がセンサで検出されるまでに振幅のみ変化して、
体動に伴う外乱と外乱を含んでいない生体情報とが混ざ
っていると見做される場合に、確率密度関数を用い、外
乱を含んでいない生体情報と外乱との統計的独立性を利
用して、外乱を含んだ生体情報から外乱を分離するため
の１つの分離マトリクスを計算する分離マトリクス計算
手段と、Ａ／Ｄ変換手段で変換されたデジタル信号に分
離マトリクスを作用させて、デジタル信号を外乱を含ま
ない生体信号と外乱を示す体動信号とに分離する信号分
離計算手段とを備える。
【００１０】請求項４では、生体情報を検出する２つ以
上のセンサから生体情報をアナログ信号として入力し、
そのアナログ信号をデジタル信号に変換するＡ／Ｄ変換
手段と、センサから体動による外乱を含んでいない生体
情報を入力し、その生体情報の確率密度関数を推定する
確率密度関数推定手段と、センサで検出された生体情報
に体動による外乱が含まれており、その外乱を含んだ生
体情報がセンサで検出されるまでに周波数と位相が線形
的に変化して、体動に伴う外乱と外乱を含んでいない生
体情報とが混ざっていると見做される場合に、確率密度
関数を用い、外乱を含んでいない生体情報と外乱との統
計的独立性を利用して、外乱を含んだ生体情報から外乱
を分離するための２つ以上の分離マトリクスを計算する
分離マトリクス計算手段と、Ａ／Ｄ変換手段で変換され
たデジタル信号に２つ以上の分離マトリクスを作用させ
て、デジタル信号を外乱を含まない生体信号と外乱を示
す体動信号とに分離する信号分離計算手段とを備える。
上のセンサから生体情報をアナログ信号として入力し、
そのアナログ信号をデジタル信号に変換するＡ／Ｄ変換
手段と、センサから体動による外乱を含んでいない生体
情報を入力し、その生体情報の確率密度関数を推定する
確率密度関数推定手段と、センサで検出された生体情報
に体動による外乱が含まれており、その外乱を含んだ生
体情報がセンサで検出されるまでに周波数と位相が線形
的に変化して、体動に伴う外乱と外乱を含んでいない生
体情報とが混ざっていると見做される場合に、確率密度
関数を用い、外乱を含んでいない生体情報と外乱との統
計的独立性を利用して、外乱を含んだ生体情報から外乱
を分離するための２つ以上の分離マトリクスを計算する
分離マトリクス計算手段と、Ａ／Ｄ変換手段で変換され
たデジタル信号に２つ以上の分離マトリクスを作用させ
て、デジタル信号を外乱を含まない生体信号と外乱を示
す体動信号とに分離する信号分離計算手段とを備える。
【００１１】請求項５では、請求項１〜４に記載した何
れかの生体信号検出装置において、信号分離計算手段で
分離された生体信号を外部に出力する信号出力手段を具
備している。
れかの生体信号検出装置において、信号分離計算手段で
分離された生体信号を外部に出力する信号出力手段を具
備している。
【００１２】請求項６では、請求項１〜５に記載した何
れかの生体信号検出装置において、センサから入力した
生体情報に体動による外乱が含まれているか否かを判定
する体動判定手段と、この体動判定手段で体動による外
乱が含まれていないと判定された場合に、Ａ／Ｄ変換手
段で変換されたデジタル信号を分離マトリクス計算手段
及び信号分離計算手段を経ずに出力する信号バイパス手
段とを備えている。
れかの生体信号検出装置において、センサから入力した
生体情報に体動による外乱が含まれているか否かを判定
する体動判定手段と、この体動判定手段で体動による外
乱が含まれていないと判定された場合に、Ａ／Ｄ変換手
段で変換されたデジタル信号を分離マトリクス計算手段
及び信号分離計算手段を経ずに出力する信号バイパス手
段とを備えている。
【００１３】請求項７では、請求項３〜５に記載した何
れかの生体信号検出装置において、センサから入力した
生体情報に体動による外乱が含まれているか否かを判定
する体動判定手段を備え、確率密度関数推定手段は、体
動判定手段で体動による外乱が含まれていないと判定さ
れた場合に、センサから入力した生体情報の確率密度関
数を推定する。
れかの生体信号検出装置において、センサから入力した
生体情報に体動による外乱が含まれているか否かを判定
する体動判定手段を備え、確率密度関数推定手段は、体
動判定手段で体動による外乱が含まれていないと判定さ
れた場合に、センサから入力した生体情報の確率密度関
数を推定する。
【００１４】請求項８では、請求項７に記載した生体信
号検出装置において、体動判定手段で体動による外乱が
含まれていないと判定された場合に、Ａ／Ｄ変換手段で
変換されたデジタル信号を分離マトリクス計算手段及び
信号分離計算手段を経ずに出力する信号バイパス手段を
備えている。
号検出装置において、体動判定手段で体動による外乱が
含まれていないと判定された場合に、Ａ／Ｄ変換手段で
変換されたデジタル信号を分離マトリクス計算手段及び
信号分離計算手段を経ずに出力する信号バイパス手段を
備えている。
【００１５】請求項９では、請求項１〜８に記載した生
体信号検出装置は、被験者の脈波（生体情報）から外乱
を含まない生体信号を検出するものである。
体信号検出装置は、被験者の脈波（生体情報）から外乱
を含まない生体信号を検出するものである。
【００１６】また、本発明の非観血血圧計は、上記の課
題を解決する技術的手段として以下の構成を有してい
る。請求項１０では、被験者の身体表面に保持され、心
拍に同期して得られる脈動成分を検出してアナログ信号
として出力する少なくとも２つ以上のセンサと、センサ
からアナログ信号を入力し、そのアナログ信号をデジタ
ル信号に変換するＡ／Ｄ変換手段と、センサで検出され
た脈動成分に体動による外乱が含まれており、その外乱
を含んだ脈動成分がセンサで検出されるまでに周波数と
位相が線形的に変化して、体動に伴う外乱と外乱を含ん
でいない心拍に同期した脈動成分とが混ざっていると見
做される場合に、外乱を含んでいない脈動成分と外乱と
の統計的独立性を利用して、外乱を含んだ脈動成分から
外乱を分離するための複数の分離マトリクスを計算する
分離マトリクス計算手段と、Ａ／Ｄ変換手段で変換され
たデジタル信号に複数の分離マトリクスを作用させて、
デジタル信号を、外乱を含まない心拍に同期した脈動成
分と外乱を示す体動成分とに分離する信号分離計算手段
と、この信号分離計算手段から取り込まれる外乱を含ま
ない心拍に同期した脈動成分を元に、連続的に血圧値を
計算する血圧値計算手段とを有している。
題を解決する技術的手段として以下の構成を有してい
る。請求項１０では、被験者の身体表面に保持され、心
拍に同期して得られる脈動成分を検出してアナログ信号
として出力する少なくとも２つ以上のセンサと、センサ
からアナログ信号を入力し、そのアナログ信号をデジタ
ル信号に変換するＡ／Ｄ変換手段と、センサで検出され
た脈動成分に体動による外乱が含まれており、その外乱
を含んだ脈動成分がセンサで検出されるまでに周波数と
位相が線形的に変化して、体動に伴う外乱と外乱を含ん
でいない心拍に同期した脈動成分とが混ざっていると見
做される場合に、外乱を含んでいない脈動成分と外乱と
の統計的独立性を利用して、外乱を含んだ脈動成分から
外乱を分離するための複数の分離マトリクスを計算する
分離マトリクス計算手段と、Ａ／Ｄ変換手段で変換され
たデジタル信号に複数の分離マトリクスを作用させて、
デジタル信号を、外乱を含まない心拍に同期した脈動成
分と外乱を示す体動成分とに分離する信号分離計算手段
と、この信号分離計算手段から取り込まれる外乱を含ま
ない心拍に同期した脈動成分を元に、連続的に血圧値を
計算する血圧値計算手段とを有している。
【００１７】請求項１１では、請求項１０に記載した非
観血血圧計において、別途測定された被験者の血圧値を
基準血圧として入力及び記憶する入力・記憶手段と、信
号分離計算手段によって得られた外乱を含まない心拍に
同期した脈動成分と基準血圧とを対応させるための係数
を求める係数計算手段とを有し、血圧値計算手段は、信
号分離計算手段から取り込まれる外乱を含まない心拍に
同期した脈動成分と係数計算手段によって算出された係
数とから連続的に血圧値を計算することを特徴とする。
観血血圧計において、別途測定された被験者の血圧値を
基準血圧として入力及び記憶する入力・記憶手段と、信
号分離計算手段によって得られた外乱を含まない心拍に
同期した脈動成分と基準血圧とを対応させるための係数
を求める係数計算手段とを有し、血圧値計算手段は、信
号分離計算手段から取り込まれる外乱を含まない心拍に
同期した脈動成分と係数計算手段によって算出された係
数とから連続的に血圧値を計算することを特徴とする。
【００１８】請求項１２では、請求項１０に記載した非
観血血圧計において、別途測定された被験者の血圧値及
び外乱を含まない心拍に同期した脈動成分の伝播時間を
それぞれ基準血圧及び基準脈波伝播時間として入力及び
記憶する入力・記憶手段と、基準血圧と基準脈波伝播時
間とを対応させるための係数を求める係数計算手段と、
信号分離計算手段によって得られた外乱を含まない心拍
に同期した脈動成分の伝播時間を算出する脈波伝播時間
計算手段と、血圧値計算手段は、脈波伝播時間計算手段
から取り込まれる脈動成分の伝播時間と係数計算手段に
よって算出された係数とから連続的に血圧値を計算する
ことを特徴とする。
観血血圧計において、別途測定された被験者の血圧値及
び外乱を含まない心拍に同期した脈動成分の伝播時間を
それぞれ基準血圧及び基準脈波伝播時間として入力及び
記憶する入力・記憶手段と、基準血圧と基準脈波伝播時
間とを対応させるための係数を求める係数計算手段と、
信号分離計算手段によって得られた外乱を含まない心拍
に同期した脈動成分の伝播時間を算出する脈波伝播時間
計算手段と、血圧値計算手段は、脈波伝播時間計算手段
から取り込まれる脈動成分の伝播時間と係数計算手段に
よって算出された係数とから連続的に血圧値を計算する
ことを特徴とする。
【００１９】請求項１３では、請求項１０〜１２に記載
した非観血血圧計において、被験者の身体に装着される
カフを有し、このカフによりセンサを被験者の拡張期圧
以下の圧力で保持することを特徴とする。
した非観血血圧計において、被験者の身体に装着される
カフを有し、このカフによりセンサを被験者の拡張期圧
以下の圧力で保持することを特徴とする。
【００２０】請求項１４では、請求項１０〜１３に記載
した非観血血圧計において、２つ以上のセンサは、心拍
に同期して得られる脈動成分を同一血管系上から検出で
きる位置に配置されることを特徴とする。
した非観血血圧計において、２つ以上のセンサは、心拍
に同期して得られる脈動成分を同一血管系上から検出で
きる位置に配置されることを特徴とする。
【００２１】請求項１５では、請求項１４に記載した非
観血血圧計において、２つ以上のセンサは、同一血管系
上に沿って直線的に配置されることを特徴とする。
観血血圧計において、２つ以上のセンサは、同一血管系
上に沿って直線的に配置されることを特徴とする。
【００２２】請求項１６では、請求項１４に記載した非
観血血圧計において、２つ以上のセンサは、同一血管系
と直交する周回方向に配置されることを特徴とする。
観血血圧計において、２つ以上のセンサは、同一血管系
と直交する周回方向に配置されることを特徴とする。
【００２３】
【作用】センサで検出された生体情報に体動による外乱
が含まれている場合は、その外乱を含んだ生体情報がセ
ンサで検出されるまでに振幅のみが変化して、体動に伴
う外乱と外乱を含んでいない生体情報とが混ざっている
と見做される場合と、線形システムを通ることにより周
波数と位相とが変化して、体動に伴う外乱と外乱を含ん
でいない生体情報とが混ざっていると見做される場合が
考えられる。
が含まれている場合は、その外乱を含んだ生体情報がセ
ンサで検出されるまでに振幅のみが変化して、体動に伴
う外乱と外乱を含んでいない生体情報とが混ざっている
と見做される場合と、線形システムを通ることにより周
波数と位相とが変化して、体動に伴う外乱と外乱を含ん
でいない生体情報とが混ざっていると見做される場合が
考えられる。
【００２４】このような統計的に独立な複数の信号が振
幅のみの変化で加法的に混合した場合、及び線形システ
ムを通って加法的に混合した場合は、逆マトリクスを推
定するアルゴリズムからセンサで検出された生体情報の
逆マトリクスを推定し、この逆マトリクスを生体情報に
作用させることで、元の信号を回復することが原理的に
可能である（アルゴリズムの詳細は、例えば“ An info
rmation-maximizationapproach to blind separation a
nd blind deconvolution,”Anthony J.Bell and Terren
ce J.Sejnowsky,Technical Report no.INC-9501,Februa
ry 1995,Institute for Neural Computation,UCSD,San
Diego,CA 92093-0523を参照）。
幅のみの変化で加法的に混合した場合、及び線形システ
ムを通って加法的に混合した場合は、逆マトリクスを推
定するアルゴリズムからセンサで検出された生体情報の
逆マトリクスを推定し、この逆マトリクスを生体情報に
作用させることで、元の信号を回復することが原理的に
可能である（アルゴリズムの詳細は、例えば“ An info
rmation-maximizationapproach to blind separation a
nd blind deconvolution,”Anthony J.Bell and Terren
ce J.Sejnowsky,Technical Report no.INC-9501,Februa
ry 1995,Institute for Neural Computation,UCSD,San
Diego,CA 92093-0523を参照）。
【００２５】そこで、センサを２つとした場合に、体動
が加わっていない生体情報（生体信号）をｓ1(n)、体動
による外乱（体動信号）をｓ2(n)とし、センサ信号をｘ
1(n)、ｘ2(n)とすると（n:デジタル信号のインデック
ス、n=0,1,2,3 …）、振幅のみが変化して加わった場合
は、 ｘ1(n)＝ａ11×ｓ1(n)＋ａ12×ｓ2(n)…………………………(1) ｘ2(n)＝ａ21×ｓ1(n)＋ａ22×ｓ2(n)…………………………(2) と記述できる。
が加わっていない生体情報（生体信号）をｓ1(n)、体動
による外乱（体動信号）をｓ2(n)とし、センサ信号をｘ
1(n)、ｘ2(n)とすると（n:デジタル信号のインデック
ス、n=0,1,2,3 …）、振幅のみが変化して加わった場合
は、 ｘ1(n)＝ａ11×ｓ1(n)＋ａ12×ｓ2(n)…………………………(1) ｘ2(n)＝ａ21×ｓ1(n)＋ａ22×ｓ2(n)…………………………(2) と記述できる。
【００２６】また、線形システムを通って変形して加わ
った場合は、 ｘ1(n)＝ａ0,11×ｓ1(n)＋ａ0,12×ｓ2(n) ＋ａ1,11×ｓ1(n-1)＋ａ1,12×ｓ2(n-1)＋… ＋ａL,11×ｓ1(n-L)＋ａL,12×ｓ2(n-L)……………(3) ｘ2(n)＝ａ0,21×ｓ1(n)＋ａ0,22×ｓ2(n) ＋ａ1,21×ｓ1(n-1)＋ａ1,22×ｓ2(n-1)＋… ＋ａL,21×ｓ1(n-L)＋ａL,22×ｓ2(n-L)……………(4) と記述できる。但し、L は１以上の正の整数である。
った場合は、 ｘ1(n)＝ａ0,11×ｓ1(n)＋ａ0,12×ｓ2(n) ＋ａ1,11×ｓ1(n-1)＋ａ1,12×ｓ2(n-1)＋… ＋ａL,11×ｓ1(n-L)＋ａL,12×ｓ2(n-L)……………(3) ｘ2(n)＝ａ0,21×ｓ1(n)＋ａ0,22×ｓ2(n) ＋ａ1,21×ｓ1(n-1)＋ａ1,22×ｓ2(n-1)＋… ＋ａL,21×ｓ1(n-L)＋ａL,22×ｓ2(n-L)……………(4) と記述できる。但し、L は１以上の正の整数である。
【００２７】上記の(1) 及び(2) の場合は、マトリクス
Ａ（Ａ11＝ａ11、Ａ12＝ａ12、Ａ21＝ａ21、Ａ22＝
ａ22）の逆マトリクスを、｛ｘ1(n)、ｘ2(n)｝t （t ：
転置）に作用させれば、体動のない生体信号と体動信号
とが得られる。また、(3) 及び(4) の場合は、マトリク
スＡ0(Ａ011=ａ1,11、Ａ012=ａ0,12、Ａ021=ａ0,21、Ａ
022=ａ0,22) 、Ａ1(Ａ111=ａ1,11、Ａ112=ａ1,12、Ａ
121=ａ1, 21、Ａ122=ａ1,22) 、……、ＡL ( ＡL11=ａ
L,11、ＡL12=ａL,12、ＡL21=ａL,21、ＡL22=ａL,22) で
決まるフィルタマトリクスの逆フィルタマトリクスを作
用させれば良い。逆フィルタマトリクスは、２つ以上の
マトリクスから成る。
Ａ（Ａ11＝ａ11、Ａ12＝ａ12、Ａ21＝ａ21、Ａ22＝
ａ22）の逆マトリクスを、｛ｘ1(n)、ｘ2(n)｝t （t ：
転置）に作用させれば、体動のない生体信号と体動信号
とが得られる。また、(3) 及び(4) の場合は、マトリク
スＡ0(Ａ011=ａ1,11、Ａ012=ａ0,12、Ａ021=ａ0,21、Ａ
022=ａ0,22) 、Ａ1(Ａ111=ａ1,11、Ａ112=ａ1,12、Ａ
121=ａ1, 21、Ａ122=ａ1,22) 、……、ＡL ( ＡL11=ａ
L,11、ＡL12=ａL,12、ＡL21=ａL,21、ＡL22=ａL,22) で
決まるフィルタマトリクスの逆フィルタマトリクスを作
用させれば良い。逆フィルタマトリクスは、２つ以上の
マトリクスから成る。
【００２８】なお、逆マトリクスを推定する際に、ｇ
(x)=１／｛１＋exp(-x) ｝等の非線形関数を使用する
が、この非線形関数を微分した関数は、求めたい生体信
号の確率密度関数に一致するものが信号の分離の観点か
ら見て最も望ましい。しかし、ｇ(x)=１／｛１＋exp(-
x) ｝等の非線形関数でも信号の分離は達成される。
(x)=１／｛１＋exp(-x) ｝等の非線形関数を使用する
が、この非線形関数を微分した関数は、求めたい生体信
号の確率密度関数に一致するものが信号の分離の観点か
ら見て最も望ましい。しかし、ｇ(x)=１／｛１＋exp(-
x) ｝等の非線形関数でも信号の分離は達成される。
【００２９】上記の作用により、体動による外乱成分を
含まない心拍に同期した脈波信号が得られれば、脈波の
特徴変化から、例えば下記の方法により連続的に血圧を
推定することができる。 ａ）被測定部における脈波の波形自体が心臓の拍動に基
づく血圧値に依存することを利用して、脈波の波形変化
から血圧を計算する方法。 ｂ）末梢血管側では、脈波の伝播時刻は心臓の拍動に比
べて遅れが生じる（脈波伝播時間）。一定距離の脈波伝
播時間に対する脈波伝播速度は、血管の容積弾性率の関
数として表され、血圧が上がると血管の容積弾性率は増
加し、血管壁が硬くなって伝播速度が速くなる。従っ
て、体動による外乱成分を含まない心拍に同期した脈波
信号が得られれば、脈波伝播時間から血圧変動を求める
ことができる。
含まない心拍に同期した脈波信号が得られれば、脈波の
特徴変化から、例えば下記の方法により連続的に血圧を
推定することができる。 ａ）被測定部における脈波の波形自体が心臓の拍動に基
づく血圧値に依存することを利用して、脈波の波形変化
から血圧を計算する方法。 ｂ）末梢血管側では、脈波の伝播時刻は心臓の拍動に比
べて遅れが生じる（脈波伝播時間）。一定距離の脈波伝
播時間に対する脈波伝播速度は、血管の容積弾性率の関
数として表され、血圧が上がると血管の容積弾性率は増
加し、血管壁が硬くなって伝播速度が速くなる。従っ
て、体動による外乱成分を含まない心拍に同期した脈波
信号が得られれば、脈波伝播時間から血圧変動を求める
ことができる。
【００３０】
【発明の実施の形態】次に、本発明の生体信号検出装置
の実施形態を図面に基づいて説明する。 （第１実施例）図１は生体信号検出装置の構成を示すブ
ロック図である。この生体信号検出装置１は、２つのセ
ンサ２で検出される生体情報から体動による外乱を分離
して、外乱を含まない生体信号を出力する装置で、以下
に詳述するＡ／Ｄ変換器３、分離マトリクス計算手段
４、信号分離計算手段５、及び外部出力手段６より構成
される。
の実施形態を図面に基づいて説明する。 （第１実施例）図１は生体信号検出装置の構成を示すブ
ロック図である。この生体信号検出装置１は、２つのセ
ンサ２で検出される生体情報から体動による外乱を分離
して、外乱を含まない生体信号を出力する装置で、以下
に詳述するＡ／Ｄ変換器３、分離マトリクス計算手段
４、信号分離計算手段５、及び外部出力手段６より構成
される。
【００３１】センサ２は、生体情報として人体の脈波を
計測するもので、例えば発光素子と受光素子を備える周
知の光学式脈波センサである。このセンサ２は、例えば
人体の指等に装着され、発光素子から指に向かって光が
照射されると、光の一部が指の内部を通る毛細動脈に当
たって毛細動脈を流れる血液中のヘモグロビンに吸収さ
れ、残りの光が毛細動脈で反射して散乱し、一部が受光
素子に入射する。この時、血液の脈動により、ヘモグロ
ビンの量が波動的に変化するので、ヘモグロビンに吸収
される光の量も波動的に変化する。その結果、毛細動脈
で反射して受光素子で検出される受光量が変化し、その
受光量の変化を脈波情報（例えば電圧信号）としてＡ／
Ｄ変換器３に出力している。
計測するもので、例えば発光素子と受光素子を備える周
知の光学式脈波センサである。このセンサ２は、例えば
人体の指等に装着され、発光素子から指に向かって光が
照射されると、光の一部が指の内部を通る毛細動脈に当
たって毛細動脈を流れる血液中のヘモグロビンに吸収さ
れ、残りの光が毛細動脈で反射して散乱し、一部が受光
素子に入射する。この時、血液の脈動により、ヘモグロ
ビンの量が波動的に変化するので、ヘモグロビンに吸収
される光の量も波動的に変化する。その結果、毛細動脈
で反射して受光素子で検出される受光量が変化し、その
受光量の変化を脈波情報（例えば電圧信号）としてＡ／
Ｄ変換器３に出力している。
【００３２】Ａ／Ｄ変換器３は、２つのセンサ２から信
号を入力し、例えば１００Ｈｚのサンプリング周波数で
アナログ信号をデジタル信号ｘ1(n)、ｘ2(n)、(n=0,1,
2,,,,) に変換して分離マトリクス計算手段４へ送る。
分離マトリクス計算手段４は、ＣＰＵとメモリを有し、
Ａ／Ｄ変換器３から送られる現在時点のデジタル信号ｘ
1(n)、ｘ2(n)を用いてメモリにある現在時点の２×２分
離マトリクスＷ(n) をΔＷ(n) だけ修正し、Ｗ(n＋1)＝
Ｗ(n) ＋ΔＷ(n) としてメモリの値を更新する。なお、
本装置が起動した時点で、メモリのＷ(0) は２×２の単
位マトリクスにセットされる。２×２分離マトリクスの
修正分ΔＷ(n) は、以下の式で計算される。 ΔＷ(n) ＝[ Ｗ(n) t ]-1＋［Ｉ−２×ｇ｛Ｗ(n) ×ｘ
(n) ｝］×ｘ(n) t ｇ(x)=１／｛１＋exp(-x) ｝ 但し、ｘ(n) ＝｛ｘ1(n)、ｘ2(n)｝t 、（t ：転置） また、Ｉは１を成分とする二次元のベクトルである。Ｉ
＝｛１、１｝t
号を入力し、例えば１００Ｈｚのサンプリング周波数で
アナログ信号をデジタル信号ｘ1(n)、ｘ2(n)、(n=0,1,
2,,,,) に変換して分離マトリクス計算手段４へ送る。
分離マトリクス計算手段４は、ＣＰＵとメモリを有し、
Ａ／Ｄ変換器３から送られる現在時点のデジタル信号ｘ
1(n)、ｘ2(n)を用いてメモリにある現在時点の２×２分
離マトリクスＷ(n) をΔＷ(n) だけ修正し、Ｗ(n＋1)＝
Ｗ(n) ＋ΔＷ(n) としてメモリの値を更新する。なお、
本装置が起動した時点で、メモリのＷ(0) は２×２の単
位マトリクスにセットされる。２×２分離マトリクスの
修正分ΔＷ(n) は、以下の式で計算される。 ΔＷ(n) ＝[ Ｗ(n) t ]-1＋［Ｉ−２×ｇ｛Ｗ(n) ×ｘ
(n) ｝］×ｘ(n) t ｇ(x)=１／｛１＋exp(-x) ｝ 但し、ｘ(n) ＝｛ｘ1(n)、ｘ2(n)｝t 、（t ：転置） また、Ｉは１を成分とする二次元のベクトルである。Ｉ
＝｛１、１｝t
【００３３】マトリクスＷ(n) t の逆マトリクスは、例
えばＬＵ分解法（「数値計算法」長嶋秀世著：槙書店を
参照）などの標準的な数値計算法により求める。また、
関数ｇの値も数値計算法により求める。このようにして
修正されたある時点の分離マトリクスＷ(n＋1)が、一瞬
前の分離マトリクスＷ(n) の代わりにメモリに貯えられ
ると共に、元の分離マトリクスＷ(n) とＡ／Ｄ変換器３
で変換されたデジタル信号ｘ(n) が信号分離計算手段５
へ送られる。信号分離計算手段５は、ＣＰＵとメモリを
有しており、分離マトリクス計算手段４から送られた分
離マトリクスＷ(n) とデジタル信号ｘ(n) との積を過去
に逆上って加えることで、体動と分離された生体信号と
体動による外乱を示す体動信号とから成るベクトルｙ
(n) ＝｛ｙ1(n), ｙ2(n)｝t を求め、外部出力手段６へ
送る。
えばＬＵ分解法（「数値計算法」長嶋秀世著：槙書店を
参照）などの標準的な数値計算法により求める。また、
関数ｇの値も数値計算法により求める。このようにして
修正されたある時点の分離マトリクスＷ(n＋1)が、一瞬
前の分離マトリクスＷ(n) の代わりにメモリに貯えられ
ると共に、元の分離マトリクスＷ(n) とＡ／Ｄ変換器３
で変換されたデジタル信号ｘ(n) が信号分離計算手段５
へ送られる。信号分離計算手段５は、ＣＰＵとメモリを
有しており、分離マトリクス計算手段４から送られた分
離マトリクスＷ(n) とデジタル信号ｘ(n) との積を過去
に逆上って加えることで、体動と分離された生体信号と
体動による外乱を示す体動信号とから成るベクトルｙ
(n) ＝｛ｙ1(n), ｙ2(n)｝t を求め、外部出力手段６へ
送る。
【００３４】外部出力手段６は、サンプリング周波数１
００Ｈｚに同期して、体動による外乱が分離された生体
信号と体動信号を外部へ出力するモニタ装置である。以
上により、センサ２で検出される生体情報に体動による
外乱が加わっている場合でも、その外乱を含んだ生体情
報がセンサ２で検出されるまでに振幅のみ変化して、体
動に伴う外乱と外乱を含んでいない生体情報とが混ざっ
ていると見做される場合は、分離マトリクス計算手段４
で計算された１つの分離マトリクスを用いることで生体
情報から外乱を分離することができる。なお、外乱を含
んだ生体情報がセンサ２で検出されるまでに周波数と位
相が線形的に変化して、体動に伴う外乱と外乱を含んで
いない生体情報とが混ざっていると見做される場合に
は、分離マトリクス計算手段４で２つ以上の分離マトリ
クスを計算し、この２つ以上の分離マトリクスを用いて
生体情報から外乱を分離することができる。
００Ｈｚに同期して、体動による外乱が分離された生体
信号と体動信号を外部へ出力するモニタ装置である。以
上により、センサ２で検出される生体情報に体動による
外乱が加わっている場合でも、その外乱を含んだ生体情
報がセンサ２で検出されるまでに振幅のみ変化して、体
動に伴う外乱と外乱を含んでいない生体情報とが混ざっ
ていると見做される場合は、分離マトリクス計算手段４
で計算された１つの分離マトリクスを用いることで生体
情報から外乱を分離することができる。なお、外乱を含
んだ生体情報がセンサ２で検出されるまでに周波数と位
相が線形的に変化して、体動に伴う外乱と外乱を含んで
いない生体情報とが混ざっていると見做される場合に
は、分離マトリクス計算手段４で２つ以上の分離マトリ
クスを計算し、この２つ以上の分離マトリクスを用いて
生体情報から外乱を分離することができる。
【００３５】（第２実施例）図２は生体信号検出装置１
の構成を示すブロック図である。本実施例は、第１実施
例の構成に加えて、下述の設定モードに応じてＡ／Ｄ変
換器３の出力を切り替えるモード制御手段７と、体動に
よる外乱を含んでいない生体情報の確率密度関数（ヒス
トグラム）を推定する信号確率密度関数推定手段８を有
している。設定モードは、例えば外部のボタン操作によ
って選択できる訓練モードと分離モードがあり、何方か
一方のモードを選択すると、そのモードに対応したモー
ド制御信号（訓練モード信号と分離モード信号）が出力
される。但し、訓練モードは、被験者が静止した状態
（体動が無い状態）で実施される。
の構成を示すブロック図である。本実施例は、第１実施
例の構成に加えて、下述の設定モードに応じてＡ／Ｄ変
換器３の出力を切り替えるモード制御手段７と、体動に
よる外乱を含んでいない生体情報の確率密度関数（ヒス
トグラム）を推定する信号確率密度関数推定手段８を有
している。設定モードは、例えば外部のボタン操作によ
って選択できる訓練モードと分離モードがあり、何方か
一方のモードを選択すると、そのモードに対応したモー
ド制御信号（訓練モード信号と分離モード信号）が出力
される。但し、訓練モードは、被験者が静止した状態
（体動が無い状態）で実施される。
【００３６】モード制御手段７は、本発明の体動判定手
段でもあり、訓練モード信号を受け取ると、Ａ／Ｄ変換
器３の出力を信号確率密度関数推定手段８へ送り、分離
モード信号を受け取ると、Ａ／Ｄ変換器３の出力を分離
マトリクス計算手段４へ送る。信号確率密度関数推定手
段８は、ＣＰＵとメモリを有し、訓練モード時にＡ／Ｄ
変換器３から順次送られてくるデジタル信号、例えばｘ
1(n)の振幅のヒストグラムを作成する。ヒストグラムの
作成は、ｘ1(n)の可能な最大振幅Ｍと可能な最小振幅ｍ
との間を例えば１００分割し、ｘ1(n)の振幅がある分割
区間に入る時に、その分割区間のカウンタを「１」増や
すことで行う。このヒストグラムの作成は、モード制御
信号が分離モード信号に変わるまで行われる。また、訓
練モードから分離モードに変わると、推定したヒストグ
ラム（折れ線近似した関数）を分離マトリクス計算手段
４へ送る。
段でもあり、訓練モード信号を受け取ると、Ａ／Ｄ変換
器３の出力を信号確率密度関数推定手段８へ送り、分離
モード信号を受け取ると、Ａ／Ｄ変換器３の出力を分離
マトリクス計算手段４へ送る。信号確率密度関数推定手
段８は、ＣＰＵとメモリを有し、訓練モード時にＡ／Ｄ
変換器３から順次送られてくるデジタル信号、例えばｘ
1(n)の振幅のヒストグラムを作成する。ヒストグラムの
作成は、ｘ1(n)の可能な最大振幅Ｍと可能な最小振幅ｍ
との間を例えば１００分割し、ｘ1(n)の振幅がある分割
区間に入る時に、その分割区間のカウンタを「１」増や
すことで行う。このヒストグラムの作成は、モード制御
信号が分離モード信号に変わるまで行われる。また、訓
練モードから分離モードに変わると、推定したヒストグ
ラム（折れ線近似した関数）を分離マトリクス計算手段
４へ送る。
【００３７】分離マトリクス計算手段４では、第１実施
例における関数ｇとして、信号確率密度関数推定手段８
より送られたヒストグラム（折れ線近似した関数）を積
分した関数を使用する。その他の動作は、第１実施例と
同じである。本実施例では、被験者が静止した状態で検
出される生体情報からヒストグラムを作成し、このヒス
トグラムを使用して分離マトリクスを計算するので、体
動による外乱を精度良く分離できる。なお、この実施例
は、外乱を含んだ生体情報がセンサ２で検出されるまで
に振幅のみ変化して、体動に伴う外乱と外乱を含んでい
ない生体情報とが混ざっていると見做される場合、及び
外乱を含んだ生体情報がセンサ２で検出されるまでに周
波数と位相が線形的に変化して、体動に伴う外乱と外乱
を含んでいない生体情報とが混ざっていると見做される
場合のどちらにも適用できることは言うまでもない。
例における関数ｇとして、信号確率密度関数推定手段８
より送られたヒストグラム（折れ線近似した関数）を積
分した関数を使用する。その他の動作は、第１実施例と
同じである。本実施例では、被験者が静止した状態で検
出される生体情報からヒストグラムを作成し、このヒス
トグラムを使用して分離マトリクスを計算するので、体
動による外乱を精度良く分離できる。なお、この実施例
は、外乱を含んだ生体情報がセンサ２で検出されるまで
に振幅のみ変化して、体動に伴う外乱と外乱を含んでい
ない生体情報とが混ざっていると見做される場合、及び
外乱を含んだ生体情報がセンサ２で検出されるまでに周
波数と位相が線形的に変化して、体動に伴う外乱と外乱
を含んでいない生体情報とが混ざっていると見做される
場合のどちらにも適用できることは言うまでもない。
【００３８】（第３実施例）図３は生体信号検出装置１
の構成を示すブロック図である。本実施例は、上記の第
２実施例で説明した構成において、訓練モードが選択さ
れている時（モード制御手段７に訓練モード信号が入力
される時）は、被験者が静止状態であるため、基本的に
体動による外乱が加わることがない。そこで、訓練モー
ド時には、図３に示すように、Ａ／Ｄ変換器３の出力を
モード制御手段７から直接外部出力手段６へ送るバイパ
ス信号回路９を設けても良い。
の構成を示すブロック図である。本実施例は、上記の第
２実施例で説明した構成において、訓練モードが選択さ
れている時（モード制御手段７に訓練モード信号が入力
される時）は、被験者が静止状態であるため、基本的に
体動による外乱が加わることがない。そこで、訓練モー
ド時には、図３に示すように、Ａ／Ｄ変換器３の出力を
モード制御手段７から直接外部出力手段６へ送るバイパ
ス信号回路９を設けても良い。
【００３９】（第４実施例）図４は生体信号検出装置１
の構成を示すブロック図である。本実施例は、図１に示
した第１実施例の構成に加えて、センサ２で検出された
生体情報に体動による外乱が含まれているか否かを判定
する体動判定手段１０を設け、この体動判定手段１０で
体動による外乱が含まれていないと判定された場合は、
Ａ／Ｄ変換器３の出力を体動判定手段１０から直接外部
出力手段６へ送るバイパス信号回路９を設けても良い。
この場合、体動判定手段１０として、例えば体動の有無
を検出できる加速度センサを使用することができる。あ
るいは、第２実施例で説明したように、被験者が静止状
態の時に出力される信号（第２実施例では訓練モード信
号）を検出できる信号検出手段を設けても良い。
の構成を示すブロック図である。本実施例は、図１に示
した第１実施例の構成に加えて、センサ２で検出された
生体情報に体動による外乱が含まれているか否かを判定
する体動判定手段１０を設け、この体動判定手段１０で
体動による外乱が含まれていないと判定された場合は、
Ａ／Ｄ変換器３の出力を体動判定手段１０から直接外部
出力手段６へ送るバイパス信号回路９を設けても良い。
この場合、体動判定手段１０として、例えば体動の有無
を検出できる加速度センサを使用することができる。あ
るいは、第２実施例で説明したように、被験者が静止状
態の時に出力される信号（第２実施例では訓練モード信
号）を検出できる信号検出手段を設けても良い。
【００４０】次に、本発明の非観血血圧計の実施形態を
図面に基づいて説明する。 （第１実施例）図５は非観血血圧計の構成を示すブロッ
ク図である。 （非観血血圧計の構成）本実施例の非観血血圧計１１
は、図５に示す様に、被験者の脈波を検出する２つのセ
ンサ２、一定圧で２つのセンサ２を保持するカフ１２、
センサ信号（アナログ信号）をデジタル信号に変換する
Ａ／Ｄ変換器３、体動による外乱を分離する外乱分離装
置１３（後述する）、別途測定した被験者の血圧値を基
準血圧として入力・記憶する入力・記憶手段１４、外乱
を含まない心拍に同期した脈波信号と基準血圧とを対応
させるための係数を求める係数計算手段１５、算出され
た係数と外乱を含まない心拍に同期した脈波信号から連
続的に血圧値を算出する血圧値計算手段１６、算出され
た連続血圧を記録する記録装置１７、算出された連続血
圧を外部へ出力する出力装置１８から成る。
図面に基づいて説明する。 （第１実施例）図５は非観血血圧計の構成を示すブロッ
ク図である。 （非観血血圧計の構成）本実施例の非観血血圧計１１
は、図５に示す様に、被験者の脈波を検出する２つのセ
ンサ２、一定圧で２つのセンサ２を保持するカフ１２、
センサ信号（アナログ信号）をデジタル信号に変換する
Ａ／Ｄ変換器３、体動による外乱を分離する外乱分離装
置１３（後述する）、別途測定した被験者の血圧値を基
準血圧として入力・記憶する入力・記憶手段１４、外乱
を含まない心拍に同期した脈波信号と基準血圧とを対応
させるための係数を求める係数計算手段１５、算出され
た係数と外乱を含まない心拍に同期した脈波信号から連
続的に血圧値を算出する血圧値計算手段１６、算出され
た連続血圧を記録する記録装置１７、算出された連続血
圧を外部へ出力する出力装置１８から成る。
【００４１】センサ２は、例えば、上記の生体信号検出
装置１で説明した光学式脈波センサであり、カフ１２に
よって被験者の腕等に密着して保持される。なお、セン
サ出力は、増幅器によって増幅されても良く、センサ２
及び増幅器とも良好な線形性を有しているものとする。
また、センサ出力は、互いに統一要件を有していること
が、外乱分離処理を行う上で望ましい。即ち、複数のセ
ンサ２の位置を離し過ぎると、異なる外乱要因が一部の
センサ出力にのみ支配的に現れ、血圧算出精度が低下す
る。極端な例としては、２つのセンサ２をそれぞれ右腕
先と左腕先に配置して、右腕のみを軽く振った場合、右
腕に付けたセンサ２には体動の影響が大きく現れるが、
左腕に付けたセンサ２には体動の影響は現れにくい。従
って、同一外乱要因が信号に重畳し易い同一血管系での
脈波測定が可能となるように個々のセンサ２を配置する
ことが望ましい。
装置１で説明した光学式脈波センサであり、カフ１２に
よって被験者の腕等に密着して保持される。なお、セン
サ出力は、増幅器によって増幅されても良く、センサ２
及び増幅器とも良好な線形性を有しているものとする。
また、センサ出力は、互いに統一要件を有していること
が、外乱分離処理を行う上で望ましい。即ち、複数のセ
ンサ２の位置を離し過ぎると、異なる外乱要因が一部の
センサ出力にのみ支配的に現れ、血圧算出精度が低下す
る。極端な例としては、２つのセンサ２をそれぞれ右腕
先と左腕先に配置して、右腕のみを軽く振った場合、右
腕に付けたセンサ２には体動の影響が大きく現れるが、
左腕に付けたセンサ２には体動の影響は現れにくい。従
って、同一外乱要因が信号に重畳し易い同一血管系での
脈波測定が可能となるように個々のセンサ２を配置する
ことが望ましい。
【００４２】例えば、２つのセンサ２を腕の長さ方向に
互いの中心を３cm離して配置し、軽く腕を振った場合、
互いの位相にずれが生じるが、体動による外乱もきれい
に一致した波形が得られる。即ち、波形形状に関して、
統一性を得ることができる。また、２つのセンサ２を腕
の周方向に互いの中心を３cm離して配置し、軽く腕を振
った場合、互いの位相はきれいに一致するが、わずかと
は言え外乱の振幅が異なる。この場合、互いの位相に関
して、統一性を得るとができる。カフ１２は、２つのセ
ンサ２を被験者の拡張期圧より低い一定圧（例えば２０
mmHg）で保持することができる。
互いの中心を３cm離して配置し、軽く腕を振った場合、
互いの位相にずれが生じるが、体動による外乱もきれい
に一致した波形が得られる。即ち、波形形状に関して、
統一性を得ることができる。また、２つのセンサ２を腕
の周方向に互いの中心を３cm離して配置し、軽く腕を振
った場合、互いの位相はきれいに一致するが、わずかと
は言え外乱の振幅が異なる。この場合、互いの位相に関
して、統一性を得るとができる。カフ１２は、２つのセ
ンサ２を被験者の拡張期圧より低い一定圧（例えば２０
mmHg）で保持することができる。
【００４３】Ａ／Ｄ変換器３は、２つのセンサ２から信
号を入力し、例えば１００Ｈｚのサンプリング周波数で
アナログ信号をデジタル信号ｘ1(n)、ｘ2(n)、(n=0,1,
2,,,,) に変換して外乱分離装置１３へ送る。外乱分離
装置１３は、Ａ／Ｄ変換器３でデジタル変換された信号
ｘ1(n)、ｘ2(n)から体動による外乱ｙ2(n)｛≒ｓ2(n)｝
を分離して、外乱を含まない心臓の拍動に基づく脈波信
号ｙ1(n)｛≒ｓ1(n)｝を回復する装置で、図６に示す様
に、分離マトリクス計算手段４と信号分離計算手段５よ
り構成される。この分離マトリクス計算手段４と信号分
離計算手段５の機能及び作用は、上述した生体信号検出
装置１の実施例と同じであり、その説明は省略する。
号を入力し、例えば１００Ｈｚのサンプリング周波数で
アナログ信号をデジタル信号ｘ1(n)、ｘ2(n)、(n=0,1,
2,,,,) に変換して外乱分離装置１３へ送る。外乱分離
装置１３は、Ａ／Ｄ変換器３でデジタル変換された信号
ｘ1(n)、ｘ2(n)から体動による外乱ｙ2(n)｛≒ｓ2(n)｝
を分離して、外乱を含まない心臓の拍動に基づく脈波信
号ｙ1(n)｛≒ｓ1(n)｝を回復する装置で、図６に示す様
に、分離マトリクス計算手段４と信号分離計算手段５よ
り構成される。この分離マトリクス計算手段４と信号分
離計算手段５の機能及び作用は、上述した生体信号検出
装置１の実施例と同じであり、その説明は省略する。
【００４４】（非観血血圧計１１の作用）外乱分離装置
１３により回復された心臓の拍動に基づく脈波信号ｙ1
(n)を基に、以下の経緯で連続血圧値が算出される。ま
ず、係数計算手段１５により、分離後の脈波信号ｙ1(n)
の例えば１０拍分をサンプリングして加算平均し、基準
脈波Ａを合成する。この基準脈波Ａからピーク値Ａpea
k、ボトム値Ａbtm 、直流成分を平均して平均値Ａmean
を求める（図７参照）。また、別途、被験者の収縮期圧
ＳＹＳ、拡張期圧ＤＩＡ、及び平均動脈圧ＭＥＡＮＰを
測定し、校正用の基準値（基準血圧値）として入力・記
憶する。
１３により回復された心臓の拍動に基づく脈波信号ｙ1
(n)を基に、以下の経緯で連続血圧値が算出される。ま
ず、係数計算手段１５により、分離後の脈波信号ｙ1(n)
の例えば１０拍分をサンプリングして加算平均し、基準
脈波Ａを合成する。この基準脈波Ａからピーク値Ａpea
k、ボトム値Ａbtm 、直流成分を平均して平均値Ａmean
を求める（図７参照）。また、別途、被験者の収縮期圧
ＳＹＳ、拡張期圧ＤＩＡ、及び平均動脈圧ＭＥＡＮＰを
測定し、校正用の基準値（基準血圧値）として入力・記
憶する。
【００４５】次に、ピーク値Ａpeakとボトム値Ａbtm と
の差と、収縮期圧ＳＹＳと拡張期圧ＤＩＡとの差を対応
づける係数Ｃを求める。以降、血圧値計算手段１６によ
って、時々刻々得られる体動の分離された現在時点の脈
波のデジタル値ｙ1(t)から、以下の式を用いて連続血圧
値ＢＰを算出することができる。 ＢＰ＝｛ｙ1(t)−Ａmean｝×Ｃ＋ＭＥＡＮＰ 算出された連続血圧値ＢＰは、記録装置１７に記録さ
れ、必要に応じて出力装置１８から外部へ出力される。
の差と、収縮期圧ＳＹＳと拡張期圧ＤＩＡとの差を対応
づける係数Ｃを求める。以降、血圧値計算手段１６によ
って、時々刻々得られる体動の分離された現在時点の脈
波のデジタル値ｙ1(t)から、以下の式を用いて連続血圧
値ＢＰを算出することができる。 ＢＰ＝｛ｙ1(t)−Ａmean｝×Ｃ＋ＭＥＡＮＰ 算出された連続血圧値ＢＰは、記録装置１７に記録さ
れ、必要に応じて出力装置１８から外部へ出力される。
【００４６】（第２実施例）図８は非観血血圧計１１の
構成を示すブロック図である。 （非観血血圧計１１の構成）本実施例の非観血血圧計１
１は、図８に示す様に、被験者の脈波を検出する２つの
センサ２、一定圧で２つのセンサ２を保持するカフ１
２、センサ信号（アナログ信号）をデジタル信号に変換
するＡ／Ｄ変換器３、体動による外乱を分離する外乱分
離装置１３、別途測定した被験者の血圧値（収縮期圧及
び拡張期圧）と脈波伝播時間（外乱を含まない心拍に同
期した脈動成分の伝播時間）とを校正用の基準値として
入力・記憶する入力・記憶手段１９、その基準値（血圧
値と脈波伝播時間）から両者を対応させる係数を算出す
る係数計算手段２０、外乱が分離された心臓の拍動に基
づく脈波信号から脈波伝播時間を求める脈波伝播時間計
算手段２１、算出された脈波伝播時間と前記係数とから
連続的に血圧値を算出する血圧値計算手段２２、算出さ
れた連続血圧を記録する記録装置１７、算出された連続
血圧を外部へ出力する出力装置１８から成る。
構成を示すブロック図である。 （非観血血圧計１１の構成）本実施例の非観血血圧計１
１は、図８に示す様に、被験者の脈波を検出する２つの
センサ２、一定圧で２つのセンサ２を保持するカフ１
２、センサ信号（アナログ信号）をデジタル信号に変換
するＡ／Ｄ変換器３、体動による外乱を分離する外乱分
離装置１３、別途測定した被験者の血圧値（収縮期圧及
び拡張期圧）と脈波伝播時間（外乱を含まない心拍に同
期した脈動成分の伝播時間）とを校正用の基準値として
入力・記憶する入力・記憶手段１９、その基準値（血圧
値と脈波伝播時間）から両者を対応させる係数を算出す
る係数計算手段２０、外乱が分離された心臓の拍動に基
づく脈波信号から脈波伝播時間を求める脈波伝播時間計
算手段２１、算出された脈波伝播時間と前記係数とから
連続的に血圧値を算出する血圧値計算手段２２、算出さ
れた連続血圧を記録する記録装置１７、算出された連続
血圧を外部へ出力する出力装置１８から成る。
【００４７】（非観血血圧計１１の作用）外乱分離装置
１３までの作用は、第１実施例と同じである。外乱分離
装置１３により回復された心臓の拍動に基づく脈波信号
ｙ1(n)を基に、以下の経緯で連続血圧値が算出される。
まず、別途、例えば安静時と運動負荷時に、それぞれに
おける脈波伝播時間と収縮期圧及び拡張期圧を測定し、
これらの測定値を校正用の基準値として入力・記憶手段
１９に入力する。
１３までの作用は、第１実施例と同じである。外乱分離
装置１３により回復された心臓の拍動に基づく脈波信号
ｙ1(n)を基に、以下の経緯で連続血圧値が算出される。
まず、別途、例えば安静時と運動負荷時に、それぞれに
おける脈波伝播時間と収縮期圧及び拡張期圧を測定し、
これらの測定値を校正用の基準値として入力・記憶手段
１９に入力する。
【００４８】次に、係数計算手段２０により、脈波伝播
時間から血圧値を算出するために必要な係数を求める。
即ち、それぞの基準値を測定した時の血圧値と脈波伝播
時間をそれぞれＰ1 、Ｐ2 、Ｔ1 、Ｔ2 、被験者によっ
て異なる固有の定数をα、βとすれば、次式が成立す
る。 Ｐ1 ＝α×Ｔ1 ＋β Ｐ2 ＝α×Ｔ2 ＋β
時間から血圧値を算出するために必要な係数を求める。
即ち、それぞの基準値を測定した時の血圧値と脈波伝播
時間をそれぞれＰ1 、Ｐ2 、Ｔ1 、Ｔ2 、被験者によっ
て異なる固有の定数をα、βとすれば、次式が成立す
る。 Ｐ1 ＝α×Ｔ1 ＋β Ｐ2 ＝α×Ｔ2 ＋β
【００４９】与式からα、βを求めてしまえば、以降、
時々刻々得られる体動の分離された脈波波形から、脈波
伝播時間計算手段２１により脈波伝播時間Ｔを求め、血
圧値計算手段２２により連続血圧値ＢＰを算出すること
ができる。 ＢＰ＝α×Ｔ＋β 算出された連続血圧値ＢＰは、記録装置１７に記録さ
れ、必要に応じて出力装置１８から外部へ出力される。
なお、脈波伝播時間を求めるには、例えば心電電極で計
測される心電波形のＲ波を基に脈波の特異点の遅れ時間
を求めれば良い（図９参照）。
時々刻々得られる体動の分離された脈波波形から、脈波
伝播時間計算手段２１により脈波伝播時間Ｔを求め、血
圧値計算手段２２により連続血圧値ＢＰを算出すること
ができる。 ＢＰ＝α×Ｔ＋β 算出された連続血圧値ＢＰは、記録装置１７に記録さ
れ、必要に応じて出力装置１８から外部へ出力される。
なお、脈波伝播時間を求めるには、例えば心電電極で計
測される心電波形のＲ波を基に脈波の特異点の遅れ時間
を求めれば良い（図９参照）。
【００５０】（実施例の効果）第１実施例及び第２実施
例に記載した非観血血圧計１１では、センサ２で検出さ
れた脈波情報から体動による外乱を分離して心拍に同期
した脈動成分を取り出すことができるので、頻繁に体動
が発生する自由行動下でも連続的に血圧測定を行うこと
ができる。また、カフ１２による締め付けを繰り返す必
要がないため、被験者に大きな負担を与えず、連続的な
血圧測定を無侵襲で行うことができる。更に、カフ圧を
制御するために必要なコンプレッサ等の機械装置や、こ
れを駆動するための大容量バッテリー等の大きく重い部
品を必要としないため、小型軽量化が容易であり、携帯
可能な血圧計として汲み上げることもできる。
例に記載した非観血血圧計１１では、センサ２で検出さ
れた脈波情報から体動による外乱を分離して心拍に同期
した脈動成分を取り出すことができるので、頻繁に体動
が発生する自由行動下でも連続的に血圧測定を行うこと
ができる。また、カフ１２による締め付けを繰り返す必
要がないため、被験者に大きな負担を与えず、連続的な
血圧測定を無侵襲で行うことができる。更に、カフ圧を
制御するために必要なコンプレッサ等の機械装置や、こ
れを駆動するための大容量バッテリー等の大きく重い部
品を必要としないため、小型軽量化が容易であり、携帯
可能な血圧計として汲み上げることもできる。
【００５１】校正用の基準値として被験者の血圧値を別
途測定する時は、センサ２を保持しているカフ１２を使
用して測定できることは言うまでもない。但し、上記の
実施例では、カフ１２によってセンサ２を保持する例を
示しているが、必ずしもカフ１２を使用する必要はな
く、被験者の拡張期圧以下の一定圧でセンサ２を保持で
きるものであれば、例えばゴム製のバンドや保持力を調
節できるベルト等を使用しても良い。
途測定する時は、センサ２を保持しているカフ１２を使
用して測定できることは言うまでもない。但し、上記の
実施例では、カフ１２によってセンサ２を保持する例を
示しているが、必ずしもカフ１２を使用する必要はな
く、被験者の拡張期圧以下の一定圧でセンサ２を保持で
きるものであれば、例えばゴム製のバンドや保持力を調
節できるベルト等を使用しても良い。
【図１】生体信号検出装置の構成を示すブロック図であ
る（第１実施例）。
る（第１実施例）。
【図２】生体信号検出装置の構成を示すブロック図であ
る（第２実施例）。
る（第２実施例）。
【図３】生体信号検出装置の構成を示すブロック図であ
る（第３実施例）。
る（第３実施例）。
【図４】生体信号検出装置の構成を示すブロック図であ
る（第４実施例）。
る（第４実施例）。
【図５】非観血血圧計の構成を示すブロック図である
（第１実施例）。
（第１実施例）。
【図６】外乱分離装置の構成を示すブロック図である。
【図７】第１実施例における血圧測定の方法を示す説明
図である。
図である。
【図８】非観血血圧計の構成を示すブロック図である
（第２実施例）。
（第２実施例）。
【図９】心電波形と脈波波形との相関を表すグラフであ
る。
る。
１ 生体信号検出装置 ２ センサ ３ Ａ／Ｄ変換器（Ａ／Ｄ変換手段） ４ 分離マトリクス計算手段 ５ 信号分離計算手段 ６ 外部出力手段（信号出力手段） ７ モード制御手段（体動判定手段） ８ 信号確率密度関数推定手段（確率密度関数推定手
段） ９ バイパス信号回路（信号バイパス手段） １０ 体動判定手段 １１ 非観血血圧計 １２ カフ １４ 入力・記憶手段（第１実施例） １５ 係数計算手段（第１実施例） １６ 血圧値計算手段（第１実施例） １９ 入力・記憶手段（第２実施例） ２０ 係数計算手段（第２実施例） ２１ 脈波伝播時間計算手段（第２実施例） ２２ 血圧値計算手段（第２実施例）
段） ９ バイパス信号回路（信号バイパス手段） １０ 体動判定手段 １１ 非観血血圧計 １２ カフ １４ 入力・記憶手段（第１実施例） １５ 係数計算手段（第１実施例） １６ 血圧値計算手段（第１実施例） １９ 入力・記憶手段（第２実施例） ２０ 係数計算手段（第２実施例） ２１ 脈波伝播時間計算手段（第２実施例） ２２ 血圧値計算手段（第２実施例）
Claims (16)
- 【請求項１】生体情報を検出する２つ以上のセンサから
前記生体情報をアナログ信号として入力し、そのアナロ
グ信号をデジタル信号に変換するＡ／Ｄ変換手段と、 前記センサで検出された生体情報に体動による外乱が含
まれており、その外乱を含んだ生体情報が前記センサで
検出されるまでに振幅のみ変化して、体動に伴う外乱と
外乱を含んでいない生体情報とが混ざっていると見做さ
れる場合に、前記外乱を含んでいない生体情報と前記外
乱との統計的独立性を利用して、前記外乱を含んだ生体
情報から外乱を分離するための１つの分離マトリクスを
計算する分離マトリクス計算手段と、 前記Ａ／Ｄ変換手段で変換されたデジタル信号に前記分
離マトリクスを作用させて、前記デジタル信号を、外乱
を含まない生体信号と外乱を示す体動信号とに分離する
信号分離計算手段とを備える生体信号検出装置。 - 【請求項２】生体情報を検出する２つ以上のセンサから
前記生体情報をアナログ信号として入力し、そのアナロ
グ信号をデジタル信号に変換するＡ／Ｄ変換手段と、 前記センサで検出された生体情報に体動による外乱が含
まれており、その外乱を含んだ生体情報が前記センサで
検出されるまでに周波数と位相が線形的に変化して、体
動に伴う外乱と外乱を含んでいない生体情報とが混ざっ
ていると見做される場合に、前記外乱を含んでいない生
体情報と前記外乱との統計的独立性を利用して、前記外
乱を含んだ生体情報から外乱を分離するための２つ以上
の分離マトリクスを計算する分離マトリクス計算手段
と、 前記Ａ／Ｄ変換手段で変換されたデジタル信号に前記２
つ以上の分離マトリクスを作用させて、前記デジタル信
号を、外乱を含まない生体信号と外乱を示す体動信号と
に分離する信号分離計算手段とを備える生体信号検出装
置。 - 【請求項３】生体情報を検出する２つ以上のセンサから
前記生体情報をアナログ信号として入力し、そのアナロ
グ信号をデジタル信号に変換するＡ／Ｄ変換手段と、 前記センサから体動による外乱を含んでいない生体情報
を入力し、その生体情報の確率密度関数を推定する確率
密度関数推定手段と、 前記センサで検出された生体情報に体動による外乱が含
まれており、その外乱を含んだ生体情報が前記センサで
検出されるまでに振幅のみ変化して、体動に伴う外乱と
外乱を含んでいない生体情報とが混ざっていると見做さ
れる場合に、前記確率密度関数を用い、前記外乱を含ん
でいない生体情報と前記外乱との統計的独立性を利用し
て、前記外乱を含んだ生体情報から外乱を分離するため
の１つの分離マトリクスを計算する分離マトリクス計算
手段と、 前記Ａ／Ｄ変換手段で変換されたデジタル信号に前記分
離マトリクスを作用させて、前記デジタル信号を、外乱
を含まない生体信号と外乱を示す体動信号とに分離する
信号分離計算手段とを備える生体信号検出装置。 - 【請求項４】生体情報を検出する２つ以上のセンサから
前記生体情報をアナログ信号として入力し、そのアナロ
グ信号をデジタル信号に変換するＡ／Ｄ変換手段と、 前記センサから体動による外乱を含んでいない生体情報
を入力し、その生体情報の確率密度関数を推定する確率
密度関数推定手段と、 前記センサで検出された生体情報に体動による外乱が含
まれており、その外乱を含んだ生体情報が前記センサで
検出されるまでに周波数と位相が線形的に変化して、体
動に伴う外乱と外乱を含んでいない生体情報とが混ざっ
ていると見做される場合に、前記確率密度関数を用い、
前記外乱を含んでいない生体情報と前記外乱との統計的
独立性を利用して、前記外乱を含んだ生体情報から外乱
を分離するための２つ以上の分離マトリクスを計算する
分離マトリクス計算手段と、 前記Ａ／Ｄ変換手段で変換されたデジタル信号に前記２
つ以上の分離マトリクスを作用させて、前記デジタル信
号を、外乱を含まない生体信号と外乱を示す体動信号と
に分離する信号分離計算手段とを備える生体信号検出装
置。 - 【請求項５】請求項１〜４に記載した何れかの生体信号
検出装置において、 前記信号分離計算手段で分離された前記生体信号を外部
に出力する信号出力手段を具備していることを特徴とす
る生体信号検出装置。 - 【請求項６】請求項１〜５に記載した何れかの生体信号
検出装置において、 前記センサから入力した生体情報に体動による外乱が含
まれているか否かを判定する体動判定手段と、 この体動判定手段で体動による外乱が含まれていないと
判定された場合に、前記Ａ／Ｄ変換手段で変換されたデ
ジタル信号を前記分離マトリクス計算手段及び前記信号
分離計算手段を経ずに出力する信号バイパス手段とを備
えていることを特徴とする生体信号検出装置。 - 【請求項７】請求項３〜５に記載した何れかの生体信号
検出装置において、 前記センサから入力した生体情報に体動による外乱が含
まれているか否かを判定する体動判定手段を備え、 前記確率密度関数推定手段は、前記体動判定手段で体動
による外乱が含まれていないと判定された場合に、前記
センサから入力した生体情報の確率密度関数を推定する
ことを特徴とする生体信号検出装置。 - 【請求項８】請求項７に記載した生体信号検出装置にお
いて、 前記体動判定手段で体動による外乱が含まれていないと
判定された場合に、前記Ａ／Ｄ変換手段で変換されたデ
ジタル信号を前記分離マトリクス計算手段及び前記信号
分離計算手段を経ずに出力する信号バイパス手段を備え
ていることを特徴とする生体信号検出装置。 - 【請求項９】請求項１〜８に記載した生体信号検出装置
において、 前記生体情報は、被験者の脈波であることを特徴とする
生体信号検出装置。 - 【請求項１０】被験者の身体表面に保持され、心拍に同
期して得られる脈動成分を検出してアナログ信号として
出力する少なくとも２つ以上のセンサと、 前記センサからアナログ信号を入力し、そのアナログ信
号をデジタル信号に変換するＡ／Ｄ変換手段と、 前記センサで検出された脈動成分に体動による外乱が含
まれており、その外乱を含んだ脈動成分が前記センサで
検出されるまでに周波数と位相が線形的に変化して、体
動に伴う外乱と外乱を含んでいない心拍に同期した脈動
成分とが混ざっていると見做される場合に、前記外乱を
含んでいない脈動成分と前記外乱との統計的独立性を利
用して、前記外乱を含んだ脈動成分から外乱を分離する
ための複数の分離マトリクスを計算する分離マトリクス
計算手段と、 前記Ａ／Ｄ変換手段で変換されたデジタル信号に前記複
数の分離マトリクスを作用させて、前記デジタル信号
を、外乱を含まない心拍に同期した脈動成分と外乱を示
す体動成分とに分離する信号分離計算手段と、 この信号分離計算手段から取り込まれる外乱を含まない
心拍に同期した脈動成分を元に、連続的に血圧値を計算
する血圧値計算手段とを有することを特徴とする非観血
血圧計。 - 【請求項１１】請求項１０に記載した非観血血圧計にお
いて、 別途測定された前記被験者の血圧値を基準血圧として入
力及び記憶する入力・記憶手段と、 前記信号分離計算手段によって得られた外乱を含まない
心拍に同期した脈動成分と前記基準血圧とを対応させる
ための係数を求める係数計算手段とを有し、 前記血圧値計算手段は、前記信号分離計算手段から取り
込まれる外乱を含まない心拍に同期した脈動成分と前記
係数計算手段によって算出された係数とから連続的に血
圧値を計算することを特徴とする非観血血圧計。 - 【請求項１２】請求項１０に記載した非観血血圧計にお
いて、 別途測定された前記被験者の血圧値及び外乱を含まない
心拍に同期した脈動成分の伝播時間をそれぞれ基準血圧
及び基準脈波伝播時間として入力及び記憶する入力・記
憶手段と、 前記基準血圧と前記基準脈波伝播時間とを対応させるた
めの係数を求める係数計算手段と、 前記信号分離計算手段によって得られた外乱を含まない
心拍に同期した脈動成分の伝播時間を算出する脈波伝播
時間計算手段と、 前記血圧値計算手段は、前記脈波伝播時間計算手段から
取り込まれる脈動成分の伝播時間と前記係数計算手段に
よって算出された係数とから連続的に血圧値を計算する
ことを特徴とする非観血血圧計。 - 【請求項１３】請求項１０〜１２に記載した非観血血圧
計において、 前記被験者の身体に装着されるカフを有し、 このカフにより前記センサを前記被験者の拡張期圧以下
の圧力で保持することを特徴とする非観血血圧計。 - 【請求項１４】請求項１０〜１３に記載した非観血血圧
計において、 前記２つ以上のセンサは、心拍に同期して得られる脈動
成分を同一血管系上から検出できる位置に配置されるこ
とを特徴とする非観血血圧計。 - 【請求項１５】請求項１４に記載した非観血血圧計にお
いて、 前記２つ以上のセンサは、前記同一血管系上に沿って直
線的に配置されることを特徴とする非観血血圧計。 - 【請求項１６】請求項１４に記載した非観血血圧計にお
いて、 前記２つ以上のセンサは、前記同一血管系と直交する周
回方向に配置されることを特徴とする非観血血圧計。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000270294A JP4415466B2 (ja) | 1999-12-27 | 2000-09-06 | 生体信号検出装置及び非観血血圧計 |
| US09/745,839 US6620104B2 (en) | 1999-12-27 | 2000-12-26 | Vital signal detecting apparatus and noninvasive blood pressure meter |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP36875199 | 1999-12-27 | ||
| JP11-368751 | 1999-12-27 | ||
| JP2000270294A JP4415466B2 (ja) | 1999-12-27 | 2000-09-06 | 生体信号検出装置及び非観血血圧計 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001245859A true JP2001245859A (ja) | 2001-09-11 |
| JP4415466B2 JP4415466B2 (ja) | 2010-02-17 |
Family
ID=26582029
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000270294A Expired - Fee Related JP4415466B2 (ja) | 1999-12-27 | 2000-09-06 | 生体信号検出装置及び非観血血圧計 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6620104B2 (ja) |
| JP (1) | JP4415466B2 (ja) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006149882A (ja) * | 2004-11-30 | 2006-06-15 | Sumitomo Osaka Cement Co Ltd | 状態解析装置及びソフトウエアプログラム |
| JP2007167369A (ja) * | 2005-12-22 | 2007-07-05 | Nippon Telegr & Teleph Corp <Ntt> | 生体センサ |
| JP2012500081A (ja) * | 2008-08-19 | 2012-01-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 患者血圧モニタリング |
| JP2013027550A (ja) * | 2011-07-28 | 2013-02-07 | Seiko Epson Corp | バイタルサイン計測装置、バイタルサイン計測プログラム、および記録媒体 |
| JP2014000105A (ja) * | 2012-06-15 | 2014-01-09 | Univ Of Tokyo | 非侵襲的連続血圧モニタリング方法及び装置 |
| JP2016152918A (ja) * | 2014-02-27 | 2016-08-25 | 京セラ株式会社 | センサ装置およびセンサ装置の駆動方法 |
| JP2017104326A (ja) * | 2015-12-10 | 2017-06-15 | 日本電信電話株式会社 | 脈波伝播時間変化推定方法、脈波伝播時間変化推定装置、及びプログラム |
| JP6257015B1 (ja) * | 2017-08-02 | 2018-01-10 | 株式会社Ａｒｂｌｅｔ | 開発支援サーバ、開発支援方法及び開発支援プログラム |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SE0102919D0 (sv) * | 2001-08-31 | 2001-08-31 | St Jude Medical | Medical device |
| US6954494B2 (en) * | 2001-10-25 | 2005-10-11 | Siemens Corporate Research, Inc. | Online blind source separation |
| JP4620959B2 (ja) * | 2004-03-26 | 2011-01-26 | キヤノン株式会社 | 生体情報モニタ装置 |
| US7434724B2 (en) * | 2006-12-22 | 2008-10-14 | Welch Allyn, Inc. | Dynamic barcode for displaying medical data |
| WO2016145522A1 (en) | 2015-03-13 | 2016-09-22 | Xinde Li | Adaptive methods and systems for detecting signals from interference contaminated data |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CH678690A5 (ja) * | 1989-03-08 | 1991-10-31 | Asulab Sa | |
| JPH0624525B2 (ja) | 1991-07-04 | 1994-04-06 | 日本光電工業株式会社 | 連続型非観血血圧測定装置 |
| JPH05154136A (ja) | 1991-11-29 | 1993-06-22 | Shimadzu Corp | 光生体計測装置 |
| JPH05184547A (ja) | 1992-01-10 | 1993-07-27 | Omron Corp | 電子血圧計 |
| JPH05207978A (ja) | 1992-01-24 | 1993-08-20 | Nec Corp | 無拘束運動負荷時血圧変動測定装置 |
| US5289823A (en) * | 1992-05-12 | 1994-03-01 | Colin Electronics Co., Ltd. | Non-invasive aortic blood flow sensor and method for non-invasively measuring aortic blood flow |
| US5795300A (en) * | 1994-06-01 | 1998-08-18 | Advanced Body Metrics Corporation | Heart pulse monitor |
| JP3318727B2 (ja) | 1994-06-06 | 2002-08-26 | 日本光電工業株式会社 | 脈波伝播時間方式血圧計 |
| US6331162B1 (en) * | 1999-02-01 | 2001-12-18 | Gary F. Mitchell | Pulse wave velocity measuring device |
| JP2001170016A (ja) * | 1999-12-16 | 2001-06-26 | Seiko Instruments Inc | 脈波検出装置 |
-
2000
- 2000-09-06 JP JP2000270294A patent/JP4415466B2/ja not_active Expired - Fee Related
- 2000-12-26 US US09/745,839 patent/US6620104B2/en not_active Expired - Fee Related
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006149882A (ja) * | 2004-11-30 | 2006-06-15 | Sumitomo Osaka Cement Co Ltd | 状態解析装置及びソフトウエアプログラム |
| JP4588423B2 (ja) * | 2004-11-30 | 2010-12-01 | 住友大阪セメント株式会社 | 状態解析装置及びソフトウエアプログラム |
| JP2007167369A (ja) * | 2005-12-22 | 2007-07-05 | Nippon Telegr & Teleph Corp <Ntt> | 生体センサ |
| JP2012500081A (ja) * | 2008-08-19 | 2012-01-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 患者血圧モニタリング |
| JP2013027550A (ja) * | 2011-07-28 | 2013-02-07 | Seiko Epson Corp | バイタルサイン計測装置、バイタルサイン計測プログラム、および記録媒体 |
| JP2014000105A (ja) * | 2012-06-15 | 2014-01-09 | Univ Of Tokyo | 非侵襲的連続血圧モニタリング方法及び装置 |
| JP2016152918A (ja) * | 2014-02-27 | 2016-08-25 | 京セラ株式会社 | センサ装置およびセンサ装置の駆動方法 |
| JP2017104326A (ja) * | 2015-12-10 | 2017-06-15 | 日本電信電話株式会社 | 脈波伝播時間変化推定方法、脈波伝播時間変化推定装置、及びプログラム |
| JP6257015B1 (ja) * | 2017-08-02 | 2018-01-10 | 株式会社Ａｒｂｌｅｔ | 開発支援サーバ、開発支援方法及び開発支援プログラム |
| JP2019028880A (ja) * | 2017-08-02 | 2019-02-21 | 株式会社Ａｒｂｌｅｔ | 開発支援サーバ、開発支援方法及び開発支援プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| US20010016691A1 (en) | 2001-08-23 |
| US6620104B2 (en) | 2003-09-16 |
| JP4415466B2 (ja) | 2010-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2764702B2 (ja) | 血圧監視装置 | |
| JP4855721B2 (ja) | 血圧測定装置 | |
| US6083171A (en) | Blood pressure monitoring apparatus | |
| US7137955B2 (en) | Methods and systems for distal recording of phonocardiographic signals | |
| US5682898A (en) | Respiration rate measuring apparatus | |
| US9131859B2 (en) | Blood pressure measurement apparatus, recording medium that records blood pressure derivation program, and blood pressure derivation method | |
| JP2001245859A (ja) | 生体信号検出装置及び非観血血圧計 | |
| US6582374B2 (en) | Automatic blood-pressure measuring apparatus | |
| US8282567B2 (en) | Method and system for determination of pulse rate | |
| JPH07327940A (ja) | 脈波伝播時間方式血圧計 | |
| WO2018043692A1 (ja) | 血圧測定装置、血圧測定方法及び血圧測定プログラムを記録した記録媒体 | |
| WO2004103185A1 (ja) | 超音波診断装置 | |
| JPH07136136A (ja) | 連続血圧モニタ装置 | |
| WO2006024871A1 (en) | Methods and apparatus for the measurement of blood pressure | |
| US6129677A (en) | Non-invasive heart monitoring apparatus and method | |
| CN104274165B (zh) | 确定设备和确定方法 | |
| JP3218890B2 (ja) | 血圧計測装置 | |
| JP3023505B2 (ja) | 血圧監視装置 | |
| JP2002320593A (ja) | 動脈血圧測定方法および動脈血圧測定装置 | |
| JP3023731B2 (ja) | 心音解析装置 | |
| JP3533122B2 (ja) | 脈波モニタ | |
| US11850039B2 (en) | Apparatus and method for determining a change in left ventricular twist of a subject's heart | |
| JPH0866377A (ja) | 多機能血圧計 | |
| KR102347155B1 (ko) | 두 개의 센서를 이용하는 실시간 혈압 측정 디바이스 및 그를 이용한 실시간 혈압 측정 방법 | |
| JP3028393B2 (ja) | 非観血血圧計 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061017 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20091014 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20091104 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20091117 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121204 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131204 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |