(57)【要約】【課題】被験者の心拍情報に基づいて精度よく睡眠段階を判定できる睡眠判定装置を提供すること。【解決手段】被験者の心拍を取得し（Ｓ１０）、取得された心拍に基づいて心拍ごとの心拍周期を演算し（Ｓ１２）、連続して演算された二つの値を座標値とする座標点を演算してリターンマップを生成し（Ｓ１４）、そのリターンマップの特徴量を演算し（Ｓ１６）、特徴量に基づいて被験者の睡眠段階を判定する（Ｓ１８）。これにより、睡眠段階に応じて変化する心拍周期の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。【選択図】図４
2010148605.tif 000002
【特許請求の範囲】
【請求項１】
被験者の心拍を取得する心拍取得手段と、
前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期、その心拍周期の前後の差分値、又はその差分値の前後の差分値について、少なくとも一つを演算する差分演算手段と、
前記差分演算手段により演算された心拍周期及び差分値について、連続して演算された二つの値を座標値とする座標点を繰り返し演算してリターンマップを生成するリターンマップ生成手段と、
前記リターンマップ生成手段により生成されたリターンマップの特徴量を演算する特徴量演算手段と、
前記特徴量演算手段により演算された特徴量に基づいて前記被験者の睡眠段階を判定する睡眠判定手段と、
を備えた睡眠判定装置。
【請求項２】
前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、
前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、
前記特徴量演算手段は、前記リターンマップ生成手段により演算される座標点の変動量を前記特徴量として演算し、
前記睡眠判定手段は、前記特徴量演算手段により演算された前記座標点の変動量のバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定する、
請求項１に記載の睡眠判定装置。
【請求項３】
前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、
前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、
前記特徴量演算手段は、前記リターンマップ生成手段により演算される座標点の変動量を前記特徴量として演算し、
前記睡眠判定手段は、前記特徴量演算手段により演算された前記座標点の変動量のヒストグラムにおける歪度が小さいほど睡眠状態が深いとして睡眠段階を判定する、
請求項１に記載の睡眠判定装置。
【請求項４】
前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、
前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、
前記特徴量演算手段は、前記リターンマップ生成手段により演算される複数の座標点における集積度を前記特徴量として演算し、
前記睡眠判定手段は、前記特徴量演算手段により演算された前記座標点の集積度が高いほど睡眠状態が深いとして睡眠段階を判定する、
請求項１に記載の睡眠判定装置。
【請求項５】
前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、
前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、
前記特徴量演算手段は、前記リターンマップ生成手段により連続して演算される三つの座標点により形成される三角形の重心位置の変動量を前記特徴量として繰り返し演算し、
前記睡眠判定手段は、前記特徴量演算手段により演算された前記三角形の重心位置の変動量が小さいほど睡眠状態が深いとして睡眠段階を判定する、
請求項１に記載の睡眠判定装置。
【請求項６】
前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期の前後の差分値を演算し、
前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの差分値を座標値とする座標点を繰り返し演算し、
前記特徴量演算手段は、前記リターンマップ生成手段により連続して演算される三つの座標点により形成される三角形の頂点の変動量を前記特徴量として演算し、
前記睡眠判定手段は、前記特徴量演算手段により演算された前記変動量が小さいほど睡眠状態が深いとして睡眠段階を判定する、
請求項１に記載の睡眠判定装置。
【請求項７】
前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期の前後の差分値を演算し、
前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの差分値を座標値とする座標点を繰り返し演算し、
前記特徴量演算手段は、前記リターンマップ生成手段により連続して演算される三つの座標点により形成される三角形の辺の傾きを前記特徴量として演算し、
前記睡眠判定手段は、前記特徴量演算手段により演算された前記三角形の辺の傾きのバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定する、
請求項１に記載の睡眠判定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被験者の睡眠段階を判定する睡眠判定装置に関するものである。
【背景技術】
【０００２】
従来、被験者の睡眠段階を判定するものとして、特開２００５−１５２３１０号公報に記載されるように、エアマットに設けた通気部を流れる流体の流速、流量センサで測定される生体の心拍変動情報及び体動情報に基づいて睡眠の段階を推定するものが知られている。この推定装置は、エアマットで横たわる被験者を無拘束状態で睡眠段階を推定しようとするものである。
【特許文献１】特開２００５−１５２３１０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、このような睡眠段階の推定装置にあっては、睡眠段階の推定にあたりエアマットなどの特別な器具が必要となり、体動情報の取得も必要となる。エアマットなどの特別な装置を用いて睡眠段階を判定できない場合や体動情報が得られない場合など、心拍情報だけで被験者の睡眠段階をできるだけ正確に判定できる装置の開発が望まれている。
【０００４】
そこで本発明は、このような問題点を解決するためになされたものであって、被験者の心拍情報に基づいて精度よく睡眠段階を判定できる睡眠判定装置を提供することを目的とする。
【課題を解決するための手段】
【０００５】
すなわち、本発明に係る睡眠判定装置は、被験者の心拍を取得する心拍取得手段と、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期、その心拍周期の前後の差分値、又はその差分値の前後の差分値について、少なくとも一つを演算する差分演算手段と、前記差分演算手段により演算された心拍周期及び差分値について、連続して演算された二つの値を座標値とする座標点を繰り返し演算してリターンマップを生成するリターンマップ生成手段と、前記リターンマップ生成手段により生成されたリターンマップの特徴量を演算する特徴量演算手段と、前記特徴量演算手段により演算された特徴量に基づいて前記被験者の睡眠段階を判定する睡眠判定手段とを備えて構成されている。
【０００６】
この発明によれば、被験者の心拍周期及びその差分値を座標値とする座標点を心拍ごとに繰り返し演算してリターンマップを生成し、そのリターンマップから生成される特徴量に基づいて被験者の睡眠段階を判定することにより、睡眠段階に応じて変化する心拍周期又はその差分値の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。
【０００７】
また本発明に係る睡眠判定装置において、前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、前記特徴量演算手段は、前記リターンマップ生成手段により演算される座標点の変動量を前記特徴量として演算し、前記睡眠判定手段は、前記特徴量演算手段により演算された前記座標点の変動量のバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定することが好ましい。
【０００８】
この発明によれば、リターンマップの座標点の変動量を演算し、その座標点の変動量のバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期のリターンマップの座標点の変動量を用いて、被験者の睡眠段階を精度よく判定することができる。
【０００９】
また本発明に係る睡眠判定装置において、前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、前記特徴量演算手段は、前記リターンマップ生成手段により演算される座標点の変動量を前記特徴量として演算し、前記睡眠判定手段は、前記特徴量演算手段により演算された前記座標点の変動量のヒストグラムにおける分散が小さいほど睡眠状態が深いとして睡眠段階を判定することが好ましい。
【００１０】
この発明によれば、リターンマップの座標点の変動量を演算し、その座標点の変動量のヒストグラムにおける歪度が小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期のリターンマップの座標点の変動量のヒストグラムにおける歪度を用いて、被験者の睡眠段階を精度よく判定することができる。
【００１１】
また本発明に係る睡眠判定装置において、前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、前記特徴量演算手段は、前記リターンマップ生成手段により演算される複数の座標点における集積度を前記特徴量として演算し、前記睡眠判定手段は、前記特徴量演算手段により演算された前記座標点の集積度が高いほど睡眠状態が深いとして睡眠段階を判定することが好ましい。
【００１２】
この発明によれば、リターンマップの複数の座標点の集積度を演算し、その座標点の集積度が高いほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する座標点の集積度を用いて被験者の睡眠段階を精度よく判定することができる。
【００１３】
また本発明に係る睡眠判定装置において、前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期を演算し、前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、前記特徴量演算手段は、前記リターンマップ生成手段により連続して演算される三つの座標点により形成される三角形の重心位置の変動量を前記特徴量として繰り返し演算し、前記睡眠判定手段は、前記特徴量演算手段により演算された前記三角形の重心位置の変動量が小さいほど睡眠状態が深いとして睡眠段階を判定することが好ましい。
【００１４】
この発明によれば、リターンマップ生成手段により連続して演算される三つの座標点によって形成される三角形の重心位置の変動量を演算し、その演算された三角形の重心位置の変動量が小さいほど睡眠状態が深いとして睡眠段階を判定する。これにより、睡眠段階によって変化するリターンマップ上の多角形の重心位置の変動量を用いて、被験者の睡眠段階を精度よく判定することができる。
【００１５】
また本発明に係る睡眠判定装置において、前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期の前後の差分値を演算し、前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの差分値を座標値とする座標点を繰り返し演算し、前記特徴量演算手段は、前記リターンマップ生成手段により連続して演算される三つの座標点により形成される三角形の頂点の変動量を前記特徴量として演算し、前記睡眠判定手段は、前記特徴量演算手段により演算された前記変動量が小さいほど睡眠状態が深いとして睡眠段階を判定することが好ましい。
【００１６】
この発明によれば、心拍に基づいて心拍ごとの心拍周期の前後の差分値を演算し、連続して演算された二つの差分値を座標値とする座標点を繰り返し演算してリターンマップを生成し、連続して演算される三つの座標点により形成される三角形の頂点の変動量を特徴量として演算し、演算された変動量が小さいほど睡眠状態が深いとして睡眠段階を判定する。これにより、睡眠段階に応じて変化するリターンマップに形成される三角形の頂点の変動量を用いて、被験者の睡眠段階を精度よく判定することができる。
【００１７】
また本発明に係る睡眠判定装置において、前記差分演算手段は、前記心拍取得手段により取得された心拍に基づいて心拍ごとの心拍周期の前後の差分値を演算し、前記リターンマップ生成手段は、前記差分演算手段により連続して演算された二つの差分値を座標値とする座標点を繰り返し演算し、前記特徴量演算手段は、前記リターンマップ生成手段により連続して演算される三つの座標点により形成される三角形の辺の傾きを前記特徴量として演算し、前記睡眠判定手段は、前記特徴量演算手段により演算された前記三角形の辺の傾きのバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定することが好ましい。
【００１８】
この発明によれば、心拍に基づいて心拍ごとの心拍周期の前後の差分値を演算し、連続して演算された二つの差分値を座標値とする座標点を繰り返し演算してリターンマップを生成し、連続して演算される三つの座標点により形成される三角形の辺の傾きを特徴量として演算し、演算された三角形の辺の傾きのバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定する。これにより、睡眠段階に応じて変化するリターンマップに形成される三角形の辺の傾きのバラツキ度合いを用いて、被験者の睡眠段階を精度よく判定することができる。
【発明の効果】
【００１９】
本発明によれば、被験者の心拍情報に基づいて精度よく被験者の睡眠段階を判定することができる。
【発明を実施するための最良の形態】
【００２０】
以下、添付図面を参照して本発明の実施の形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。
（第一実施形態）
【００２１】
図１は本発明の第一実施形態に係る睡眠判定装置の概略構成図である。
【００２２】
図１に示すように、本実施形態に係る睡眠判定装置１は、被験者の睡眠状態を判定する装置であって、例えば車両に設置され、駐車場等で安静中の被験者である運転者の睡眠段階、すなわち睡眠の深さレベル（睡眠レベル）の判定に適用されるものである。
【００２３】
この睡眠判定装置１は、例えば心拍検出部２、ＥＣＵ（ElectronicControl Unit）３及び出力部４を備えて構成される。心拍検出部２は、被験者の心拍を取得する心拍取得手段として機能するものであって、例えば心電図検査器に用いられる電極部のように心拍に対応した微弱な電流を検知するものが用いられる。なお、この心拍検出部２としては、被験者の心電図波形などの心拍状態を検出できるものであればいずれのものを用いてもよい。また、車両の運転者の心拍状態を検出する場合、例えば心拍検出部２が運転シートに設置されることが好ましい。なお、ここでいう心拍とは、心拍に関連する脈拍を含むものである。例えば、心拍検出部２としては、被験者の脈拍を心拍として検出するものであってもよい。
【００２４】
ＥＣＵ３は、睡眠判定装置１全体の制御を行うものであって、例えばＣＰＵ、ＲＯＭ、ＲＡＭを含むコンピュータを主体として構成されている。ＥＣＵ３は、差分演算部３１、リターンマップ生成部３２、特徴量演算部３３、睡眠判定部３４、データベース３５を備えている。
【００２５】
差分演算部３１は、心拍検出部２から出力される心拍信号に基づいて心拍ごとの心拍周期を演算する差分演算手段と機能するものである。
【００２６】
例えば、図２に示すように、心拍検出部２から所定の周期でＲ波を生ずる心拍信号が入力される。この心拍信号を信号処理することによってＲ波が抽出され、Ｒ波とＲ波の時間Ｘ（ＲＲＩ）が一つの心拍周期として演算され、連続して心拍周期Ｘ０、Ｘ１、…、Ｘｎが演算される。
【００２７】
図１において、リターンマップ生成部３２は、差分演算部３１により演算された心拍周期について、連続して演算された二つの値を座標値とする座標点を繰り返し演算してリターンマップを生成するものであって、リターンマップ生成手段として機能するものである。ここでリターンマップとは、二つの値の座標点を打点した地図のことである。
【００２８】
図３は、連続して得られた心拍周期Ｘｋ−１、Ｘｋを座標値（パラメータ）とする座標点からなる心拍周期（０次差分値）のリターンマップの一例を示したものである。例えば、Ｘｋ−１を横軸の座標値とし、Ｘｋを縦軸の座標値とする座標点ｐをプロットして心拍周期のリターンマップが生成される。このとき、予め設定された所定時間（例えば、計測時間３０秒）における心拍周期についてリターンマップを生成すればよい。
【００２９】
図３において、座標点ｐ、ｐの間の線分Ｌは、座標点ｐの推移を示すものである。この線分Ｌの長さを演算することにより、リターンマップにおける座標点の変動量を取得することができる。
【００３０】
なお、このリターンマップ生成部３２において、連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算するものであれば、その座標点をリターンマップとしてプロットしない場合もある。例えば、連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算した座標点のデータが実質的にリターンマップとして機能する場合、複数の座標軸を有するリターンマップの形で生成しなくてもよい。
【００３１】
図１において、特徴量演算部３３は、リターンマップ生成部３２により生成されたリターンマップの特徴量を演算するものであって、特徴量演算手段と機能するものである。リターンマップの特徴量として、例えばリターンマップを構成する座標点の変動量（ステップサイズ）が演算される。具体的には、図３において、連続して演算される座標点ｐ、ｐの間の距離、すなわち線分Ｌの長さが座標点の変動量として演算される。
【００３２】
図１において、睡眠判定部３４は、リターンマップの特徴量に基づいて被験者の睡眠段階を判定するものであって、睡眠判定手段として機能するものである。例えば、睡眠判定部３４は、特徴量演算部３３により演算された座標点の変動量のバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定する。
【００３３】
具体的には、現在から過去における所定の時間（例えば３０秒）以内の座標点の変動量又は現在から過去における所定のデータ数の（例えば３０個）の座標点の変動量についてバラツキを演算し、その変動量のバラツキに応じて睡眠段階Ｓ０〜Ｓ４を判定する。この睡眠段階の判定は、新たな座標点のデータが演算されるごとに行うことが好ましい。睡眠段階Ｓ０は覚醒状態の段階であり、Ｓ１、Ｓ２、Ｓ３、Ｓ４の順に睡眠が深い状態になっていく。
【００３４】
座標点の変動量のバラツキは、例えば座標点の変動量のデータについて、偏差や分散を算出して演算すればよい。また、睡眠段階Ｓ０〜Ｓ４の判定は、座標点の変動量のバラツキの閾値を睡眠段階ごとに予め設定しておき、その閾値との大小関係に応じて行えばよい。
【００３５】
図１において、データベース３５は、睡眠判定に用いられるデータを記録するものであり、例えばリターンマップの座標点データ、特徴量データ、睡眠判定に用いられる閾値データを記録している。出力部４は、ＥＣＵ３により判定された被験者の睡眠段階を出力するものであり、例えば表示モニタなどが用いられる。
【００３６】
次に本実施形態に係る睡眠判定装置の動作について説明する。
【００３７】
図４は、本実施形態に係る睡眠判定装置の動作を示すフローチャートである。この図４の制御処理は、例えばＥＣＵ３によって所定の周期を繰り返し実行される。
【００３８】
まず、図１において、心拍検出部２により被験者の心拍信号が検出される。そして、この心拍信号に基づいて図４の睡眠判定処理により被験者の睡眠段階の判定が行われる。
【００３９】
図４のステップＳ１０（以下、単に「Ｓ１０」という。ステップＳ１２以下についても同じ。）に示すように、心拍信号の読み込み処理が行われる。心拍信号の読み込み処理は、心拍検出部２により出力される被験者の心拍信号を読み込む処理である。心拍信号は、図２に示すように心拍に応じて繰り返しＲ波を生じる波形で入力される。
【００４０】
そして、図４のＳ１２に移行し、差分演算処理が行われる。差分演算処理は、心拍信号に基づいて心拍ごとの心拍周期を演算する処理である。例えば、心拍周期は、心拍信号からＲ波を検出し、繰り返し検出されるＲ波とＲ波の間の時間として算出される。
【００４１】
そして、Ｓ１４に移行し、リターンマップ生成処理が行われる。リターンマップ生成処理は、Ｓ１２の差分演算処理により連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算してリターンマップを生成する処理である。
【００４２】
例えば、図３に示すように、連続して得られた心拍周期Ｘｋ−１、Ｘｋを座標値（パラメータ）とする座標点ｐが繰り返し演算される。すなわち、Ｘｋ−１を横軸の座標値とし、Ｘｋを縦軸の座標値とする座標点ｐ（Ｘｋ−１、Ｘｋ）をプロットして心拍周期のリターンマップが生成される。このとき、予め設定された所定時間（例えば、計測時間３０秒）以内における心拍周期についてリターンマップが生成される。
【００４３】
なお、このリターンマップ生成処理において、連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算するものであれば、その座標点をリターンマップとしてプロットしない場合もある。例えば、連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算した座標点のデータが実質的にリターンマップとして機能する場合、複数の座標軸を有するリターンマップの形で生成しなくてもよい。
【００４４】
そして、図４のＳ１６に移行し、特徴量演算処理が行われる。特徴量演算処理は、Ｓ１４のリターンマップ生成処理により生成されたリターンマップの特徴量を演算する処理である。リターンマップの特徴量として、例えばリターンマップを構成する座標点の変動量（ステップサイズ）が演算される。具体的には、図３において、連続して演算される座標点ｐ、ｐの間の距離、すなわち線分Ｌの長さが座標点の変動量として演算される。演算された座標点の変動量は、データベース３５に記録される。座標点の変動量は、心拍信号の入力されるごとに演算され、随時データベースに記録される。
【００４５】
そして、図４のＳ１８に移行し、睡眠段階判定処理が行われる。睡眠段階判断処理は、Ｓ１６にて演算されたリターンマップの特徴量に基づいて被験者の睡眠段階を判定する処理である。例えば、Ｓ１６にて演算された座標点の変動量のバラツキが小さいほど睡眠状態が深いとして睡眠段階が判定される。
【００４６】
具体的に説明すると、今回演算された座標点の変動量と過去に演算された所定回数の座標点の変動量（例えば過去３０回分の変動量）について、座標点の変動量のバラツキを示す値として標準偏差又は分散値が演算され、その標準偏差又は分散値に応じて睡眠段階が判定される。
【００４７】
睡眠段階の判定しきい値として、睡眠段階Ｓ０と睡眠段階Ｓ１を判別するためのしきい値、睡眠段階Ｓ１と睡眠段階Ｓ２を判別するためのしきい値、睡眠段階Ｓ２と睡眠段階Ｓ３を判別するためのしきい値、睡眠段階Ｓ３と睡眠段階Ｓ４を判別するためのしきい値がそれぞれ設定される。そして、座標点の変動量のバラツキを示す値とこれらの判別のしきい値とが比較され、被験者の睡眠段階がＳ０〜Ｓ４のいずれであるかが判定される。判定された睡眠段階の結果は、出力部４を通じて出力される。そして、Ｓ１８の睡眠段階判定処理を終えたら、一連の制御処理を終了する。
【００４８】
以上のように、本実施形態に係る睡眠判定装置によれば、被験者の心拍周期を座標値とする座標点を心拍ごとに繰り返し演算してリターンマップを生成し、そのリターンマップの特徴量に基づいて被験者の睡眠段階を判定することにより、睡眠段階に応じて変化する心拍周期の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。
【００４９】
また、本実施形態に係る睡眠判定装置においては、リターンマップの座標点の変動量を演算し、その座標点の変動量のバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期のリターンマップの座標点の変動量を用いて、被験者の睡眠段階を精度よく判定することができる。
【００５０】
図５〜９は、被験者の心拍周期のリターンマップ（０次リターンマップ：心拍周期の０次差分のリターンマップ）を睡眠段階Ｓ０〜Ｓ４ごとに示したものである。図５は覚醒状態の睡眠段階Ｓ０のリターンマップ、図６は睡眠段階Ｓ１のリターンマップ、図７は睡眠段階Ｓ２のリターンマップ、図８は睡眠段階Ｓ３のリターンマップ、図９は睡眠段階Ｓ４のリターンマップである。この図５〜９において、横軸のＸｋ−１は心拍周期の前回値であり、縦軸のＸｋは心拍周期の今回値である。また、図５〜９において、座標点のみならず、座標点から次回値の座標点への線分も図示している。
【００５１】
図５〜９に示されるように、睡眠段階がＳ０からＳ４へ深くなるに連れて、座標点から次回値の座標点への線分の長さ、すなわち座標点の変動量が一定の値に近づく傾向がある。つまり、睡眠状態がＳ０からＳ４へ深くなるほど、座標点の変動量のバラツキが小さくなることが分かる。
【００５２】
このため、被験者の心拍周期についてリターンマップを生成し座標点の変動量のバラツキ度合いを演算することにより、その座標点の変動量のバラツキ度合いに応じて睡眠段階を精度よく判定することが可能となる。
【００５３】
また、本実施形態に係る睡眠判定装置においては、被験者の心拍状態のみに基づいて判定するため、心電図や呼吸のデータがなくても被験者の睡眠段階を判定が可能である。
【００５４】
さらに、心拍信号を受けるごとに心拍周期を更新してリアルタイムで被験者の睡眠段階を判定することができる。
（第二実施形態）
【００５５】
次に本発明の第二実施形態に係る睡眠判定装置について説明する。
【００５６】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置とほぼ同様な構成を備え、図４に示す第一実施形態に係る睡眠判定装置とほぼ同様な睡眠判定処理を実行するものであって、座標点の変動量のヒストグラムにおける歪度に基づいて睡眠段階の判定を行う点で、座標点の変動量のバラツキに基づいて睡眠段階の判定を行う第一実施形態に係る睡眠判定装置と異なっている。
【００５７】
すなわち、本実施形態に係る睡眠判定装置は、被験者の心拍に基づいて心拍ごとの心拍周期を演算し、連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、座標点の変動量（ステップサイズ）を特徴量として演算し、座標点の変動量のヒストグラムにおける歪度が小さいほど睡眠状態が深いとして睡眠段階を判定するものである。
【００５８】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置と同様に、心拍検出部２、ＥＣＵ３及び出力部４を備えて構成される。心拍検出部２及び出力部４は、第一実施形態に係る睡眠判定装置と同様なものが用いられる。本実施形態に係る睡眠判定装置のＥＣＵ３は、特徴量演算部３３と睡眠判定部３４の機能が第一実施形態に係る睡眠判定装置のものと異なっている。
【００５９】
特徴量演算部３３は、リターンマップ生成部３２により生成されたリターンマップの特徴量を演算するものであって、特徴量演算手段と機能するものである。この特徴量演算部３３は、リターンマップの特徴量として、座標点の変動量の分布特性を演算する。
【００６０】
例えば、リターンマップを構成する座標点の変動量（ステップサイズ）が演算され、現在から過去における所定の時間（例えば３０秒）以内に演算される座標点の変動量についてヒストグラムが生成され、そのヒストグラムにおける分布特性をとらえた統計量が演算される。
【００６１】
具体的に示すと、図１０に示すように、所定時間に演算される座標点の変動量ごとの頻度を示すヒストグラムＨが生成され、このヒストグラムＨの統計量が演算される。この統計量は、ヒストグラムの非対称性の度合いを示すものであり、左右対称の場合にはゼロとなり、左右対称性が低いほど大きな値となる。具体的な統計量の値は、例えば統計学の公知の手法により算出すればよい。
【００６２】
睡眠判定部３４は、リターンマップの特徴量に基づいて被験者の睡眠段階を判定するものであって、リターンマップの座標点の変動量におけるヒストグラム統計量に基づいて睡眠段階を判定する。
【００６３】
例えば、睡眠判定部３４は、特徴量演算部３３により演算された座標点の変動量におけるヒストグラムの統計量が小さいほど睡眠状態が深いとして睡眠段階を判定する。
【００６４】
具体的には、睡眠段階Ｓ０、Ｓ１を判別する統計量しきい値、睡眠段階Ｓ１、Ｓ２を判別する統計量しきい値、睡眠段階Ｓ２、Ｓ３を判別する統計量しきい値、睡眠段階Ｓ３、Ｓ４を判別する統計量しきい値を予め設定しておき、特徴量演算部３３により演算された統計量の値を統計量しきい値を比較することにより、被験者の睡眠段階が判定される。
【００６５】
次に、本実施形態に係る睡眠判定装置の動作を説明する。
【００６６】
まず、図１において、心拍検出部２により被験者の心拍信号が検出される。そして、この心拍信号に基づいて図４の睡眠判定処理により被験者の睡眠段階の判定が行われる。
【００６７】
そして、図４に示すように、心拍信号の読み込み処理、差分演算処理、リターンマップ生成処理が行われる（Ｓ１０、Ｓ１２、Ｓ１４）。これらの処理は、第一実施形態に説明した処理と同様に行われる。
【００６８】
そして、Ｓ１６に移行し、特徴量演算処理が行われる。特徴量演算処理は、Ｓ１４のリターンマップ生成処理により生成されたリターンマップの特徴量を演算する処理である。リターンマップの特徴量として、座標点の変動量のヒストグラムにおける統計量が演算される。例えば、リターンマップを構成する座標点の変動量（ステップサイズ）が演算され、今回演算された座標点の変動量と過去に演算された所定回数の座標点の変動量（例えば過去３０回分の変動量）について、図１０に示すようなヒストグラムが生成され、そのヒストグラムの統計量が演算される。
【００６９】
そして、図４のＳ１８に移行し、睡眠段階判定処理が行われる。睡眠段階判断処理は、Ｓ１６にて演算されたリターンマップの特徴量に基づいて被験者の睡眠段階を判定する処理である。例えば、Ｓ１６にて演算された座標点の変動量のヒストグラムの統計量が小さいほど睡眠状態が深いとして睡眠段階が判定される。
【００７０】
具体的に説明すると、演算されたヒストグラムの統計量の値と予め設定された統計量しきい値が比較される。そして、演算されたヒストグラムの統計量の値が睡眠段階Ｓ０、Ｓ１を判別する統計量しきい値より大きい場合は、被験者の睡眠段階がＳ０であると判定される。
【００７１】
演算されたヒストグラムの統計量の値が睡眠段階Ｓ０、Ｓ１を判別する統計量しきい値以下であり睡眠段階Ｓ１、Ｓ２を判別する統計量しきい値より大きい場合は、被験者の睡眠段階がＳ１であると判定される。演算されたヒストグラムの統計量の値が睡眠段階Ｓ１、Ｓ２を判別する統計量しきい値以下であり睡眠段階Ｓ２、Ｓ３を判別する統計量しきい値より大きい場合は、被験者の睡眠段階がＳ２であると判定される。
【００７２】
演算されたヒストグラムの統計量の値が睡眠段階Ｓ２、Ｓ３を判別する統計量しきい値以下であり睡眠段階Ｓ３、Ｓ４を判別する統計量しきい値より大きい場合は、被験者の睡眠段階がＳ３であると判定される。演算されたヒストグラムの統計量の値が睡眠段階Ｓ３、Ｓ４を判別する統計量しきい値以下である場合は、被験者の睡眠段階がＳ４であると判定される。そして、Ｓ１８の睡眠段階判定処理を終えたら、一連の制御処理を終了する。
【００７３】
以上のように、本実施形態に係る睡眠判定装置によれば、第一実施形態に係る睡眠判定装置と同様に、被験者の心拍周期を座標値とする座標点を心拍ごとに繰り返し演算してリターンマップを生成し、そのリターンマップの特徴量に基づいて被験者の睡眠段階を判定することにより、睡眠段階に応じて変化する心拍周期の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。
【００７４】
また、本実施形態に係る睡眠判定装置においては、リターンマップの座標点の変動量の分布状態（ヒストグラム）の統計量を演算し、その統計量が小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期のリターンマップの座標点の変動状態を用いて、被験者の睡眠段階を精度よく判定することができる。
【００７５】
図１１は、睡眠段階ごとのリターンマップの座標点の変動量のヒストグラムを示したものである。図１１（ａ）は睡眠段階Ｓ０におけるリターンマップの座標点の変動量のヒストグラム、図１１（ｂ）は睡眠段階Ｓ１におけるリターンマップの座標点の変動量のヒストグラム、図１１（ｃ）は睡眠段階Ｓ２におけるリターンマップの座標点の変動量のヒストグラム、図１１（ｄ）は睡眠段階Ｓ３におけるリターンマップの座標点の変動量のヒストグラム、図１１（ｅ）は睡眠段階Ｓ４におけるリターンマップの座標点の変動量のヒストグラムである。
【００７６】
この図１１（ａ）〜（ｅ）のヒストグラムを見ると、睡眠段階Ｓ０からＳ４へ行くほどヒストグラムの左右対称性が高くなることが分かる。このため、ヒストグラムの非対称性を示す指標となるヒストグラムの統計量に基づいて睡眠段階の判定を行うことにより、精度の高い睡眠段階の判定が行える。
【００７７】
なお、本実施形態では、リターンマップの座標点の変動量における分布状態の統計量に基づいて被験者の睡眠段階を判定する場合について説明したが、リターンマップの座標点の変動量における分布状態の尖度に基づいて被験者の睡眠段階を判定してもよい。
【００７８】
図１１を見ると、睡眠段階Ｓ０からＳ４へ行くほど変動量の分布状態（ヒストグラム）の尖度が高くなることが分かる。このため、ヒストグラムの尖度が高いほど被験者の睡眠状態が深いとして睡眠段階の判定を行うことができる。
（第三実施形態）
【００７９】
次に本発明の第三実施形態に係る睡眠判定装置について説明する。
【００８０】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置とほぼ同様な構成を備え、図４に示す第一実施形態に係る睡眠判定装置とほぼ同様な睡眠判定処理を実行するものであって、リターンマップの座標点の集積度に基づいて睡眠段階の判定を行う点で、座標点の変動量のバラツキに基づいて睡眠段階の判定を行う第一実施形態に係る睡眠判定装置と異なっている。
【００８１】
すなわち、本実施形態に係る睡眠判定装置は、被験者の心拍に基づいて心拍ごとの心拍周期を演算し、連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、座標点の集積度を特徴量として演算し、その座標点の集積度が高いほど睡眠状態が深いとして睡眠段階を判定するものである。
【００８２】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置と同様に、心拍検出部２、ＥＣＵ３及び出力部４を備えて構成される。心拍検出部２及び出力部４は、第一実施形態に係る睡眠判定装置と同様なものが用いられる。本実施形態に係る睡眠判定装置のＥＣＵ３は、特徴量演算部３３と睡眠判定部３４の機能が第一実施形態に係る睡眠判定装置のものと異なっている。
【００８３】
特徴量演算部３３は、リターンマップ生成部３２により生成されたリターンマップの特徴量を演算するものであって、特徴量演算手段と機能するものである。この特徴量演算部３３は、リターンマップの特徴量として、座標点の集積度を演算する。
【００８４】
例えば、図１２に示すように、リターンマップを構成する複数の座標点ｐが演算され、それらの座標点ｐについて重心位置Ｇが演算され、その重心位置Ｇを中心とする所定の半径の円Ｃの内部に含まれる座標点ｐの数が座標点の集積度として演算される。円Ｃの半径は、睡眠段階を判別しやすい大きさに設定され、実際の座標点データの分布状態に応じて設定することが好ましい。
【００８５】
なお、座標点の集積度の度合いを算出することができれば、座標点の集積度の算出の仕方は上述したものに限られるものではない。例えば、所定の数の座標点を含む円の半径に応じて集積度を算出するなどしてもよい。
【００８６】
睡眠判定部３４は、リターンマップの特徴量に基づいて被験者の睡眠段階を判定するものであって、リターンマップの座標点の集積度が高いほど睡眠状態が深いとして被験者の睡眠段階を判定する。
【００８７】
例えば、睡眠段階Ｓ０、Ｓ１を判別する集積しきい値、睡眠段階Ｓ１、Ｓ２を判別する集積しきい値、睡眠段階Ｓ２、Ｓ３を判別する集積しきい値、睡眠段階Ｓ３、Ｓ４を判別する集積しきい値が予め設定され、特徴量演算部３３により演算された集積度の値と集積しきい値を比較することにより、被験者の睡眠段階が判定される。
【００８８】
次に、本実施形態に係る睡眠判定装置の動作を説明する。
【００８９】
まず、図１において、心拍検出部２により被験者の心拍信号が検出される。そして、この心拍信号に基づいて図４の睡眠判定処理により被験者の睡眠段階の判定が行われる。
【００９０】
そして、図４に示すように、心拍信号の読み込み処理、差分演算処理、リターンマップ生成処理が行われる（Ｓ１０、Ｓ１２、Ｓ１４）。これらの処理は、第一実施形態に説明した処理と同様に行われる。
【００９１】
そして、Ｓ１６に移行し、特徴量演算処理が行われる。特徴量演算処理は、Ｓ１４のリターンマップ生成処理により生成されたリターンマップの特徴量を演算する処理である。リターンマップの特徴量として、座標点の集積度が演算される。
【００９２】
例えば、図１２に示すように、リターンマップを構成する複数の座標点ｐが演算され、それらの座標点ｐについて重心位置Ｇが演算され、その重心位置Ｇを中心とする所定の半径の円Ｃの内部に含まれる座標点ｐの数が座標点の集積度として演算される。このとき、座標点ｐとしては、例えば所定時間内に演算されたもの又は所定回数演算されたものが用いられる。
【００９３】
なお、座標点の集積度の度合いを算出することができれば、座標点の集積度の算出の仕方は上述したものに限られるものではない。例えば、所定の数の座標点を含む円の半径に応じて集積度を算出するなどしてもよい。
【００９４】
そして、図４のＳ１８に移行し、睡眠段階判定処理が行われる。睡眠段階判断処理は、Ｓ１６にて演算されたリターンマップの特徴量に基づいて被験者の睡眠段階を判定する処理である。例えば、Ｓ１６にて演算された座標点の集積度が高いほど睡眠状態が深いとして睡眠段階が判定される。
【００９５】
具体的には、図１２に示すように、重心位置Ｇを中心とする所定の半径の円Ｃの内部に含まれる座標点ｐの数と集積しきい値を比較し、被験者の睡眠段階が判定される。そして、Ｓ１８の睡眠段階判定処理を終えたら、一連の制御処理を終了する。
【００９６】
以上のように、本実施形態に係る睡眠判定装置によれば、第一実施形態に係る睡眠判定装置と同様に、被験者の心拍周期を座標値とする座標点を心拍ごとに繰り返し演算してリターンマップを生成し、そのリターンマップの特徴量に基づいて被験者の睡眠段階を判定することにより、睡眠段階に応じて変化する心拍周期の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。
【００９７】
また、本実施形態に係る睡眠判定装置においては、リターンマップの座標点の集積度を演算し、その集積度が小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期のリターンマップの座標点の変動状態を用いて、被験者の睡眠段階を精度よく判定することができる。
【００９８】
すなわち、図５〜９のリターンマップを見ると、睡眠状態がＳ０からＳ４へ深くなるに連れてリターンマップの座標点の集積度が徐々に高くなることが分かる。このため、リターンマップの座標点の集積度に基づいて睡眠段階を判定することにより、被験者の睡眠段階を精度よく判定することができる。
（第四実施形態）
【００９９】
次に本発明の第四実施形態に係る睡眠判定装置について説明する。
【０１００】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置とほぼ同様な構成を備え、図４に示す第一実施形態に係る睡眠判定装置とほぼ同様な睡眠判定処理を実行するものであって、リターンマップにおいて連続する三つの座標点により形成される三角形の重心位置の変動量に基づいて睡眠段階の判定を行う点で、座標点の変動量のバラツキに基づいて睡眠段階の判定を行う第一実施形態に係る睡眠判定装置と異なっている。
【０１０１】
すなわち、本実施形態に係る睡眠判定装置は、被験者の心拍に基づいて心拍ごとの心拍周期を演算し、連続して演算された二つの心拍周期の値を座標値とする座標点を繰り返し演算し、連続する三つの座標点により形成される三角形の重心位置の変動量を特徴量として演算し、その重心位置の変動量が小さいほど睡眠状態が深いとして睡眠段階を判定するものである。
【０１０２】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置と同様に、心拍検出部２、ＥＣＵ３及び出力部４を備えて構成される。心拍検出部２及び出力部４は、第一実施形態に係る睡眠判定装置と同様なものが用いられる。本実施形態に係る睡眠判定装置のＥＣＵ３は、特徴量演算部３３と睡眠判定部３４の機能が第一実施形態に係る睡眠判定装置のものと異なっている。
【０１０３】
特徴量演算部３３は、リターンマップ生成部３２により生成されたリターンマップの特徴量を演算するものであって、特徴量演算手段と機能するものである。この特徴量演算部３３は、リターンマップの特徴量として、リターンマップにおいて連続する三つの座標点により形成される三角形の重心位置の変動量を演算する。
【０１０４】
例えば、図１３に示すように、リターンマップにおいて連続して演算される三つの座標点ｐ１〜ｐ３の座標値に基づいてその三つの座標点ｐ１〜ｐ３により形成される三角形の重心座標が演算される。そして、次の座標点ｐ４が演算されると、三つの座標点ｐ２〜ｐ４の座標値に基づいて三つの座標点ｐ２〜ｐ４により形成される三角形の重心座標が演算される。そして、三角形の重心位置が演算されるごとに前回演算された重心位置からの変動量が特徴量として演算される。
【０１０５】
睡眠判定部３４は、リターンマップの特徴量に基づいて被験者の睡眠段階を判定するものであって、リターンマップにおける三つの座標点により形成される三角形の重心位置の変動量が小さいほど睡眠状態が深いとして睡眠段階を判定するものである。
【０１０６】
次に、本実施形態に係る睡眠判定装置の動作を説明する。
【０１０７】
まず、図１において、心拍検出部２により被験者の心拍信号が検出される。そして、この心拍信号に基づいて図４の睡眠判定処理により被験者の睡眠段階の判定が行われる。
【０１０８】
そして、図４に示すように、心拍信号の読み込み処理、差分演算処理、リターンマップ生成処理が行われる（Ｓ１０、Ｓ１２、Ｓ１４）。これらの処理は、第一実施形態に説明した処理と同様に行われる。
【０１０９】
そして、Ｓ１６に移行し、特徴量演算処理が行われる。特徴量演算処理は、Ｓ１４のリターンマップ生成処理により生成されたリターンマップの特徴量を演算する処理である。リターンマップの特徴量として、リターンマップにおいて連続する三つの座標点により形成される三角形の重心位置の変動量が演算される。
【０１１０】
例えば、図１３に示すように、リターンマップにおいて連続して演算される三つの座標点ｐ１〜ｐ３の座標値に基づいてその三つの座標点ｐ１〜ｐ３により形成される三角形の重心座標が演算される。そして、次の座標点ｐ４が演算されると、三つの座標点ｐ２〜ｐ４の座標値に基づいて三つの座標点ｐ２〜ｐ４により形成される三角形の重心座標が演算される。このように、三角形の重心位置が演算されるごとに前回演算された重心位置からの変動量が特徴量として演算される。
【０１１１】
そして、図４のＳ１８に移行し、睡眠段階判定処理が行われる。睡眠段階判断処理は、Ｓ１６にて演算されたリターンマップの特徴量に基づいて被験者の睡眠段階を判定する処理である。例えば、Ｓ１６にて演算された三角形の重心位置の変動量が小さいほど睡眠状態が深いとして睡眠段階が判定される。
【０１１２】
例えば、所定時間内の三角形の重心位置の変動量又は所定回数の三角形の重心位置の変動量をヒストグラムとし、三角形の重心位置の変動量の分布状態に基づいてその変動量が小さいほど睡眠状態が深いと判定する。また、所定時間内の三角形の重心位置の変動量又は所定回数の三角形の重心位置の変動量について平均値を算出し、その平均値を予め設定されるしきい値と比較して睡眠段階を判定してもよい。そして、Ｓ１８の睡眠段階判定処理を終えたら、一連の制御処理を終了する。
【０１１３】
以上のように、本実施形態に係る睡眠判定装置によれば、第一実施形態に係る睡眠判定装置と同様に、被験者の心拍周期を座標値とする座標点を心拍ごとに繰り返し演算してリターンマップを生成し、そのリターンマップの特徴量に基づいて被験者の睡眠段階を判定することにより、睡眠段階に応じて変化する心拍周期の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。
【０１１４】
また、本実施形態に係る睡眠判定装置においては、リターンマップにおいて連続する三つの座標点により形成される三角形の重心位置の変動量を演算し、その変動量が小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期のリターンマップの座標点の変動状態を用いて、被験者の睡眠段階を精度よく判定することができる。
【０１１５】
すなわち、図５〜９のリターンマップを見ると、睡眠状態がＳ０からＳ４へ深くなるに連れてリターンマップの座標点により形成される三角形の重心位置の変動量が徐々に小さくなることが分かる。このため、リターンマップの座標点により形成される三角形の重心位置の変動量に基づいて睡眠段階を判定することにより、被験者の睡眠段階を精度よく判定することができる。
（第五実施形態）
【０１１６】
次に本発明の第五実施形態に係る睡眠判定装置について説明する。
【０１１７】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置とほぼ同様な構成を備え、図４に示す第一実施形態に係る睡眠判定装置とほぼ同様な睡眠判定処理を実行するものであって、被験者の心拍周期の差分値を演算し、連続して演算された二つの差分値を座標値とする座標点を繰り返し演算し、連続して演算される三つの座標点により形成される三角形の頂点の変動量を特徴量として演算し、演算された三角形の頂点の変動量が小さいほど睡眠状態が深いとして睡眠段階を判定するものである。
【０１１８】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置と同様に、心拍検出部２、ＥＣＵ３及び出力部４を備えて構成される。心拍検出部２及び出力部４は、第一実施形態に係る睡眠判定装置と同様なものが用いられる。本実施形態に係る睡眠判定装置のＥＣＵ３は、差分演算部３１、リターンマップ生成部３２、特徴量演算部３３、睡眠判定部３４及びデータベース３５を備えている。
【０１１９】
差分演算部３１は、心拍検出部２から出力される心拍信号に基づいて心拍ごとの心拍周期の前後の差分値を演算する差分演算手段と機能するものである。
【０１２０】
例えば、図２に示すように、心拍検出部２から所定の周期でＲ波を生ずる心拍信号が入力される。この心拍信号を信号処理することによってＲ波が抽出され、Ｒ波とＲ波の時間Ｘ（ＲＲＩ）が一つの心拍周期として演算され、連続して心拍周期Ｘ０、Ｘ１、…、Ｘｎが演算される。
【０１２１】
そして、図１４に示すように、心拍周期の前後の差分値Ｘ０（１）、Ｘ１（１）…が演算される。心拍周期の差分値Ｘｋ（１）は、次の式（１）により算出される。
【０１２２】
Ｘｋ（１）＝Ｘｋ＋１−Ｘｋ  …  （１）
【０１２３】
Ｘｋ＋１、Ｘｋは、連続して演算される心拍周期の値である。このＸｋ＋１、Ｘｋの差として、Ｘｋ（１）が演算される。
【０１２４】
図１において、リターンマップ生成部３２は、差分演算部３１により演算された心拍周期の差分値について、連続して演算された二つの値を座標値とする座標点を繰り返し演算してリターンマップを生成するものであって、リターンマップ生成手段として機能するものである。
【０１２５】
図１５は、連続して得られた心拍周期Ｘｋ−１（１）、Ｘｋ（１）を座標値（パラメータ）とする座標点からなる心拍周期の差分値（１次差分値）のリターンマップの一例を示したものである。例えば、Ｘｋ−１（１）を横軸の座標値とし、Ｘｋ（１）を縦軸の座標値とする座標点ｐをプロットして心拍周期のリターンマップが生成される。このとき、予め設定された所定時間（例えば、計測時間３０秒）における心拍周期の差分値についてリターンマップを生成すればよい。図１５において、座標点ｐ、ｐの間の線分Ｌは、座標点ｐの推移を示すものである。
【０１２６】
なお、このリターンマップ生成部３２において、連続して演算された二つの心拍周期の差分値を座標値とする座標点を繰り返し演算するものであれば、その座標点をリターンマップとしてプロットしない場合もある。例えば、連続して演算された二つの差分値を座標値とする座標点を繰り返し演算した座標点のデータが実質的にリターンマップとして機能する場合、複数の座標軸を有するリターンマップの形で生成しなくてもよい。
【０１２７】
図１において、特徴量演算部３３は、リターンマップ生成部３２により生成されたリターンマップの特徴量を演算するものであって、特徴量演算手段と機能するものである。リターンマップの特徴量として、リターンマップにおいて連続して演算される三つの座標点により形成される三角形の頂点の変動量が演算される。
【０１２８】
例えば、三つの座標点のｐｋ、ｐｋ＋１、ｐｋ＋２を頂点とする三角形において、各頂点の座標値をｐｋ（Ｘｋ、Ｘｋ＋１）、ｐｋ＋１（Ｘｋ＋１、Ｘｋ＋２）、ｐｋ＋２（Ｘｋ＋２、Ｘｋ＋３）とすると、三角形の頂点の変動量Ｚは、次の式（２）により演算される。
【数１】
2010148605.tif 000003
【０１２９】
この式（２）において、ｉは演算時の現時点を示す数値である。この式（２）は、３０個の座標値を演算した際の変動量Ｚを算出した場合の数式を示したものである。なお、三角形の頂点の変動量の演算の仕方は、上述したものに限られるものではない。
【０１３０】
睡眠判定部３４は、リターンマップの特徴量に基づいて被験者の睡眠段階を判定するものであって、睡眠判定手段として機能するものである。例えば、睡眠判定部３４は、特徴量演算部３３により演算された三角形の頂点の変動量が小さいほど睡眠状態が深いとして睡眠段階を判定する。例えば、式（２）により算出された変動量Ｚに応じて被験者の睡眠段階が判定される。このとき、睡眠段階を判定するための判定しきい値を予め設定しておくことが好ましい。
【０１３１】
また、睡眠判定の手法として、三角形の頂点の変動量のヒストグラムを生成し、そのヒストグラムの位置と形状によって睡眠段階の判定を行ってもよい。この場合、三角形の頂点の変動量の分布状態を加味して睡眠段階の判定を行えるため、より精度のよい睡眠段階の判定が行える。
【０１３２】
次に、本実施形態に係る睡眠判定装置の動作を説明する。
【０１３３】
まず、図１において、心拍検出部２により被験者の心拍信号が検出される。そして、この心拍信号に基づいて図４の睡眠判定処理により被験者の睡眠段階の判定が行われる。
【０１３４】
そして、図４に示すように、心拍信号の読み込み処理が行われる（Ｓ１０）。これらの処理は、第一実施形態に説明した処理と同様に行われる。
【０１３５】
そして、Ｓ１２に移行し、差分演算処理が行われる。差分演算処理は、心拍信号に基づいて心拍ごとの心拍周期の差分値を演算する処理である。例えば、心拍信号からＲ波を検出し、繰り返し検出されるＲ波とＲ波の間の時間として心拍周期を算出し、連続して算出される心拍周期の差として差分値が演算される。
【０１３６】
そして、Ｓ１４に移行し、リターンマップ生成処理が行われる。リターンマップ生成処理は、Ｓ１２の差分演算処理により連続して演算された二つの心拍周期の差分値を座標値とする座標点を繰り返し演算してリターンマップを生成する処理である。
【０１３７】
そして、Ｓ１６に移行し、特徴量演算処理が行われる。特徴量演算処理は、Ｓ１４のリターンマップ生成処理により生成されたリターンマップの特徴量を演算する処理である。リターンマップの特徴量として、リターンマップにおいて連続して演算される三つの座標点により形成される三角形の頂点の変動量が演算される。演算された座標点の変動量は、データベース３５に記録される。三角形の頂点の変動量は、心拍信号の入力されるごとに演算され、随時データベースに記録される。
【０１３８】
そして、Ｓ１８に移行し、睡眠段階判定処理が行われる。睡眠段階判断処理は、Ｓ１６にて演算されたリターンマップの特徴量に基づいて被験者の睡眠段階を判定する処理である。例えば、Ｓ１６にて演算された三角形の頂点の変動量が小さいほど睡眠状態が深いとして睡眠段階が判定される。そして、Ｓ１８の睡眠段階判定処理を終えたら、一連の制御処理を終了する。
【０１３９】
以上のように、本実施形態に係る睡眠判定装置によれば、被験者の心拍周期の差分値を座標値とする座標点を心拍ごとに繰り返し演算してリターンマップを生成し、そのリターンマップの特徴量に基づいて被験者の睡眠段階を判定することにより、睡眠段階に応じて変化する心拍周期の差分値の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。
【０１４０】
また、本実施形態に係る睡眠判定装置においては、リターンマップにおいて連続して演算される三つの座標点により形成される三角形の頂点の変動量を演算し、その変動量が小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期の差分値のリターンマップの座標点の変動状態を用いて、被験者の睡眠段階を精度よく判定することができる。
【０１４１】
図１６〜２０は、被験者の心拍周期の差分値のリターンマップ（１次リターンマップ：心拍周期の１次差分のリターンマップ）を睡眠段階Ｓ０〜Ｓ４ごとに示したものである。図１６は覚醒状態の睡眠段階Ｓ０のリターンマップ、図１７は睡眠段階Ｓ１のリターンマップ、図１８は睡眠段階Ｓ２のリターンマップ、図１９は睡眠段階Ｓ３のリターンマップ、図２０は睡眠段階Ｓ４のリターンマップである。この図１６〜２０において、横軸のＸｋ−１（１）は心拍周期の差分値の前回値であり、縦軸のＸｋ（１）は心拍周期の差分値の今回値である。また、図１６〜２０において、座標点のみならず、座標点から次回値の座標点への線分も図示している。
【０１４２】
図１６〜２０に示されるように、睡眠段階がＳ０からＳ４へ深くなるに連れて、連続する三つの座標点により形成される三角形の頂点が所定の箇所に集まってくる。つまり、睡眠状態がＳ０からＳ４へ深くなるほど、三角形の頂点の変動量が小さくなることが分かる。
【０１４３】
このため、被験者の心拍周期の差分値についてリターンマップを生成し座標点により形成される三角形の頂点の変動量を演算することにより、その変動量に応じて睡眠段階を精度よく判定することが可能となる。
【０１４４】
また、本実施形態に係る睡眠判定装置においては、被験者の心拍状態のみに基づいて判定するため、心電図や呼吸のデータがなくても被験者の睡眠段階を判定が可能である。
【０１４５】
さらに、心拍信号を受けるごとに心拍周期の差分値を更新してリアルタイムで被験者の睡眠段階を判定することができる。
（第六実施形態）
【０１４６】
次に本発明の第六実施形態に係る睡眠判定装置について説明する。
【０１４７】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置とほぼ同様な構成を備え、図４に示す第一実施形態に係る睡眠判定装置とほぼ同様な睡眠判定処理を実行するものであって、被験者の心拍周期の差分値を演算し、連続して演算された二つの差分値を座標値とする座標点を繰り返し演算し、連続して演算される三つの座標点により形成される三角形の辺の傾きのバラツキを特徴量として演算し、演算された三角形の辺の傾きのバラツキが小さいほど睡眠状態が深いとして睡眠段階を判定するものである。
【０１４８】
本実施形態に係る睡眠判定装置は、図１に示す第一実施形態に係る睡眠判定装置と同様に、心拍検出部２、ＥＣＵ３及び出力部４を備えて構成される。心拍検出部２及び出力部４は、第一実施形態に係る睡眠判定装置と同様なものが用いられる。本実施形態に係る睡眠判定装置のＥＣＵ３は、差分演算部３１、リターンマップ生成部３２、特徴量演算部３３、睡眠判定部３４及びデータベース３５を備えている。
【０１４９】
差分演算部３１及びリターンマップ生成部３２は、第五実施形態に係る睡眠判定装置と同様なものが用いられる。
【０１５０】
特徴量演算部３３は、リターンマップ生成部３２により生成されたリターンマップの特徴量を演算するものであって、特徴量演算手段と機能するものである。リターンマップの特徴量として、リターンマップにおいて連続して演算される三つの座標点により形成される三角形の辺の傾きのバラツキ量が演算される。
【０１５１】
例えば、リターンマップの座標点が演算されるごとに座標点から座標点への傾きを演算していき、図２１に示すように、所定時間又は所定演算回数における傾きデータのヒストグラムを生成する。このヒストグラムにおいて三峰性が表れるので、この三峰性が際立つ度合いに応じて三角形の辺の傾きのバラツキ量を演算すればよい。この三峰性が際立つ度合いが大きいほど三角形の辺の傾きのバラツキ量が小さくなることとなる。なお、三角形の辺の傾きのバラツキ量の演算の仕方は、上述したものに限られるものではない。
【０１５２】
睡眠判定部３４は、リターンマップの特徴量に基づいて被験者の睡眠段階を判定するものであって、睡眠判定手段として機能するものである。例えば、睡眠判定部３４は、特徴量演算部３３により演算された三角形の辺の傾きのバラツキ量が小さいほど睡眠状態が深いとして睡眠段階を判定する。例えば、睡眠段階を判定するための判定しきい値を予め設定しておき、そのしきい値と三角形の辺の傾きのバラツキ量を大小比較することにより、睡眠段階を判定すればよい。
【０１５３】
次に、本実施形態に係る睡眠判定装置の動作を説明する。
【０１５４】
まず、図１において、心拍検出部２により被験者の心拍信号が検出される。そして、この心拍信号に基づいて図４の睡眠判定処理により被験者の睡眠段階の判定が行われる。
【０１５５】
そして、図４に示すように、心拍信号の読み込み処理が行われる（Ｓ１０）。これらの処理は、第一実施形態に説明した処理と同様に行われる。
【０１５６】
そして、差分演算処理、リターンマップ生成処理が順次行われる（Ｓ１２、Ｓ１４）。この差分演算処理及びリターンマップ生成処理は、第五実施形態に係る睡眠判定装置と同様に行えばよい。
【０１５７】
そして、Ｓ１６に移行し、特徴量演算処理が行われる。特徴量演算処理は、Ｓ１４のリターンマップ生成処理により生成されたリターンマップの特徴量を演算する処理である。リターンマップの特徴量として、リターンマップにおいて連続して演算される三つの座標点により形成される三角形の辺の傾きのバラツキ量が演算される。演算された三角形の辺の傾きのバラツキ量は、データベース３５に記録される。三角形の辺の傾きのバラツキ量は、心拍信号の入力されるごとに演算され、随時データベースに記録される。
【０１５８】
そして、Ｓ１８に移行し、睡眠段階判定処理が行われる。睡眠段階判断処理は、Ｓ１６にて演算されたリターンマップの特徴量に基づいて被験者の睡眠段階を判定する処理である。例えば、Ｓ１６にて演算された三角形の辺の傾きのバラツキ量が小さいほど睡眠状態が深いとして睡眠段階が判定される。そして、Ｓ１８の睡眠段階判定処理を終えたら、一連の制御処理を終了する。
【０１５９】
以上のように、本実施形態に係る睡眠判定装置によれば、被験者の心拍周期の差分値を座標値とする座標点を心拍ごとに繰り返し演算してリターンマップを生成し、そのリターンマップの特徴量に基づいて被験者の睡眠段階を判定することにより、睡眠段階に応じて変化する心拍周期の差分値の特徴量を用いて被験者の睡眠段階を精度よく判定することができる。
【０１６０】
また、本実施形態に係る睡眠判定装置においては、リターンマップにおいて連続して演算される三つの座標点により形成される三角形の辺の傾きのバラツキ量を演算し、そのバラツキ量が小さいほど睡眠状態が深いとして睡眠段階を判定することにより、睡眠段階によって変化する心拍周期の差分値のリターンマップの座標点の変動状態を用いて、被験者の睡眠段階を精度よく判定することができる。
【０１６１】
すなわち、図１６〜２０に示されるように、睡眠段階がＳ０からＳ４へ深くなるに連れて、連続する三つの座標点により形成される三角形の辺の傾きのバラツキが小さくなることが分かる。
【０１６２】
このため、被験者の心拍周期の差分値についてリターンマップを生成し座標点により形成される三角形の辺の傾きのバラツキ量を演算することにより、そのバラツキ量に応じて睡眠段階を精度よく判定することが可能となる。
【０１６３】
なお、上述した各実施形態は本発明に係る睡眠判定装置の実施形態を説明したものであり、本発明に係る睡眠判定装置は本実施形態に記載したものに限定されるものではない。本発明に係る睡眠判定装置は、各請求項に記載した要旨を変更しないように実施形態に係る睡眠判定装置を変形し、又は他のものに適用したものであってもよい。
【０１６４】
例えば、実施形態においては、被験者が車両の運転者である場合について説明したが、被験者が車両の運転者以外の者であってもよく、車両に乗車していないものであってもよい。
【図面の簡単な説明】
【０１６５】
【図１】本発明の第一実施形態に係る睡眠判定装置の構成概要図である。
【図２】心拍信号の説明図である。
【図３】図１の睡眠判定装置により生成されるリターンマップの説明図である。
【図４】図１の睡眠判定装置の動作を示すフローチャートである。
【図５】睡眠段階Ｓ０における心拍周期のリターンマップの説明図である。
【図６】睡眠段階Ｓ１における心拍周期のリターンマップの説明図である。
【図７】睡眠段階Ｓ２における心拍周期のリターンマップの説明図である。
【図８】睡眠段階Ｓ３における心拍周期のリターンマップの説明図である。
【図９】睡眠段階Ｓ４における心拍周期のリターンマップの説明図である。
【図１０】本発明の第二実施形態に係る睡眠判定装置により生成されるヒストグラムの説明図である。
【図１１】睡眠段階ごとのヒストグラムについての説明図である。
【図１２】本発明の第三実施形態に係る睡眠判定装置の睡眠段階判定の説明図である。
【図１３】本発明の第四実施形態に係る睡眠判定装置の睡眠段階判定の説明図である。
【図１４】本発明の第五実施形態に係る睡眠判定装置における心拍周期の差分値演算の説明図である。
【図１５】本発明の第五実施形態に係る睡眠判定装置により生成されるリターンマップの説明図である。
【図１６】睡眠段階Ｓ０における心拍周期の差分値のリターンマップの説明図である。
【図１７】睡眠段階Ｓ１における心拍周期の差分値のリターンマップの説明図である。
【図１８】睡眠段階Ｓ２における心拍周期の差分値のリターンマップの説明図である。
【図１９】睡眠段階Ｓ３における心拍周期の差分値のリターンマップの説明図である。
【図２０】睡眠段階Ｓ４における心拍周期の差分値のリターンマップの説明図である。
【図２１】本発明の第六実施形態に係る睡眠判定装置におけるリターンマップの座標点から座標点の傾きについてヒストグラムの説明図である。
【符号の説明】
【０１６６】
１…睡眠判定装置、２…心拍検出部、３…ＥＣＵ、４…出力部、３１…差分演算部、３２…リターンマップ生成部、３３…ヒストグラム生成部、３４…睡眠判定部、３５…データベース。
【図１】
2010148605.tif 000004
【図２】
2010148605.tif 000005
【図３】
2010148605.tif 000006
【図４】
2010148605.tif 000007
【図５】
2010148605.tif 000008
【図６】
2010148605.tif 000009
【図７】
2010148605.tif 000010
【図８】
2010148605.tif 000011
【図９】
2010148605.tif 000012
【図１０】
2010148605.tif 000013
【図１１】
2010148605.tif 000014
【図１２】
2010148605.tif 000015
【図１３】
2010148605.tif 000016
【図１４】
2010148605.tif 000017
【図１５】
2010148605.tif 000018
【図１６】
2010148605.tif 000019
【図１７】
2010148605.tif 000020
【図１８】
2010148605.tif 000021
【図１９】
2010148605.tif 000022
【図２０】
2010148605.tif 000023
【図２１】
2010148605.tif 000024