JP5333427B2

JP5333427B2 - 心拍数検出装置、心拍数検出方法、およびプログラム

Info

Publication number: JP5333427B2
Application number: JP2010282946A
Authority: JP
Inventors: 理敏関根
Original assignee: Oki Electric Industry Co Ltd
Current assignee: Oki Electric Industry Co Ltd
Priority date: 2010-12-20
Filing date: 2010-12-20
Publication date: 2013-11-06
Anticipated expiration: 2030-12-20
Also published as: JP2012130391A

Description

本発明は、心拍数検出装置、心拍数検出方法、およびプログラムに関し、特にドップラーセンサの出力信号に基づいて心拍数を検出する心拍数検出装置、心拍数検出方法、およびプログラムに関する。

近年、ドップラーセンサを用いて心拍数を検知するシステムに関する研究が進められている。ドップラーセンサは、検知対象に対して放射する電磁波又は超音波である放射波の周波数と、放射波が検知対象により反射した反射波の周波数との差分の周波数を有するドップラーセンサ出力信号を出力する。従って、ドップラーセンサを用いた心拍数検出装置は、被験者に接触することなく心拍数を検出することができるという利点がある。例えば特許文献１には、このドップラーセンサ出力信号のうち、心拍数の周波数帯域（０．８〜２．５Ｈｚ）を通過帯域とするフィルタを通した波形のうち最大振幅の周波数から心拍数を検出する装置が開示されている。

特開２０００−８３９２７号公報

しかし、ドップラーセンサ出力信号には、心拍による影響のみならず、呼吸による身体の動き、身体の揺れ、その他の環境雑音などが含まれている。このため、心拍数を正確に抽出するのが困難であるという問題があった。
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、より精度の高い心拍数を検出することが可能な、新規かつ改良された心拍数検出装置、心拍数検出方法、およびプログラムを提供することにある。

上記課題を解決するために、本発明のある観点によれば、被験者に対して放射した電磁波又は超音波である放射波の周波数と、上記放射波が上記被験者により反射した反射波の周波数との差分の周波数を有するドップラーセンサ出力信号を取得するデータ取得部と、上記ドップラーセンサ出力信号を複数の周波数における成分信号に分解し、複数の上記成分信号についてそれぞれ信号強度の時間変動を解析する時間周波数解析部と、複数の上記成分信号についてそれぞれ信号強度の自己相関を解析する自己相関解析部と、上記複数の成分信号の中から自己相関関数に基づいて周期性の最も高い成分信号を選択し、選択された成分信号の上記自己相関の解析結果に基づいて上記被験者の心拍数を検出する心拍数検出部と、を有することを特徴とする、心拍数検出装置が提供される。

かかる構成によれば、ドップラーセンサ出力信号を複数の成分信号に分け、この成分信号のうち周期性の高い成分信号に基づいて心拍数を検出することができる。これにより、分解しないそのままのドップラーセンサ出力信号の自己相関を解析する場合と比較して精度の高い心拍数を検出することができる。

また、上記ドップラーセンサ出力信号から、心拍に起因した上記被験者の身体の動きの速度変化に応じた周波数成分を含む信号を抽出するフィルタ部、をさらに有してもよい。

また、上記被験者の身体の動きの位置変化に応じた周波数成分を含む信号を除去するフィルタ部、をさらに有してもよい。

また、上記心拍数検出部は、自己相関関数のピークのうち最も振幅の大きい第一ピークの値と２番目に振幅の大きい第二ピークの値との比に基づいて、周期性が最も高い成分信号を選択してもよい。

また、上記心拍数検出部は、上記第一ピークの時間を心拍間隔として上記心拍数を検出してもよい。

また、上記課題を解決するために、本発明の別の観点によれば、被験者に対して電磁波又は超音波である放射波を放射するステップと、上記放射波が上記被験者により反射した反射波を検出するステップと、上記放射波の周波数と上記反射波の周波数との差分の周波数を有するドップラーセンサ出力信号を取得するステップと、上記ドップラーセンサ出力信号を複数の周波数における成分信号に分解するステップと、複数の上記成分信号についてそれぞれ信号強度の時間変動を解析するステップと、複数の上記成分信号についてそれぞれ上記信号強度の自己相関を解析するステップと、複数の上記成分信号の中から最も自己相関関数の周期性が高い成分信号を選択するステップと、上記選択された成分信号の上記自己相関の解析結果に基づいて、上記被験者の心拍数を検出するステップと、を含む、心拍数検出方法が提供される。

また、上記課題を解決するために、本発明の別の観点によれば、コンピュータを、被験者に対して放射した電磁波又は超音波である放射波の周波数と、上記放射波が上記被験者により反射した反射波の周波数との差分の周波数を有するドップラーセンサ出力信号を取得するデータ取得部と、上記ドップラーセンサ出力信号を複数の周波数における成分信号に分解し、複数の上記成分信号についてそれぞれ信号強度の時間変動を解析する時間周波数解析部と、複数の上記成分信号についてそれぞれ信号強度の自己相関を解析する自己相関解析部と、上記複数の成分信号の中から最も自己相関関数の周期性が高い成分信号を選択し、選択された成分信号の上記自己相関の解析結果に基づいて上記被験者の心拍数を検出する心拍数検出部とを有することを特徴とする、心拍数検出装置として機能させるためのプログラムが提供される。

以上説明したように本発明によれば、より精度の高い心拍数を検出することができる。

本発明の一実施形態に係る心拍数検出装置の概要について説明するための説明図である。同実施形態に係る心拍数検出装置の構成を示すブロック図である。同実施形態に係る心拍数検出装置の心拍数抽出部１０６の詳細な構成を示すブロック図である。同実施形態に係る心拍数検出装置の動作の一例を示すフローチャートである。ドップラーセンサ出力信号から心拍に依存した低周波数成分を含む成分を抽出した信号の一例を示すグラフである。ドップラーセンサ出力信号から心拍に依存した高周波数成分を抽出した信号の一例を示すグラフである。５Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。１０Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。１５Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。２０Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。図７のウェーブレット係数の自己相関関数を示すグラフである。図８のウェーブレット係数の自己相関関数を示すグラフである。図９のウェーブレット係数の自己相関関数を示すグラフである。図１０のウェーブレット係数の自己相関関数を示すグラフである。周期性を評価するための評価関数について説明するための説明図である。同実施形態による精度向上の効果を説明するための説明図である。

以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。

なお、説明は以下の順序で行うものとする。
１．概要
２．構成
３．動作例
４．効果の例

＜１．概要＞
まず初めに、本発明の一実施形態に係る心拍数検出装置の概要について図１を参照しながら説明する。図１は、本発明の一実施形態に係る心拍検出装置の概要について説明するための説明図である。

図１に示されるように、本実施形態に係る心拍数検出装置１００は、被験者Ｍの胸部付近に放射波Ｔｘを放射し、放射波Ｔｘが被験者Ｍに反射した反射波Ｒｘを取得するドップラーセンサを有する。このドップラーセンサは、放射波Ｔｘの周波数と反射波Ｒｘの周波数との差分の周波数を有するドップラーセンサ出力信号を出力する。心拍数検出装置１００は、このドップラーセンサ出力信号の解析結果に基づいて被験者Ｍの心拍数を検出する。

心拍数検出装置１００は、電磁波又は超音波である放射波Ｔｘを用いるため、被験者Ｍに接触することなく心拍数を検出することができる。このため、心拍数検出装置１００は、ヘルスケアの分野における利用が期待されている。例えば心拍数検出装置１００は、高齢者などの健康状態を遠隔監視するシステムなどに用いることができる。

なお、ドップラーセンサは、対象物（ここでは被験者Ｍ）の移動速度に比例して反射波Ｒｘの周波数がシフトする現象（ドップラー効果）を利用したセンサである。しかし、ここで取得されるドップラーセンサ出力信号には、検出対象である心拍による影響だけでなく、例えば呼吸による身体の動き、身体の揺れ、およびその他の環境雑音等が含まれる。精度の高い心拍数を検出するためには、これらの心拍以外の要因によるドップラーセンサ出力信号への影響を取り除き、心拍による周波数シフトを検出することが重要である。

被験者Ｍが動くとき、ドップラーセンサから観測されるドップラーセンサ出力信号Ｖ_ｄは、以下の式で表される。

なお、振幅Ｋ_ｄは、被験者Ｍの身体の動きの変位に依存する。また、周波数ｆ_ｄは、被験者Ｍの身体の動きの速度に依存する。１分間の心拍数が６０である場合には、１秒の周期で被験者Ｍの身体の動きが生じる。このため、ドップラーセンサ出力信号にもこの動きに起因する振幅Ｋ_ｄの変化に依存した成分が含まれている。ところが、この１秒前後の周期を持つ信号には、呼吸による影響や身体の揺れによる影響が含まれやすい。このため、正確に心拍数を検出することが困難であった。

そこで、本実施形態に係る心拍数検出装置１００は、心拍に起因した身体の動きの変位に応じた周波数成分ではなく、身体の動きによる速度の大小、つまりｆ_ｄの大小変化に依存した周波数成分を観測する。このとき、ｆ_ｄの大小変化に依存した成分を含む信号の周波数は、個人差がある。また、同一人物であっても放射波Ｔｘの放射される位置や放射波Ｔｘに対する身体の向きが異なると、ｆ_ｄの大小変化に依存した成分を含む信号の周波数、およびその信号の振幅の大きさは異なる。このため、心拍数検出装置１００は、ドップラーセンサ出力信号を複数の周波数成分に分解して得られる複数の成分信号から最も周期性の高い成分信号を選択し、選択された成分信号を用いて心拍数の検出を行う。

かかる心拍数検出装置１００の機能を実現するための機能構成と、動作の一例について次に説明される。

＜２．構成＞
ここで、図２および図３を参照しながら、本実施形態に係る心拍数検出装置１００の機能構成の一例について説明する。図２は、本実施形態に係る心拍数検出装置の構成を示すブロック図である。図３は、本実施形態に係る心拍数検出装置の心拍数抽出部１０６の詳細な構成を示すブロック図である。

心拍数検出装置１００は、ドップラーセンサ１０１と、増幅器１０２と、アナログフィルタ１０３と、Ａ／Ｄ（Ａｎａｌｏｇ／Ｄｉｇｉｔａｌ）変換器１０４と、デジタルフィルタ１０５と、心拍数抽出部１０６と、心拍数表示部１０７とを主に有する。

ドップラーセンサ１０１は、上述の通り、対象物（ここでは被験者Ｍ）の移動速度に比例して反射波Ｒｘの周波数がシフトする現象（ドップラー効果）を利用したセンサである。ドップラーセンサ１０１は、放射波Ｔｘが被験者Ｍに反射した反射波Ｒｘを検出する。そしてドップラーセンサ１０１は、放射波Ｔｘの周波数ｆ_ｔと反射波Ｒｘの周波数ｆ_ｒとの差分の周波数ｆ_ｄ＝ｆ_ｔ−ｆ_ｒを有するドップラーセンサ出力信号を出力する。このドップラーセンサ出力信号の電圧Ｖｄは、既出の通り下記の数式１で示される。

ここで、振幅Ｋ_ｄは、電波受信強度に依存する値となり、電波受信強度Ｐ_ｒは、一般的に以下の数式２で表されるレーダ方程式に従う。ここで、Ｐ_ｔは送信電力、Ｇ_ｔは送信アンテナ利得、σは当該場所における散乱断面積、Ａ_ｒは受信アンテナの有効開口面積、Ｒはセンサと認識対象間の距離である。

なお、ここで用いられる放射波Ｔｘは、電磁波又は超音波などであってよい。ドップラーセンサ１０１は、一般的なドップラーセンサモジュールが用いられてよい。そして、ドップラーセンサ１０１は、取得したドップラーセンサ出力信号を増幅器１０２に入力する。

増幅器１０２は、ドップラーセンサ１０１から入力されたドップラーセンサ出力信号を増幅する機能を有する。増幅器１０２は、一般的な増幅回路が用いられてよい。このとき増幅器１０２は、アナログフィルタ１０３およびデジタルフィルタ１０５を通した後のドップラーセンサ出力信号の振幅が観察に適した振幅となるように増幅率が決められる。増幅器１０２は、増幅処理を施した後のドップラーセンサ出力信号をアナログフィルタ１０３に入力する。

アナログフィルタ１０３は、増幅器１０２から入力されたドップラーセンサ出力信号から雑音等による影響の大きい周波数成分を除去し、心拍の影響を含む周波数成分を抽出するフィルタである。なお、雑音等の影響が大きい周波数成分を除去することにより、心拍による動きの変位に起因する周波数成分も除去される。アナログフィルタ１０３は、例えばローパスフィルタ、ハイパスフィルタ、またはバンドパスフィルタ等のいずれかであってよい。或いはアナログフィルタ１０３は、ローパスフィルタ、ハイパスフィルタ、またはバンドパスフィルタ等が組み合わせて用いられてもよい。アナログフィルタ１０３は、フィルタリング処理後のドップラーセンサ出力信号をＡ／Ｄ変換器１０４に入力する。

Ａ／Ｄ変換器１０４は、アナログフィルタ１０３から入力されたアナログ信号であるドップラーセンサ出力信号をデジタル信号に変換する機能を有する。Ａ／Ｄ変換器１０４は、デジタル信号に変換されたドップラーセンサ出力信号をデジタルフィルタ１０５に入力する。

デジタルフィルタ１０５は、Ａ／Ｄ変換器１０４によりデジタル化されたドップラーセンサ出力信号にフィルタリング処理を行うフィルタである。デジタルフィルタ１０５は、アナログフィルタ１０３では除ききれなかった周波数成分を除去する機能を有する。デジタルフィルタ１０５は、フィルタリング処理後のドップラーセンサ出力信号を心拍数抽出部１０６に入力する。

心拍数抽出部１０６は、入力されたドップラーセンサ出力信号に基づいて心拍数を抽出する機能を有する。図３を参照すると、心拍数抽出部１０６は、ドップラーセンサ出力信号を取得するデータ取得部１０６１と、取得したドップラーセンサ出力信号をウェーブレット変換により複数の周波数における成分信号に分解し、成分信号の信号強度の時間変動を解析するウェーブレット解析部１０６２と、ウェーブレット係数の自己相関を解析する自己相関解析部１０６３と、自己相関の解析結果に基づいて心拍数を検出する心拍数検出部１０６４として主に機能する。

ここでウェーブレット解析部１０６２は、データ取得部が取得したドップラーセンサ出力信号を複数の周波数における成分信号に分解し、この成分信号について信号強度の時間変動を解析する時間周波数解析部の一例である。ウェーブレット解析部１０６２はウェーブレット変換を用いて、ドップラーセンサ出力信号を複数の擬似周波数における成分信号に分解する。本実施形態においては、５Ｈｚ、１０Ｈｚ、１５Ｈｚ、２０Ｈｚの擬似周波数についてウェーブレット解析が実施される。

また、自己相関解析部１０６３は、ウェーブレット変換後の各成分信号の自己相関関数からピークを抽出し、抽出したピークに基づいて複数の成分信号から最も周期性の高い成分信号を選択する。心拍数検出部１０６４は、自己相関解析部により選択された成分信号を用いて心拍数を検出する。

心拍数表示部１０７は、心拍数抽出部１０６により抽出された心拍数の情報を表示する機能を有する。心拍数表示部１０７は、例えば心拍数抽出部１０６により抽出された心拍数の情報を含む表示画面を生成する表示画面生成部と、生成された表示画面を表示する表示装置とを含んでもよい。表示装置は、例えば液晶ディスプレイ（ＬＣＤ）装置、ＯＬＥＤ（ＯｒｇａｎｉｃＬｉｇｈｔＥｍｉｔｔｉｎｇＤｉｏｄｅ）装置、ＣＲＴ（ＣａｔｈｏｄｅＲａｙＴｕｂｅ）ディスプレイ装置、およびランプなどの表示装置であってよい。表示装置がランプである場合には、例えば心拍数抽出部１０６が抽出した心拍数に基づいて、心拍と同じ間隔でランプを点滅させてもよい。

以上、本実施形態に係る心拍数検出装置１００の機能の一例を示した。上記の各構成要素は、汎用的な部材や回路を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。また、各構成要素の機能を、ＣＰＵ（ＣｅｎｔｒａｌＰｒｏｃｅｓｓｉｎｇＵｎｉｔ）などの演算装置がこれらの機能を実現する処理手順を記述した制御プログラムを記憶したＲＯＭ（ＲｅａｄＯｎｌｙＭｅｍｏｒｙ）やＲＡＭ（ＲａｎｄｏｍＡｃｃｅｓｓＭｅｍｏｒｙ）などの記憶媒体から制御プログラムを読出し、そのプログラムを解釈して実行することにより行ってもよい。従って、本実施形態を実施する時々の技術レベルに応じて、適宜、利用する構成を変更することが可能である。

なお、上述のような本実施形態に係る心拍数検出装置１００の各機能を実現するためのコンピュータプログラムを作製し、パーソナルコンピュータ等に実装することが可能である。また、このようなコンピュータプログラムが格納された、コンピュータで読み取り可能な記録媒体も提供することができる。記録媒体は、例えば、磁気ディスク、光ディスク、光磁気ディスク、フラッシュメモリなどである。また、上記のコンピュータプログラムは、記録媒体を用いずに、例えばネットワークを介して配信してもよい。

＜３．動作＞
次に、図４〜図１４を参照しながら、本実施形態に係る心拍数検出装置１００の動作について説明する。図４は、本実施形態に係る心拍数検出装置の動作の一例を示すフローチャートである。図５は、ドップラーセンサ出力信号の一例を示すグラフである。図６は、図５のドップラーセンサ出力信号から高周波数成分を抽出した信号の一例を示すグラフである。図７は、５Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。図８は、１０Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。図９は、１５Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。図１０は、２０Ｈｚにおけるウェーブレット係数の時間変化を示すグラフである。図１１は、図７のウェーブレット係数の自己相関関数を示すグラフである。図１２は、図８のウェーブレット係数の自己相関関数を示すグラフである。図１３は、図９のウェーブレット係数の自己相関関数を示すグラフである。図１４は、図１０のウェーブレット係数の自己相関関数を示すグラフである。

まず図４を参照すると、本実施形態に係る心拍数検出装置１００の動作の一例が示されている。ここに示される動作は、主に１．ドップラーセンサ出力信号を取得し、心拍の影響を含む高周波成分を抽出する第１段階（Ｓ１０１〜Ｓ１０５）、２．複数の周波数における成分信号を抽出する第２段階（Ｓ１０７）、３．成分信号の自己相関関数に基づいて心拍数を検出する第３段階（Ｓ１０９〜Ｓ１１５）、の３つに区分される。

まず、第１段階について説明する。心拍数検出装置１００のドップラーセンサから放射される放射波Ｔｘの放射方向が被験者Ｍの方向となるように心拍数検出装置１００の設置位置または被験者Ｍの位置を調整する。そして、心拍数検出装置１００は、被験者Ｍの胸部付近に放射波Ｔｘを放射する（Ｓ１０１）。このとき、放射波Ｔｘの放射位置は、検出したい動きの影響を観測しやすい部位とすることが好ましい。本実施形態においては、心拍の影響は胸部付近において観測しやすいため、胸部に放射波Ｔｘを放射する。また、観測時において被験者Ｍの身体の動きは少ない方が心拍数の検出精度は高まる。

そして、ドップラーセンサ１０１は、放射波Ｔｘと、放射波Ｔｘが被験者Ｍにより反射した反射波Ｒｘとに基づいて、放射波Ｔｘの周波数ｆｔと反射波Ｒｘの周波数ｆｒとの差分の周波数ｆｄを有するドップラーセンサ出力信号を出力する（Ｓ１０３）。このとき出力されるドップラーセンサ出力信号は、例えば図５に示されるように、呼吸による信号強度の変動と心拍による信号強度の変動とを含む。図５は、ドップラーセンサ１０１が出力したドップラーセンサ出力信号を増幅器１０２で増幅した後、電源ノイズをカットした後の波形を示したものである。

次に、心拍数検出装置１００は、増幅器１０２、アナログフィルタ１０３、Ａ／Ｄ変換器１０４、およびデジタルフィルタ１０５を用いて、取得されたドップラーセンサ出力信号に各種の信号処理を行う（Ｓ１０５）。ここで行われるフィルタリング処理は、例えば環境雑音、身体の揺れ、および呼吸などの心拍による影響とは無関係な成分を除き、心拍による身体の動きの速度の変化による周波数成分を抽出するために行われる。なお、本実施形態においては、心拍による身体の位置変化に応じた周波数成分もドップラーセンサ出力信号から除く。心拍による身体の位置変化に応じた周波数成分は、概ね１秒前後の心拍間隔の周期を有する。かかる周波数成分は、身体のゆれなどによって生じる周波数成分と重なりやすい。このため、本実施形態においては、１秒前後の周期を有する信号を除き、これより高い周波数を有する信号を抽出する。

なお、図６に、ドップラーセンサ出力信号を増幅器１０２で増幅した後、カットオフ周波数が２Ｈｚのハイパスフィルタ及びカットオフ周波数２０Ｈｚのローパスフィルタを通した後の信号が示される。図５に存在する低周波の信号成分がカットされている。

次に、第２段階について説明する。第１段階の処理によって、各種の信号処理が施されたドップラーセンサ出力信号を用いて、心拍数抽出部１０６は、各周波数成分に分解した成分信号の解析を行う。具体的には、データ取得部１０６１はデジタルフィルタ１０５から入力されたドップラーセンサ出力信号を取得する。そして、ウェーブレット解析部１０６２は、取得されたドップラーセンサ出力信号を連続ウェーブレット変換して各擬似周波数の成分信号に分解して解析する。

第１段階の処理によって、ドップラーセンサ出力信号は、主な雑音成分等が取り除かれ、心拍による身体の動きの速度変化の影響による周波数成分を含む信号となった。しかし、このドップラーセンサ出力信号の中にも雑音となる身体の動きの影響による周波数成分が含まれる場合がある。また、心拍による身体の動きの速度変化の影響による周波数成分の周波数の高さやその振幅の大きさには個人差がある。さらに、同一の被験者Ｍであっても、ドップラーセンサの放射波Ｔｘが放射される位置又は放射波Ｔｘの放射方向に対する身体の向きが異なると、心拍による身体の動きの速度変化の影響による周波数成分の周波数の高さやその振幅の大きさは異なる。このように、心拍による身体の動きの速度変化の影響による周波数成分の周波数は被験者Ｍおよびドップラーセンサのセンシング環境によって異なる。このため、本実施形態に係る心拍数検出装置１００は、心拍による身体の動きの速度変化の影響による周波数成分を含む信号を複数の周波数成分に分解する。具体的には、ウェーブレット解析部１０６２は、心拍による身体の動きの速度変化の影響による周波数成分を含むドップラーセンサ出力信号を連続ウェーブレット変換により複数の擬似周波数における成分信号に分解する（Ｓ１０７）。図７〜図１０には、上記ドップラーセンサ出力信号に連続ウェーブレット変換を行ったときの５Ｈｚ、１０Ｈｚ、１５Ｈｚ、２０Ｈｚの各擬似周波数におけるウェーブレット係数が示される。ウェーブレット係数は、信号強度に比例する値となる。

次に、第３段階について説明する。第２段階により、心拍による身体の動きの速度変化の影響による周波数成分を含むドップラーセンサ出力信号は、複数の周波数における成分信号に分解された。第３段階は、これら複数の成分信号の中から、心拍による身体の動きの速度変化の影響による周期的な振幅変動を含む成分信号を選択し、選択した成分信号から心拍数を検出する段階である。

心拍による身体の動きに依存したドップラーセンサ出力信号の振幅に比べ、一般的に呼吸や身体の揺れ等に依存した振幅の方が大きい。このため、元の信号の自己相関からは心拍の周期を明確に抽出することが困難な場合もある。このため、本実施形態においては、第２段階において求めた、複数の擬似周波数におけるそれぞれのウェーブレット係数を求めて、それぞれの周波数成分における自己相関関数をとることにより、周期性が低く振幅の大きい周波数成分の影響を分離することができる。

一般的に、人間の心拍の時間間隔には上限と下限が存在する。このため、自己相関解析部１０６３は、この上限と下限との間の範囲内における自己相関関数を図１１〜図１４に示すように求める（Ｓ１０９）。本実施形態においては、これらの自己相関関数の中から最も周期性が高いものを選択し、選択した自己相関関数に基づいて心拍数を検出する。例えば、ここでは周期性の判断に、自己相関関数の最も大きいピークと２番目に大きいピークの値の比が用いられる。

図１５に示されるように、自己相関関数には、一般的に複数のピークが存在する。最も大きい振幅を有するピークを第一ピークと称し、２番目に大きい振幅を有するピークを第二ピークと称する。このとき、自己相関解析部１０６３は、第一ピークと第二ピークの比を評価関数として周期性の評価を実行し（Ｓ１１１）、複数の成分信号についての自己相関関数の中から心拍による身体の動きの速度変化の影響による周波数成分を最も含む成分信号の自己相関関数を選択する。すなわち、評価関数は次のように表される。
評価関数＝（第一ピークの振幅の大きさ）÷（第二ピークの振幅の大きさ）

ここで、心拍数検出部１０６４は、自己相関解析部１０６３により選択された、評価関数の最も高い擬似周波数のウェーブレット係数の自己相関関数に基づいて、選択された自己相関関数の第一ピークの時間を心拍間隔とする（Ｓ１１３）。例えば、図１１〜図１４の例においては、５Ｈｚ、１０Ｈｚ、１５Ｈｚ、２０Ｈｚのそれぞれの擬似周波数における評価関数の評価値は、それぞれ１．４０、２．０３、２．８４、３．３６となる。すなわち、最も評価値の高い２０Ｈｚにおける自己相関関数を用いて心拍数は求められる。なお、若し自己相関関数のピークが１つしかない場合には、身体の揺れなどの影響を含む低い擬似周波数帯における信号成分の自己相関関数である場合がある。このため、その擬似周波数におけるウェーブレット係数の自己相関関数については評価の対象外とする。心拍数検出部１０６４は、自己相関関数のピークの時間を心拍間隔として、単位時間（例えば１分間）当たりの心拍数を求める（Ｓ１１５）。心拍数検出部１０６４は、求めた心拍数の情報を心拍数表示部１０７に入力する。心拍数表示部１０７は、入力された情報に基づいて心拍数の情報を提供する。

＜４．効果の例＞
以上説明してきたように、本実施形態に係る心拍数検出装置１００によれば、取得したドップラーセンサ出力信号から雑音成分をフィルタリングにより取り除き、さらに複数の擬似周波数における成分信号に分解する。そして複数の成分信号の中から、心拍による身体の動きの速度変化による周波数成分を含み、最も周期性が高い成分信号が選択される。心拍数検出装置１００は、選択された成分信号の自己相関関数のピーク値に基づいて心拍間隔を抽出する。そして心拍数検出装置１００は、この心拍間隔から単位時間当たりの心拍数を検出することができる。

複数の成分信号への分解前のドップラーセンサ出力信号は、心拍による影響以外の雑音、例えば身体の揺れなどの影響による周波数成分を含んでいる場合もある。ここで身体の揺れなどの影響による周波数成分は、心拍による影響と比べて一般的に振幅が大きい。このため、このドップラーセンサ出力信号の自己相関を取ったとしても雑音による影響を受けて周期性の検出は困難である場合がある。また、心拍による身体の動きの速度変化の影響による周波数成分の周波数の高さやその振幅の大きさは、個人差やセンシング環境によって異なる。そこで本実施形態においては、ドップラーセンサ出力信号を各周波数における成分信号に分解し、この成分信号毎に自己相関を解析する。また、この複数の成分信号の中から、心拍による影響を含み周期性の高い成分信号を選択し、選択された成分信号の自己相関関数のピーク値に基づいて心拍数が検出される。かかる構成により、心拍数検出装置１００は、より精度の高い心拍数の情報を提供することができる。

なお、心拍数検出の精度について検証した実験結果が図１６に示される。ここでは、三人の被験者について、本実施形態で検出した心拍数の精度と比較例の条件で検出した心拍数の精度とが示される。なお、比較例は、適切な振幅に増幅処理を行った後のドップラーセンサ出力信号の自己相関関数から心拍数を検出する方法が用いられた。なお、エラー率は、本実施形態と比較例のそれぞれについて、脈波計から求めた実際の心拍数の値との誤差率を示す。この図１６から、被験者Ａ、被験者Ｂ、および被験者Ｃのいずれについても、本実施形態の方法を用いて検出した心拍数が比較例の方法により検出した心拍数よりも精度が向上していることがわかる。

以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。

例えば、上記実施形態では、心拍を検出する装置について説明したが、本発明はかかる例に限定されない。例えば、呼吸、貧乏ゆすり、震え、筋音など他の生体情報を測定することとしてもよい。フィルタの通過帯域の周波数などのパラメータを変更することにより、周期性を有するものの検知に広く応用できることはいうまでもない。

また、上記実施形態では、求められた心拍数をそのままの値で提供することとしたが、本発明はかかる例に限定されない。例えば、所定期間の心拍数の中央値をとることにより、異常値の影響を低減するなど様々な精度向上のための手法が用いられてよい。

また、上記実施形態では、自己相関関数の第一ピークの値と第二ピークの値との比に基づいて周期性の高い成分信号が選択されたが、本発明はかかる例に限定されない。例えば、心拍間隔の範囲における第一ピークの値と、自己相関関数の時間が０のときの値との比に基づいて成分信号が選択されてもよい。その他、自己相関関数の周期性を評価するための手法が用いられてよい。

尚、本明細書において、フローチャートに記述されたステップは、記載された順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくとも、並列的に又は個別的に実行される処理をも含む。また時系列的に処理されるステップでも、場合によっては適宜順序を変更することが可能であることは言うまでもない。

１００心拍数検出装置
１０１ドップラーセンサ
１０２増幅器
１０３アナログフィルタ
１０４Ａ／Ｄ変換器
１０５デジタルフィルタ
１０６心拍数抽出部
１０７心拍数表示部

Claims

被験者に対して放射した電磁波又は超音波である放射波の周波数と、前記放射波が前記被験者により反射した反射波の周波数との差分の周波数を有するドップラーセンサ出力信号を取得するデータ取得部と、
前記ドップラーセンサ出力信号を複数の周波数における成分信号に分解し、複数の前記成分信号についてそれぞれ信号強度の時間変動を解析する時間周波数解析部と、
複数の前記成分信号についてそれぞれ信号強度の自己相関を解析する自己相関解析部と、
前記複数の成分信号の中から自己相関関数に基づいて周期性の最も高い成分信号を選択し、選択された成分信号の前記自己相関の解析結果に基づいて前記被験者の心拍数を検出する心拍数検出部と、
を備え、
前記心拍数検出部は、自己相関関数のピークのうち最も振幅の大きい第一ピークの値と２番目に振幅の大きい第二ピークの値との比に基づいて、周期性が最も高い成分信号を選択する、心拍数検出装置。
前記ドップラーセンサ出力信号から、心拍に起因した前記被験者の身体の動きの速度変化に応じた周波数成分を含む信号を抽出するフィルタ部、
をさらに備えることを特徴とする、請求項１に記載の心拍数検出装置。
前記ドップラーセンサ出力信号から、心拍に起因した前記被験者の身体の動きの位置変化に応じた周波数成分を含む信号を除去するフィルタ部、
をさらに備えることを特徴とする、請求項１に記載の心拍数検出装置。
前記心拍数検出部は、前記第一ピークの時間を心拍間隔として前記心拍数を検出する、請求項１〜３のいずれか一項に記載の心拍数検出装置。
被験者に対して電磁波又は超音波である放射波を放射するステップと、
前記放射波が前記被験者により反射した反射波を検出するステップと、
前記放射波の周波数と前記反射波の周波数との差分の周波数を有するドップラーセンサ出力信号を取得するステップと、
前記ドップラーセンサ出力信号を複数の周波数における成分信号に分解するステップと、
複数の前記成分信号についてそれぞれ信号強度の時間変動を解析するステップと、
複数の前記成分信号についてそれぞれ前記信号強度の自己相関を解析するステップと、
複数の前記成分信号の中から自己相関関数のピークのうち最も振幅の大きい第一ピークの値と２番目に振幅の大きい第二ピークの値との比に基づいて、周期性が最も高い成分信号を選択するステップと、
前記選択された成分信号の前記自己相関の解析結果に基づいて、前記被験者の心拍数を検出するステップと、
を含む、心拍数検出方法。
コンピュータを、
被験者に対して放射した電磁波又は超音波である放射波の周波数と、前記放射波が前記被験者により反射した反射波の周波数との差分の周波数を有するドップラーセンサ出力信号を取得するデータ取得部と、
前記ドップラーセンサ出力信号を複数の周波数における成分信号に分解し、複数の前記成分信号についてそれぞれ信号強度の時間変動を解析する時間周波数解析部と、
複数の前記成分信号についてそれぞれ信号強度の自己相関を解析する自己相関解析部と、
前記複数の成分信号の中から最も自己相関関数の周期性が高い成分信号を選択し、選択された成分信号の前記自己相関の解析結果に基づいて前記被験者の心拍数を検出する心拍数検出部と、
を備え、
前記心拍数検出部は、自己相関関数のピークのうち最も振幅の大きい第一ピークの値と２番目に振幅の大きい第二ピークの値との比に基づいて、周期性が最も高い成分信号を選択する、心拍数検出装置として機能させるためのプログラム。