(57)【要約】      （修正有）【課題】ウルトラディアンリズムのサイクル毎におけるδ波のパワーを推定することにより、使用者にとって少ない負担で使用者の睡眠の質を推定できる装置を提供する。【解決手段】睡眠状態推定装置１は、使用者の睡眠時の生体情報を取得する生体情報取得部３と、生体情報に基づいて、使用者の睡眠段階を推定する睡眠段階推定部１０と、睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割部１１と、睡眠段階の推定結果に基づいて、サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定部１２と、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定するパワー推定部１３と、サイクル毎におけるパワーの推定結果に基づいて、使用者の睡眠の質を推定する質推定部１４と、を備える。【選択図】図１
2015097713.tif 000002
【特許請求の範囲】
【請求項１】
使用者の睡眠の質を推定する睡眠状態推定装置であって、
前記使用者の睡眠時の生体情報を取得する生体情報取得部と、
前記生体情報に基づいて、前記使用者の睡眠段階を推定する睡眠段階推定部と、
前記睡眠段階の推定結果に基づいて、前記使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割部と、
前記睡眠段階の推定結果に基づいて、前記サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定部と、
前記サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、前記サイクル毎における前記使用者のδ波のパワーを推定するパワー推定部と、
前記サイクル毎における前記パワーの推定結果に基づいて、前記使用者の睡眠の質を推定する質推定部と、
を備え、
前記パワー推定部は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーが前記パワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーが前記パワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、前記サイクル毎におけるノンレム睡眠の継続時間の推定結果から前記サイクル毎における前記使用者のδ波のパワーを推定する、睡眠状態推定装置。
【請求項２】
前記パワー推定部は、前記使用者の入眠からの前記サイクルの順番に応じて、それぞれ設定された前記δ波のパワー最大値を用いて、当該サイクルにおける前記使用者のδ波のパワーを推定する、請求項１に記載の睡眠状態推定装置。
【請求項３】
前記パワー推定部は、前記使用者の入眠からの前記サイクルの順番に応じて、それぞれ設定された前記開始時定数を用いて、当該サイクルにおける前記使用者のδ波のパワーを推定する、請求項１又は２に記載の睡眠状態推定装置。
【請求項４】
前記パワー推定部は、前記使用者の入眠からの前記サイクルの順番に応じて、それぞれ設定された前記終了時定数を用いて、当該サイクルにおける前記使用者のδ波のパワーを推定する、請求項１〜３のうち何れか一項に記載の睡眠状態推定装置。
【請求項５】
使用者の睡眠の質を推定する睡眠状態推定方法であって、
前記使用者の睡眠時の生体情報を取得する生体情報取得工程と、
前記生体情報に基づいて、前記使用者の睡眠段階を推定する睡眠段階推定工程と、
前記睡眠段階の推定結果に基づいて、前記使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割工程と、
前記睡眠段階の推定結果に基づいて、前記サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定工程と、
前記サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、前記サイクル毎における前記使用者のδ波のパワーを推定するパワー推定工程と、
前記サイクル毎における前記パワーの推定結果に基づいて、前記使用者の睡眠の質を推定する質推定工程と、
を備え、
前記パワー推定工程は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーが前記パワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーが前記パワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、前記サイクル毎におけるノンレム睡眠の継続時間の推定結果から前記サイクル毎における前記使用者のδ波のパワーを推定する、睡眠状態推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、使用者の睡眠の質を推定する睡眠状態推定装置及び睡眠状態推定方法に関する。
【背景技術】
【０００２】
睡眠の質は、睡眠時の脳波と相関があることが知られている。脳波は、α波、β波、γ波、δ波に分類され、δ波は眠りの深い睡眠段階において現われる傾向がある。具体的に、睡眠段階は、脳が覚醒状態にあるレム睡眠と脳が休息している状態にあるノンレム睡眠とに分けられ、ノンレム睡眠においてδ波の出現確率が高くなる。そこで、睡眠時におけるδ波等の脳波の測定結果に基づいて睡眠の質を推定することが従来から行われている。
【０００３】
例えば、特開２０１３−１２１４８９号公報には、使用者の前額部に配置される複数個の電極を備え、これらの電極により測定した脳波に基づいて使用者の睡眠状態の検出を行う睡眠状態検出装置が示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−１２１４８９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上述した従来の睡眠状態検出装置のように使用者の頭部へ電極を配置することは、使用者にとって負担であり、使用者の睡眠を妨げるおそれがある。このため、従来の睡眠状態検出装置を日常的に使用することは現実的ではないという問題があった。
【０００６】
そこで、本発明は、ウルトラディアンリズムのサイクル毎におけるδ波のパワーを推定することにより、使用者にとって少ない負担で使用者の睡眠の質を推定できる睡眠状態推定装置及び睡眠状態推定方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明者等は、複数人の被験者の睡眠時の脳波（単位：μＶ）測定を十数日にわたって行う中で、睡眠時の脳活動の休息を意味するδ波が約９０分のウルトラディアンリズムのサイクルに基づいて変化することに着目した。本発明者等は、研究の結果、ウルトラディアンリズムのサイクル毎におけるノンレム睡眠時のδ波のパワー（単位：μＶ２／Ｈｚ）が人によらず同様のダイナミクスを示すことを新たに見いだし、このダイナミクスを利用してサイクル毎のノンレム睡眠の継続時間からδ波のパワーを推定することを想到した。
【０００８】
すなわち、本発明の一態様は、使用者の睡眠の質を推定する睡眠状態推定装置であって、使用者の睡眠時の生体情報を取得する生体情報取得部と、生体情報に基づいて、使用者の睡眠段階を推定する睡眠段階推定部と、睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割部と、睡眠段階の推定結果に基づいて、サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定部と、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定するパワー推定部と、サイクル毎におけるパワーの推定結果に基づいて、使用者の睡眠の質を推定する質推定部と、を備え、パワー推定部は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定する。なお、上記生体情報に脳波は含まれない。
【０００９】
この睡眠状態推定装置は、睡眠時の呼吸、心電（心拍）、血圧、皮膚電気活動などの生体情報に基づいて使用者の睡眠段階を推定することで、当該睡眠段階の推定結果から睡眠時間をウルトラディアンリズムのサイクル毎に分け、当該睡眠段階の推定結果からサイクル毎におけるノンレム睡眠の継続時間を推定することができる。そして、この睡眠状態推定装置は、サイクル毎におけるノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらず同様の傾向を示すことから、この傾向を利用してサイクル毎のノンレム睡眠の継続時間からサイクル毎のδ波のパワーを推定することができる。すなわち、この睡眠状態推定装置は、サイクル毎におけるノンレム睡眠の開始からδ波のパワーが増加してδ波のパワー最大値に至り、その後にδ波のパワーが減少してノンレム睡眠が終了するという傾向を利用して、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定することができる。従って、この睡眠状態推定装置は、睡眠時の使用者の頭部に配置する電極等を用いてδ波を直接測定することなく、使用者にとって少ない負担で睡眠時のδ波のパワーを推定することができ、δ波のパワーの推定結果に基づいて使用者の睡眠の質を推定することができる。
【００１０】
一実施形態において、上記パワー推定部は、使用者の入眠からのサイクルの順番に応じて、それぞれ設定されたδ波のパワー最大値を用いて、当該サイクルにおける使用者のδ波のパワーを推定してもよい。
この実施形態に係る睡眠状態推定装置は、使用者の入眠からのサイクルの順番により、当該サイクルにおけるノンレム睡眠の継続時間とδ波のパワーの関係性が変化することから、サイクルの順番に応じてそれぞれ設定されたδ波のパワー最大値を用いることで、δ波のパワー最大値として全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーの推定精度の向上を図ることができる。
【００１１】
一実施形態において、上記パワー推定部は、使用者の入眠からのサイクルの順番に応じて、それぞれ設定された開始時定数を用いて、当該サイクルにおける使用者のδ波のパワーを推定してもよい。
この実施形態に係る睡眠状態推定装置は、使用者の入眠からのサイクルの順番により、当該サイクルにおけるノンレム睡眠の継続時間とδ波のパワーの関係性が変化することから、サイクルの順番に応じてそれぞれ設定された開始時定数を用いることで、開始時定数として全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーの推定精度の向上を図ることができる。
【００１２】
一実施形態において、上記パワー推定部は、使用者の入眠からのサイクルの順番に応じて、それぞれ設定された終了時定数を用いて、当該サイクルにおける使用者のδ波のパワーを推定してもよい。
この実施形態に係る睡眠状態推定装置は、使用者の入眠からのサイクルの順番により、当該サイクルにおけるノンレム睡眠の継続時間とδ波のパワーの関係性が変化することから、サイクルの順番に応じてそれぞれ設定された終了時定数を用いることで終了時定数として全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーの推定精度の向上を図ることができる。
【００１３】
本発明の他の態様は、使用者の睡眠の質を推定する睡眠状態推定方法であって、使用者の睡眠時の生体情報を取得する生体情報取得工程と、生体情報に基づいて、使用者の睡眠段階を推定する睡眠段階推定工程と、睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割工程と、睡眠段階の推定結果に基づいて、サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定工程と、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定するパワー推定工程と、サイクル毎におけるパワーの推定結果に基づいて、使用者の睡眠の質を推定する質推定工程と、を備え、パワー推定工程は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定する。
【００１４】
この睡眠状態推定方法によれば、睡眠時の呼吸、心電（心拍）、血圧、皮膚電気活動などの生体情報に基づいて使用者の睡眠段階を推定することで、当該睡眠段階の推定結果から睡眠時間をウルトラディアンリズムのサイクル毎に分け、当該睡眠段階の推定結果からサイクル毎におけるノンレム睡眠の継続時間を推定することができる。そして、この睡眠状態推定方法では、サイクル毎におけるノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらず同様のダイナミクスを示すことから、このダイナミクスを利用してサイクル毎のノンレム睡眠の継続時間からサイクル毎のδ波のパワーを推定することができる。すなわち、この睡眠状態推定方法では、サイクル毎におけるノンレム睡眠の開始からδ波のパワーが増加してδ波のパワー最大値に至り、その後にδ波のパワーが減少してノンレム睡眠が終了するという傾向を利用して、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定することができる。従って、この睡眠状態推定方法によれば、睡眠時の使用者の頭部に配置する電極等を用いてδ波を直接測定することなく、使用者にとって少ない負担で睡眠時のδ波のパワーを推定することができ、δ波のパワーの推定結果に基づいて使用者の睡眠の質を推定することができる。
【発明の効果】
【００１５】
本発明の一態様に係る睡眠状態推定装置及び睡眠状態推定方法によれば、ウルトラディアンリズムのサイクル毎におけるδ波のパワーを推定することにより、使用者にとって少ない負担で使用者の睡眠の質を推定できる。
【図面の簡単な説明】
【００１６】
【図１】本実施形態に係る睡眠状態推定装置を示すブロック図である。
【図２】使用者の睡眠段階の推移の一例を示す図である。
【図３】パワー推定部の入力及び出力を説明するためのブロック図である。
【図４】（ａ）サイクル毎におけるδ波のパワーの実測データを時間軸について正規化したグラフである。（ｂ）図４（ａ）と異なる被験者について、サイクル毎におけるδ波のパワーの実測データを時間軸について正規化したグラフである。
【図５】（ａ）サイクル１におけるδ波のパワー実測値を示すグラフである。（ｂ）図５（ａ）に示すフィルタ後のサイクル１のグラフである。
【図６】（ａ）サイクル毎におけるδ波のパワーの実測値にフィルタを掛けたグラフである。（ｂ）サイクル毎のδ波のパワーモデルを示すグラフである。
【図７】本実施形態に係る睡眠状態推定方法を示すフローチャートである。
【発明を実施するための形態】
【００１７】
以下、本発明の好適な実施形態について、図面を参照して詳細に説明する。
【００１８】
図１に示す睡眠状態推定装置１は、睡眠時の使用者の生体情報を利用してウルトラディアンリズムのサイクル毎におけるδ波（脳波の一種）のパワーを推定することで、睡眠の質の推定を行うものである。使用者の生体情報には、呼吸、心電（心拍）、血圧、体温、皮膚電気活動に関する情報などが含まれる。なお、睡眠状態推定装置１は、使用者の脳波を直接測定するものではないため、使用者の生体情報に脳波は含まれない。
【００１９】
ウルトラディアンリズムとは、人体の生体リズムの一つであり、睡眠時におけるレム睡眠とノンレム睡眠のサイクルである。レム睡眠とは脳が覚醒状態にある浅い睡眠段階を意味し、ノンレム睡眠とは脳が休息している状態にある深い睡眠段階を意味する。ウルトラディアンリズムの周期は、個人差もあるが、一般的に約９０分である。例えば、使用者の一晩の睡眠時間が６時間（３６０分）である場合には、９０分周期のウルトラディアンリズムのサイクルが一晩で４回繰り返される。以後の説明において、特に断りがない限り、サイクルという言葉はウルトラディアンリズムのサイクルを意味する。
【００２０】
δ波は、約１〜３Ｈｚの周波数を持つ脳波の一種であり、単位は電圧と同じμＶである。δ波は、脳波のうち最もゆっくりした波であり、深い睡眠段階（ノンレム睡眠）において出現確率が高いことが知られている。
【００２１】
ウルトラディアンリズムのサイクル毎におけるδ波（脳波の一種）のパワーとは、サイクル毎におけるδ波のパワーの時間積分であるパワースペクトルに相当する。このパワーの単位は、μＶ２／Ｈｚである。
【００２２】
［睡眠状態推定装置の構成］
本実施形態における睡眠状態推定装置１は、装置全体を統括的に制御する演算部２を備えている。演算部２は、ＣＰＵ[Central  Processing  Unit]、ＲＯＭ[Read  Only  Memory]、ＲＡＭ[Random  Access  Memory]などを含む電子制御ユニットである。この演算部２は、生体情報取得部３及び表示部４と接続されている。
【００２３】
生体情報取得部３は、睡眠時の使用者の生体情報を取得可能なバイタルセンサから構成されている。具体的に、生体情報取得部３は、リストバンド型や貼付け型のウェアラブルセンサを含んでもよい。ウェアラブルセンサは、使用者の負担軽減の観点からワイヤレスとすることができる。ウェアラブルセンサは、使用者の呼吸、心電（心拍）、血圧、体温、皮膚電気活動等の生体情報を検出する。なお、本実施形態に係る生体情報取得部３は、使用者の脳波を直接測定する脳波計を含まない。
【００２４】
表示部４は、睡眠状態推定装置１における各種の検出結果及び推定結果を表示するためのディスプレイである。表示部４に表示される内容としては、使用者の入眠から起床までの睡眠時間、中途覚醒の回数及び合計時間、寝返り回数、レム睡眠の合計時間、ノンレム睡眠の合計時間、ウルトラディアンリズムのサイクル数等が挙げられる。
【００２５】
続いて、演算部２の構成について説明する。演算部２は、睡眠段階推定部１０、睡眠時間分割部１１、ノンレム睡眠時間推定部１２、パワー推定部１３、質推定部１４、及び記憶部１５を有している。記憶部１５は、過去の睡眠段階の比率や過去の睡眠時間などの睡眠履歴を個人別に記憶しているデータベースである。また、記憶部１５は、後述するサイクル毎のδ波のパワーモデル等の推定に用いるデータを記憶している。
【００２６】
睡眠段階推定部１０は、生体情報取得部３が取得した睡眠時の使用者の生体情報に基づいて、現在の使用者の睡眠段階の推定を行う。睡眠段階は、レム睡眠とノンレム睡眠に大きく分けられる。睡眠段階推定部１０は、ノンレム睡眠を更に複数の段階に分けて推定を行ってもよい。
【００２７】
睡眠段階推定部１０は、例えば、使用者の生体情報から特徴量を算出し、当該特徴量に基づいて睡眠段階を推定する。特徴量とは、例えば生体情報としての心拍の検出データから抽出される単位時間当たりの心拍数、心拍の周期のゆらぎ、心拍の強さなどである。睡眠段階推定部１０は、使用者の生体情報から様々な観点の特徴量を算出し、これらの特徴量に基づいて使用者の睡眠段階を多元的に推定する。
【００２８】
睡眠時間分割部１１は、睡眠段階推定部１０による睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分割する。ここで、図２は、使用者の睡眠段階の推定結果の一例を示す図である。
【００２９】
図２に示されるように、使用者は、入眠によってレム睡眠に入った後、周期的にレム睡眠とノンレム睡眠とを繰り返す。睡眠時間分割部１１は、例えば、睡眠段階の推定結果からレム睡眠とノンレム睡眠との繰り返しを検出することで、睡眠時間をウルトラディアンリズムのサイクル毎に分ける。使用者の入眠から起床までの睡眠時間は、３〜４のサイクルに分けられることが通常である。図２では、使用者の入眠から起床までの睡眠時間が４つのサイクルに分けられている。
【００３０】
以下、図２に示されるように、使用者の入眠から１番目のサイクルをサイクル１、２番目のサイクルをサイクル２、３番目のサイクルをサイクル３、４番目のサイクルをサイクル４として説明に用いる。
【００３１】
なお、睡眠時間分割部１１は、ウルトラディアンリズムのサイクルを９０分と仮定して睡眠時間を９０分毎に分割した後、レム睡眠とノンレム睡眠との繰り返しを検出して睡眠時間の分割位置を補正する態様であってもよい。
【００３２】
図１に示すノンレム睡眠時間推定部１２は、睡眠段階推定部１０による睡眠段階の推定結果に基づいて、睡眠時間分割部１１によって分けられたサイクル毎におけるノンレム睡眠の継続時間を推定する。具体的に、図２において、ノンレム睡眠時間推定部１２は、睡眠段階の推定結果からサイクル１におけるノンレム睡眠の継続時間Ｔ１を推定する。同様に、ノンレム睡眠時間推定部１２は、サイクル２におけるノンレム睡眠の継続時間Ｔ２を推定し、サイクル３におけるノンレム睡眠の継続時間Ｔ３を推定する。また、ノンレム睡眠時間推定部１２は、サイクル４におけるノンレム睡眠の継続時間Ｔ４を推定する。
【００３３】
図１及び図３に示すパワー推定部１３は、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定する。パワー推定部１３は、予め設定されたδ波のパワーモデルを利用して、サイクル毎における使用者のδ波のパワーを推定する。
【００３４】
以下、δ波のパワーモデルについて説明する。δ波のパワーモデルとは、サイクル毎におけるノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらず共通する傾向を示すことに着目して作成されたものである。このようなδ波のパワーモデルを用いることで、あるサイクルにおけるノンレム睡眠の継続時間から当該サイクルにおけるδ波のパワーを推定することができる。δ波のパワーモデルのデータは、例えば記憶部１５に記憶されている。
【００３５】
また、δ波のパワーモデルは、使用者の入眠からのサイクルの順番に応じて複数のモデルが予め設定されている。これは、サイクルの順番が同じ場合にノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらずに同様のダイナミクスを示すためである。なお、図３には、サイクル１〜４の各々のδ波のパワーモデルが示されているが、サイクル５以降のパワーモデルについて予め用意していてもよい。
【００３６】
サイクル１〜４のδ波のパワーモデルは、下記の式（１）〜（４）に示すモデルパラメータとして表すことができる。式（１）は、サイクル１に対応するモデルパラメータであり、式（２）はサイクル１に対応するモデルパラメータである。式（３）は、サイクル３に対応するモデルパラメータであり、式（４）はサイクル４に対応するモデルパラメータである。
【数１】
2015097713.tif 000003
【００３７】
上記の式（１）において、Ｐ１は、サイクル１におけるδ波のパワーである。Ｐmax１は、サイクル１において予め設定されたδ波のパワー最大値である。τｓ１は、サイクル１において予め設定された開始時定数である。τｅ１は、サイクル１において予め設定された終了時定数である。Ｔ１は、サイクル１におけるノンレム睡眠の継続時間である。ｔ１は、サイクル１が始まってからの時間である。
【００３８】
同様に、式（２）〜式（４）において、Ｐ２〜Ｐ４は、サイクル２〜４のそれぞれにおけるδ波のパワーである。Ｐmax２〜Ｐmax４は、サイクル２〜４のそれぞれにおいて予め設定されたδ波のパワー最大値である。τｓ２〜τｓ４は、サイクル２〜４のそれぞれにおいて予め設定された開始時定数である。τｅ２〜τｅ４は、サイクル２〜４のそれぞれにおいて予め設定された終了時定数である。Ｔ２〜Ｔ４は、サイクル２〜４のそれぞれにおけるノンレム睡眠の継続時間である。ｔ2〜ｔ4は、サイクル２〜４のそれぞれにおいてサイクルが始まってからの時間である。
【００３９】
以下、サイクル１におけるδ波のパワー最大値Ｐmax１、開始時定数τｓ１、終了時定数τｅ１について説明する。
【００４０】
まず、図４（ａ）は、サイクル毎におけるδ波のパワーの実測値をサイクルの時間的長さについて正規化したグラフである。図４（ｂ）は、図４（ａ）と異なる被験者について、サイクル毎におけるδ波のパワーの実測データを時間軸について正規化したグラフである。
図４（ａ）及び図４（ｂ）において、縦軸はδ波のパワー、横軸は時間を示している。図４（ａ）及び図４（ｂ）に示すδ波のパワーの実測値は、被験者のδ波を脳波計で直接測定して得たものである。サイクル毎の時間的長さのバラツキを無くすため、時間軸について正規化を行った後、サイクル毎におけるδ波のパワーの実測値を重ねて表している。
【００４１】
図４（ａ）及び図４（ｂ）に示されるように、δ波のパワーのダイナミクスは被験者ごとに差異はあるがパワーの最大値（上記の式（１）〜式（４）におけるＰmax１〜Ｐmax４）がサイクル１からサイクル４にかけて徐々に小さくなること、δ波のパワーの立ち上がり（上記の式（１）〜式（４）におけるτｓ１〜τｓ４）がサイクル１からサイクル４にかけて遅くなることに関して共通していることが分かる。このような特徴は、ここでは示していない他の被験者でも確認されたことから、δ波のパワー最大値とδ波のパワーの立ち上がり等をパラメータとする式（１）〜式（４）により、ある程度の精度でδ波のダイナミクスを推定できると考えられる。
【００４２】
図５（ａ）は、サイクル１におけるδ波のパワーの実測値を示すグラフである。なお、図５（ａ）に示す実測値は、図４（ａ）又は図４（ｂ）に示すサイクル１のδ波のパワーの実測値と同一ではない。図５（ａ）において、サイクル１の実測値をフィルタにより鈍らせたグラフを太線で示す。
【００４３】
図５（ｂ）に、図５（ａ）に示すフィルタ後のサイクル１のグラフを示す。図５（ｂ）では、図５（ａ）よりも縦軸の寸法を拡大している。図５（ｂ）において、予め設定されたδ波のパワー最大値Ｐmax１を示す。
【００４４】
まず、δ波のパワー最大値Ｐmax１について説明する。δ波のパワー最大値Ｐmax１は、例えば、実験により統計的に求められた値であり、睡眠状態推定装置１に予め設定された設定値である。δ波のパワー最大値Ｐmax１は、例えば、被験者ごとに十数日以上にわたって測定された睡眠時のサイクル１におけるδ波のパワーの実測値の最大値の平均とすることができる。同様に、他のサイクル２〜４におけるδ波のパワー最大値Ｐmax２〜Ｐmax４も、被験者ごとに測定されたサイクル２〜４におけるδ波のパワーの実測値の最大値の平均とすることができる。
【００４５】
なお、δ波のパワー最大値Ｐmax１は、予め設定された複数の設定値の中から性別や年齢等に応じて使用者が選択する態様であってもよい。また、δ波のパワー最大値Ｐmax１は、使用者の過去の睡眠履歴等の個人データを利用して予め設定された値を次第に補正する態様であってもよい。サイクル２〜４におけるδ波のパワー最大値Ｐmax２〜Ｐmax４においても同様である。
【００４６】
ここで、図４（ａ）及び図４（ｂ）に示されるように、使用者の入眠からのサイクルの順番が早いほどδ波のパワーが大きく現われ、サイクルの順番が遅いほどδ波のパワーが小さくなることが経験的に認められる。このようにサイクルの順番に応じて、ノンレム睡眠の継続時間とδ波のパワーとの関係性が変化することから、δ波のパワー最大値Ｐmax１〜Ｐmax４をそれぞれ適切な設定値とすることができる。例えば、δ波のパワー最大値Ｐmax１〜Ｐmax４は、サイクルの順番が遅いほど小さな値とすることができる。この場合、サイクルの順番に応じてそれぞれ設定されたδ波のパワー最大値Ｐmax１〜Ｐmax４を用いることで、δ波のパワー最大値として全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーＰ１〜Ｐ４の推定精度を向上させることができる。
【００４７】
なお、必ずしもδ波のパワー最大値Ｐmax１〜Ｐmax４の全てを異なる値とする必要はない。例えば、サイクル１，２におけるδ波のパワーのグラフは、サイクル３，４と比べて、統計的に同様のダイナミクスが認められるため、δ波のパワーの最大値Ｐmax１及びδ波のパワー最大値Ｐmax２を同じ値としてもよい。同様に、δ波のパワーの最大値Ｐmax３及びδ波のパワー最大値Ｐmax４を同じ値としてもよい。また、δ波のパワー最大値Ｐmax１〜Ｐmax４の全てを同じ値とすることもできる。
【００４８】
次に、開始時定数τｓ１及び終了時定数τｅ１について説明する。開始時定数τｓ１及び終了時定数τｅ１は、δ波のパワー最大値Ｐmax１と同様に、睡眠状態推定装置１に予め設定された設定値である。
【００４９】
開始時定数τｓ１は、図５（ｂ）において、予め設定されたδ波のパワー最大値Ｐmax１に１−ｅｘｐ−１を乗じた値（Ｐmax１×（１−ｅｘｐ−１））に対応する線Ｌｐとフィルタ後のサイクル１のグラフとの交点のうち、サイクル１の開始側の交点Ｃ１に対応する時間に相当する。開始時定数τｓ１は、サイクル１において、ノンレム睡眠の開始からδ波のパワーがパワー最大値Ｐmax１に至るまでの時間に関する統計結果に基づいて予め設定された値である。具体的には、開始時定数τｓ１は、複数人の被験者から十数日以上にわたって測定された睡眠時のサイクル１において、ノンレム睡眠の開始からδ波のパワーがパワー最大値Ｐmax１に至るまでの時間の平均とすることができる。
【００５０】
一方、終了時定数τｅ１は、図５（ｂ）において、線Ｌｐとフィルタ後のサイクル１のグラフとの交点のうち、サイクル１の終了側の交点Ｃ２に対応する時間に相当する。終了時定数τｅ１は、サイクル１において、δ波のパワーがパワー最大値Ｐmax１から減少してノンレム睡眠の終了に至るまでの時間に関する統計結果に基づいて予め設定された値である。具体的には、終了時定数τｅ１は、複数人の被験者から十数日以上にわたって測定された睡眠時のサイクル１において、ノンレム睡眠の開始からδ波のパワーがパワー最大値Ｐmax１に至るまでの時間の平均とすることができる。
【００５１】
なお、開始時定数τｓ１及び終了時定数τｅ１も、予め設定された複数の設定値の中から性別や年齢等に応じて使用者が選択する態様であってもよい。また、開始時定数τｓ１及び終了時定数τｅ１は、使用者の過去の睡眠履歴等の個人データを利用して予め設定された値を次第に補正する態様であってもよい。以上、サイクル１における開始時定数τｓ１及び終了時定数τｅ１について説明したが、他のサイクル２〜４における開始時定数τｓ２〜τｓ４及び終了時定数τｅ２〜τｅ４についても同様である。
【００５２】
また、開始時定数τｓ１〜τｓ４及び終了時定数τｅ１〜τｅ４においても、使用者の入眠からのサイクルの順番に応じてそれぞれ適切な設定値とすることで、全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーＰ１〜Ｐ４の推定精度を向上させることができる。
【００５３】
例えば、図４に示されるように、サイクル１，２においてδ波のパワーが最大値付近に至るまでの時間は、サイクル３，４においてδ波のパワーが最大値付近に至るまでの時間より短いことが経験的に認められる。このため、開始時定数τｓ１，τｓ２を開始時定数τｓ３，τｓ４より小さな値として予め設定することができる。なお、必ずしも開始時定数τｓ１〜τｓ４の全てを異なる値とする必要はない。例えば、開始時定数τｓ１，τｓ２を同じ値としてもよく、開始時定数τｓ３，τｓ４を同じ値としてもよい。また、開始時定数τｓ１〜τｓ４の全てを同じ値とすることもできる。
【００５４】
同様に、図４に示されるように、サイクル１，２においてδ波のパワーが減少を始めてノンレム睡眠が終了するまでの時間は、サイクル３，４においてδ波のパワーが減少を始めてノンレム睡眠が終了するまでの時間より長いことが経験的に認められる。このため、終了時定数τｅ１，τｅ２を終了時定数τｅ３，τｅ４より大きな値として予め設定することができる。なお、必ずしも終了時定数τｅ１〜τｅ４の全てを異なる値とする必要はない。例えば、終了時定数τｅ１，τｅ２を同じ値としてもよく、開始時定数τｅ３，τｅ４を同じ値としてもよい。また、終了時定数τｅ１〜τｅ４の全てを同じ値とすることもできる。
【００５５】
パワー推定部１３は、ノンレム睡眠時間推定部１２におけるサイクル毎のノンレム睡眠の継続時間Ｔ１〜Ｔ４と、上述したδ波のパワー最大値Ｐmax１〜Ｐmax４、開始時定数τｓ１〜τｓ４、及び終了時定数τｅ１〜τｅ４とを用いて、上記の式（１）〜式（４）からサイクル毎のδ波のパワーＰ１〜Ｐ４を推定する。
【００５６】
ここで、図６（ａ）は、サイクル毎におけるδ波のパワーの実測値にフィルタを掛けたグラフである。図６（ｂ）は、サイクル毎のδ波のパワーモデルを示すグラフである。図６（ａ）及び図６（ｂ）に示されるように、上記式（１）〜式（４）をモデルパラメータとするサイクル毎のδ波のパワーモデルは、実測値と近い傾向を再現している。パワー推定部１３は、これらのδ波のパワーモデルを用いることにより、サイクル毎のノンレム睡眠の継続時間からサイクル毎のδ波のパワーを十分な精度で推定することができる。
【００５７】
続いて、図１及び図３に示す質推定部１４は、パワー推定部１３によるサイクル毎のδ波のパワーに基づいて、使用者の睡眠の質を推定する。質推定部１４は、例えば、サイクル毎のδ波のパワーＰ１〜Ｐ４の積算値を用いて使用者の睡眠の質を段階的に推定する。
【００５８】
また、質推定部１４は、使用者の生体情報から得られる使用者の睡眠時間及び寝返り回数等の情報を睡眠の質の推定に用いることができる。更に、質推定部１４は、睡眠段階の推定結果から得られる中途覚醒の回数、中途覚醒の合計時間、睡眠時間に占めるレム睡眠の合計時間、睡眠時間に含まれるウルトラディアンリズムのサイクル数等に基づいて、使用者の睡眠の質を推定することができる。質推定部１４は、例えば、睡眠時間におけるレム睡眠とノンレム睡眠の比率に応じて、当該比率が予め設定された最適値に近いほど睡眠の質が高いと推定する。
【００５９】
［睡眠状態推定方法］
次に、本実施形態に係る睡眠状態推定装置１を用いた睡眠状態推定方法について図面を参照して説明する。図７は、本実施形態に係る睡眠状態推定方法を示すフローチャートである。
【００６０】
図７に示されるように、睡眠状態推定装置１は、ステップＳ１において生体情報取得部３による使用者の生体情報の取得を行う。生体情報取得部３は、使用者の各種の生体情報を所定のサンプリング周期で取得する（生体情報取得工程）。
【００６１】
次に、ステップＳ２において、睡眠段階推定部１０は、生体情報取得部３が取得した睡眠時の使用者の生体情報に基づいて現在の使用者の睡眠段階の推定を行う（睡眠段階推定工程）。その後、ステップＳ３において、睡眠段階推定部１０は、使用者の生体情報に基づいて使用者が起床したか否かの判定を行う。使用者が起床したか否かの判定には、周知の様々な方法を採用することができる。睡眠段階推定部１０は、使用者が起床していないと判定した場合、再びステップＳ１に戻り、新たな生体情報に基づいて睡眠段階の推定を行う。睡眠段階推定部１０は、使用者が起床したと判定した場合、ステップＳ４に移行する。
【００６２】
ステップＳ４において、睡眠時間分割部１１は、睡眠段階推定部１０による睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分割する（睡眠時間分割工程）。
【００６３】
その後、ステップＳ５において、ノンレム睡眠時間推定部１２は、睡眠段階推定部１０による睡眠段階の推定結果に基づいて、睡眠時間分割部１１によって分けられたサイクル毎におけるノンレム睡眠の継続時間を推定する（ノンレム睡眠時間推定工程）。
【００６４】
続いて、ステップＳ６において、パワー推定部１３は、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定する（パワー推定工程）。パワー推定部１３は、予め設定されたδ波のパワーモデルを利用して、サイクル毎における使用者のδ波のパワーを推定する。
【００６５】
その後、ステップＳ７において、質推定部１４は、パワー推定部１３によるサイクル毎のδ波のパワーに基づいて、使用者の睡眠の質を推定する（質推定工程）。δ波のパワーに基づく睡眠の質の推定には、周知の方法を採用することができる。
【００６６】
以上説明した睡眠状態推定装置１によれば、睡眠時の呼吸、心電（心拍）、血圧、皮膚電気活動などの生体情報に基づいて使用者の睡眠段階を推定することで、当該睡眠段階の推定結果から睡眠時間をウルトラディアンリズムのサイクル毎に分け、当該睡眠段階の推定結果からサイクル毎におけるノンレム睡眠の継続時間を推定することができる。そして、この睡眠状態推定装置１は、サイクル毎におけるノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらず同様の傾向を示すことから、この傾向を利用してサイクル毎のノンレム睡眠の継続時間からサイクル毎のδ波のパワーを推定することができる。すなわち、この睡眠状態推定装置１は、サイクル毎におけるノンレム睡眠の開始からδ波のパワーが増加してδ波のパワー最大値に至り、その後にδ波のパワーが減少してノンレム睡眠が終了するという傾向を利用して、予め設定されたδ波のパワー最大値Ｐmax１〜Ｐmax４、開始時定数τｓ１〜τｓ４、及び終了時定数τｅ１〜τｅ４を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果Ｔ１〜Ｔ４からサイクル毎における使用者のδ波のパワーを推定することができる。従って、この睡眠状態推定装置１は、睡眠時の使用者の頭部に配置する電極等を用いてδ波を直接測定することなく、使用者にとって少ない負担で睡眠時のδ波のパワーを推定することができ、δ波のパワーの推定結果に基づいて使用者の睡眠の質を推定することができる。
【００６７】
以上、本発明の好適な実施形態について説明したが、本発明は上述した実施形態に限定されるものではない。例えば、上述した実施形態に係る睡眠状態推定方法では、使用者が起床した後に、睡眠時間をウルトラディアンリズムのサイクル毎に分割して、サイクル毎のδ波のパワーを推定したが、入眠からの睡眠時間をサイクルに分割できる場合には睡眠中であってもサイクルにおけるδ波のパワーを推定してもよい。また、睡眠時間のサイクル毎の分割数は、４である場合に限られず、１以上であれば本発明を適用可能である。
【符号の説明】
【００６８】
１…睡眠状態推定装置  ２…演算部  ３…生体情報取得部  …表示部  １０…睡眠段階推定部  １１…睡眠時間分割部  １２…ノンレム睡眠時間推定部  １３…パワー推定部  １４…質推定部  １５…記憶部  Ｐ１-Ｐ４…δ波のパワー  Ｐmax１-Ｐmax４…δ波のパワー最大値  Ｔ１-Ｔ４…継続時間  τｓ１-τｓ４…開始時定数  τｅ１-τｅ４…終了時定数
【図１】
2015097713.tif 000004
【図２】
2015097713.tif 000005
【図３】
2015097713.tif 000006
【図５】
2015097713.tif 000007
【図６】
2015097713.tif 000008
【図７】
2015097713.tif 000009
【図４】
2015097713.tif 000010