(57)【特許請求の範囲】
【請求項１】
生体の振動に応じてドップラーセンサにより検出されたビート信号を一次元信号に変換する信号変換部と、
前記一次元信号の周波数を推定する周波数推定部と、
推定された前記一次元信号の周波数に対応する周期に相当する区間長を有する基準信号を前記一次元信号から切り出し、前記基準信号と前記周波数を有する正弦波との相関係数を算出することにより、前記一次元信号の基準時刻を推定する基準時刻推定部と、
前記基準時刻を用いて前記ビート信号の振動情報を特定する振動情報特定部と、
特定された前記振動情報に基づいて前記生体の睡眠深度を推定する睡眠深度推定部と、
を備える、睡眠深度推定装置。
【請求項２】
前記基準時刻推定部は、前記基準信号と前記正弦波との相関係数が最大となるときの前記正弦波の位相に基づいて前記基準時刻を推定する、請求項１に記載の睡眠深度推定装置。
【請求項３】
前記基準時刻推定部は、前記基準時刻を推定する際に用いる前記基準信号の区間において過去に推定された複数の基準時刻を抽出し、抽出された前記複数の基準時刻の分布に基づいて前記基準時刻を特定する、請求項１または２に記載の睡眠深度推定装置。
【請求項４】
前記振動情報特定部は、前記基準時刻における二次元ベクトルで表現された前記ビート信号のスカラー量を、前記ビート信号の振幅として特定する、請求項１〜３のいずれか１項に記載の睡眠深度推定装置。
【請求項５】
前記振動情報特定部は、前記基準時刻における前記一次元信号の振幅を、前記ビート信号の振幅として特定する、請求項１〜３のいずれか１項に記載の睡眠深度推定装置。
【請求項６】
前記振動情報特定部は、連続する２つの前記基準時刻の間隔を、前記ビート信号の周期として特定する、請求項１〜５のいずれか１項に記載の睡眠深度推定装置。
【請求項７】
前記睡眠深度推定部は、前記生体の体動の振動情報に基づいて前記生体の睡眠深度を推定する、請求項１〜６のいずれか１項に記載の睡眠深度推定装置。
【請求項８】
前記睡眠深度推定部は、前記生体の睡眠の継続時間に応じて前記生体の睡眠深度の推定条件を変更する、請求項１〜７のいずれか１項に記載の睡眠深度推定装置。
【請求項９】
前記ビート信号からカットオフ周波数より大きい周波数の成分を減少させるフィルタ処理を行うフィルタ部をさらに備える、請求項１〜８のいずれか１項に記載の睡眠深度推定装置。
【請求項１０】
前記フィルタ部の前記カットオフ周波数は１．５Ｈｚ以上である、請求項９に記載の睡眠深度推定装置。
【請求項１１】
前記生体の振動は、前記生体の呼吸による振動を含む、請求項１〜１０のいずれか１項に記載の睡眠深度推定装置。
【請求項１２】
生体の振動に応じてドップラーセンサにより検出されたビート信号を一次元信号に変換するステップと、
前記一次元信号の周波数を推定するステップと、
推定された前記一次元信号の周波数に対応する周期に相当する区間長を有する基準信号を前記一次元信号から切り出し、前記基準信号と前記周波数を有する正弦波との相関係数を算出することにより、前記基準信号の基準時刻を推定するステップと、
前記基準時刻を用いて前記ビート信号の振幅および周期を含む振動情報を特定するステップと、
特定された前記振動情報に基づいて前記生体の睡眠深度を推定するステップと、
を含む、睡眠深度推定方法。
【請求項１３】
コンピュータを、
生体の振動に応じてドップラーセンサにより検出されたビート信号を一次元信号に変換する信号変換部と、
前記一次元信号の周波数を推定する周波数推定部と、
推定された前記一次元信号の周波数に対応する周期に相当する区間長を有する基準信号を前記一次元信号から切り出し、前記基準信号と前記周波数を有する正弦波との相関係数を算出することにより、前記基準信号の基準時刻を推定する基準時刻推定部と、
前記基準時刻を用いて前記ビート信号の振幅および周期を含む振動情報を特定する振動情報特定部と、
特定された前記振動情報に基づいて前記生体の睡眠深度を推定する睡眠深度推定部と、
として機能させるための、プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、睡眠深度推定装置、睡眠深度推定方法、およびプログラムに関する。
【背景技術】
【０００２】
近年、被験者の睡眠状態等の健康状態を判定するために、被験者の呼吸や体動などの振動を非接触で検出し、検出された振動情報に基づき睡眠深度を推定する技術が開発されている。例えば、下記特許文献１には、被験者が使用する寝具に備えつけられたエアパッドを介して被験者の振動を計測し、計測された信号から被験者の振動の周期を算出することにより被験者の睡眠深度を推定する技術が開示されている。また、下記特許文献２には、ドップラーセンサにより取得された被験者の振動を示すＩＱ信号（Ｉ：Ｉｎ−ｐｈａｓｅ、Ｑ：Ｑｕａｄｒａｔｕｒｅ−ｐｈａｓｅ）について、ＩＱ平面と呼ばれる二次元平面上における信号のノルムから振動の一周期ごとの波形の境界を特定することにより振動の周期を算出し、算出された振動の周期に基づいて被験者の睡眠深度を判定する技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００６−２６３０３２号公報
【特許文献２】特開２０１４−０１４７０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、上記特許文献１においては、振動波形に含まれる微細な揺らぎをローパスフィルタにより減少させているため、微細な揺らぎを含む振動の波形は正確に抽出されていない。よって、呼吸の周期の乱れなどを含む睡眠状態の変化を捉えることは困難である。また、上記特許文献２においては、振動を示すＩＱ信号にノイズが重畳している場合、または、振動が連続的で切れ目なく生じている場合、振動の周期の区切りとなる振動の停止区間が検出されにくくなるので、推定される振動の周期の精度は低くなり、それに応じて推定される睡眠深度の精度も低くなる。そのため、微細な揺らぎを含む振動の周期を精度高く推定し、被験者の睡眠深度の推定精度を向上させることが望まれる。
【０００５】
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、被験者の睡眠深度をより高い精度で推定することが可能な、新規かつ改良された睡眠深度推定装置を提供することにある。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、本発明のある観点によれば、生体の振動に応じてドップラーセンサにより検出されたビート信号を一次元信号に変換する信号変換部と、上記一次元信号の周波数を推定する周波数推定部と、推定された上記一次元信号の周波数に対応する周期に相当する区間長を有する基準信号を上記一次元信号から切り出し、上記基準信号と上記周波数を有する正弦波との相関係数を算出することにより、上記一次元信号の基準時刻を推定する基準時刻推定部と、上記基準時刻を用いて上記ビート信号の振動情報を特定する振動情報特定部と、特定された上記振動情報に基づいて上記生体の睡眠深度を推定する睡眠深度推定部と、を備える、睡眠深度推定装置が提供される。
【０００７】
上記基準時刻推定部は、上記基準信号と上記正弦波との相関係数が最大となるときの上記正弦波の位相に基づいて上記基準時刻を推定してもよい。
【０００８】
上記基準時刻推定部は、上記基準時刻を推定する際に用いる上記基準信号の区間において過去に推定された複数の基準時刻を抽出し、抽出された上記複数の基準時刻の分布に基づいて上記基準時刻を特定してもよい。
【０００９】
上記振動情報特定部は、上記基準時刻における二次元ベクトルで表現された上記ビート信号のスカラー量を、上記ビート信号の振幅として特定してもよい。
【００１０】
上記振動情報特定部は、上記基準時刻における上記一次元信号の振幅を、上記ビート信号の振幅として特定してもよい。
【００１１】
上記振動情報特定部は、連続する２つの上記基準時刻の間隔を、上記ビート信号の周期として特定してもよい。
【００１２】
上記睡眠深度推定部は、上記生体の体動の振動情報に基づいて上記生体の睡眠深度を推定してもよい。
【００１３】
上記睡眠深度推定部は、上記生体の睡眠の継続時間に応じて上記生体の睡眠深度の推定条件を変更してもよい。
【００１４】
上記睡眠深度推定装置は、上記ビート信号からカットオフ周波数より大きい周波数の成分を減少させるフィルタ処理を行うフィルタ部をさらに備えてもよい。
【００１５】
上記フィルタ部の上記カットオフ周波数は１．５Ｈｚ以上であってもよい。
【００１６】
上記生体の振動は、上記生体の呼吸による振動を含んでもよい。
【００１７】
また、上記課題を解決するために、本発明の別の観点によれば、生体の振動に応じてドップラーセンサにより検出されたビート信号を一次元信号に変換するステップと、上記一次元信号の周波数を推定するステップと、推定された上記一次元信号の周波数に対応する周期に相当する区間長を有する基準信号を上記一次元信号から切り出し、上記基準信号と上記周波数を有する正弦波との相関係数を算出することにより、上記基準信号の基準時刻を推定するステップと、上記基準時刻を用いて上記ビート信号の振幅および周期を含む振動情報を特定するステップと、特定された上記振動情報に基づいて上記生体の睡眠深度を推定するステップと、を含む、睡眠深度推定方法が提供される。
【００１８】
また、上記課題を解決するために、本発明の別の観点によれば、コンピュータを、生体の振動に応じてドップラーセンサにより検出されたビート信号を一次元信号に変換する信号変換部と、上記一次元信号の周波数を推定する周波数推定部と、推定された上記一次元信号の周波数に対応する周期に相当する区間長を有する基準信号を上記一次元信号から切り出し、上記基準信号と上記周波数を有する正弦波との相関係数を算出することにより、上記基準信号の基準時刻を推定する基準時刻推定部と、上記基準時刻を用いて上記ビート信号の振幅および周期を含む振動情報を特定する振動情報特定部と、特定された上記振動情報に基づいて上記生体の睡眠深度を推定する睡眠深度推定部と、として機能させるための、プログラムが提供される。
【発明の効果】
【００１９】
以上説明したように本発明によれば、被験者の睡眠深度をより高い精度で推定することが可能である。
【図面の簡単な説明】
【００２０】
【図１】本発明の一実施形態に係る睡眠深度推定システムの概要を示す図である。
【図２】本発明の一実施形態に係る睡眠深度推定装置の構成例を示すブロック図である。
【図３】通過帯域の異なるフィルタによりフィルタ処理されたビート信号から変換された一次元信号の波形を比較する例を示す図である。
【図４】ビート信号のＩＱ平面における軌跡の一例を示す図である。
【図５】ビート信号のＩＱ平面における軌跡について第２変換部による変換処理を実施する例を示す図である。
【図６】本発明の一実施形態に係る基準時刻推定部による基準時刻の推定方法の一例を示す図である。
【図７】本発明の一実施形態に係る基準時刻推定部による基準時刻の特定方法の一例を示す図である。
【図８】本発明の一実施形態に係る深度判別部による睡眠深度の判別フローを示すフローチャートの一例を示す図である。
【図９】睡眠深度の推定結果を出力したグラフの一例を示す図である。
【図１０】本発明の一実施形態に係る睡眠深度推定装置の動作例を示すフローチャートである。
【図１１】本発明の実施形態に係る睡眠深度推定装置のハードウェア構成例を示すブロック図である。
【発明を実施するための形態】
【００２１】
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００２２】
＜１．睡眠深度推定システムの構成例＞
図１は、本発明の一実施形態に係る睡眠深度推定システム１の概要を示す図である。図１を参照すると、睡眠深度推定システム１は、ドップラーセンサ２、および睡眠深度推定装置１０を有する。
【００２３】
ドップラーセンサ２は、図１に示すように、例えば部屋の天井等に設置され、検出エリアである部屋の屋内に向けて、光、電磁波、または音波等の放射波を放射し、物体（例えば図１においては人物Ｐ）により反射された反射波を受信する。このとき、反射波の周波数は、物体の振動等の運動により生じるドップラー効果により、放射波の周波数から変化する。ドップラーセンサ２は、放射波の周波数と反射波の周波数との差分の周波数を有するビート信号を生成する。ドップラーセンサ２は直交検波方式を搭載してもよく、その場合、ドップラーセンサ２は、余弦波成分（Ｉ成分）および正弦波成分（Ｑ成分）の２種類のビート信号を生成する。ドップラーセンサ２は、生成したビート信号を睡眠深度推定装置１０に出力する。
【００２４】
なお、ドップラーセンサ２は、物体の振動が検出可能であれば任意の位置に設置可能である。また、図１に示した例では、ドップラーセンサ２は放射波の送信部と反射波の受信部とが一体となって構成されているが、ドップラーセンサ２は、送信部と受信部とを分離した構成により実現されてもよい。また、ドップラーセンサ２が放射する放射波は、物体の振動によりドップラー効果を生じさせることが可能であれば、任意の周波数帯の波であってもよい。例えば、呼吸の振動をドップラーセンサにより検知する場合、２．４ＧＨｚ〜１０．５ＧＨｚ帯などの、呼吸の振動の変位と比較して波長が大きい波が用いられることが多い。しかし、本実施形態においては、呼吸等の振動の幅広い変位を対象として、かつ波形に含まれる微細なゆらぎを検出するために、２４ＧＨｚ帯などの準ミリ波やミリ波の帯域の波が用いられることが好ましい。
【００２５】
睡眠深度推定装置１０は、ドップラーセンサ２から出力されたビート信号から被験者（例えば人物Ｐ）の振動情報を特定し、特定された振動情報に基づいて被験者の睡眠深度を推定する装置である。睡眠深度推定装置１０は、例えば、ネットワーク上の一または複数の情報処理装置によって実現されてもよい。より具体的には、睡眠深度推定装置１０は、サーバやＰＣ（Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ）等により実現されてもよい。ここで、本明細書において振動情報とは、例えば物体の振動波形から得られる振幅や周期を意味する。例えば、睡眠深度推定装置１０は、人物Ｐの呼吸や体動等の振動から得られるビート信号から人物Ｐの振動情報を特定し、特定された振動情報を用いて人物Ｐの睡眠深度を推定することが可能である。睡眠深度推定装置１０は、例えば、不図示の通信部を介して、推定した睡眠深度を有線または無線により外部装置等に送信することが可能である。これにより、例えば、上記外部装置において、睡眠深度について解析を行うことが可能となる。睡眠深度の解析により、例えば、人物Ｐの睡眠状態や睡眠時のストレス等を把握することが可能となる。なお、睡眠深度を用いた解析については、上記のような外部装置ではなく、睡眠深度推定装置１０本体において行うことも可能である。
【００２６】
ここで、人物Ｐの振動から得られるビート信号には、人物Ｐの様々な運動に起因する信号が含まれ得る。例えば、人物Ｐの呼吸の振動のみならず、人物Ｐの体動や心拍などの振動もビート信号に含まれ得る。さらに、ドップラーセンサ２が生成したビート信号には、周囲の環境に起因するノイズ等も含まれ得る。そのため、例えば呼吸の振動を検出対象とする場合、呼吸の波形情報のみをビート信号から検出することが望ましい。
【００２７】
例えば、特開２００６−２６３０３２号公報においては、呼吸の周期を検出するために、振動信号に対して、呼吸の周波数帯に相当する０．１５〜０．３０Ｈｚのバンドパスフィルタによるフィルタ処理が施されている。しかしその場合、呼吸の一周期ごとに微細な揺らぎが生じている場合において、その揺らぎがバンドパスフィルタにより低減されてしまうため、微細な揺らぎを伴う振動を表現する波形を抽出することが困難となる。そのため、呼吸の周期の乱れなどを含む睡眠状態の変化等を検知することは困難である。
【００２８】
また、特開２０１４−０１４７０８号公報においては、ＩＱ平面上におけるビート信号のノルムを算出し、ノルムが略ゼロとなる時刻を呼吸の吸い込みと吐き出しの切り替わり時として検出することにより、呼吸の周期が推定されている。しかし、例えばビート信号にノイズが多く重畳している場合や、呼吸の動作が切れ目なく連続的である場合、ノルムが略ゼロとなる時刻の検出が困難となる。そのため、呼吸の周期を誤って推定してしまう可能性がある。それに応じて、被験者の睡眠深度の変化等も誤って推定されてしまうおそれがある。
【００２９】
そこで、上記事情を一着眼点にして、睡眠深度推定装置１０を創作するに至った。本発明の一実施形態に係る睡眠深度推定装置１０は、微細な揺らぎを伴う被験者の呼吸の振動状態を精度高く推定することが可能であるので、睡眠深度の推定精度をさらに高めることができる。以下、一実施形態に係る睡眠深度推定装置１０の構成についてより詳細に説明する。
【００３０】
＜２．睡眠深度推定装置の構成例＞
図２は、本発明の一実施形態に係る睡眠深度推定装置１０の構成例を示すブロック図である。図２を参照すると、睡眠深度推定装置１０は、呼吸推定部１００、体動推定部２００、および睡眠深度推定部３００を含む。図示された睡眠深度推定装置１０は、微細な揺らぎを伴う被験者の呼吸の振動情報、および被験者の体動の振動情報を推定し、推定された各振動情報に基づいて被験者の睡眠深度を推定することができる。以下、各構成の機能について説明する。
【００３１】
［２−１．呼吸推定部］
呼吸推定部１００は、ビート信号取得部１１０、フィルタ部１２０、信号変換部１３０、周波数推定部１４０、基準時刻推定部１５０、および振動情報特定部１６０を含む。呼吸推定部１００は、取得したビート信号を一次元信号に変換し、変換した一次元信号の周波数を推定し、推定された周波数からビート信号の周期を算出するための一次元信号における基準位置を推定する。そして、呼吸推定部１００は、推定された基準位置を用いて、呼吸の振動の周期、または振幅を含む振動情報を特定する。特定された振動情報は、睡眠深度推定部３００へ出力される。以下、呼吸推定部１００に含まれる各機能の構成について説明する。
【００３２】
（ビート信号取得部）
ビート信号取得部１１０は、ドップラーセンサ２から出力されるビート信号Ｄ（ｔ）を取得する。ドップラーセンサ２から出力されるビート信号Ｄ（ｔ）はＩ成分とＱ成分の２波の成分を有し、振幅をＡ（ｔ）、波長をλ、時刻ｔにおけるドップラーセンサ２と対象物体（例えば図１に示された人物Ｐ）との距離をｄ（ｔ）、初期位相をφ０、直流成分をＯ、ノイズ成分をｗとすると、ビート信号Ｄ（ｔ）は、下記の数式１のように表現される。
【００３３】
【数１】
0006515670.tif 000002
【００３４】
ビート信号取得部１１０は、取得したビート信号Ｄ（ｔ）をフィルタ部１２０に出力する。また、ビート信号取得部１１０は、取得したビート信号Ｄ（ｔ）を時系列に不図示の記憶部に記憶してもよい。睡眠時における呼吸の基本周波数は、約０．１５〜０．３０Ｈｚと言われている。しかし、実際の呼吸の振動を表す信号は正弦波ではなく、複雑な波形を有している。そのため、微細な揺らぎを含む呼吸の周期を正確に抽出するために、サンプリング周波数は、少なくとも呼吸の基本周波数である０．３０Ｈｚの５倍以上である１．５Ｈｚ以上の周波数成分であることが好ましい。より具体的には、サンプリング周波数は、３Ｈｚ以上であることがさらに好ましい。ちなみに、折り返し歪みを抑制するためには、信号のサンプリングの前にナイキスト周波数以下の周波数成分を得るためのローパスフィルタを適用する必要がある。これは、信号からナイキスト周波数以上の成分をある程度除去しないと、波形の歪みにより適切な推定が難しくなるためである。いずれにせよ、呼吸の振動波形の抽出のためには、サンプリング周波数が１．５Ｈｚ以上である必要がある。そして、より正確に呼吸の周期を抽出するためには、サンプリング周波数を、３Ｈｚ〜２０Ｈｚ程度に設定することが好ましい。これよりさらに高い周波数とすると、呼吸の周期の推定の精度が改善する可能性もある一方で、帯域が広がるに従い信号に重畳される雑音レベルも多くなる。そのため、サンプリング周波数を２０Ｈｚより高く設定しても、呼吸の周期の推定精度の改善以上に、処理量の増加によるデメリットが大きくなると考えられる。ただし、呼吸の周期の推定にとって有効である周波数成分がサンプリングできていれば、サンプリング周波数を高く設定することは、処理量が増加する以外に、周期の推定精度に影響を及ぼさない。なお、本実施形態においては、サンプリング周波数は２０Ｈｚである。
【００３５】
（フィルタ部）
フィルタ部１２０は、ビート信号取得部１１０が取得したビート信号Ｄ（ｔ）に含まれる直流成分Ｏなどの低周波成分や、ノイズ成分ｗなどの高周波成分を減少または除去するフィルタ処理を実施し、フィルタ処理されたビート信号Ｄ（ｔ）を信号変換部１３０に出力する。例えば、フィルタ部１２０は、ビート信号Ｄ（ｔ）に含まれる直流成分Ｏを減少または除去するハイパスフィルタを備えてもよい。ビート信号Ｄ（ｔ）の振幅が微弱である場合において、ビート信号Ｄ（ｔ）に直流成分Ｏが含まれたまま不図示の増幅器等により増幅させると、直流成分Ｏも増幅される。そうすると、ビート信号Ｄ（ｔ）のうち呼吸の振動に寄与する特徴を示す信号の波形のみを抽出することが困難である。そこで、ハイパスフィルタによって直流成分Ｏを減少または除去することにより、ビート信号Ｄ（ｔ）のうち、呼吸の振動に寄与する特徴を示す信号の波形を適切に増幅させることが可能となる。上記のようなハイパスフィルタは、例えば睡眠深度推定装置１０に備えられるフィルタ回路等のアナログフィルタによって実現されてもよいし、ＩＩＲ（Ｉｎｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ：無限インパルス応答）フィルタ等のデジタルフィルタによって実現されてもよい。また、フィルタ部１２０は、直流成分Ｏを特定して除去するように減算回路を配置するフィルタ回路によって直流成分Ｏを減少または除去してもよい。
【００３６】
また、フィルタ部１２０は、ビート信号Ｄ（ｔ）に含まれるノイズ成分ｗを減少または除去するローパスフィルタを備えてもよい。ローパスフィルタは、例えばカットオフ周波数が１．５Ｈｚ以上であるローパスフィルタであってもよい。一般的に、睡眠時における呼吸の基本周波数は、上述したように、約０．１５〜０．３０Ｈｚと言われている。しかし、被験者の呼吸は、被験者の睡眠状態によって、揺らぎが生じる。例えば、被験者の睡眠が浅い状態である場合、または発熱等により被験者の呼吸が不安定な場合は、呼吸に揺らぎが生じ得る。呼吸の揺らぎに起因する振動は、呼吸の基本周波数よりも高い周波数成分を有する。このような呼吸の揺らぎを正確に抽出するために、ローパスフィルタのカットオフ周波数を、呼吸の基本周波数よりも高い値（例えば５〜１０倍以上）に設定することが好ましい。
【００３７】
図３は、通過帯域の異なるフィルタによりフィルタ処理されたビート信号Ｄ（ｔ）から変換された一次元信号ｒ（ｔ）の波形を比較する例を示す図である。波形Ｗ０は、フィルタ処理がされていないビート信号Ｄ（ｔ）から変換された一次元信号ｒ（ｔ）の波形である。波形Ｗ１は、通過帯域が０．１５〜０．３０Ｈｚであるバンドパスフィルタによるフィルタ処理がされたビート信号Ｄ（ｔ）から変換された一次元信号ｒ（ｔ）の波形である。波形Ｗ２は、通過帯域が３．０Ｈｚ以下であるローパスフィルタによるフィルタ処理がされたビート信号Ｄ（ｔ）から変換された一次元信号ｒ（ｔ）の波形である。なお、波形Ｗ０〜Ｗ２で示される信号の位相は同一である。
【００３８】
領域Ｒ１において３つの波形を比較してみると、波形Ｗ２は波形Ｗ０と略同一の形状を示しているのに対し、波形Ｗ１は波形Ｗ０とは異なる形状を示している。また、波形Ｗ０の下限のピーク位置を示す時刻Ｔ１の位置について、波形Ｗ２の下限のピーク位置を示す時刻は時刻Ｔ１と略同一であるのに対し、波形Ｗ１の下限のピーク位置を示す時刻は時刻Ｔ１と乖離していることが分かる。つまり、通過帯域が０．１５〜０．３０Ｈｚであるバンドパスフィルタによるフィルタ処理がされたビート信号Ｄ（ｔ）から変換された一次元信号ｒ（ｔ）の波形Ｗ１は、単純な正弦波に近い波形であり、実際の呼吸の振動を正確に表現していないと言える。この場合、例えば、呼吸の振動に乱れが生じ得るレム睡眠状態を、波形Ｗ１を用いて推定することは困難である。一方、通過帯域が３．０Ｈｚ以下であるローパスフィルタによるフィルタ処理がされたビート信号Ｄ（ｔ）から変換された一次元信号ｒ（ｔ）の波形Ｗ２は、波形Ｗ０からノイズが取り除かれた呼吸の振動の波形を示している。よって、カットオフ周波数を上記した１．５Ｈｚ以上の値である３Ｈｚとするローパスフィルタは、ノイズを取り除きつつ、微細な揺らぎの成分を含む呼吸の振動を抽出することが可能である。しかし、上記のようなローパスフィルタから抽出された波形Ｗ２は単純な正弦波形ではないため、例えば特開２００６−２６３０３２号公報に開示された手法では、呼吸の振動の周期を推定することは困難である。そのため、本発明は、後述する基準時刻推定部１５０等によって、微細な揺らぎの成分を含む呼吸の振動の周期を推定する。
【００３９】
なお、フィルタ部１２０を実現する上記のようなローパスフィルタのカットオフ周波数は、上述の通り１．５Ｈｚ以上が好ましく、より具体的には、３Ｈｚとすることが好ましいが、呼吸に含まれるより微細な揺らぎを含む振動を抽出したい場合、カットオフ周波数は３Ｈｚを超えてもよい。例えば、上記のカットオフ周波数は、１０Ｈｚ〜２０Ｈｚの間に設定されてもよい。これにより、呼吸に含まれる微細な揺らぎを含む振動を示す信号の周期をより正確に推定することが可能である。ただし、カットオフ周波数が２０Ｈｚを超えると、呼吸の振動には微細な揺らぎのみならずノイズが多く含まれ得る。そのため、フィルタ部１２０を実現するローパスフィルタのカットオフ周波数は２０Ｈｚ以下であることが好ましい。
【００４０】
（信号変換部）
信号変換部１３０は、フィルタ部１２０から出力されたビート信号Ｄ（ｔ）を一次元信号ｒ（ｔ）に変換する。例えば、信号変換部１３０は、ビート信号Ｄ（ｔ）のＩＱ平面における分布に基づいてビート信号Ｄ（ｔ）を一次元の第１候補信号Ｒ１（ｔ）に変換する、第１変換部を有してもよい。また、信号変換部１３０は、ビート信号Ｄ（ｔ）のＩＱ平面における位置の時間変化に基づいてビート信号Ｄ（ｔ）を一次元の第２候補信号Ｒ２（ｔ）に変換する、第２変換部を有してもよい。さらに、信号変換部１３０は、上記の２つの変換部を有する場合、第１候補信号Ｒ１（ｔ）と第２候補信号Ｒ２（ｔ）の少なくともいずれかを用いることにより１次元信号ｒ（ｔ）を決定する信号決定部を有してもよい。以下、上述した各機能部について説明する。
【００４１】
−第１変換部
第１変換部は、例えば、ビート信号Ｄ（ｔ）を表現する二次元ベクトルと、ビート信号Ｄ（ｔ）の共分散行列の最大固有値に対応する固有ベクトルとの内積を算出することにより、ビート信号Ｄ（ｔ）を第１候補信号Ｒ１（ｔ）に変換してもよい。つまり、第１変換部は、ビート信号Ｄ（ｔ）を最大固有値に対応する主成分方向に射影した第１候補信号Ｒ１（ｔ）を得ることができる。具体的には、第１候補信号Ｒ１（ｔ）は、下記の数式２により算出されてもよい。
【００４２】
【数２】
0006515670.tif 000003
【００４３】
ここで、ｂ（ｔ）はビート信号Ｄ（ｔ）のＩ成分とＱ成分を表現する二次元ベクトル、ｐ（ｔ）は、ビート信号Ｄ（ｔ）の共分散行列の最大固有値に対応する固有ベクトルである。なお、ビート信号Ｄ（ｔ）を表現する二次元ベクトルの原点は、例えばビート信号Ｄ（ｔ）の分布に基づいて推定されてもよい。例えば、過去数秒〜数十秒におけるビート信号Ｄ（ｔ）の軌跡の、ＩＱ平面の原点からの最大位置および最小位置の中点が、原点であってもよい。また、最小二乗法等によりビート信号Ｄ（ｔ）の軌跡を円や楕円として近似し、その近似した図形の中心点が原点であってもよい。また、フィルタ部１２０においてハイパスフィルタにより処理されたビート信号Ｄ（ｔ）はＩＱ平面の原点に近接するので、ＩＱ平面の原点が二次元ベクトルの原点であってもよい。
【００４４】
図４は、ビート信号Ｄ（ｔ）のＩＱ平面における軌跡の一例を示す図である。図４に示されるグラフＧ１を参照すると、ＩＱ平面上においてビート信号Ｄ（ｔ）の軌跡Ｔｒ１が示されている。例えば、呼吸の振動の変位がドップラーセンサから放射される放射波の波長の１／１０程度であるなど比較的小さい場合、ビート信号Ｄ（ｔ）は、軌跡Ｔｒ１のような扁平な８の字の分布となることが多い。このような場合、例えば軌跡Ｔｒ１について主成分分析を行った場合、軸Ｐｃ１が示す第１主成分に対応する最大固有値は、他の主成分と比較して大きな値を示す。第１変換部は、ビート信号Ｄ（ｔ）を第１主成分に射影することにより、第１候補信号Ｒ１（ｔ）を得ることができる。
【００４５】
なお、本実施形態に係る第１変換部は、軌跡Ｔｒ１について主成分分析等を行うことにより得られる固有値および固有ベクトルを用いてビート信号Ｄ（ｔ）を第１候補信号Ｒ１（ｔ）に変換する。しかし、例えば、軌跡Ｔｒ１のような扁平な８の字の分布について、上記扁平方向の軸方向を抽出し、上記軸方向における信号の動きを抽出することが可能であれば、第１変換部は主成分分析以外の方法を用いてビート信号Ｄ（ｔ）を第１候補信号Ｒ１（ｔ）に変換してもよい。
【００４６】
−第２変換部
第２変換部は、例えば、ビート信号Ｄ（ｔ）のＩＱ平面における分布に基づいて推定される分布の中心からのビート信号Ｄ（ｔ）の位置までの距離（振幅に相当）と、分布の中心を基準とするビート信号Ｄ（ｔ）の回転角度の変化量との積を算出することにより、ビート信号Ｄ（ｔ）を第２候補信号Ｒ２（ｔ）に変換してもよい。例えば、Ｒ２（ｔ）は、下記の数式３における連続関数（本実施形態においては、シグモイド関数）を利用した関数を用いて算出されてもよい。
【００４７】
【数３】
0006515670.tif 000004
【００４８】
ここで、Ａｍｐ（ｔ）はビート信号Ｄ（ｔ）の振幅、θ（ｔ）はビート信号Ｄ（ｔ）のＩＱ平面上における角度であり、θ’（ｔ）は、角度θ（ｔ）の単位時間の角度の変化量、または角度θ（ｔ）の時間微分値である。また、ａ、ｂ、およびｃはいずれも定数であり、自由に設定可能である。
【００４９】
図５は、ビート信号Ｄ（ｔ）のＩＱ平面における軌跡について第２変換部による変換処理を実施する例を示す図である。図５に示されるグラフＧ２を参照すると、軌跡Ｔｒ２上におけるビート信号Ｄ（ｔ）について、振幅Ａｍｐ（ｔ）および角度θ（ｔ）が示されている。振幅Ａｍｐ（ｔ）および角度θ（ｔ）の変化量θ’（ｔ）の積は、図５に示されている領域Ａｖ１の面積速度に相当する。この面積速度は、呼吸の振動の単位時間当たりの大きさに相当する。また、変化量θ’（ｔ）は、呼吸の振動の方向（息の吸い込み、吐き出し）により正負の値を取るため、呼吸の振動の方向を表現できる。そのため、振幅Ａｍｐ（ｔ）と変化量θ’（ｔ）の積により、呼吸運動の状態を的確に表現することが可能となる。
【００５０】
なお、振幅Ａｍｐ（ｔ）は、ビート信号Ｄ（ｔ）のＩＱ平面における分布に基づいて推定される分布の中心からのビート信号Ｄ（ｔ）の位置までの距離に相当し、また、角度θ（ｔ）は、分布の中心を基準とするビート信号Ｄ（ｔ）の回転角度に相当する。ここで、分布の中心とは、例えば図５に示した中心点Ｃ１を指す。ここで、ビート信号Ｄ（ｔ）の分布の中心は、ビート信号Ｄ（ｔ）の分布に基づいて推定されてもよい。例えば、過去数秒〜数十秒におけるビート信号Ｄ（ｔ）の分布について、ＩＱ平面の原点からの最大位置と最小位置との中間位置や、分布の中央に相当する位置が、ビート信号Ｄ（ｔ）の分布の中心であってもよい。また、最小二乗法等によりビート信号Ｄ（ｔ）の軌跡を円や楕円として近似し、その近似した図形の中心点が分布の中心であってもよい。また、フィルタ部１２０においてハイパスフィルタにより処理されたビート信号Ｄ（ｔ）はＩＱ平面の原点に近接するように分布するので、ＩＱ平面の原点がビート信号Ｄ（ｔ）の分布の中心であってもよい。
【００５１】
−信号決定部
信号決定部は、例えば、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）の一方を一次元信号ｒ（ｔ）と決定してもよい。また、信号決定部は、第１候補信号Ｒ１（ｔ）と第２候補信号Ｒ２（ｔ）とを統合させた信号を一次元信号ｒ（ｔ）と決定してもよい。
【００５２】
図４に示した呼吸の振動の変位が小さい場合のビート信号Ｄ（ｔ）の軌跡Ｔｒ１においては、分布の中心からの角度θ（ｔ）の変化量θ’（ｔ）は、図４に示したビート信号Ｄ（ｔ）の軌跡Ｔｒ２の場合と比較して、ノイズや量子化誤差の影響により精度が低くなる。つまり、呼吸の振動が小さい場合は、第２候補信号Ｒ２（ｔ）よりも、第１候補信号Ｒ１（ｔ）の方が、呼吸の振動の波形を正確に表現し得る。一方で、図５に示した呼吸の振動の変位が大きい場合のビート信号Ｄ（ｔ）の軌跡Ｔｒ２においては、第１候補信号Ｒ１（ｔ）よりも、第２候補信号Ｒ２（ｔ）の方が、呼吸の振動の波形を正確に表現し得る。これは、主成分方向に射影される第１候補信号Ｒ１（ｔ）の波は呼吸の周期よりも短い周期（高周波）の振動を含んでしまうが、第２候補信号Ｒ２（ｔ）は、ビート信号Ｄ（ｔ）の逐次的な変化を追いかけることにより得られるので、ビート信号Ｄ（ｔ）の分布に起因する高周波成分を含まないためである。よって、呼吸の振動の変位が大きい場合は、第１候補信号Ｒ１（ｔ）よりも、第２候補信号Ｒ２（ｔ）の方が、呼吸の振動の波形を正確に表現し得る。
【００５３】
信号決定部は、第１候補信号Ｒ１（ｔ）または第２候補信号Ｒ２（ｔ）についてそれぞれ比較し、その比較結果に基づいて出力する一次元信号ｒ（ｔ）を決定してもよい。信号決定部は、例えば、どちらの候補信号を優先的に用いるかを決定するためのパラメータとして、評価値ｓを算出してもよい。評価値ｓは、例えば各候補信号の周波数の比較結果や、ビート信号Ｄ（ｔ）の角度の変化量の偏差等の指標値から決定されてもよい。例えば、信号決定部は、上記のような指標値等について多次元のマッピングを行い、得られた分布に基づいて評価値ｓを算出してもよい。また、信号決定部は、機械学習等のデータ処理手法を用いて、上記のような指標値を入力することにより得られる出力結果を評価値ｓとして取得してもよい。
【００５４】
また、信号決定部は、例えば、一次元信号ｒ（ｔ）として一時刻前に一方の候補信号を選択したあとに他方の候補信号に切り替える場合、または、双方の候補信号を統合する場合、双方の候補信号の位相を整合させる必要がある。例えば、第１候補信号Ｒ１（ｔ）の位相φ１と第２候補信号Ｒ２（ｔ）の位相φ２は必ずしも一致しているとは限らない。そのため、信号決定部は、一方の候補信号の位相を他方の候補信号の位相に整合させる。これにより、一次元信号ｒ（ｔ）の信号の切り替え処理や、双方の候補信号の統合処理において、各候補信号の位相の不整合による波形の乱れを防ぐことができる。なお、本実施形態では、第１候補信号Ｒ１（ｔ）の位相を第２候補信号Ｒ２（ｔ）の位相に整合させる処理が実施される。例えば、一方の候補信号をヒルベルト変換により解析信号に変換し、解析信号と他方の候補信号とを積算させて得られる２つの信号の位相差を算出することにより、位相の整合処理を実施することが可能である。
【００５５】
さらに、信号決定部は、上記の位相の整合処理を実施した場合、双方の候補信号を統合する処理を実施してもよい。例えば、信号決定部は、上記のような評価値ｓを重みとして用いて算出された第１候補信号Ｒ１ａ（ｔ）と第２候補信号Ｒ２（ｔ）との加重和を、一次元信号ｒ（ｔ）として決定することができる。
【００５６】
以上より、信号変換部１３０は、呼吸の振動の変位の大小にかかわらず適切な変換手法を選択することができる。信号変換部１３０によりビート信号Ｄ（ｔ）から変換された一次元信号ｒ（ｔ）は、周波数推定部１４０に出力される。
【００５７】
（周波数推定部）
周波数推定部１４０は、信号変換部１３０から出力された一次元信号ｒ（ｔ）の有する周波数を推定する。例えば、周波数推定部１４０は、呼吸周期の大きな揺らぎに合わせるために、一次元信号ｒ（ｔ）と参照信号との位相差、および位相差の時間変化に基づいて一次元信号ｒ（ｔ）の周波数を推定する、第１周波数推定部を有してもよい。また、周波数推定部１４０は、第１周波数推定部により推定された推定周波数に対応する周期に相当する区間長を有する比較信号を一次元信号ｒ（ｔ）から切り出し、比較信号と上記の区間長を有する一次元信号ｒ（ｔ）との相関係数を算出することにより一次元信号ｒ（ｔ）の周波数を推定する、第２周波数推定部を有してもよい。以下、第１周波数推定部により推定された一次元信号ｒ（ｔ）の周波数を第１推定周波数と、また、第２周波数推定部により推定された一次元信号ｒ（ｔ）の周波数を第２推定周波数と呼称する。
【００５８】
−第１周波数推定部
第１の周波数推定部は、例えば、時刻ｔｋにおける一次元信号ｒ（ｔｋ）と、参照信号とを積算し、積算された信号の低周波成分を、ローパスフィルタ等を用いて抽出することにより、一次元信号ｒ（ｔk）と参照信号との位相差ψ（ｔk）を算出する。ここでは参照信号として、位相変化による周期推定を容易とするため、位相とその時間変化をパラメータにもつ正弦波、または余弦波が用いられる。その場合、参照信号は、時刻ｔk−１に推定された周波数ｆＡ（ｔk−１）を有する正弦波、あるいは余弦波の信号である。そして、第１の周波数推定部は、位相差ψ（ｔk）と、時刻ｔk−１に算出された位相差ψ（ｔk−１）との変化分を算出する。位相差の時間変化分ψ（ｔk）−ψ（ｔk−１）は、一次元信号ｒ（ｔ）の周波数の変化によって生じたものと仮定して、第１周波数推定部は、参照信号の周波数ｆＡ（ｔk−１）に、位相差の時間変化分ψ（ｔk）−ψ（ｔk−１）に相当する周波数を加算して、第１推定周波数ｆＡ（ｔk）を出力する。
【００５９】
なお、上記位相差の時間変化分が所定の閾値以上であった場合、周波数のフィードバック値が大きく見積もられやすいため、第１推定周波数ｆＡ（ｔk）が振動を伴って収束しないといった問題がある。そのため、参照信号の周波数ｆＡ（ｔk−１）に加算する周波数の値は、位相差の時間変化分ψ（ｔk）−ψ（ｔk−１）に相当する周波数に、１未満の係数を積算した値でもよい。これにより、比較的短い時間でｆＡ（ｔｋ）が収束することが可能である。
【００６０】
また、算出された第１推定周波数ｆＡ（ｔk）が所定の周波数帯域に存在しない場合、参照信号の周波数ｆＡ（ｔk−１）を、第１推定周波数ｆＡ（ｔk）として出力してもよい。ここで所定の周波数帯域とは、例えば、本実施形態においては、呼吸運動の周期に相当する周波数帯域を意味する。
【００６１】
−第２の周波数推定部
第２の周波数推定部は、入力された一次元信号ｒ（ｔ）の相関分析を行うことにより、第１の周波数推定部により得られた第１推定周波数ｆＡ（ｔk）より局所的な変化に対して精度の高い第２推定周波数ｆＢ（ｔk）を推定する。まず、第２周波数推定部は、一次元信号ｒ（ｔ）から、時刻ｔk−ＴＡから時刻ｔkまでの区間を切り出し、これを比較信号として扱う。ＴＡは、先に推定された第１の周波数ｆＡ（ｔk−１）に対応する周期である。次に、第２周波数推定部は、切り出された比較信号を用いて、相関係数を算出する。相関係数の算出範囲は、ｔk−２ＴＡからｔk−ＴＡの範囲を含むことが好ましい。そして第２周波数推定部は、上記の算出範囲において、区間長ＴＡを有する一次元信号ｒ（ｔ）と、比較信号との相関係数Ｒ（τ）が最大となるタイムラグτｍａｘを算出する。このとき、第２周波数推定部は、τ＝τｍａｘのときの相関係数Ｒ（τｍａｘ）を不図示の記憶部に記憶してもよい。第２周波数推定部は、相関係数算出範囲において算出される相関係数が最大となるときのタイムラグτｍａｘから、第２推定周波数ｆＢ（ｔｋ）を算出する。例えば、第２推定周波数ｆＢ（ｔｋ）は、タイムラグτｍａｘの逆数であってもよい。
【００６２】
−判定部
周波数推定部１４０は、上述した２つの推定された周波数について、いずれの周波数を用いるかを判定する判定部を有してもよい。例えば、判定部は、２つの推定された周波数の差が所定の閾値以内であるかに基づいて、いずれの周波数を推定周波数ｆＥ（ｔ）として用いるかを判定してもよい。また、第２周波数推定部で算出されたタイムラグτｍａｘの値に基づいて、第１推定周波数ｆＡ（ｔk）と第２推定周波数ｆＢ（ｔｋ）のいずれの周波数を用いるかを判定してもよい。
【００６３】
以上より、周波数推定部１４０は、おおよその周波数を推定できる第１周波数推定部と、イレギュラーな変動には追従することが困難であるが、波形の局所的な変化を検出することにより周波数を推定する第２周波数推定部とのどちらを用いるかを判定部により判定することにより、微細な揺らぎを多く含む波形の周波数を、１回の呼吸ごとに高い精度で推定することができる。周波数推定部１４０により推定された推定周波数ｆＥ（ｔ）は、基準時刻推定部１５０に出力される。
【００６４】
（基準時刻推定部）
基準時刻推定部１５０は、推定周波数ｆＥ（ｔ）に基づいて一次元信号ｒ（ｔ）の基準時刻を推定し、推定した基準時刻を、ビート信号Ｄ（ｔ）または一次元信号ｒ（ｔ）の少なくともいずれかとともに振動情報特定部１６０に出力する。具体的には、基準時刻推定部１５０は、推定周波数ｆＥ（ｔ）に対応する周期に相当する区間長を有する基準信号を一次元信号ｒ（ｔ）から切り出し、基準信号と推定周波数ｆＥ（ｔ）を有する正弦波との相関係数を算出することにより、一次元信号ｒ（ｔ）の基準時刻を推定する。ここで基準時刻とは、本明細書においては、一次元信号ｒ（ｔ）の一周期ごとの波形について、各波形の類似する波形のうち、波形のピークを示す時刻などの共通する特徴を持った時刻を意味する。また、基準時刻は、例えば波形の極小値を示す位置に相当する時刻でもよいし、ｒ（ｔ）＝０となるときの時刻でもよい。
【００６５】
図６は、本発明の一実施形態に係る基準時刻推定部１５０による基準時刻の推定方法の一例を示す図である。図６に示されるグラフＧ３には、一次元信号ＯＳ１、および正弦波ＲＳ１が示されている。ここで、基準時刻推定部１５０は、時刻ｔｋにおける一次元信号ｒ（ｔｋ）の基準時刻ＲＴ１を推定する。まず、基準時刻推定部１５０は、推定周波数ｆＥ（ｔｋ）に対応する周期ＴＥ（ｔｋ）に相当する区間長を有する基準信号ＢＳ１（破線部）を一次元信号ＯＳ１から切り出す。このとき、区間の終端はｔｋである。次に、基準時刻推定部１５０は、推定周波数ｆＥ（ｔｋ）を有する正弦波ＲＳ１と基準信号ＢＳ１との相関係数を、上記区間において算出する。例えば、基準時刻推定部１５０は、上記区間において正弦波ＲＳ１の時刻ｔｋにおける任意の初期位相をφ０（ｔｋ）とし、時刻ｔｋにおける位相を変化させ、正弦波ＲＳ１と基準信号ＢＳ１の相関係数が最大となるときの正弦波ＲＳ１の位相を特定する。このときの正弦波ＲＳ１の位相は、例えば正弦波ＲＳ１を解析信号に変換することにより、容易に特定可能である。ここで特定される位相をφＲｍａｘ（ｔｋ）とすると、基準時刻ＲＴ１は、ｔｋ−ＴＥ（φ（ｔｋ）／２π）となる。ただし、φ（ｔｋ）＝φＲｍａｘ（ｔｋ）−φ０（ｔｋ）である。例えば、初期位相φ０（ｔｋ）＝０の場合は、グラフＧ３に示されるように、正弦波ＲＳ１のピークにあたる位相を示す時刻が基準時刻ＲＴ１となる。これにより、一次元信号ｒ（ｔ）の基準時刻ＲＴ１を得ることができる。なお、基準時刻推定部１５０は、推定した基準時刻ＲＴ１を不図示の記憶部に記憶してもよい。
【００６６】
また、基準時刻推定部１５０は、基準時刻を推定する際に用いた上記のような基準信号の区間ＴＥにおいて過去に推定された複数の基準時刻を抽出し、抽出された複数の基準時刻の分布に基づいて基準時刻を特定してもよい。例えば、基準時刻推定部１５０は、基準時刻の推定を実施した時刻ｔｋを終点とする区間ＴＥの範囲内において過去に推定された基準時刻を不図示の記憶部から抽出し、各基準時刻の分布を取得してもよい。抽出された各基準時刻の分布は、単位区間で仕切られるヒストグラムにおいて表現されてもよい。
【００６７】
図７は、本発明の一実施形態に係る基準時刻推定部１５０による基準時刻の特定方法の一例を示す図である。図７に示されているグラフＧ４には、一次元信号ＯＳ１が示されている。ここで、時刻ｔｋ−１において推定された基準時刻ＲＴ３が、ｔｋ、ｔｋ−２、およびｔｋ−３において推定された基準時刻ＲＴ２とずれることが生じる。つまり、基準時刻推定部１５０において推定される基準時刻は必ずしも正確ではなく、誤差を含むことも考えられる。そのため、基準時刻推定部１５０は、推定した複数の基準時刻の分布から、基準時刻を統計的に特定してもよい。
【００６８】
具体的には、まず、基準時刻推定部１５０は、ｔｋを終点とする基準信号の区間ＴＥ内において推定された基準時刻の情報を、不図示の記憶部から抽出する。例えば、基準時刻推定部１５０は、ｔｋ−ＴＥからｔｋまでの区間において推定処理が実施された時刻ｔｋ、ｔｋ−１、ｔｋ−２、ｔｋ−３、・・・における各基準時刻の情報を記憶部から抽出し、抽出された各基準時刻に基づくヒストグラムＨ１を作成する。そして、基準時刻推定部１５０は、作成したヒストグラムＨ１をもとに、時刻ｔｋにおける基準時刻を特定してもよい。例えば、基準時刻推定部１５０は、ヒストグラムＨ１の最頻値Ｍ１を示す基準時刻を、時刻ｔｋにおける基準時刻として特定してもよい。また、基準時刻推定部１５０は、区間ＴＥ内において推定された基準時刻の平均値を算出し、平均値を中心とし、任意の分散値を有するガウス関数とヒストグラムＨ１とを乗算することにより得られる新たなヒストグラムの最頻値における基準時刻を、時刻ｔｋにおける基準時刻として特定してもよい。また、基準時刻推定部１５０は、最頻値に限定されず、平均値および中間値などの確率分布を代表する値を用いて、時刻ｔｋにおける基準時刻を特定してもよい。
【００６９】
（振動情報特定部）
振動情報特定部１６０は、推定された基準時刻を用いてビート信号Ｄ（ｔ）（または一次元信号ｒ（ｔ））の振動情報を特定する。ここで、本明細書において振動情報は、例えば、ビート信号Ｄ（ｔ）の周期や振幅を含む。
【００７０】
例えば、振動情報特定部１６０は、基準時刻推定部１５０で推定された連続する２つの基準時刻の間隔をビート信号Ｄ（ｔ）の周期として特定する。これにより、呼吸の振動の周期を、一周期ごとに特定することが可能である。
【００７１】
また、振動情報特定部１６０は、基準時刻推定部１５０で推定された基準時刻における呼吸の振動を表す信号の振幅を、ビート信号Ｄ（ｔ）の振幅として特定してもよい。例えば、振動情報特定部１６０は、基準時刻における二次元ベクトルで表現されたビート信号Ｄ（ｔ）のスカラー量を、ビート信号Ｄ（ｔ）の振幅として特定してもよい。ここで、スカラー量は、例えば、図５で示した中心Ｃ１とビート信号Ｄ（ｔ）までの距離に相当する振幅Ａｍｐ（ｔ）の値であってもよい。
【００７２】
また、振動情報特定部１６０は、基準時刻における一次元信号ｒ（ｔ）の振幅を、ビート信号Ｄ（ｔ）の振幅として特定してもよい。例えば、信号変換部１３０以降の処理においてビート信号Ｄ（ｔ）が保持されず、一次元信号ｒ（ｔ）の情報が保持されている場合、振動情報特定部１６０は、一次元信号ｒ（ｔ）の振幅から、ビート信号Ｄ（ｔ）の振幅を特定してもよい。このとき、振動情報特定部１６０は、一次元信号ｒ（ｔ）をヒルベルト変換により解析信号に変換し、解析信号と、解析信号の共役信号との積により、ビート信号Ｄ（ｔ）の振幅を特定してもよい。これにより、一次元信号ｒ（ｔ）から振幅に係る振動情報を特定することが可能である。また、振動情報特定部１６０は、特定された振幅を対数変換してもよい。これにより、ドップラーセンサと被験者との距離の変化による振幅の変化の影響（Ｓ／Ｎ比や量子化誤差など）を低減させることが可能である。
【００７３】
呼吸推定部１００は、振動情報特定部１６０により特定された振動情報を、睡眠深度推定部３００へ出力する。ここで、特定された振動情報は呼吸一周期ごとに変化する情報であるため、例えば、振動情報は、ドップラーセンサ２のサンプリング周波数（３Ｈｚ以上）よりも低い頻度で睡眠深度推定部３００へ出力されてもよい。本実施形態においては、呼吸推定部１００は、１Ｈｚの頻度で振動情報を睡眠深度推定部３００へ出力する。
【００７４】
以上より、呼吸推定部１００は、ドップラーセンサ２から出力されたビート信号Ｄ（ｔ）から、被験者の呼吸の振動に係る振動情報（周期や振幅）を推定することができる。特に、呼吸推定部１００は、微細な揺らぎを含む呼吸の振動を表す波形の基準位置を推定することにより、呼吸の振動の周期をより正確に推定することが可能である。
【００７５】
［２−２．体動推定部］
体動推定部２００は、ドップラーセンサ２から出力されたビート信号Ｄ（ｔ）から被験者の体動に係る振動情報を推定し、推定した振動情報を睡眠深度推定部３００へ出力する。ここで、体動とは、例えば被験者の寝返りや手足の動きなどを意味する。また、体動に係る振動情報は、例えば、体動の振幅量を含む。
【００７６】
体動推定部２００は、まず、ドップラーセンサ２から出力されたビート信号Ｄ（ｔ）を取得する。なお、体動推定部２００は、呼吸推定部１００においてビート信号Ｄ（ｔ）を取得したビート信号取得部１１０から取得してもよい。
【００７７】
次に、体動推定部２００は、取得したビート信号Ｄ（ｔ）について、ハイパスフィルタやバンドパスフィルタによりフィルタ処理を行う。例えば、ドップラーセンサ２が２４ＧＨｚの放射波を用いて睡眠中の寝返り等の体動を検出する場合、体動推定部２００は、通過帯域が５〜２０Ｈｚのバンドパスフィルタを用いることにより、体動を検出することができる。この場合、体動推定部２００は、バンドパスフィルタによるフィルタ処理がされたビート信号Ｄ（ｔ）のノルムを、体動の振幅量として抽出してもよい。
【００７８】
また、体動推定部２００は、バンドパスフィルタによる処理がされたビート信号Ｄ（ｔ）について、バンドパスフィルタの通過帯域の上限カットオフ周波数よりも低いサブサンプリング処理を実施してもよい。例えば、体動推定部２００は、上限カットオフ周波数が２０Ｈｚであるバンドパスフィルタを用いたフィルタ処理がされたビート信号Ｄ（ｔ）について、５Ｈｚのサブサンプリング処理を実施してもよい。この場合、ナイキスト周波数である２．５Ｈｚよりも大きい周波数成分が２．５Ｈｚ以下の成分として折り返される。しかし、体動推定部２００はビート信号Ｄ（ｔ）の振幅のみを抽出するので、折り返しにより発生する雑音の影響は受けない。そのため、体動推定部２００は、より低いサブサンプリング周波数を用いてサブサンプリング処理を実施することにより、体動の振動情報の推定処理の負荷を低減させることができる。
【００７９】
また、体動推定部２００は、抽出した振幅の単位時間当たりの積算値を振幅量として推定してもよい。また、体動推定部２００は、上記の積算値を対数変換したものを振幅量として推定してもよい。上記の対数変換を実施することにより、体動推定部２００は、微細な体動の差異を強調することができる。信号の強度はセンサと被験者との距離の二乗に反比例して変化するが、対数変換を実施することにより、信号の強度は、上記距離と線形に変化する。これにより、体動推定部２００は、センサと被験者との距離による振幅量への影響を低下させることができる。例えば、８畳間の家屋において睡眠深度推定システム１を適用させる場合、被験者とセンサの距離は約２．５〜５ｍとなるため、対数変換により、信号の強度の減衰率を５０％程度に抑えることが可能である。
【００８０】
体動推定部２００により推定された振幅量を含む振動情報は、睡眠深度推定部３００へ出力される。ここで、推定された振動情報は時間変化する値であるが、例えば、体動推定部２００は、呼吸の振動の周波数よりも高い頻度で振動情報を睡眠深度推定部３００へ出力してもよい。より具体的には、体動推定部２００は、呼吸推定部１００の出力頻度である１Ｈｚの頻度で振動情報を睡眠深度推定部３００へ出力してもよい。このとき、体動推定部２００による出力のタイミングが、呼吸推定部１００による出力のタイミングと同期してもよい。これにより、後述する睡眠深度推定部３００は、同一のタイミングに推定された呼吸および体動の振動情報に基づいて、被験者の睡眠深度を推定することが可能となる。
【００８１】
［２−３．睡眠深度推定部］
睡眠深度推定部３００は、特徴量算出部３１０、モデル選択部３２０、および深度判別部３３０を含む。睡眠深度推定部３００は、呼吸推定部１００から取得した被験者の呼吸の振動情報、または体動推定部２００から取得した被験者の体動の振動情報の少なくともいずれかを用いて、被験者の睡眠深度を推定する。この場合、睡眠深度は所定の区間ごとに推定され、上記の所定の区間の大きさは、特に限定されないが、本実施形態においては、Rechtschaffen，  A．  and  Kales，  A．“A  Manual  of  Standardized  Terminology，  Techniques  and  Scoring  System  for  Sleep  Stage  of  Human  Subjects，  Washington  Public  Health  Service  U.S.  Government  Printing  Office，  Washington，  D.C.(1968)に記載された、国際的な睡眠深度判定基準である「Ｒ＆Ｋ法」にのっとり、３０秒とする。この３０秒の区間をエポックと呼称する。
【００８２】
（睡眠深度について）
まず、睡眠深度推定部３００が推定する睡眠深度について説明する。睡眠深度とは、生体の睡眠の状態により分別される睡眠の深さを意味する。睡眠段階は、例えば、上述したＲ＆Ｋ法によれば、覚醒、レム睡眠、およびノンレム睡眠の状態に区別可能である。さらに、ノンレム睡眠の状態については、例えば、軽睡眠および深睡眠の状態に区別可能である。本実施形態においては、睡眠深度は４段階に分けて推定するが、他の実施形態においては睡眠深度の段階数について特に限定されない。例えば、他の実施形態において睡眠深度は、覚醒と睡眠の２段階であってもよいし、ノンレム睡眠をさらに細分化したものであってもよい。これらの睡眠深度に応じて、睡眠中の被験者の脳の活動状態や自律神経の活性度、および骨格筋の緊張状態もしくは弛緩状態が変化する。これらの変化は、例えば、被験者の体動の大きさ、呼吸の周期の速さや揺らぎ、または呼吸の振幅の揺らぎ等の特徴に現れる。そのため、睡眠深度推定部３００は、上記の特徴について解析を行うことにより、被験者の睡眠深度を推定することができる。以下、睡眠深度推定部３００に含まれる各機能の構成について説明する。
【００８３】
（特徴量算出部）
特徴量算出部３１０は、被験者の睡眠深度を推定するために用いる特徴量を算出する。ここで、特徴量とは例えば以下の特徴量を含み得る。
【００８４】
−体動に関する特徴量ｂｍ
体動に関する特徴量ｂｍは、被験者の寝返りなどの体動の大きさを示す特徴量である。ｂｍは、睡眠または覚醒の状態、および、レム睡眠またはノンレム睡眠の状態を判別するために用いられる。例えば、ｂｍの値が低い場合は、体動が小さいので、被験者は睡眠状態であると考えられる。特徴量算出部３１０は、１秒ごとに１つ得られる体動の振動情報に含まれる振幅量から、ｂｍを算出する。具体的には、特徴量算出部３１０は、１エポックに含まれる３０個の振幅量の平均値を当該エポックにおけるｂｍとして算出してもよい。
【００８５】
−呼吸の振幅に関する特徴量ｃｖ
呼吸の振幅に関する特徴量ｃｖは、呼吸における吸い込みや吐き出しの大きさの揺らぎを示す特徴量である。ｃｖは、睡眠または覚醒の状態、および、軽睡眠または深睡眠の状態を判別するために用いられる。例えば、ｃｖの値が低い場合は、呼吸が安定であるので、被験者は睡眠状態であると考えられる。特徴量算出部３１０は、１秒ごとに１つ得られる呼吸の振動情報に含まれる振幅から、ｃｖを算出する。具体的には、特徴量算出部３１０は、１エポックに含まれる３０個の振幅の平均値及び標準偏差を算出し、標準偏差を平均値で割ることにより得られる値をｃｖとして算出してもよい。
【００８６】
−呼吸の周期の長さに関する特徴量ｒｅｓｐ
呼吸の周期の長さに関する特徴量ｒｅｓｐは、呼吸の周期の長短を示す特徴量である。ｒｅｓｐは、軽睡眠または深睡眠の状態を判別するために用いられる。例えば、ｒｅｓｐの値が低い場合は、呼吸の周期が短いので、被験者はノンレム睡眠のうち睡眠が浅いとされる軽睡眠であると考えられる。特徴量算出部３１０は、１秒ごとに１つ得られる呼吸の振動情報に含まれる周期から、ｒｅｓｐを算出する。具体的には、特徴量算出部３１０は、１エポックに含まれる３０個の周期の平均値を算出し、平均値を係数αで割ることにより得られる値をｒｅｓｐとして算出してもよい。ここで、係数αは、被験者の呼吸の周期に合わせて変動する値でもよいし、固定値であってもよい。例えば、係数αは、被験者の睡眠時における呼吸の周波数の平均値であってもよい。これにより、個人差によるｒｅｓｐの値のばらつきを抑制することが可能である。また、係数αを固定値として用いる場合、係数αは、０．２〜０．５程度の値としてもよい。
【００８７】
−呼吸の周期の揺らぎに関する特徴量ｒｅｓｐ＿ａｖｅ
呼吸の周期の揺らぎに関する特徴量ｒｅｓｐ＿ａｖｅは、呼吸の周期の揺らぎの大きさを示す特徴量である。ｒｅｓｐ＿ａｖｅは、レム睡眠またはノンレム睡眠の状態を判別するために用いられる。例えば、ｒｅｓｐ＿ａｖｅの値が低い場合は、呼吸の周期の揺らぎが小さいので、被験者はノンレム睡眠の状態であると考えられる。特徴量算出部３１０は、あるエポックにおいて算出されたｒｅｓｐと１エポック前において算出されたｒｅｓｐとを比較することにより算出する。具体的には、まず、特徴量算出部３１０は、あるエポックにおいて算出されたｒｅｓｐと１エポック前において算出されたｒｅｓｐとの差分の二乗を算出し、連続する所定数のエポックにおいて算出された上記の差分の二乗の平均値をｒｅｓｐ＿ａｖｅとして算出してもよい。上記の所定数は、過去数分程度の呼吸の周期の揺らぎを捉えるために、例えば１０であってもよい。
【００８８】
以上の４つの特徴量は、深度判別部３３０に出力される。なお、例えば、ｂｍ、ｃｖ、ｒｅｓｐ、およびｒｅｓｐ＿ａｖｅについて、算術平均を用いて平均値が算出されているが、例えば、算出される平均値は、移動平均や指数平滑化平均であってもよい。また、平均値の代わりに、中央値や最頻値などが用いられてもよい。また、ｃｖの算出に用いられる標準偏差の代わりに、振幅値の最大および最小の差、振幅値の第３四分位点と第１四分位点の差などが用いられてもよい。また、ｒｅｓｐは、１エポックに含まれる複数の周期について、周波数スペクトルのパワーや振幅を重みとした加重平均である平均パワー周波数の逆数による求まる値であってもよい。また、ｒｅｓｐ＿ａｖｅは、１エポックに含まれる複数の周期について、周波数スペクトルのパワーや振幅を重みとした加重分散から求まる値であってもよい。
【００８９】
一方、以下で説明する２つの特徴量は、深度判別部３３０において睡眠深度を判別するために用いられるモデルをモデル選択部３２０において決定するために用いられる。ここで、レム睡眠の状態、軽睡眠の状態、および深睡眠の状態を、睡眠状態と総称する。
【００９０】
−睡眠継続時間ｓｌｐ＿ｅｔ
睡眠継続時間ｓｌｐ＿ｅｔは、過去の所定期間において、被験者の睡眠状態の継続時間を示す特徴量である。ｓｌｐ＿ｅｔは、覚醒状態（または被験者の睡眠深度を推定していない状態）から睡眠状態であると推定された時点を起点とする、睡眠状態が継続している時間を指す。なお、所定期間とは、被験者の睡眠状態に基づいて決定してもよいし、固定値であってもよい。例えば、所定期間は、４時間＝１４４００秒であってもよい。ここで、例えば睡眠状態が途切れて覚醒状態に変化した場合、ｓｌｐ＿ｅｔは、上記起点から覚醒状態に変化した時刻までの時間としてもよい。また、覚醒状態に変化した後に第１の所定時間以内に再び睡眠状態に変化した場合、覚醒状態であった時間は、睡眠継続時間としてｓｌｐ＿ｅｔに加算されてもよい。第１の所定時間は特に限定されないが、覚醒状態が直後の睡眠状態に影響を及ぼさない程度の時間であることが好ましく、例えば、１時間であってもよい。また、覚醒状態から睡眠状態に変化した後に第２の所定時間以内に再び覚醒状態に変化した場合、ｓｌｐ＿ｅｔの値は０にリセットされてもよい。第２の所定時間は特に限定されないが、一時的な眠気による睡眠状態と判別するのに必要な時間であることが望ましく、例えば、１０分であってもよい。第１の所定時間および第２の所定時間は、被験者の睡眠の特性や被験者の健康状態等に応じて、自由に変更可能である。
【００９１】
−睡眠深度積算時間ｓｕｍ＿ｒｅｍ、ｓｕｍ＿ｌｉｇｈｔ、ｓｕｍ＿ｄｅｅｐ
睡眠深度積算時間ｓｕｍ＿ｒｅｍ、ｓｕｍ＿ｌｉｇｈｔ、およびｓｕｍ＿ｄｅｅｐ（以下、３つの睡眠深度積算時間を総称する場合は、ｓｕｍｓと呼称する）は、過去の所定期間において発生した、被験者のレム睡眠、軽睡眠、および深睡眠の状態の各々の積算時間を示す特徴量である。ｓｕｍｓの各々の値は、１エポック単位で判定される各睡眠状態に基づいて決定されてもよい。例えば、あるエポックにおいてレム睡眠であると推定された場合、ｓｕｍ＿ｒｅｍに３０秒が加算される。
【００９２】
特徴量算出部３１０は、深度判別部３３０による睡眠深度の特定結果のフィードバックを受けることにより、ｓｌｐ＿ｅｔおよびｓｕｍｓを更新する。そして、特徴量算出部３１０は、更新されたｓｌｐ＿ｅｔまたはｓｕｍｓをモデル選択部３２０へ出力する。なお、詳しくは後述するが、深度判別部３３０による睡眠深度の特定処理を実施する前段階として深度判別部３３０において用いられるモデルについて学習を行う場合は、特徴量算出部３１０は、ＰＳＧ（ＰｏｌｙＳｏｍｎｏＧｒａｐｈｙ：睡眠ポリグラフ検査）などの他の睡眠深度測定手法により得られる睡眠深度の正解値に基づいてｓｌｐ＿ｅｔおよびｓｕｍｓを更新してもよい。更新されたｓｌｐ＿ｅｔまたはｓｕｍｓがモデル選択部３２０へ出力されることにより、ｓｌｐ＿ｅｔまたはｓｕｍｓに基づいて選択された他のモデルについての学習も行うことができる。また、ｓｌｐ＿ｅｔおよびｓｕｍｓは、睡眠深度の特定結果のフィードバックと、睡眠深度の正解値との二つの手法を並行して用いることにより更新されてもよい。この場合、二つの手法により更新されたｓｌｐ＿ｅｔおよびｓｕｍｓの値が二つの手法の間で乖離している場合、一方の手法により得られた値を他方の手法により得られた値に反映させてもよい。また、二つの手法の間において得られた値が乖離している場合、その二つの値の差分をモデル選択部３２０へ反映させてもよい。モデル選択部３２０は、上記の差分を用いて、後述する判別モデルに含まれる各パラメータの値を最適化することができる。
【００９３】
なお、ｓｌｐ＿ｅｔおよびｓｕｍｓの単位は特に限定されないが、本実施形態においては、ｓｌｐ＿ｅｔおよびｓｕｍｓの単位は、秒とする。
【００９４】
（モデル選択部）
モデル選択部３２０は、被験者の睡眠の継続時間や各睡眠深度における積算時間に応じて、睡眠深度の推定条件を変更する。例えば、モデル選択部３２０は、深度判別部３３０による睡眠深度の特定処理に用いられる判別モデルを、ｓｌｐ＿ｅｔまたはｓｕｍｓの少なくともいずれかに基づいて決定してもよい。判別モデルは、ｂｍ、ｃｖ、ｒｅｓｐ、およびｒｅｓｐ＿ａｖｅを独立変数として入力し、各段階の睡眠深度を推定するための目的変数を得るのに用いられる判別行列Ｍおよび補正ベクトルｂを有する。判別行列Ｍおよび補正ベクトルｂのパラメータは、例えば、上述したＰＳＧなどの他の睡眠深度測定手法により得られる睡眠深度の正解値を教師データとする教師付学習や半教師付学習等の手法を用いて決定されてもよい。また、判別行列Ｍおよび補正ベクトルｂのパラメータは、教師なし学習に基づくクラス分けが実施された後に、教師付学習により最適化されてもよい。各パラメータの最適化に用いられる教師付学習のアルゴリズムは、公知のアルゴリズムであってもよい。
【００９５】
ここで、睡眠の深さに基づく呼吸や体動の変化の傾向は、睡眠継続時間に依存すると言われている。そのため、モデル選択部３２０は、睡眠継続時間を示すｓｌｐ＿ｅｔや、各睡眠深度における睡眠深度積算時間ｓｕｍｓに基づいて、最適と考えられる判別モデルを選択し、深度判別部３３０に出力する。例えば、モデル選択部３２０は、睡眠の初期、中期および後記に合わせて、３種類の判別モデルを備えてもよい。具体的には、睡眠の初期に用いられる判別モデル、睡眠の中期に用いられる判別モデル、および睡眠の後期に用いられる判別モデルのインデックスをそれぞれ１、２、および３とすると、モデル１、モデル２、およびモデル３において用いられる判別行列Ｍおよび補正ベクトルｂは、以下の数式４〜６のように表される。
【００９６】
【数４】
0006515670.tif 000005
【００９７】
【数５】
0006515670.tif 000006
【００９８】
【数６】
0006515670.tif 000007
【００９９】
モデル選択部３２０は、例えば、ｓｌｐ＿ｅｔと所定の閾値を比較することにより判別モデルを選択してもよい。例えば、ｓｌｐ＿ｅｔ＜７０００の場合、モデル選択部３２０はモデル１を選択し、７０００≦ｓｌｐ＿ｅｔ＜１４０００の場合、モデル選択部３２０はモデル２を選択し、１４０００≦ｓｌｐ＿ｅｔの場合、モデル選択部３２０はモデル３を選択してもよい。これにより、睡眠継続時間に応じて、判別モデルを選択することが可能である。
【０１００】
また、モデル選択部３２０は、ｓｕｍｓと所定の閾値を比較することにより判別モデルを選択してもよい。例えば、ｓｕｍ＿ｒｅｍ＋ｓｕｍ＿ｌｉｇｈｔ＋ｓｕｍ＿ｄｅｅｐ＜７０００である場合、モデル選択部３２０はモデル１を選択してもよい。また、ｓｕｍ＿ｒｅｍ＋ｓｕｍ＿ｌｉｇｈｔ＋ｓｕｍ＿ｄｅｅｐ≧７０００であって、（ｓｕｍ＿ｒｅｍ＋ｓｕｍ＿ｌｉｇｈｔ）×０．７５＜ｓｕｍ＿ｄｅｅｐである場合、モデル選択部３２０はモデル２を選択してもよい。また、ｓｕｍ＿ｒｅｍ＋ｓｕｍ＿ｌｉｇｈｔ＋ｓｕｍ＿ｄｅｅｐ≧７０００であって、（ｓｕｍ＿ｒｅｍ＋ｓｕｍ＿ｌｉｇｈｔ）×０．７５≧ｓｕｍ＿ｄｅｅｐである場合、モデル選択部３２０はモデル３を選択してもよい。これにより、各睡眠深度の積算時間に応じて、判別モデルを選択することが可能である。つまり、被験者の睡眠の深さに応じて判別モデルを選択することが可能となる。
【０１０１】
モデル選択部３２０は、選択したモデルを深度判別部３３０に出力する。なお、数式４〜６に示した各判別モデルの判別行列Ｍおよび補正ベクトルｂに含まれるパラメータの値は、あくまで本実施形態において用いられる値の一例であり、パラメータの値は、上述したように、教師付学習等のアルゴリズムを用いて適宜変更され得る。なお、本実施形態においては、判別モデルの数は３としたが、他の実施形態において判別モデルの数は特に限定されない。また、本実施形態においては、判別モデルは睡眠継続時間や各睡眠深度の積算時間に応じて選択されたが、例えば、被験者の健康状態、疾病の種類、年齢層、または性別等に応じてモデルが構築されてもよい。また、判別モデルのパラメータは、被験者ごとに決定されてもよい。また、本実施形態においては、モデル選択部３２０は、ｓｌｐ＿ｅｔまたはｓｕｍｓのいずれか１つを用いて判別モデルを選択したが、ｓｌｐ＿ｅｔとｓｕｍｓの双方に基づいて判別モデルを選択してもよい。例えば、モデル選択部３２０は、ニューラルネットワーク等の機械学習を用いて、選択すべき判別モデルを決定してもよい。
【０１０２】
（深度判別部）
深度判別部３３０は、特徴量算出部３１０から出力された特徴量、およびモデル選択部３２０により選択された判別モデルを用いて、被験者の睡眠深度を判別する。例えば、深度判別部３３０は、ｂｍ、ｃｖ、ｒｅｓｐ、およびｒｅｓｐ＿ａｖｅの各特徴量を独立変数として、判別モデルに含まれる判別行列Ｍおよび補正ベクトルｂを用いて、線形判別分析を行うことにより、被験者の睡眠深度を判別してもよい。
【０１０３】
まず、深度判別部３３０は、上記の４つの特徴量を含む特徴ベクトルｖｔを、以下の数式７のように定義する。
【０１０４】
【数７】
0006515670.tif 000008
【０１０５】
ただし、添え字ｔは、時刻を表す。また、睡眠深度を判別するのに用いる判別ベクトルｄｔを、以下の数式８のように定義する。
【０１０６】
【数８】
0006515670.tif 000009
【０１０７】
ここで、ＷｔとＳｔは、覚醒状態と睡眠状態との判別に寄与する変数であり、ＲｔとＮＲｔは、レム睡眠とノンレム睡眠との判別に寄与する変数であり、ＤｔとＬｔは、深睡眠と軽睡眠との判別に寄与する変数である。モデル選択部３２０において選択されたモデルのインデックスを示す数字をｍｏｄｅｌ＿ｉｄｘとすると、判別ベクトルｄｔは、以下の数式９を用いて算出される。
【０１０８】
【数９】
0006515670.tif 000010
【０１０９】
次に、深度判別部３３０は、算出された判別ベクトルｄｔに基づいて、被験者の睡眠深度を判別する。具体的には、深度判別部３３０は、判別ベクトルｄｔに含まれる各変数の値を比較することにより、被験者の４段階の睡眠深度を判別する。
【０１１０】
図８は、本発明の一実施形態に係る深度判別部３３０による睡眠深度の判別フローを示すフローチャートの一例を示す図である。まず、深度判別部３３０は、被験者が覚醒しているか、または睡眠しているかを判別する（Ｓ４０１）。具体的には、深度判別部３３０は、覚醒状態に寄与する変数Ｗｔと、睡眠状態に寄与する変数Ｓｔとを比較することにより、被験者が覚醒状態であるか、または睡眠状態であるかを判別する。ここで、ＷｔおよびＳｔは、体動に関する特徴量ｂｍ、および呼吸の振幅に関する特徴量ｃｖを用いて算出される変数である。例えば、覚醒状態は、睡眠状態と明確に被験者の状態が異なるため、覚醒状態においては、被験者の体動や呼吸が大きくなることが知られている。そのため、深度判別部３３０は、体動の振幅が大きく、呼吸の振幅の揺らぎが大きい場合、被験者は覚醒状態であると判別することができる。ここで、Ｗｔ＞Ｓｔである場合（ＹＥＳ）、深度判別部３３０は、被験者は覚醒状態であると判別する（Ｓ４０２）。一方、ステップＳ４０１においてＷｔ≦Ｓｔである場合（ＮＯ）、深度判別部３３０は、被験者が睡眠状態であるとして、次のステップＳ４０３へ進む。
【０１１１】
深度判別部３３０は、次に、被験者がレム睡眠状態であるか、またはノンレム睡眠状態であるかを判別する（Ｓ４０３）。具体的には、深度判別部３３０は、レム睡眠状態に寄与する変数Ｒｔと、ノンレム睡眠状態に寄与する変数ＮＲｔとを比較することにより、被験者がレム睡眠状態であるか、またはノンレム睡眠状態であるかを判別する。ここで、ＲｔおよびＮＲｔは、体動に関する特徴量ｂｍ、および呼吸の周期の揺らぎに関する特徴量ｒｅｓｐ＿ａｖｅを用いて算出される変数である。例えば、レム睡眠状態は、被験者の自律神経がノンレム睡眠状態よりも乱れているため、レム睡眠状態においては、被験者の体動が大きくなったり、呼吸間隔が乱れるということが知られている。そのため、深度判別部３３０は、体動の振幅および呼吸の周期の揺らぎ大きい場合、被験者はレム睡眠状態であると判別することができる。ここで、Ｒｔ＞ＮＲｔである場合（ＹＥＳ）、深度判別部３３０は、被験者はレム睡眠状態であると判別する（Ｓ４０４）。一方、ステップＳ４０３においてＲｔ≦ＮＲｔである場合（ＮＯ）、深度判別部３３０は、被験者がノンレム睡眠状態であるとして、次のステップＳ４０５へ進む。
【０１１２】
深度判別部３３０は、次に、被験者が軽睡眠状態であるか、または深睡眠状態であるかを判別する（Ｓ４０５）。具体的には、深度判別部３３０は、軽睡眠状態に寄与する変数Ｌｔと、深睡眠状態に寄与するＤｔとを比較することにより、被験者が軽睡眠状態であるか、または深睡眠状態であるかを判別する。ここで、ＬｔおよびＤｔは、呼吸の振幅に関する特徴量ｃｖ、および呼吸の周期の長さに関する特徴量ｒｅｓｐを用いて算出される変数である。例えば、深睡眠においては呼吸が安定しているため、呼吸の周期および振幅の揺らぎは小さくなる。そのため、深度判別部３３０は、呼吸の周期および振幅の揺らぎが小さい場合、被験者は深睡眠状態であると判別することができる。ここで、Ｌｔ＞Ｄｔである場合（ＹＥＳ）、深度判別部３３０は、被験者は軽睡眠状態であると判別する（Ｓ４０６）。一方、ステップＳ４０５においてＬｔ≦Ｄｔである場合（ＮＯ）、深度判別部３３０は、被験者が深睡眠状態であると判別する（Ｓ４０７）。以上のように、深度判別部３３０は、特徴量算出部３１０において算出された、被験者の呼吸や体動から得られる特徴量、およびモデル選択部３２０において選択されたモデルを用いることにより、被験者の睡眠状態を時系列に判別することができる。
【０１１３】
図９は、睡眠深度の推定結果を出力したグラフの一例を示す図である。図９には、睡眠深度の推定結果の時系列変化を示すグラフＰｈ１が示されている。グラフＰｈ１は被験者の睡眠深度の変化を、覚醒状態、レム睡眠状態、軽睡眠状態、および深睡眠状態の４段階において示している。睡眠深度推定部３００は、深度判別部３３０において推定された被験者の睡眠状態について、例えば１エポック（３０秒）ごとに出力してもよい。また、睡眠深度推定部３００は、各エポックにおいて推定された被験者の睡眠状態を不図示の記憶部に記憶させ、事後的に出力してもよい。
【０１１４】
なお、本実施形態において睡眠深度推定部３００は、４段階の睡眠深度を推定するが、睡眠深度の段階数については特に限定されない。例えば、睡眠深度推定部３００は、覚醒状態と睡眠状態の２段階の睡眠深度を推定してもよいし、覚醒状態、レム睡眠状態、および４つのノンレム睡眠状態の、計６段階の睡眠深度を推定してもよい。この場合、睡眠深度推定部３００は、段階数等に応じて、睡眠深度を推定するのに用いる判別ベクトルの変数の数を設定してもよい。これにより、多段階の睡眠深度を推定することが可能である。
【０１１５】
以上のような睡眠深度推定装置１０の構成によって、微細な揺らぎを伴う被験者の呼吸および体動の振動情報をドップラーセンサ２により取得されたビート信号Ｄ（ｔ）から推定し、さらに得られた振動情報をもとに、被験者の睡眠深度を推定することが可能である。特に、睡眠深度推定装置１０は、例えば、レム睡眠状態等の睡眠状態にみられる呼吸の乱れにより変化する呼吸の周期の揺らぎを推定できるため、被験者の睡眠深度をより精度高く推定することが可能となる。
【０１１６】
＜３．睡眠深度推定装置の動作例＞
次に、本発明の一実施形態に係る睡眠深度推定装置１０の動作例について説明する。図１０は、本発明の一実施形態に係る睡眠深度推定装置１０の動作例を示すフローチャートである。
【０１１７】
まず、呼吸推定部１００のビート信号取得部１１０、および体動推定部２００は、ドップラーセンサ２から出力されたビート信号Ｄ（ｔ）を取得する（Ｓ５０１）。呼吸推定部１００は、取得したビート信号Ｄ（ｔ）から、被験者の呼吸振幅や呼吸周期などの振動情報を推定する（Ｓ５０２）。また、体動推定部２００は、取得したビート信号Ｄ（ｔ）から、被験者の振幅量などの体動に関する振動情報を推定する（Ｓ５０３）。ステップＳ５０１〜Ｓ５０３については、呼吸および体動に関する振動情報の推定処理を継続する場合、これらの処理は繰り返し実施される。
【０１１８】
次に、特徴量算出部３１０は、呼吸および体動に関する振動情報から、睡眠深度を推定するのに用いる特徴量を算出する（Ｓ５０４）。続いて、モデル選択部３２０は、睡眠継続時間および各睡眠深度の積算時間に基づいて、睡眠深度を推定するのに用いる判別モデルを選択する（Ｓ５０５）。そして、深度推定部は、上記の特徴量および判別モデルを用いて、被験者の睡眠深度を判別する（Ｓ５０６）。
【０１１９】
＜４．ハードウェア構成例＞
以上、本発明の一実施形態に係る睡眠深度推定装置１０の動作例について説明した。上述した睡眠深度推定装置１０の情報処理は、ソフトウェアと、睡眠深度推定装置１０との協働により実現される。以下では、本発明の実施形態に係る睡眠深度推定装置１０のハードウェア構成を説明する。
【０１２０】
図１１は、本発明の実施形態に係る睡眠深度推定装置１０のハードウェア構成例を示すブロック図である。図１１を参照すると、睡眠深度推定装置１０は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）９０１と、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）９０２と、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）９０３と、ホストバス９０４と、を備える。また、睡眠深度推定装置１０は、ブリッジ９０５と、外部バス９０６と、インタフェース９０７と、入力装置９０８と、出力装置９０９と、ストレージ装置９１０と、ドライブ９１１と、ネットワークインタフェース９１２と、を備える。
【０１２１】
ＣＰＵ９０１は、演算処理装置および制御装置として機能し、各種プログラムに従って睡眠深度推定装置１０内の動作全般を制御する。また、ＣＰＵ９０１は、マイクロプロセッサであってもよい。なお、ＣＰＵ９０１は、睡眠深度推定装置１０内の動作全般またはその一部を制御する。例えば、ＣＰＵ９０１は、ＲＯＭ９０２は、ＣＰＵ９０１が使用するプログラムや演算パラメータ等を記憶する。ＲＡＭ９０３は、ＣＰＵ９０１の実行において使用するプログラムや、その実行において適宜変化するパラメータ等を一時記憶する。これらはＣＰＵバスなどから構成されるホストバス９０４により相互に接続されている。
【０１２２】
ホストバス９０４は、ブリッジ９０５を介して、ＰＣＩ（Ｐｅｒｉｐｈｅｒａｌ  Ｃｏｍｐｏｎｅｎｔ  Ｉｎｔｅｒｃｏｎｎｅｃｔ／Ｉｎｔｅｒｆａｃｅ）バスなどの外部バス９０６に接続されている。なお、必ずしもホストバス９０４、ブリッジ９０５および外部バス９０６を分離構成する必要はなく、１つのバスにこれらの機能を実装してもよい。
【０１２３】
入力装置９０８は、マウス、キーボード、タッチパネル、ボタン、マイクロフォン、スイッチおよびレバーなどユーザが情報を入力するための入力手段と、ユーザによる入力に基づいて入力信号を生成し、ＣＰＵ９０１に出力する入力制御回路などから構成されている。睡眠深度推定装置１０のユーザは、当該入力装置９０８を操作することにより、睡眠深度推定装置１０に対して各種のデータを入力したり処理動作を指示したりすることができる。
【０１２４】
出力装置９０９は、例えば、ＣＲＴディスプレイ装置、液晶ディスプレイ（ＬＣＤ）装置、ＯＬＥＤ装置およびランプなどの表示装置を含む。さらに、出力装置９０９は、スピーカ及びヘッドホンなどの音声出力装置を含む。出力装置９０９は、例えば、再生されたコンテンツを出力する。具体的には、表示装置は再生された映像データ等の各種情報をテキストまたはイメージで表示する。一方、音声出力装置は、再生された音声データや表示装置に表示されたテキストデータ等を音声に変換して出力する。
【０１２５】
ストレージ装置９１０は、本発明の実施形態に係る睡眠深度推定装置１０におけるデータ格納用の装置である。ストレージ装置９１０は、記憶媒体、記憶媒体にデータを記録する記録装置、記憶媒体からデータを読み出す読み出し装置および記憶媒体に記憶されたデータを削除する削除装置などを含んでも良い。ストレージ装置は、例えば、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｃ  Ｄｒｉｖｅ）やＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）で構成される。このストレージ装置９１０は、ＣＰＵ９０１が実行するプログラムや各種データを格納する。なお、ストレージ装置９１０は、不図示の記憶部の機能を実現する。
【０１２６】
ドライブ９１１は、記憶媒体用リーダライタであり、睡眠深度推定装置１０に内蔵、あるいは外付けされる。ドライブ９１１は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記憶媒体９６に記録されている情報を読みだして、ＲＡＭ９０３に出力する。また、ドライブ９１１は、リムーバブル記憶媒体９６に情報を書き込むこともできる。
【０１２７】
ネットワークインタフェース９１２は、例えば、他の装置に接続するための通信デバイス等で構成された通信インタフェースである。また、ネットワークインタフェース９１２は、無線ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）対応通信装置、ＬＴＥ（Ｌｏｎｇ  Ｔｅｒｍ  Ｅｖｏｌｕｔｉｏｎ）対応通信装置、またはブルートゥース通信装置であってもよい。また、ネットワークインタフェース９１２は、有線による通信を行うワイヤー通信装置であってもよい。なお、ネットワークインタフェース９１２は、不図示の通信部の機能を実現する。
【０１２８】
以上、睡眠深度推定装置１０のハードウェア構成の一例を示した。上記の各構成要素は、汎用的な部材を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。かかる構成は、実施する時々の技術レベルに応じて適宜変更され得る。
【０１２９】
＜５．まとめ＞
ここまで、図１〜図１１を用いて、本発明の一実施形態について説明した。本発明の一実施形態によれば、睡眠深度推定装置１０は、取得したビート信号から被験者の呼吸の振動情報、および体動の振動情報を推定し、推定された各々の振動情報を用いて被験者の睡眠深度を推定する。特に、睡眠深度推定装置１０は、被験者の呼吸に含まれる微細な揺らぎを伴う振動について、正弦波との相関を算出することにより、呼吸の振動の周期を精度高く推定することが可能である。これにより、睡眠深度推定装置１０は、呼吸の周期の乱れを生じさせるレム睡眠など、呼吸の振動について異なる特徴を有する各段階の睡眠深度を推定することが可能となる。
【０１３０】
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。
【０１３１】
例えば、上記実施形態では、睡眠深度推定装置１０は、ドップラーセンサ２から出力されるビート信号をリアルタイムで取得し、並行して解析する動作を実施していたが、本発明はかかる例に限定されない。例えば、他の実施形態に係る睡眠深度推定装置は、一度取得した信号を記憶部に保存し、推定処理を実施したい際に記憶部から再度信号を取り出して解析を行ってもよい。これにより、過去に取得したビート信号についてまとめて解析を行うことが可能となる。
【０１３２】
また、上記実施形態では、被験者がドップラーセンサ２による検知が可能なエリアに存在していることが想定されていたが、本発明はかかる例に限定されない。例えば、被験者が椅子やベッドなどの位置にいる場合にのみドップラーセンサが被験者を検知することができるように、ドップラーセンサが設置されてもよい。これにより、解析対象となる被験者を限定することが可能となる。また、睡眠深度推定装置は、被験者の有無や被験者による歩行などの大きな活動を検知する手段を組み合わせてもよい。これにより、被験者が睡眠しているなど限られた状態にある場合にのみ推定処理を実施することが可能となる。
【０１３３】
また、睡眠深度推定装置１０および睡眠深度推定システム１は、他の装置やシステム等に組み込まれて構成されてもよい。例えば、睡眠深度推定装置またはドップラーセンサは、テレビ等の家電機器や空調機、医療介護機器、およびベッドなどの寝具に組み込まれてもよい。さらに、睡眠深度推定装置は、推定した睡眠深度に係る情報を、上記のような他の装置の制御のために用いてもよい。例えば、睡眠深度推定装置は被験者の睡眠深度に応じて、空調機の空調の温度や風量などを調節してもよいし、テレビの電源を制御してもよい。
【０１３４】
また、上記実施形態では、睡眠深度推定装置１０は被験者の睡眠深度を推定したが、本発明は係る例に限定されない。例えば、睡眠深度推定装置は、被験者の心的ストレスを推定してもよい。この場合、例えば他の手段を用いて推定された被験者の心的ストレスを教師データとして判別学習を実施することにより、睡眠深度推定装置はビート信号を用いて被験者の心的ストレスを推定することが可能となる。なお、上述した他の手段とは、例えば、唾液に含まれるアミラーゼの検査や、心拍の周波比（高周波成分および低周波成分の比率）を用いる手段が含まれる。また、睡眠深度推定装置は、人間のみならず、動物などの生体に関する睡眠深度を推定してもよい。
【０１３５】
また、上記実施形態では、ドップラーセンサ２と睡眠深度推定装置１０とが睡眠深度推定システム１を構成していたが、本発明は係る例に限定されない。例えば、ドップラーセンサ２と睡眠深度推定装置１０とが一体となった睡眠深度推定装置または睡眠深度推定システムが提供されてもよい。
【０１３６】
また、本明細書の睡眠深度推定装置１０の処理における各ステップは、必ずしもフローチャートとして記載された順序に沿って時系列に処理する必要はない。例えば、睡眠深度推定装置１０の処理における各ステップは、フローチャートとして記載した順序と異なる順序で処理されても、並列的に処理されてもよい。
【０１３７】
また、睡眠深度推定装置１０に内蔵されるＣＰＵ、ＲＯＭおよびＲＡＭなどのハードウェアを、上記の睡眠深度推定装置１０の各構成と同等の機能を発揮させるためのコンピュータプログラムも作成可能である。また、該コンピュータプログラムを記憶させた記憶媒体も提供される。
【符号の説明】
【０１３８】
１      睡眠深度推定システム
２      ドップラーセンサ
１０    睡眠深度推定装置
１００  呼吸推定部
１１０  ビート信号取得部
１２０  フィルタ部
１３０  信号変換部
１４０  周波数推定部
１５０  基準時刻推定部
１６０  振動情報特定部
２００  体動推定部
３００  睡眠深度推定部
３１０  特徴量算出部
３２０  モデル選択部
３３０  深度判別部
【図１】
0006515670.tif 000011
【図２】
0006515670.tif 000012
【図３】
0006515670.tif 000013
【図４】
0006515670.tif 000014
【図５】
0006515670.tif 000015
【図６】
0006515670.tif 000016
【図７】
0006515670.tif 000017
【図８】
0006515670.tif 000018
【図９】
0006515670.tif 000019
【図１０】
0006515670.tif 000020
【図１１】
0006515670.tif 000021