拍動検出装置
【課題】 拍動検出装置の検出性能を高めること。
【解決手段】 拍動検出装置１００は、脈波センサー１０と、フィルター部３０と、周波数解析部５０と、を含み、周波数解析部５０は、フィルター後信号の周波数スペクトルに基づいて、拍動成分と体動信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行して拍動呈示スペクトルを特定する第１処理部６０と、脈波信号の周波数スペクトルに基づいて、拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部７０と、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部８０と、を含む。
【解決手段】 拍動検出装置１００は、脈波センサー１０と、フィルター部３０と、周波数解析部５０と、を含み、周波数解析部５０は、フィルター後信号の周波数スペクトルに基づいて、拍動成分と体動信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行して拍動呈示スペクトルを特定する第１処理部６０と、脈波信号の周波数スペクトルに基づいて、拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部７０と、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部８０と、を含む。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、拍動検出装置等に関する。
【背景技術】
【０００２】
拍動検出装置は、人体の心拍に由来する拍動を検出するための装置であって、例えば、腕、手のひら、手指などに装着される脈波センサーからの信号（脈波信号）から、人体の体動の影響により発生する信号成分（体動影響信号）を雑音として除去し、心拍に由来する信号（拍動信号）を検出する装置である。
【０００３】
人の指や手首に装着するタイプの脈波センサーは、例えば、特許文献１〜特許文献３に記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００５−１９８８２９号公報
【特許文献２】特開２００７−５４４７１号公報
【特許文献３】特開２００５−１３１４２６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
拍動信号成分を明確に峻別するに際し、脈波センサーから出力される脈波信号に含まれる拍動信号成分の信号レベルが著しく低い場合には、信号レベルが十分な場合に比べる、峻別処理の難易度が高くなる。
【０００６】
本発明の少なくとも一つの態様によれば、例えば、拍動検出装置の検出性能を高めることができる。
【課題を解決するための手段】
【０００７】
（１）本発明の拍動検出装置の一態様は、被検体の拍動に由来する拍動信号を検出する拍動検出装置であって、前記拍動信号と、前記被検体の体動に由来する体動信号を含むノイズ信号とが混在した脈波信号を検出して出力する脈波センサーと、前記脈波信号から、前記ノイズ信号の少なくとも一部を除去するフィルター部と、前記フィルター部から出力されるフィルター後信号、または前記脈波センサーから出力される前記脈波信号に基づいて、所定時間毎に前記周波数解析処理を行って、前記拍動信号を示す拍動呈示スペクトルを特定する周波数解析部と、を含み、前記周波数解析部は、前記フィルター後信号の周波数スペクトルに基づいて、前記拍動信号と前記ノイズ信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行して前記拍動呈示スペクトルを特定する第１処理部と、前記脈波信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部と、前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部と、を含む。
【０００８】
本態様では、周波数解析部は、第１処理部、第２処理部ならびに判断部を有する。第１処理部は、脈波センサーによって検出された脈波信号をフィルター部によってフィルタリングして得られるフィルター後信号に基づいて、拍動信号を示す拍動呈示スペクトルを特定する。第１処理部は、拍動信号と体動信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行する。
【０００９】
本態様では、例えば、脈波センサーから出力される脈波信号に含まれる拍動信号成分の信号レベルが著しく低い場合においても、第２処理部によって拍動呈示スペクトルの捕捉を試みることができ、これによって拍動検出の可能性が向上する。この第２処理部は、フィルター前の脈波信号に基づいて拍動呈示スペクトルの捕捉を試みる。
【００１０】
フィルタリング処理を行った場合、残留した体動信号による体動スペクトルと、拍動呈示スペクトルとの大きさが同程度になることがあり、この場合、単純にスペクトルの大きさを比較するだけでは、拍動呈示スペクトルがどれなのかを判定することが難しい。この場合であっても、フィルター前の脈波信号に基づいた拍動呈示スペクトルの捕捉処理を試みることによって、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低い拍動信号成分を検出することができる場合がある。
【００１１】
但し、フィルター前の脈波信号は、ノイズ対策がまったく行われていない原信号であり、例えば、過大なノイズしか存在しない信号である場合もあり得る。そこで、本態様では、第２処理部による拍動呈示スペクトルの捕捉処理を試みる前に、判断部が、拍動呈示スペクトル捕捉処理の実行の可否を判断することとした。したがって、本態様によれば、拍動検出装置の検出性能を高めることができる。
【００１２】
（２）本発明の拍動検出装置の他の態様では、前記判断部は、前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、を有し、前記評価部による評価に基づいて、前記第２処理部による前記拍動呈示スペクトル捕捉処理の実行の可否を判断する。
【００１３】
本態様では、判断部は、脈波信号のきれいさに関する判断指標を取得する判断指標取得部と、判断指標を用いて、脈波信号に含まれるノイズ量の程度を評価する評価部とを有する。判断部は、脈波信号のノイズ状態を判断するための判断指標を取得し、その判断指標を用いて、脈波信号に含まれるノイズ量の程度を評価する。そして、その評価に基づいて、第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する。
【００１４】
判断部は、例えば、評価部によって、脈波信号に混入されるノイズ量（例えば外乱ノイズ成分のノイズ量）が所定の基準以下であると判断され場合に、第２処理部による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。換言すれば、脈波信号が、許容レベル以上にきれいであることを第２処理部による拍動呈示スペクトル捕捉処理の実行の条件としてもよい。
【００１５】
第２処理部による救済を試みると、ノイズを拍動信号成分と誤検出する可能性が高い。本態様によれば、例えば、外乱ノイズが多い状態では拍動呈示スペクトルの捕捉処理が実行されない。本態様によれば、脈波信号のノイズ状態を慎重に判断しつつ、第２処理部による拍動呈示スペクトルの捕捉処理の可否を判断することができる。
【００１６】
（３）本発明の拍動検出装置の他の態様では、前記判断部は、前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、を有し、前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、前記評価部による評価が許容範囲にあることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断する。
【００１７】
本態様では、周波数解析部による前回の周波数解析によって、拍動呈示スペクトルが検出されていること、ならびに、評価部による、脈波信号に含まれるノイズ量の程度の評価が許容範囲にあることを条件として、第２処理部による拍動呈示スペクトルの捕捉処理を実行する。
【００１８】
拍動呈示スペクトルは、所定時間前の過去の拍動信号成分と現在の拍動信号成分との相関性に基づいて捕捉することができる。よって、拍動呈示スペクトルの捕捉処理を実行する条件として、前回の周波数解析処理、すなわち、所定時間前の周波数解析処理において拍動信号の検出に成功していることを条件としてもよい。本態様では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【００１９】
（４）本発明の拍動検出装置の他の態様では、前記判断部は、前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、前記被検体の運動状態を判断する運動状態判断部と、を有し、前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、前記評価部による評価が許容範囲にあり、かつ、前記運動状態判断部によって前記被検体が定常状態にあると判断されることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断する。
【００２０】
被検体の運動状態が急激に変化している場合には、脈拍等の変化が急である。
【００２１】
よって、本態様では、よって、拍動呈示スペクトルの捕捉処理を実行する条件として、運動状態判断部によって被検体が定常状態にあると判断されることを条件としてもよい。被検体が定常状態である場合とは、例えば、被検体が安定した運動状態や静養状態にある場合をいう。
【００２２】
被検体の運動状態は、例えば、体動センサーから得られる体動信号のスペクトルが、所定の時間にわたって大きく変化しないという観点から判断することができる（但し、これに限定されるものではない）。本態様では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【００２３】
（５）本発明の拍動検出装置の他の態様では、前記判断部は、前記第１処理部が拍動呈示スペクトルを特定することができなかった場合に、前記判断部は、前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する。
【００２４】
本態様では、判断部による拍動呈示スペクトル捕捉処理の実行の可否判断は、第１処理部が、フィルター後信号に基づいて拍動呈示スペクトルを特定することができなかった場合に行われる。
【００２５】
第１処理部と第２処理部とを並列に動作させることも可能であるが、第１処理部によって拍動呈示スペクトルが特定できた場合には、第２処理部による拍動呈示スペクトルの捕捉は不要である。そこで、第２処理部の無駄な動作をなくすために、本態様では、第１処理部による拍動呈示スペクトルの特定の失敗を条件として、判断部による拍動呈示スペクトルの捕捉処理の可否判断を開始させることとした。
【００２６】
（６）本発明の拍動検出装置の他の態様では、前記被検体の体動に伴う体動信号を検出して出力する体動センサーと、前記体動センサーから出力される体動信号に基づいて周波数分解処理を行う周波数分解部と、を有し、前記周波数分解部による周波数分解結果に基づいて推定される、前記被検体の体動由来の体動成分を体動基線とし、異なる周波数のｍ本（ｍは１以上の整数）の前記体動基線のうちの第ｎ番目（１≦ｎ≦ｍ）の体動基線を第ｎ体動基線としたとき、前記第２処理部は、周波数解析処理に使用される前記脈波信号の周波数スペクトルにおける前記被検体の拍動由来の拍動成分と推定される第１拍動基線と、前記所定時間前の周波数解析処理に使用された前記脈波信号または前記フィルター後信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルとして特定された第２拍動基線との相関性を重視した第１捕捉処理を実行して、前記拍動呈示スペクトルの捕捉を試みる第１捕捉処理部と、前記脈波信号に、前記第１処理部による前記拍動／ノイズ分離処理と同じ処理を施した後、前記第１拍動基線と、前記第ｎ体動基線との区別の可否に基づいて第２捕捉処理を実行して前記拍動呈示スペクトルの捕捉を試みる第２捕捉処理部と、を有し、前記判断部は、前記拍動呈示スペクトル捕捉処理を実行可と判断した場合、さらに、前記評価部による評価が第１の許容レベルであるか、前記第１許容レベルよりも低い第２許容レベルであるかを判断し、前記判断部によって前記第１許容レベルであると判断されたときは、前記第１捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行し、前記第２許容レベルであると判断されたときは、前記第２捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行する。
【００２７】
本態様では、第２処理部の構成と動作の一例が示される。第２処理部は、第１捕捉処理部および第２捕捉処理部を有する。判断部は、拍動呈示スペクトル捕捉処理を実行可と判断した場合、さらに、評価部による評価（ノイズ量の程度に基づく評価）が第１の許容レベルであるか、第１許容レベルよりも低い第２許容レベルであるかを判断する。
【００２８】
第１捕捉処理部は、判断結果が第１許容レベルのとき（例えば、ノイズ量の程度が小さく、脈波信号がきれいであるとき）に動作する。第２捕捉処理部は、判断結果が、第１許容レベルよりも低い第２許容レベルのとき（例えば、ノイズ量の程度が第１許容レベルを基準とする場合よりも大きく、脈波信号のきれいさが中程度のとき）に動作する。
【００２９】
そこで、第１捕捉処理部は、第１捕捉処理を実行する。この第１捕捉処理では、第１拍動基線と第２拍動基線との相関性に基づく拍動呈示スペクトルの捕捉処理を実行する。第１拍動基線は、今回の周波数解析処理に使用される、「脈波信号」の周波数スペクトルにおける被検体の拍動由来の拍動成分と推定されるスペクトルである。第２拍動基線とは、前回の周波数解析処理に使用された、「脈波信号またはフィルター後信号」の周波数スペクトルに基づいて拍動呈示スペクトルとして特定された拍動基線である。
【００３０】
第１捕捉処理部は、第１拍動基線と第２拍動基線との相関を重視して判断したとき、相関性が有る場合には、第１拍動基線を拍動呈示スペクトルと判断する。
【００３１】
また、第２捕捉処理部は、例えば、脈波信号が中程度である場合に動作することができる。この場合には、拍動呈示スペクトルの捕捉のためには、拍動基線と外乱ノイズとの分離が必要となる。
【００３２】
第２捕捉処理部は、第１処理部で行われる拍動／ノイズ分離処理と同じ処理を実行してもよい。また、第２捕捉処理部は、さらに、今回の検出時の周波数スペクトルにおける第１拍動基線と、第ｎ体動基線との区別の可否に基づいて拍動呈示スペクトルの捕捉を試みてもよい。
【００３３】
ここで、体動基線（広義）は、体動センサーから出力される体動信号の、現在ならびに過去の周波数分解結果に基づいて、体動成分と推定されるスペクトルの総称である。外乱ノイズがある程度存在する場合、例えば、被検体が運動している場合等では、体動基線には、異なる周波数のｍ本（ｍは１以上の整数）の体動基線が含まれることが多い。
【００３４】
第ｎ体動基線は、異なる周波数のｍ本（ｎ、ｍは１以上の整数）の体動基線のうち、一番周波数の低い基線を第１番目としそのｎ倍の周波数に相当する基線を第ｎ体動基線と定義する。例えば、第１波（基本波）の拍動基線が、第ｎ体動基線（１≦ｎ≦ｍ）の位置とは異なる位置にある場合には、第２捕捉処理部は、第１拍動基線が拍動呈示スペクトルであると判断する。この場合、例えば、ｎが２〜４である場合の体動基線は、ｎが１である場合の体動基線よりも信号値が小さく、外乱ノイズ中に潜んでいる場合があるため、第ｎ体動基線（例えばｎ＝１〜４）との区別を試みることは、拍動呈示スペクトルの捕捉精度を高める点で有効である。
【００３５】
なお、第２捕捉処理においては、第１拍動基線と第２拍動基線との相関性を考慮してもよい。例えば、今回検出時（現在の検出期間における検出時）の周波数スペクトルにおける第１拍動基線と、前回検出時（所定時間前の検出持）の周波数スペクトルにおける第２拍動基線との間の周波数軸上における距離の観点からの判断を行ってもよい。
【００３６】
本態様では、脈波信号に含まれるノイズ量の程度に応じて、拍動呈示スペクトルの捕捉処理の内容を切り替えることから、脈波信号の状態に応じた適切な拍動呈示スペクトルの捕捉処理を実行することができる。
【００３７】
（７）本発明の拍動検出装置の他の態様では、前記第１捕捉処理部は、前記第１拍動基線の周波数軸上における第１位置と、前記第２拍動基線の周波数軸上における第２位置との間の距離によって、前記相関性を判断する。
【００３８】
本態様では、第１捕捉処理部の動作の一例が示される。今回の検出（現在の検出期間における検出）に係る第１拍動基線の周波数軸上の第１位置と、前回の検出（所定時間前の検出）に係る第２拍動基線の周波数軸上における第２位置との間の距離を、相関性の判断の基準としてもよい。例えば、周波数軸上における距離が、例えば、距離２以内であるとき、つまり、周波数軸上において、第２拍動基線の第２位置を基準として、第１拍動基線が２本のスペクトル以内で隣り合っているとき、第１捕捉処理部は、相関ありと判断して、第１拍動基線を拍動呈示スペクトルと判断してもよい。
【００３９】
（８）本発明の拍動検出装置の他の態様では、前記第２捕捉処理部は、前記第１拍動基線が前記第ｎ体動基線と重なっていないという第１条件と、前記第１拍動基線が前記第２拍動基線に対して、周波数軸上で所定の距離以内であるという第２条件と、を満足する場合に、前記第１拍動基線を、拍動呈示スペクトルとして特定する。
【００４０】
本態様では、第２捕捉処理部は、第１条件の他に、第２条件が満足されるか否かの判断も実行し、より慎重に拍動呈示スペクトルを捕捉する。
【００４１】
例えば、第１拍動基線が、体動センサー出力の周波数分解結果に基づいて特定される第ｎ体動基線に重なっていないことを第１条件とし、第１拍動基線が、周波数軸上で、第２拍動基線に対して所定距離内にあることを、第２条件としてもよい。
【００４２】
「所定距離内」は、例えば、距離２以内とすることができる。距離２とは、第１拍動基線が、第２拍動基線に対して、１本のスペクトルを隔てて隣り合っていることである。距離２以内とは、距離０、距離１または距離２のいずれかであることを意味する。すなわち、本態様では、第２捕捉処理部は、周波数軸上でスペクトル（基線）が３本分離れるということが発生しないものとして、拍動呈示スペクトルの呈示処理を行う。上述の第１条件と第２条件を満足する場合には、第２捕捉処理部は、第１拍動基線を拍動呈示スペクトルとして捕捉する。これによって、例えば、目立つノイズに引きずられることなく、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような拍動呈示スペクトルを特定できる可能性が高まる。
【００４３】
このように、本発明の少なくとも一つの態様によれば、例えば、拍動検出装置の検出性能を高めることができる。
【図面の簡単な説明】
【００４４】
【図１】本発明の拍動検出装置の一例の構成を示す図
【図２】図２（Ａ）および図２（Ｂ）は、第１捕捉処理部による第１捕捉処理の一例を示す図
【図３】図３（Ａ）および図３（Ｂ）は、第２捕捉処理部による第２捕捉処理の一例を示す図
【図４】本発明の拍動検出装置の他の例の構成を示す図
【図５】図５（Ａ）〜図５（Ｃ）は、脈波信号に含まれるノイズ量の程度の判断動作を説明するための脈波信号の波形および周波数スペクトルの一例を示す図
【図６】第２処理部による、拍動呈示スペクトルの捕捉処理（第２処理）の手順の一例を示すフローチャート
【図７】図７（Ａ）および図７（Ｂ）は、第２捕捉処理の適用例を示す図
【図８】図８（Ａ）〜図８（Ｃ）は、第１捕捉処理の適用例を示す図
【図９】図９（Ａ）〜図９（Ｄ）は、外乱ノイズの混入が著しい例を示す図
【図１０】腕時計型の脈拍計測装置を被検体の手首に装着し、被検体が歩行（ウォーキング）を行ったときにおける、検出された脈拍数の変化を示す図
【図１１】図１１（Ａ）および図１１（Ｂ）は、拍動検出装置の、被検体への装着例を示す図
【図１２】図１２（Ａ）および図１２（Ｂ）は、拍動信号成分が強い場合におけるフィルタリングの効果を示す図
【図１３】図１３（Ａ）および図１３（Ｂ）は、拍動信号成分が弱い場合におけるフィルタリングの効果を示す図、
【図１４】図１４（Ａ）および図１４（Ｂ）は、被検体が安定運動している最中に、非周期的な外乱ノイズが混入した場合の例を示す図
【発明を実施するための形態】
【００４５】
本実施形態の説明に入る前に、フィルター後信号に基づいて拍動呈示スペクトルを特定するのが困難である場合について考察する。フィルター後信号とは、脈波センサーから出力される脈波信号に、フィルター部によるフィルタリング処理を施すことによって得られる信号である。フィルタリング処理は、脈波信号に含まれるノイズを除去あるいは抑制するための信号処理である。
【００４６】
図１２（Ａ）および図１２（Ｂ）は、拍動信号成分が強い場合におけるフィルタリングの効果を示す図である。図１２（Ａ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図１２（Ｂ）の１６秒分フィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【００４７】
図１２（Ａ）および図１２（Ｂ）の例は、拍動信号成分が強い場合の例である。図１２（Ａ）では、拍動信号のスペクトルＳ２０の信号値は、比較的大きなノイズ成分のスペクトルＮ２１〜Ｎ２３と比べて十分に大きい。したがって、脈波信号とノイズとを区別するためのフィルタリング処理を実行すると、図１２（Ｂ）に示すように、ノイズは十分に低減されてＳ／Ｎが向上し、拍動信号のスペクトルＳ３０をノイズと区別することが容易となる。
【００４８】
ここで、図１３（Ａ）および図１３（Ｂ）を参照する。図１３（Ａ）および図１３（Ｂ）は、拍動信号成分が弱い場合におけるフィルタリングの効果を示す図である。図１３（Ａ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図１３（Ｂ）の上側には１６秒分のフィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【００４９】
図１３（Ａ）では、ノイズのスペクトルＮ４０，Ｎ４１の信号値が大きく、これに対して、拍動呈示スペクトルＳ４０の信号値は小さく、拍動呈示スペクトルＳ４０はノイズに拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低い状態となっている。図１３（Ｂ）の例では、フィルタリングによるノイズ除去の効果が顕著に表れ、フィルタリング後のＳ／Ｎは、フィルター前のＳ／Ｎよりも格段に向上している。
【００５０】
しかしながら、フィルタリング処理を行った場合、残留した体動信号による体動スペクトルと、拍動呈示スペクトルとの大きさが同程度になることがあり、この場合、単純にスペクトルの大きさを比較するだけでは、拍動呈示スペクトルがどれなのかを判定することが難しい。
【００５１】
図１４（Ａ）および図１４（Ｂ）は、被検体が安定運動している最中に、非周期的な外乱ノイズ（例えば脈波センサーを物体に衝突させた場合の衝撃ノイズ）が混入した場合の例を示す図である。図１４（Ａ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図１４（Ｂ）の上側には１６秒分のフィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【００５２】
図１４（Ａ）では、拍動呈示スペクトルＳ６０の周波数よりも低い周波数帯域において、大きな外乱ノイズのスペクトルＮ６０が存在している。図１４（Ｂ）では、フィルタリングの効果によって、外乱ノイズのスペクトルＮ６０は抑圧されるが、拍動呈示スペクトルＳ６０の信号値も低下しており、結果的に、拍動呈示スペクトルＳ６０の特定が困難である。
【００５３】
以下に説明する本実施形態の拍動検出装置では、このような状況であっても、フィルター前の脈波信号に基づく周波数解析を実行して、拍動呈示スペクトルの捕捉を試みる。ここで、図１３（Ａ）および図１３（Ｂ）に示す例に着目する。
【００５４】
図１３（Ｂ）の例では、拍動呈示スペクトルＳ５０が、ノイズのスペクトルＮ５０と同程度の大きさである。また、拍動呈示スペクトルＳ５０と近い周波数帯に外乱ノイズのスペクトルも多数存在しており、よって、図１２の場合と比較すると拍動呈示スペクトル特定処理の難易度が高い。
【００５５】
これに対し、図１３（Ａ）の例では、Ｎ４０、Ｎ４１といった体動呈示スペクトルが存在するものの、これらとは離れた周波数帯に、拍動呈示スペクトルＳ４０が存在する。また、更には拍動呈示スペクトルＳ４０は、隣り合う他のスペクトルと比較すると相対的には大きいスペクトルとなっている。よって、脈波信号ｄが図１３（Ａ）のような信号状態であることを検出した上であれば、Ｓ４０が確かに拍動呈示スペクトルであると容易に判定可能と言える。
【００５６】
つまり、例えば、所定時間前の過去と現在の拍動信号成分の相関性等に着目すれば、拍動呈示スペクトルＳ４０が存在する周波数域から離れた周波数域における大きなノイズのスペクトルに惑わされることなく、拍動呈示スペクトルを特定する可能性があるといえる。
【００５７】
そこで、以下の本実施形態では、拍動検出装置に、フィルター前の脈波信号に基づく周波数解析の機能を付与し、新規な検出アルゴリズムを用いて拍動呈示スペクトルの捕捉を行えるようにした。
【００５８】
但し、フィルター前の脈波信号は、ノイズフィルタリング処理を行っていない原信号であり、例えば、拍動信号成分と比較して大きなノイズが含まれる場合があり得る。この場合に該当するのが、図１４に示す例である。この場合には、拍動呈示スペクトルの捕捉処理を行わないのが好ましい。
【００５９】
このような観点から、以下に述べる実施形態においては、判断部を設けて、拍動呈示スペクトルの捕捉処理を試みる前に、判断部が、拍動呈示スペクトル捕捉処理の実行の可否を判断することとした。
【００６０】
判断部は、例えば、脈波信号に含まれるノイズが多い状態である場合や、被検体（人ならびに動物を含む）の運動状態が急激に変化している場合、所定時間前の過去の検出時点において拍動呈示スペクトルの検出に失敗している場合等においては、拍動呈示スペクトルの捕捉処理は不可と判断してもよい（但し、上述の場合は一例であり、これらの場合に限定されるものではない）。
【００６１】
以下、図１〜図３を参照して、具体的に説明する。なお、以下に説明する本実施形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また、本実施形態で説明される構成の全てが、本発明の必須構成要件であるとは限らない。
【００６２】
（第１実施形態）
図１は、本発明の拍動検出装置の一例の構成を示す図である。図１に示される拍動検出装置１００は、被検体（人や動物を含む）の拍動に由来する拍動信号、拍動信号に対応する心拍等の生体情報等を検出するセンサー装置の一種である。
【００６３】
ここで、拍動とは、医学的には心臓のみならず内臓一般の周期的な収縮、弛緩が繰り返された場合に起こる運動のことをいう。ここでは、心臓が周期的に血液を送るポンプとしての動きを拍動と呼ぶ。なお、心拍数とは、１分間の心臓の拍動の数をいう。また、脈拍数は、末梢血管における脈動の数をいう。心臓が血液を送り出す際に、動脈に脈動が生じるので、この回数を数えたものを脈拍数あるいは単に脈拍と呼ぶ。腕で脈を計測する限りは、医学的には心拍数とは呼ばずに脈拍数と呼ぶのが通常である。また、以下の説明では、体動という用語が使用される。体動とは、広い意味では、体を動かすことすべてを意味する（広義の体動）。体動センサーから出力される体動信号は、広義の体動に対応する信号であり、体動センサー信号と言い換えることができる。但し、体動信号の中には、被検体の周期的な体動（例えば、歩行・ジョギングなどに伴う定常的、周期的な腕（脈拍計の装着部位近辺）の振り等）に伴う周期的な体動信号が含まれる。この周期性をもつ体動信号は、狭義の体動（周期性のある体動）に対応する信号である。なお、例えば、広義の体動信号を「体動センサー信号」という場合、狭義の体動信号を、単に「体動信号」といってもよい。
【００６４】
（拍動検出装置の構成例）
図１に示される拍動検出装置１００は、脈波センサー１０と、体動センサー（加速度センサー等）２０と、フィルター部（前処理部）３０と、周波数解析部５０と、脈拍数算出部９０と、表示処理部９２と、表示部９４と、を有する。
【００６５】
脈波センサー１０は、拍動信号と、被検体の体動に由来する体動ノイズ信号を含むノイズ信号とが混在した脈波信号ｄを検出して出力する。脈波センサー１０は、例えば、光電脈波センサー、およびその原理に基づく脈波センサーとすることができる。
【００６６】
体動センサー２０は、被検体の体動、例えば、歩行やジョギングなどに伴う定常的、周期的な腕（脈拍計の装着部位近辺）の動きを検出するセンサーであり、例えば、加速度センサーやジャイロセンサーを含むことができる。
【００６７】
フィルター部３０は、脈波信号ｄおよび体動信号（体動センサー２０から出力される信号）ｆを入力として、フィルタリング処理を実行して、脈波信号ｄに含まれるノイズを最小化する。フィルター部３０は、特に限定されるものではないが、適応フィルター（図１では不図示、図５の参照符号３２）を含んでもよい。
【００６８】
適応フィルターは、入力信号に基づいて、伝達関数を自己適応的に変化させることができるデジタルフィルターである。適応フィルターは一般にデジタル信号処理を行うデジタルフィルターとして実装される。フィルターの所望の性能（例えば、入力信号に含まれるノイズ成分を最小化するといった性能）を維持するために、所定のアルゴリズム（例えば最適化アルゴリズム）が使用される。所定のアルゴリズムには、フィルタリング後の信号に基づいて得られる信号がフィードバックされ、そのアルゴリズムによる適応処理によってフィルター係数が適応的に変化し、その結果として、適応フィルターの周波数応答特性が変化する。適応フィルターを使用することによって、例えば、周期性のある信号と非周期的なノイズとを分離することができる。
【００６９】
周波数解析部５０は、フィルター部３０から出力されるフィルター後信号ｅ、または脈波センサー１０から出力される脈波信号ｄに基づいて、所定時間毎（例えば４秒毎）に周波数解析処理を行って、拍動信号を示す拍動呈示スペクトルを特定する。ここで、フィルター後信号ｅは、脈波信号ｄをフィルタリングして得られる信号である。
【００７０】
周波数解析部５０は、高速フーリエ変換（ＦＦＴ）によって実時間上の信号を周波数軸上の信号に変換する周波数分解部４０（第１周波数分解部４０ａ、第２周波数分解部４０ｂならびに第３周波数分解部４０ｃを含む）と、フィルター後信号ｅのＦＦＴ結果ｇ１ならびに体動信号ｆのＦＦＴ結果ｇ２に基づいて、拍動呈示スペクトルを特定する第１処理部（後処理部）６０と、脈波信号ｄのＦＦＴ結果ｇ３に基づいて、拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部（捕捉処理部）７０と、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部８０と、を含む。
【００７１】
なお、第１周波数分解部４０ａは、フィルター後信号ｅにＦＦＴ（高速フーリエ変換）を施す。第２周波数分解部４０ｂは、体動信号ｆにＦＦＴを施す。また、第３周波数分解部４０ｃは、脈波信号ｄにＦＦＴを施す。
【００７２】
第１処理部（後処理部）６０は、ピーク順ソート部６２と、相関判断部６４と、拍動／ノイズ分離部６６と、拍動呈示スペクトル特定部６８と、を含むことができる。
【００７３】
第２処理部（捕捉処理部）７０は、相関性に基づく捕捉処理を実行する第１捕捉処理部７２と、第２捕捉処理部７４（図１において破線で囲んで示すブロック）とを含む。第２捕捉処理部７４は、第１処理部６０の構成要素であるピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６に、距離と区別性に基づく捕捉部７６を追加して構成される。ピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６と同等の機能をもつブロックを新たに設けてもよい。但し、この場合、同じ機能をもつブロックが重複する。よって本実施形態では、第１処理部６０の構成要素であるピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６を、第２捕捉処理部７４の構成要素としても活用している。
【００７４】
第１処理部６０による第１処理（後処理）は、脈波信号ｄに含まれる拍動信号成分の信号値（振幅値）が大きい場合に有効である。第２処理部７０による第２処理（拍動呈示スペクトルの捕捉処理）は、例えば、第１処理部６０による拍動呈示スペクトルの特定に失敗した場合に行うことができる（但し、これに限定されるものではない）。
【００７５】
判断部８０は、周波数解析部５０に含まれる第３周波数分解部４０ｃが生成した周波数スペクトルｇ３に基づいて、第２処理部７０による拍動呈示スペクトルの捕捉処理の可否を判断する。この判断は、例えば、脈波信号ｄに含まれるノイズ量の程度の評価に基づいて行うことができる。判断部８０による判断動作については、第２実施形態にて詳述する。
【００７６】
例えば、判断部８０が、脈波信号ｄに含まれるノイズ量の程度が大きい（脈波信号ｄのきれいさの程度がきれい：目立つノイズ少）と判断したとき、第２処理部７０による拍動呈示スペクトルの捕捉処理が実行される。この場合、例えば、第１捕捉処理部７２が捕捉処理を実行する。また、脈波信号ｄに含まれるノイズ量の程度が中程度（脈波信号ｄのきれいさの程度がまあまあ：外乱中程度）の場合には、第２捕捉処理部７４が捕捉処理を実行する。第１捕捉処理ならびに第２捕捉処理の内容の一例については、図２および図３を用いて後述する。
【００７７】
第１処理部６０または第２処理部７０が拍動呈示スペクトルの特定、捕捉に成功すると、脈拍数算出部９０は、被検体の脈拍数を算出する。脈拍数は、例えば、拍動呈示スペクトルの周波数に基づいて算出する。求められた脈拍数の情報は、脈拍数算出部９０から表示処理部９２を経由して表示部９４に供給され、例えば、脈波数を示す数値が、表示部９４によって表示される。なお、脈拍数ではなく、検出した脈の、時間軸上における変化を信号波形やグラフの形式で表示してもよい。
【００７８】
（図１の構成によって得られる効果等）
第１処理部６０は、フィルター後信号ｅのＦＦＴ結果ｇ１および体動信号ｆのＦＦＴ結果ｇ２に基づいて、拍動信号を示す拍動呈示スペクトルを特定する。図１の拍動検出装置は、第２処理部７０を備える。第２処理部７０は、脈波信号ｄに外乱ノイズが多く含まれ、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような場合であっても、第２処理部７０によって拍動呈示スペクトルの捕捉を試みる。これによって拍動検出の可能性が向上する。この第２処理部７０は、フィルター前の脈波信号ｄに基づいて拍動呈示スペクトルの捕捉を試みる。フィルター前の脈波信号ｄに基づいた拍動呈示スペクトルの捕捉処理を試みることによって、拍動信号成分を検出することができる可能性（拍動呈示スペクトルの捕捉処理による救済の効果）が生じる。
【００７９】
但し、フィルター前の脈波信号ｄは、ノイズ対策がまったく行われていない原信号であり、例えば、過大なノイズしか存在しない信号である場合もあり得る。
【００８０】
そこで、図１の拍動検出装置１００では、第２処理部７０による拍動呈示スペクトルの捕捉処理を試みる前に、判断部８０が、拍動呈示スペクトル捕捉処理の実行の可否を判断することとした。判断部８０は、例えば、脈波信号ｄに含まれるノイズが多い状態である場合や、被検体（人ならびに動物を含む）の運動状態が急激に変化している場合、所定時間前の過去の検出時点において拍動呈示スペクトルの検出に失敗している場合等においては、第２処理部７０による拍動呈示スペクトルの捕捉処理は不可と判断する。具体例については後述する。
【００８１】
したがって、図１の拍動検出装置１００によれば、拍動検出性能を高めることができる。また、例えば、脈波センサーを手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着するタイプの拍動検出装置において、拍動信号を表す拍動呈示スペクトルの特定性能を向上させることができる。
【００８２】
（第１処理部による第１処理の内容）
第１処理部６０による第１処理は、フィルター後信号ｅのＦＦＴ結果ｇ１ならびに体動信号ｆのＦＦＴ結果ｇ２に基づいて実行することができる。第１処理は、脈波信号ｄに含まれる拍動信号成分の信号値（振幅値）が大きい場合に有効である。
【００８３】
第１処理は、例えば、以下の手順により行われる。ピーク順ソート部６２は、フィルター後信号ｅ１のＦＦＴ結果ｇ１に基づいて、スペクトルを信号値が大きい順にソーティングする。例えば、上位１０本のスペクトルのいずれかが拍動信号を示す拍動呈示スペクトルである可能性が高い。
【００８４】
次に、相関判断部６４が、相関判断を実行する。周波数解析部５０による周波数解析が４秒毎に行われる場合、相関判断部６４は、所定時間前の過去（つまり直近の過去）である４秒前の検出時に検出された拍動スペクトルを基準として、今回の検出（現在の検出期間における検出）に係る上位１０本のスペクトルについて相関の有無を判断する。これにより、拍動呈示スペクトルと推定される１本の候補が抽出される。
【００８５】
次に、拍動／ノイズ分離部６６は、抽出された候補について、体動を示す体動呈示スペクト等のノイズとの分離が可能かを判断する。例えば、抽出された候補の周波数が、体動スペクトルの周波数とは明らかに異なる場合には分離可能と判断される。
【００８６】
また、希に、抽出された候補の周波数と体動スペクトルの周波数とが一致する場合もあるが、この場合でも、例えば、所定時間前の過去および現在の体動信号の相関性等を考慮して、その候補を、体動呈示スペクトルと区別できる場合がある。
【００８７】
拍動呈示スペクトル特定部６８は、体動呈示スペクトルとは異なることが明らかとなった候補を、拍動呈示スペクトルとして特定することができる。
【００８８】
（第２処理部による第２処理の内容）
上述のとおり、第２処理部７０による第２処理（拍動呈示スペクトルの捕捉処理）は、例えば、第１処理部６０による拍動呈示スペクトルの特定ができなかった場合に行ってもよい。すなわち、判断部８０による拍動呈示スペクトル捕捉処理の実行の可否判断は、第１処理部６０が、拍動呈示スペクトルを特定することができなかった場合に行ってもよい。
【００８９】
第１処理部６０と第２処理部７０とを並列に動作させることも可能であるが、第１処理部６０によって拍動呈示スペクトルが特定できた場合には、第２処理部７０による拍動呈示スペクトルの捕捉は不要である。第２処理部７０の無駄な動作をなくすという観点から、第１処理部６０による拍動呈示スペクトルの特定ができないときに、判断部８０による拍動呈示スペクトルの捕捉処理の可否判断を開始させるのが好ましい。
【００９０】
第２処理部７０による第２処理（拍動呈示スペクトルの捕捉処理）は、第１捕捉処理部７２による第１捕捉処理と、第２捕捉処理部７４による第２捕捉処理と、に大別される。上述のとおり、第１捕捉処理部７２は、例えば、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルである場合（例えば、ノイズ量の程度が小さく、脈波信号がきれいであるとき）に動作する。また、第２捕捉処理部７４は、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルよりも低い第２許容レベルである場合（例えば、ノイズ量の程度が第１許容レベルを基準とする場合よりも大きく、脈波信号のきれいさが中程度のとき）に動作する。
【００９１】
脈波信号に含まれるノイズ量の程度に応じて、拍動呈示スペクトルの捕捉処理の内容を切り替えることによって、脈波信号の状態に応じた適切な拍動呈示スペクトルの捕捉処理を実行することができる。
【００９２】
なお、以下の説明では、「拍動基線」ならびに「体動基線」という用語を使用する場合がある。拍動基線は、周波数分解部４０の周波数分解処理によって得られる周波数スペクトルにおける、被検体の拍動由来の拍動成分と推定されるスペクトルである。体動基線は、周波数分解部４０の周波数分解処理によって得られる周波数スペクトルにおける、被検体の体動由来の体動成分と推定されるスペクトルである。
【００９３】
また、「拍動基線（広義）」には、今回の（つまり、現在の検出期間での）周波数解析処理に使用される脈波信号の周波数スペクトルにおける、被検体の拍動由来の成分と推定される「第１拍動基線」と、前回の（つまり所定時間前の）周波数解析処理に使用された脈波信号またはフィルター後信号の周波数スペクトルに基づいて、拍動呈示スペクトルとして特定された「第２拍動基線」とが含まれる。
【００９４】
第１拍動基線は、今回の検出処理に係る拍動基線であり、「現在の拍動基線あるいは今回の拍動基線」と称する場合がある。また、第２拍動基線は、前回（所定時間前の過去）の検出処理に係る拍動基線であり、「所定時間前の過去の拍動基線あるいは前回の拍動基線」と称する場合がある。これらは相関性の有無の判断において使用される。
【００９５】
また、体動基線（広義）は、体動センサー２０から出力される体動信号ｆの、現在ならびに過去の周波数分解結果に基づいて、体動成分と推定されるスペクトルの総称である。外乱ノイズがある程度存在する場合、例えば、被検体が運動している場合等では、体動基線には、異なる周波数のｍ本（ｍは１以上の整数）の体動基線が含まれることが多い。
【００９６】
第ｎ体動基線は、異なる周波数のｍ本（ｎ、ｍは１以上の整数）の体動基線のうち、一番周波数の低い基線を第１番目としそのｎ倍の周波数に相当する基線を第ｎ体動基線と定義する。例えば、第１波（基本波）の拍動基線が、第ｎ体動基線（１≦ｎ≦ｍ）の位置とは異なる位置にある場合には、第２処理部７０は、第１拍動基線が拍動呈示スペクトルであると判断する。この場合、例えば、ｎが２〜４である場合の体動基線は、ｎが１である場合の体動基線よりも信号値が小さく、外乱ノイズ中に潜んでいる場合があるため、第ｎ体動基線（例えばｎ＝１〜４）との区別を試みることは、拍動呈示スペクトルの捕捉精度を高める点で有効である。
【００９７】
（第１捕捉処理）
まず、第１捕捉処理部７２による第１捕捉処理について、図２（Ａ）および図２（Ｂ）を用いて説明する。図２（Ａ）および図２（Ｂ）は、第１捕捉処理部による第１捕捉処理の一例を示す図である。図２（Ａ）および図２（Ｂ）には、０〜４Ｈｚまでの脈波信号をＦＦＴして得られる周波数スペクトルが示されている。
【００９８】
上述のとおり、第１捕捉処理部７２は、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルである場合（例えば、ノイズ量の程度が小さく、脈波信号がきれいであるとき）に動作する。ここで、脈波信号が「きれい」である場合とは、例えば、広範囲の周波数域にわたってスペクトルのレベルが総じて小さく、その中に、信号値の大きな目立つスペクトル(一般的には、検出の目的とするスペクトルの候補となり得る)が含まれており、その目立つスペクトルの数がそれほど多くはない、というような場合である。一般的には、このような場合、目立つ複数本のスペクトルを抽出し、過去の拍動成分との相関性や、ノイズとの峻別性等に基づく判断を行うことが有効とされている。
【００９９】
但し、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に著しく低い場合等には、一見、信号がきれいな状態であっても、第１処理部６０による拍動／ノイズ分離処理を含む信号処理（参照符号６２〜６６のブロックによってなされる一連の信号処理）では、拍動信号を特定できない場合がある。
【０１００】
また、一見、信号がきれいに見える場合としては、拍動信号が存在せず、全部がノイズであるという場合がある。つまり、スペクトル全体の形状からみると、一見、きれいな信号に見えるが、実際は、目的とする信号成分はノイズに埋もれてほとんど認識することができない状態（つまり、全部がノイズであって目的とする拍動信号がない状態）である場合がある。このような場合、検出の目的である拍動信号が無いのであるから、外乱ノイズと拍動信号との分離を試みても効果がないと言える。
【０１０１】
そこで、第１捕捉処理部７２では、第１捕捉処理を実行する。この第１捕捉処理において、拍動呈示スペクトルの捕捉のためには、上述した第１拍動基線（今回の検出に係る拍動基線）と第２拍動基線（所定時間前の過去の検出に係る拍動基線）との間の相関性に着目することが有効である。
【０１０２】
上述のとおり、第１拍動基線は、今回の周波数解析処理に使用される、「脈波信号」の周波数スペクトルにおける被検体の拍動由来の拍動成分と推定されるスペクトルである。第２拍動基線とは、前回の周波数解析処理に使用された、「脈波信号またはフィルター後信号」の周波数スペクトルに基づいて、拍動呈示スペクトルとして特定された拍動基線である。
【０１０３】
前回検出時の周波数スペクトルにおける第１拍動基線と、今回検出時の周波数スペクトルにおける第２拍動基線との相関を重視して判断したとき、相関性が有るとみることができれば、第１捕捉処理部７２は、第１拍動基線は、検出目的の拍動呈示スペクトルであると判断する。よって、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような場合であっても、第１拍動基線を拍動呈示スペクトルとして捕捉することができる。この場合、第１拍動基線と第２拍動基線との間の相関によって判断していることから、仮に比較的大きなノイズがあったとしても、拍動呈示スペクトルを捕捉することができる。
【０１０４】
また、外乱ノイズ（非周期的ノイズ）との分離は困難であっても、体動（周期的な腕振り等）に由来する周期的な体動信号との区別は可能であることが多い。よって、相関性有りと判断され、かつ、周期的な体動信号と区別できることを、拍動呈示スペクトルの捕捉の条件としてもよい。図２（Ａ）および図２（Ｂ）の例では、この判断基準を採用して、拍動呈示スペクトルを捕捉する。
【０１０５】
図２（Ａ）において、Ｘ１が体動基線であり、Ｘ２が拍動基線であるとする。拍動基線Ｘ２は、体動基線Ｘ１よりも十分に小さく、ノイズに埋もれた状態となっている。
【０１０６】
この場合において、図２（Ｂ）に示すように、まず、０〜４Ｈｚの周波数スペクトルの中から、信号値の大きいもの順に、例えば１０本のスペクトルを抽出する。ここでは、信号値の大きいスペクトルから順に、スペクトルＬ１〜Ｌ１０とする。抽出されたスペクトルの中に、拍動呈示スペクトルが存在する可能が高い。但し、信号値に着目した判断を実行すると、信号値の大きなノイズのスペクトルを拍動呈示スペクトルとして誤検出する可能性が高い。よって、所定時間前の過去の検出に係る拍動基線（第２拍動基線）と、今回の検出に係る第１拍動基線との相関性を重視して、拍動呈示スペクトルを特定する。
【０１０７】
図２（Ｂ）において、第２拍動基線（前回の拍動基線）はＬ０である。ここでは、実際の拍動呈示スペクトルはＬ６であるとする。Ｌ１〜Ｌ１０が、拍動呈示スペクトルの候補であり、各スペクトルＬ１〜Ｌ１０が、今回の検出に係る第１拍動基線（一応、拍動成分と推定できる基線）となる。
【０１０８】
次に、所定時間前の過去（４秒前）に拍動基線と推定され、かつ拍動呈示スペクトルとして特定された拍動基線（第２拍動基線）Ｌ０と、各スペクトルＬ１〜Ｌ１０との周波数軸上における距離を算出する。
【０１０９】
第１拍動基線Ｌ０に最も近い第１拍動基線はＬ６である。ここで、周波数軸上におけるＬ０とＬ６の距離に着目する。Ｌ０とＬ６の距離は「２」である。「距離２」とは、周波数軸上で、第２拍動基線Ｌ０の位置（第２位置）を基準として、第１拍動基線Ｌ６の位置（第１位置）が１本のスペクトルを隔てて隣り合っている場合の距離である。例えば、「基準のスペクトルに対して距離２以内のスペクトル」というとき、基準のスペクトルに対して、２本のスペクトル以内で、低周波数側（左）あるいは高周波数側（右）において隣り合っているスペクトル、という意味である。第１捕捉処理部７２は、基準となる第２拍動基線Ｌ０と、今回の検出に係る第１拍動基線Ｌ６の距離が「距離２以内（＝距離０，距離１，距離２）」であるとき、第１拍動基線Ｌ６を拍動呈示スペクトルとしてもよい。
【０１１０】
例えば、スペクトル１本分が、脈拍数３．７５拍分に相当する。被検体である人が定常的に同じようなリズムで運動しているとき、４秒間に、脈拍数が“２０”以上変化することは極めて希である。ここで、周波数軸上での距離が３であるということ（周波数軸上で、スペクトル（基線）が３本分離れること）は、被検体の脈波数が“１０”（上記の脈拍数“２０”の半分）変化するのに相当する。
【０１１１】
本実施形態では、第１捕捉処理部７２は、４秒間における脈拍数の“１０”の変化は発生しないものとして、判断処理を実行する。上記の脈拍数２０に対して半分の値１０を基準としたのは、判定処理の安全度を加味したためである。
【０１１２】
したがって、図２（Ｂ）の例では、第１捕捉処理部７２は、ピーク順が６番のスペクトルである第１拍動基線Ｌ６を、拍動呈示スペクトルとして特定する。
【０１１３】
以上の手順は、第１捕捉処理の一例であり、変形が可能である。例えば、上述の手順によれば、信号ピークの上位１０本の全部について、第２拍動基線Ｌ２との間の距離判断やスペクトルの大きさの比に基づく判断が必要であるが、次に示す判断を行うことで、第１拍動基線を減らすことができる。すなわち、判断処理量を減らすことができる。
【０１１４】
ここで、体動センサー２０から出力される体動信号ｆのＦＦＴ結果に基づいて、スペクトルＬ１が体動基線（第ｎ体動基線）と推定できる場合、このスペクトルＬ１と、このスペクトルＬ１の、低周波数側（左）ならびに高周波数側（右）に隣り合う各１本のスペクトルＬ２，Ｌ３は、第１拍動基線（今回の検出に係る、拍動成分と推定される基線）から除外することができる。３本のスペクトルＬ１〜Ｌ３は、体動成分である可能性が極めて高いからである。
【０１１５】
これによって、第１拍動基線は、Ｌ４〜Ｌ１０の７本となる。この場合、第１拍動基線の数が減少した分、第２拍動基線Ｌ２との間の距離判断のための、第１捕捉処理部７２の負担が軽減されることになる。なお、後述する図６のステップＳＴ１０，ＳＴ１１では、体動基線（第ｎ体動基線）と、その体動基線の左右に位置する各１本のスペクトルを第１拍動基線から除外した後、第２拍動基線との間の距離判定を実行している。
【０１１６】
（第２捕捉処理）
第２捕捉処理部７４は、上述のとおり、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルよりも低い第２許容レベルである場合、例えば、ノイズ量の程度が第１許容レベルを基準とする場合よりも大きく、脈波信号のきれいさが中程度（まあまあきれい）のときに動作する。
【０１１７】
脈波信号ｄには、ノイズ少（きれい）の場合と比べてかなり目立つ非周期的なノイズ（外乱ノイズ）が多く含まれている。よって、拍動呈示スペクトルの捕捉のためには、拍動基線とその他の基線（外乱ノイズや体動ノイズに起因してスペクトル値が大きくなっている）とを区別するために別の判断処理、すなわち第２捕捉処理が必要となる。
【０１１８】
第２捕捉処理部７４は、第１処理部６０で行われる拍動／ノイズ分離処理（図１の例における、参照符号６２〜６６のブロックによって行われる処理）と同じ処理を実行することができる。但し、これだけでは、拍動信号と外乱ノイズとを十分には区別できない可能性がある。
【０１１９】
そこで、第２捕捉処理部は、さらに、今回の検出時の周波数スペクトルにおける第１拍動基線と、第ｎ体動基線との区別の可否に基づいて拍動呈示スペクトルの捕捉を試みる。
【０１２０】
上述のとおり、第ｎ体動基線は、異なる周波数のｍ本（ｎ、ｍは１以上の整数）の体動基線のうち、一番周波数の低い基線を第１番目としそのｎ倍の周波数に相当する基線を第ｎ体動基線と定義する。例えば、第１波（基本波）の拍動基線が、第ｎ体動基線（１≦ｎ≦ｍ）の位置とは異なる位置にある場合には、第２捕捉処理部は、第１拍動基線が拍動呈示スペクトルであると判断する。この場合、例えば、ｎが２〜４である場合の体動基線は、ｎが１である場合の体動基線よりも信号値が小さく、外乱ノイズ中に潜んでいる場合があるため、第ｎ体動基線（例えばｎ＝１〜４）との区別を試みることは、拍動呈示スペクトルの捕捉精度を高める点で有効である。
【０１２１】
このように、第２捕捉処理部７４による第２捕捉処理は、例えば、第１処理部で行われる拍動／ノイズ分離処理と同じ処理を、フィルター前の脈波信号のＦＦＴ結果に対して実行し、さらに、体動基線（第ｎ体動基線）との区別の可否を慎重に判断する処理ということができる。
【０１２２】
なお、第２捕捉処理においては、第１拍動基線と第２拍動基線との相関性も考慮してもよい。例えば、今回検出時の周波数スペクトルにおける第１拍動基線と、第２拍動基線との間の周波数軸上における距離の観点からの判断を、拍動呈示スペクトルの捕捉の条件としてもよい。
【０１２３】
例えば、今回検出に係る第１拍動基線が、体動センサー出力の周波数分解結果に基づいて特定される体動基線（基本波成分ならびに高調波成分を含む）に重なっていないことを第１条件とし、今回検出に係る第１拍動基線が、周波数軸上で、第２拍動基線に対して所定距離内にあることを第２条件としてもよい。
【０１２４】
本実施形態では、第２捕捉処理部７４は、周波数軸上でスペクトル（基線）が３本分離れるということが発生しないものとして、拍動呈示スペクトルの呈示処理を行う。上述の第１条件と第２条件を満足する場合には、第２捕捉処理部は、第１拍動基線を拍動呈示スペクトルとして捕捉する。これによって、例えば、目立つノイズに引きずられることなく、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような拍動呈示スペクトルを特定できる可能性が高まる。
【０１２５】
図３（Ａ）および図３（Ｂ）は、第２捕捉処理部による第２捕捉処理の一例を示す図である。図３（Ａ）および図３（Ｂ）には、０〜４Ｈｚまでの脈波信号をＦＦＴして得られる周波数スペクトルが示されている。図３（Ａ）は第２捕捉処理によって救済が可能である例を示し、図３（Ｂ）は第２捕捉処理によっても救済することができない例を示す。
【０１２６】
第２捕捉処理部７４は、フィルター前の脈波信号ｄについて、ピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６による拍動／ノイズ分離処理を実行する。この後、捕捉部７６が、体動基線（第ｎ体動基線）との峻別（区別）の可否を判断する。捕捉部７６は、さらに、第１拍動基線と第２拍動基線との相関性を考慮してもよい。例えば、今回検出時の周波数スペクトルにおける第１拍動基線と、第２拍動基線との間の周波数軸上における距離の観点からの判断を行ってもよい。
【０１２７】
図３（Ａ）の例では、第２拍動基線がＬ１３であり、今回の検出に係る第１拍動基線がＬ１４である。また、第ｎ体動基線（ｎ＝２）をＬ１５とし、第ｎ体動基線（ｎ＝３）をＬ１６とする。
【０１２８】
図３（Ａ）の例では、第１拍動基線Ｌ１４は、第ｎ体動基線（ｎ＝２，３）Ｌ１５，Ｌ１６と重なってはいない。また、第２拍動基線Ｌ１３と、今回の検出に係る第１拍動基線Ｌ１４との間の、周波数軸上での距離は「距離２」である。よって、第１拍動基線Ｌ１４は、第２拍動基線Ｌ１３と相関があると判断できる。よって、第１拍動基線Ｌ１４を、拍動呈示スペクトルとして捕捉することができる。
【０１２９】
図３（Ｂ）の例では、第２拍動基線がＬ１７であり、今回の検出に係る第１拍動基線がＬ１８である。また、第ｎ体動基線（例えばｎ＝１）をＬ１９とする。
【０１３０】
図３（Ｂ）の例では、第１拍動基線Ｌ１８が、第ｎ体動基線（例えばｎ＝１）Ｌ１９と重なっている。つまり、第１拍動基線Ｌ１８は、第ｎ体動基線（例えばｎ＝１）Ｌ１９と区別がつかない。上述の第１条件を満足しないことから、第１拍動基線Ｌ１８を、拍動呈示スペクトルとして捕捉することができない。
【０１３１】
また、図３（Ｂ）の例では、第１拍動基線Ｌ１８は、第２拍動基線Ｌ１７に対する周波数軸上での距離が３である。すなわち、第１拍動基線Ｌ８は、２本のスペクトルを隔てて、第２拍動基線Ｌ２に隣り合っている。よって、第２条件（距離２以内である）を満足せず、この点でも、第１拍動基線Ｌ１８を、拍動呈示スペクトルとして捕捉することができない。この場合には、拍動呈示スペクトルを検出できない。
【０１３２】
（第２実施形態）
本実施形態では、拍動検出装置の構成と動作をより具体的に説明する。本実施形態では
特に、判断部８０の構成と動作について詳述する。
【０１３３】
図４は、本発明の拍動検出装置の他の例の構成を示す図である。図４において、図１と共通する部分には同じ参照符号を付し、その説明を省略している。
【０１３４】
図４に示される拍動検出装置１００は、脈波センサー１０と、脈波信号蓄積部（４秒分の脈波信号ｄのデータを蓄積する第１バッファメモリー１３および１６秒分の脈波信号ｄのデータを蓄積する第２バッファメモリー１５を有する）１２と、適応フィルター３２および体動成分除去フィルター３４を含むフィルター部３０と、体動センサー（加速度センサーやジャイロセンサー等）２０と、体動信号蓄積部２２と、周波数解析部５０と、脈拍数算出部９０と、表示処理部９２と、表示部９４と、を有する。
【０１３５】
上述のとおり、脈波センサー１０は、例えば、光電脈波センサー及びその原理に基づく脈波センサーである。脈波センサー１０は、拍動信号と、被検体（人や動物）の体動に由来する体動ノイズ信号を含むノイズ信号とが混在した脈波信号ｄを出力する。
【０１３６】
ここで、脈波信号ｄは、例えば、拍動信号成分（定常的成分あるいは周期的成分）と、体動ノイズ成分（定常的あるいは周期的成分）と、外乱ノイズ成分（衝撃ノイズ等の、非定常的あるいは非周期的成分）とを含む。
【０１３７】
脈波センサー１０から出力される脈波信号ｄの、４秒分の信号が、第１バッファメモリー１３に蓄積される。４秒分の脈波信号ｄは、４秒周期で、第２バッファメモリー１５に転送される。第２バッファメモリー１５はＦＩＦＯ（ファーストイン・ファーストアウト）メモリーであり、１６秒分の脈波信号は、４秒分ずつ更新される。１６秒分の脈波信号を蓄積するのは、周波数解析によって拍動成分を特定するとき、ある程度の時間幅で信号の推移を観測し、脈拍の時間的変化を含めて慎重に検討する必要があるからである。
【０１３８】
フィルター部３０は、入力信号に、定常的（周期的）な周波数成分とその他の非定常的（非周期的）な成分が含まれるときに、それらを分離して出力することのできる適応フィルターの一種である。フィルター部３０に含まれる適応フィルター３２によって、脈波信号ｄに含まれる定常的（周期的）な周波数成分とその他の非定常的（非周期的）な成分とを分離することができる。また、体動成分除去フィルター３４によって、例えば、脈波信号ｄに含まれる目立つ体動成分を除去・抑制することができる。
【０１３９】
周波数解析部５０は、信号分配部３９と、第１周波数分解部４０ａ〜第３周波数分解部４０ｃと、タイミング部４３と、第１処理部６０と、第２処理部７０（第１捕捉処理部７２および第２捕捉処理部７４とを含む）と、判断部８０を含む。
【０１４０】
信号分配部３９は、フィルター後信号ｅを第１周波数分解部４０ａに供給し、体動信号ｆを第２周波数分解部４０ｂに供給し、フィルター前の脈波信号ｄを第３周波数分解部４０ｃに供給する。
【０１４１】
タイミング部４３は、第１処理部６０ならびに第２処理部７０の動作タイミングを制御する。なお、第１処理部６０ならびに第２処理部７０の動作は、前掲の実施形態にて説明したとおりである。
【０１４２】
判断部８０は、周波数解析部（ここでは、第３周波数分解部４０ｃ）が生成した周波数スペクトルに基づいて、脈波信号ｄのきれいさに関する判断指標を取得する判断指標取得部８１と、判断指標を用いて、脈波信号ｄに含まれるノイズ量の程度を評価する評価部８３と、被検体の運動状態を判断する運動状態判断部８５と、前回の拍動検出の結果を格納している拍動検出履歴格納部８７と、捕捉処理の実行可否判断部８９と、を有する。
【０１４３】
捕捉処理の実行可否判断部８９は、例えば、評価部８３による評価に基づいて、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断してもよい（第１判断例）。
【０１４４】
また、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、評価部８３による評価が許容範囲にあることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい（第２判断例）。
【０１４５】
また、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、評価部８３による評価が許容範囲にあり、かつ、運動状態判断部８５によって被検体が定常状態にあると判断されることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。
【０１４６】
なお、第１処理部６０と第２処理部７０とを並列に動作させることも可能であるが、第１処理部６０によって拍動呈示スペクトルが特定できた場合には、第２処理部７０による拍動呈示スペクトルの捕捉は不要である。そこで、第２処理部７０の無駄な動作をなくすために、第１処理部６０によって拍動呈示スペクトルの特定ができなかったときに、判断部８０による拍動呈示スペクトルの捕捉処理の可否判断を開始させてもよい。
【０１４７】
（第１判断例）
判断部８０に含まれる捕捉処理の実行可否判断部８９は、脈波信号ｄに含まれるノイズ量の程度の評価に基づいて、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断する。なお、脈波信号ｄに含まれるノイズ量の程度の判断指標の具体例、ならびに評価の例については、図５を用いて後述する。
【０１４８】
捕捉処理の実行可否判断部８９は、脈波信号ｄに含まれるノイズ量の程度の評価を参照して、例えば、脈波信号に混入されるノイズ量（例えば外乱ノイズ成分のノイズ量）が所定の基準以下であると判断される場合に、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。換言すれば、脈波信号ｄが、許容レベル以上にきれいであることを第２処理部７０による拍動呈示スペクトル捕捉処理の実行の条件としてもよい。
【０１４９】
第１判断例によれば、外乱ノイズ（非周期的ノイズ）が多い状態では、このような場合には、拍動呈示スペクトルの捕捉処理が実行されない。よって、拍動呈示スペクトルの誤検出の可能性が低減される。第１判断例によれば、脈波信号ｄのノイズ状態を慎重に判断しつつ、第２処理部７０による拍動呈示スペクトルの捕捉処理の可否を判断してもよい。
【０１５０】
（第２判断例）
判断部８０に含まれる捕捉処理の実行可否判断部８９は、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、評価部８３による評価が許容範囲にあることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。
【０１５１】
第２判断例では、周波数解析部（ここでは、第１周波数分解部４０ａまたは第３周波数分解部４０ｃ）による前回の周波数解析によって、拍動呈示スペクトルが検出（特定）されていること、ならびに、評価部８３による、脈波信号ｄに含まれるノイズ量の程度の評価が許容範囲にあることを条件として、第２処理部７０による拍動呈示スペクトルの捕捉処理を実行する。
【０１５２】
上述の第１判断例では、ノイズについての基準によって、第２処理部７０による拍動呈示スペクトルの捕捉処理の可否を判断していた。第２判断例では、前回の拍動信号の検出に成功しているという条件を満足するか否かも判断する。前回の拍動検出の結果は、拍動検出履歴格納部８７に格納されている。捕捉処理の実行可否判断部８９は、拍動検出履歴格納部８７から前回の結果を読み出して参照することによって、前回の拍動信号の検出の成功／不成功を確認することができる。
【０１５３】
上述のとおり、周波数分解部４０（４０ａ〜４０ｃ）は、一定時間毎（例えば４秒毎）に周波数解析処理を行って拍動呈示スペクトルを特定する。第２判断例では、前回の検出に成功していない場合には、第２処理部７０による拍動呈示スペクトルの捕捉処理は実行されない。
【０１５４】
上述のとおり、拍動呈示スペクトルの検出には、所定時間前の過去の拍動信号成分と現在の拍動信号成分との相関性に着目することが有効である。よって、前回の拍動信号の検出に成功しているか否かを条件としてもよい。第２判断例では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【０１５５】
（第３判断例）
第３判断例では、判断部８０に含まれる捕捉処理の実行可否判断部８９は、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、評価部８３による評価が許容範囲にあり、かつ、運動状態判断部８５によって被検体が定常状態にあると判断されることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。
【０１５６】
被検体の運動状態が急激に変化している場合には、脈拍等の変化が急である。よって、運動状態判断部８５によって被検体が定常状態にあると判断される、という条件によって、第２処理部７０による拍動呈示スペクトルの捕捉処理の実行の可否を判断してもよい。
【０１５７】
被検体が定常状態にある場合とは、例えば、被検体が一定の速度で歩行している場合のように、体動センサー２０から得られる体動信号の状態が安定している状態、あるいは、被検体が安静状態にある場合等である。例えば、被検体が静止している状態から歩行状態に移行したとき、当初は加速過程にあることから、体動信号のスペクトル（体動スペクトル）は一定しない。被検体の歩行状態が安定する（例えば、一定の速度で歩行する状態となる）と、体動信号のスペクトルが安定し、時間が経過しても、体動スペクトルに大きな変化が生じない状態となる。
【０１５８】
運動状態判断部８５は、例えば、２０秒間（１回の測定が４秒であるとすると、その５倍の期間）にわたって、体動スペクトルに大きな変化が見られない場合に、被検体が定常状態にあると判断する。但し、これは一例であり、これに限定されるものではない。例えば、現在ならびに過去の脈波信号のスペクトルの安定性を判断基準とすることもできる。また、後述するｒ５，ｒ１０という指標によっても、被検体の運動状態を判断することができる。
【０１５９】
運動状態判断部８５による被検体の運動状態の判断結果（定常状態／非定常状態の判断結果）は、捕捉処理の実行可否判断部８９に供給される。第３判断例では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【０１６０】
（脈波信号に含まれるノイズ量の程度（ノイズ量の程度）の判断について）
以下、脈波信号ｄに含まれるノイズ量の程度（ノイズ量の程度）の判断の例について説明する。図５（Ａ）〜図５（Ｃ）は、脈波信号に含まれるノイズ量の程度の判断動作を説明するための脈波信号の波形および周波数スペクトルの一例を示す図である。
【０１６１】
図５（Ａ）〜図５（Ｃ）において、上側には、１６秒間のＦＦＴ前の脈波信号ｄの信号波形が示されている。横軸は時間を示し、縦軸は信号の振幅を示す。また、下側には、０から４Ｈｚの周波数帯域における周波数スペクトルが示されている。横軸は周波数を示し、縦軸はスペクトル値を示す。
【０１６２】
ここで、図５（Ａ）は、ノイズ少（脈波信号のきれいさの程度がきれい）の場合における脈波信号ｄの波形と周波数スペクトルを示し、図５（Ｂ）は、ノイズが中程度（脈波信号のきれいさの程度がまあまあきれい）の場合における脈波信号の波形と周波数スペクトルを示しており、図５（Ｃ）は、脈波信号ｄに多くのノイズが含まれる場合（脈波信号のきれいさの程度がノイジー）における脈波信号ｄの波形と周波数スペクトルを示している。図５（Ａ）〜図５（Ｃ）の各々の比較から明らかなように、脈波信号ｄの波形と周波数スペクトルとは密接に関連しており、脈波信号の波形に対応して、周波数スペクトルの分布状態やスペクトル値が変化する。よって、ＦＦＴによって得られる周波数スペクトルに基づいて、脈波信号ｄにする重畳するノイズの状態（ノイズ量の程度）を推定することが可能である。
【０１６３】
図５の例では、ノイズ量の程度の推定のための指標として、例えば、主要な周波数スペクトルのスペクトル値の比（つまり、基線の高さの比）を用いることができる。具体的には、ｒ５およびｒ１０という指標を用いる（ただし、一例であり、他の統計的指標、例えば、標準偏差等を用いてもよい）。ここで、ｒ５とは、１６秒分の脈波信号の周波数スペクトルの中から、ピーク値の大きさの順に５本のスペクトルを並べたとき（つまり、ソーティングしたとき）、第１番目のスペクトルのスペクトル値（パワー）を分母とし、第５番目のスペクトルのスペクトル値（パワー）を分子とすることによって得られる指標である。
【０１６４】
また、ｒ１０とは、１６秒分の脈波信号の周波数スペクトルの中から、ピーク値の大きさの順に１０本のスペクトルを並べたとき（つまりソーティングしたとき）、第１番目のスペクトルのスペクトル値（パワー）を分母とし、第１０番目のスペクトルのスペクトル値（パワー）を分子とすることによって得られる指標である。指標ｒ５，ｒ１０は、判断部８０に含まれる判断指標取得部８１によって取得される。
【０１６５】
ここでは、一例として、ｒ５＜０．５かつｒ１０＜０．２のときをノイズ少（きれい）とし、ｒ５＞０．７かつｒ１０＞０．５のときをノイズ多（ノイジー）とし、上記いずれでもない場合をノイズが中程度（まあまあ）とする。
【０１６６】
図５（Ａ）の例では、ｒ５＝０．１４かつｒ１０＝０．０８であることから、ノイズ少（きれい）と判断される。また、図５（Ｂ）の例では、ｒ５＝０．５６かつｒ１０＝０．３５であることから、ノイズが中程度（まあまあ）と判断される。図５（Ｃ）の例では、ｒ５＝０．８２かつｒ１０＝０．６２であることから、ノイズ多（ノイジー）と判断される。
【０１６７】
このような判断処理が、判断部８０に含まれる評価部８３によって実行される。なお、被検体の運動状態も、上述の指標ｒ５、ｒ１０を用いて推定することもできる。被検体の運動状態によって脈波信号ｄの波形が変化すると、その変化は、周波数スペクトルの変化となって現れ、周波数スペクトルの変化は、指標ｒ５，ｒ１０に反映されるからである。被検体の運動状態は、上述のとおり、運動状態判断部８５によって実行される。
【０１６８】
（拍動呈示スペクトルの捕捉処理のフローの例について）
図６は、第２処理部による、拍動呈示スペクトルの捕捉処理（第２処理）の手順の一例を示すフローチャートである。
【０１６９】
上述のとおり、例えば、第１処理部６０による拍動呈示スペクトルの特定（第１処理）に失敗した場合に、第２処理部７０による拍動呈示スペクトルの捕捉処理（第２処理）を実行することができる。
【０１７０】
まず、判断部８０は、前回拍動呈示スペクトル（拍動基線）を特定できたかを判断する（ステップＳＴ１）。特定できたとき（ステップＳＴ１：Ｙ）のときは、被検体が定常状態であるかを判断する（ステップＳＴ２）。定常状態であるときは（ステップＳＴ２：Ｙ）のときは、脈波信号に含まれるノイズ量の程度（脈波信号のきれいさの程度）を判断する（ステップＳＴ３）。具体的には、きれいさの程度が、きれい（外乱ノイズ少：第１許容レベル）、まあまあ（外乱ノイズ中程：第２許容レベル）、ノイジー（許容レベル未満）のいずれであるかを判断する。
【０１７１】
ステップＳＴ１でＮのとき、あるいは、ステップＳＴ２でＮのときは、拍動呈示スペクトルの検出失敗とする（ステップＳＴ９）。また、ステップＳＴ３で、ノイジーと判断されたときも、拍動呈示スペクトルの検出失敗とする（ステップＳＴ９）。
【０１７２】
ステップＳＴ３で、外乱ノイズ少（脈波信号のきれいさの程度がきれい：第１許容レベル）と判断された場合には、第１捕捉処理部７２が、第１捕捉処理（図６の左側において、破線で囲んで示されるステップＳＴ１０〜ＳＴ１６）を実行する。ステップＳＴ３で、外乱ノイズ中程度（脈波信号のレベルがまあまあきれい：第２許容レベル）と判断された場合には、第２捕捉処理部７４が、第２捕捉処理（図６の右下において、破線で囲んで示されるステップＳＴ４〜ＳＴ７）を実行する。
【０１７３】
第１捕捉処理が実行される場合、まず、第１捕捉処理部７２は、ピーク順ソート部６２でソートされた結果の順位が１０位以内のスペクトルを抽出し、その１０本の中から、体動基線（複数本あるときは、好ましくは各線）と重なっていない、あるいは、各体動基線の、低周波数側（左）あるいは高周波数側（右）に隣り合う１本を除くという条件を満足するスペクトルをリストアップする（ステップＳＴ１０）。リストアップされたスペクトルを、第１拍動基線（今回の検出に係る拍動基線）とすることができる。
【０１７４】
次に、第１捕捉処理部７２は、リストアップした第１拍動基線毎に、第２拍動基線（前回の検出に係る拍動基線）との間の周波数軸上における距離を算出する（ステップＳＴ１１）。
【０１７５】
ステップＳＴ１２では、距離０の第１拍動基線があるかを検出し、有れば（ステップＳＴ１２：Ｙ）、拍動呈示スペクトルの検出に成功する（ステップＳＴ１３）。
【０１７６】
ステップＳＴ１４では、距離１の第１拍動基線があるかを検出する。第１拍動基線が有る場合（ステップＳＴ１４：Ｙ）であって、第２拍動基線の左右の両側に各１本が存在する場合には、信号値が大きい方の第１拍動基線を選択する（ステップＳＴ１５）。これによって、拍動呈示スペクトルの検出に成功する（ステップＳＴ１３）。
【０１７７】
ステップＳＴ１６では、距離２の第１拍動基線があるかを検出し、有れば（ステップＳＴ１６：Ｙ）、拍動呈示スペクトルの検出に成功する（ステップＳＴ１３）。
【０１７８】
ステップＳＴ１２，ＳＴ１４，ＳＴ１６のいずれかにおいてＮの場合は、拍動呈示スペクトルを検出せずと判断する。（ステップＳＴ１７）。
【０１７９】
次に、第１捕捉処理の手順について説明する。上述のとおり、ステップＳＴ３で、脈波信号に含まれるノイズ量の程度が中程度（脈波信号のきれいさの程度がまあまあきれい：第２許容レベル）と判断された場合には、第２捕捉処理部７４が、第２捕捉処理（図６の右下において、破線で囲んで示されるステップＳＴ４〜ＳＴ７）を実行する。
【０１８０】
第２捕捉処理部７４は、まず、ステップＳＴ４にて、フィルター後信号の最大スペクトル値（最大信号値）が所定値（ここでは相対値１０００）以下であるかを判定する。このステップＳＴ４は必須ではないが、実行するのが好ましい。第２捕捉処理の趣旨は、例えば、総じて信号値が大きくなく、かつ全体のスペクトルの形状として凹凸が少ない場合に、信号レベルが相対的に小さい拍動呈示スペクトルを見つけ出すということにある。よって、例えば、最大信号値が１０００を超えるスペクトルが含まれる場合には、上述の趣旨に反することになる。つまり、ステップＳＴ４でＮの場合には、拍動呈示スペクトルの検出に失敗する（ステップＳＴ９）。
【０１８１】
ステップＳＴ４でＹの場合には、フィルター前の脈波信号に対して第１処理（第１処理部６０で行われる拍動／ノイズ分離のための処理と同等の処理：図１の参照符号６２〜６４のブロックで実行される処理）を実行する（ステップＳＴ５）。
【０１８２】
次に、拍動基線（拍動成分と推定される第１拍動基線）が存在するか否かを判定する（ステップＳＴ６）。ステップＳＴ６でＮの場合は、拍動呈示スペクトルの検出に失敗する（ステップＳＴ９）。
【０１８３】
ステップＳＴ６でＹの場合は、ステップＳＴ７の処理を実行する。ステップＳＴ７では、今回の検出に係る拍動基線（第１拍動基線）が、第ｎ体動基線（例えば、ｎ＝１〜４の各基線）に重なっていないかを判断し（第１条件の判断）、さらに、第１拍動基線は第２拍動基線に対して２本以内で隣り合っている（換言すれば、３本以上離れていない）か、を判断する（第２条件の判断）。これらの２つの条件を満たす場合（Ｙの場合）、拍動呈示スペクトルの検出に成功する（ステップＳＴ８）。Ｎの場合は、拍動呈示スペクトルの検出に失敗する（ステップＳＴ９）。
【０１８４】
（第２捕捉処理の適用例）
図７（Ａ）および図７（Ｂ）は、第２捕捉処理の適用例を示す図である。図７（Ａ）および図７（Ｂ）の例では、体動ノイズが多く混入し、外乱ノイズの混入は少ない。
【０１８５】
図７（Ａ）の上側には、１６秒分のフィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図７（Ｂ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図７（Ｂ）の例では、脈波信号に含まれるノイズ量の程度は中程度であり、脈波信号ｄのきれいさの程度は、「まあまあきれい」と評価される。
【０１８６】
図７（Ａ）の例では、フィルター後信号に基づく周波数分解によって得られるスペクトルに含まれる拍動基線Ｓ２は小さく、体動ノイズＮ３と区別する困難性は、拍動基線Ｓ２が大きい場合に比べて高くなる。しかし、図７（Ｂ）に示すように、フィルター前の脈波信号周波数分解によって得られるスペクトルに含まれる拍動基線Ｓ１（第１拍動基線Ｓ１）はかなり大きく、かつ、拍動基線Ｓ１の周波数軸上における位置は、過大な体動ノイズＮ１，Ｎ２の位置から離れている。
【０１８７】
図７（Ｂ）の例では、第２捕捉処理部７４による第２捕捉処理が実行される。すなわち、フィルター前の脈波信号に対して、フィルター後信号に対する拍動／ノイズ分離処理と同様の信号処理を施し、第１拍動基線Ｓ１を、複数本の第ｎ体動基線Ｎ１，Ｎ２（ｎ＝１，２に相当）から分離する。図７（Ｂ）の場合、第１拍動基線Ｓ１は、第ｎ体動基線Ｎ１，Ｎ２と重なっておらず、区別することができる。よって、第１拍動基線Ｓ１を、拍動呈示スペクトルとして捕捉（特定）することができる。また、前回の検出時において、脈波信号に基づいて第２拍動基線の特定に成功しているのならば、その第２拍動基線との相関性を併せて判断すると、判断の精度が向上する。
【０１８８】
（第１捕捉処理の適用例）
図８（Ａ）〜図８（Ｃ）は、第１捕捉処理の適用例を示す図である。図８（Ａ）〜図８（Ｃ）は、体動ノイズの混入は多くはないが、外乱ノイズの混入がある例を示す。この場合、脈波信号に含まれるノイズ量の程度は見かけ上「小さい」、つまり、脈波信号ｄのきれいさの程度は見かけ上「きれい」と評価される。
【０１８９】
図８（Ａ）の上側には、前回の１６秒分のフィルター後信号が示されており、下側には前回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図８（Ｂ）の上側には、今回の１６秒分のフィルター後信号が示されており、下側には今回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図８（Ｃ）の上側には、今回の１６秒分の脈波信号が示されており、下側には今回の脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【０１９０】
図８（Ａ）の例では、前回のフィルター後信号に基づく周波数分解によって得られるスペクトルに含まれる拍動基線（第２拍動基線）Ｓ３の信号値は小さいものの、周囲のノイズよりも信号値は大きく、第１処理部６０による第１処理によって、拍動呈示スペクトルが特定されている。
【０１９１】
但し、図８（Ｂ）の例では、今回のフィルター後信号に基づく周波数分解によって得られるスペクトルでは、拍動基線（第１拍動基線）の信号レベルは、他の基線の信号レベルよりも相対的に低い。
【０１９２】
この場合には、第１捕捉処理部７２による第１捕捉処理が実行される。第１捕捉処理を実行した結果を示すのが図８（Ｃ）の例である。図８（Ｃ）の例では、上述のとおり、脈波信号に含まれるノイズ量の程度は、目立つ体動基線等がないことから、一応、小さい（脈波信号ｄのきれいさの程度はきれい）と判断される。但し、外乱ノイズが多いため、第１拍動基線Ｓ４を検出する難易度は、外乱ノイズが少ない場合に比べて高くなる。
【０１９３】
この場合には、例えば、第１捕捉処理部７２が、第２拍動基線Ｓ３と、今回の検出に係る第１拍動基線Ｓ４との相関性を重視した判断を実行し、その相関性に基づいて、第１拍動基線Ｓ４を拍動呈示スペクトルとして特定する。このとき、第１拍動基線Ｓ４が、体動基線（例えば複数本の第ｎ体動基線）と重なっていないという観点からの判断を併用することによって、判断の精度が向上する。
【０１９４】
図９（Ａ）〜図９（Ｄ）は、外乱ノイズの混入が著しい例（脈波信号に含まれるノイズ量の程度が大きい、すなわち脈波信号のきれいさの程度がノイジーである例）を示す図である。例えば、拍動検出装置を、周囲の物体に衝突させたことによって衝撃ノイズが混入したような場合が該当する。
【０１９５】
図９（Ａ）の上側には、前回の１６秒分のフィルター後信号が示されており、下側には前回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図９（Ｂ）の上側には、前回の１６秒分の脈波信号が示されており、下側には前回の脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。また、図９（Ｃ）の上側には、今回の１６秒分のフィルター後信号が示されており、下側には今回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図９（Ｄ）の上側には、今回の１６秒分の脈波信号が示されており、下側には今回の脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【０１９６】
図９（Ａ）〜図９（Ｄ）に示されるいずれの例でも、目立つ外乱ノイズばかりであり、拍動基線Ｓは見あたらない。拍動基線Ｓが仮に存在したとしても、その拍動基線Ｓの信号レベルは、他の基線の信号レベルよりも相対的に低い。よって、この場合には、判断部８０は、脈波信号がノイジーな状態であると判断する。この結果、拍動呈示スペクトルの捕捉処理は実行されない。
【０１９７】
（脈拍計測の結果の一例）
図１０は、腕時計型の脈拍計測装置を被検体の手首に装着し、被検体が歩行（ウォーキング）を行ったときにおける、検出された脈拍数の変化を示す図である。図１０の例では、計測開始（ｔ＝０）から６０秒間の間、被検体（ユーザー）は安静な状態を保ち、６０秒が経過した時点から歩行を開始している。
【０１９８】
計測開始後６０秒〜１２０秒までの期間Ｔ１では、被検体の運動開始に伴って脈拍数が、かなり急激に上昇している。この期間Ｔ１において、脈拍検出に失敗する可能性が高いと言える。すなわち、第１処理部６０による、フィルター後信号に基づく周波数解析処理では、脈拍数の変化が大きい期間Ｔ１において、周波数スペクトル値が小さくなり、かつ個々の周波数スペクトルが分散する結果となり、脈拍数の変化が少ない期間に比べて、特徴的な拍動基線（拍動呈示スペクトル）の検出の難易度が高くなる。
【０１９９】
しかし、本実施形態の拍動検出装置によれば、拍動呈示スペクトルの捕捉処理が実行されることから、期間Ｔ１においても、正しい脈拍数を安定的に計測することが可能である。
【０２００】
（第３実施形態）
図１１（Ａ）および図１１（Ｂ）は、拍動検出装置の、被検体への装着例を示す図である。
【０２０１】
図１１（Ａ）の例は、腕時計型の拍動検出装置の例である。脈波センサー１０および表示部９４を含むベース部４００は、リストバンド３００によって、被検体（ユーザー）の左手首２００に装着されている。
【０２０２】
図１１（Ｂ）の例は、指装着型の拍動検出装置の例である。被検体の指先に挿入するためのリング状のガイド３０２の底部に、脈波センサー１０が設けられている。
【０２０３】
本発明はいずれのタイプの拍動検出装置にも適用可能であるが、腕時計型の拍動検出装置（図１１（Ａ）の例）に適用するのが、より好ましい。図１１（Ａ）の例では、脈波センサー１０が、手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着される。このため、脈波センサー１０から出力される脈波信号の振幅が総じて小さくなる傾向にある。本発明を適用することによって、腕時計型の拍動検出装置における拍動呈示スペクトルの特定性能が向上する。
【０２０４】
このように、本発明の少なくとも一つの実施形態によれば、例えば、拍動検出装置の検出性能を高めることができる。また、例えば、脈波センサーを手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着するタイプの拍動検出装置において、拍動信号を示す拍動呈示スペクトルの特定性能を向上させることができる。
【０２０５】
なお、上記のように本実施形態について詳細に説明したが、本発明の新規事項及び効果から実体的に逸脱しない多くの変形が可能であることは当業者には容易に理解できるであろう。従って、このような変形例はすべて本発明の範囲に含まれるものとする。また、明細書又は図面において、少なくとも一度、より広義又は同義な異なる用語と共に記載された用語は、明細書又は図面のいかなる箇所においても、その異なる用語に置き換えることができる。
【符号の説明】
【０２０６】
１０ 脈波センサー、１１ 体動センサー（加速度センサーやジャイロセンサー等）、
１２ 脈波信号蓄積部、１７ 適応フィルター切り替え部、３０ フィルター部、
４０（４０ａ〜４０ｃ） 周波数分解部、
６０ 第１処理部（後処理部）、６２ ピーク順ソート部、
６４ 相関判定部、６６ 拍動／体動分離部（拍動／ノイズ分離部）、
７０ 第２処理部（捕捉処理部）、７２ 第１捕捉処理部、７４ 第２捕捉処理部、
７６ 捕捉部、８０ 判断部、９０ 脈拍数算出部、９２ 表示処理部、
９４ 表示部、１００ 拍動検出装置
【技術分野】
【０００１】
本発明は、拍動検出装置等に関する。
【背景技術】
【０００２】
拍動検出装置は、人体の心拍に由来する拍動を検出するための装置であって、例えば、腕、手のひら、手指などに装着される脈波センサーからの信号（脈波信号）から、人体の体動の影響により発生する信号成分（体動影響信号）を雑音として除去し、心拍に由来する信号（拍動信号）を検出する装置である。
【０００３】
人の指や手首に装着するタイプの脈波センサーは、例えば、特許文献１〜特許文献３に記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００５−１９８８２９号公報
【特許文献２】特開２００７−５４４７１号公報
【特許文献３】特開２００５−１３１４２６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
拍動信号成分を明確に峻別するに際し、脈波センサーから出力される脈波信号に含まれる拍動信号成分の信号レベルが著しく低い場合には、信号レベルが十分な場合に比べる、峻別処理の難易度が高くなる。
【０００６】
本発明の少なくとも一つの態様によれば、例えば、拍動検出装置の検出性能を高めることができる。
【課題を解決するための手段】
【０００７】
（１）本発明の拍動検出装置の一態様は、被検体の拍動に由来する拍動信号を検出する拍動検出装置であって、前記拍動信号と、前記被検体の体動に由来する体動信号を含むノイズ信号とが混在した脈波信号を検出して出力する脈波センサーと、前記脈波信号から、前記ノイズ信号の少なくとも一部を除去するフィルター部と、前記フィルター部から出力されるフィルター後信号、または前記脈波センサーから出力される前記脈波信号に基づいて、所定時間毎に前記周波数解析処理を行って、前記拍動信号を示す拍動呈示スペクトルを特定する周波数解析部と、を含み、前記周波数解析部は、前記フィルター後信号の周波数スペクトルに基づいて、前記拍動信号と前記ノイズ信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行して前記拍動呈示スペクトルを特定する第１処理部と、前記脈波信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部と、前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部と、を含む。
【０００８】
本態様では、周波数解析部は、第１処理部、第２処理部ならびに判断部を有する。第１処理部は、脈波センサーによって検出された脈波信号をフィルター部によってフィルタリングして得られるフィルター後信号に基づいて、拍動信号を示す拍動呈示スペクトルを特定する。第１処理部は、拍動信号と体動信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行する。
【０００９】
本態様では、例えば、脈波センサーから出力される脈波信号に含まれる拍動信号成分の信号レベルが著しく低い場合においても、第２処理部によって拍動呈示スペクトルの捕捉を試みることができ、これによって拍動検出の可能性が向上する。この第２処理部は、フィルター前の脈波信号に基づいて拍動呈示スペクトルの捕捉を試みる。
【００１０】
フィルタリング処理を行った場合、残留した体動信号による体動スペクトルと、拍動呈示スペクトルとの大きさが同程度になることがあり、この場合、単純にスペクトルの大きさを比較するだけでは、拍動呈示スペクトルがどれなのかを判定することが難しい。この場合であっても、フィルター前の脈波信号に基づいた拍動呈示スペクトルの捕捉処理を試みることによって、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低い拍動信号成分を検出することができる場合がある。
【００１１】
但し、フィルター前の脈波信号は、ノイズ対策がまったく行われていない原信号であり、例えば、過大なノイズしか存在しない信号である場合もあり得る。そこで、本態様では、第２処理部による拍動呈示スペクトルの捕捉処理を試みる前に、判断部が、拍動呈示スペクトル捕捉処理の実行の可否を判断することとした。したがって、本態様によれば、拍動検出装置の検出性能を高めることができる。
【００１２】
（２）本発明の拍動検出装置の他の態様では、前記判断部は、前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、を有し、前記評価部による評価に基づいて、前記第２処理部による前記拍動呈示スペクトル捕捉処理の実行の可否を判断する。
【００１３】
本態様では、判断部は、脈波信号のきれいさに関する判断指標を取得する判断指標取得部と、判断指標を用いて、脈波信号に含まれるノイズ量の程度を評価する評価部とを有する。判断部は、脈波信号のノイズ状態を判断するための判断指標を取得し、その判断指標を用いて、脈波信号に含まれるノイズ量の程度を評価する。そして、その評価に基づいて、第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する。
【００１４】
判断部は、例えば、評価部によって、脈波信号に混入されるノイズ量（例えば外乱ノイズ成分のノイズ量）が所定の基準以下であると判断され場合に、第２処理部による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。換言すれば、脈波信号が、許容レベル以上にきれいであることを第２処理部による拍動呈示スペクトル捕捉処理の実行の条件としてもよい。
【００１５】
第２処理部による救済を試みると、ノイズを拍動信号成分と誤検出する可能性が高い。本態様によれば、例えば、外乱ノイズが多い状態では拍動呈示スペクトルの捕捉処理が実行されない。本態様によれば、脈波信号のノイズ状態を慎重に判断しつつ、第２処理部による拍動呈示スペクトルの捕捉処理の可否を判断することができる。
【００１６】
（３）本発明の拍動検出装置の他の態様では、前記判断部は、前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、を有し、前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、前記評価部による評価が許容範囲にあることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断する。
【００１７】
本態様では、周波数解析部による前回の周波数解析によって、拍動呈示スペクトルが検出されていること、ならびに、評価部による、脈波信号に含まれるノイズ量の程度の評価が許容範囲にあることを条件として、第２処理部による拍動呈示スペクトルの捕捉処理を実行する。
【００１８】
拍動呈示スペクトルは、所定時間前の過去の拍動信号成分と現在の拍動信号成分との相関性に基づいて捕捉することができる。よって、拍動呈示スペクトルの捕捉処理を実行する条件として、前回の周波数解析処理、すなわち、所定時間前の周波数解析処理において拍動信号の検出に成功していることを条件としてもよい。本態様では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【００１９】
（４）本発明の拍動検出装置の他の態様では、前記判断部は、前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、前記被検体の運動状態を判断する運動状態判断部と、を有し、前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、前記評価部による評価が許容範囲にあり、かつ、前記運動状態判断部によって前記被検体が定常状態にあると判断されることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断する。
【００２０】
被検体の運動状態が急激に変化している場合には、脈拍等の変化が急である。
【００２１】
よって、本態様では、よって、拍動呈示スペクトルの捕捉処理を実行する条件として、運動状態判断部によって被検体が定常状態にあると判断されることを条件としてもよい。被検体が定常状態である場合とは、例えば、被検体が安定した運動状態や静養状態にある場合をいう。
【００２２】
被検体の運動状態は、例えば、体動センサーから得られる体動信号のスペクトルが、所定の時間にわたって大きく変化しないという観点から判断することができる（但し、これに限定されるものではない）。本態様では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【００２３】
（５）本発明の拍動検出装置の他の態様では、前記判断部は、前記第１処理部が拍動呈示スペクトルを特定することができなかった場合に、前記判断部は、前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する。
【００２４】
本態様では、判断部による拍動呈示スペクトル捕捉処理の実行の可否判断は、第１処理部が、フィルター後信号に基づいて拍動呈示スペクトルを特定することができなかった場合に行われる。
【００２５】
第１処理部と第２処理部とを並列に動作させることも可能であるが、第１処理部によって拍動呈示スペクトルが特定できた場合には、第２処理部による拍動呈示スペクトルの捕捉は不要である。そこで、第２処理部の無駄な動作をなくすために、本態様では、第１処理部による拍動呈示スペクトルの特定の失敗を条件として、判断部による拍動呈示スペクトルの捕捉処理の可否判断を開始させることとした。
【００２６】
（６）本発明の拍動検出装置の他の態様では、前記被検体の体動に伴う体動信号を検出して出力する体動センサーと、前記体動センサーから出力される体動信号に基づいて周波数分解処理を行う周波数分解部と、を有し、前記周波数分解部による周波数分解結果に基づいて推定される、前記被検体の体動由来の体動成分を体動基線とし、異なる周波数のｍ本（ｍは１以上の整数）の前記体動基線のうちの第ｎ番目（１≦ｎ≦ｍ）の体動基線を第ｎ体動基線としたとき、前記第２処理部は、周波数解析処理に使用される前記脈波信号の周波数スペクトルにおける前記被検体の拍動由来の拍動成分と推定される第１拍動基線と、前記所定時間前の周波数解析処理に使用された前記脈波信号または前記フィルター後信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルとして特定された第２拍動基線との相関性を重視した第１捕捉処理を実行して、前記拍動呈示スペクトルの捕捉を試みる第１捕捉処理部と、前記脈波信号に、前記第１処理部による前記拍動／ノイズ分離処理と同じ処理を施した後、前記第１拍動基線と、前記第ｎ体動基線との区別の可否に基づいて第２捕捉処理を実行して前記拍動呈示スペクトルの捕捉を試みる第２捕捉処理部と、を有し、前記判断部は、前記拍動呈示スペクトル捕捉処理を実行可と判断した場合、さらに、前記評価部による評価が第１の許容レベルであるか、前記第１許容レベルよりも低い第２許容レベルであるかを判断し、前記判断部によって前記第１許容レベルであると判断されたときは、前記第１捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行し、前記第２許容レベルであると判断されたときは、前記第２捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行する。
【００２７】
本態様では、第２処理部の構成と動作の一例が示される。第２処理部は、第１捕捉処理部および第２捕捉処理部を有する。判断部は、拍動呈示スペクトル捕捉処理を実行可と判断した場合、さらに、評価部による評価（ノイズ量の程度に基づく評価）が第１の許容レベルであるか、第１許容レベルよりも低い第２許容レベルであるかを判断する。
【００２８】
第１捕捉処理部は、判断結果が第１許容レベルのとき（例えば、ノイズ量の程度が小さく、脈波信号がきれいであるとき）に動作する。第２捕捉処理部は、判断結果が、第１許容レベルよりも低い第２許容レベルのとき（例えば、ノイズ量の程度が第１許容レベルを基準とする場合よりも大きく、脈波信号のきれいさが中程度のとき）に動作する。
【００２９】
そこで、第１捕捉処理部は、第１捕捉処理を実行する。この第１捕捉処理では、第１拍動基線と第２拍動基線との相関性に基づく拍動呈示スペクトルの捕捉処理を実行する。第１拍動基線は、今回の周波数解析処理に使用される、「脈波信号」の周波数スペクトルにおける被検体の拍動由来の拍動成分と推定されるスペクトルである。第２拍動基線とは、前回の周波数解析処理に使用された、「脈波信号またはフィルター後信号」の周波数スペクトルに基づいて拍動呈示スペクトルとして特定された拍動基線である。
【００３０】
第１捕捉処理部は、第１拍動基線と第２拍動基線との相関を重視して判断したとき、相関性が有る場合には、第１拍動基線を拍動呈示スペクトルと判断する。
【００３１】
また、第２捕捉処理部は、例えば、脈波信号が中程度である場合に動作することができる。この場合には、拍動呈示スペクトルの捕捉のためには、拍動基線と外乱ノイズとの分離が必要となる。
【００３２】
第２捕捉処理部は、第１処理部で行われる拍動／ノイズ分離処理と同じ処理を実行してもよい。また、第２捕捉処理部は、さらに、今回の検出時の周波数スペクトルにおける第１拍動基線と、第ｎ体動基線との区別の可否に基づいて拍動呈示スペクトルの捕捉を試みてもよい。
【００３３】
ここで、体動基線（広義）は、体動センサーから出力される体動信号の、現在ならびに過去の周波数分解結果に基づいて、体動成分と推定されるスペクトルの総称である。外乱ノイズがある程度存在する場合、例えば、被検体が運動している場合等では、体動基線には、異なる周波数のｍ本（ｍは１以上の整数）の体動基線が含まれることが多い。
【００３４】
第ｎ体動基線は、異なる周波数のｍ本（ｎ、ｍは１以上の整数）の体動基線のうち、一番周波数の低い基線を第１番目としそのｎ倍の周波数に相当する基線を第ｎ体動基線と定義する。例えば、第１波（基本波）の拍動基線が、第ｎ体動基線（１≦ｎ≦ｍ）の位置とは異なる位置にある場合には、第２捕捉処理部は、第１拍動基線が拍動呈示スペクトルであると判断する。この場合、例えば、ｎが２〜４である場合の体動基線は、ｎが１である場合の体動基線よりも信号値が小さく、外乱ノイズ中に潜んでいる場合があるため、第ｎ体動基線（例えばｎ＝１〜４）との区別を試みることは、拍動呈示スペクトルの捕捉精度を高める点で有効である。
【００３５】
なお、第２捕捉処理においては、第１拍動基線と第２拍動基線との相関性を考慮してもよい。例えば、今回検出時（現在の検出期間における検出時）の周波数スペクトルにおける第１拍動基線と、前回検出時（所定時間前の検出持）の周波数スペクトルにおける第２拍動基線との間の周波数軸上における距離の観点からの判断を行ってもよい。
【００３６】
本態様では、脈波信号に含まれるノイズ量の程度に応じて、拍動呈示スペクトルの捕捉処理の内容を切り替えることから、脈波信号の状態に応じた適切な拍動呈示スペクトルの捕捉処理を実行することができる。
【００３７】
（７）本発明の拍動検出装置の他の態様では、前記第１捕捉処理部は、前記第１拍動基線の周波数軸上における第１位置と、前記第２拍動基線の周波数軸上における第２位置との間の距離によって、前記相関性を判断する。
【００３８】
本態様では、第１捕捉処理部の動作の一例が示される。今回の検出（現在の検出期間における検出）に係る第１拍動基線の周波数軸上の第１位置と、前回の検出（所定時間前の検出）に係る第２拍動基線の周波数軸上における第２位置との間の距離を、相関性の判断の基準としてもよい。例えば、周波数軸上における距離が、例えば、距離２以内であるとき、つまり、周波数軸上において、第２拍動基線の第２位置を基準として、第１拍動基線が２本のスペクトル以内で隣り合っているとき、第１捕捉処理部は、相関ありと判断して、第１拍動基線を拍動呈示スペクトルと判断してもよい。
【００３９】
（８）本発明の拍動検出装置の他の態様では、前記第２捕捉処理部は、前記第１拍動基線が前記第ｎ体動基線と重なっていないという第１条件と、前記第１拍動基線が前記第２拍動基線に対して、周波数軸上で所定の距離以内であるという第２条件と、を満足する場合に、前記第１拍動基線を、拍動呈示スペクトルとして特定する。
【００４０】
本態様では、第２捕捉処理部は、第１条件の他に、第２条件が満足されるか否かの判断も実行し、より慎重に拍動呈示スペクトルを捕捉する。
【００４１】
例えば、第１拍動基線が、体動センサー出力の周波数分解結果に基づいて特定される第ｎ体動基線に重なっていないことを第１条件とし、第１拍動基線が、周波数軸上で、第２拍動基線に対して所定距離内にあることを、第２条件としてもよい。
【００４２】
「所定距離内」は、例えば、距離２以内とすることができる。距離２とは、第１拍動基線が、第２拍動基線に対して、１本のスペクトルを隔てて隣り合っていることである。距離２以内とは、距離０、距離１または距離２のいずれかであることを意味する。すなわち、本態様では、第２捕捉処理部は、周波数軸上でスペクトル（基線）が３本分離れるということが発生しないものとして、拍動呈示スペクトルの呈示処理を行う。上述の第１条件と第２条件を満足する場合には、第２捕捉処理部は、第１拍動基線を拍動呈示スペクトルとして捕捉する。これによって、例えば、目立つノイズに引きずられることなく、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような拍動呈示スペクトルを特定できる可能性が高まる。
【００４３】
このように、本発明の少なくとも一つの態様によれば、例えば、拍動検出装置の検出性能を高めることができる。
【図面の簡単な説明】
【００４４】
【図１】本発明の拍動検出装置の一例の構成を示す図
【図２】図２（Ａ）および図２（Ｂ）は、第１捕捉処理部による第１捕捉処理の一例を示す図
【図３】図３（Ａ）および図３（Ｂ）は、第２捕捉処理部による第２捕捉処理の一例を示す図
【図４】本発明の拍動検出装置の他の例の構成を示す図
【図５】図５（Ａ）〜図５（Ｃ）は、脈波信号に含まれるノイズ量の程度の判断動作を説明するための脈波信号の波形および周波数スペクトルの一例を示す図
【図６】第２処理部による、拍動呈示スペクトルの捕捉処理（第２処理）の手順の一例を示すフローチャート
【図７】図７（Ａ）および図７（Ｂ）は、第２捕捉処理の適用例を示す図
【図８】図８（Ａ）〜図８（Ｃ）は、第１捕捉処理の適用例を示す図
【図９】図９（Ａ）〜図９（Ｄ）は、外乱ノイズの混入が著しい例を示す図
【図１０】腕時計型の脈拍計測装置を被検体の手首に装着し、被検体が歩行（ウォーキング）を行ったときにおける、検出された脈拍数の変化を示す図
【図１１】図１１（Ａ）および図１１（Ｂ）は、拍動検出装置の、被検体への装着例を示す図
【図１２】図１２（Ａ）および図１２（Ｂ）は、拍動信号成分が強い場合におけるフィルタリングの効果を示す図
【図１３】図１３（Ａ）および図１３（Ｂ）は、拍動信号成分が弱い場合におけるフィルタリングの効果を示す図、
【図１４】図１４（Ａ）および図１４（Ｂ）は、被検体が安定運動している最中に、非周期的な外乱ノイズが混入した場合の例を示す図
【発明を実施するための形態】
【００４５】
本実施形態の説明に入る前に、フィルター後信号に基づいて拍動呈示スペクトルを特定するのが困難である場合について考察する。フィルター後信号とは、脈波センサーから出力される脈波信号に、フィルター部によるフィルタリング処理を施すことによって得られる信号である。フィルタリング処理は、脈波信号に含まれるノイズを除去あるいは抑制するための信号処理である。
【００４６】
図１２（Ａ）および図１２（Ｂ）は、拍動信号成分が強い場合におけるフィルタリングの効果を示す図である。図１２（Ａ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図１２（Ｂ）の１６秒分フィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【００４７】
図１２（Ａ）および図１２（Ｂ）の例は、拍動信号成分が強い場合の例である。図１２（Ａ）では、拍動信号のスペクトルＳ２０の信号値は、比較的大きなノイズ成分のスペクトルＮ２１〜Ｎ２３と比べて十分に大きい。したがって、脈波信号とノイズとを区別するためのフィルタリング処理を実行すると、図１２（Ｂ）に示すように、ノイズは十分に低減されてＳ／Ｎが向上し、拍動信号のスペクトルＳ３０をノイズと区別することが容易となる。
【００４８】
ここで、図１３（Ａ）および図１３（Ｂ）を参照する。図１３（Ａ）および図１３（Ｂ）は、拍動信号成分が弱い場合におけるフィルタリングの効果を示す図である。図１３（Ａ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図１３（Ｂ）の上側には１６秒分のフィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【００４９】
図１３（Ａ）では、ノイズのスペクトルＮ４０，Ｎ４１の信号値が大きく、これに対して、拍動呈示スペクトルＳ４０の信号値は小さく、拍動呈示スペクトルＳ４０はノイズに拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低い状態となっている。図１３（Ｂ）の例では、フィルタリングによるノイズ除去の効果が顕著に表れ、フィルタリング後のＳ／Ｎは、フィルター前のＳ／Ｎよりも格段に向上している。
【００５０】
しかしながら、フィルタリング処理を行った場合、残留した体動信号による体動スペクトルと、拍動呈示スペクトルとの大きさが同程度になることがあり、この場合、単純にスペクトルの大きさを比較するだけでは、拍動呈示スペクトルがどれなのかを判定することが難しい。
【００５１】
図１４（Ａ）および図１４（Ｂ）は、被検体が安定運動している最中に、非周期的な外乱ノイズ（例えば脈波センサーを物体に衝突させた場合の衝撃ノイズ）が混入した場合の例を示す図である。図１４（Ａ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図１４（Ｂ）の上側には１６秒分のフィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【００５２】
図１４（Ａ）では、拍動呈示スペクトルＳ６０の周波数よりも低い周波数帯域において、大きな外乱ノイズのスペクトルＮ６０が存在している。図１４（Ｂ）では、フィルタリングの効果によって、外乱ノイズのスペクトルＮ６０は抑圧されるが、拍動呈示スペクトルＳ６０の信号値も低下しており、結果的に、拍動呈示スペクトルＳ６０の特定が困難である。
【００５３】
以下に説明する本実施形態の拍動検出装置では、このような状況であっても、フィルター前の脈波信号に基づく周波数解析を実行して、拍動呈示スペクトルの捕捉を試みる。ここで、図１３（Ａ）および図１３（Ｂ）に示す例に着目する。
【００５４】
図１３（Ｂ）の例では、拍動呈示スペクトルＳ５０が、ノイズのスペクトルＮ５０と同程度の大きさである。また、拍動呈示スペクトルＳ５０と近い周波数帯に外乱ノイズのスペクトルも多数存在しており、よって、図１２の場合と比較すると拍動呈示スペクトル特定処理の難易度が高い。
【００５５】
これに対し、図１３（Ａ）の例では、Ｎ４０、Ｎ４１といった体動呈示スペクトルが存在するものの、これらとは離れた周波数帯に、拍動呈示スペクトルＳ４０が存在する。また、更には拍動呈示スペクトルＳ４０は、隣り合う他のスペクトルと比較すると相対的には大きいスペクトルとなっている。よって、脈波信号ｄが図１３（Ａ）のような信号状態であることを検出した上であれば、Ｓ４０が確かに拍動呈示スペクトルであると容易に判定可能と言える。
【００５６】
つまり、例えば、所定時間前の過去と現在の拍動信号成分の相関性等に着目すれば、拍動呈示スペクトルＳ４０が存在する周波数域から離れた周波数域における大きなノイズのスペクトルに惑わされることなく、拍動呈示スペクトルを特定する可能性があるといえる。
【００５７】
そこで、以下の本実施形態では、拍動検出装置に、フィルター前の脈波信号に基づく周波数解析の機能を付与し、新規な検出アルゴリズムを用いて拍動呈示スペクトルの捕捉を行えるようにした。
【００５８】
但し、フィルター前の脈波信号は、ノイズフィルタリング処理を行っていない原信号であり、例えば、拍動信号成分と比較して大きなノイズが含まれる場合があり得る。この場合に該当するのが、図１４に示す例である。この場合には、拍動呈示スペクトルの捕捉処理を行わないのが好ましい。
【００５９】
このような観点から、以下に述べる実施形態においては、判断部を設けて、拍動呈示スペクトルの捕捉処理を試みる前に、判断部が、拍動呈示スペクトル捕捉処理の実行の可否を判断することとした。
【００６０】
判断部は、例えば、脈波信号に含まれるノイズが多い状態である場合や、被検体（人ならびに動物を含む）の運動状態が急激に変化している場合、所定時間前の過去の検出時点において拍動呈示スペクトルの検出に失敗している場合等においては、拍動呈示スペクトルの捕捉処理は不可と判断してもよい（但し、上述の場合は一例であり、これらの場合に限定されるものではない）。
【００６１】
以下、図１〜図３を参照して、具体的に説明する。なお、以下に説明する本実施形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また、本実施形態で説明される構成の全てが、本発明の必須構成要件であるとは限らない。
【００６２】
（第１実施形態）
図１は、本発明の拍動検出装置の一例の構成を示す図である。図１に示される拍動検出装置１００は、被検体（人や動物を含む）の拍動に由来する拍動信号、拍動信号に対応する心拍等の生体情報等を検出するセンサー装置の一種である。
【００６３】
ここで、拍動とは、医学的には心臓のみならず内臓一般の周期的な収縮、弛緩が繰り返された場合に起こる運動のことをいう。ここでは、心臓が周期的に血液を送るポンプとしての動きを拍動と呼ぶ。なお、心拍数とは、１分間の心臓の拍動の数をいう。また、脈拍数は、末梢血管における脈動の数をいう。心臓が血液を送り出す際に、動脈に脈動が生じるので、この回数を数えたものを脈拍数あるいは単に脈拍と呼ぶ。腕で脈を計測する限りは、医学的には心拍数とは呼ばずに脈拍数と呼ぶのが通常である。また、以下の説明では、体動という用語が使用される。体動とは、広い意味では、体を動かすことすべてを意味する（広義の体動）。体動センサーから出力される体動信号は、広義の体動に対応する信号であり、体動センサー信号と言い換えることができる。但し、体動信号の中には、被検体の周期的な体動（例えば、歩行・ジョギングなどに伴う定常的、周期的な腕（脈拍計の装着部位近辺）の振り等）に伴う周期的な体動信号が含まれる。この周期性をもつ体動信号は、狭義の体動（周期性のある体動）に対応する信号である。なお、例えば、広義の体動信号を「体動センサー信号」という場合、狭義の体動信号を、単に「体動信号」といってもよい。
【００６４】
（拍動検出装置の構成例）
図１に示される拍動検出装置１００は、脈波センサー１０と、体動センサー（加速度センサー等）２０と、フィルター部（前処理部）３０と、周波数解析部５０と、脈拍数算出部９０と、表示処理部９２と、表示部９４と、を有する。
【００６５】
脈波センサー１０は、拍動信号と、被検体の体動に由来する体動ノイズ信号を含むノイズ信号とが混在した脈波信号ｄを検出して出力する。脈波センサー１０は、例えば、光電脈波センサー、およびその原理に基づく脈波センサーとすることができる。
【００６６】
体動センサー２０は、被検体の体動、例えば、歩行やジョギングなどに伴う定常的、周期的な腕（脈拍計の装着部位近辺）の動きを検出するセンサーであり、例えば、加速度センサーやジャイロセンサーを含むことができる。
【００６７】
フィルター部３０は、脈波信号ｄおよび体動信号（体動センサー２０から出力される信号）ｆを入力として、フィルタリング処理を実行して、脈波信号ｄに含まれるノイズを最小化する。フィルター部３０は、特に限定されるものではないが、適応フィルター（図１では不図示、図５の参照符号３２）を含んでもよい。
【００６８】
適応フィルターは、入力信号に基づいて、伝達関数を自己適応的に変化させることができるデジタルフィルターである。適応フィルターは一般にデジタル信号処理を行うデジタルフィルターとして実装される。フィルターの所望の性能（例えば、入力信号に含まれるノイズ成分を最小化するといった性能）を維持するために、所定のアルゴリズム（例えば最適化アルゴリズム）が使用される。所定のアルゴリズムには、フィルタリング後の信号に基づいて得られる信号がフィードバックされ、そのアルゴリズムによる適応処理によってフィルター係数が適応的に変化し、その結果として、適応フィルターの周波数応答特性が変化する。適応フィルターを使用することによって、例えば、周期性のある信号と非周期的なノイズとを分離することができる。
【００６９】
周波数解析部５０は、フィルター部３０から出力されるフィルター後信号ｅ、または脈波センサー１０から出力される脈波信号ｄに基づいて、所定時間毎（例えば４秒毎）に周波数解析処理を行って、拍動信号を示す拍動呈示スペクトルを特定する。ここで、フィルター後信号ｅは、脈波信号ｄをフィルタリングして得られる信号である。
【００７０】
周波数解析部５０は、高速フーリエ変換（ＦＦＴ）によって実時間上の信号を周波数軸上の信号に変換する周波数分解部４０（第１周波数分解部４０ａ、第２周波数分解部４０ｂならびに第３周波数分解部４０ｃを含む）と、フィルター後信号ｅのＦＦＴ結果ｇ１ならびに体動信号ｆのＦＦＴ結果ｇ２に基づいて、拍動呈示スペクトルを特定する第１処理部（後処理部）６０と、脈波信号ｄのＦＦＴ結果ｇ３に基づいて、拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部（捕捉処理部）７０と、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部８０と、を含む。
【００７１】
なお、第１周波数分解部４０ａは、フィルター後信号ｅにＦＦＴ（高速フーリエ変換）を施す。第２周波数分解部４０ｂは、体動信号ｆにＦＦＴを施す。また、第３周波数分解部４０ｃは、脈波信号ｄにＦＦＴを施す。
【００７２】
第１処理部（後処理部）６０は、ピーク順ソート部６２と、相関判断部６４と、拍動／ノイズ分離部６６と、拍動呈示スペクトル特定部６８と、を含むことができる。
【００７３】
第２処理部（捕捉処理部）７０は、相関性に基づく捕捉処理を実行する第１捕捉処理部７２と、第２捕捉処理部７４（図１において破線で囲んで示すブロック）とを含む。第２捕捉処理部７４は、第１処理部６０の構成要素であるピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６に、距離と区別性に基づく捕捉部７６を追加して構成される。ピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６と同等の機能をもつブロックを新たに設けてもよい。但し、この場合、同じ機能をもつブロックが重複する。よって本実施形態では、第１処理部６０の構成要素であるピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６を、第２捕捉処理部７４の構成要素としても活用している。
【００７４】
第１処理部６０による第１処理（後処理）は、脈波信号ｄに含まれる拍動信号成分の信号値（振幅値）が大きい場合に有効である。第２処理部７０による第２処理（拍動呈示スペクトルの捕捉処理）は、例えば、第１処理部６０による拍動呈示スペクトルの特定に失敗した場合に行うことができる（但し、これに限定されるものではない）。
【００７５】
判断部８０は、周波数解析部５０に含まれる第３周波数分解部４０ｃが生成した周波数スペクトルｇ３に基づいて、第２処理部７０による拍動呈示スペクトルの捕捉処理の可否を判断する。この判断は、例えば、脈波信号ｄに含まれるノイズ量の程度の評価に基づいて行うことができる。判断部８０による判断動作については、第２実施形態にて詳述する。
【００７６】
例えば、判断部８０が、脈波信号ｄに含まれるノイズ量の程度が大きい（脈波信号ｄのきれいさの程度がきれい：目立つノイズ少）と判断したとき、第２処理部７０による拍動呈示スペクトルの捕捉処理が実行される。この場合、例えば、第１捕捉処理部７２が捕捉処理を実行する。また、脈波信号ｄに含まれるノイズ量の程度が中程度（脈波信号ｄのきれいさの程度がまあまあ：外乱中程度）の場合には、第２捕捉処理部７４が捕捉処理を実行する。第１捕捉処理ならびに第２捕捉処理の内容の一例については、図２および図３を用いて後述する。
【００７７】
第１処理部６０または第２処理部７０が拍動呈示スペクトルの特定、捕捉に成功すると、脈拍数算出部９０は、被検体の脈拍数を算出する。脈拍数は、例えば、拍動呈示スペクトルの周波数に基づいて算出する。求められた脈拍数の情報は、脈拍数算出部９０から表示処理部９２を経由して表示部９４に供給され、例えば、脈波数を示す数値が、表示部９４によって表示される。なお、脈拍数ではなく、検出した脈の、時間軸上における変化を信号波形やグラフの形式で表示してもよい。
【００７８】
（図１の構成によって得られる効果等）
第１処理部６０は、フィルター後信号ｅのＦＦＴ結果ｇ１および体動信号ｆのＦＦＴ結果ｇ２に基づいて、拍動信号を示す拍動呈示スペクトルを特定する。図１の拍動検出装置は、第２処理部７０を備える。第２処理部７０は、脈波信号ｄに外乱ノイズが多く含まれ、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような場合であっても、第２処理部７０によって拍動呈示スペクトルの捕捉を試みる。これによって拍動検出の可能性が向上する。この第２処理部７０は、フィルター前の脈波信号ｄに基づいて拍動呈示スペクトルの捕捉を試みる。フィルター前の脈波信号ｄに基づいた拍動呈示スペクトルの捕捉処理を試みることによって、拍動信号成分を検出することができる可能性（拍動呈示スペクトルの捕捉処理による救済の効果）が生じる。
【００７９】
但し、フィルター前の脈波信号ｄは、ノイズ対策がまったく行われていない原信号であり、例えば、過大なノイズしか存在しない信号である場合もあり得る。
【００８０】
そこで、図１の拍動検出装置１００では、第２処理部７０による拍動呈示スペクトルの捕捉処理を試みる前に、判断部８０が、拍動呈示スペクトル捕捉処理の実行の可否を判断することとした。判断部８０は、例えば、脈波信号ｄに含まれるノイズが多い状態である場合や、被検体（人ならびに動物を含む）の運動状態が急激に変化している場合、所定時間前の過去の検出時点において拍動呈示スペクトルの検出に失敗している場合等においては、第２処理部７０による拍動呈示スペクトルの捕捉処理は不可と判断する。具体例については後述する。
【００８１】
したがって、図１の拍動検出装置１００によれば、拍動検出性能を高めることができる。また、例えば、脈波センサーを手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着するタイプの拍動検出装置において、拍動信号を表す拍動呈示スペクトルの特定性能を向上させることができる。
【００８２】
（第１処理部による第１処理の内容）
第１処理部６０による第１処理は、フィルター後信号ｅのＦＦＴ結果ｇ１ならびに体動信号ｆのＦＦＴ結果ｇ２に基づいて実行することができる。第１処理は、脈波信号ｄに含まれる拍動信号成分の信号値（振幅値）が大きい場合に有効である。
【００８３】
第１処理は、例えば、以下の手順により行われる。ピーク順ソート部６２は、フィルター後信号ｅ１のＦＦＴ結果ｇ１に基づいて、スペクトルを信号値が大きい順にソーティングする。例えば、上位１０本のスペクトルのいずれかが拍動信号を示す拍動呈示スペクトルである可能性が高い。
【００８４】
次に、相関判断部６４が、相関判断を実行する。周波数解析部５０による周波数解析が４秒毎に行われる場合、相関判断部６４は、所定時間前の過去（つまり直近の過去）である４秒前の検出時に検出された拍動スペクトルを基準として、今回の検出（現在の検出期間における検出）に係る上位１０本のスペクトルについて相関の有無を判断する。これにより、拍動呈示スペクトルと推定される１本の候補が抽出される。
【００８５】
次に、拍動／ノイズ分離部６６は、抽出された候補について、体動を示す体動呈示スペクト等のノイズとの分離が可能かを判断する。例えば、抽出された候補の周波数が、体動スペクトルの周波数とは明らかに異なる場合には分離可能と判断される。
【００８６】
また、希に、抽出された候補の周波数と体動スペクトルの周波数とが一致する場合もあるが、この場合でも、例えば、所定時間前の過去および現在の体動信号の相関性等を考慮して、その候補を、体動呈示スペクトルと区別できる場合がある。
【００８７】
拍動呈示スペクトル特定部６８は、体動呈示スペクトルとは異なることが明らかとなった候補を、拍動呈示スペクトルとして特定することができる。
【００８８】
（第２処理部による第２処理の内容）
上述のとおり、第２処理部７０による第２処理（拍動呈示スペクトルの捕捉処理）は、例えば、第１処理部６０による拍動呈示スペクトルの特定ができなかった場合に行ってもよい。すなわち、判断部８０による拍動呈示スペクトル捕捉処理の実行の可否判断は、第１処理部６０が、拍動呈示スペクトルを特定することができなかった場合に行ってもよい。
【００８９】
第１処理部６０と第２処理部７０とを並列に動作させることも可能であるが、第１処理部６０によって拍動呈示スペクトルが特定できた場合には、第２処理部７０による拍動呈示スペクトルの捕捉は不要である。第２処理部７０の無駄な動作をなくすという観点から、第１処理部６０による拍動呈示スペクトルの特定ができないときに、判断部８０による拍動呈示スペクトルの捕捉処理の可否判断を開始させるのが好ましい。
【００９０】
第２処理部７０による第２処理（拍動呈示スペクトルの捕捉処理）は、第１捕捉処理部７２による第１捕捉処理と、第２捕捉処理部７４による第２捕捉処理と、に大別される。上述のとおり、第１捕捉処理部７２は、例えば、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルである場合（例えば、ノイズ量の程度が小さく、脈波信号がきれいであるとき）に動作する。また、第２捕捉処理部７４は、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルよりも低い第２許容レベルである場合（例えば、ノイズ量の程度が第１許容レベルを基準とする場合よりも大きく、脈波信号のきれいさが中程度のとき）に動作する。
【００９１】
脈波信号に含まれるノイズ量の程度に応じて、拍動呈示スペクトルの捕捉処理の内容を切り替えることによって、脈波信号の状態に応じた適切な拍動呈示スペクトルの捕捉処理を実行することができる。
【００９２】
なお、以下の説明では、「拍動基線」ならびに「体動基線」という用語を使用する場合がある。拍動基線は、周波数分解部４０の周波数分解処理によって得られる周波数スペクトルにおける、被検体の拍動由来の拍動成分と推定されるスペクトルである。体動基線は、周波数分解部４０の周波数分解処理によって得られる周波数スペクトルにおける、被検体の体動由来の体動成分と推定されるスペクトルである。
【００９３】
また、「拍動基線（広義）」には、今回の（つまり、現在の検出期間での）周波数解析処理に使用される脈波信号の周波数スペクトルにおける、被検体の拍動由来の成分と推定される「第１拍動基線」と、前回の（つまり所定時間前の）周波数解析処理に使用された脈波信号またはフィルター後信号の周波数スペクトルに基づいて、拍動呈示スペクトルとして特定された「第２拍動基線」とが含まれる。
【００９４】
第１拍動基線は、今回の検出処理に係る拍動基線であり、「現在の拍動基線あるいは今回の拍動基線」と称する場合がある。また、第２拍動基線は、前回（所定時間前の過去）の検出処理に係る拍動基線であり、「所定時間前の過去の拍動基線あるいは前回の拍動基線」と称する場合がある。これらは相関性の有無の判断において使用される。
【００９５】
また、体動基線（広義）は、体動センサー２０から出力される体動信号ｆの、現在ならびに過去の周波数分解結果に基づいて、体動成分と推定されるスペクトルの総称である。外乱ノイズがある程度存在する場合、例えば、被検体が運動している場合等では、体動基線には、異なる周波数のｍ本（ｍは１以上の整数）の体動基線が含まれることが多い。
【００９６】
第ｎ体動基線は、異なる周波数のｍ本（ｎ、ｍは１以上の整数）の体動基線のうち、一番周波数の低い基線を第１番目としそのｎ倍の周波数に相当する基線を第ｎ体動基線と定義する。例えば、第１波（基本波）の拍動基線が、第ｎ体動基線（１≦ｎ≦ｍ）の位置とは異なる位置にある場合には、第２処理部７０は、第１拍動基線が拍動呈示スペクトルであると判断する。この場合、例えば、ｎが２〜４である場合の体動基線は、ｎが１である場合の体動基線よりも信号値が小さく、外乱ノイズ中に潜んでいる場合があるため、第ｎ体動基線（例えばｎ＝１〜４）との区別を試みることは、拍動呈示スペクトルの捕捉精度を高める点で有効である。
【００９７】
（第１捕捉処理）
まず、第１捕捉処理部７２による第１捕捉処理について、図２（Ａ）および図２（Ｂ）を用いて説明する。図２（Ａ）および図２（Ｂ）は、第１捕捉処理部による第１捕捉処理の一例を示す図である。図２（Ａ）および図２（Ｂ）には、０〜４Ｈｚまでの脈波信号をＦＦＴして得られる周波数スペクトルが示されている。
【００９８】
上述のとおり、第１捕捉処理部７２は、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルである場合（例えば、ノイズ量の程度が小さく、脈波信号がきれいであるとき）に動作する。ここで、脈波信号が「きれい」である場合とは、例えば、広範囲の周波数域にわたってスペクトルのレベルが総じて小さく、その中に、信号値の大きな目立つスペクトル(一般的には、検出の目的とするスペクトルの候補となり得る)が含まれており、その目立つスペクトルの数がそれほど多くはない、というような場合である。一般的には、このような場合、目立つ複数本のスペクトルを抽出し、過去の拍動成分との相関性や、ノイズとの峻別性等に基づく判断を行うことが有効とされている。
【００９９】
但し、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に著しく低い場合等には、一見、信号がきれいな状態であっても、第１処理部６０による拍動／ノイズ分離処理を含む信号処理（参照符号６２〜６６のブロックによってなされる一連の信号処理）では、拍動信号を特定できない場合がある。
【０１００】
また、一見、信号がきれいに見える場合としては、拍動信号が存在せず、全部がノイズであるという場合がある。つまり、スペクトル全体の形状からみると、一見、きれいな信号に見えるが、実際は、目的とする信号成分はノイズに埋もれてほとんど認識することができない状態（つまり、全部がノイズであって目的とする拍動信号がない状態）である場合がある。このような場合、検出の目的である拍動信号が無いのであるから、外乱ノイズと拍動信号との分離を試みても効果がないと言える。
【０１０１】
そこで、第１捕捉処理部７２では、第１捕捉処理を実行する。この第１捕捉処理において、拍動呈示スペクトルの捕捉のためには、上述した第１拍動基線（今回の検出に係る拍動基線）と第２拍動基線（所定時間前の過去の検出に係る拍動基線）との間の相関性に着目することが有効である。
【０１０２】
上述のとおり、第１拍動基線は、今回の周波数解析処理に使用される、「脈波信号」の周波数スペクトルにおける被検体の拍動由来の拍動成分と推定されるスペクトルである。第２拍動基線とは、前回の周波数解析処理に使用された、「脈波信号またはフィルター後信号」の周波数スペクトルに基づいて、拍動呈示スペクトルとして特定された拍動基線である。
【０１０３】
前回検出時の周波数スペクトルにおける第１拍動基線と、今回検出時の周波数スペクトルにおける第２拍動基線との相関を重視して判断したとき、相関性が有るとみることができれば、第１捕捉処理部７２は、第１拍動基線は、検出目的の拍動呈示スペクトルであると判断する。よって、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような場合であっても、第１拍動基線を拍動呈示スペクトルとして捕捉することができる。この場合、第１拍動基線と第２拍動基線との間の相関によって判断していることから、仮に比較的大きなノイズがあったとしても、拍動呈示スペクトルを捕捉することができる。
【０１０４】
また、外乱ノイズ（非周期的ノイズ）との分離は困難であっても、体動（周期的な腕振り等）に由来する周期的な体動信号との区別は可能であることが多い。よって、相関性有りと判断され、かつ、周期的な体動信号と区別できることを、拍動呈示スペクトルの捕捉の条件としてもよい。図２（Ａ）および図２（Ｂ）の例では、この判断基準を採用して、拍動呈示スペクトルを捕捉する。
【０１０５】
図２（Ａ）において、Ｘ１が体動基線であり、Ｘ２が拍動基線であるとする。拍動基線Ｘ２は、体動基線Ｘ１よりも十分に小さく、ノイズに埋もれた状態となっている。
【０１０６】
この場合において、図２（Ｂ）に示すように、まず、０〜４Ｈｚの周波数スペクトルの中から、信号値の大きいもの順に、例えば１０本のスペクトルを抽出する。ここでは、信号値の大きいスペクトルから順に、スペクトルＬ１〜Ｌ１０とする。抽出されたスペクトルの中に、拍動呈示スペクトルが存在する可能が高い。但し、信号値に着目した判断を実行すると、信号値の大きなノイズのスペクトルを拍動呈示スペクトルとして誤検出する可能性が高い。よって、所定時間前の過去の検出に係る拍動基線（第２拍動基線）と、今回の検出に係る第１拍動基線との相関性を重視して、拍動呈示スペクトルを特定する。
【０１０７】
図２（Ｂ）において、第２拍動基線（前回の拍動基線）はＬ０である。ここでは、実際の拍動呈示スペクトルはＬ６であるとする。Ｌ１〜Ｌ１０が、拍動呈示スペクトルの候補であり、各スペクトルＬ１〜Ｌ１０が、今回の検出に係る第１拍動基線（一応、拍動成分と推定できる基線）となる。
【０１０８】
次に、所定時間前の過去（４秒前）に拍動基線と推定され、かつ拍動呈示スペクトルとして特定された拍動基線（第２拍動基線）Ｌ０と、各スペクトルＬ１〜Ｌ１０との周波数軸上における距離を算出する。
【０１０９】
第１拍動基線Ｌ０に最も近い第１拍動基線はＬ６である。ここで、周波数軸上におけるＬ０とＬ６の距離に着目する。Ｌ０とＬ６の距離は「２」である。「距離２」とは、周波数軸上で、第２拍動基線Ｌ０の位置（第２位置）を基準として、第１拍動基線Ｌ６の位置（第１位置）が１本のスペクトルを隔てて隣り合っている場合の距離である。例えば、「基準のスペクトルに対して距離２以内のスペクトル」というとき、基準のスペクトルに対して、２本のスペクトル以内で、低周波数側（左）あるいは高周波数側（右）において隣り合っているスペクトル、という意味である。第１捕捉処理部７２は、基準となる第２拍動基線Ｌ０と、今回の検出に係る第１拍動基線Ｌ６の距離が「距離２以内（＝距離０，距離１，距離２）」であるとき、第１拍動基線Ｌ６を拍動呈示スペクトルとしてもよい。
【０１１０】
例えば、スペクトル１本分が、脈拍数３．７５拍分に相当する。被検体である人が定常的に同じようなリズムで運動しているとき、４秒間に、脈拍数が“２０”以上変化することは極めて希である。ここで、周波数軸上での距離が３であるということ（周波数軸上で、スペクトル（基線）が３本分離れること）は、被検体の脈波数が“１０”（上記の脈拍数“２０”の半分）変化するのに相当する。
【０１１１】
本実施形態では、第１捕捉処理部７２は、４秒間における脈拍数の“１０”の変化は発生しないものとして、判断処理を実行する。上記の脈拍数２０に対して半分の値１０を基準としたのは、判定処理の安全度を加味したためである。
【０１１２】
したがって、図２（Ｂ）の例では、第１捕捉処理部７２は、ピーク順が６番のスペクトルである第１拍動基線Ｌ６を、拍動呈示スペクトルとして特定する。
【０１１３】
以上の手順は、第１捕捉処理の一例であり、変形が可能である。例えば、上述の手順によれば、信号ピークの上位１０本の全部について、第２拍動基線Ｌ２との間の距離判断やスペクトルの大きさの比に基づく判断が必要であるが、次に示す判断を行うことで、第１拍動基線を減らすことができる。すなわち、判断処理量を減らすことができる。
【０１１４】
ここで、体動センサー２０から出力される体動信号ｆのＦＦＴ結果に基づいて、スペクトルＬ１が体動基線（第ｎ体動基線）と推定できる場合、このスペクトルＬ１と、このスペクトルＬ１の、低周波数側（左）ならびに高周波数側（右）に隣り合う各１本のスペクトルＬ２，Ｌ３は、第１拍動基線（今回の検出に係る、拍動成分と推定される基線）から除外することができる。３本のスペクトルＬ１〜Ｌ３は、体動成分である可能性が極めて高いからである。
【０１１５】
これによって、第１拍動基線は、Ｌ４〜Ｌ１０の７本となる。この場合、第１拍動基線の数が減少した分、第２拍動基線Ｌ２との間の距離判断のための、第１捕捉処理部７２の負担が軽減されることになる。なお、後述する図６のステップＳＴ１０，ＳＴ１１では、体動基線（第ｎ体動基線）と、その体動基線の左右に位置する各１本のスペクトルを第１拍動基線から除外した後、第２拍動基線との間の距離判定を実行している。
【０１１６】
（第２捕捉処理）
第２捕捉処理部７４は、上述のとおり、脈波信号ｄに含まれるノイズ量の程度に基づく評価が第１許容レベルよりも低い第２許容レベルである場合、例えば、ノイズ量の程度が第１許容レベルを基準とする場合よりも大きく、脈波信号のきれいさが中程度（まあまあきれい）のときに動作する。
【０１１７】
脈波信号ｄには、ノイズ少（きれい）の場合と比べてかなり目立つ非周期的なノイズ（外乱ノイズ）が多く含まれている。よって、拍動呈示スペクトルの捕捉のためには、拍動基線とその他の基線（外乱ノイズや体動ノイズに起因してスペクトル値が大きくなっている）とを区別するために別の判断処理、すなわち第２捕捉処理が必要となる。
【０１１８】
第２捕捉処理部７４は、第１処理部６０で行われる拍動／ノイズ分離処理（図１の例における、参照符号６２〜６６のブロックによって行われる処理）と同じ処理を実行することができる。但し、これだけでは、拍動信号と外乱ノイズとを十分には区別できない可能性がある。
【０１１９】
そこで、第２捕捉処理部は、さらに、今回の検出時の周波数スペクトルにおける第１拍動基線と、第ｎ体動基線との区別の可否に基づいて拍動呈示スペクトルの捕捉を試みる。
【０１２０】
上述のとおり、第ｎ体動基線は、異なる周波数のｍ本（ｎ、ｍは１以上の整数）の体動基線のうち、一番周波数の低い基線を第１番目としそのｎ倍の周波数に相当する基線を第ｎ体動基線と定義する。例えば、第１波（基本波）の拍動基線が、第ｎ体動基線（１≦ｎ≦ｍ）の位置とは異なる位置にある場合には、第２捕捉処理部は、第１拍動基線が拍動呈示スペクトルであると判断する。この場合、例えば、ｎが２〜４である場合の体動基線は、ｎが１である場合の体動基線よりも信号値が小さく、外乱ノイズ中に潜んでいる場合があるため、第ｎ体動基線（例えばｎ＝１〜４）との区別を試みることは、拍動呈示スペクトルの捕捉精度を高める点で有効である。
【０１２１】
このように、第２捕捉処理部７４による第２捕捉処理は、例えば、第１処理部で行われる拍動／ノイズ分離処理と同じ処理を、フィルター前の脈波信号のＦＦＴ結果に対して実行し、さらに、体動基線（第ｎ体動基線）との区別の可否を慎重に判断する処理ということができる。
【０１２２】
なお、第２捕捉処理においては、第１拍動基線と第２拍動基線との相関性も考慮してもよい。例えば、今回検出時の周波数スペクトルにおける第１拍動基線と、第２拍動基線との間の周波数軸上における距離の観点からの判断を、拍動呈示スペクトルの捕捉の条件としてもよい。
【０１２３】
例えば、今回検出に係る第１拍動基線が、体動センサー出力の周波数分解結果に基づいて特定される体動基線（基本波成分ならびに高調波成分を含む）に重なっていないことを第１条件とし、今回検出に係る第１拍動基線が、周波数軸上で、第２拍動基線に対して所定距離内にあることを第２条件としてもよい。
【０１２４】
本実施形態では、第２捕捉処理部７４は、周波数軸上でスペクトル（基線）が３本分離れるということが発生しないものとして、拍動呈示スペクトルの呈示処理を行う。上述の第１条件と第２条件を満足する場合には、第２捕捉処理部は、第１拍動基線を拍動呈示スペクトルとして捕捉する。これによって、例えば、目立つノイズに引きずられることなく、拍動信号成分の信号レベルが他の信号成分の信号レベルと比べて相対的に低いような拍動呈示スペクトルを特定できる可能性が高まる。
【０１２５】
図３（Ａ）および図３（Ｂ）は、第２捕捉処理部による第２捕捉処理の一例を示す図である。図３（Ａ）および図３（Ｂ）には、０〜４Ｈｚまでの脈波信号をＦＦＴして得られる周波数スペクトルが示されている。図３（Ａ）は第２捕捉処理によって救済が可能である例を示し、図３（Ｂ）は第２捕捉処理によっても救済することができない例を示す。
【０１２６】
第２捕捉処理部７４は、フィルター前の脈波信号ｄについて、ピーク順ソート部６２、相関判断部６４ならびに拍動／ノイズ分離部６６による拍動／ノイズ分離処理を実行する。この後、捕捉部７６が、体動基線（第ｎ体動基線）との峻別（区別）の可否を判断する。捕捉部７６は、さらに、第１拍動基線と第２拍動基線との相関性を考慮してもよい。例えば、今回検出時の周波数スペクトルにおける第１拍動基線と、第２拍動基線との間の周波数軸上における距離の観点からの判断を行ってもよい。
【０１２７】
図３（Ａ）の例では、第２拍動基線がＬ１３であり、今回の検出に係る第１拍動基線がＬ１４である。また、第ｎ体動基線（ｎ＝２）をＬ１５とし、第ｎ体動基線（ｎ＝３）をＬ１６とする。
【０１２８】
図３（Ａ）の例では、第１拍動基線Ｌ１４は、第ｎ体動基線（ｎ＝２，３）Ｌ１５，Ｌ１６と重なってはいない。また、第２拍動基線Ｌ１３と、今回の検出に係る第１拍動基線Ｌ１４との間の、周波数軸上での距離は「距離２」である。よって、第１拍動基線Ｌ１４は、第２拍動基線Ｌ１３と相関があると判断できる。よって、第１拍動基線Ｌ１４を、拍動呈示スペクトルとして捕捉することができる。
【０１２９】
図３（Ｂ）の例では、第２拍動基線がＬ１７であり、今回の検出に係る第１拍動基線がＬ１８である。また、第ｎ体動基線（例えばｎ＝１）をＬ１９とする。
【０１３０】
図３（Ｂ）の例では、第１拍動基線Ｌ１８が、第ｎ体動基線（例えばｎ＝１）Ｌ１９と重なっている。つまり、第１拍動基線Ｌ１８は、第ｎ体動基線（例えばｎ＝１）Ｌ１９と区別がつかない。上述の第１条件を満足しないことから、第１拍動基線Ｌ１８を、拍動呈示スペクトルとして捕捉することができない。
【０１３１】
また、図３（Ｂ）の例では、第１拍動基線Ｌ１８は、第２拍動基線Ｌ１７に対する周波数軸上での距離が３である。すなわち、第１拍動基線Ｌ８は、２本のスペクトルを隔てて、第２拍動基線Ｌ２に隣り合っている。よって、第２条件（距離２以内である）を満足せず、この点でも、第１拍動基線Ｌ１８を、拍動呈示スペクトルとして捕捉することができない。この場合には、拍動呈示スペクトルを検出できない。
【０１３２】
（第２実施形態）
本実施形態では、拍動検出装置の構成と動作をより具体的に説明する。本実施形態では
特に、判断部８０の構成と動作について詳述する。
【０１３３】
図４は、本発明の拍動検出装置の他の例の構成を示す図である。図４において、図１と共通する部分には同じ参照符号を付し、その説明を省略している。
【０１３４】
図４に示される拍動検出装置１００は、脈波センサー１０と、脈波信号蓄積部（４秒分の脈波信号ｄのデータを蓄積する第１バッファメモリー１３および１６秒分の脈波信号ｄのデータを蓄積する第２バッファメモリー１５を有する）１２と、適応フィルター３２および体動成分除去フィルター３４を含むフィルター部３０と、体動センサー（加速度センサーやジャイロセンサー等）２０と、体動信号蓄積部２２と、周波数解析部５０と、脈拍数算出部９０と、表示処理部９２と、表示部９４と、を有する。
【０１３５】
上述のとおり、脈波センサー１０は、例えば、光電脈波センサー及びその原理に基づく脈波センサーである。脈波センサー１０は、拍動信号と、被検体（人や動物）の体動に由来する体動ノイズ信号を含むノイズ信号とが混在した脈波信号ｄを出力する。
【０１３６】
ここで、脈波信号ｄは、例えば、拍動信号成分（定常的成分あるいは周期的成分）と、体動ノイズ成分（定常的あるいは周期的成分）と、外乱ノイズ成分（衝撃ノイズ等の、非定常的あるいは非周期的成分）とを含む。
【０１３７】
脈波センサー１０から出力される脈波信号ｄの、４秒分の信号が、第１バッファメモリー１３に蓄積される。４秒分の脈波信号ｄは、４秒周期で、第２バッファメモリー１５に転送される。第２バッファメモリー１５はＦＩＦＯ（ファーストイン・ファーストアウト）メモリーであり、１６秒分の脈波信号は、４秒分ずつ更新される。１６秒分の脈波信号を蓄積するのは、周波数解析によって拍動成分を特定するとき、ある程度の時間幅で信号の推移を観測し、脈拍の時間的変化を含めて慎重に検討する必要があるからである。
【０１３８】
フィルター部３０は、入力信号に、定常的（周期的）な周波数成分とその他の非定常的（非周期的）な成分が含まれるときに、それらを分離して出力することのできる適応フィルターの一種である。フィルター部３０に含まれる適応フィルター３２によって、脈波信号ｄに含まれる定常的（周期的）な周波数成分とその他の非定常的（非周期的）な成分とを分離することができる。また、体動成分除去フィルター３４によって、例えば、脈波信号ｄに含まれる目立つ体動成分を除去・抑制することができる。
【０１３９】
周波数解析部５０は、信号分配部３９と、第１周波数分解部４０ａ〜第３周波数分解部４０ｃと、タイミング部４３と、第１処理部６０と、第２処理部７０（第１捕捉処理部７２および第２捕捉処理部７４とを含む）と、判断部８０を含む。
【０１４０】
信号分配部３９は、フィルター後信号ｅを第１周波数分解部４０ａに供給し、体動信号ｆを第２周波数分解部４０ｂに供給し、フィルター前の脈波信号ｄを第３周波数分解部４０ｃに供給する。
【０１４１】
タイミング部４３は、第１処理部６０ならびに第２処理部７０の動作タイミングを制御する。なお、第１処理部６０ならびに第２処理部７０の動作は、前掲の実施形態にて説明したとおりである。
【０１４２】
判断部８０は、周波数解析部（ここでは、第３周波数分解部４０ｃ）が生成した周波数スペクトルに基づいて、脈波信号ｄのきれいさに関する判断指標を取得する判断指標取得部８１と、判断指標を用いて、脈波信号ｄに含まれるノイズ量の程度を評価する評価部８３と、被検体の運動状態を判断する運動状態判断部８５と、前回の拍動検出の結果を格納している拍動検出履歴格納部８７と、捕捉処理の実行可否判断部８９と、を有する。
【０１４３】
捕捉処理の実行可否判断部８９は、例えば、評価部８３による評価に基づいて、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断してもよい（第１判断例）。
【０１４４】
また、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、評価部８３による評価が許容範囲にあることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい（第２判断例）。
【０１４５】
また、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、評価部８３による評価が許容範囲にあり、かつ、運動状態判断部８５によって被検体が定常状態にあると判断されることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。
【０１４６】
なお、第１処理部６０と第２処理部７０とを並列に動作させることも可能であるが、第１処理部６０によって拍動呈示スペクトルが特定できた場合には、第２処理部７０による拍動呈示スペクトルの捕捉は不要である。そこで、第２処理部７０の無駄な動作をなくすために、第１処理部６０によって拍動呈示スペクトルの特定ができなかったときに、判断部８０による拍動呈示スペクトルの捕捉処理の可否判断を開始させてもよい。
【０１４７】
（第１判断例）
判断部８０に含まれる捕捉処理の実行可否判断部８９は、脈波信号ｄに含まれるノイズ量の程度の評価に基づいて、第２処理部７０による拍動呈示スペクトル捕捉処理の実行の可否を判断する。なお、脈波信号ｄに含まれるノイズ量の程度の判断指標の具体例、ならびに評価の例については、図５を用いて後述する。
【０１４８】
捕捉処理の実行可否判断部８９は、脈波信号ｄに含まれるノイズ量の程度の評価を参照して、例えば、脈波信号に混入されるノイズ量（例えば外乱ノイズ成分のノイズ量）が所定の基準以下であると判断される場合に、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。換言すれば、脈波信号ｄが、許容レベル以上にきれいであることを第２処理部７０による拍動呈示スペクトル捕捉処理の実行の条件としてもよい。
【０１４９】
第１判断例によれば、外乱ノイズ（非周期的ノイズ）が多い状態では、このような場合には、拍動呈示スペクトルの捕捉処理が実行されない。よって、拍動呈示スペクトルの誤検出の可能性が低減される。第１判断例によれば、脈波信号ｄのノイズ状態を慎重に判断しつつ、第２処理部７０による拍動呈示スペクトルの捕捉処理の可否を判断してもよい。
【０１５０】
（第２判断例）
判断部８０に含まれる捕捉処理の実行可否判断部８９は、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、評価部８３による評価が許容範囲にあることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。
【０１５１】
第２判断例では、周波数解析部（ここでは、第１周波数分解部４０ａまたは第３周波数分解部４０ｃ）による前回の周波数解析によって、拍動呈示スペクトルが検出（特定）されていること、ならびに、評価部８３による、脈波信号ｄに含まれるノイズ量の程度の評価が許容範囲にあることを条件として、第２処理部７０による拍動呈示スペクトルの捕捉処理を実行する。
【０１５２】
上述の第１判断例では、ノイズについての基準によって、第２処理部７０による拍動呈示スペクトルの捕捉処理の可否を判断していた。第２判断例では、前回の拍動信号の検出に成功しているという条件を満足するか否かも判断する。前回の拍動検出の結果は、拍動検出履歴格納部８７に格納されている。捕捉処理の実行可否判断部８９は、拍動検出履歴格納部８７から前回の結果を読み出して参照することによって、前回の拍動信号の検出の成功／不成功を確認することができる。
【０１５３】
上述のとおり、周波数分解部４０（４０ａ〜４０ｃ）は、一定時間毎（例えば４秒毎）に周波数解析処理を行って拍動呈示スペクトルを特定する。第２判断例では、前回の検出に成功していない場合には、第２処理部７０による拍動呈示スペクトルの捕捉処理は実行されない。
【０１５４】
上述のとおり、拍動呈示スペクトルの検出には、所定時間前の過去の拍動信号成分と現在の拍動信号成分との相関性に着目することが有効である。よって、前回の拍動信号の検出に成功しているか否かを条件としてもよい。第２判断例では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【０１５５】
（第３判断例）
第３判断例では、判断部８０に含まれる捕捉処理の実行可否判断部８９は、前回の周波数解析処理において、被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、評価部８３による評価が許容範囲にあり、かつ、運動状態判断部８５によって被検体が定常状態にあると判断されることを条件として、第２処理部７０による拍動呈示スペクトル捕捉処理の実行可と判断してもよい。
【０１５６】
被検体の運動状態が急激に変化している場合には、脈拍等の変化が急である。よって、運動状態判断部８５によって被検体が定常状態にあると判断される、という条件によって、第２処理部７０による拍動呈示スペクトルの捕捉処理の実行の可否を判断してもよい。
【０１５７】
被検体が定常状態にある場合とは、例えば、被検体が一定の速度で歩行している場合のように、体動センサー２０から得られる体動信号の状態が安定している状態、あるいは、被検体が安静状態にある場合等である。例えば、被検体が静止している状態から歩行状態に移行したとき、当初は加速過程にあることから、体動信号のスペクトル（体動スペクトル）は一定しない。被検体の歩行状態が安定する（例えば、一定の速度で歩行する状態となる）と、体動信号のスペクトルが安定し、時間が経過しても、体動スペクトルに大きな変化が生じない状態となる。
【０１５８】
運動状態判断部８５は、例えば、２０秒間（１回の測定が４秒であるとすると、その５倍の期間）にわたって、体動スペクトルに大きな変化が見られない場合に、被検体が定常状態にあると判断する。但し、これは一例であり、これに限定されるものではない。例えば、現在ならびに過去の脈波信号のスペクトルの安定性を判断基準とすることもできる。また、後述するｒ５，ｒ１０という指標によっても、被検体の運動状態を判断することができる。
【０１５９】
運動状態判断部８５による被検体の運動状態の判断結果（定常状態／非定常状態の判断結果）は、捕捉処理の実行可否判断部８９に供給される。第３判断例では、より慎重に拍動呈示スペクトル捕捉処理の可否を判断することができる。
【０１６０】
（脈波信号に含まれるノイズ量の程度（ノイズ量の程度）の判断について）
以下、脈波信号ｄに含まれるノイズ量の程度（ノイズ量の程度）の判断の例について説明する。図５（Ａ）〜図５（Ｃ）は、脈波信号に含まれるノイズ量の程度の判断動作を説明するための脈波信号の波形および周波数スペクトルの一例を示す図である。
【０１６１】
図５（Ａ）〜図５（Ｃ）において、上側には、１６秒間のＦＦＴ前の脈波信号ｄの信号波形が示されている。横軸は時間を示し、縦軸は信号の振幅を示す。また、下側には、０から４Ｈｚの周波数帯域における周波数スペクトルが示されている。横軸は周波数を示し、縦軸はスペクトル値を示す。
【０１６２】
ここで、図５（Ａ）は、ノイズ少（脈波信号のきれいさの程度がきれい）の場合における脈波信号ｄの波形と周波数スペクトルを示し、図５（Ｂ）は、ノイズが中程度（脈波信号のきれいさの程度がまあまあきれい）の場合における脈波信号の波形と周波数スペクトルを示しており、図５（Ｃ）は、脈波信号ｄに多くのノイズが含まれる場合（脈波信号のきれいさの程度がノイジー）における脈波信号ｄの波形と周波数スペクトルを示している。図５（Ａ）〜図５（Ｃ）の各々の比較から明らかなように、脈波信号ｄの波形と周波数スペクトルとは密接に関連しており、脈波信号の波形に対応して、周波数スペクトルの分布状態やスペクトル値が変化する。よって、ＦＦＴによって得られる周波数スペクトルに基づいて、脈波信号ｄにする重畳するノイズの状態（ノイズ量の程度）を推定することが可能である。
【０１６３】
図５の例では、ノイズ量の程度の推定のための指標として、例えば、主要な周波数スペクトルのスペクトル値の比（つまり、基線の高さの比）を用いることができる。具体的には、ｒ５およびｒ１０という指標を用いる（ただし、一例であり、他の統計的指標、例えば、標準偏差等を用いてもよい）。ここで、ｒ５とは、１６秒分の脈波信号の周波数スペクトルの中から、ピーク値の大きさの順に５本のスペクトルを並べたとき（つまり、ソーティングしたとき）、第１番目のスペクトルのスペクトル値（パワー）を分母とし、第５番目のスペクトルのスペクトル値（パワー）を分子とすることによって得られる指標である。
【０１６４】
また、ｒ１０とは、１６秒分の脈波信号の周波数スペクトルの中から、ピーク値の大きさの順に１０本のスペクトルを並べたとき（つまりソーティングしたとき）、第１番目のスペクトルのスペクトル値（パワー）を分母とし、第１０番目のスペクトルのスペクトル値（パワー）を分子とすることによって得られる指標である。指標ｒ５，ｒ１０は、判断部８０に含まれる判断指標取得部８１によって取得される。
【０１６５】
ここでは、一例として、ｒ５＜０．５かつｒ１０＜０．２のときをノイズ少（きれい）とし、ｒ５＞０．７かつｒ１０＞０．５のときをノイズ多（ノイジー）とし、上記いずれでもない場合をノイズが中程度（まあまあ）とする。
【０１６６】
図５（Ａ）の例では、ｒ５＝０．１４かつｒ１０＝０．０８であることから、ノイズ少（きれい）と判断される。また、図５（Ｂ）の例では、ｒ５＝０．５６かつｒ１０＝０．３５であることから、ノイズが中程度（まあまあ）と判断される。図５（Ｃ）の例では、ｒ５＝０．８２かつｒ１０＝０．６２であることから、ノイズ多（ノイジー）と判断される。
【０１６７】
このような判断処理が、判断部８０に含まれる評価部８３によって実行される。なお、被検体の運動状態も、上述の指標ｒ５、ｒ１０を用いて推定することもできる。被検体の運動状態によって脈波信号ｄの波形が変化すると、その変化は、周波数スペクトルの変化となって現れ、周波数スペクトルの変化は、指標ｒ５，ｒ１０に反映されるからである。被検体の運動状態は、上述のとおり、運動状態判断部８５によって実行される。
【０１６８】
（拍動呈示スペクトルの捕捉処理のフローの例について）
図６は、第２処理部による、拍動呈示スペクトルの捕捉処理（第２処理）の手順の一例を示すフローチャートである。
【０１６９】
上述のとおり、例えば、第１処理部６０による拍動呈示スペクトルの特定（第１処理）に失敗した場合に、第２処理部７０による拍動呈示スペクトルの捕捉処理（第２処理）を実行することができる。
【０１７０】
まず、判断部８０は、前回拍動呈示スペクトル（拍動基線）を特定できたかを判断する（ステップＳＴ１）。特定できたとき（ステップＳＴ１：Ｙ）のときは、被検体が定常状態であるかを判断する（ステップＳＴ２）。定常状態であるときは（ステップＳＴ２：Ｙ）のときは、脈波信号に含まれるノイズ量の程度（脈波信号のきれいさの程度）を判断する（ステップＳＴ３）。具体的には、きれいさの程度が、きれい（外乱ノイズ少：第１許容レベル）、まあまあ（外乱ノイズ中程：第２許容レベル）、ノイジー（許容レベル未満）のいずれであるかを判断する。
【０１７１】
ステップＳＴ１でＮのとき、あるいは、ステップＳＴ２でＮのときは、拍動呈示スペクトルの検出失敗とする（ステップＳＴ９）。また、ステップＳＴ３で、ノイジーと判断されたときも、拍動呈示スペクトルの検出失敗とする（ステップＳＴ９）。
【０１７２】
ステップＳＴ３で、外乱ノイズ少（脈波信号のきれいさの程度がきれい：第１許容レベル）と判断された場合には、第１捕捉処理部７２が、第１捕捉処理（図６の左側において、破線で囲んで示されるステップＳＴ１０〜ＳＴ１６）を実行する。ステップＳＴ３で、外乱ノイズ中程度（脈波信号のレベルがまあまあきれい：第２許容レベル）と判断された場合には、第２捕捉処理部７４が、第２捕捉処理（図６の右下において、破線で囲んで示されるステップＳＴ４〜ＳＴ７）を実行する。
【０１７３】
第１捕捉処理が実行される場合、まず、第１捕捉処理部７２は、ピーク順ソート部６２でソートされた結果の順位が１０位以内のスペクトルを抽出し、その１０本の中から、体動基線（複数本あるときは、好ましくは各線）と重なっていない、あるいは、各体動基線の、低周波数側（左）あるいは高周波数側（右）に隣り合う１本を除くという条件を満足するスペクトルをリストアップする（ステップＳＴ１０）。リストアップされたスペクトルを、第１拍動基線（今回の検出に係る拍動基線）とすることができる。
【０１７４】
次に、第１捕捉処理部７２は、リストアップした第１拍動基線毎に、第２拍動基線（前回の検出に係る拍動基線）との間の周波数軸上における距離を算出する（ステップＳＴ１１）。
【０１７５】
ステップＳＴ１２では、距離０の第１拍動基線があるかを検出し、有れば（ステップＳＴ１２：Ｙ）、拍動呈示スペクトルの検出に成功する（ステップＳＴ１３）。
【０１７６】
ステップＳＴ１４では、距離１の第１拍動基線があるかを検出する。第１拍動基線が有る場合（ステップＳＴ１４：Ｙ）であって、第２拍動基線の左右の両側に各１本が存在する場合には、信号値が大きい方の第１拍動基線を選択する（ステップＳＴ１５）。これによって、拍動呈示スペクトルの検出に成功する（ステップＳＴ１３）。
【０１７７】
ステップＳＴ１６では、距離２の第１拍動基線があるかを検出し、有れば（ステップＳＴ１６：Ｙ）、拍動呈示スペクトルの検出に成功する（ステップＳＴ１３）。
【０１７８】
ステップＳＴ１２，ＳＴ１４，ＳＴ１６のいずれかにおいてＮの場合は、拍動呈示スペクトルを検出せずと判断する。（ステップＳＴ１７）。
【０１７９】
次に、第１捕捉処理の手順について説明する。上述のとおり、ステップＳＴ３で、脈波信号に含まれるノイズ量の程度が中程度（脈波信号のきれいさの程度がまあまあきれい：第２許容レベル）と判断された場合には、第２捕捉処理部７４が、第２捕捉処理（図６の右下において、破線で囲んで示されるステップＳＴ４〜ＳＴ７）を実行する。
【０１８０】
第２捕捉処理部７４は、まず、ステップＳＴ４にて、フィルター後信号の最大スペクトル値（最大信号値）が所定値（ここでは相対値１０００）以下であるかを判定する。このステップＳＴ４は必須ではないが、実行するのが好ましい。第２捕捉処理の趣旨は、例えば、総じて信号値が大きくなく、かつ全体のスペクトルの形状として凹凸が少ない場合に、信号レベルが相対的に小さい拍動呈示スペクトルを見つけ出すということにある。よって、例えば、最大信号値が１０００を超えるスペクトルが含まれる場合には、上述の趣旨に反することになる。つまり、ステップＳＴ４でＮの場合には、拍動呈示スペクトルの検出に失敗する（ステップＳＴ９）。
【０１８１】
ステップＳＴ４でＹの場合には、フィルター前の脈波信号に対して第１処理（第１処理部６０で行われる拍動／ノイズ分離のための処理と同等の処理：図１の参照符号６２〜６４のブロックで実行される処理）を実行する（ステップＳＴ５）。
【０１８２】
次に、拍動基線（拍動成分と推定される第１拍動基線）が存在するか否かを判定する（ステップＳＴ６）。ステップＳＴ６でＮの場合は、拍動呈示スペクトルの検出に失敗する（ステップＳＴ９）。
【０１８３】
ステップＳＴ６でＹの場合は、ステップＳＴ７の処理を実行する。ステップＳＴ７では、今回の検出に係る拍動基線（第１拍動基線）が、第ｎ体動基線（例えば、ｎ＝１〜４の各基線）に重なっていないかを判断し（第１条件の判断）、さらに、第１拍動基線は第２拍動基線に対して２本以内で隣り合っている（換言すれば、３本以上離れていない）か、を判断する（第２条件の判断）。これらの２つの条件を満たす場合（Ｙの場合）、拍動呈示スペクトルの検出に成功する（ステップＳＴ８）。Ｎの場合は、拍動呈示スペクトルの検出に失敗する（ステップＳＴ９）。
【０１８４】
（第２捕捉処理の適用例）
図７（Ａ）および図７（Ｂ）は、第２捕捉処理の適用例を示す図である。図７（Ａ）および図７（Ｂ）の例では、体動ノイズが多く混入し、外乱ノイズの混入は少ない。
【０１８５】
図７（Ａ）の上側には、１６秒分のフィルター後信号が示されており、下側にはフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図７（Ｂ）の上側には、１６秒分の脈波信号が示されており、下側には脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図７（Ｂ）の例では、脈波信号に含まれるノイズ量の程度は中程度であり、脈波信号ｄのきれいさの程度は、「まあまあきれい」と評価される。
【０１８６】
図７（Ａ）の例では、フィルター後信号に基づく周波数分解によって得られるスペクトルに含まれる拍動基線Ｓ２は小さく、体動ノイズＮ３と区別する困難性は、拍動基線Ｓ２が大きい場合に比べて高くなる。しかし、図７（Ｂ）に示すように、フィルター前の脈波信号周波数分解によって得られるスペクトルに含まれる拍動基線Ｓ１（第１拍動基線Ｓ１）はかなり大きく、かつ、拍動基線Ｓ１の周波数軸上における位置は、過大な体動ノイズＮ１，Ｎ２の位置から離れている。
【０１８７】
図７（Ｂ）の例では、第２捕捉処理部７４による第２捕捉処理が実行される。すなわち、フィルター前の脈波信号に対して、フィルター後信号に対する拍動／ノイズ分離処理と同様の信号処理を施し、第１拍動基線Ｓ１を、複数本の第ｎ体動基線Ｎ１，Ｎ２（ｎ＝１，２に相当）から分離する。図７（Ｂ）の場合、第１拍動基線Ｓ１は、第ｎ体動基線Ｎ１，Ｎ２と重なっておらず、区別することができる。よって、第１拍動基線Ｓ１を、拍動呈示スペクトルとして捕捉（特定）することができる。また、前回の検出時において、脈波信号に基づいて第２拍動基線の特定に成功しているのならば、その第２拍動基線との相関性を併せて判断すると、判断の精度が向上する。
【０１８８】
（第１捕捉処理の適用例）
図８（Ａ）〜図８（Ｃ）は、第１捕捉処理の適用例を示す図である。図８（Ａ）〜図８（Ｃ）は、体動ノイズの混入は多くはないが、外乱ノイズの混入がある例を示す。この場合、脈波信号に含まれるノイズ量の程度は見かけ上「小さい」、つまり、脈波信号ｄのきれいさの程度は見かけ上「きれい」と評価される。
【０１８９】
図８（Ａ）の上側には、前回の１６秒分のフィルター後信号が示されており、下側には前回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図８（Ｂ）の上側には、今回の１６秒分のフィルター後信号が示されており、下側には今回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図８（Ｃ）の上側には、今回の１６秒分の脈波信号が示されており、下側には今回の脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【０１９０】
図８（Ａ）の例では、前回のフィルター後信号に基づく周波数分解によって得られるスペクトルに含まれる拍動基線（第２拍動基線）Ｓ３の信号値は小さいものの、周囲のノイズよりも信号値は大きく、第１処理部６０による第１処理によって、拍動呈示スペクトルが特定されている。
【０１９１】
但し、図８（Ｂ）の例では、今回のフィルター後信号に基づく周波数分解によって得られるスペクトルでは、拍動基線（第１拍動基線）の信号レベルは、他の基線の信号レベルよりも相対的に低い。
【０１９２】
この場合には、第１捕捉処理部７２による第１捕捉処理が実行される。第１捕捉処理を実行した結果を示すのが図８（Ｃ）の例である。図８（Ｃ）の例では、上述のとおり、脈波信号に含まれるノイズ量の程度は、目立つ体動基線等がないことから、一応、小さい（脈波信号ｄのきれいさの程度はきれい）と判断される。但し、外乱ノイズが多いため、第１拍動基線Ｓ４を検出する難易度は、外乱ノイズが少ない場合に比べて高くなる。
【０１９３】
この場合には、例えば、第１捕捉処理部７２が、第２拍動基線Ｓ３と、今回の検出に係る第１拍動基線Ｓ４との相関性を重視した判断を実行し、その相関性に基づいて、第１拍動基線Ｓ４を拍動呈示スペクトルとして特定する。このとき、第１拍動基線Ｓ４が、体動基線（例えば複数本の第ｎ体動基線）と重なっていないという観点からの判断を併用することによって、判断の精度が向上する。
【０１９４】
図９（Ａ）〜図９（Ｄ）は、外乱ノイズの混入が著しい例（脈波信号に含まれるノイズ量の程度が大きい、すなわち脈波信号のきれいさの程度がノイジーである例）を示す図である。例えば、拍動検出装置を、周囲の物体に衝突させたことによって衝撃ノイズが混入したような場合が該当する。
【０１９５】
図９（Ａ）の上側には、前回の１６秒分のフィルター後信号が示されており、下側には前回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図９（Ｂ）の上側には、前回の１６秒分の脈波信号が示されており、下側には前回の脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。また、図９（Ｃ）の上側には、今回の１６秒分のフィルター後信号が示されており、下側には今回のフィルター後信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。図９（Ｄ）の上側には、今回の１６秒分の脈波信号が示されており、下側には今回の脈波信号のＦＦＴ結果（０〜４Ｈｚまでの周波数スペクトル）が示されている。
【０１９６】
図９（Ａ）〜図９（Ｄ）に示されるいずれの例でも、目立つ外乱ノイズばかりであり、拍動基線Ｓは見あたらない。拍動基線Ｓが仮に存在したとしても、その拍動基線Ｓの信号レベルは、他の基線の信号レベルよりも相対的に低い。よって、この場合には、判断部８０は、脈波信号がノイジーな状態であると判断する。この結果、拍動呈示スペクトルの捕捉処理は実行されない。
【０１９７】
（脈拍計測の結果の一例）
図１０は、腕時計型の脈拍計測装置を被検体の手首に装着し、被検体が歩行（ウォーキング）を行ったときにおける、検出された脈拍数の変化を示す図である。図１０の例では、計測開始（ｔ＝０）から６０秒間の間、被検体（ユーザー）は安静な状態を保ち、６０秒が経過した時点から歩行を開始している。
【０１９８】
計測開始後６０秒〜１２０秒までの期間Ｔ１では、被検体の運動開始に伴って脈拍数が、かなり急激に上昇している。この期間Ｔ１において、脈拍検出に失敗する可能性が高いと言える。すなわち、第１処理部６０による、フィルター後信号に基づく周波数解析処理では、脈拍数の変化が大きい期間Ｔ１において、周波数スペクトル値が小さくなり、かつ個々の周波数スペクトルが分散する結果となり、脈拍数の変化が少ない期間に比べて、特徴的な拍動基線（拍動呈示スペクトル）の検出の難易度が高くなる。
【０１９９】
しかし、本実施形態の拍動検出装置によれば、拍動呈示スペクトルの捕捉処理が実行されることから、期間Ｔ１においても、正しい脈拍数を安定的に計測することが可能である。
【０２００】
（第３実施形態）
図１１（Ａ）および図１１（Ｂ）は、拍動検出装置の、被検体への装着例を示す図である。
【０２０１】
図１１（Ａ）の例は、腕時計型の拍動検出装置の例である。脈波センサー１０および表示部９４を含むベース部４００は、リストバンド３００によって、被検体（ユーザー）の左手首２００に装着されている。
【０２０２】
図１１（Ｂ）の例は、指装着型の拍動検出装置の例である。被検体の指先に挿入するためのリング状のガイド３０２の底部に、脈波センサー１０が設けられている。
【０２０３】
本発明はいずれのタイプの拍動検出装置にも適用可能であるが、腕時計型の拍動検出装置（図１１（Ａ）の例）に適用するのが、より好ましい。図１１（Ａ）の例では、脈波センサー１０が、手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着される。このため、脈波センサー１０から出力される脈波信号の振幅が総じて小さくなる傾向にある。本発明を適用することによって、腕時計型の拍動検出装置における拍動呈示スペクトルの特定性能が向上する。
【０２０４】
このように、本発明の少なくとも一つの実施形態によれば、例えば、拍動検出装置の検出性能を高めることができる。また、例えば、脈波センサーを手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着するタイプの拍動検出装置において、拍動信号を示す拍動呈示スペクトルの特定性能を向上させることができる。
【０２０５】
なお、上記のように本実施形態について詳細に説明したが、本発明の新規事項及び効果から実体的に逸脱しない多くの変形が可能であることは当業者には容易に理解できるであろう。従って、このような変形例はすべて本発明の範囲に含まれるものとする。また、明細書又は図面において、少なくとも一度、より広義又は同義な異なる用語と共に記載された用語は、明細書又は図面のいかなる箇所においても、その異なる用語に置き換えることができる。
【符号の説明】
【０２０６】
１０ 脈波センサー、１１ 体動センサー（加速度センサーやジャイロセンサー等）、
１２ 脈波信号蓄積部、１７ 適応フィルター切り替え部、３０ フィルター部、
４０（４０ａ〜４０ｃ） 周波数分解部、
６０ 第１処理部（後処理部）、６２ ピーク順ソート部、
６４ 相関判定部、６６ 拍動／体動分離部（拍動／ノイズ分離部）、
７０ 第２処理部（捕捉処理部）、７２ 第１捕捉処理部、７４ 第２捕捉処理部、
７６ 捕捉部、８０ 判断部、９０ 脈拍数算出部、９２ 表示処理部、
９４ 表示部、１００ 拍動検出装置
【特許請求の範囲】
【請求項１】
被検体の拍動に由来する拍動信号を検出する拍動検出装置であって、
前記拍動信号と、前記被検体の体動に由来する体動信号を含むノイズ信号とが混在した脈波信号を検出して出力する脈波センサーと、
前記脈波信号から、前記ノイズ信号の少なくとも一部を除去するフィルター部と、
前記フィルター部から出力されるフィルター後信号、または前記脈波センサーから出力される前記脈波信号に基づいて、所定時間毎に前記周波数解析処理を行って、前記拍動信号を示す拍動呈示スペクトルを特定する周波数解析部と、を含み、
前記周波数解析部は、
前記フィルター後信号の周波数スペクトルに基づいて、前記拍動信号と前記ノイズ信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行して前記拍動呈示スペクトルを特定する第１処理部と、
前記脈波信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部と、
前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部と、を含むことを特徴とする拍動検出装置。
【請求項２】
請求項１記載の拍動検出装置であって、
前記判断部は、
前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、
前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、を有し、
前記評価部による評価に基づいて、前記第２処理部による前記拍動呈示スペクトル捕捉処理の実行の可否を判断することを特徴とする拍動検出装置。
【請求項３】
請求項１記載の拍動検出装置であって、
前記判断部は、
前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、
前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、
前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、を有し、
前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、前記評価部による評価が許容範囲にあることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断することを特徴とする拍動検出装置。
【請求項４】
請求項１記載の拍動検出装置であって、
前記判断部は、
前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、
前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、
前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、
前記被検体の運動状態を判断する運動状態判断部と、を有し、
前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、前記評価部による評価が許容範囲にあり、かつ、前記運動状態判断部によって前記被検体が定常状態にあると判断されることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断することを特徴とする拍動検出装置。
【請求項５】
請求項１記載の拍動検出装置であって、
前記判断部は、前記第１処理部が拍動呈示スペクトルを特定することができなかった場合に、前記判断部は、前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断することを特徴とする拍動検出装置。
【請求項６】
請求項５記載の拍動検出装置であって、
前記被検体の体動に伴う体動信号を検出して出力する体動センサーと、
前記体動センサーから出力される体動信号に基づいて周波数分解処理を行う周波数分解部と、を有し、
前記周波数分解部による周波数分解結果に基づいて推定される、前記被検体の体動由来の体動成分を体動基線とし、異なる周波数のｍ本（ｍは１以上の整数）の前記体動基線のうちの第ｎ番目（１≦ｎ≦ｍ）の体動基線を第ｎ体動基線としたとき、
前記第２処理部は、周波数解析処理に使用される前記脈波信号の周波数スペクトルにおける前記被検体の拍動由来の拍動成分と推定される第１拍動基線と、前記所定時間前の周波数解析処理に使用された前記脈波信号または前記フィルター後信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルとして特定された第２拍動基線との相関性を重視した第１捕捉処理を実行して、前記拍動呈示スペクトルの捕捉を試みる第１捕捉処理部と、
前記脈波信号に、前記第１処理部による前記拍動／ノイズ分離処理と同じ処理を施した後、前記第１拍動基線と、前記第ｎ体動基線との区別の可否に基づいて第２捕捉処理を実行して前記拍動呈示スペクトルの捕捉を試みる第２捕捉処理部と、を有し、
前記判断部は、前記拍動呈示スペクトル捕捉処理を実行可と判断した場合、さらに、前記評価部による評価が第１の許容レベルであるか、前記第１許容レベルよりも低い第２許容レベルであるかを判断し、
前記判断部によって前記第１許容レベルであると判断されたときは、前記第１捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行し、前記第２許容レベルであると判断されたときは、前記第２捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行することを特徴とする拍動検出装置。
【請求項７】
請求項６記載の拍動検出装置であって、
前記第１捕捉処理部は、前記第１拍動基線の周波数軸上における第１位置と、前記第２拍動基線の周波数軸上における第２位置との間の距離によって、前記相関性を判断することを特徴とする拍動検出装置。
【請求項８】
請求項６記載の拍動検出装置であって、
前記第２捕捉処理部は、
前記第１拍動基線が前記第ｎ体動基線と重なっていないという第１条件と、前記第１拍動基線が前記第２拍動基線に対して、周波数軸上で所定の距離以内であるという第２条件と、を満足する場合に、前記第１拍動基線を、拍動呈示スペクトルとして特定することを特徴とする拍動検出装置。
【請求項１】
被検体の拍動に由来する拍動信号を検出する拍動検出装置であって、
前記拍動信号と、前記被検体の体動に由来する体動信号を含むノイズ信号とが混在した脈波信号を検出して出力する脈波センサーと、
前記脈波信号から、前記ノイズ信号の少なくとも一部を除去するフィルター部と、
前記フィルター部から出力されるフィルター後信号、または前記脈波センサーから出力される前記脈波信号に基づいて、所定時間毎に前記周波数解析処理を行って、前記拍動信号を示す拍動呈示スペクトルを特定する周波数解析部と、を含み、
前記周波数解析部は、
前記フィルター後信号の周波数スペクトルに基づいて、前記拍動信号と前記ノイズ信号とを分離する拍動／ノイズ分離処理を含む信号処理を実行して前記拍動呈示スペクトルを特定する第１処理部と、
前記脈波信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルの捕捉を試みる拍動呈示スペクトル捕捉処理を実行する第２処理部と、
前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断する判断部と、を含むことを特徴とする拍動検出装置。
【請求項２】
請求項１記載の拍動検出装置であって、
前記判断部は、
前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、
前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、を有し、
前記評価部による評価に基づいて、前記第２処理部による前記拍動呈示スペクトル捕捉処理の実行の可否を判断することを特徴とする拍動検出装置。
【請求項３】
請求項１記載の拍動検出装置であって、
前記判断部は、
前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、
前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、
前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、を有し、
前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、かつ、前記評価部による評価が許容範囲にあることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断することを特徴とする拍動検出装置。
【請求項４】
請求項１記載の拍動検出装置であって、
前記判断部は、
前記周波数解析部が生成した周波数スペクトルに基づいて、前記脈波信号に含まれるノイズ量の程度に関する判断指標を取得する判断指標取得部と、
前記判断指標を用いて、前記脈波信号に含まれるノイズ量の程度を評価する評価部と、
前記所定時間前の拍動検出の結果を格納している拍動検出履歴格納部と、
前記被検体の運動状態を判断する運動状態判断部と、を有し、
前記所定時間前の周波数解析処理において、前記被検体の拍動由来の拍動成分と推定される拍動基線を特定することができ、前記評価部による評価が許容範囲にあり、かつ、前記運動状態判断部によって前記被検体が定常状態にあると判断されることを条件として、前記第２処理部による前記拍動呈示スペクトル捕捉処理を実行可と判断することを特徴とする拍動検出装置。
【請求項５】
請求項１記載の拍動検出装置であって、
前記判断部は、前記第１処理部が拍動呈示スペクトルを特定することができなかった場合に、前記判断部は、前記第２処理部による拍動呈示スペクトル捕捉処理の実行の可否を判断することを特徴とする拍動検出装置。
【請求項６】
請求項５記載の拍動検出装置であって、
前記被検体の体動に伴う体動信号を検出して出力する体動センサーと、
前記体動センサーから出力される体動信号に基づいて周波数分解処理を行う周波数分解部と、を有し、
前記周波数分解部による周波数分解結果に基づいて推定される、前記被検体の体動由来の体動成分を体動基線とし、異なる周波数のｍ本（ｍは１以上の整数）の前記体動基線のうちの第ｎ番目（１≦ｎ≦ｍ）の体動基線を第ｎ体動基線としたとき、
前記第２処理部は、周波数解析処理に使用される前記脈波信号の周波数スペクトルにおける前記被検体の拍動由来の拍動成分と推定される第１拍動基線と、前記所定時間前の周波数解析処理に使用された前記脈波信号または前記フィルター後信号の周波数スペクトルに基づいて、前記拍動呈示スペクトルとして特定された第２拍動基線との相関性を重視した第１捕捉処理を実行して、前記拍動呈示スペクトルの捕捉を試みる第１捕捉処理部と、
前記脈波信号に、前記第１処理部による前記拍動／ノイズ分離処理と同じ処理を施した後、前記第１拍動基線と、前記第ｎ体動基線との区別の可否に基づいて第２捕捉処理を実行して前記拍動呈示スペクトルの捕捉を試みる第２捕捉処理部と、を有し、
前記判断部は、前記拍動呈示スペクトル捕捉処理を実行可と判断した場合、さらに、前記評価部による評価が第１の許容レベルであるか、前記第１許容レベルよりも低い第２許容レベルであるかを判断し、
前記判断部によって前記第１許容レベルであると判断されたときは、前記第１捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行し、前記第２許容レベルであると判断されたときは、前記第２捕捉処理部が前記拍動呈示スペクトル捕捉処理を実行することを特徴とする拍動検出装置。
【請求項７】
請求項６記載の拍動検出装置であって、
前記第１捕捉処理部は、前記第１拍動基線の周波数軸上における第１位置と、前記第２拍動基線の周波数軸上における第２位置との間の距離によって、前記相関性を判断することを特徴とする拍動検出装置。
【請求項８】
請求項６記載の拍動検出装置であって、
前記第２捕捉処理部は、
前記第１拍動基線が前記第ｎ体動基線と重なっていないという第１条件と、前記第１拍動基線が前記第２拍動基線に対して、周波数軸上で所定の距離以内であるという第２条件と、を満足する場合に、前記第１拍動基線を、拍動呈示スペクトルとして特定することを特徴とする拍動検出装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【公開番号】特開２０１２−１７０７０１（Ｐ２０１２−１７０７０１Ａ）
【公開日】平成２４年９月１０日（２０１２．９．１０）
【国際特許分類】
【出願番号】特願２０１１−３６８００（Ｐ２０１１−３６８００）
【出願日】平成２３年２月２３日（２０１１．２．２３）
【出願人】（０００００２３６９）セイコーエプソン株式会社 (51,324)
【Ｆターム（参考）】
【公開日】平成２４年９月１０日（２０１２．９．１０）
【国際特許分類】
【出願日】平成２３年２月２３日（２０１１．２．２３）
【出願人】（０００００２３６９）セイコーエプソン株式会社 (51,324)
【Ｆターム（参考）】
[ Back to top ]