(57)【特許請求の範囲】
【請求項１】
生体の脈波を時間領域で示す時間領域脈波信号を取得する第１信号取得部と、
前記生体の体動を時間領域で示す時間領域体動信号を取得する第２信号取得部と、
前記時間領域脈波信号と前記時間領域体動信号のそれぞれを周波数領域に変換して、周波数領域脈波信号と周波数領域体動信号をそれぞれ生成し、前記周波数領域脈波信号と周波数領域体動信号に基づき前記生体の脈波の周波数を推定する周波数領域分析部と、
前記時間領域脈波信号と、前記生体の脈波の周波数に基づき、前記生体の生体情報を算出する時間領域分析部と、
を備える測定装置。
【請求項２】
前記周波数領域分析部は、
前記時間領域脈波信号をダウンサンプリングし、ダウンサンプリング後の信号を周波数領域に変換することにより、前記周波数領域脈波信号を取得する第１変換部と、
前記時間領域体動信号をダウンサンプリングし、ダウンサンプリング後の信号を周波数領域に変換することにより、前記周波数領域体動信号を取得する第２変換部と、
前記周波数領域脈波信号を、前記周波数領域体動信号に基づき加工し、周波数領域加工信号を取得する周波数スペクトル演算部と、
前記周波数領域加工信号の周波数領域でのピークを検出する周波数ピーク検出部と、
前記周波数領域でのピークの位置に基づき前記生体の脈波の周波数を推定する周波数推定部と、
を備える
請求項１に記載の測定装置。
【請求項３】
前記時間領域分析部は、前記生体の脈波の周波数に基づき、前記時間領域脈波信号をフィルタリングし、フィルタリング後の信号に基づき、前記生体の生体情報を算出する
請求項１または２に記載の測定装置。
【請求項４】
前記時間領域分析部は、
前記生体の脈波の周波数に基づき、前記時間領域脈波信号をフィルタリングする脈拍成分抽出部と、
前記脈拍成分抽出部によりフィルタリングされた信号の時間領域でのピークを検出する時間ピーク検出部と、
前記時間ピーク検出部により検出された、前記時間領域でのピークの位置に基づき、前記生体の前記生体情報を取得する生体情報取得部
を備える請求項１ないし３のいずれか一項に記載の測定装置。
【請求項５】
前記生体情報は心拍数に関する情報である
請求項１ないし４のいずれか一項に記載の測定装置。
【請求項６】
前記時間領域脈波信号から、予め定められた周波数の信号を除去する第１ノイズ除去部と、
前記時間領域体動信号から、予め定められた周波数の信号を除去する第２ノイズ除去部と、
をさらに備える請求項１ないし５のいずれか一項に記載の測定装置。
【請求項７】
前記時間領域脈波信号から、直流成分の信号を除去する第１直流除去部と、
前記時間領域体動信号から、直流成分の信号を除去する第２直流除去部と、
をさらに備える請求項１ないし６のいずれか一項に記載の測定装置。
【請求項８】
前記時間領域脈波信号を取得する第３〜第ｎ（ｎは３以上の整数）の信号取得部をさらに備え、
前記時間領域分析部は、前記第１信号取得部で取得された前記時間領域脈波信号と、前記生体の脈波の周波数と、前記第３〜第ｎの信号取得部で取得された前記時間領域脈波信号とに基づき、前記生体の生体情報を算出する
請求項１ないし７のいずれか一項に記載の測定装置。
【請求項９】
生体の脈波を時間領域で示す時間領域脈波信号を取得する第１信号取得ステップと、
前記生体の体動を時間領域で示す時間領域体動信号を取得する第２信号取得ステップと、
前記時間領域脈波信号と前記時間領域体動信号のそれぞれを周波数領域に変換して、周波数領域脈波信号と周波数領域体動信号をそれぞれ生成し、前記周波数領域脈波信号と周波数領域体動信号に基づき前記生体の脈波の周波数を推定する周波数分析ステップと、
前記時間領域脈波信号と、前記生体の脈波の周波数に基づき、前記生体の生体情報を算出する時間領域分析ステップと、
をコンピュータが実行する測定方法。
【請求項１０】
生体の脈波を時間領域で示す時間領域脈波信号を取得する第１信号取得ステップと、
前記生体の体動を時間領域で示す時間領域体動信号を取得する第２信号取得ステップと、
前記時間領域脈波信号と前記時間領域体動信号のそれぞれを周波数領域に変換して、周波数領域脈波信号と周波数領域体動信号をそれぞれ生成し、前記周波数領域脈波信号と周波数領域体動信号に基づき前記生体の脈波の周波数を推定する周波数分析ステップと、
前記時間領域脈波信号と、前記生体の脈波の周波数に基づき、前記生体の生体情報を算出する時間領域分析ステップと、
をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、測定装置、測定方法およびプログラムに関する。
【背景技術】
【０００２】
近年、健康志向の高まりを受けて、ウェアラブル端末にて生体情報を測定することが普及しつつある。特に、反射型光電脈波センサを用いた腕時計型の脈拍測定装置は、測定が比較的容易であり、ユーザに不快感を与えないという利点がある。
【０００３】
反射型光電脈波センサは、動脈に照射した光の反射光を利用して、脈拍数を測定する。動脈中のヘモグロビンは光を吸収する性質があるため、反射光は拍動時の血管容積変化によるヘモグロビン量の変化に応じて変動する。この反射光の変動により、脈拍数を測定する。
【０００４】
しかし、反射型光電脈波センサは、センサと生体との接触状態や、生体内の血流の変化により、その測定波形に大きな乱れが生じる事が知られている。そのため、運動時など、被験者の体動が大きい場合は、これらの雑音を除去する処理を行う必要がある。
【０００５】
雑音を除去する処理は、処理演算器の大型化、および消費電力量増加に対応するためのバッテリの大型化につながり、小型化を求めるウェアラブルデバイスにとって問題が生ずる。また、生産コストも高くなる。ゆえに、測定装置の処理量を減らすことが求められる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第２８１６９４４号
【特許文献２】特許第３５８４１４３号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明の実施形態は、測定装置の精度を保ちつつ、処理量を抑えることを目的とする。
【課題を解決するための手段】
【０００８】
本発明の一態様としての測定装置は、生体の脈波信号を取得する第１信号取得部と、前記生体の体動信号を取得する第２信号取得部と、前記脈波信号と前記体動信号のそれぞれを周波数領域に変換して周波数領域信号をそれぞれ生成し、前記周波数領域信号に基づき前記生体の脈波の周波数を推定する周波数分析部と、前記脈波信号と、前記周波数に基づき、前記生体の生体情報を算出する時間領域分析部と、を備える。
【図面の簡単な説明】
【０００９】
【図１】第１の実施形態に係る測定装置の概略構成を示すブロック図。
【図２】脈波の一例を示す図。
【図３】脈波信号と体動信号並びに脈波信号スペクトルと体動信号スペクトルの一例を示す図。
【図４】スペクトルサブトラクションの一例を示す図。
【図５】脈拍成分抽出部による処理前と処理後の脈波信号の一例を示す図。
【図６】第１の実施形態に係る測定装置の処理結果の一例を示す図。
【図７】第１の実施形態に係る測定装置の全体処理の概略フローチャート。
【図８】脈波信号処理部による処理のフローチャート。
【図９】体動信号処理部による処理のフローチャート。
【図１０】周波数領域分析部による処理のフローチャート。
【図１１】時間領域分析部による処理のフローチャート。
【図１２】第２の実施形態に係る測定装置の概略構成を示すブロック図。
【図１３】本実施形態に係る測定装置を実現したハードウェア構成例を示すブロック図。
【発明を実施するための形態】
【００１０】
以下、図面を参照しながら、本発明の実施形態について説明する。
【００１１】
（第１の実施形態）
図１は、第１の実施形態に係る測定装置の概略構成を示すブロック図である。第１の実施形態に係る測定装置は、測定した生体の脈波信号から、生体情報を取得する。測定装置は、生体の一部に装着される事を想定する。
【００１２】
第１の実施形態に係る測定装置は、脈波信号処理部１００と、体動信号処理部２００と、周波数領域分析部３００と、時間領域分析部４００と、出力部５０１と、を備える。
【００１３】
以下、各部について説明する。
【００１４】
脈波信号処理部１００は、生体の脈波信号を図示しないセンサから取得する。脈波信号処理部１００は、脈波信号取得部１０１と、第１直流除去部１０２と、第１ノイズ除去部１０３と、を備える。
【００１５】
体動信号処理部２００は、生体の体動などから生ずる測定装置自体の動きに関する信号を、図示しないセンサから取得する。体動信号処理部２００は、体動信号取得部２０１と、第２直流除去部２０２と、第２ノイズ除去部２０３と、を備える。
【００１６】
周波数領域分析部３００は、脈波信号と体動信号を周波数領域に変換し、前記脈波信号の周波数を算出する。周波数領域分析部３００は、第１周波数領域変換部３０１と、第１周波数領域変換部３０２と、第２周波数領域変換部３１１と、第２周波数領域変換部３１２と、周波数スペクトル演算部３２１と、周波数ピーク検出部３２２と、周波数推定部３２３と、を備える。
【００１７】
時間領域分析部は、脈波信号と周波数に基づき、脈拍数などの生体情報を算出する。時間領域分析部４００は、脈拍成分抽出部４０１と、時間ピーク検出部４０２と、生体情報取得部４０３と、を備える。
【００１８】
以下、各部の詳細について説明する。
【００１９】
脈波信号取得部１０１は、生体の脈波信号を取得する。脈波とは、心臓の収縮により血液が大動脈に押し出されることで発生する血管内の圧力変化が、抹消方向に伝達する時の波動のことである。
【００２０】
なお、以下の説明において、脈波信号は、反射型光電脈波センサにより、拍動時の血管容積変化に伴うヘモグロビン量の変化に基づき、測定されることを想定する。
【００２１】
第１直流除去部１０２は、脈波信号に対し、直流成分の信号を除去する。反射型光電脈波センサを用いた脈波信号には、脈動により交流成分として検出される脈波信号以外に、吸収および反射された光により検出される直流成分が含まれる。直流成分は、所定の周波数帯域成分のみを通過させるハイパスフィルタなどにて、取り除くことができる。
【００２２】
なお、図１において、第１直流除去部１０２は第１ノイズ除去部の手前に図示されているが、後にあってもよい。なお、第１直流除去部１０２は、省略可能である。
【００２３】
第１ノイズ除去部１０３は、脈波信号に対し、予め定められた一定範囲外の信号を除去する。例えば、ヒトの正常な１分あたりの脈拍数は４０から２２０ｂｐｍ程度であるため、脈波の周波数は約０．６６Ｈｚから約３．７Ｈｚの範囲を取る。ゆえに、第１ノイズ除去部１０３は、当該範囲外の周波数成分の信号をノイズとみなし除去してもよい。この例に限らず、当該範囲は任意に定めてよい。
【００２４】
また、第１ノイズ除去部１０３は、第１周波数領域変換部３０２が行う後述する処理にて、折り返し歪み（エイリシアシング）の発生を防ぐために、予め原因となる周波数成分の信号を取り除いてもよい。
【００２５】
第１ノイズ除去部１０３による除去は、バンドパスフィルタなどにて実現することができる。また、折り返し歪みは、ローパスフィルタにて除去することができる。
【００２６】
なお、第１ノイズ除去部１０３は、省略可能である。
【００２７】
体動信号取得部２０１は、生体の動きなどから生ずる測定装置自体の動きに関する信号を取得する。以降、この信号を体動信号と称する。生体の動きは、センサと生体との接触状態および生体内の血流状態を変化させ、測定波形に大きな乱れ（体動雑音）を生じさせる。例えば、測定装置を付けた腕を振ると、生体と測定装置のずれおよび接触圧力の変化等から生ずるノイズ、測定部位の血流変化から生ずるノイズなどが、脈波信号に含まれる。体動信号は、これらのノイズを脈波信号から除去するために用いられる。
【００２８】
図２は、脈波信号取得部１０１で測定された脈波信号の一例を示す図である。図２（Ａ）は、安静時の脈波信号である。図２（Ｂ）は運動時の脈波信号である。脈波信号には脈波と体動雑音、またその他のノイズが含まれる。
【００２９】
図２（Ａ）の脈波信号には、体動雑音はほとんど含まれておらず、脈波信号は脈波とほぼ一致する。一方、図２（Ｂ）の脈波信号は、体動雑音の影響により、乱れた波形となる。体動雑音は、脈波よりも大きいため、体動雑音が含まれる脈波信号では、脈波の波長等を測定することはできない。そのため、体動雑音を除去する必要がある。
【００３０】
なお、以下の説明では、体動信号は、例えば、加速度センサや角速度センサ（ジャイロセンサ）など、測定部位や測定装置自体の動きを検知するセンサにより、測定されることを想定する。
【００３１】
第２直流除去部２０２は、体動信号取得部２０１から取得した体動信号に対し、直流成分の信号を除去する。役割および動作は、第１直流除去部１０２と同様である。
【００３２】
第２ノイズ除去部２０３は、体動信号から予め定められた一定範囲外の信号を除去する。役割および動作は、第１ノイズ除去部１０３と同様である。
【００３３】
第１周波数変換部３０１は、脈波信号からサンプル信号を時間方向で間引くダウンサンプリングを行う。ピーク間隔を高精度に算出するには、サンプリング数は多い方が望ましい。しかし、第１周波数領域変換部３０２および第２周波数領域変換部３１２における処理負荷が高くなり、測定装置の大型化、消費電力量の増加につながる。そこで、負荷が高い周波数領域での処理を抑え、第１周波数領域変換部３０２および第２周波数領域変換部３１２の分解能を小さくするために、ダウンサンプリングを行い、精度は別の機能で担保する。ダウンサンプリングにより、第１周波数領域変換部３０２および第２周波数領域変換部３１２は、後述するＦＦＴの処理量を抑えることができる。
【００３４】
第１周波数領域変換部３０２および第２周波数領域変換部３１２は、入力されたサンプリングデータに基づき、ＦＦＴ（高速フーリエ変換）を行い、スペクトルを算出する。
【００３５】
図３は、脈波信号と体動信号並びに脈波信号のスペクトルと体動信号のスペクトルの一例を示す図である。図３（Ａ）の上が脈波信号、下が体動信号である。図３（Ｂ）は、図３（Ａ）の点線枠で囲まれた範囲におけるスペクトルである。図３（Ｂ）の上が脈波信号のスペクトル、下が体動信号のスペクトルである。測定された脈波信号には、体動雑音が含まれている。このデータにおける本来の脈波の周波数は２．０Ｈｚ近傍であるのだが、体動雑音により、２．３Ｈｚ近傍の周波数が最大となっている。
【００３６】
周波数スペクトル演算部３２１は、第１周波数領域変換部３０２で得られた脈波信号スペクトルを加工し、より正しい脈波信号スペクトルに近づける。加工方法は、例えば、第１周波数領域変換部３０２で得られた脈波信号スペクトルから、第２周波数領域変換部３１２で得られた体動信号スペクトルを差し引くスペクトルサブトラクションがある。
【００３７】
図４は、スペクトルサブトラクションの一例を示す図である。図４（Ａ）は脈波信号のスペクトル、図４（Ｂ）は体動信号のスペクトルを示す。図４（Ａ）および４（Ｂ）のグラフの極大値にはバツ印が付けられており、グラフの最大値には円で囲まれたバツ印が付けられている。
【００３８】
図４（Ａ）および４（Ｂ）のグラフをスペクトルサブトラクションすると、図４（Ｃ）のグラフが得られる。図４（Ａ）では、周波数が２．３Ｈｚ近傍にある極大値が最大値であったが、図４（Ｃ）の当該極大値は、体動信号の出力値の約１０ｄＢ程度、減算されていることが分かる。これにより、図４（Ｃ）の最大値は、２．０Ｈｚ近傍における極大値となる。
【００３９】
周波数ピーク検出部３２２は、スペクトルのピーク（極値）を算出する。ピークの算出は、例えば、スペクトルを周波数で微分し、微分値が０の時点を極値と判定することで求めることができる。
【００４０】
周波数推定部３２３は、脈拍信号のスペクトルから、周波数を推定する。ダウンサンプリングにより、周波数ピーク検出部が算出した極値は、複数あることが想定される。周波数推定部３２３は、これらの極値から脈波の周波数を推定する。推定方法は、ピークの最大値を検出結果としてもよい。また、予め過去のデータ等に基づく評価関数を生成し、評価関数からの予測値に最も近い値を、脈波の周波数としてもよい。
【００４１】
脈拍成分抽出部４０１は、周波数推定部３２３が算出した推定周波数を取得し、周波数領域分析部３００の処理がされていない脈波に対する適応フィルタとして機能する。
【００４２】
適応フィルタは、フィルタの特性を決定するフィルタ係数が可変なために、通過帯域を変更することができるフィルタである。一般には、入力信号に対する出力信号をフィードバックさせて、誤差をより小さくするようにフィルタ係数を決定する。ここでは、出力信号をフィードバックさせるのではなく、脈波信号の推定周波数を、適応フィルタの参照値として用いる。フィルタ係数を求める適応アルゴリズムは、任意に定めてよい。
【００４３】
推定周波数は、第１周波数変換部３０１および第２周波数変換部３１１のダウンサンプリングにより精度が高くない。そのため、適応フィルタを通過することができる周波数は余裕を持たせた方がよい。
【００４４】
図５は、脈拍成分抽出部４０１による処理前と処理後の脈波信号の一例を示す図である。図５（Ａ）は処理前の脈波信号、図５（Ｂ）は処理後の脈波信号である。図５（Ａ）では、第１直流除去部１０２と第１ノイズ除去部１０３によるノイズ除去が行われているものの、時間領域において脈波信号のピークを判別することはできない。図５（Ｂ）では、振幅の変動はあるものの、時間領域での脈波のピークの判別が可能な状態にまでなっていることが分かる。
【００４５】
時間ピーク検出部４０２は、フィルタリング処理した脈波信号からピークを検出する。ピークは、予め定められた一定長の時間間隔（検出時間窓幅）における脈波の最大値とする。検出時間窓の幅は、予め定めておいてもよいし、周波数領域分析部３００が算出した推定周波数を基づき決定してもよい。また、時間ピーク検出部４０２が、前回の検出結果に基づき、逐次更新してもよい。最大値は、例えば、脈波を時間で微分し、微分値が０の時点を極値とし、極値の最大値とするといった方法で求めることができる。
【００４６】
生体情報取得部４０３は、時間ピーク検出部４０２が算出したピークに基づき、脈波信号から脈拍数などの生体情報を算出する。生体情報取得部４０３は、脈拍数を、１拍ごとに逐次算出してもよい。具体的には、時間ピーク検出部４０２が算出したピークと、一つ前のピークとの時間間隔を求める。
【００４７】
脈拍数は、拍動ごとに値が異なる。このことをゆらぎと称する。自律神経解析では、脈拍、心拍のゆらぎを解析する。自律神経を失調すると、このゆらぎが無くなることが知られている。そのため、脈拍数を１拍ごとに測定し、ゆらぎを求める必要がある。
【００４８】
周波数領域分析部３００で推定した周波数は、所定時間における平均周波数であるため、ゆらぎを把握することができないが、時間領域分析部４００は、１拍ごとに脈拍数を測定できるため、ゆらぎを把握することができる。
【００４９】
なお、ここでは取得する生体情報を脈拍数と想定しているが、フィルタリング処理した脈波信号から算出できる生体情報であればよい。
【００５０】
出力部５０１は、算出した脈拍数などの生体情報を出力する。出力方法は、測定装置の図示しない画面に表示してもよい。また、無線通信機能を持たせ、別の通信装置などに、データ送信などによる出力を行ってもよい。
【００５１】
図６は、第１の実施形態に係る測定装置の処理結果の一例を示す図である。図６（Ａ）は、胸につけた当該測定装置とは異なる測定装置による生体の心拍数を示す。図６（Ｂ）は、腕につけた当該測定装置の生体の脈拍数を示す。図６（Ａ）では、心拍数は平均値として平滑化され、曲線で表されている。図６（Ｂ）では、１拍ごとに判定した脈拍数がプロットされている。
【００５２】
図６では、測定時間で測定者の行動を変化させている。開始から３００秒までは、測定者は座った状態（Ｓｉｔｔｉｎｇ）であり、以降３００秒ごとに、歩いている状態（Ｗａｌｋｉｎｇ）、早足の状態（Ｊｏｇｇｉｎｇ）、走っている状態（Ｒｕｎｎｉｎｇ）、座った状態（Ｓｉｔｔｉｎｇ）と変化させている。両図を比較すると、体動の激しいＲｕｎｎｉｇなどの状態でも、ほぼ等しい値を取ることが分かる。これにより、体動雑音の影響を受ける腕型の測定装置にもかかわらず、本装置が高精度に測定できることが分かる。
【００５３】
次に、第１の実施形態の処理フローについて具体的に説明する。
【００５４】
図７は、第１の実施形態に係る測定装置の全体処理の概略フローチャートである。当該処理は、測定装置の電源起動時、またはユーザから操作開始の指示などのタイミングで開始されるものとする。
【００５５】
処理が開始されると、脈波信号処理部１００と体動信号処理部２００は別々に処理を開始する（Ｓ１０１ＡとＳ１０１Ｂ）。脈波信号処理部１００からは脈波信号が、体動信号処理部２００からは体動信号が出力される。
【００５６】
周波数領域分析部３００は、脈波信号処理部１００からの脈波信号と、体動信号処理部２００からの体動信号の両方を取得した後に、これらの信号に基づき、脈波信号の推定周波数を算出する（Ｓ１０２）。時間領域分析部４００は、周波数領域分析部３００からの推定周波数に基づき、脈拍の周期を拍動ごとに求め、１拍あたりの脈拍数を算出する（Ｓ１０３）。出力部５０１は、時間領域分析部４００が算出した脈拍数を出力する（Ｓ１０４）。以上が、概略処理のフローとなる。次に各部の処理の詳細を説明する。
【００５７】
図８は、脈波信号処理部１００による処理のフローチャートである。脈波信号取得部１０１は、脈波信号を取得する（Ｓ２０１）。取得された脈波信号は、第１直流除去部１０２であるバンドパスフィルタ等を介して、直流成分が除去される（Ｓ２０２）。さらに、抽出された脈波信号は、第１ノイズ除去部１０３であるバンドパスフィルタ等を介して、通過可能範囲外の周波数成分が除去される（Ｓ２０３）。以上が、脈波信号処理部１００の処理のフローとなる。
【００５８】
図９は、体動信号処理部２００による処理のフローチャートである。体動信号取得部２０１は、体動信号を取得する（Ｓ３０１）。取得された体動信号は、脈波信号と同様に、第２直流除去部２０２を介して、直流成分が除去され（Ｓ３０２）、第２ノイズ除去部を介して、通過可能範囲外の周波数成分が除去される（Ｓ３０３）。以上が、体動信号処理部２００の処理のフローとなる。
【００５９】
図１０は、周波数領域処理部３００による処理のフローチャートである。周波数領域処理部３００は、脈波信号と体動信号の処理を別々に開始する。
【００６０】
第１周波数変換部３０１は、取得した脈波信号に対しダウンサンプリングを行う（Ｓ４０１Ａ）。サンプリングする周期は、予め定められているものとする。サンプリングされた脈波信号は、第１周波数領域変換部３０２に送られる。
【００６１】
第１周波数領域変換部３０２は、ＦＦＴを行い、サンプリングされた脈波信号を周波数スペクトルに変換する（Ｓ４０２Ａ）。サンプリングされた脈波信号は、周波数スペクトル演算部３１１に送られる。
【００６２】
第２周波数変換部３１１は、第１周波数変換部３０１と同様、取得した体動信号に対しダウンサンプリングを行う（Ｓ４０１Ｂ）。サンプリングする周期は、予め定められているものとする。サンプリングされた脈波信号は、第２周波数領域変換部３１２に送られる。
【００６３】
第２周波数領域変換部３１２は、第１周波数領域変換部３０２と同様、ＦＦＴを行い、サンプリングされた脈波信号を周波数スペクトルに変換する（Ｓ４０２Ｂ）。サンプリングされた脈波信号は、周波数スペクトル演算部３１１に送られる。
【００６４】
周波数スペクトル演算部３２１は、スペクトルサブトラクションにより、脈波信号の周波数スペクトルから、体動信号の周波数スペクトルを減算する（Ｓ４０３）。減算された脈波スペクトルは、周波数ピーク検出部３２２に送られる。
【００６５】
周波数ピーク検出部３２２は、予め定められた手法により、ピークを算出する（Ｓ４０４）。算出したピークは複数あってもよい。算出結果は、周波数推定部３２３に送られる。
【００６６】
周波数推定部３２３は、予め定められた手法により、脈波の周波数を推定する（Ｓ４０５）。以上が、周波数領域処理部３００の処理のフローとなる。
【００６７】
図１１は、時間領域分析部４００による処理のフローチャートである。
【００６８】
脈拍成分抽出部４０１は、脈波信号処理部１００からの脈波信号に対し、周波数領域分析部３００からの推定周波数に基づき、フィルタリングを行う（Ｓ５０１）。適応フィルタのアルゴリズムは予め定められているものとする。フィルタリングされた脈波信号は、時間ピーク検出部４０２に送られる。
【００６９】
時間ピーク検出部４０２は、フィルタリング処理した脈波信号からピークを検出する（Ｓ５０２）。ピーク検出方法は予め定められているものとする。算出されたピークに関する情報は、生体情報取得部４０３に送られる。
【００７０】
生体情報取得部４０３は、ピークに関する情報または脈波信号に基づき、生体情報を算出する（Ｓ５０３）。算出する生体情報は、予め定められているものとする。算出する生体情報は、複数でもよい。以上が、時間領域分析部４００の処理のフローとなる。
【００７１】
以上のように、第１の実施形態によれば、分解能が小さい周波数領域処理を補助的に用いた時間領域処理を行うことにより、処理負荷の軽減と、高精度な脈波の算出を両立させることができる。これにより、測定装置の軽量化および省電力化を図ることができる。また、拍動ごとに１拍あたりの脈拍数を算出することができるため、測定データを自律神経解析などの他の解析にも活用することができる。
【００７２】
（第２の実施形態）
図１２は、第２の実施形態に係る測定装置の概略構成を示すブロック図である。第２の実施形態に係る測定装置は、第１の実施形態に対し、新たに第３信号処理部６００が追加されたものである。
【００７３】
以下、第２の実施形態について説明する。第１の実施形態と重複する説明は省略する。
【００７４】
第３信号処理部６００は、脈波信号処理部１００が処理する脈波信号（第１脈波信号）とは異なる脈波信号（第２脈波信号）に関する処理を行う。
【００７５】
第１脈波信号を測定した測定装置と、第２脈波信号を測定した装置は異なっていてもよいし、同じでもよい。測定装置が同じ場合は、測定方法などが異なることにより、測定信号が異なればよい。例えば、脈波信号を測定する反射型光電脈波センサが、赤色の光を受信して第１脈波信号を測定した場合に、第２脈波信号は、赤色以外の光を受信して測定された信号であればよい。
【００７６】
第３信号処理部６００は、第３信号取得部６０１と、第３直流除去部６０２と、第３ノイズ除去部６０３を備える。第３信号処理部６００の内部の各部の機能は、脈波信号処理部１００内部の対応する各部の機能と同一であるため、省略する。第３信号取得部６０１は、脈波信号処理部１００の脈波信号取得部１０１と、第３直流除去部６０２は、第１直流除去部１０２と、第３ノイズ除去部６０３は第１ノイズ除去部１０３と対応する。
【００７７】
第３信号処理部６００により加工された第２脈波信号は、第１脈波信号同様に、周波数領域分析部３００と時間領域分析部４００に送られ、第１の実施形態同様に、処理される。これにより、第１脈波信号に基づく生体情報と、第２脈波信号に基づく生体情報の２つを得ることができる。得られた２つの生体情報を比較することで、新たな測定情報を提供することでもできる。
【００７８】
生体情報取得部４０３は、第１脈波信号に基づく生体情報と、第２脈波信号に基づく生体情報を個別に算出してもよい。または、２つの生体情報の比較に基づく情報を算出してもよい。
【００７９】
なお、第２脈波信号は、周波数領域分析部３００に送らずに、第１脈波信号の推定情報に基づく脈拍成分抽出部４０１のフィルタリングを受けてもよい。第１脈波信号と第２脈波信号の周波数が大幅に異なる場合以外は、問題が生じる可能性は低いからである。
【００８０】
第２の実施形態における処理のフローチャートは、脈波信号処理部１００および当該内部の各部の処理が、それぞれ第３信号処理部６００および当該内部の各部に置き換えたものであるため、省略する。
【００８１】
なお、上記の第２の実施形態は、新たに第３信号処理部６００を備えた形態であるが、さらに第４から第ｎ（ｎは５以上の整数）の信号処理部を備えてもよい。その場合も、第４から第ｎの信号処理部の各部の機能および動作は、第３信号処理部６００と同じである。
【００８２】
以上のように、第２の実施形態によれば、センサまたは測定方法等の異なる２種類の生体情報を取得することができ、取得した生体情報を比較した結果に基づき、新たな測定情報を提供することができる。また、第３信号処理部６００と同様の構成の第４から第ｎの信号処理部を用いることで、複数種の生体情報を取得し、取得した複数種の生体情報を比較することで新たな測定情報を提供することができる。さらに、複数の信号処理部の出力の平均をとることで、測定精度を高めることも可能である。
【００８３】
また、上記に説明した実施形態における各処理は、ソフトウェア（プログラム）によって実現することが可能である。よって、上記に説明した実施形態における測定装置は、例えば、汎用のコンピュータ装置を基本ハードウェアとして用い、コンピュータ装置に搭載されたプロセッサにプログラムを実行させることにより実現することが可能である。
【００８４】
図１３は、本実施形態に係る測定装置を実現したハードウェア構成例を示すブロック図である。測定装置は、プロセッサ７０１、主記憶装置７０２、補助記憶装置７０３、デバイスインタフェース７０４、ネットワークインタフェース７０５、入力装置７０６、出力装置７０７を備え、これらがバス７０８を介して接続されたコンピュータ装置として実現できる。
【００８５】
プロセッサ７０１が、補助記憶装置７０３からプログラムを読み出して、主記憶装置７０２に展開して、実行することで、脈波信号処理部１００と、体動信号処理部２００と、周波数領域分析部３００と、時間領域分析部４００の機能を実現することができる。
【００８６】
本実施形態の測定装置は、当該測定装置で実行されるプログラムをコンピュータ装置に予めインストールすることで実現してもよいし、プログラムをＣＤ−ＲＯＭなどの記憶媒体に記憶して、あるいはネットワークを介して配布して、コンピュータ装置に適宜インストールすることで実現してもよい。
【００８７】
主記憶装置７０２は、プロセッサ７０１が実行する命令、および各種データ等を一時的に記憶するメモリ装置であり、ＤＲＡＭ等の揮発性メモリでも、ＭＲＡＭ等の不揮発性メモリでもよい。補助記憶装置７０３は、プログラムやデータ等を永続的に記憶する記憶装置であり、例えば、フラッシュメモリ等がある。
【００８８】
デバイスインタフェース７０４は、出力結果などを記録する外部記憶媒体と接続するＵＳＢなどのインタフェースである。外部記憶媒体は、ＨＤＤ、ＣＤ−Ｒ、ＣＤ−ＲＷ、ＤＶＤ−ＲＡＭ、ＤＶＤ−Ｒ、ＳＡＮ（Ｓｔｏｒａｇｅ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）等の任意の記録媒体でよい。また、デバイスインタフェース７０４を介して、図示しない反射型光電脈波センサなどの脈波を測定するセンサ、加速度センサや角速度センサ（ジャイロセンサ）など加速度を測定するセンサなどと接続されていてもよい。
【００８９】
ネットワークインタフェース５０４は、無線ＬＡＮなどのネットワークに接続するためのインタフェースである。出力部５０１は、ネットワークインタフェース５０４を介して、出力結果などを他の通信装置に送信してもよい。
【００９０】
ユーザは、入力装置７０６から測定開始などの情報を入力してもよい。また出力装置７０８は、画像を表示する表示装置でもよい。
【００９１】
上記に、本発明の一実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００９２】
１００  脈波信号処理部
１０１  脈波信号取得部
１０２  第１直流除去部
１０３  第１ノイズ除去部
２００  体動信号処理部
２０１  体動信号取得部
２０２  第２直流除去部
２０３  第２ノイズ除去部
３００  周波数領域分析部
３０１  第１周波数変換部
３０２  第１周波数領域変換部
３１１  第２周波数変換部
３１２  第２周波数領域変換部
３２１  周波数スペクトル演算部
３２２  周波数ピーク検出部
３２３  周波数推定部
４００  時間領域分析部
４０１  脈拍成分抽出部
４０２  時間ピーク検出部
４０３  生体情報取得部
５０１  出力部
６００  第３信号処理部
６０１  第３信号取得部
６０２  第３直流除去部
６０３  第３ノイズ除去部
７０１  プロセッサ
７０２  主記憶装置
７０３  補助記憶装置
７０４  デバイスインタフェース
７０５  ネットワークインタフェース
７０６  入力装置
７０７  出力装置
７０８  バス
【図１】
0006464004.tif 000002
【図２】
0006464004.tif 000003
【図３】
0006464004.tif 000004
【図４】
0006464004.tif 000005
【図５】
0006464004.tif 000006
【図６】
0006464004.tif 000007
【図７】
0006464004.tif 000008
【図８】
0006464004.tif 000009
【図９】
0006464004.tif 000010
【図１０】
0006464004.tif 000011
【図１１】
0006464004.tif 000012
【図１２】
0006464004.tif 000013
【図１３】
0006464004.tif 000014