(57)【要約】【課題】睡眠中の不整脈や無呼吸などの異常を検出すると共に脈波データへの体動の影響を低減し、もって、睡眠状態の判定精度を向上させること。【解決手段】被験者の脈波の一周期の時間間隔データである脈拍間隔データと被験者の体動を示す体動データに基づいて被験者の睡眠状態を判定する睡眠状態判定装置１０であって、体動判定部１９１は、体動データの変動が第１の所定の閾値より大きい場合に体動と判定し、脈拍間隔データ処理部２２は、体動判定部１９１が体動である判定した場合に体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理し、睡眠状態判定部２４は、脈拍間隔データ処理部２２によってデータ処理された一連の脈拍間隔データから取得した自律神経指標に基づいて睡眠状態を判定する。【選択図】        図１
2005279113.tif 000002
【特許請求の範囲】
【請求項１】
被験者の脈波の一周期の時間間隔を示す脈拍間隔データと該被験者の体動を示す体動データとに基づいて該被験者の睡眠状態を判定する睡眠状態判定装置であって、
前記体動データの変動量が第１の所定の閾値より大きい場合に体動ありと判定する体動判定手段と、
前記体動判定手段によって体動ありと判定された場合に該体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理する脈拍間隔データ処理手段と、
前記脈拍間隔データ処理手段によってデータ処理された一連の脈拍間隔データから取得した自律神経指標に基づいて前記睡眠状態を判定する睡眠状態判定手段と、
を備えたことを特徴とする睡眠状態判定装置。
【請求項２】
前記体動判定手段によって判定された体動の発生頻度が第２の所定の閾値以上である場合に覚醒中の体動と判定し、該体動の発生頻度が該第２の閾値未満である場合は睡眠中の体動と判定する覚醒判定手段とをさらに備え、
前記脈拍間隔データ処理手段は、前記覚醒判定手段によって睡眠中の体動と判定された場合に該睡眠中の体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理をすることを特徴とする請求項１に記載の睡眠状態判定装置。
【請求項３】
前記覚醒判定手段は、前記体動判定手段によって判定された体動の発生頻度が第２の所定の閾値以上であり、かつ、前記脈拍間隔データが過去の睡眠中の脈拍間隔データの平均値より短い場合は、覚醒中の体動と判定することを特徴とする請求項２に記載の睡眠状態判定装置。
【請求項４】
前記覚醒判定手段によって前記体動と判定された場合に該体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外する異常データ除外手段をさらに備え、
前記脈拍間隔データ処理手段は、前記異常データ除外手段によって前記体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理をすることを特徴とする請求項２または３に記載の睡眠状態判定装置。
【請求項５】
前記異常データ除外手段によって該体動と判定された体動データと並行して計測された脈拍間隔データを除外された一連の脈拍間隔データを高次の多項式で補間する脈拍間隔データ補間手段をさらに備え、
前記脈拍間隔データ処理手段は、前記脈拍間隔データ補間手段によって補間された一連の脈拍間隔データに基づいて時間的に等間隔の一連の脈拍間隔データにデータ処理をすることを特徴とする請求項４に記載の睡眠状態判定装置。
【請求項６】
前記脈拍間隔データ補間手段によって高次の多項式で補間された一連の脈拍間隔データを周波数スペクトル分布に変換し、周波数スペクトル分布に変換された一連の脈拍間隔データの複数のパワースペクトルの値から自律神経指標を取得する自律神経指標取得手段とを、さらに備え、
前記睡眠状態判定手段は、前記自律神経指標取得手段によって取得された自律神経指標に基づいて前記睡眠状態を判定することを特徴とする請求項５に記載の睡眠状態判定装置。
【請求項７】
前記睡眠状態判定手段は、前記複数のパワースペクトルの値のうち、低周波数の領域にあるパワースペクトルの値をＬＦ、高周波数の領域にあるパワースペクトルの値をＨＦとしたときに、ＬＦ／ＨＦの値が第１の判定閾値よりも小さく、かつ、ＨＦの値が第２の判定閾値２よりも大きいときは深睡眠と判定し、また、ＬＦ／ＨＦの値が第３の判定閾値よりも大きく、かつ、ＨＦの値が第４の判定閾値より小さく、かつ、ＬＦ、ＨＦの標準偏差の合計が第５の判定閾値より大きいときはレム睡眠と判定し、深睡眠およびレム睡眠以外のときは浅睡眠と前記睡眠状態を判定することを特徴とする請求項６に記載の睡眠状態判定装置。
【請求項８】
時刻を計測する時刻計測手段をさらに備え、
前記睡眠状態判定手段は、前記時刻計測手段によって計測された時刻に応じた前記被験者のサーカディアンリズムの影響を考慮に入れて前記第１の判定閾値から前記第５の判定閾値までの判定閾値を変更することを特徴とする請求項７に記載の睡眠状態判定装置。
【請求項９】
前記脈波から脈波データをサンプリングし、サンプリングした一連の脈波データを時間微分して該一連の脈波データから該直流変動成分を除去する直流変動成分除去手段と、
前記直流変動成分除去手段によって直流変動成分を除去された一連の脈波データの所定の時間間隔内の脈波データの最大値と最小値を取得し、該最大値と該最小値との間の所定の値を第３の所定の閾値とする閾値取得手段と、
前記直流変動成分除去手段によって直流変動成分を除去された一連の脈波データから前記閾値取得手段によって取得された第３の所定の閾値に一致する該一連の脈波データの値が現れた時刻を算出し、算出された時刻の間隔から脈拍間隔データを取得する脈拍間隔データ取得手段と、をさらに備え、
前記脈拍間隔データ処理手段は、前記脈拍間隔データ取得手段によって取得された脈拍間隔データから一連の脈拍間隔データを生成し、該一連の脈拍間隔データをデータ処理することを特徴とする請求項１に記載の睡眠状態判定装置。
【請求項１０】
前記脈拍間隔データ取得手段は、前記第３の所定の閾値を超える直前の前記一連の脈波データの値と時刻をｙ１とｔ１とし、該第３の所定の閾値を超えた直後の該一連の脈波データの値と時刻をｙ２とｔ２として、該第３の所定の閾値に一致する該一連の脈波データの値と時刻をｙ０とｔとしたときに、ｔ＝ｔ１＋（ｔ２−ｔ１）×（ｙ０−ｙ１）／（ｙ２−ｙ１）によって該第３の所定の閾値に一致する前記一連の脈波データの値が現れた時刻ｔを算出することを特徴とする請求項９に記載の睡眠状態判定装置。
【請求項１１】
前記体動判定手段によって体動と判定された脈拍間隔データの波形と正常な脈拍間隔データの波形との一致度を算出して該一致度が第４の所定の閾値未満である場合は、該体動と判定された脈拍間隔データを波形異常として検出する波形異常検出手段をさらに備え、
前記脈拍間隔データ処理手段は、前記波形異常検出手段によって波形異常として検出された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理をすることを特徴とする請求項１に記載の睡眠状態判定装置。
【請求項１２】
前記脈拍間隔データ処理手段は、前記覚醒判定手段が前記睡眠中の体動と判定した場合に、前記波形異常検出手段によって波形異常として検出された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理をすることを特徴とする請求項１に記載の睡眠状態判定装置。
【請求項１３】
前記波形異常検出手段は、前記覚醒判定手段によって睡眠中の体動と判定された脈拍間隔データの波形と正常な脈拍間隔データの波形との相関係数を算出して該相関係数が第４の所定の閾値未満の場合は、該脈拍間隔データを波形異常として検出することを特徴とする請求項１１または１２に記載の睡眠状態判定装置。
【請求項１４】
前記波形異常検出手段は、前記覚醒判定手段によって睡眠中の体動と判定された脈拍間隔データの波形と異常な脈拍間隔データの波形とをマッチングして、波形の識別と発生頻度を検出することを特徴とする請求項１１乃至１３のいずれか一つに記載の睡眠状態判定装置。
【請求項１５】
３軸方向の加速度データを前記体動データとして計測する体動計測手段と、
前記体動計測手段によって計測された３軸方向の加速度データを時間微分して３軸方向の加速度データの微係数を求め、該３軸方向の微係数のそれぞれの二乗和の平方根である前記体動データの変動量を算出する体動サンプリングデータ処理手段と、
をさらに備え、
前記体動判定手段は、前記体動サンプリングデータ処理手段によって算出された体動データの変動量が第１の所定の閾値より大きい場合に、前記体動と判定することを特徴とする請求項１乃至４または１１のいずれか一つに記載の睡眠状態判定装置。
【請求項１６】
被験者の脈波の一周期の時間間隔データである脈拍間隔データと該被験者の体動を示す体動データに基づいて該被験者の睡眠状態を判定する睡眠状態判定方法であって、
前記体動データの変動量が第１の所定の閾値より大きい場合に前記体動と判定する体動判定工程と、
前記体動判定工程によって前記体動と判定された場合に該体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理する脈拍間隔データ処理工程と、
前記脈拍間隔データ処理工程によってデータ処理された一連の脈拍間隔データから取得した自律神経指標に基づいて前記睡眠状態を判定する睡眠状態判定工程と、
を含んだことを特徴とする睡眠状態判定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、睡眠状態判定装置および睡眠状態判定方法に関し、特に、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる睡眠状態判定装置および睡眠状態判定方法に関するものである。
【背景技術】
【０００２】
従来、被験者の脈波の一周期のデータである脈拍間隔データと被験者の体動を示す体動データに基づいて被験者の睡眠状態を判定する睡眠状態判定装置が知られていた。かかる睡眠状態判定装置は、睡眠ポリグラフと呼ばれる脳波、眼球運動、筋電、心電などの生体信号のパターンから睡眠状態を自動的に判定する大掛かりな装置に比較して、日常生活において手軽に睡眠状態を判定することができる睡眠状態判定装置として注目されている。
【０００３】
かかる睡眠状態判定装置は、睡眠中の自律神経の活動である心拍の拍動間隔を脈波の脈拍間隔として捉え、脈拍間隔の変動から取得した自律神経指標に基づいて睡眠状態の判定を行っている。例えば、手の血管の血流変化である脈波は、心拍に同期して変動するために、脈波の脈拍間隔から心拍の拍動間隔を取得できる。例えば、特許文献１および特許文献２の従来技術では、脈波データの周波数スペクトル成分から求められた自律神経指標に基づく睡眠状態の判定が行われている。即ち、脈波データから一連の脈拍間隔データを求め、一連の脈拍間隔データを周波数スペクトル分布に変換し、周波数スペクトル分布に変換した一連の脈拍間隔データから求めた低周波数領域（０．０５〜０．１５Ｈｚ付近）と高周波数領域（０．１５〜０．４Ｈｚ付近）におけるパワースペクトルの値から自律神経指標を取得し、自律神経指標から睡眠状態の判定を行っている。また、特許文献３の従来技術では、脈波データと共に体動データとを計測して、覚醒、レム睡眠、ノンレム睡眠、中途覚醒などの睡眠状態の判定を行っている。
【０００４】
【特許文献１】特開２００２−２９１７１０公報
【特許文献２】特開平０７−１４３９７２号公報
【特許文献３】特開２００２−３４９５５公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記の特許文献１〜３の従来技術では、覚醒、レム睡眠、ノンレム睡眠、中途覚醒などの睡眠状態の判定を行うことはできるが、脈波データは、手の血管の血流変化である脈波を計測するものであるので、手足の動きなど体動の影響を受け易く睡眠状態の判定精度が低いという問題点があった。また、脈波データは、不整脈や無呼吸状態などの脈波異常の影響を受け易く睡眠状態の判定精度が低いという問題点があった。
【課題を解決するための手段】
【０００６】
上述した課題を解決し、目的を達成するため、本発明に係る睡眠状態判定装置は、被験者の脈波の一周期の時間間隔データである脈拍間隔データと該被験者の体動を示す体動データとに基づいて該被験者の睡眠状態を判定する睡眠状態判定装置であって、前記体動データの変動が第１の所定の閾値より大きい場合に前記体動と判定する体動判定手段と、  前記体動判定手段によって前記体動と判定された場合に該体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理する脈拍間隔データ処理手段と、前記脈拍間隔データ処理手段によってデータ処理された一連の脈拍間隔データから取得した自律神経指標に基づいて前記睡眠状態を判定する睡眠状態判定手段と、を備えたことを特徴とする。
【０００７】
また、本発明に係る睡眠状態判定方法は、被験者の脈波の一周期のデータである脈拍間隔データと該被験者の体動を示す体動データとに基づいて該被験者の睡眠状態を判定する睡眠状態判定方法であって、前記体動データの変動が第１の所定の閾値より大きい場合に前記体動と判定する体動判定工程と、前記体動判定工程によって前記体動と判定された場合に該体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理する脈拍間隔データ処理工程と、前記脈拍間隔データ処理工程によってデータ処理された一連脈拍間隔データから取得した自律神経指標に基づいて前記睡眠状態を判定する睡眠状態判定工程と、を含んだことを特徴とする。
【０００８】
この発明によれば、体動データの変動が第１の所定の閾値より大きい場合に体動と判定し、体動と判定された場合に体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理し、データ処理された一連脈拍間隔データから取得した自律神経指標に基づいて睡眠状態を判定することとしたので、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【発明の効果】
【０００９】
本発明によれば、体動データの変動が第１の所定の閾値より大きい場合に体動と判定し、体動と判定された場合に体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理し、データ処理された一連脈拍間隔データから取得した自律神経指標に基づいて睡眠状態を判定するよう構成したので、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができるという効果を奏する。
【発明を実施するための最良の形態】
【００１０】
以下に添付図面を参照して、この発明に係る睡眠状態判定装置の好適な実施の形態を詳細に説明する。実施の形態１では、加速度センサを用いて脈波センサへの体動の影響を低減する場合に適用し、実施の形態２では、睡眠中の不整脈や無呼吸などの異常を検出する場合に適用する。そして、実施の形態３では、体動センサとしてマットセンサを用いる場合に適用する。最後に、実施の形態４では、体動センサおよび脈波センサとしてマットセンサを用いる場合に適用する。
【００１１】
（実施の形態１）
実施の形態１では、本発明に係る睡眠状態判定装置を加速度センサを用いて脈波センサへの体動の影響を低減する場合について説明する。まず、本発明に係る睡眠状態判定装置の構成について説明する。図１は、本発明に係る睡眠状態判定装置の構成を示す機能ブロック図である。同図に示すように、睡眠状態判定装置１０は、入力部１１と、表示部１２と、記憶部１３と、体動計測部１４と、脈波計測部１５と、データ通信部１６と、時刻取得部１７と、電力供給部１８と、覚醒／睡眠判定部１９と、体動サンプリングデータ処理部２０と、脈波データサンプリング処理部２１と、脈波間隔データ処理部２２と、自律神経指標データ処理部２３と、睡眠状態判定部２４と、制御部２５とを有する。
【００１２】
ここで、図１に示す睡眠状態判定装置１０の装着の例について説明する。図２は、図１に示す睡眠状態判定装置１０の装着の一例を示す図である。また、図３は、図１に示す睡眠状態判定装置１０の装着の別の例を示す図である。図２では、脈波センサ１５１は指に、睡眠状態判定装置本体１０は腕時計のように手首に装着されているが、脈波センサ１４１は、手のひらにバンソウ膏で装着することもできる。また、図３に示すように、赤外や赤色ＬＥＤを用いた脈波センサ１５１と睡眠状態判定装置本体１０とを一体化して手首の動脈上に装着することもできる。
【００１３】
図１の説明に戻ると、入力部１１は、ユーザが電源をＯＮ／ＯＦＦする、または表示を切り替える要求や指示を行うスイッチである。また、表示部１２は、睡眠状態判定結果を表示する表示装置であり、具体的には、ＬＣＤなどである。また、記憶部１３は、脈波データ、体動データなどの計測データ、脈拍間隔データなど処理後のデータ、睡眠状態を判定する閾値などのデータを記憶する記憶部であり、具体的には、フラッシュメモリなどである。
【００１４】
体動計測部１４は、被験者の体動を示す体動データとして加速度データを計測し、データ変換をする計測部であり、加速度センサ１４１を有する。加速度センサ１４１は、３軸方向の−２ｇ〜２ｇの加速度を計測する加速度計であり、睡眠状態判定装置本体１０内に搭載されている。また、体動計測部１４は、加速度センサ１４１のアナログデータのゲイン、オフセットを調整回路で調整した後、１０ビットＡ／Ｄ変換器でデジタル量に変換して制御部２５に入力する。
【００１５】
脈波計測部１５は、被験者の脈波データを計測し、データ変換をする計測部であり、脈波センサ１５１を有する。脈波センサ１５１は、青色ＬＥＤとフォトダイオードからなり、指の皮膚表面に光を照射し、毛細血管内の血流変化により変化する反射光の変動をフォトダイオードで捉えることで脈波を計測する。また、脈波計測部１５は、脈波センサ１５１のフォトダイオードからの出力電流を電流電圧変換器で電圧に変換し、増幅器で電圧を増幅して、ハイパスフィルタ（カットオフ周波数：０．１Ｈｚ）とローパスフィルタ（カットオフ周波数：５０Ｈｚ）を施した後、１０ビットＡ／Ｄ変換器でデジタル量に変換して制御部２５に入力する。
【００１６】
データ通信部１６は、無線ＬＡＮなどでパソコンやＰＤＡ端末とデータ通信を行う通信部であり、具体的には、Ｂｌｕｅｔｏｏｔｈ（ＴＭ）などである。また、時刻計測部１７は、時刻を計測する装置であり、具体的には、リアルクロックＩＣなどである。また、電力供給部１８は、睡眠状態判定装置１０の電力を供給する電源であり、具体的には、バッテリである。
【００１７】
覚醒／睡眠判定部１９は、被験者の体動データに基づいて被験者が覚醒しているか否かを判定する判定部であり、体動判定部１９１と覚醒判定部１９２を有する。体動判定部１９１は、体動データの変動量が第１の所定の閾値より大きい場合に体動と判定する判定部である。例えば、第１の所定の閾値として体動計に使用されている微小な体動の最小値である０．０１Ｇを用いる。また、覚醒判定部１９２は、体動判定部１９１によって判定された体動の発生頻度が第２の所定の閾値以上である場合に覚醒中の体動と判定し、体動の発生頻度が第２の所定の閾値未満である場合は睡眠中の体動と判定する判定部である。また、覚醒判定部１９２は、体動判定部１９１によって判定された体動の発生頻度が第２の所定の閾値以上であり、かつ、脈拍間隔データが過去の睡眠中の脈拍間隔データの平均値より短い場合は、覚醒中の体動と判定する。例えば、第２の所定の閾値として過去の覚醒時における体動頻度から２０回／分を用いる。
【００１８】
体動サンプリングデータ処理部２０は、体動計測部１４から取得した３軸方向の加速度データを時間微分して３軸方向の加速度の微係数を求め、３軸方向の加速度のそれぞれの微係数の二乗和の平方根である体動データの変動量および脈拍間隔内の体動データの変動量の平均である体動量を求めるデータ処理部である。体動サンプリングデータ処理部２０は、体動データの変動量および体動量を体動判定のためのデータとして体動判定部１９１に提供する。
【００１９】
脈波サンプリングデータ処理部２１は、被験者の脈波から脈波データをサンプリングして脈拍間隔データを取得するデータ処理部であり、脈波間隔データ取得部２１１を有する。脈拍間隔データ取得部２１１は、被験者の脈波の一周期の時間間隔データである脈拍間隔データを取得するデータ処理部である。
【００２０】
具体的には、脈波から脈波データをサンプリングし、サンプリングした一連の脈波データを時間微分して、一連の脈波データから直流変動成分を除去する。そして、直流変動成分を除去された一連の脈波データの処理ポイントを中心とした前後役１秒の脈波データの最大値と最小値を取得し、最大値と最小値との間の所定の値を第３の所定の閾値とする。例えば、第３の所定の閾値として最大値、最小値の差を振幅として、最小値から振幅の９割の値を用いる。さらに、直流変動成分を除去された一連の脈波データから第３の所定の閾値に一致する一連の脈波データの値が現れた時刻を算出し、算出された時刻の間隔から脈拍間隔データを取得する。ここで、さらに具体的に、図１に示す脈拍間隔データ取得部２１１が脈拍間隔データを算出する算出方法について説明する。図４は、図１に示す脈拍間隔データ取得部２１１が脈拍間隔データを算出する算出方法を示す図である。
【００２１】
同図に示すように、第３の所定の閾値を超える直前の一連の脈波データの値と時刻をｙ１とｔ１とし、第３の所定の閾値を超えた直後の一連の脈波データの値と時刻をｙ２とｔ２として、第３の所定の閾値に一致する一連の脈波データの値と時刻をｙ０とｔとした場合に、次式で第３の所定の閾値に一致する一連の脈波データの値が現れた時刻ｔを算出する。
ｔ＝ｔ１＋（ｔ２−ｔ１）×（ｙ０−ｙ１）／（ｙ２−ｙ１）      （式１）
したがって、このようにして算出された第３の所定の閾値に一致する一連の脈波データの値が現れた前回の時刻をｔ０とすると、脈拍間隔データは、時間間隔（ｔ−ｔ０）で求められる。
【００２２】
図１の説明に戻ると、脈拍間隔データ処理部２２は、脈波サンプリングデータ処理部２１によって取得された脈拍間隔データから一連の脈拍間隔データ、例えば、１分間のデータセットを生成し、一連の脈拍間隔データをデータ処理する処理部であり、波形異常検出部２２１と、異常データ除外部２２２と、脈拍間隔データ補間部２２３とを有する。波形異常検出部２２１は、体動判定部１９１によって体動と判定された場合に、体動と判定された体動データの体動量が第５の所定の閾値以上の場合に、体動量が大きいと判定し、体動と判定された体動データと並行して計測された脈拍間隔データを波形異常として検出する検出部である。例えば、第５の所定の閾値として１Ｇを用いる。
【００２３】
異常データ除外部２２２は、波形異常検出部２２１によって波形異常として検出された脈拍間隔データを一連の脈拍間隔データから除外する処理部である。また、脈拍間隔データ補間部２２３は、異常データ除外部２２２によって波形異常として検出された脈拍間隔データを除外された一連の脈拍間隔データを高次の多項式で補間する処理部である。ここで、具体的に、図１に示す異常データ除外部２２２が体動と判定された脈拍間隔データを除外した一連の脈拍間隔データの一例について説明する。また、図１に示す脈拍間隔データ補間部２２３が一連の脈拍間隔データを補間する一例について説明する。図５は、図１に示す異常データ除外部２２２が体動と判定された脈拍間隔データを除外した一連の脈拍間隔データの一例を示す図である。また、図６は、図１に示す脈拍間隔データ補間部２２３が一連の脈拍間隔データを補間する一例を示す図である。
【００２４】
図５示すように、脈拍間隔データ処理部２１がデータ処理する脈拍間隔データは、元々不等間隔データである上に、体動や波形異常によって異常データが除外されている。そこで、図６に示すように、不等間隔の脈拍間隔データを補間、再サンプリングし、等間隔の脈拍間隔データを生成する。例えば、３次の多項式補間法によって補間する点の前後それぞれ３点のサンプリング点を用いて等間隔の脈拍間隔データを生成する。このようにして脈拍間隔データ処理部２２が不等間隔の一連の脈拍間隔データから等間隔の脈拍間隔データを生成した後に、自律神経指標データ処理部２３が等間隔の脈拍間隔データを周波数スペクトル分布に変換する。
【００２５】
図１の説明に戻ると、自律神経指標データ処理部２３は、睡眠状態を判定する低周波数領域（０．０５〜０，１５Ｈｚ付近）の指標ＬＦと高周波数領域（０．１５〜０．４Ｈｚ付近）の指標ＨＦという二つの自律神経指標を取得するデータ処理部であり、周波数スペクトル変換部２３１と自律神経指標取得部２３２とを有する。周波数スペクトル変換部２３１は、脈拍間隔データ処理部２２によってデータ処理された一連の脈拍間隔データをＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）法などの解析手法によって周波数スペクトル分布に変換するデータ処理部である。また、自律神経指標取得部２３２は、周波数スペクトル変換部２３１によって周波数スペクトル分布に変換された一連の脈拍間隔データの複数のパワースペクトルの値から自律神経指標ＬＦ，ＨＦを取得するデータ処理部である。具体的には、複数のパワースペクトルのピーク値とピーク値を中心として前後等間隔の１点との３点の合計値の算術平均をとってＬＦ、ＨＦとする。
【００２６】
ここで、図１に示す自律神経指標データ処理部２３がパワースペクトルの値から自律神経指標ＬＦ，ＨＦを取得する一例について説明する。図７は、図１に示す自律神経指標データ処理部２３がパワースペクトルの値から自律神経指標ＬＦ，ＨＦを取得する一例を示す図である。自律神経指標データ処理部２３は、脈拍間隔データ処理部２２によって等間隔の脈拍間隔データが生成されると、ＦＦＴ法を用いて周波数解析を行って一連の脈拍間隔データを周波数スペクトル分布に変換する。周波数解析法は、ＡＲモデル、最大エントロピー法、ウェーブレット法などを用いても良いが、データ処理の負担の軽いＦＦＴ法を用いている。
【００２７】
図１の説明に戻ると、睡眠状態判定部２４は、複数のパワースペクトルの値のうち、低周波数領域（０．０５〜０．１５Ｈｚ付近）にあるパワースペクトルの値をＬＦとし、高周波数領域（０．１５〜０．４Ｈｚ付近）にあるパワースペクトルの値をＨＦとして、ＬＦとＨＦの値によって被験者の睡眠状態を判定する判定部である。具体的には、ＬＦ／ＨＦの値が第１の判定閾値よりも小さく、かつ、ＨＦの値が第２の判定閾値２よりも大きいときは深睡眠と判定し、また、ＬＦ／ＨＦの値が第３の判定閾値よりも大きく、かつ、ＨＦの値が第４の判定閾値より小さく、かつ、ＬＦ、ＨＦの標準偏差の合計が第５の判定閾値より大きいときはレム睡眠と判定し、深睡眠およびレム睡眠以外のときは浅睡眠と睡眠状態を判定する。
【００２８】
なお、第１の判定閾値から第５の判定閾値は、例えば、被験者毎に一晩計測したＬＦ，ＨＦ，ＬＦ／ＨＦ，ＬＦ，ＨＦの標準偏差の合計のそれぞれの分布の密度の高い点を２点選び、ＬＦ／ＨＦの２点の中点を第１の判定閾値＝第３の判定閾値、ＨＦの２点の中点を第２の判定閾値＝第４の判定閾値、ＬＦ，ＨＦの標準偏差の合計の中点を第５の判定閾値として設定することができる。
【００２９】
また、睡眠状態判定部２４は、ＬＦ，ＨＦの標準偏差の合計の中点の値が第５の判定閾値よりも小さく、かつ、ＨＦの値が第２の判定閾値より大きいときに深睡眠と判定するか、または、ＬＦ／ＨＦの値が第１の判定閾値より小さく、かつ、ＨＦの値が第２の判定閾値よりも小さいときに深睡眠と睡眠状態を判定することができる。
【００３０】
また、睡眠状態判定部２４は、ＬＦ／ＨＦの値が第３の判定閾値より大きく、かつ、ＬＦ、ＨＦの標準偏差値の合計が第５の判定閾値より大きいときにレム睡眠と判定するか、または、ＬＦ／ＨＦの値が第３の判定閾値より大きく、かつ、ＬＦ，ＨＦの標準偏差の合計が第５の判定閾値よりも大きく、かつ、脈波の一周期の平均時間間隔が第６の判定閾値以下のときにレム睡眠と睡眠状態を判定することができる。
【００３１】
睡眠状態判定部２４は、時刻に応じた被験者のサーカディアンリズムの影響を考慮に入れて第１の判定閾値から前記第６の判定閾値までの判定閾値を変更することができる。
自律神経系の変動は、睡眠状態の変化と、生体の持つサーカディアンリズム（概日リズム）とによって影響を受けることが知られている。この影響を低減するために睡眠状態判定装置１０は、時刻計測部１７を内蔵し、時刻と対応付けて第１の判定閾値から第６の判定閾値までの判定閾値を変更することによりサーカディアンリズムの影響を低減し、睡眠状態の判定精度を向上することができる。
【００３２】
睡眠状態判定部２４は、以上のような判定に基づいて睡眠状態を総括するため、次の式で示すような快眠度を算出して、表示する。なお、次の式のａｂｓは、絶対値の演算記号である。
快眠度（％）＝（１−定数Ａ×ａｂｓ（深睡眠割合−標準深睡眠割合）／標準深睡眠割合
−定数Ｂ×ａｂｓ（レム睡眠割合−標準レム睡眠割合）／標準レム睡眠割合）×１００
（式２）
【００３３】
ここで、式２に示す標準深睡眠割合、標準レム睡眠割合について説明する。図８は、式２に示す標準深睡眠割合、標準レム睡眠割合を示す図である。同図に示すよう、標準深睡眠割合、標準レム睡眠割合が年齢、性別に対して統計的に知られているので、これらの値にもとづいて固定的に設定することができる。また、これらの値は、時間帯によって、または就寝後の経過時間によっても変化するので時間帯ごとに快眠度を取得して、平均を求めることができる。また、これらの値は、被験者の過去のデータの平均値を用いることができる。また、上記快眠度の算出には、体動量をパラメータとして加えることができる。
【００３４】
ここで、図１に示す睡眠状態判定部２４が判定した睡眠状態判定結果の表示の例について説明する。図９は、図１に示す睡眠状態判定部２４が判定した睡眠状態判定結果の表示の一例を示す図である。また、図１０は、図１に示す睡眠状態判定部２４が判定した睡眠状態判定結果の表示の別の例を示す図である。
【００３５】
図９に示すように、睡眠状態判定結果は、快眠度と睡眠時間で睡眠状態を総括し、深睡眠の割合、レム睡眠の割合、体動量の割合、中途覚醒の割合のそれぞれの標準的な割合からの差分を示している。すなわち、深睡眠の割合は、標準的な深睡眠の割合より５％多いことを示している。また、深睡眠の（  ）内の３５％は、睡眠時間７：３０に占める深睡眠時間の割合を示している。また、図１０に示すように、睡眠状態判定結果は、グラフで表示し、時間帯ごとに体動の頻度を示すことができる。
【００３６】
図１の説明に戻ると、制御部２５は、睡眠状態判定装置１０全体を制御する制御部であり、被験者の要求および指示を受け付けて各処理部に対する処理要求およびデータの流れを制御する。具体的には、被験者の要求を受け付けて電源のＯＮ／ＯＦＦ、睡眠状態判定機能の起動、および睡眠状態判定結果の表示などを制御する。
【００３７】
次に、図１に示す睡眠状態判定装置１０の睡眠状態判定手順について説明する。図１１は、図１に示す睡眠状態判定装置１０の睡眠状態判定手順を示すフローチャートである。同図に示すように、被験者は、睡眠前に本睡眠状態判定装置１０を装着し、入力部１１から電源および睡眠状態判定機能を起動すると（ステップＳ１１０１）、睡眠状態判定装置１０は、加速度センサおよび脈波センサによって加速度データと脈波データの測定を開始する（ステップＳ１１０２、ステップＳ１１１０）。
【００３８】
体動サンプリングデータ処理部２０は、体動計測部１４から取得した３軸方向の加速度データを時間微分して３軸方向の加速度の微係数を求め（ステップＳ１１０３）、３軸方向の加速度のそれぞれの微係数の二乗和の平方根を求める（ステップＳ１１０４）。そして、体動判定部１９１は、体動サンプリングデータ処理部２０によって求められた３軸方向の加速度の微係数の二乗和の平方根を受け取って、体動が現れるまで待つ（ステップＳ１１０５）。具体的には、体動サンプリングデータ処理部２０によって求められた３軸方向の加速度の微係数の二乗和の平方根である体動データの変動量が第１の所定の閾値より大きい場合に体動と判定する。例えば、第１の所定の閾値として体動計に使用されている微小な体動の最小値である０．０１Ｇを用いる。
【００３９】
そして、体動判定部１９１は、体動が現れた場合は（ステップ１１０５肯定）、体動が現れたことを波形異常検出部２２１に通知する。また、体動判定部１９１は、体動が現れた場合は（ステップ１１０５肯定）、覚醒判定部１９２は、設定区間、例えば、１分間の体動発生頻度を取得して（ステップ１１０６）、覚醒か否かを判定する（ステップＳ１１０７）。具体的には、覚醒判定部１９２は、体動判定部１９１によって判定された体動の発生頻度が第２の所定の閾値以上である場合に覚醒中の体動と判定し、体動の発生頻度が第２の閾値未満である場合は睡眠中の体動と判定する。例えば、第２の所定の閾値として過去の覚醒時における体動頻度から２０回／分を用いる。
【００４０】
その結果、覚醒判定部１９２が、覚醒状態と判定した場合は（ステップＳ１１０７肯定）、覚醒判定部１９２は、記憶部１３に入眠時刻、覚醒時刻、中途覚醒回数を保持させる（ステップ１１０８）。一方、覚醒判定部１９２が、睡眠中の体動と判定した場合は（ステップＳ１１０７否定）、覚醒判定部１９２は、体動と判定された体動データが睡眠中の体動であることを異常波形検出部２１１に通知する。また、覚醒判定部１９２は、記憶部１３に睡眠中の体動量を保持させる（ステップＳ１１０９）。
【００４１】
これに対して、脈波サンプリング取得部２１は、脈波から脈波データをサンプリングし、サンプリングした一連の脈波データを時間微分して、一連の脈波データから直流変動成分を除去する（ステップＳ１１１１）。そして、直流変動成分を除去された一連の脈波データの所定の時間間隔内の脈波データの最大値と最小値を取得し、最大値と最小値との間の所定の値を第３の所定の閾値として決定する（ステップＳ１１１２）。例えば、第３の所定の閾値として最大値、最小値の差を振幅として、最小値から振幅の９割の値を用いる。さらに、直流変動成分を除去された一連の脈波データから第３の所定の閾値に一致する一連の脈波データの値が現れた時刻を算出し、算出された時刻の間隔から脈拍間隔データを取得する（ステップＳ１１１３）。
【００４２】
そして、異常波形検出部２２１は、体動判定部１９１および覚醒判定部１９２からの通知によって体動から睡眠中の体動を識別し、睡眠中の体動と判定された体動データの変動量を脈拍間隔内で平均した体動量が第５の所定の閾値以上の場合に、体動量が大きいと判定し、睡眠中の体動と判定された体動データと並行して計測された脈拍間隔データを異常波形として検出する。さらに、異常データ除外部２２２は、覚醒中の体動と判定された体動データと並行して計測された脈拍間隔データと波形異常検出部２２１によって異常波形として検出された脈拍間隔データを一連の脈拍間隔データから除外する（ステップＳ１１１４）。
【００４３】
さらに、脈拍間隔データ処理部２２は、異常データ除外部２２２によって異常波形が除外された一連の脈拍間隔データを用いて１分間のデータセット生成する（ステップＳ１１１５）。そして、脈拍間隔データ補間部２２３は、異常データ除外部２２２によって異常波形が除外された一連の脈拍間隔データをスプライン補間、すなわち、高次の多項式で補間し、等間隔の脈拍間隔データを生成する（ステップＳ１１１６）。
【００４４】
さらに、周波数スペクトル変換部２３１は、脈拍間隔データ処理部２２によってデータ処理された一連の脈拍間隔データをＦＦＴ法などの周波数解析法によって周波数スペクトル分布に変換する（ステップ１１１７）。そして、自律神経指標取得部２３２は、周波数スペクトル変換部２３１によって周波数スペクトル分布に変換された一連の脈拍間隔データの複数のパワースペクトルの値から自律神経指標ＬＦ，ＨＦを取得する（ステップＳ１１１８）。具体的には、複数のパワースペクトルのピーク値とピーク値を中心として前後等間隔の１点との３点の合計値の算術平均をとってＬＦ、ＨＦとする。
【００４５】
さらに、睡眠状態判定部２４は、自律神経指標ＬＦ，ＨＦに基づいて睡眠状態を判定し、記憶部１３に保持させる（ステップＳ１１１９）。そして、睡眠状態判定部２４は、被験者の睡眠状態の総括として快眠度を算出し、記憶部１３に保持させる（ステップＳ１１２０）。快眠度の表示を覚醒判定時に自動的に行う。さらに、睡眠状態判定部２４は、被験者から睡眠状態判定結果の表示要求があったか否かを調べる（ステップＳ１１２１）。その結果、被験者から睡眠状態判定結果の表示要求があった場合は（ステップＳ１１２１肯定）、表示部１２に被験者の要求に応じて記憶部１３に保持された睡眠状態判定結果を表示させる（ステップＳ１１２２）。一方、被験者から睡眠状態判定結果の表示要求がなかった場合は（ステップＳ１１２１否定）、ステップＳ１１０２に戻り、ステップＳ１１２０までの手順を繰り返す。
【００４６】
ここで、図１１に示す自律神経指標ＬＦ，ＨＦに基づく睡眠状態判定手順をさらに詳細に説明する。図１２は、図１１に示す自律神経指標ＬＦ，ＨＦに基づく睡眠状態判定手順をさらに詳細に示すフローチャートである。同図に示すように、睡眠状態判定部２４は、睡眠状態判定のために自律神経指標ＬＦ、ＨＦの標準偏差の合計を取得し（ステップＳ１２０１）、ＬＦ／ＨＦを取得する（ステップＳ１２０２）。
【００４７】
そして、睡眠状態判定部２４は、ＬＦ／ＨＦの値が第１の判定閾値よりも小さいか否かを調べる（ステップＳ１２０３）。その結果、ＬＦ／ＨＦの値が第１の判定閾値よりも小さい場合は（ステップＳ１２０３肯定）、さらに、ＨＦの値が第２の判定閾値よりも大きいか否かを調べる（ステップＳ１２０５）。その結果、ＨＦの値が第２の判定閾値よりも大きい場合は（ステップＳ１２０５肯定）、深睡眠と判定する（ステップＳ１２０９）。
【００４８】
一方、睡眠状態判定部２４は、ＬＦ／ＨＦの値が第１の判定閾値以上である場合は（ステップＳ１２０３否定）、さらに、ＬＦ／ＨＦの値が第３の判定閾値より大きいか否かを調べる（ステップＳ１２０４）。その結果、ＬＦ／ＨＦの値が第３の判定閾値より大きい場合は（ステップＳ１２０４肯定）、さらに、ＨＦの値が第２の判定閾値よりも大きいか否かを調べる（ステップＳ１２０５）。その結果、ＨＦの値が第２の判定閾値以下である場合は（ステップＳ１２０５否定）、さらに、ＨＦの値が第４の判定閾値よりも小さいか否かを調べる（ステップＳ１２０６）。その結果、ＨＦの値が第４の判定閾値よりも小さい場合は（ステップＳ１２０６肯定）、さらに、ＬＦ、ＨＦの標準偏差の合計が第５の判定閾値より大きいか否かを調べる（ステップＳ１２０７）。その結果、ＬＦ、ＨＦの標準偏差の合計が第５の判定閾値より大きい場合は（ステップＳ１２０７肯定）、レム睡眠と睡眠状態を判定する（ステップＳ１２０８）。
【００４９】
一方、睡眠状態判定部２４は、ＬＦ／ＨＦの値が第２の判定閾値以下である場合（ステップＳ１２０４否定）、および、ＨＦが第４の判定閾値以上である場合（ステップＳ１２０６否定）、および、ＬＦ、ＨＦの標準偏差の合計が第５の判定閾値以下である場合は（ステップＳ１２０７否定）、浅睡眠と睡眠状態を判定する（ステップＳ１２１０）。
【００５０】
なお、第１の判定閾値から第５の判定閾値は、例えば、被験者毎に一晩計測したＬＦ，ＨＦ，ＬＦ／ＨＦのそれぞれの分布の密度の高い点を２点選び、ＬＦ／ＨＦの２点の中点を第１の判定閾値＝第３の判定閾値、ＨＦの２点の中点を第２の判定閾値＝第４の判定閾値、ＬＦの２点の中点を第５の判定閾値として設定することができる。
【００５１】
上述してきたように、本実施の形態１では、体動データの変動が第１の所定の閾値より大きい場合に体動と判定し、体動と判定された場合に体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理し、一連の脈拍間隔データから取得した自律神経指標に基づいて前記睡眠状態を判定することとしたので、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００５２】
また、睡眠中の体動と判定された場合に睡眠中の体動と判定された体動データと並行して計測された脈拍間隔データを一連の脈拍間隔データから除外してデータ処理し、一連の脈拍間隔データから取得された自律神経指標に基づいて睡眠状態を判定することとしたので、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００５３】
また、体動と判定された体動データと並行して計測された脈拍間隔データを除外された一連の脈拍間隔データを高次の多項式で補間することとしたので、異常データが除外されて不等間隔になった一連の脈拍間隔データを用いて周波数スペクトル分布に変換することができ、睡眠状態の判定に必要な自律神経指標を取得できる。したがって、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００５４】
また、複数のパワースペクトルの値のうち、低周波数の領域にあるパワースペクトルの値をＬＦ、高周波数の領域にあるパワースペクトルの値をＨＦとしたときに、ＬＦ／ＨＦの値が第１の判定閾値よりも小さく、かつ、ＨＦの値が第２の判定閾値２よりも大きいときは深睡眠と判定し、また、ＬＦ／ＨＦの値が第３の判定閾値よりも大きく、かつ、ＨＦの値が第４の判定閾値より小さく、かつ、ＬＦ、ＨＦの標準偏差の合計が第５の判定閾値より大きいときはレム睡眠と判定し、深睡眠およびレム睡眠以外は浅睡眠と睡眠状態を判定することとしたので、自律神経指標を用いて深睡眠、レム睡眠、浅睡眠が容易に識別できる。したがって、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００５５】
また、３軸方向の加速度データを体動データとして計測することとしたので、体動を手軽で精度よく体動を測定することができる。したがって、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００５６】
（実施の形態２）
ところで、実施の形態１では、本発明に係る睡眠状態判定装置１０を加速度センサ１４１を用いて脈波センサ１５１への体動の影響を低減する場合について説明したが、本発明はそれに限定されるものではなく、本発明に係る睡眠状態判定装置１０を睡眠中の不整脈や無呼吸などの異常を検出する場合に適用することができる。そこで、本実施の形態２では、本発明に係る睡眠状態判定装置１０を睡眠中の不整脈や無呼吸などの異常を検出する場合について説明する。なお、本実施の形態２のでは、実施の形態１と共通な部分については説明を省略する。
【００５７】
まず、本実施の形態２に係る睡眠状態判定装置の構成について説明する。図１３は、本実施の形態２に係る睡眠状態判定装置の構成を示す機能ブロック図である。同図に示すように、図１に示す実施の形態１の機能ブロック図との構成上の相違はないが、脈拍間隔データ処理部２２の波形異常検出部２２１の機能が相違する。
【００５８】
本実施の形態２に係る睡眠状態判定装置の波形異常検出部２２１は、取得した脈拍間隔データの波形と正常な脈拍間隔データの波形との一致度を算出して一致度が第４の所定の閾値未満である場合は、取得した脈拍間隔データの波形が異常であることを検出する処理部である。具体的には、脈拍間隔データの波形と正常な脈拍間隔データの波形との相関係数を算出して相関係数が第４の所定の閾値未満の場合は、脈拍間隔データの波形異常を検出する。例えば、第４の所定の閾値として０．５を用いる。また、波形異常検出部２２１は、脈拍間隔データの波形と異常な脈拍間隔データの波形とをＤＰ（Ｄｙｎａｍｉｃ  Ｐｒｏｇｒａｍｍｉｎｇ）マッチングして、波形識別と発生頻度を検出する。ここで、図１３に示す睡眠状態判定部２４が判定した睡眠状態判定結果の例について説明する。図１４は、図１３に示す睡眠状態判定部２４が判定した睡眠状態判定結果の一例を示す図である。また、図１５は、図１３に示す睡眠状態判定部２４が判定した睡眠状態判定結果の別の例を示す図である。
【００５９】
図１４に示すように、睡眠状態判定結果は、図９の睡眠状態判定結果に無呼吸の割合、不整脈の割合が追加されている。これらの結果は、波形異常検出部２２１によって検出された異常である。また、図１５に示すように、睡眠状態判定結果は、時間帯ごとに覚醒、レム睡眠、深睡眠、浅睡眠などの睡眠状態と体動、無呼吸および不整脈の頻度とをグラフで表示することができる。
【００６０】
次に、図１３に示睡眠状態判定装置の睡眠状態判定手順について説明する。図１６は、図１３に示す睡眠状態判定装置の睡眠状態判定手順を示すフローチャートである。同図に示すように、本手順は、ステップＳ１６１４〜ステップＳ１６１５だけが実施の形態１の睡眠状態判定手順と相違し、その他の手順は全く同じなので、ここではステップＳ１６１４〜ステップＳ１６１５について説明する。
【００６１】
波形異常検出部２２１は、覚醒判定部１９２によって睡眠中の体動と判定された脈拍間隔データの波形と正常な脈拍間隔データの波形との一致度を算出して一致度が第４の所定の閾値未満である場合は、脈拍間隔データの波形が異常であることを検出する（ステップＳ１６１４）。具体的には、覚醒判定部１９２によって睡眠中の体動と判定された脈拍間隔データの波形と正常な脈拍間隔データの波形との相関係数を算出して相関係数が第４の所定の閾値未満の場合は、脈拍間隔データの波形が異常であることを検出する。例えば、第４の所定の閾値として０．５を用いる。
【００６２】
そして、波形異常検出部２２１は、覚醒判定部１９２によって体動と判定された脈拍間隔データの波形と睡眠中の不整脈や無呼吸などの異常な脈拍間隔データの波形とをマッチングして、波形の識別と発生頻度を検出し、集計する（ステップＳ１６１５）。さらに、脈拍間隔データ処理部２２は、波形異常検出部２２１によって検出された波形異常の発生頻度を取得し、睡眠状態判定が可能か否かを調べる（ステップＳ１６１６）。
【００６３】
その結果、波形異常の発生頻度が多く、所定のデータ数が確保できないために睡眠状態判定ができない場合は（ステップＳ１６１６否定）、睡眠状態判定部２４は、体動または波形異常と判定する（ステップＳ１６２４）。一方、波形異常の発生頻度が少なく、所定のデータ数が確保できるので睡眠状態判定ができる場合は（ステップＳ１６１６肯定）、ステップＳ１６１５以降の手順を行う。
【００６４】
上述してきたように、本実施形態２では、覚醒判定部１９２によって睡眠中の体動と判定された脈拍間隔データの波形と正常な脈拍間隔データの波形との一致度を算出して一致度が第４の所定の閾値未満である場合は、覚醒判定部１９２によって睡眠中の体動と判定された脈拍間隔データを波形異常として検出することとしたので、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００６５】
また、覚醒判定部１９２によって睡眠中の体動と判定された脈拍間隔データの波形と正常な脈拍間隔データの波形との相関係数を算出して相関係数が第４の所定の閾値未満の場合は、覚醒判定部１９２によって睡眠中の体動と判定された脈拍間隔データを波形異常として検出することとしたので、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００６６】
また、覚醒判定部１９２によって睡眠中の体動と判定された脈拍間隔データの波形と異常な脈拍間隔データの波形とをマッチングして、波形の識別と発生頻度を検出することとしたので、脈波に対する体動の影響および不整脈や無呼吸状態などの脈波異常の影響を低減し、もって、睡眠状態の判定精度を向上させることができる。
【００６７】
（実施の形態３）
ところで、上記の実施の形態１および２では、本発明に係る睡眠状態判定装置１０を体動センサとして加速度センサ１４１を用いる場合についてしたが、本発明はそれに限定されるものではなく、本発明に係る睡眠状態判定装置を体動センサとしてマットセンサを用いる場合に適用することができる。ここでは、本発明に係る睡眠状態判定装置を体動センサとしてマットセンサを用いる場合について説明する。なお、本実施の形態３のでは、実施の形態１と共通な部分については説明を省略する。
【００６８】
まず、本実施の形態３に係る睡眠状態判定装置１０の構成について説明する。図１７は、本実施の形態３に係る睡眠状態判定装置１０の構成を示す機能ブロック図である。同図に示すように、図１に示す実施の形態１の機能ブロック図との構成上の相違は、体動計測部２６とマットセンサ２６１だけである。ここで、図１７に示す睡眠状態判定装置１０のマットセンサの配置について説明する。図１８は、図１７に示す睡眠状態判定装置１０のマットセンサの配置の一例を示す図である。
【００６９】
同図に示すように、マットセンサ２６１は、ベッド２７のマットレス２８の表面に設置され、被験者の不在、在床、体動を検出するセンサであり、具体的には、ポリフッ化ビニリデンなどの高分子圧電材料を薄膜状にし、両面に可とう性の電極膜を付着させてテープ状に形成した圧電素子である。また、体動計測部２６は、マットセンサ２６１の出力信号をフィルタ、増幅器を介してＡ／Ｄ変換器でデジタル量に変換する計測部である。
【００７０】
本発明に係る睡眠状態判定装置１０は、体動計測部２６がマットセンサ２６１の出力信号をデジタル量に変換した後は実施の形態１および２と同じ処理手順で睡眠状態を判定することができる。したがって、実施の形態１および２に係る睡眠状態判定装置１０では、体動センサとして加速度センサ１４１を使用したが、本発明はこれに限定されるものではなく、マットセンサ２６１などの他の体動センサを使用することができる。
【００７１】
（実施の形態４）
また、図１９を参照して、本実施の形態４に係る睡眠状態判定装置１０の構成について説明する。図１９は、本実施の形態４に係る睡眠状態判定装置１０の構成を示す機能ブロック図である。同図に示すように、睡眠状態判定装置１０は、マットセンサ２６１内に体動抽出部２６１１、脈波抽出部２６１２を設けた構成でも適用できる。体動抽出部２６１１は前述したようにベッド２７のマットレス２８の表面に設置された体動等を検出するセンサからの体動成分を抽出する。一方、脈波抽出部２６１２はマットセンサの出力から周波数解析等で脈波成分を抽出する。この抽出した体動成分と脈波成分はそれぞれ体動計測部２６及び脈波計測部１５に入力され、実施の形態１及び２と同様な処理手順で睡眠状態が判定される。
【産業上の利用可能性】
【００７２】
以上のように、本発明に係る睡眠状態判定装置は、日常生活で手軽に使用でき、判定精度の高い睡眠状態判定装置に有用であり、特に、携帯用の睡眠状態判定装置に適している。
【図面の簡単な説明】
【００７３】
【図１】本発明に係る睡眠状態判定装置の構成を示す機能ブロック図である。
【図２】図１に示す睡眠状態判定装置の装着の一例を示す図である。
【図３】図１に示す睡眠状態判定装置の装着の別の例を示す図である。
【図４】図１に示す脈拍間隔データ取得部が脈拍間隔データを算出する算出方法を示す図である。
【図５】図１に示す異常データ除外部が体動と判定された脈拍間隔データを除外した一連の脈拍間隔データの一例を示す図である。
【図６】図１に示す脈拍間隔データ補間部が一連の脈拍間隔データを補間する一例を示す図である。
【図７】図１に示す自律神経指標データ処理部がパワースペクトルの値から自律神経指標を取得する一例を示す図である。
【図８】式２に示す標準深睡眠割合、標準レム睡眠割合を示す図である。
【図９】図１に示す睡眠状態判定部が判定した睡眠状態判定結果の表示の一例を示す図である。
【図１０】図１に示す睡眠状態判定部が判定した睡眠状態判定結果の表示の別の例を示す図である。
【図１１】図１に示す睡眠状態判定装置の睡眠状態判定手順を示すフローチャートである。
【図１２】図１１に示す自律神経指標に基づく睡眠状態判定手順をさらに詳細に示すフローチャートである。
【図１３】本実施の形態２に係る睡眠状態判定装置の構成を示す機能ブロック図である。
【図１４】図１３に示す睡眠状態判定部が判定した睡眠状態判定結果の一例を示す図である。
【図１５】図１３に示す睡眠状態判定部が判定した睡眠状態判定結果の別の例を示す図である。
【図１６】図１３に示す睡眠状態判定装置の睡眠状態判定手順を示すフローチャートである。
【図１７】本実施の形態３に係る睡眠状態判定装置の構成を示す機能ブロック図である。
【図１８】図１７に示す睡眠状態判定装置のマットセンサの配置の一例を示す図である。
【図１９】本実施の形態４に係る睡眠状態判定装置の構成を示す機能ブロック図である。
【符号の説明】
【００７４】
１０    睡眠状態判定装置
１１    入力部
１２    表示部
１３    記憶部
１４、２６    体動計測部
１５    脈波計測部
１６    データ通信部
１７    時刻計測部
１８    電力供給部
１９    覚醒／睡眠判定部
２０    体動サンプリングデータ取得部
２１    脈波サンプリングデータ取得部
２２    脈拍間隔データ処理部
２３    自律神経指標データ処理部
２４    睡眠状態判定部
２５    制御部
２７    ベッド
２８    マットレス
１９１    体動判定部
１９２    覚醒判定部
２０１    脈波間隔データ取得部
２２１    異常データ除外部
２２２    脈拍間隔データ補間部
２２３    波形異常検出部
２３１    周波数スペクトル変換部
２３２    自律神経指標取得部
２６１    マットセンサ
２６１１  体動抽出部
２６１２  脈波抽出部
【図１】
2005279113.tif 000003
【図２】
2005279113.tif 000004
【図３】
2005279113.tif 000005
【図４】
2005279113.tif 000006
【図５】
2005279113.tif 000007
【図６】
2005279113.tif 000008
【図７】
2005279113.tif 000009
【図８】
2005279113.tif 000010
【図９】
2005279113.tif 000011
【図１０】
2005279113.tif 000012
【図１１】
2005279113.tif 000013
【図１２】
2005279113.tif 000014
【図１３】
2005279113.tif 000015
【図１４】
2005279113.tif 000016
【図１５】
2005279113.tif 000017
【図１６】
2005279113.tif 000018
【図１７】
2005279113.tif 000019
【図１８】
2005279113.tif 000020
【図１９】
2005279113.tif 000021