(57)【特許請求の範囲】
【請求項１】
所定サンプリング時間でサンプリングした一連の複数の基本脈波を２回微分し一連の複数の加速度脈波を求める加速度脈波検出手段と、
前記加速度脈波検出手段によって検出した加速度脈波から、前記複数の基本脈波の繰返し周期のばらつき程度を示すばらつき指標を算出する指標算出手段と、
前記ばらつき指標に前記繰返し周期のばらつきが加速度脈波加齢値に及ぼす影響度を補正する補正係数を予め対応させて形成した補正データと、
前記加速度脈波検出手段で求めた一連の複数の加速度脈波に対し前記指標算出手段で前記ばらつき指標を求め、求めた前記ばらつき指標に対応する前記補正係数を前記補正データから抽出する補正係数抽出手段と、
前記加速度脈波検出手段で求めた一連の複数の加速度脈波に対し仮の加速度脈波加齢値を算出するとともに前記仮の加速度脈波加齢値を前記補正データから抽出した前記補正係数で補正し加速度脈波加齢値を求める加速度脈波算出手段と、
を備えることを特徴とする脈波検出装置。
【請求項２】
前記ばらつき指標は第１のばらつき指標と第２のばらつき指標とから構成されており、
前記補正データの前記補正係数は前記第１のばらつき指標と前記第２のばらつき指標とを行列として特定され、
前記第１のばらつき指標は前記複数の基本脈波の繰返し周期の平均値である
ことを特徴とする請求項１に記載の脈波検出装置。
【請求項３】
前記第２のばらつき指標は、繰返し周期の標準偏差、平方和、分散、範囲、変動係数のいずれかである
ことを特徴とする請求項２に記載の脈波検出装置。
【請求項４】
前記加速度脈波検出手段で検出した前記一連の複数の加速度脈波の周波数分布を求める周波数分析手段をさらに備え、
前記ばらつき指標は、加速度脈波の周波数分布に係る量である
ことを特徴とする請求項１に記載の脈波検出装置。
【請求項５】
前記周波数分布に係る量は、加速度脈波の周波数分布における低周波帯域面積と高周波帯域面積との面積比である
ことを特徴とする請求項４に記載の脈波検出装置。
【請求項６】
前記周波数分布に係わる量は、加速度脈波の周波数分布における低周波帯域面積をＬＦ、中周波帯域面積をＭＦ、高周波帯域面積をＨＦとしたときに（２ＬＦ＋３ＭＦ＋ＨＦ）／６であることを特徴とする請求項４に記載の脈波検出装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、医療分野や健康管理分野などで医用電子機器として使用される脈波検出装置に関する。
【０００２】
【従来技術】
脈波とは、心臓の収縮により血流が大動脈に押し出されたときに血管内に発生する圧力変化が抹消方向に伝達するときの波動のことである。
【０００３】
従来の脈波検出装置は、光電式や圧電式等の脈波センサから脈波を測定開始から一定の時間サンプリングをして、測定結果として脈波を加速度脈波の加齢値（ＡＰＧｉｎｄｅｘ）などの波形として表示するものであった。
【０００４】
【発明が解決しようとする課題】
心臓の拍動間隔は常に一定であるというわけではなく、変化することが知られている。そして、この心臓の拍動間隔に対応して、本来変化しない脈波も約１０秒前後の周期で変化してしまうことがわかっている。
【０００５】
このために、従来の脈波計は、脈波が周期変化することに関係なくサンプリングされるため、変化する周期のどの帯域の脈波を捕らえるかにより、測定結果が測定のたびに変動しばらつきが生じる、という問題があった。また、測定結果に対するばらつきを小さなものにするためには、サンプリング時間を長く取らなければならないという問題があった。
【０００６】
そこで、本発明の目的は、上記従来技術の有する問題を解消し、短い測定時間で、ばらつきの小さい脈波を高精度に測定可能な脈波検出装置を提供することである。
【０００７】
【課題を解決するための手段】
上記目的を達成するために、本発明の脈波検出装置は、所定サンプリング時間でサンプリングした一連の複数の基本脈波を２回微分し一連の複数の加速度脈波を求める加速度脈波検出手段と、前記加速度脈波検出手段によって検出した加速度脈波から、前記複数の基本脈波の繰返し周期のばらつき程度を示すばらつき指標を算出する指標算出手段と、前記ばらつき指標に前記繰返し周期のばらつきが加速度脈波加齢値に及ぼす影響度を補正する補正係数を予め対応させて形成した補正データと、前記加速度脈波検出手段で求めた一連の複数の加速度脈波に対し前記指標算出手段で前記ばらつき指標を求め、求めた前記ばらつき指標に対応する前記補正係数を前記補正データから抽出する補正係数抽出手段と、前記加速度脈波検出手段で求めた一連の複数の加速度脈波に対し仮の加速度脈波加齢値を算出するとともに前記仮の加速度脈波加齢値を前記補正データから抽出した前記補正係数で補正し加速度脈波加齢値を求める加速度脈波算出手段と、を備えることを特徴とする。
【０００８】
また、前記ばらつき指標は第１のばらつき指標と第２のばらつき指標とから構成されており、前記補正データの前記補正係数は前記第１のばらつき指標と前記第２のばらつき指標とを行列として特定され、前記第１のばらつき指標は前記複数の基本脈波の繰返し周期の平均値であることを特徴とする。
【０００９】
また、前記第２のばらつき指標は、繰返し周期の標準偏差、平方和、分散、範囲、変動係数のいずれかであることを特徴とする。
【００１０】
また、前記加速度脈波検出手段で検出した前記一連の複数の加速度脈波の周波数分布を求める周波数分析手段をさらに備え、前記ばらつき指標は、加速度脈波の周波数分布に係る量であることを特徴とする。
【００１１】
また、前記周波数分布に係る量は、加速度脈波の周波数分布における低周波帯域面積と高周波帯域面積との面積比であることを特徴とする。
【００１２】
上述の発明において、心臓の拍動間隔の変動が加速度脈波加齢値の測定に及ぼす影響を除去するために、加速度脈波の繰返し周期のばらつきの程度を示すばらつき指標を用いる。そして、加速度脈波検出手段によって加速度脈波を検出する毎に加速度脈波の繰返し周期のばらつき指標を求め、このばらつき指標に対応する補正係数を予め形成されている補正データを参照して補正係数抽出手段によって抽出し、加速度脈波算出手段によって仮の加速度脈波加齢値を算出するとともに仮の加速度脈波加齢値を抽出した補正係数で補正する。これによって、心臓の拍動間隔の変動の影響を受けることなく加速度脈波加齢値を測定することができる。
【００１３】
【発明の実施の形態】
以下に図面を参照して、本発明に係る脈波検出装置の実施の形態について説明する。
【００１４】
図２に示すように、脈波検出装置は、装置本体１と、センサ部３と、装置本体１とセンサ部３とを接続する導線４とを有する。センサ部３は、サック状に形成され手の指５に装着して使用される。装置本体１の表面部には、データ入力モードの切り替えや電源投入等のためのスイッチ部６と計測結果等の表示をする表示部７とが設けられている。なお、図２に示す例は、一例であり、これ以外の他の構成でもよく、例えば、センサ部３を耳たぶに装着し装置本体１を手首に装着するようなタイプでもよく、また、センサ部３を血圧計のようなカフで腕に巻くようなタイプでもよい。
【００１５】
図１に、脈波検出装置のブロック図を示す。脈波検出装置は、所定サンプリング時間でサンプリングした一連の複数の基本脈波を２回微分し一連の複数の加速度脈波を求める加速度脈波検出手段１０を備え、加速度脈波検出手段１０は、ＬＥＤ等の発光部とフォトダイオード又はフォトトランジスタ等からなる受光部を備え心臓の拍動による血流量の増減を光の透過量や散乱光量の変化により一連の基本脈波を検出するセンサ部３と、センサ部３からのアナログ信号を増幅する増幅回路２２と、増幅回路２２からの基本脈波のアナログ信号を２回微分して加速度脈波信号に変換する２段に接続された微分回路２３，２４と、微分回路２４からのアナログ信号をデジタル信号に変換するＡ／Ｄ変換回路２５を有する。図８（ａ）に１個の基本脈波１０１の波形の一例を示し、図８（ｂ）に基本脈波１０１を２回微分して得られた加速度脈波１０２の波形を示す。加速度脈波に対し、後述するように基線Ｌから変曲点までの大きさａ，ｂ，ｃ，ｄが算出され、加速度脈波加齢値が算出される。
【００１６】
Ａ／Ｄ変換回路２５では、所定サンプリング時間、例えば数１０秒間に渡ってサンプリングした複数の加速度脈波から構成される図４に示すような加速度脈波信号が得られる。
【００１７】
脈波検出装置は、また、加速度脈波検出手段１０で検出した加速度脈波に対し各種制御処理をするＣＰＵ１１と、データを記憶する記憶部１２とを有する。
【００１８】
ＣＰＵ１１は、指標算出手段２６と、補正係数抽出手段１７と、加速度脈波算出手段１９とを有する。記憶部１２は、計測データ記憶部１３と関連データ記憶部１５とを有する。
【００１９】
指標算出手段２６は、加速度脈波検出手段１０によって検出した加速度脈波から、複数の基本脈波あるいは加速度脈波の繰返し周期のばらつき程度を示すばらつき指標を算出する。ここで、複数の基本脈波あるいは加速度脈波の繰返し周期は図４に示すＲで表される。心臓の拍動間隔は常に一定であるというわけではなく、繰返し周期Ｒは変化することが知られている。繰返し周期の第１のばらつき指標として、繰返し周期の平均値Ｒバーが採用され、また、第２のばらつき指標としては、例えば、繰返し周期の標準偏差σが採用される。繰返し周期の第１のばらつき指標及び第２のばらつき指標のデータは、計測毎に計測データ記憶部１３に記憶される。
【００２０】
また、計測データ記憶部１３には、加速度脈波検出手段１０によって検出された加速度脈波から加速度脈波加齢値を算出するのに要するａ，ｂ，ｃ，ｄ等の計測要素のデータも記憶される。
【００２１】
【数１】
0003781956.tif 000002
具体的には、指標算出手段２６は、サンプリングした複数の加速度脈波からなる加速度脈波信号について、図４に示すように基線Ｌからの各変曲点までの大きさ（ａ１〜ｄ１、ａ２〜ｄ２、ａ３〜ｄ３、ａｎ〜ｄｎ）を求め、各加速度脈波派の同じポイントの変曲点の大きさ同士（ａ１〜ａｎ、ｂ１〜ｂｎ、ｃ１〜ｃｎ、ｄ１〜ｄｎ）を平均した値ａバー、ｂバー、ｃバー、ｄバーを式（３）乃至（６）に従って求め、また、サンプリングした複数の加速度脈波のピーク間隔の平均時間Ｒバーを式（１）に従って求め、サンプリングした複数の加速度脈波の脈波間隔（繰返し周期Ｒ）の標準偏差σを式（２）に従って求める。
【００２２】
関連データ記憶部１５には、図３に示すように補正データが記憶されている。補正データは、第１のばらつき指標と第２のばらつき指標に、繰返し周期のばらつきが加速度脈波加齢値に及ぼす影響度を補正する補正係数を予め対応させて形成されている。第１のばらつき指標としては、例えば繰返し周期の平均値（Ｒバーｎ）が取られ、第２のばらつき指標としては例えば繰返し周期の標準偏差（σｍ）が取られる。
【００２３】
繰返し周期のばらつきが加速度脈波加齢値に及ぼす影響度は、予め実験データあるいは経験データを収集して得られ、繰返し周期にばらつきがないとした場合のあるべき加速度脈波加齢値を想定して補正データが定められる。図３において、行部は、第２のばらつき指標の一例としての繰返し周期の標準偏差（σｍ）を示し、列部は第１のばらつき指標である繰返し周期の平均値（Ｒバーｎ）を示し、行列（σｍ、Ｒバーｎ）によって補正係数（Ｓｍｎ）が定まる。
【００２４】
なお、ここでは、第２のばらつき指標として標準偏差（σｍ）をとった場合を示したが、標準偏差に限らず、例えば、繰返し周期の平方和、分散、範囲、変動係数等を採用してもよい。
【００２５】
補正係数抽出手段１７は、指標算出手段２６で求めた繰返し周期の平均値及び標準偏差に対応する補正係数Ｓを、関連データ記憶部１５に記憶された補正データから抽出する。
【００２６】
加速度脈波算出手段１９は、加速度脈波検出手段１０で求めた一連の複数の加速度脈波に対し式（７）に従い、仮の加速度脈波加齢値（ＡＰＧｉｎｄｅｘバー）を算出し、仮の加速度脈波加齢値を補正係数抽出手段１７によって抽出した補正係数Ｓで乗算し加速度脈波加齢値（ＡＰＧｉｎｄｅｘ’）を求める。ここで、仮の加速度脈波加齢値とは、補正データによって補正する前の加速度脈波加齢値をいう。
【００２７】
以上、本実施形態によれば、予め関連データ記憶部１５に補正データを記憶しておき、加速度脈波検出手段１０によって加速度脈波を検出する毎に繰返し周期の平均値及び標準偏差を算出し、この算出結果に対応する補正係数Ｓを補正係数抽出手段１７によって抽出し、加速度脈波算出手段１９によって仮の加速度脈波加齢値（ＡＰＧｉｎｄｅｘバー）を算出するとともに仮の加速度脈波加齢値を補正係数Ｓで補正するようにしたので、心臓の拍動間隔の変動の影響を受けることなく加速度脈波加齢値を測定することができる。
【００２８】
また、心臓の拍動間隔の変動の影響を小さくするためにサンプリング時間を長く取る必要がなくなり、測定時間を短くすることができる。
【００２９】
次に、図５乃至図７を参照して、本発明の第２の実施形態について説明する。
図５に、本実施形態の脈波検出装置のブロック図を示す。図１に示す脈波検出装置と同一部材には同一の符号をつけ説明を省略し、異なる部分についてのみ説明する。
【００３０】
図５に示すように、ＣＰＵ１１には、加速度脈波検出手段１０におけるＡ／Ｄ変換回路２５で検出した加速度脈波を高速フーリェ変換（ＦＦＴ）する周波数分析手段２８を備えている。周波数分析手段２８により、図６に示すような周波数分布データが得られる。
【００３１】
指標算出手段２６は、周波数分析手段２８により得た周波数分布データを図６に示すように、二つの低周波帯域ＬＨ（０．０４〜０．１４Ｈｚ）と高周波帯域ＨＦ（０．１４〜０．５０Ｈｚ）とを設け、低周波帯域ＬＨと高周波帯域ＨＦの面積比ＬＦ／ＨＦを求める。そして、繰返し周期のばらつき指標として面積比ＬＦ／ＨＦを採用する。
【００３２】
なお、繰返し周期のばらつき指標としては、面積比ＬＦ／ＨＦに限らず、繰返し周期のばらつきの指標となるものであれば、他の加速度脈波の周波数分布に係る量であってもよい。例えば、加速度脈波の周波数分布を３つの低周波帯域ＬＨ、中周波帯域ＭＦ及び高周波帯域ＨＦとに分け、各々の面積をＬＦ、ＭＦ及びＨＦとしたときに、（２ＬＦ＋３ＭＦ＋ＨＦ）／６をばらつき指標とし、図７の横軸にとってもよい。
【００３３】
関連データ記憶部１５には、図７に示すように補正データが記憶されている。図７に示す補正データは、横軸に取られた面積比ＬＦ／ＨＦと、縦軸に取られた補正係数Ｓで定まる直線状のグラフによって与えられる。
【００３４】
次に、本実施形態の作用について説明する。
指標算出手段２６は、加速度脈波検出手段１０で検出した加速度脈波に対し面積比ＬＦ／ＨＦを算出する。補正係数抽出手段１７は、指標算出手段２６で算出した面積比ＬＦ／ＨＦに対応する補正係数Ｓを、関連データ記憶部１５に記憶された図７に示す補正データから抽出する。加速度脈波算出手段１９は、加速度脈波検出手段１０で求めた一連の複数の加速度脈波に対し式（７）に従い、仮の加速度脈波加齢値（ＡＰＧｉｎｄｅｘバー）を算出し、仮の加速度脈波加齢値を補正係数抽出手段１７によって抽出した補正係数Ｓで乗算し加速度脈波加齢値（ＡＰＧｉｎｄｅｘ’）を求める。
【００３５】
本実施の形態によれば、心臓の拍動間隔の変動の影響を受けることなく加速度脈波加齢値を測定することができる。また、心臓の拍動間隔の変動の影響を小さくするためにサンプリング時間を長く取る必要がなくなり、測定時間を短くすることができる。
【００３６】
【発明の効果】
以上説明したように、本発明の構成によれば、心臓の拍動間隔の変動の影響を受けることなく高精度な加速度脈波加齢値の測定を短時間で可能にする脈波検出装置を提供することができる。
【図面の簡単な説明】
【図１】本発明に係る脈波検出装置の一実施形態の概略構成を示すブロック図。
【図２】本発明に係る脈波検出装置の装置構成の概略構成を示す図。
【図３】第１のばらつき指標としての繰返し周期の平均値と第２のばらつき指標としての標準偏差との行列で指定される補正データを示す図。
【図４】サンプリングされた一連の複数の加速度脈波を示す波形図。
【図５】本発明に係る脈波検出装置の他の実施形態の概略構成を示すブロック図。
【図６】ＦＦＴ処理によって得られた周波数スペクトル分布を示す図。
【図７】ばらつき指標として、加速度脈波の周波数分布における低周波帯域面積と高周波帯域面積との面積比を採用した場合の補正データを示す図。
【図８】基本脈波の信号波形を示す図（ａ）と、（ａ）に示す信号波形を２回微分して得られた加速度脈波の信号波形を示す図（ｂ）。
【符号の説明】
１  装置本体
３  センサ部
１０  加速度脈波検出手段
１１  ＣＰＵ
１３  計測データ記憶部
１５  関連データ記憶部（補正データ）
１７  補正係数抽出手段
１９  加速度脈波算出手段
２６  指標算出手段
２８  周波数分析手段
【図１】
0003781956.tif 000003
【図２】
0003781956.tif 000004
【図３】
0003781956.tif 000005
【図４】
0003781956.tif 000006
【図５】
0003781956.tif 000007
【図６】
0003781956.tif 000008
【図７】
0003781956.tif 000009
【図８】
0003781956.tif 000010