JP2015097713A - 睡眠状態推定装置及び睡眠状態推定方法 - Google Patents

睡眠状態推定装置及び睡眠状態推定方法 Download PDF

Info

Publication number
JP2015097713A
JP2015097713A JP2013239900A JP2013239900A JP2015097713A JP 2015097713 A JP2015097713 A JP 2015097713A JP 2013239900 A JP2013239900 A JP 2013239900A JP 2013239900 A JP2013239900 A JP 2013239900A JP 2015097713 A JP2015097713 A JP 2015097713A
Authority
JP
Japan
Prior art keywords
sleep
power
cycle
user
wave
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2013239900A
Other languages
English (en)
Other versions
JP6121886B2 (ja
Inventor
和英 重藤
Kazuhide Shigefuji
和英 重藤
光之 中尾
Mitsuyuki Nakao
光之 中尾
彰洋 辛島
Akihiro Karashima
彰洋 辛島
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tohoku University NUC
Toyota Motor Corp
Original Assignee
Tohoku University NUC
Toyota Motor Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tohoku University NUC, Toyota Motor Corp filed Critical Tohoku University NUC
Priority to JP2013239900A priority Critical patent/JP6121886B2/ja
Publication of JP2015097713A publication Critical patent/JP2015097713A/ja
Application granted granted Critical
Publication of JP6121886B2 publication Critical patent/JP6121886B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

【課題】ウルトラディアンリズムのサイクル毎におけるδ波のパワーを推定することにより、使用者にとって少ない負担で使用者の睡眠の質を推定できる装置を提供する。
【解決手段】睡眠状態推定装置１は、使用者の睡眠時の生体情報を取得する生体情報取得部３と、生体情報に基づいて、使用者の睡眠段階を推定する睡眠段階推定部１０と、睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割部１１と、睡眠段階の推定結果に基づいて、サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定部１２と、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定するパワー推定部１３と、サイクル毎におけるパワーの推定結果に基づいて、使用者の睡眠の質を推定する質推定部１４と、を備える。
【選択図】図１

Description

本発明は、使用者の睡眠の質を推定する睡眠状態推定装置及び睡眠状態推定方法に関する。
睡眠の質は、睡眠時の脳波と相関があることが知られている。脳波は、α波、β波、γ波、δ波に分類され、δ波は眠りの深い睡眠段階において現われる傾向がある。具体的に、睡眠段階は、脳が覚醒状態にあるレム睡眠と脳が休息している状態にあるノンレム睡眠とに分けられ、ノンレム睡眠においてδ波の出現確率が高くなる。そこで、睡眠時におけるδ波等の脳波の測定結果に基づいて睡眠の質を推定することが従来から行われている。
例えば、特開２０１３−１２１４８９号公報には、使用者の前額部に配置される複数個の電極を備え、これらの電極により測定した脳波に基づいて使用者の睡眠状態の検出を行う睡眠状態検出装置が示されている。
特開２０１３−１２１４８９号公報
しかしながら、上述した従来の睡眠状態検出装置のように使用者の頭部へ電極を配置することは、使用者にとって負担であり、使用者の睡眠を妨げるおそれがある。このため、従来の睡眠状態検出装置を日常的に使用することは現実的ではないという問題があった。
そこで、本発明は、ウルトラディアンリズムのサイクル毎におけるδ波のパワーを推定することにより、使用者にとって少ない負担で使用者の睡眠の質を推定できる睡眠状態推定装置及び睡眠状態推定方法を提供することを目的とする。
本発明者等は、複数人の被験者の睡眠時の脳波（単位：μＶ）測定を十数日にわたって行う中で、睡眠時の脳活動の休息を意味するδ波が約９０分のウルトラディアンリズムのサイクルに基づいて変化することに着目した。本発明者等は、研究の結果、ウルトラディアンリズムのサイクル毎におけるノンレム睡眠時のδ波のパワー（単位：μＶ２／Ｈｚ）が人によらず同様のダイナミクスを示すことを新たに見いだし、このダイナミクスを利用してサイクル毎のノンレム睡眠の継続時間からδ波のパワーを推定することを想到した。
すなわち、本発明の一態様は、使用者の睡眠の質を推定する睡眠状態推定装置であって、使用者の睡眠時の生体情報を取得する生体情報取得部と、生体情報に基づいて、使用者の睡眠段階を推定する睡眠段階推定部と、睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割部と、睡眠段階の推定結果に基づいて、サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定部と、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定するパワー推定部と、サイクル毎におけるパワーの推定結果に基づいて、使用者の睡眠の質を推定する質推定部と、を備え、パワー推定部は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定する。なお、上記生体情報に脳波は含まれない。
この睡眠状態推定装置は、睡眠時の呼吸、心電（心拍）、血圧、皮膚電気活動などの生体情報に基づいて使用者の睡眠段階を推定することで、当該睡眠段階の推定結果から睡眠時間をウルトラディアンリズムのサイクル毎に分け、当該睡眠段階の推定結果からサイクル毎におけるノンレム睡眠の継続時間を推定することができる。そして、この睡眠状態推定装置は、サイクル毎におけるノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらず同様の傾向を示すことから、この傾向を利用してサイクル毎のノンレム睡眠の継続時間からサイクル毎のδ波のパワーを推定することができる。すなわち、この睡眠状態推定装置は、サイクル毎におけるノンレム睡眠の開始からδ波のパワーが増加してδ波のパワー最大値に至り、その後にδ波のパワーが減少してノンレム睡眠が終了するという傾向を利用して、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定することができる。従って、この睡眠状態推定装置は、睡眠時の使用者の頭部に配置する電極等を用いてδ波を直接測定することなく、使用者にとって少ない負担で睡眠時のδ波のパワーを推定することができ、δ波のパワーの推定結果に基づいて使用者の睡眠の質を推定することができる。
一実施形態において、上記パワー推定部は、使用者の入眠からのサイクルの順番に応じて、それぞれ設定されたδ波のパワー最大値を用いて、当該サイクルにおける使用者のδ波のパワーを推定してもよい。
この実施形態に係る睡眠状態推定装置は、使用者の入眠からのサイクルの順番により、当該サイクルにおけるノンレム睡眠の継続時間とδ波のパワーの関係性が変化することから、サイクルの順番に応じてそれぞれ設定されたδ波のパワー最大値を用いることで、δ波のパワー最大値として全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーの推定精度の向上を図ることができる。
一実施形態において、上記パワー推定部は、使用者の入眠からのサイクルの順番に応じて、それぞれ設定された開始時定数を用いて、当該サイクルにおける使用者のδ波のパワーを推定してもよい。
この実施形態に係る睡眠状態推定装置は、使用者の入眠からのサイクルの順番により、当該サイクルにおけるノンレム睡眠の継続時間とδ波のパワーの関係性が変化することから、サイクルの順番に応じてそれぞれ設定された開始時定数を用いることで、開始時定数として全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーの推定精度の向上を図ることができる。
一実施形態において、上記パワー推定部は、使用者の入眠からのサイクルの順番に応じて、それぞれ設定された終了時定数を用いて、当該サイクルにおける使用者のδ波のパワーを推定してもよい。
この実施形態に係る睡眠状態推定装置は、使用者の入眠からのサイクルの順番により、当該サイクルにおけるノンレム睡眠の継続時間とδ波のパワーの関係性が変化することから、サイクルの順番に応じてそれぞれ設定された終了時定数を用いることで終了時定数として全てのサイクルで一定値を用いる場合と比べて、サイクル毎におけるδ波のパワーの推定精度の向上を図ることができる。
本発明の他の態様は、使用者の睡眠の質を推定する睡眠状態推定方法であって、使用者の睡眠時の生体情報を取得する生体情報取得工程と、生体情報に基づいて、使用者の睡眠段階を推定する睡眠段階推定工程と、睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割工程と、睡眠段階の推定結果に基づいて、サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定工程と、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定するパワー推定工程と、サイクル毎におけるパワーの推定結果に基づいて、使用者の睡眠の質を推定する質推定工程と、を備え、パワー推定工程は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定する。
この睡眠状態推定方法によれば、睡眠時の呼吸、心電（心拍）、血圧、皮膚電気活動などの生体情報に基づいて使用者の睡眠段階を推定することで、当該睡眠段階の推定結果から睡眠時間をウルトラディアンリズムのサイクル毎に分け、当該睡眠段階の推定結果からサイクル毎におけるノンレム睡眠の継続時間を推定することができる。そして、この睡眠状態推定方法では、サイクル毎におけるノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらず同様のダイナミクスを示すことから、このダイナミクスを利用してサイクル毎のノンレム睡眠の継続時間からサイクル毎のδ波のパワーを推定することができる。すなわち、この睡眠状態推定方法では、サイクル毎におけるノンレム睡眠の開始からδ波のパワーが増加してδ波のパワー最大値に至り、その後にδ波のパワーが減少してノンレム睡眠が終了するという傾向を利用して、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーがパワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーがパワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、サイクル毎におけるノンレム睡眠の継続時間の推定結果からサイクル毎における使用者のδ波のパワーを推定することができる。従って、この睡眠状態推定方法によれば、睡眠時の使用者の頭部に配置する電極等を用いてδ波を直接測定することなく、使用者にとって少ない負担で睡眠時のδ波のパワーを推定することができ、δ波のパワーの推定結果に基づいて使用者の睡眠の質を推定することができる。
本発明の一態様に係る睡眠状態推定装置及び睡眠状態推定方法によれば、ウルトラディアンリズムのサイクル毎におけるδ波のパワーを推定することにより、使用者にとって少ない負担で使用者の睡眠の質を推定できる。
本実施形態に係る睡眠状態推定装置を示すブロック図である。 使用者の睡眠段階の推移の一例を示す図である。 パワー推定部の入力及び出力を説明するためのブロック図である。 （ａ）サイクル毎におけるδ波のパワーの実測データを時間軸について正規化したグラフである。（ｂ）図４（ａ）と異なる被験者について、サイクル毎におけるδ波のパワーの実測データを時間軸について正規化したグラフである。 （ａ）サイクル１におけるδ波のパワー実測値を示すグラフである。（ｂ）図５（ａ）に示すフィルタ後のサイクル１のグラフである。 （ａ）サイクル毎におけるδ波のパワーの実測値にフィルタを掛けたグラフである。（ｂ）サイクル毎のδ波のパワーモデルを示すグラフである。 本実施形態に係る睡眠状態推定方法を示すフローチャートである。
以下、本発明の好適な実施形態について、図面を参照して詳細に説明する。
図１に示す睡眠状態推定装置１は、睡眠時の使用者の生体情報を利用してウルトラディアンリズムのサイクル毎におけるδ波（脳波の一種）のパワーを推定することで、睡眠の質の推定を行うものである。使用者の生体情報には、呼吸、心電（心拍）、血圧、体温、皮膚電気活動に関する情報などが含まれる。なお、睡眠状態推定装置１は、使用者の脳波を直接測定するものではないため、使用者の生体情報に脳波は含まれない。
ウルトラディアンリズムとは、人体の生体リズムの一つであり、睡眠時におけるレム睡眠とノンレム睡眠のサイクルである。レム睡眠とは脳が覚醒状態にある浅い睡眠段階を意味し、ノンレム睡眠とは脳が休息している状態にある深い睡眠段階を意味する。ウルトラディアンリズムの周期は、個人差もあるが、一般的に約９０分である。例えば、使用者の一晩の睡眠時間が６時間（３６０分）である場合には、９０分周期のウルトラディアンリズムのサイクルが一晩で４回繰り返される。以後の説明において、特に断りがない限り、サイクルという言葉はウルトラディアンリズムのサイクルを意味する。
δ波は、約１〜３Ｈｚの周波数を持つ脳波の一種であり、単位は電圧と同じμＶである。δ波は、脳波のうち最もゆっくりした波であり、深い睡眠段階（ノンレム睡眠）において出現確率が高いことが知られている。
ウルトラディアンリズムのサイクル毎におけるδ波（脳波の一種）のパワーとは、サイクル毎におけるδ波のパワーの時間積分であるパワースペクトルに相当する。このパワーの単位は、μＶ２／Ｈｚである。
［睡眠状態推定装置の構成］
本実施形態における睡眠状態推定装置１は、装置全体を統括的に制御する演算部２を備えている。演算部２は、ＣＰＵ[Central Processing Unit]、ＲＯＭ[Read Only Memory]、ＲＡＭ[Random Access Memory]などを含む電子制御ユニットである。この演算部２は、生体情報取得部３及び表示部４と接続されている。
生体情報取得部３は、睡眠時の使用者の生体情報を取得可能なバイタルセンサから構成されている。具体的に、生体情報取得部３は、リストバンド型や貼付け型のウェアラブルセンサを含んでもよい。ウェアラブルセンサは、使用者の負担軽減の観点からワイヤレスとすることができる。ウェアラブルセンサは、使用者の呼吸、心電（心拍）、血圧、体温、皮膚電気活動等の生体情報を検出する。なお、本実施形態に係る生体情報取得部３は、使用者の脳波を直接測定する脳波計を含まない。
表示部４は、睡眠状態推定装置１における各種の検出結果及び推定結果を表示するためのディスプレイである。表示部４に表示される内容としては、使用者の入眠から起床までの睡眠時間、中途覚醒の回数及び合計時間、寝返り回数、レム睡眠の合計時間、ノンレム睡眠の合計時間、ウルトラディアンリズムのサイクル数等が挙げられる。
続いて、演算部２の構成について説明する。演算部２は、睡眠段階推定部１０、睡眠時間分割部１１、ノンレム睡眠時間推定部１２、パワー推定部１３、質推定部１４、及び記憶部１５を有している。記憶部１５は、過去の睡眠段階の比率や過去の睡眠時間などの睡眠履歴を個人別に記憶しているデータベースである。また、記憶部１５は、後述するサイクル毎のδ波のパワーモデル等の推定に用いるデータを記憶している。
睡眠段階推定部１０は、生体情報取得部３が取得した睡眠時の使用者の生体情報に基づいて、現在の使用者の睡眠段階の推定を行う。睡眠段階は、レム睡眠とノンレム睡眠に大きく分けられる。睡眠段階推定部１０は、ノンレム睡眠を更に複数の段階に分けて推定を行ってもよい。
睡眠段階推定部１０は、例えば、使用者の生体情報から特徴量を算出し、当該特徴量に基づいて睡眠段階を推定する。特徴量とは、例えば生体情報としての心拍の検出データから抽出される単位時間当たりの心拍数、心拍の周期のゆらぎ、心拍の強さなどである。睡眠段階推定部１０は、使用者の生体情報から様々な観点の特徴量を算出し、これらの特徴量に基づいて使用者の睡眠段階を多元的に推定する。
睡眠時間分割部１１は、睡眠段階推定部１０による睡眠段階の推定結果に基づいて、使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分割する。ここで、図２は、使用者の睡眠段階の推定結果の一例を示す図である。
図２に示されるように、使用者は、入眠によってレム睡眠に入った後、周期的にレム睡眠とノンレム睡眠とを繰り返す。睡眠時間分割部１１は、例えば、睡眠段階の推定結果からレム睡眠とノンレム睡眠との繰り返しを検出することで、睡眠時間をウルトラディアンリズムのサイクル毎に分ける。使用者の入眠から起床までの睡眠時間は、３〜４のサイクルに分けられることが通常である。図２では、使用者の入眠から起床までの睡眠時間が４つのサイクルに分けられている。
以下、図２に示されるように、使用者の入眠から１番目のサイクルをサイクル１、２番目のサイクルをサイクル２、３番目のサイクルをサイクル３、４番目のサイクルをサイクル４として説明に用いる。
なお、睡眠時間分割部１１は、ウルトラディアンリズムのサイクルを９０分と仮定して睡眠時間を９０分毎に分割した後、レム睡眠とノンレム睡眠との繰り返しを検出して睡眠時間の分割位置を補正する態様であってもよい。
図１に示すノンレム睡眠時間推定部１２は、睡眠段階推定部１０による睡眠段階の推定結果に基づいて、睡眠時間分割部１１によって分けられたサイクル毎におけるノンレム睡眠の継続時間を推定する。具体的に、図２において、ノンレム睡眠時間推定部１２は、睡眠段階の推定結果からサイクル１におけるノンレム睡眠の継続時間Ｔ１を推定する。同様に、ノンレム睡眠時間推定部１２は、サイクル２におけるノンレム睡眠の継続時間Ｔ２を推定し、サイクル３におけるノンレム睡眠の継続時間Ｔ３を推定する。また、ノンレム睡眠時間推定部１２は、サイクル４におけるノンレム睡眠の継続時間Ｔ４を推定する。
図１及び図３に示すパワー推定部１３は、サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、サイクル毎における使用者のδ波のパワーを推定する。パワー推定部１３は、予め設定されたδ波のパワーモデルを利用して、サイクル毎における使用者のδ波のパワーを推定する。
以下、δ波のパワーモデルについて説明する。δ波のパワーモデルとは、サイクル毎におけるノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらず共通する傾向を示すことに着目して作成されたものである。このようなδ波のパワーモデルを用いることで、あるサイクルにおけるノンレム睡眠の継続時間から当該サイクルにおけるδ波のパワーを推定することができる。δ波のパワーモデルのデータは、例えば記憶部１５に記憶されている。
また、δ波のパワーモデルは、使用者の入眠からのサイクルの順番に応じて複数のモデルが予め設定されている。これは、サイクルの順番が同じ場合にノンレム睡眠の継続時間とδ波のパワーとの関係性が人によらずに同様のダイナミクスを示すためである。なお、図３には、サイクル１〜４の各々のδ波のパワーモデルが示されているが、サイクル５以降のパワーモデルについて予め用意していてもよい。
サイクル１〜４のδ波のパワーモデルは、下記の式（１）〜（４）に示すモデルパラメータとして表すことができる。式（１）は、サイクル１に対応するモデルパラメータであり、式（２）はサイクル１に対応するモデルパラメータである。式（３）は、サイクル３に対応するモデルパラメータであり、式（４）はサイクル４に対応するモデルパラメータである。
上記の式（１）において、Ｐ１は、サイクル１におけるδ波のパワーである。Ｐ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４…δ波のパワー最大値 Ｔ１-Ｔ４…継続時間 τｓ１-τｓ４…開始時定数 τｅ１-τｅ４…終了時定数

Claims (5)

  1. 使用者の睡眠の質を推定する睡眠状態推定装置であって、
    前記使用者の睡眠時の生体情報を取得する生体情報取得部と、
    前記生体情報に基づいて、前記使用者の睡眠段階を推定する睡眠段階推定部と、
    前記睡眠段階の推定結果に基づいて、前記使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割部と、
    前記睡眠段階の推定結果に基づいて、前記サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定部と、
    前記サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、前記サイクル毎における前記使用者のδ波のパワーを推定するパワー推定部と、
    前記サイクル毎における前記パワーの推定結果に基づいて、前記使用者の睡眠の質を推定する質推定部と、
    を備え、
    前記パワー推定部は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーが前記パワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーが前記パワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、前記サイクル毎におけるノンレム睡眠の継続時間の推定結果から前記サイクル毎における前記使用者のδ波のパワーを推定する、睡眠状態推定装置。
  2. 前記パワー推定部は、前記使用者の入眠からの前記サイクルの順番に応じて、それぞれ設定された前記δ波のパワー最大値を用いて、当該サイクルにおける前記使用者のδ波のパワーを推定する、請求項１に記載の睡眠状態推定装置。
  3. 前記パワー推定部は、前記使用者の入眠からの前記サイクルの順番に応じて、それぞれ設定された前記開始時定数を用いて、当該サイクルにおける前記使用者のδ波のパワーを推定する、請求項１又は２に記載の睡眠状態推定装置。
  4. 前記パワー推定部は、前記使用者の入眠からの前記サイクルの順番に応じて、それぞれ設定された前記終了時定数を用いて、当該サイクルにおける前記使用者のδ波のパワーを推定する、請求項１〜３のうち何れか一項に記載の睡眠状態推定装置。
  5. 使用者の睡眠の質を推定する睡眠状態推定方法であって、
    前記使用者の睡眠時の生体情報を取得する生体情報取得工程と、
    前記生体情報に基づいて、前記使用者の睡眠段階を推定する睡眠段階推定工程と、
    前記睡眠段階の推定結果に基づいて、前記使用者の入眠からの睡眠時間をウルトラディアンリズムのサイクル毎に分ける睡眠時間分割工程と、
    前記睡眠段階の推定結果に基づいて、前記サイクル毎におけるノンレム睡眠の継続時間を推定するノンレム睡眠時間推定工程と、
    前記サイクル毎におけるノンレム睡眠の継続時間の推定結果に基づいて、前記サイクル毎における前記使用者のδ波のパワーを推定するパワー推定工程と、
    前記サイクル毎における前記パワーの推定結果に基づいて、前記使用者の睡眠の質を推定する質推定工程と、
    を備え、
    前記パワー推定工程は、予め設定されたδ波のパワー最大値と、ノンレム睡眠の開始からδ波のパワーが前記パワー最大値に至るまでの時間に関して予め設定された開始時定数と、δ波のパワーが前記パワー最大値から減少してノンレム睡眠の終了に至るまでの時間に関して予め設定された終了時定数と、を用いて、前記サイクル毎におけるノンレム睡眠の継続時間の推定結果から前記サイクル毎における前記使用者のδ波のパワーを推定する、睡眠状態推定方法。
JP2013239900A 2013-11-20 2013-11-20 睡眠状態推定装置及び睡眠状態推定方法 Active JP6121886B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013239900A JP6121886B2 (ja) 2013-11-20 2013-11-20 睡眠状態推定装置及び睡眠状態推定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013239900A JP6121886B2 (ja) 2013-11-20 2013-11-20 睡眠状態推定装置及び睡眠状態推定方法

Publications (2)

Publication Number Publication Date
JP2015097713A true JP2015097713A (ja) 2015-05-28
JP6121886B2 JP6121886B2 (ja) 2017-04-26

Family

ID=53374767

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013239900A Active JP6121886B2 (ja) 2013-11-20 2013-11-20 睡眠状態推定装置及び睡眠状態推定方法

Country Status (1)

Country Link
JP (1) JP6121886B2 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111657860A (zh) * 2020-05-14 2020-09-15 清华大学 一种基于对抗学习的睡眠阶段识别方法及系统
CN113576410A (zh) * 2021-07-20 2021-11-02 电子科技大学 一种入睡过程动态连续分析方法
US11166664B2 (en) 2016-03-28 2021-11-09 Fujitsu Limited Information processing device, information processing method, and recording medium
WO2022025678A1 (en) * 2020-07-31 2022-02-03 Samsung Electronics Co., Ltd. Electronic device for evaluating sleep quality and method for operation in the electronic device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120253220A1 (en) * 2011-03-31 2012-10-04 Toyota Motor Engineering And Manufacturing North America Methods and Systems for Determining Optimum Wake Time
WO2012176790A1 (ja) * 2011-06-21 2012-12-27 スリープウェル株式会社 精神疾患解析装置、精神疾患解析方法、およびプログラム
JP2013220315A (ja) * 2012-04-19 2013-10-28 Toyota Motor Corp 睡眠品質推定装置、睡眠品質推定方法及び睡眠品質推定用プログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120253220A1 (en) * 2011-03-31 2012-10-04 Toyota Motor Engineering And Manufacturing North America Methods and Systems for Determining Optimum Wake Time
WO2012176790A1 (ja) * 2011-06-21 2012-12-27 スリープウェル株式会社 精神疾患解析装置、精神疾患解析方法、およびプログラム
JP2013220315A (ja) * 2012-04-19 2013-10-28 Toyota Motor Corp 睡眠品質推定装置、睡眠品質推定方法及び睡眠品質推定用プログラム

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11166664B2 (en) 2016-03-28 2021-11-09 Fujitsu Limited Information processing device, information processing method, and recording medium
CN111657860A (zh) * 2020-05-14 2020-09-15 清华大学 一种基于对抗学习的睡眠阶段识别方法及系统
CN111657860B (zh) * 2020-05-14 2021-07-20 清华大学 一种基于对抗学习的睡眠阶段识别方法及系统
WO2022025678A1 (en) * 2020-07-31 2022-02-03 Samsung Electronics Co., Ltd. Electronic device for evaluating sleep quality and method for operation in the electronic device
CN113576410A (zh) * 2021-07-20 2021-11-02 电子科技大学 一种入睡过程动态连续分析方法
CN113576410B (zh) * 2021-07-20 2022-09-02 电子科技大学 一种入睡过程动态连续分析方法

Also Published As

Publication number Publication date
JP6121886B2 (ja) 2017-04-26

Similar Documents

Publication Publication Date Title
EP2747638B1 (en) Mental balance or imbalance estimation system and method
EP3324842B1 (en) Device and method for assessing the level of consciousness, pain and nociception during wakefulness, sedation and general anaesthesia
JP5904275B2 (ja) 保有活動量算出装置、保有活動量算出方法及び保有活動量算出システム
JP4469747B2 (ja) 目覚まし装置、目覚まし方法および目覚ましプログラム
WO2014083657A1 (ja) 睡眠制御装置及び睡眠制御方法
JP2018524080A (ja) 被検者の生理学的状態を監視する装置及び方法
CN110662487B (zh) 用于确定用户的情绪状态的设备、系统和方法
TW201538127A (zh) 睡眠檢測系統和方法
JP5696501B2 (ja) 覚醒時データ生成装置、覚醒時データ生成方法、覚醒時データ生成プログラム及び覚醒度判定装置
JP5632986B2 (ja) 睡眠段階判定装置及び睡眠段階判定方法
CN104615851B (zh) 一种睡眠监控方法及终端
JP6121886B2 (ja) 睡眠状態推定装置及び睡眠状態推定方法
US20080071152A1 (en) Autonomic-nervous state judging device, autonomic-nervous state judging method, and computer program product
JP5874489B2 (ja) 睡眠状態判定装置及び睡眠状態判定方法
JP2004254827A (ja) 睡眠状態判定装置
JP6040874B2 (ja) 睡眠段階推定装置
US20040230398A1 (en) Sleep analyzer and program product for giving sleep analysis function to computer
JP6304050B2 (ja) 生体状態推定装置
EP2191772B1 (en) Measurement of responsiveness of a subject
CN115171850A (zh) 一种睡眠方案生成方法、装置、终端设备及存储介质
Okada et al. Wearable ECG recorder with acceleration sensors for monitoring daily stress: Office work simulation study
Klein et al. Sleep stages classification using vital signals recordings
JP3293213B2 (ja) 生体リズム曲線測定装置
US20080214908A1 (en) Method and Apparatus For Monitoring the Sedation Level of a Sedated Patient
JPH053877A (ja) 生体リズム曲線測定装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160112

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20161011

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20161012

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20161109

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170307

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170330

R151 Written notification of patent or utility model registration

Ref document number: 6121886

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250