(57)【特許請求の範囲】
【請求項１】
被検体の拍動に由来する拍動信号を検出する拍動検出装置であって、
前記拍動信号と、前記被検体の体動に由来する体動ノイズ信号を含むノイズ信号とが混在した脈波信号を出力する脈波センサーと、
前記脈波信号をフィルタリングするフィルターであって、周波数応答特性を自己適応させる適応フィルターと、前記適応フィルターが継続動作している期間における第１時点において、前記適応フィルターの状態を、前記第１時点よりも前の時点であり、かつ前記適応フィルターが動作を開始した時点を含む第２時点の状態に再構成させる適応フィルター再構成部と、を含む脈波信号フィルタリング部と、
前記脈波信号フィルタリング部から出力される信号に基づいて、所定時間毎に周波数解析処理を行って、前記拍動信号を表す拍動呈示スペクトルを検出する周波数解析部と、を含むことを特徴とする拍動検出装置。
【請求項２】
請求項１記載の拍動検出装置であって、
前記適応フィルターが動作を開始した時点あるいは前記適応フィルターが再構成された直近の時点から、第１時間が経過した第３時点までを、前記適応フィルターの再構成が不要な第１期間とし、前記第３時点から、第２時間が経過した第４時点までを前記適応フィルターの再構成が可能な第２期間としたとき、
前記第２時点は、前記第１期間において設定され、前記第１時点は、前記第２期間において設定されることを特徴とする拍動検出装置。
【請求項３】
請求項２記載の拍動検出装置であって、
前記第１時点は、前記適応フィルターが動作を開始した時点を基準として、所定の周期で設定されることを特徴とする拍動検出装置。
【請求項４】
請求項２または請求項３記載の拍動検出装置であって、
前記適応フィルター再構成部は、前記第１時点よりも前の直近の時点における前記周波数解析処理によって、前記拍動呈示スペクトルが検出されていることを条件として、前記適応フィルターを再構成させることを特徴とする拍動検出装置。
【請求項５】
請求項２または請求項３記載の拍動検出装置であって、
前記脈波信号の周波数解析結果に基づいて、前記ノイズ信号のノイズ量の程度を判定する信号評価部を有し、
前記適応フィルター再構成部は、前記第１時点前の直近の時点における、前記周波数解析処理によって、前記拍動呈示スペクトルが検出されており、かつ、前記信号評価部によって、前記ノイズ量が所定の基準以下であると判定されたことを条件として、前記適応フィルターを再構成させることを特徴とする拍動検出装置。
【請求項６】
請求項２または請求項３記載の拍動検出装置であって、
前記脈波信号の周波数解析結果に基づいて、前記被検体の運動状態を判別する信号評価部を有し、
前記適応フィルター再構成部は、前記第１時点前の直近の時点における、前記周波数解析処理によって、前記拍動呈示スペクトルが検出されており、かつ、前記信号評価部によって、前記被検体の運動状態が安定状態であると判定されたことを条件として、前記適応フィルターを再構成させることを特徴とする拍動検出装置。
【請求項７】
請求項２または請求項３記載の拍動検出装置であって、
前記脈波信号の周波数解析結果に基づいて、前記ノイズ信号のノイズ量の程度ならびに前記被検体の運動状態を判別する信号評価部を有し、
前記適応フィルター再構成部は、前記第１時点前の直近の時点における、前記周波数解析処理によって、前記拍動呈示スペクトルが検出されており、前記信号評価部によって、前記ノイズ量が所定の基準以下であると判定され、かつ、前記被検体の運動状態が安定状態であると判定されたことを条件として、前記適応フィルターを再構成させることを特徴とする拍動検出装置。
【請求項８】
請求項１〜請求項７のいずれかに記載の拍動検出装置であって、
前記脈波信号フィルタリング部は、
前記脈波信号を、所定時間遅延させる遅延処理部と、
前記適応フィルターの係数を更新するフィルター係数更新部と、
減算部と、を有し、
前記適応フィルターは、前記遅延処理部の出力信号の中から、自己相関性が高い第１信号を出力し、
前記減算部は、前記脈波信号から前記第１信号を減算して、前記第１信号よりも自己相関性が低い第２信号を生成し、前記第２信号を前記フィルター係数更新部に供給し、
前記フィルター係数更新部は、前記第２信号が抑制されるように、前記適応フィルターの係数を更新し、
前記適応フィルター再構成部は、前記第１時点において前記適応フィルターを再構成させる場合に、前記適応フィルターの係数を、前記適応フィルターが動作を開始した時点における係数に設定することを特徴とする拍動検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、拍動検出装置等に関する。
【背景技術】
【０００２】
拍動検出装置は、人体の心拍に由来する拍動を検出するための装置であって、例えば、腕、手のひら、手指などに装着される脈波センサーからの信号（脈波信号）から、人体の体動の影響により発生する信号成分（体動影響信号）を雑音として除去し、心拍に由来する信号（拍動信号）を検出する装置である。
【０００３】
人の指や手首に装着するタイプの脈波センサーは、例えば、特許文献１に記載されている。また、脈波センサーから出力される脈波信号に含まれるノイズ成分を、フィルタリングによって除去する技術は、例えば、特許文献２に記載されている。特許文献２に記載される技術では、複数のバンドパスフィルターの中から、その時点における脈拍を示す周波数に近い周波数の信号を通過させるバンドパスフィルターを選択する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００５−１９８８２９号公報
【特許文献２】特開２００７−５４４７１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１に記載される技術では、脈波センサーは、人の指、手のひら、手首等に装着される。指、手のひら、手首等は、身体の中でも特に動きの多い部位である。よって、それらの部位に脈波センサーを装着した場合には、脈波センサーから出力される脈波信号には、多くのノイズが混在する。例えば、装着者が、手、あるいは手の周辺の部位を動かすと、心拍による血流とは独立して別の血流の変化が発生し、脈波センサーがキャッチする信号にノイズが混入する。
【０００６】
また、例えば、装着者が、手を物や装着者自身の別の身体部位へぶつけることによって、心拍による血流変動とは独立して別の血流の変化が発生する。脈波センサーはこの血流変化をキャッチすることから、脈波センサー出力信号にノイズが混入する。
【０００７】
特に、手首型の拍動検出装置（人等の手首に装着するタイプの拍動検出装置）は、指型の拍動検出装置（人等の指に装着するタイプの拍動検出装置）に比べて、得られる脈波信号が弱いことから、十分なノイズ対策が必要となる。
【０００８】
例えば、手首には尺骨、擁骨（トウコツ）などの骨と腱、筋肉が集まっており、指、手、手首等を動かすことによって、腱や筋肉の形状が大きく変化する。この際、血流の変化が生じる。動脈や静脈の血液の流れを観察すると、動脈は、静脈よりも心拍による血流変化がより鮮明に現れ、よって、脈波信号にも、心拍のリズムがより明確に現れる。静脈の血流では、心拍の動きは見にくい。しかしながら、手首外側の皮下組織（浅い部位）には、動脈血管が少ない。つまり、手首には動脈が少なく、また、動脈はより深い位置にある場合が多い。よって、脈波センサーで血流の変化を捕捉しようとしたときに、心拍による血流変化よりも、外的な要因による血流変化のほうが支配的となり易い。そのため、手の動きや手周辺への衝撃が入ると、心拍による血流変化が見えにくくなる可能性が高い。
【０００９】
また、特許文献２に記載される技術では、複数のバンドパスフィルターの中から、その時点における脈拍を示す周波数に近い周波数の信号を通過させるバンドパスフィルターを選択する。但し、このような処理を実現するためには、例えば、ソフトウエア上にて、多数の判断処理を実施する必要があり、処理上の負担が増加し、また、処理に要する時間が長くなるのは否めない。このことは、拍動検出装置の大型化、消費電力の増大につながる。例えば、腕時計程度の大きさの拍動検出装置を実現するためには、装置の処理負担や消費電力を低減しつつ、効果的なノイズ対策を実行することが重要である。
【００１０】
本発明の少なくとも一つの態様によれば、例えば、装置の処理負担を低減しつつ、効果的なノイズ対策を実行できる拍動検出装置を実現することができる。また、例えば、脈波センサーを手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着するタイプの拍動検出装置において、拍動信号を表す拍動呈示スペクトルの特定性能を向上させることができる。
【課題を解決するための手段】
【００１１】
（１）本発明の拍動検出装置の一態様は、被検体の拍動に由来する拍動信号を検出する拍動検出装置であって、前記拍動信号と、前記被検体の体動に由来する体動ノイズ信号を含むノイズ信号とが混在した脈波信号を出力する脈波センサーと、前記脈波信号をフィルタリングするフィルターであって、周波数応答特性を自己適応させる適応フィルターと、前記適応フィルターが継続動作している期間における第１時点において、前記適応フィルターの状態を、前記第１時点よりも前の時点であり、かつ前記適応フィルターが動作を開始した時点を含む第２時点の状態に再構成させる適応フィルター再構成部と、を含む脈波信号フィルタリング部と、前記脈波信号フィルタリング部から出力される信号に基づいて、所定時間毎に周波数解析処理を行って、前記拍動信号を表す拍動呈示スペクトルを検出する周波数解析部と、を含む。
【００１２】
適応フィルターは、入力信号に基づいて、伝達関数を自己適応的に変化させることができるデジタルフィルターである。すなわち、適応フィルターは、周波数応答特性を自己適応させるフィルターである。適応フィルターは一般にデジタル信号処理を行うデジタルフィルターとして実装される。フィルターの所望の性能（例えば、入力信号に含まれるノイズ成分を最小化するといった性能）を維持するために、所定のアルゴリズム（例えば最適化アルゴリズム）が使用される。所定のアルゴリズムには、フィルタリング後の信号に基づいて得られる信号がフィードバックされ、そのアルゴリズムによる適応処理によってフィルター係数が適応的に変化し、その結果として、適応フィルターの周波数応答特性が変化する。適応フィルターが継続的に動作している期間では、アルゴリズムによる適応処理が継続的（断続的）に行われ、適応フィルターの所望の性能が維持される。
【００１３】
但し、適応フィルターが、ある程度の時間にわたって継続的に動作すると、適応フィルターの性能が、所望のレベルよりも低下する場合がある。これは、脈波センサーから出力される脈波信号に、拍動信号（定常的成分）と、被検体の体動に由来する体動ノイズ信号を含むノイズ信号（非定常的成分）とが混在していることに起因する。つまり、適応フィルターの性能は、拍動信号（定常的成分）に追従して変化するが、その一方で、ノイズ信号にも追従する可能性がある。拍動信号ならびにノイズ成分のレベル（信号振幅）は、時間経過と共に変動し、また、例えば、突発的かつ非定常的な変動もあり得る。よって、例えば、突発的に拍動成分のレベルが低下し、同時に、ノイズ成分のレベルが上昇した場合には、適応フィルターの性能は、ノイズ成分に追従して変化する。適応フィルターは、入力信号の中から、自己相関性の高い信号を選別して出力することから、時間経過と共に、少しずつノイズへの追従が進行し、徐々にノイズ信号も通過させるようになる。この場合には、ある程度の時間が経過すると、適応フィルターの性能はかなり低下して、正常な状態への復帰がむずかしくなる。
【００１４】
例えば、腕時計型の拍動検出装置は、被検体の手首に装着して使用され、例えば、被検体の経時的な動作状態の変化を記録するために、長時間にわたって継続的に使用される場合がある。よって、長時間にわたって継続的に使用しても、適応フィルターの性能が低下しないようにすることが好ましい。
【００１５】
そこで、本態様では、適応フィルター再構成部を設けて、適切なタイミングで、適応フィルターを再構成させるようにした。すなわち、適応フィルター再構成部は、脈波信号フィルタリング部が継続的に動作している期間における第１時点において、適応フィルターの状態を、第１時点よりも前の時点であり、かつ適応フィルターが動作を開始した時点を含む第２時点の状態に再構成させる。
【００１６】
ここで、フィルターの再構成には、例えば、フィルター係数の初期化（フィルター係数の値を、動作開始時点のフィルター係数値に戻すこと）や、フィルター係数の値を、フィルターの動作期間中における過去の時点（例えば、好ましいフィルター係数が得られていた時点）のフィルター係数値に設定することが含まれる。すなわち、第２時点は、動作開始時点を含む、第１時点よりも前の時点である。また、例えば、一対の適応フィルターを用意し、一方のフィルターから他方のフィルターに切り換えるような場合もフィルターの再構成に含まれる。
【００１７】
仮に、適応フィルターの継続的な使用によってフィルターの性能が低下したとしても（例えば、ノイズ成分への追従が支配的になることによって、拍動信号に追従しないような場合が生じたとしても）、適応フィルターの再構成を実行することによって、適応フィルターは、再構成後に、新たに適応処理を開始することができ、この結果、拍動信号を捕捉するように自己適応がすすむ。よって、適応フィルターは、拍動信号に追従する正常な状態に復帰することができる。
【００１８】
本態様では、適応フィルターを使用することから、複数のバンドパスフィルターを使用する必要がない。また、フィルターの再構成は、フィルター係数の初期化等によって実現でき、装置の処理負担の増加は少なく、また、装置の消費電力の増加も特に問題とならない。
【００１９】
よって、本態様によれば、例えば、装置の処理負担を低減しつつ、効果的なノイズ対策を実行できる拍動検出装置（例えば、長時間の使用に耐える拍動検出装置）を実現することができる。また、例えば、脈波センサーを手首外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着するタイプの拍動検出装置において、拍動信号を表す拍動呈示スペクトルの特定性能を向上させることが可能となる。
【００２０】
（２）本発明の拍動検出装置の他の態様では、前記適応フィルターが動作を開始した時点あるいは前記適応フィルターが再構成された直近の時点から、第１時間が経過した第３時点までを、前記適応フィルターの再構成が不要な第１期間とし、前記第３時点から、第２時間が経過した第４時点までを前記適応フィルターの再構成可能な第２期間としたとき、前記第２時点は、前記第１期間において設定され、前記第１時点は、前記第２期間において設定される。
【００２１】
本態様では、適応フィルターの再構成が不要な第１期間と、適応フィルターの再構成が可能な第２期間とを設定し、適応フィルターの再構成は、第２期間において行う。また、適応フィルターの再構成によって、適応フィルターの状態は、過去の第１時点の状態に変化するが、この第１時点は、第１期間において設定される。適応フィルターは、再構成後に、好ましい状態から、新たに適応処理を開始することができ、この結果、拍動信号を捕捉するように自己適応がすすむ。よって、適応フィルターは、適切な性能を、継続的に維持することができる。
【００２２】
なお、下記の（４）〜（７）の態様では、適応フィルターの再構成の可否を、所定の条件によって判定することとしている。この場合、第２期間（適応フィルターの再構成可能期間）において、所定の条件が充足されない場合があり得る。この場合には、第２期間の終点である第４時点に達した時点で、所定条件の充足／非充足に関係なく、強制的に適応フィルターを再構成することができる。このようにすれば、適応フィルターの再構成ができない期間が長く続くというような、好ましくない状況が生じない。
【００２３】
（３）本発明の拍動検出装置の他の態様では、前記第１時点は、前記適応フィルターが動作を開始した時点を基準として、所定の周期で設定される。
【００２４】
本態様では、適応フィルターの再構成を、適応フィルター（脈波信号フィルタリング部）の動作開始時点を基準として、所定の周期で周期的（定期的）に実行する。本態様では、適応フィルターの再構成のタイミングは、例えば、タイマーによって簡単に管理することができ、実現が容易である。
【００２５】
（４）本発明の拍動検出装置の他の態様では、前記適応フィルター再構成部は、前記第１時点よりも前の直近の時点における前記周波数解析処理によって、前記拍動呈示スペクトルが検出されていることを条件として、前記適応フィルターを再構成させる。
【００２６】
適応フィルターの再構成は、適応フィルターの再構成が可能な第２期間において行われるが、本態様では、フィルターの再構成は第２期間において、いつでも行えるというものではなく、所定の条件を満たす場合に行うことができる。すなわち、本態様では、脈波周波数解析部による前回の周波数解析によって、拍動呈示スペクトルが検出されていることが、適応フィルター再構成の条件となる。
【００２７】
すなわち、脈波周波数解析部は、脈波信号フィルタリング部から出力されるフィルタリング後の信号に基づいて、所定時間毎（例えば４秒毎）に周波数解析処理を行って拍動呈示スペクトルを特定する。適応フィルター再構成部は、第１時点においてフィルター再構成処理を実行しようとする場合、第１時点前の直近の時点における周波数解析（つまり、第１時点を基準とした場合における前回の周波数解析処理）において、拍動呈示スペクトル（拍動の周期及びその信号強度を示す周波数スペクトル）の検出に成功しているかを確認し、その条件を満たすときに、第１時点において適応フィルターの再構成を実行する。
【００２８】
適応フィルターを再構成すると、適応フィルターは過去の第２時点の状態に戻って、新たに適応動作を開始する。この新たな動作の開始時点において、例えば、脈波信号に混入しているノイズが少なく、脈波信号の波形のきれいさの程度が許容範囲である場合には、適応フィルターは、適切な自己適応が可能である。一方、脈波信号に混入しているノイズが多く、脈波信号の波形のきれいさの程度が、許容範囲外である場合には、適応フィルターはノイズに追従する可能性が高い。適応フィルターの再構成は、その時点で捕捉している自己相関性が高いと考えられる信号への追従を断つ処理でもあるため、再構成の直後における拍動呈示スペクトルの捕捉に失敗すると、例えば、第１時点の直前まで検出できていた拍動信号成分をカットしてしまうという結果となり、かえって悪い結果となる。
【００２９】
ここで、前回の周波数解析にて、拍動信号成分の検出に成功している場合には、第１時点における脈波信号のきれいさの程度が、許容範囲内であると推定することができる。よって、本態様では、前回の周波数解析によって拍動呈示スペクトルが検出されていることを、適応フィルター再構成の条件とした。
【００３０】
（５）本発明の拍動検出装置の他の態様では、前記脈波信号の周波数解析結果に基づいて、前記ノイズ信号のノイズ量の程度を判定する信号評価部を有し、前記適応フィルター再構成部は、前記第１時点前の直近の時点における、前記周波数解析処理によって、前記拍動呈示スペクトルが検出されており、かつ、前記信号評価部によって、前記ノイズ量が所定の基準以下であると判定されたことを条件として、前記適応フィルターを再構成させる。
【００３１】
本態様では、上記（４）の態様における条件に加えて、信号評価部によって、脈波信号に混入されるノイズ量（外乱ノイズ成分のノイズ量）が所定の基準以下であると判定されたことを、フィルター再構成の条件とした。信号評価部によるノイズ量の判定によって、脈波信号のきれいさの程度を、より正確に推定することができる。
【００３２】
つまり、外乱ノイズの少ない状態のときにフィルターを再構成することによって、再構成後の適応フィルターによって適切な自己適応が実行され、拍動成分をノイズ成分と区別することができる。よって、適応フィルターの再構成の可否を、より的確に判断することができる。
【００３３】
（６）本発明の拍動検出装置の他の態様では、前記脈波信号の周波数解析結果に基づいて、前記被検体の運動状態を判別する信号評価部を有し、前記適応フィルター再構成部は、前記第１時点前の直近の時点における、前記周波数解析処理によって、前記拍動呈示スペクトルが検出されており、かつ、前記信号評価部によって、前記被検体の運動状態が安定状態であると判定されたことを条件として、前記適応フィルターを再構成させる。
【００３４】
被検体（人や動物）の運動状態は、絶えず変化する可能性がある。定常的な運動状態では、拍動信号も安定しているが、運動状態が急に変化したとき（運動状態の変化の過渡期）においては、拍動信号の振幅や周波数が不安定化し、例えば突発的に大きく変化する場合もあり得る。このような場合に、適応フィルターの再構成を実行すると、拍動信号成分のみならず、ノイズ成分に対しても適応処理の追従をする可能性が高くなる。
【００３５】
よって、本態様では、上述の（４）の態様の条件に、さらに、信号評価部による周波数解析の結果として、被検体の運動状態が安定状態であることが確認されるという条件を追加（加重）する。よって、適応フィルターの再構成の可否を、より的確に判断することができる。
【００３６】
例えば、被検体である人の運動状態が、歩行中、ジョギング中などの定常的な運動をしている状態、もしくは安静状態のいずれかに該当しない状態である場合（例えば球技をしている、体操している状態等の場合）には、適応フィルターの再構成が実行されない。よって、拍動検出の失敗や誤検出の発生を低減することができる。
【００３７】
（７）本発明の拍動検出装置の他の態様では、前記脈波信号の周波数解析結果に基づいて、前記ノイズ信号のノイズ量の程度ならびに前記被検体の運動状態を判別する信号評価部を有し、前記適応フィルター再構成部は、前記第１時点前の直近の時点における、前記周波数解析処理によって、前記拍動呈示スペクトルが検出されており、前記信号評価部によって、前記ノイズ量が所定の基準以下であると判定され、かつ、前記被検体の運動状態が安定状態であると判定されたことを条件として、前記適応フィルターを再構成させる。
【００３８】
本態様では、適応フィルターの再構成の条件として、上述の（６）の態様の条件に、さらに、上述の（５）の条件を追加（加重）する。これによって、適応フィルターの再構成の可否を、さらに的確に判断することができる。
【００３９】
（８）本発明の拍動検出装置の他の態様では、前記脈波信号フィルタリング部は、前記脈波信号を、所定時間遅延させる遅延処理部と、前記適応フィルターの係数を更新するフィルター係数更新部と、減算部と、を有し、前記適応フィルターは、前記遅延処理部の出力信号の中から、自己相関性が高い第１信号を出力し、前記減算部は、前記脈波信号から前記第１信号を減算して、前記第１信号よりも自己相関性が低い第２信号を生成し、前記第２信号を前記フィルター係数更新部に供給し、前記フィルター係数更新部は、前記第２信号が抑制されるように、前記適応フィルターの係数を更新し、前記適応フィルター再構成部は、前記第１時点において前記適応フィルターを再構成させる場合に、前記適応フィルターの係数を、前記適応フィルターが動作を開始した時点における係数に設定する。
【００４０】
本態様では、脈波信号フィルタリング部は、遅延処理部と、適応フィルターと、フィルター係数更新部と、減算部とを有する。フィルター係数更新部は、減算部から出力される第２信号（相関性が低い信号であり、エラー信号と呼ばれることがある）が抑制されるように（例えば最小化されるように）、フィルター係数を更新する。また、本態様では、適応フィルター再構成部は、適応フィルターの再構成を実行させる場合には、フィルター係数を初期化する。すなわち、フィルター係数を、適応フィルターの動作開始時点におけるフィルター係数値に設定する（例えば、フィルター係数値をオールゼロとする）。
【００４１】
本態様によれば、脈波信号フィルタリング部に含まれる適応フィルターを、適切なタイミングで再初期化することができる。よって、長時間の連続計測を行う場合であっても、適応フィルターの性能を、適切なレベルに維持することができる。これにより拍動検出の失敗や誤検出の発生を低減することができる。
【００４２】
このように、上述の本発明の態様の少なくとも一つによれば、例えば、脈波信号フィルタリング部を構成する適応フィルターを適切なタイミングで再構成することによって、適応フィルターの性能を適切なレベルに維持することができる。また、例えば、適応フィルターを有効活用することから、フィルター構成を簡素化でき、また、適応フィルターの再構成を行うことは比較的容易であり、したがって、装置の処理負担を低減しつつ、効果的なノイズ対策を実行できる拍動検出装置を実現することができる。
【図面の簡単な説明】
【００４３】
【図１】本発明の拍動検出装置の一例の構成を示す図
【図２】図２（Ａ）〜図２（Ｃ）は、脈波センサーならびに拍動検出装置の構成の一例を示す図
【図３】図３（Ａ）〜図３（Ｃ）は、脈波信号のきれいさの程度（ノイズ量の程度）の判定動作を説明するための、ＦＦＴ前の信号の波形および周波数スペクトルの一例を示す図
【図４】被検体の運動状態の判定処理の一例を示すフローチャート
【図５】適応フィルターの再構成が実行されるタイミングの一例を示す図
【図６】適応フィルターの再構成処理の手順例を示すフローチャート
【図７】適応フィルタリング処理（前処理）ならびに周波数解析処理（後処理）の手順の一例を示すフローチャート
【図８】図８（Ａ）〜図８（Ｃ）は、実際の処理における脈波信号波形と、ＦＦＴによって得られる周波数スペクトルの一例を示す図
【図９】図９（Ａ）および図９（Ｂ）は、適応フィルターの再構成処理を実行する例と、しない場合の例（比較例）の各々におけるＳ／Ｎ指標（ＳＮ３）の、時間経過に対応した変動を調べた結果を示す図
【発明を実施するための形態】
【００４４】
以下、本実施形態について説明する。なお、以下に説明する本実施形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また、本実施形態で説明される構成の全てが、本発明の必須構成要件であるとは限らない。
【００４５】
（第１実施形態）
（全体構成例）
図１は、本発明の拍動検出装置の一例の構成を示す図である。図１に示される拍動検出装置１００は、被検体（人や動物を含む）の拍動に由来する拍動信号、拍動信号に対応する心拍等の生体情報等を検出するセンサー装置の一種である。
【００４６】
ここで、拍動とは、医学的には心臓のみならず内臓一般の周期的な収縮、弛緩が繰り返された場合に起こる運動のことをいう。ここでは、心臓が周期的に血液を送るポンプとしての動きを拍動と呼ぶ。なお、心拍数とは、１分間の心臓の拍動の数をいう。また、脈拍数は、末梢血管における脈動の数をいう。心臓が血液を送り出す際に、動脈に脈動が生じるので、この回数を数えたものを脈拍数あるいは単に脈拍と呼ぶ。腕で脈を計測する限りは、医学的には心拍数とは呼ばずに脈拍数と呼ぶのが通常である。また、以下の説明では、体動という用語が使用される。体動とは、広い意味では、体を動かすことすべてを意味する。また、被検体の定常的、周期的な体動を狭義の体動といってもよい。狭義の体動は、例えば、歩行・ジョギングなどに伴う定常的、周期的な腕（脈拍計の装着部位近辺）の動きである。
【００４７】
図１に示される拍動検出装置１００は、脈波センサー１０と、脈波信号蓄積部（４秒分の脈波信号ｄのデータを蓄積する第１バッファメモリー１３および１６秒分の脈波信号ｄのデータを蓄積する第２バッファメモリー１５を有する）１２と、適応フィルター再構成部１７を含む脈波信号フィルタリング部２００と、体動センサー（加速度センサーやジャイロセンサー等）１１と、体動ノイズ除去処理部３００と、脈波周波数解析部４００と、検出結果等を表示する表示部（液晶パネル等を含む）６００と、を有する。
【００４８】
脈波センサー１０は、例えば、光電脈波センサー及びその原理に基づく脈波センサーである。脈波センサー１０は、拍動信号と、被検体（人や動物）の体動に由来する体動ノイズ信号を含むノイズ信号とが混在した脈波信号ｄを出力する。
【００４９】
ここで、脈波信号ｄは、例えば、拍動成分信号（定常的成分あるいは周期的成分）と、体動ノイズ成分（定常的あるいは周期的成分）と、外乱ノイズ成分（衝撃ノイズ等の、非定常的あるいは非周期的成分）とを含む。
【００５０】
脈波センサー１０から出力される脈波信号ｄの、４秒分の信号が、第１バッファメモリー１３に蓄積される。４秒分の脈波信号ｄは、４秒周期で、第２バッファメモリー１５に転送される。第２バッファメモリー１５はＦＩＦＯ（ファーストイン・ファーストアウト）メモリーであり、１６秒分の脈波信号は、４秒分ずつ更新される。１６秒分の脈波信号を蓄積するのは、周波数解析によって拍動成分を特定するとき、ある程度の時間幅で信号の推移を観測し、相関の有無等を慎重に検討する必要があるからである。
【００５１】
脈波信号フィルタリング部２００は、入力信号に、定常的な周波数成分とその他の非定常な成分が含まれるときに、それらを分離して出力することのできる適応フィルターの一種である。
【００５２】
脈波信号フィルタリング部２００は、脈波信号を、所定時間（ここでは１サンプル時間）遅延する遅延処理部２０１と、遅延処理部２０１の出力信号の中から、自己相関性が高い第１信号ｙを選別して出力する適応フィルター２０２と、適応フィルター２０２の係数を更新するフィルター係数更新部（ｎＬＭＳフィルター係数更新部）２１０と、脈波信号ｄから第１信号ｙを減算して、第１信号ｙよりも自己相関性が低い第２信号ｅ（＝ｄ−ｙ）を生成し、第２信号ｅをフィルター係数更新部（ｎＬＭＳフィルター係数更新部）２１０に供給（フィードバック）する減算部２０４と、を有している。
【００５３】
なお、遅延処理部２０１と、適応フィルター２０２と、フィルター係数更新部（ｎＬＭＳフィルター係数更新部）２１０と、を含む機能ブロックは、適応線スペクトル強調器と呼ばれることがある。
【００５４】
ｎＬＭＳフィルター係数更新部２１０は、第２信号ｅの値が抑制されるように（例えば、最小化されるように）、適応フィルターの係数（ここでは、正規化最小平均自乗係数（ｎＬＭＳ係数））を、適応的に更新する。
【００５５】
また、第１信号ｙおよび第２信号ｅ（＝ｄ−ｙ）の各々に、ゲイン係数ｈ１，ｈ２の各々を乗算し（この処理は、ゲイン調整部２０６に含まれる係数乗算器２０７ａ，２０７ｂによって実行される）、その後、両信号は、合算器２０８によって合算される。例えば、ｈ１≧１．０に設定され、ｈ２＜１．０に設定される。これによれば、衝撃による影響を軽減するとともに、拍動成分や体動成分の急変化への追従性を高めることができる。例えば、安静中の脈拍が６０だった人が速いぺースのジョギングを開始し、脈拍が１５０まで急上昇していくようなとき、適応フィルター２０２の追従性が、脈拍の上昇より遅れてしまうと、適応フィルター２０２は、急上昇中の心拍成分の信号を遮断してしまう可能性がある。これを排除することができる。
【００５６】
適応フィルター２０２は、入力信号に基づいて、伝達関数を自己適応的に変化させることができる。適応フィルター２０２は、一般にデジタル信号処理を行うデジタルフィルターとして実装される。フィルターの所望の性能（例えば、入力信号に含まれるノイズ成分を最小化するといった性能）を維持するために、所定のアルゴリズム（例えば最適化アルゴリズム）が使用される。所定のアルゴリズムには、フィルタリング後の信号に基づいて得られる信号がフィードバックされ、そのアルゴリズムによる適応処理によってフィルター係数が適応的に変化し、その結果として、適応フィルターの周波数応答特性が変化する。適応フィルターが継続的に動作している期間では、アルゴリズムによる適応処理が継続的（断続的）に行われ、適応フィルター２０２の所望の性能が維持される。
【００５７】
但し、適応フィルター２０２が、ある程度の時間にわたって継続的に動作すると、適応フィルター２０２の性能が、所望のレベルよりも低下する場合がある。これは、脈波センサー１０から出力される脈波信号ｄに、拍動信号（定常的成分）と、被検体の体動に由来する体動ノイズ信号を含むノイズ信号（非定常的成分）とが混在していることに起因する。つまり、適応フィルター２０２の性能は、拍動信号（定常的成分）に追従して変化するが、その一方で、ノイズ信号にも追従する可能性がある。拍動信号ならびにノイズ成分のレベル（信号振幅）は、時間経過と共に変動し、また、例えば、突発的かつ非定常的な変動もあり得る。よって、例えば、突発的に拍動成分のレベルが低下し、同時に、ノイズ成分のレベルが上昇した場合には、適応フィルター２０２の性能は、ノイズ成分に追従して変化する。適応フィルター２０２は、基本的には、入力信号の中から、自己相関性の高い信号を選別して出力することから、時間経過と共に、少しずつノイズへの追従が進行し、徐々にノイズ信号も通過させるようになる。この場合、ある程度の時間が経過すると、適応フィルター２０２の性能はかなり低下して、正常な状態への復帰がむずかしくなることがある。
【００５８】
そこで、本実施形態では、適応フィルター再構成部１７が、適切なタイミングで、適応フィルター２０２の再構成を実行する。図１に示されるように、適応フィルター再構成部１７は、ＦＦＴ（高速フーリエ変換）を実行する高速フーリエ変換部２１２と、必要に応じて設けることができる高速フーリエ変換部２１３，２１５と、信号評価部２１４と、タイミング判定部（タイマー）２１６と、フィルター再構成処理部２１８と、を有する。
【００５９】
高速フーリエ変換部２１２は、脈波信号ｄを高速フーリエ変換する。なお、高速フーリエ変換部２１３、２１５が設けられる場合、高速フーリエ変換部２１３は、体動センサー１１から出力される加速度信号（Ｘ方向成分）を高速フーリエ変換し、高速フーリエ変換部２１５は、体動センサー１１から出力される加速度信号（Ｙ方向成分）を高速フーリエ変換する。
【００６０】
信号評価部２１４は、適応フィルター２０２の再構成の可否を判定するための基礎情報としての、脈波信号ｄのきれいさの程度（すなわちノイズ量の程度）を判定する、信号のきれいさ判定部（ノイズ量の程度判定部）２１７と、被検体の運動状態を判定する運動状態判定部２１９と、を有する。信号のきれいさ判定部（ノイズ量の程度判定部）２１７ならびに運動状態判定部２１９の動作については、図２および図３を参照して後述する。
【００６１】
タイミング判定部（タイマー）２１６は、適応フィルター２０２の再構成のタイミングを決定する上で必要となるタイミング情報を出力する。フィルター再構成処理部２１８は、時期的な条件を満足し、かつ、脈波信号のきれいさの程度等に関係する、所定の実体的条件を満足する場合に、適応フィルター２０２を再構成する。
【００６２】
フィルター再構成処理部２１８には、脈波周波数解析部４００に含まれる拍動成分特定部４０８から、拍動呈示スペクトル（拍動成分）の検出ができたか否かの情報が、拍動成分の検出処理が行われる毎に供給される。フィルター再構成処理部２１８には、メモリー（不図示）が設けられており、拍動呈示スペクトルの検出の成功／失敗の情報が、履歴情報として、そのメモリーに蓄積される。フィルター再構成処理部２１８は、例えば、タイミング判定部（タイマー）２１６からのタイミング情報や、前回の検出における拍動呈示スペクトルの検出の成功／失敗の情報等に基づいて、適応フィルター２０２の再構成の可否を判定し、可能と判定したときに、適応フィルター２０２を再構成する（例えば、フィルター係数をオールゼロとして初期化する）。
【００６３】
例えば、フィルター再構成処理部２１８は、脈波信号フィルタリング部２００が継続的に動作している期間における第１時点において、適応フィルター２０２の状態を、第１時点よりも前の時点であり、かつ適応フィルター２０２が動作を開始した時点を含む第２時点の状態に変化させる。ここで、適応フィルター２０２の再構成には、例えば、フィルター係数の初期化（フィルター係数の値を、動作開始時点のフィルター係数値に戻すこと）や、フィルター係数の値を、フィルターの動作期間中における過去の時点（例えば、好ましいフィルター係数が得られていた時点）のフィルター係数値に設定することが含まれる。すなわち、第２時点は、動作開始時点を含む、第１時点よりも前の時点である。また、例えば、一対の適応フィルターを用意し、一方のフィルターから他方のフィルターに切り換えるような場合もフィルターの再構成に含まれる。
【００６４】
適応フィルター２０２の再構成を実行することによって、適応フィルターを継続使用することで、仮に、フィルターの性能が低下したとしても（例えば、ノイズ成分への追従が支配的になることによって、拍動信号に追従しないような場合が生じたとしても）、適応フィルターは、再構成後に、新たに適応処理を開始することができ、この結果、拍動信号を捕捉するように自己適応がすすむ。よって、適応フィルター２０２は、拍動信号に追従する正常な状態に復帰することができる。なお、詳細については後述する。
【００６５】
また、図１において、体動センサー１１は、被検体の体動、例えば、歩行やジョギングなどに伴う定常的、周期的な腕（脈拍計の装着部位近辺）の動きを検出するセンサーであり、例えば、加速度センサーやジャイロセンサーを含むことができる。
【００６６】
また、体動ノイズ除去処理部３００は、体動ノイズ除去用の、第１適応フィルター３０２および第２適応フィルター３０４を有する。体動ノイズ（体動ノイズ成分）は、脈波信号（より正確には、合算部２０８の出力信号）に含まれる、人間等の定常的な運動や動作（体動）に起因して生じた血管の容積変化を示すノイズ成分である。腕や指に装着する脈拍計の場合、歩行中やジョギング中における腕振りの影響で、その腕振りのリズムに合わせて血管に容積変化が生じる。人間が定常的な動作をすることにより、脈波信号は、その動作の周波数を持った成分信号となる。体動ノイズ成分は、脈波センサー装着部位近辺に装着した加速度センサーが出力する信号の波形と相関性が高いことがわかっている。
【００６７】
また、脈波周波数解析部４００は、体動ノイズ除去後の脈波信号が入力される高速フーリエ変換部４０２と、体動センサー１１からの加速度信号（Ｘ軸方向成分）が入力される高速フーリエ変換部４０４と、体動センサー１１からの加速度信号（Ｙ軸方向成分）が入力される高速フーリエ変換部４０６と、拍動成分特定部（拍動呈示スペクトル特定部）４０８と、脈拍数算出部４１０と、を有し、必要に応じて、さらに運動状態判定部４１２を有することができる。
【００６８】
拍動成分特定部（拍動呈示スペクトル特定部）４０８は、ＦＦＴ後の、１６秒分の脈波信号について４秒毎に周波数解析を行い、スペクトル値やスペクトルの分布等に基づいて、過去に得られた拍動成分との相関性等を検討し、拍動呈示スペクトルを特定する。拍動呈示スペクトルとは、一定期間の拍動成分信号のＦＦＴ結果として得られる周波数スペクトルのうち、拍動の周期およびその信号強度を示す周波数スペクトルである。なお、体動呈示スペクトルとは、一定期間の体動ノイズ成分信号のＦＦＴ結果として得られるスペクトルのうち、体動（例えば歩行中の腕振）の周期およびその信号強度を示す周波数スペクトルである。
【００６９】
脈拍数算出部４１０は、脈拍数を算出する。周波数軸上における拍動呈示スペクトルの位置（周波数）が定まれば、そのスペクトルの位置に対応して、脈拍数が一義的に定まる。運動状態判定部４１２は、脈波信号や体動信号に基づいて、被検体の運動状態（定常運動状態、非定常運動状態、安静状態等）を判定する。表示部６００には、検出された脈拍数や拍動を示す波形、検出された運動状態、被検体の消費カロリー、現在の時刻等を表示することができる。
【００７０】
図１に示される拍動検出装置１００では、適応フィルター２０２の再構成が適切なタイミングで実行される。よって、適応フィルター２０２の継続的な使用によって、仮に、フィルターの性能が低下したとしても（例えば、ノイズ成分への追従が支配的になることによって、拍動信号に追従しないような場合が生じたとしても）、適応フィルター２０２は、再構成後に、新たに適応処理を開始することができ、この結果、拍動信号を捕捉するように自己適応がすすむ。よって、適応フィルター２０２は、拍動信号に追従する正常な状態に復帰することができる。
【００７１】
図１の拍動検出装置１００では、適応フィルター２０２を使用することから、複数のバンドパスフィルターを使用する必要がない。また、フィルターの再構成は、フィルター係数の初期化等によって実現でき、装置の処理負担の増加は少なく、また、装置の消費電力の増加も特に問題とならない。
【００７２】
よって、本実施形態によれば、例えば、装置の処理負担を低減しつつ、効果的なノイズ対策を実行できる拍動検出装置１００（例えば、長時間の使用に耐える拍動検出装置１００）を実現することができる。また、例えば、脈波センサー１０を手首の外側（腕時計の裏蓋面と接触する部位）など、脈波信号を取得しにくい部位に装着するタイプの拍動検出装置において、拍動呈示スペクトルの特定性能を向上させることが可能となる。
【００７３】
（腕時計型脈拍計の例）
図２（Ａ）〜図２（Ｃ）は、脈波センサーならびに拍動検出装置の構成の一例を示す図である。図２（Ａ）は、脈波センサー１０の一例の断面構造を示している。図２（Ａ）に示される脈波センサー１０は、基板３の裏面（下面）に設けられている発光素子１と、透明カバー（光透過性の材料で構成される接触部材）２と、光反射ドーム（光反射部）４と、基板３の表面（上面）に設けられる受光部５と、を有する。発光素子１から出射された光Ｒ１は、被検出部位ＳＡ（ここでは手首とする）における血管（生体情報源）Ｏに到達し、反射される。血管Ｏの容積は、拍動に伴って周期的に変動することから、反射光Ｒ１’の強度は、拍動に対応して周期的に変動する。反射光Ｒ１’は、反射ドーム４で反射された後、受光部５に入射する。受光部５は、入射光を電気信号に変換する。
【００７４】
図２（Ｂ）は、脈波センサー１０における回路構成の一例を示している。発光部１の発光は、制御回路１６１によって制御される。また、受光部５から出力される受光信号は、増幅回路１６２によって増幅された後、Ａ／Ｄ変換回路１６３によってデジタル信号に変換される。このようにして、脈波信号ｄが得られる。
【００７５】
図２（Ｃ）は、腕時計型の脈拍計（本実施形態の拍動検出装置を内蔵する）の使用例を示している。脈拍計は、リストバンド５００と、表示部６００と、本実施形態の拍動検出装置１００（内蔵されている）と、を有している。図２（Ｃ）の例では、脈拍計（すなわち拍動検出装置１００）は、被検体である人（ユーザー）の左手首に装着されている。
【００７６】
（脈波信号のきれいさの程度（ノイズ量の程度）の判定と、被検体の運動状態の判定について）
適応フィルター２０２の再構成の可否を判定する際に、脈波信号のきれいさの程度（ノイズ量の程度）の判定や、被検体の運動状態の判定等が必要となる場合がある。以下、これらの判定動作について、図３および図４を参照して説明する。
【００７７】
まず、脈波信号のきれいさの程度（ノイズ量の程度）の判定について説明する。図３（Ａ）〜図３（Ｃ）は、脈波信号のきれいさの程度（ノイズ量の程度）の判定動作を説明するための、ＦＦＴ前の信号の波形および周波数スペクトルの一例を示す図である。
【００７８】
図３（Ａ）〜図３（Ｃ）において、上側には、１６秒間のＦＦＴ前の脈波信号ｄの信号波形が示されている。横軸は時間を示し、縦軸は信号の振幅を示す。また、下側には、０から４Ｈｚの周波数帯域における周波数スペクトルが示されている。横軸は周波数を示し、縦軸はスペクトル値を示す。
【００７９】
ここで、図３（Ａ）は、ノイズ少（きれい）の場合における脈波信号ｄの波形と周波数スペクトルを示し、図３（Ｂ）は、ノイズが中程度（まあまあ）の場合における脈波信号の波形と周波数スペクトルを示しており、図３（Ｃ）は、脈波信号ｄに多くのノイズが含まれる場合（ノイジー）における脈波信号ｄの波形と周波数スペクトルを示している。図３（Ａ）〜図３（Ｃ）の各々の比較から明らかなように、脈波信号ｄの波形と周波数スペクトルとは密接に関連しており、脈波信号の波形に対応して、周波数スペクトルの分布状態やスペクトル値が変化する。よって、ＦＦＴによって得られる周波数スペクトルに基づいて、脈波信号ｄにする重畳するノイズの状態（ノイズ量の程度）を推定することが可能である。
【００８０】
本実施形態では、ノイズ量の程度の推定のための指標として、主要な周波数スペクトルのスペクトル値の比（つまり、基線の高さの比）を用いる。具体的には、ｒ５およびｒ１０という指標を用いる（ただし、一例であり、他の統計的指標、例えば、標準偏差等を用いてもよい）。ここで、ｒ５とは、１６秒分の脈波信号の周波数スペクトルの中から、ピーク値の大きさの順に５本のスペクトルを並べたとき（つまり、ソーティングしたとき）、第１番目のスペクトルのスペクトル値（パワー）を分母とし、第５番目のスペクトルのスペクトル値（パワー）を分子とすることによって得られる指標である。
【００８１】
また、ｒ１０とは、１６秒分の脈波信号の周波数スペクトルの中から、ピーク値の大きさの順に１０本のスペクトルを並べたとき（つまりソーティングしたとき）、第１番目のスペクトルのスペクトル値（パワー）を分母とし、第１０番目のスペクトルのスペクトル値（パワー）を分子とすることによって得られる指標である。
【００８２】
ここでは、一例として、ｒ５＜０．５かつｒ１０＜０．２のときをノイズ少（きれい）とし、ｒ５＞０．７かつｒ１０＜０．５のときをノイズ多（ノイジー）とし、上記いずれでもない場合をノイズが中程度（まあまあ）とする。
【００８３】
図３（Ａ）の例では、ｒ５＝０．１４かつｒ１０＝０．０８であることから、ノイズ少（きれい）と判定される。また、図３（Ｂ）の例では、ｒ５＝０．５６かつｒ１０＝０．３５であることから、ノイズが中程度（まあまあ）と判定される。図３（Ｃ）の例では、ｒ５＝０．８２かつｒ１０＝０．６２であることから、ノイズ多（ノイジー）と判定される。
【００８４】
このような判定処理が、図１に示される信号評価部２１４に含まれる信号のきれいさ判定部（ノイズ量の程度判定部）２１７によって実行される。
【００８５】
また、被検体の運動状態も、上述の指標ｒ５、ｒ１０を用いて推定することができる。被検体の運動状態によって脈波信号ｄの波形が変化すると、その変化は、周波数スペクトルの変化となって現れ、周波数スペクトルの変化は、指標ｒ５，ｒ１０に反映されるからである。
【００８６】
被検体の運動状態は、図１に示される信号評価部２１４に含まれる運動状態判定部２１９によって実行される。また、後処理を実行する脈波周波数解析部４００にも、運動状態判定部４１２が含まれており、この運動状態判定部４１２も同様の処理を実行して、被検体の運動状態を判定する。
【００８７】
図４は、被検体の運動状態の判定処理の一例を示すフローチャートである。まず、１６秒間に対応する２５６サンプルの脈波信号データを取得する（ステップＳＴ１）。次に、脈波信号データにＦＦＴ処理を施す（ステップＳＴ２）。
【００８８】
続いて、例えば、信号評価部２１４に含まれる、信号のきれいさ判定部（ノイズ量の程度判定部）２１７が、周波数スペクトルのソーティング処理を実行（ステップＳＴ３）する。ここで、ソーティング処理は、例えば、複数の周波数スペクトルを、スペクトル値が大きい順に並べる処理である。次に、指標ｒ５，ｒ１０を算出する（ステップＳＴ４）。
【００８９】
次に、信号のきれいさ判定部（ノイズ量の程度判定部）２１７は、ｒ５＞０．７かつｒ１０＜０．５であるか、すなわち、ノイズ多（ノイジー）の状態であるかを判定する（ステップＳＴ５）。そうであるならば（Ｙの場合）、非定常的なノイズが多い状態であるため、非定常的運動中と判定する（ステップＳＴ６）。
【００９０】
ステップＳＴ５において、Ｎの場合は、次に、信号のきれいさ判定部（ノイズ量の程度判定部）２１７は、ｒ５＜０．５かつｒ１０＜０．２であるか、すなわち、ノイズ少（きれい）の状態であるかを判定する（ステップＳＴ７）。そうであるならば（Ｙの場合）、次に、ステップＳＴ３で周波数スペクトルがソーティングされたスペクトルソート結果上位１０本以内に、体動（例えば、腕振りのような周期性のある体動）に相当する周波数スペクトルがあるかを判定する（ステップＳＴ８）。この判定に際して、体動センサー１１から得られる体動信号を参照して、体動スペクトルを特定することによって、判定の精度が向上する。なお、この場合には、図１において、破線で示される構成（体動センサー１１から出力されるＸ，Ｙの各方向の加速度信号を高速フーリエ変換部２１３，２１５に供給する構成）を使用することができる。
【００９１】
ステップＳＴ８にて、Ｙの場合は、ノイズ少（きれい）の状態で、かつ体動スペクトルが特定できる状態であるため、定常的運動中（例えば、同じ速度でウォーキング中）と判定する（ステップＳＴ９）。ステップＳＴ８にて、Ｎの場合は、ノイズ少（きれい）の状態で、かつ体動スペクトルが見あたらない状態であるため、安静状態と判定する（ステップＳＴ１０）。
【００９２】
また、ステップＳＴ７において、Ｎの場合は、次に、スペクトルソート結果上位１０本以内に、体動（例えば、腕振りのような周期性のある体動）に相当する周波数スペクトルがあるかを判定する（ステップＳＴ１１）。この判定に際して、体動センサー１１から得られる体動信号を参照して、体動スペクトルを特定することによって、判定の精度が向上する。なお、この場合には、図１において、破線で示される構成（体動センサー１１から出力されるＸ，Ｙの各方向の加速度信号を高速フーリエ変換部２１３，２１５に供給する構成）を使用することができる。
【００９３】
ステップＳＴ１１において、Ｙである場合には、ノイズ量が中程度（まあまあ）であり、かつ、体動スペクトルが特定できる状態であることから、定常運動中、但し非定常要素あり（例えば、同じ速度でウォーキング中であるが、腕時計を見るといった非定常動作が行われているような状態）と判定する（ステップＳＴ１２）。
【００９４】
また、ステップＳＴ１１において、Ｎである場合には、ノイズ量が中程度（まあまあ）であり、かつ、体動スペクトルが見あたらない状態であることから、安静中、但し非定常要素あり（例えば、ベッドに横たわっているが、寝返りをうつといった非定常的動作をしているような状態）と判定する（ステップＳＴ１３）。
【００９５】
このようにして、被検体の運動状態を推定することができる。但し、上述の例は一例であり、これに限定されるものではない。
【００９６】
（適応フィルターの再構成タイミングに関する時期的条件について）
図５は、適応フィルターの再構成が実行されるタイミングの一例を示す図である。例えば、図２に示されるような腕時計型の拍動検出装置は、被検体の手首に装着して使用され、例えば、被検体の経時的な動作状態の変化を記録するために、長時間にわたって継続的に使用される場合がある。よって、長時間にわたって継続的に使用しても、適応フィルター２０２の性能が低下しないようにすることが好ましい。よって、適応フィルター２０２が継続的に動作している期間中における適切なタイミングで、適応フィルター２０２を再構成する。
【００９７】
図５において、適応フィルターが動作を開始した時点（時刻ｔ１）から、第１時間（例えば１０分）が経過した第３時点（時刻ｔ２）までを、適応フィルター２０２の再構成が不要な第１期間ＴＡとし、第３時点（時刻ｔ２）から、第２時間（例えば５分）が経過した第４時点（時刻ｔ４）までを適応フィルターの再構成可能な第２期間ＴＣとする。なお、時刻ｔ１から時刻ｔ４までの期間ＴＢは、適応フィルター２０２を、再構成することなく使用可能な最大の期間である。
【００９８】
適応フィルター２０２が再構成されるタイミング（第１時点）は、第２期間ＴＣにおいて設定される。適応フィルター２０２の状態は、再構成の結果として、過去の第２時点の状態に戻るが、この第２時点は、第１期間ＴＡにおいて設定される。
【００９９】
適応フィルター再構成部１７は、時刻ｔ２からフィルター再構成の可否を、他の条件（後述）を考慮しつつ判定し、再構成可と判定すると、適応フィルター２０２の再構成処理を実行する。図５の例では、時刻ｔ３において、第１回目の適応フィルター２０２の再構成が実行されている。時刻ｔ３は、フィルター再構成の可否を判断するための所定の条件（判断要件）を満足する最初のタイミングである。適応フィルター２０２は、再構成後に、好ましい状態から、新たに適応処理を開始することができ、この結果、拍動信号を捕捉するように自己適応がすすむ。よって、適応フィルター２０２は、適切な性能を、継続的に維持することができる。
【０１００】
なお、適応フィルターの再構成の可否を、所定の条件（時期的条件以外の条件：図６を参照して後述する）によって判定した結果、所定条件が満足されず、その状態が、第２期間（適応フィルターの再構成可能期間）ＴＣにわたって継続する場合がある。この場合には、第２期間ＴＣの終点である第４時点（時刻ｔ４）を限界時点とし、この第４時点（時刻ｔ４）に達した時点で、所定条件の充足／非充足に関係なく、強制的に適応フィルター２０２を再構成することができる。このようにすれば、適応フィルター２０２の再構成ができない期間が長く続くというような、好ましくない状況が生じない。以上の動作が、その後、繰り返される。
【０１０１】
次に、適応フィルターが再構成された直近の時点（ここでは時刻ｔ３）から、第１期間ＴＡが開始される。時刻ｔ５から時刻ｔ７までが第２期間ＴＣとなる。なお、時刻ｔ６は、第２回目のフィルター再構成タイミングである。このような動作が、時刻ｔ７以降の、適応フィルター２０２の継続動作期間（ＴＤ）において、繰り返し実行される。適応フィルター２０２の動作は、時刻ｔ８に終了する。
【０１０２】
また、適応フィルター２０２の再構成タイミングである第１時点（図５の例では時刻ｔ３）は、適応フィルター２０２の動作開始時点（時刻ｔ０）を基準として、周期的（定期的）に設定することもできる。この場合、適応フィルターの再構成のタイミングは、例えば、タイミング判定部（タイマー）２１６（図１参照）によって簡単に管理することができ、実現が容易である。
【０１０３】
（フィルター再構成処理の具体的処理手順の例について）
図６は、適応フィルターの再構成処理の手順例を示すフローチャートである。まず、１６秒間の、２５６サンプルに相当する脈波信号データ、ならびに体動信号データが取得される（ステップＳＴ１０１）。次に、時期的条件（図５を用いて説明した条件）が満足されるか否かが判定される（ステップＳＴ１０２）。
【０１０４】
ステップＳＴ１０２において、時期的条件が満足される場合（ステップＳＴ１０２：Ｙ）には、次に、現在が、限界時点（図５における時刻ｔ４や時刻ｔ７）に達しているか否か、つまり、適応フィルター２０２の、累積使用時間が、第２期間ＴＣに相当する時間以上であるかが判定される（ステップＳＴ１０３）。限界時点に達している場合（ステップＳＴ１０３：Ｙ）には、適応フィルター再構成部１７は、適応フィルター２０２を再構成する（ステップＳＴ１０５）。
【０１０５】
ステップＳＴ１０３において、Ｎｏの場合には、ステップＳＴ１０４にすすみ、時期的条件以外の所定条件が満足されるか否かが判定される。所定条件は様々な条件が考えられるが、例えば、条件Ｉを満足する場合、条件ＩおよびIIを満足する場合、条件ＩおよびIIIを満足する場合、条件Ｉ，IIおよびIIIの全部を満足する場合等を考えることができる。ステップＳＴ１０４で所定条件を満足すれば（ステップＳＴ１０４：Ｙ）、適応フィルター再構成部１７は、適応フィルター２０２を再構成する（ステップＳＴ１０５）。以下、所定条件Ｉ〜IIIの内容と、各条件を組み合わせて判定する複数の例について説明する。
【０１０６】
（ステップＳＴ１０４における判定例Ａ：条件Ｉで判定）
適応フィルター再構成部１７は、第１時点前の直近の時点における、脈波周波数解析部４００（図１参照）による周波数解析によって、拍動呈示スペクトルが検出されていることを条件として（これが条件Ｉである）、適応フィルター２０２の再構成を実行する。
【０１０７】
適応フィルター２０２の再構成は、適応フィルター２０２の再構成が可能な第２期間ＴＣにおいて行われるが、適応フィルター２０２の再構成は第２期間ＴＣにおいて、いつでも行えるというものではなく、所定の条件を満たす場合に行うことができる。ここでは、脈波周波数解析部４００による前回の周波数解析によって、拍動呈示スペクトルが検出されていることが、適応フィルター再構成の条件（条件Ｉ）となる。
【０１０８】
すなわち、脈波周波数解析部４００は、脈波信号フィルタリング部から出力されるフィルタリング後の信号に基づいて、所定時間毎（例えば４秒毎）に周波数解析処理を行って拍動呈示スペクトルを特定する。適応フィルター再構成部１７は、第１時点においてフィルター再構成処理を実行しようとする場合、第１時点前の直近の時点における周波数解析（つまり、第１時点を基準とした場合における前回の周波数解析処理）において、拍動呈示スペクトル（拍動に対応する拍動基線の周波数スペクトル）の検出に成功しているかを確認し、その条件（条件Ｉ）を満たすときに、第１時点において適応フィルターの再構成を実行する。
【０１０９】
適応フィルター２０２を再構成すると、適応フィルター２０２は過去の第２時点の状態に戻って、新たに適応動作を開始する。この新たな動作の開始時点において、例えば、脈波信号ｄに混入しているノイズが少なく、脈波信号の波形のきれいさの程度が許容範囲である場合には、適応フィルターは、適切な自己適応が可能である。
【０１１０】
（ステップＳＴ１０４における判定例Ｂ：条件Ｉおよび条件IIで判定）
一方、脈波信号ｄに混入しているノイズが多く、脈波信号の波形のきれいさの程度が、許容範囲外である場合には、適応フィルター２０２はノイズに追従する可能性が高い。適応フィルター２０２の再構成は、その時点で捕捉している自己相関性が高いと考えられる信号への追従を断つ処理でもあるため、再構成の直後における拍動呈示スペクトルの捕捉に失敗すると、例えば、第１時点の直前まで検出できていた拍動信号成分をカットしてしまうという結果となり、かえって悪い結果となる場合がある。
【０１１１】
適応フィルター再構成部１７は、第１時点前の直近の時点における、脈波周波数解析部４００による周波数解析によって、拍動呈示スペクトルが検出されており（条件Ｉを満足）、かつ、信号評価部２１４によって、ノイズ量が所定の基準以下であると判定されたこと（条件II）が満足される場合に、適応フィルター２０２を再構成する。ノイズ量の判定にあたっては、図３（Ａ）〜図３（Ｃ）を用いて説明した判定方法を利用することができる。
【０１１２】
つまり、この判定例Ｂでは、条件Ｉに加えて、信号評価部２１４によって、脈波信号ｄに混入されるノイズ量（外乱ノイズ成分のノイズ量）が所定の基準以下であると判定されること（条件II）も満足することが要求される。信号評価部２１４によるノイズ量の判定によって、脈波信号ｄのきれいさの程度を、より正確に推定することができる。
【０１１３】
この例では、外乱ノイズの少ない状態のときにフィルターを再構成することによって、再構成後の適応フィルターによって適切な自己適応が実行され、拍動成分をノイズ成分と区別することができる。よって、適応フィルターの再構成の可否を、より的確に判断することができる。
【０１１４】
（ステップＳＴ１０４における判定例Ｃ：条件Ｉおよび条件IIIで判定）
この例では、適応フィルター再構成部１７は、第１時点前の直近の時点における、脈波周波数解析部４００による周波数解析によって、拍動呈示スペクトルが検出されており（条件Ｉが満足される）、かつ、運動状態判定部４１２（図１参照）によって、被検体の運動状態が安定状態である場合（つまり、条件IIIも満足される場合）に、適応フィルター２０２を再構成する。
【０１１５】
被検体（人や動物）の運動状態は、絶えず変化する可能性がある。定常的な運動状態では、拍動信号も安定しているが、運動状態が急に変化したとき（運動状態の変化の過渡期）においては、拍動信号の振幅や周波数が不安定化し、例えば突発的に大きく変化する場合もあり得る。このような場合に、適応フィルターの再構成を実行すると、拍動信号成分ではなく、ノイズ成分に追従する可能性が高くなる。
【０１１６】
よって、判定例Ｃでは、上述の条件Ｉに、さらに、被検体の運動状態が安定状態であることが確認されるという条件（条件III）を追加（加重）する。よって、適応フィルター２０２の再構成の可否を、より的確に判断することができる。被検体の運動状態の判定にあたっては、図４を用いて説明した判定方法を利用することができる。
【０１１７】
例えば、被検体である人の運動状態が、歩行中、ジョギング中などの定常的な運動をしている状態、もしくは安静状態のいずれかに該当しない状態である場合（例えば球技をしている、体操している状態等の場合）には、適応フィルター２０２の再構成が実行されない。よって、拍動検出の失敗や誤検出の発生を低減することができる。
【０１１８】
（ステップＳＴ１０４における判定例Ｄ：条件Ｉ、条件IIおよび条件IIIで判定）
本例では、適応フィルター再構成部１７は、第１時点前の直近の時点における、脈波周波数解析部による周波数解析によって、拍動呈示スペクトルが検出されており（条件Ｉを満足する）、信号評価部２１４によって、ノイズ量が所定の基準以下であると判定され（条件IIを満足する）、かつ、被検体の運動状態が、定常的運動状態および安静状態を含む安定状態である場合（つまり、条件IIIを満足する場合）に、適応フィルター２０２を再構成する。
【０１１９】
判定例Ｄでは、適応フィルター２０２の再構成の条件として、条件Ｉ，IIならびにIIIをすべて満たすか否かを判定する。よって、適応フィルター２０２の再構成の可否を、さらに的確に判断することができる。
【０１２０】
次に、適応フィルタリング処理（周波数解析前において、脈波信号フィルタリング部２００によって実行される前処理）と、脈波周波数解析部４００によって行われる周波数解析処理（後処理）の手順について説明する。
【０１２１】
図７は、適応フィルタリング処理（前処理）ならびに周波数解析処理（後処理）の手順の一例を示すフローチャートである。前処理では、まず、遅延処理部２０１にて、脈波信号ｄを遅延させる（ステップＳＴ１０６）。次に、処理したデータ数が、全サンプル数に到達しているかを判定し（ステップＳＴ１０７）、到達していない場合（ステップＳＴ１０７：Ｙ）には、適応フィルター２０２による処理を行って、自己相関性が高い信号である第１信号（定常的成分）ｙを出力する（ステップＳＴ１０８）。続いて、脈波信号ｄから第１信号ｙを減算することによって、自己相関性が低い信号ｅを出力する（ステップＳＴ１０９）。
【０１２２】
次に、サンプリング区間内の脈波信号に過大振幅が検出されたか否かを判定する（ステップＳＴ１１０）。過大振幅が検出された場合（ステップＳＴ１１０：Ｙ）には、衝撃ノイズが加わった可能性が高いことから、衝撃検出モードのゲイン係数（例えば、ｈ１＝１．２，ｈ２＝０）が選択され、選択されたゲイン係数に基づいて出力値が算出される（ステップＳＴ１１１）。
【０１２３】
ステップＳＴ１１０において、過大振幅が検出されない場合（ステップＳＴ１１０：Ｎ）には、最後に過大振幅が検出されてから、一定期間、過大振幅が検出されていないかを判定し（ステップＳＴ１１２）、検出されていないならば（ステップＳＴ１１２：Ｙ）、適応フィルター係数更新処理を実行し（ステップＳＴ１１３）、次に、通常モードのゲイン係数（例えば、ｈ１＝１．０，ｈ２＝０．５）が選択され、選択されたゲイン係数に基づいて出力値が算出される（ステップＳＴ１１４）。
【０１２４】
次に、後処理の手順について説明する。まず、フィルタリング後の信号と、加速度信号等（体動信号）とが取得される（ステップＳＴ１１５）。次に、体動ノイズ除去処理部３００（図１参照）による、体動除去適応フィルタリング処理が実行される（ステップＳＴ１１６）。次に、脈波周波数解析部４００によって、周波数解析処理が実行される（ステップＳＴ１１７）。次に、脈拍数算出部４１０によって、脈拍数が算出される（ステップＳＴ１１８）。
【０１２５】
（実際の処理例）
図８（Ａ）〜図８（Ｃ）は、実際の処理における脈波信号波形と、ＦＦＴによって得られる周波数スペクトルの一例を示す図である。図８では、１６秒分の脈波信号の波形（上段）と、そのＦＦＴ結果（下段）を示している。なお、この例で取り扱っている脈波信号は、歩行やジョギング等の周期的運動を行っていないときに得られた脈波信号である。
【０１２６】
図８（Ａ）の上段の波形では、１６秒分の信号の終端部分付近において、外乱ノイズ成分と思われる信号の混入が見られる（点線で囲んで示される部分Ａ）。下段のＦＦＴ結果においては、脈波信号には、検出対象成分である拍動信号成分ＭＰ１ａと、外乱ノイズ信号成分ＳＰ１ａとが混在している。
【０１２７】
次に、図８（Ｂ）を参照する。図８（Ｂ）は、図８（Ａ）の上段に示される波形をもつ脈波信号を、フィルター再構成処理が実施されない適応フィルターにてフィルタリングして得られる脈波信号と、そのＦＦＴ結果を示している。ここでは、適応フィルターは、計測開始から４０分程度、継続的に動作している状態（フィルター係数を適応的に更新しながら使用してきた状態）となっている（但し、フィルター再構成処理は実施されていない）。
【０１２８】
図８（Ｂ）から明らかなように、外乱ノイズ成分（ノイズ信号成分）ＳＰ１ｂは、カットしきれずに残存している。たまたま、図８（Ｂ）の例では、拍動呈示スペクトルＭＰ１ｂのスペクトル値が最も大きいため、各スペクトルのピーク値に基づくソーティング処理等によって、信号ＭＰ１ｂを、拍動呈示スペクトルであると判定することは可能である。但し、拍動呈示スペクトルＭＰ１ｂが小さい場合には、外乱ノイズ成分（ノイズ信号成分）ＳＰ１ｂと区別することが困難となり、拍動成分の検出が失敗に終わる可能性が高くなる。
【０１２９】
次に、図８（Ｃ）を参照する。図８（Ｃ）は、図８（Ａ）の上段に示される波形をもつ脈波信号を、適応フィルター２０２（フィルター再構成処理を例えば１０分間隔で実施した場合）にてフィルタリングして得られる脈波信号と、そのＦＦＴ結果を示している。ここでは、適応フィルター２０２は、計測開始から４０分程度、継続的に動作している状態（フィルター係数を適応的に更新しながら使用してきた状態）となっている（フィルター再構成処理は、例えば１０分間隔で実施されたものとする）。
【０１３０】
図８（Ｃ）から明らかなように、外乱ノイズ成分ＳＰ１ｃは、適応フィルタリング処理によって十分に抑制され、その一方、拍動呈示スペクトルＭＰ１ｂは強調されており、よって、拍動成分を、確実に特定することができる。
【０１３１】
フィルタの性能は、Ｓ／Ｎ指標を利用して客観的に評価することが可能である。以下、Ｓ／Ｎ指標に基づく評価によって、本実施形態の効果を実証する。なお、以下の説明では、Ｓ／Ｎ指標として、ＳＮ３と呼ばれる評価値を利用する。ここで、ＳＮ３は、以下の算出式によって表すことができる。
ＳＮ３＝（拍動呈示スペクトルとその左右１本のスペクトル値の合計）／（全周波数０〜４Ｈｚにおけるスペクトル値の合計）（単位：％）
図８（Ａ）に示される脈波信号（原信号）におけるＳ／Ｎ指標（ＳＮ３）は、１３．８％であった。図８（Ｂ）に示される脈波信号（適応フィルターの再構成無し）におけるＳ／Ｎ指標（ＳＮ３）は、１３．３％であり、原信号に比べて、やや低下している。これに対して、図８（Ｃ）に示される脈波信号（適応フィルターの再構成有り）におけるＳ／Ｎ指標（ＳＮ３）は、１８．７％であり、Ｓ／Ｎは大幅に改善されている。このＳ／Ｎ指標の改善の事実によって、適応フィルター２０２を再構成することによる効果が、数値的に確認できた。
【０１３２】
（第２実施形態）
図９（Ａ）および図９（Ｂ）は、適応フィルターの再構成処理を実行する例と、しない場合の例（比較例）の各々におけるＳ／Ｎ指標（ＳＮ３）の、時間経過に対応した変動を調べた結果を示す図である。
【０１３３】
ここでは、被検体である人が腕時計型の脈拍計（すなわち拍動検出装置１００）を装着し、被検体が２０分間にわたって歩行し、脈波信号を適応フィルターによってフィルタリングし、そして、フィルタリング後の信号について、Ｓ／Ｎ指標（ＳＮ３）を算出し、その算出結果を、時間軸上でプロットした。図９（Ａ）および図９（Ｂ）において、横軸が時間軸であり、縦軸がＳ／Ｎ指標であるＳＮ３の値（％）を示す軸である。図９（Ａ）は、被検体であるユーザーが歩行を開始してから、１０分経過するまでのＳＮ３の経時的変化が示されている。また、図９（Ｂ）は、１５分経過時点から、２０分経過するまでのＳＮ３の経時的変化が示されている。また、本実施形態に対応するデータは、黒塗りの四角形で示され、比較例のデータは、白抜きの菱形で示されている。
【０１３４】
計測開始から１０分経過頃の時点まで（図９（Ａ））においては、本実施形態と比較例との間には、ＳＮ３の値には差異は、ほとんど現れない。しかしながら、図９（Ｂ）から明らかなように、１５分経過頃までには差異がはっきりと現れる。
【０１３５】
図９（Ａ）および図９（Ｂ）によって、脈波信号フィルタリング部（適応線スペクトル強調器）を構成する適応フィルター２０２を、適切なタイミングで再構成することによって、適応フィルター２０２の性能を適切なレベルに維持することができることが裏付けされた。よって、長時間の連続計測を行う場合であっても、適応フィルター２０２の性能を、適切なレベルに維持することができ、拍動検出の失敗や誤検出の発生を低減することができる。
【０１３６】
以上説明したように、本発明の少なくとも一つの実施形態によれば、脈波信号フィルタリング部を構成する適応フィルターを適切なタイミングで再構成することによって、適応フィルターの性能を適切なレベルに維持することができる。よって、例えば、長時間の連続計測を行う場合であっても、適応フィルターの性能を、適切なレベルに維持することができる。これにより拍動検出の失敗や誤検出の発生を低減することができる。
【０１３７】
また、例えば、適応フィルターを有効活用することから、フィルター構成を簡素化でき、また、適応フィルターの再構成を行うことは容易であり、よって、装置の処理負担を低減しつつ、効果的なノイズ対策を実行できる拍動検出装置を実現することができる。
【０１３８】
なお、上記のように本実施形態について詳細に説明したが、本発明の新規事項及び効果から実体的に逸脱しない多くの変形が可能であることは当業者には容易に理解できるであろう。従って、このような変形例はすべて本発明の範囲に含まれるものとする。例えば、明細書又は図面において、少なくとも一度、より広義又は同義な異なる用語と共に記載された用語は、明細書又は図面のいかなる箇所においても、その異なる用語に置き換えることができる。
【符号の説明】
【０１３９】
１０  脈波センサー、１１  体動センサー（加速度センサーやジャイロセンサー等）、
１２  脈波信号蓄積部、１７  フィルター再構成部、
２００  脈波信号フィルタリング部、
３００  体動ノイズ除去処理部、４００  脈波周波数解析部、６００  表示部
【図１】
0005640819.tif 000002
【図２】
0005640819.tif 000003
【図３】
0005640819.tif 000004
【図４】
0005640819.tif 000005
【図５】
0005640819.tif 000006
【図６】
0005640819.tif 000007
【図７】
0005640819.tif 000008
【図８】
0005640819.tif 000009
【図９】
0005640819.tif 000010