(57)【要約】  在宅で睡眠状態についての測定を行うために、使用者の身体に接している圧力検出部から圧力波形を得る圧力波形取得部１０と、圧力波形に所定の処理を行って拍動波形を抽出する拍動抽出部２２と、拍動間隔を算出する拍動間隔算出部２４と、前記使用者の睡眠の深度の段階である睡眠ステージを拍動間隔から判定する睡眠ステージ判定部２６とを備える睡眠状態モニタリング装置を用い、また、コンピュータネットワークを通じて睡眠中の拍動間隔のデータを収集して健康管理を行う。
2009150744.tif 000002
【特許請求の範囲】
【請求項１】
使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、
該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、
該拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、
前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定部と
を備える睡眠状態モニタリング装置。
【請求項２】
前記圧力波形に所定の処理を行って呼吸波形を抽出する呼吸波形抽出部と、
該呼吸波形における各回の呼吸について、直前の呼吸からの時間間隔である呼吸間隔を算出する呼吸間隔算出部と、
使用者が睡眠しているか覚醒しているかを該呼吸間隔から判定する覚醒判定部と
をさらに備える請求項１に記載の睡眠状態モニタリング装置。
【請求項３】
前記覚醒判定部が、前記拍動間隔から拍動間隔の変動係数を算出し、該拍動間隔の変動係数によって睡眠か覚醒かを判定するものである、請求項２に記載の睡眠状態モニタリング装置。
【請求項４】
前記覚醒判定部が、前記呼吸間隔から呼吸間隔の変動係数を算出し、該呼吸間隔の変動係数によって睡眠か覚醒かを判定する、請求項２に記載の睡眠状態モニタリング装置。
【請求項５】
前記拍動間隔算出部は、微分処理とピーク判定とを組み合わせて前記拍動波形から各拍動のタイミングを示す時刻を特定する微分ピーク処理部を有する、請求項１に記載の睡眠状態モニタリング装置。
【請求項６】
前記睡眠ステージ判定部は、
拍動間隔の時系列データから拍動間隔のゆらぎについてのパワースペクトルを算出するパワースペクトル算出部と、
該パワースペクトルの示すパターンを判定対象として睡眠ステージを判定するスペクトルパターン解析部と
を備えている、請求項１に記載の睡眠状態モニタリング装置。
【請求項７】
前記スペクトルパターン解析部が、ニューラルネットワークの学習によるパターン判定部を備える、請求項６に記載の睡眠状態モニタリング装置。
【請求項８】
前記睡眠ステージ判定部は、拍動間隔の時系列データに含まれる拍動間隔のゆらぎについての低周波成分と高周波成分とを算出して該低周波成分と該高周波成分との比を計算するＬＦ／ＨＦ比算出部を備え、該比の値によって睡眠ステージを判定するものである、請求項１に記載の睡眠状態モニタリング装置。
【請求項９】
前記圧力波形取得部が、使用者の身体を支える寝具と身体との間に配置される空気圧センサーを含む、請求項１に記載の睡眠状態モニタリング装置。
【請求項１０】
使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得ステップと、
該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出ステップと、
該拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出ステップと、
前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定ステップと
を含む睡眠状態モニタリング方法
【請求項１１】
前記圧力波形に所定の処理を行って呼吸波形を抽出する呼吸波形抽出ステップと、
該呼吸波形における各回の呼吸について、直前の呼吸からの時間間隔である呼吸間隔を算出する呼吸間隔算出ステップと、
使用者が睡眠しているか覚醒しているかを呼吸間隔から判定する覚醒判定ステップと
をさらに含む請求項１０に記載の装置。
【請求項１２】
前記覚醒判定ステップが、前記拍動間隔から拍動間隔の変動係数を算出し、該拍動間隔の変動係数によって睡眠か覚醒かを判定するものである、請求項１１に記載の睡眠状態モニタリング方法。
【請求項１３】
前記覚醒判定ステップが、前記呼吸間隔から呼吸間隔の変動係数を算出し、該呼吸間隔の変動係数によって睡眠か覚醒かを判定するものである、請求項１１に記載の睡眠状態モニタリング方法。
【請求項１４】
前記拍動間隔算出ステップは、微分処理とピーク判定とを組み合わせて前記拍動波形から各拍動のタイミングを示す時刻を特定する微分ピーク処理ステップを含む、請求項１０に記載の睡眠状態モニタリング方法。
【請求項１５】
前記睡眠ステージ判定ステップは、
拍動間隔の時系列データから拍動間隔のゆらぎについてのパワースペクトルを算出するパワースペクトル算出ステップと、
該パワースペクトルの示すパターンを判定対象として睡眠ステージを判定するスペクトルパターン解析ステップと
を含む、請求項１０に記載の睡眠状態モニタリング方法。
【請求項１６】
前記スペクトルパターン解析ステップが、学習を行わせたニューラルネットワークにパターンを判定させるステップを含む、請求項１５に記載の装置。
【請求項１７】
前記睡眠ステージ判定ステップは、拍動間隔の時系列データに含まれる拍動間隔のゆらぎについての低周波成分と高周波成分とを算出して該低周波成分と該高周波成分との比を計算するＬＦ／ＨＦ比算出ステップを含み、該比の値によって睡眠ステージを判定するステップである、請求項１０に記載の装置。
【請求項１８】
前記圧力波形取得ステップが、使用者の身体を支える寝具と身体との間に配置される空気圧センサーによって圧力波形を取得するステップを含む、請求項１０に記載の睡眠状態モニタリング方法。
【請求項１９】
使用者の身体のいずれかの部分に直接または着衣を介して接しており電気信号として圧力波形を得る圧力検出部に接続されたコンピュータに、請求項１０〜１８に記載のいずれかの方法を実行させるためのコンピュータプログラム。
【請求項２０】
使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、コンピュータネットワークを通じてデータを送信するデータ送信部とを備えたコンピュータネットワークに接続されたモニタリング部と、
該コンピュータネットワークを通じて該装置から少なくとも拍動間隔を含む時系列データを複数回受信し、受信した時系列データを使用者に対応させて該時系列データを記録するための記憶部と
を有する睡眠状態モニタリングシステムであって、前記コンピュータネットワークに接続されたいずれかのコンピュータが、複数回にわたって受信した前記時系列データを使用者に対応付けて時間に対するグラフとして表示し、過去のある期間における該使用者の拍動間隔のゆらぎの変化を提示するデータ表示部を備えている、睡眠状態モニタリングシステム。
【請求項２１】
前記コンピュータネットワークに接続されたいずれかのコンピュータが、過去のある期間における該使用者の拍動間隔のデータにおいて、拍動間隔のゆらぎを示す数値についての前記期間より短い短期間における変動を検知してアラーム信号を出力するゆらぎ変化探知部を有している請求項２０に記載の睡眠状態モニタリングシステム。
【請求項２２】
使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、コンピュータネットワークを通じてデータを送信するデータ送信部とを備え、コンピュータネットワークに接続されているモニタリング部に該コンピュータネットワークを通じて接続可能にされたコンピュータを、請求項２０または２１のいずれかに記載のコンピュータとして動作させるためのコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、睡眠状態モニタリング装置、モニタリングシステムおよびコンピュータプログラムに関する。
【背景技術】
【０００２】
現在、人の健康管理を行ったり心身へのストレスを適切に管理するための手法のひとつとして、睡眠の質を管理することが求められている。睡眠の質を具体的に評価するには、睡眠を６つのカテゴリ、すなわち、覚醒、ＲＥＭ（Ｒａｐｉｄ  Ｅｙｅ  Ｍｏｖｅｍｅｎｔｓ）睡眠、ステージ１（Ｓ１）〜ステージ４（Ｓ４）に分類して評価することが行われており、このためには、通常、電気的計測が利用されている。なお、覚醒、ＲＥＭ睡眠、Ｓ１〜Ｓ４、は、この順に睡眠の深度が大きくなるが、通常の睡眠の状態の推移では、睡眠開始直後にＲＥＭ睡眠になることはなく、一旦深い眠り（Ｓ３またはＳ４）になってから覚醒に向かって眠りが浅くなる際にＲＥＭ睡眠の状態になる。この測定では、人（使用者または被験者）の身体のいくつかの部位に電極（皿電極）を装着して、例えば、頭皮からの脳波、アゴからの筋電位、眼の周辺からの眼球電位、および心電図を利用する必要がある。脳波から睡眠状態を判定する技術が、Ｒｅｃｈｔｓｃｈａｆｆｅｎ−Ｋａｌｅｓ（Ｒ−Ｋ）の方法として知られている。また、心電図を利用する技術は、例えば特許文献１（特開２００３−２２５２１１号公報）に開示されている。これらの測定手法による睡眠状態の評価は精度が高いという特徴を有するものの、正しい評価結果を得るためには、専門的知識を有する専門家または医療従事者がデータの解析を行う必要がある。また、電気的測定のために電極を測定時に装着するため、使用者に対する拘束性が高く、睡眠中の使用者に対して精神的な負担を与えてしまう。このため、これらの従来の技術は、睡眠の状態に影響すると共に、長期間にわたる継続的な測定を行うことも難しい。【０００３】
また、睡眠状態の評価を上記以外の生体情報によって実現する手法について、多くの報告がなされている。この手法は、体振動の情報の採り方によりいくつかのバリエーションがあるが、その一つとして、体振動の大きさを指標として定義する方法（非特許文献１）が知られている。【０００４】
さらに、自律神経成分比の研究によって睡眠ステージをごく大きく分類することが報告されている。非特許文献３には、心電図波形に基づいて自律神経成分比を求めて、浅い睡眠（Ｓ１またはＳ２）と、深い睡眠（Ｓ３またはＳ４）との間で判定を行うことが開示されている。【０００５】
そして、 非特許文献５には、心疾患における心拍のスペクトルの変化について開示されており、また、非特許文献６には、心拍変動のパワースペクトルを応用して、糖尿病による初期段階の自律神経機能不全を検出しうることについて開示されている。
【特許文献１】特開２００３−２２５２１１号公報
【非特許文献１】渡辺崇士、渡辺嘉二郎：  就寝時無拘束計測生体データによる睡眠段階の推定、計測自動制御学会論文集、３８（７），  ５８１−５８９，  （２００２）
【非特許文献２】早野順一郎他：  心拍変動と自律神経機能、生物物理、２８，  １９８―２０２，  （１９８８）
【非特許文献３】Ｂａｈａｒａｖ  Ａ，  ｅｔ  ａｌ．  ：  Ｆｌｕｃｔｕａｔｉｏｎｓ  ｉｎ  ａｕｔｏｎｏｍｉｃ  ｎｅｒｖｏｕｓ  ａｃｔｉｖｉｔｙ  ｄｕｒｉｎｇ  ｓｌｅｅｐ  ｄｉｓｐｌａｙｅｄ  ｂｙ  ｐｏｗｅｒ  ｓｐｅｃｔｒｕｍ  ａｎａｌｙｓｉｓ  ｏｆ  ｈｅａｒｔ  ｒａｔｅ  ｖａｒｉａｂｉｌｉｔｙ，  Ｎｅｕｒｏｌｏｇｙ，  ４５，  １１８３―１１８７，  （１９９５）
【非特許文献４】Ｔａｓｋ  Ｆｏｒｃｅ  ｏｆ  ｔｈｅ  Ｅｕｒｏｐｅａｎ  Ｓｏｃｉｅｔｙ  ｏｆ  Ｃａｒｄｉｏｌｏｇｙ  ａｎｄ  ｔｈｅ  Ｎｏｒｔｈ  Ａｍｅｒｉｃａｎ  Ｓｏｃｉｅｔｙ  ｏｆ  Ｐａｃｉｎｇ  ａｎｄ  Ｅｌｅｃｔｒｏｐｈｙｓｉｏｌｏｇｙ：  Ｈｅａｒｔ  Ｒａｔｅ  Ｖａｒｉａｂｉｌｉｔｙ（ＨＲＶ）：Ｓｔａｎｄａｒｄｓ  ｏｆ  Ｍｅａｓｕｒｅｍｅｎｔ，  Ｐｈｙｓｉｏｌｏｇｉｃａｌ  Ｉｎｔｅｒｐｒｅｔａｔｉｏｎ，  ａｎｄ  Ｃｌｉｎｉｃａｌ  Ｕｓｅ，  Ｃｉｒｃｕｌａｔｉｏｎ，  ９３（５），  １０４３−１０６５，  （１９９６）
【非特許文献５】Ｓａｕｌ  ＪＰ，  Ａｒａｉ  Ｙ，  Ｂｅｒｇｅｒ  ＲＤ．  ｅｔ  ａｌ.：  Ａｓｓｅｓｓｍｅｎｔ  ｏｆ  Ａｕｔｏｎｏｍｉｃ  ｒｅｇｕｌａｔｉｏｎ  ｉｎ  ｃｈｒｏｎｉｃ  ｃｏｎｇｅｓｔｉｖｅ  ｈｅａｒｔ  ｆａｉｌｕｒｅ  ｂｙ  ｈｅａｒｔ  ｒａｔｅ  ｓｐｅｃｔｒａｌ  ａｎａｌｙｓｉｓ，  Ａｍｅｒｉｃａｎ  Ｊｏｕｒｎａｌ  ｏｆ  Ｃａｒｄｉｏｌｏｇｙ，  ６１，  １２９２−１２９９，  （１９８８）
【非特許文献６】Ｙａｍａｓａｋｉ，  Ｙ，  Ｕｅｄａ，  Ｎ．，  Ｋｉｓｈｉｍｏｔｏ，  Ｍ  ｅｔ  ａｌ．：  Ａｓｓｅｓｓｍｅｎｔ  ｏｆ  ｅａｒｌｙ  ｓｔａｇｅ  ａｕｔｏｎｏｍｉｃ  ｎｅｒｖｅ  ｄｙｓｆｕｎｃｔｉｏｎ  ｉｎ  ｄｉａｂｅｔｉｃ  ｓｕｂｊｅｃｔｓ；  ａｐｐｌｉｃａｔｉｏｎ  ｏｆ  ｐｏｗｅｒ  ｓｐｅｃｔｒａｌ  ａｎａｌｙｓｉｓ  ｏｆ  ｈｅａｒｔ  ｒａｔｅ  ｖａｒｉａｂｉｌｉｔｙ，  Ｄｉａｂｅｔｅｓ  Ｒｅｓｅａｒｃｈ，  １７（２）,  ７３−８０，  （１９９１）
【発明の開示】
【発明が解決しようとする課題】
【０００６】
専門家または医療従事者の助けを必要とせずに使用者が在宅のまま睡眠の質を簡単にチェックできるような装置や方法があれば、自己の健康状態を把握することが可能となり、健康管理が容易になる。また、睡眠の質の低下が手軽に検出できれば、専門家または医療従事者が診断するための的確な情報を与えることができる。また、非特許文献３の手法では、心電図を測定する必要があるばかりでなく、睡眠ステージを細分するような判定は難しく、また、睡眠か覚醒かを判定することもできない。さらに、在宅のままの使用者の拍動間隔のゆらぎのデータを収集して使用者の健康状態を監視することも容易に行うことはできない。本発明は、上記いくつかの課題を解決することを課題とする。
【課題を解決するための手段】
【０００７】
本発明においては、使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、該拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定部とを備える睡眠状態モニタリング装置が提供される。また本発明の睡眠状態モニタリングは、方法として、あるいは、圧力検出部に接続されたコンピュータにおいてコンピュータプログラムとして実施することもできる。【０００８】
本明細書において、睡眠状態とは、覚醒、ＲＥＭ睡眠、ステージ１（Ｓ１）〜ステージ４（Ｓ４）という覚醒と睡眠の全ての状態を指すために用い、睡眠ステージとは、睡眠状態から覚醒を除いたもの、すなわち、ＲＥＭ睡眠、Ｓ１〜Ｓ４を指すために用いる。また、拍動とは、心臓の拍動およびそれに影響されて生じる振動、例えば、動脈の脈動による身体表面の振動を含む。【０００９】
本発明においては、使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、コンピュータネットワークを通じてデータを送信するデータ送信部とを備えたネットワークに接続されたモニタリング部と、該コンピュータネットワークを通じて該装置から少なくとも拍動間隔を含む時系列データを複数回受信し、受信した時系列データを使用者に対応させて該時系列データを記録するための記憶部とを有する睡眠状態モニタリングシステムであって、前記コンピュータネットワークに接続されたいずれかのコンピュータが、複数回にわたって受信した前記時系列データを使用者に対応付けて時間に対するグラフとして表示し、過去のある期間における該使用者の拍動間隔のゆらぎの変化を提示するデータ表示部を備えている、睡眠状態モニタリングシステムが提供される。この睡眠状態モニタリングシステムは、モニタリング部にコンピュータネットワークを通じて接続可能にされているコンピュータにおいてプログラムとして実施することもできる。
【発明の効果】
【００１０】
本発明のいずれかの実施形態によれば、使用者を拘束せず睡眠状態にも影響を与えることなく、覚醒しているか睡眠しているか、あるいは、使用者が睡眠しているのがどの睡眠ステージであるかなどの睡眠状態を監視することが可能になる。また、本発明のいずれかの実施形態によれば、拍動間隔のゆらぎによって使用者の健康状態を監視することが可能になる。
【図面の簡単な説明】
【００１１】
【図１】本発明の睡眠状態モニタリング装置の実施形態を示すブロック図。
【図２】本発明の睡眠状態モニタリング方法の実施形態を示すフローチャート。
【図３】本発明の睡眠状態モニタリング方法の実施形態において、睡眠か覚醒かを判定する処理を示すフローチャート。
【図４】本発明の睡眠状態モニタリング方法の実施形態において、睡眠ステージを判定する処理を示すフローチャート。
【図５】本発明の睡眠状態モニタリング方法の実施形態において、スペクトルパターン法によって睡眠ステージを判定する処理を示すフローチャート。
【図６】本発明の睡眠状態モニタリング方法の実施形態において、自律神経成分比法によって睡眠ステージを判定する処理を示すフローチャート。
【図７】本発明の睡眠状態モニタリングシステムの実施形態を示すブロック図。
【図８】本発明の実施例においてＲ波を特定する処理を行った結果を示すコンピュータ画面の例。
【図９】本発明の実施例において被験者から得た体振動からＲＲ間隔のパワースペクトルを得た結果を示すコンピュータ画面の例であり、覚醒状態（図９（Ａ））、Ｓ１〜Ｓ４（図９（Ｂ）〜図９（Ｅ））およびＲＥＭ睡眠（図９（Ｆ））の場合のパワースペクトルを示す。
【図１０】本発明の実施例において被験者から得た体振動によって睡眠状態を判定した結果を示すコンピュータ画面の例。
【図１１】本発明の実施例において被験者から得た体振動によって睡眠状態を自動判定した結果と、従来の手法によって判定した結果を対比して示すコンピュータ画面の例。
【符号の説明】
【００１２】
１  睡眠状態モニタリング装置
２  モニタリング部
１０  圧力波形取得部
１２  空気圧センサー
１４  トランスデューサ部
１６  Ａ／Ｄ変換機
２０、５０  コンピュータ
２４  拍動間隔算出部
２６  睡眠ステージ判定部
３２  呼吸波形抽出部
３４  呼吸間隔算出部
３６  覚醒判定部
４２  睡眠状態判定部
４６  データ送信部
６  睡眠状態モニタリングシステム
６０  コンピュータネットワーク
７０  記憶部
８０  データ生成部
【発明を実施するための最良の形態】
【００１３】
１．睡眠状態の判定
以下図面を参照して、本発明の実施の形態について説明する。本発明の実施の形態においては、体振動波形から睡眠状態を判定する。体振動波形の取得は、図１に示したように、圧力波形取得部１０とコンピュータ２０とを備える睡眠状態モニタリング装置１によって行う。圧力波形取得部１０は、空気圧センサー１２、圧力を電気信号に変換するトランスデューサ部１４、および、その電気信号をアナログ／デジタル変換するＡ／Ｄコンバータ１６を含んでおり、本発明においては、これらを集合的に圧力波形取得部１０と呼ぶ。空気圧センサー１２は、例えばベッド１００に配置されていて、その上に使用者が横たわることができるように作られている。例えば、気密性の高い袋状のエンベロープにスポンジ状の詰め物が封入されていて、全体として１ｃｍ〜３ｃｍ程度の厚み、３０ｃｍ四方の正方形形状となるクッション状になっていて、使用者がその上に身体を横たえても睡眠が可能なような柔軟性や肌触りとなるように作られている。そして、そのエンベロープの周縁部の一箇所には、内部の空気と連通して圧力を伝播させることができるチューブが接続されていて、そのチューブがトランスデューサ部１４に延びている。トランスデューサ部１４には、そのチューブが接続される気密室があり、その気密室の一つの側面には開口部が設けられて、その開口部を塞ぐようにセラミックス製の圧電素子が配置されている。その圧電素子には、空気圧センサー１２の内部の圧力波が伝達してアナログ値の圧電信号が生成され、その圧電信号は圧力波信号となって適当な増幅回路によって増幅される。増幅された圧力波信号は、Ａ／Ｄ変換機１６に入力され、デジタル化された圧力波信号が出力される。【００１４】
デジタル化された圧力波信号は、コンピュータ２０に入力される。コンピュータ２０は、プログラムされた命令に従って、所定の処理を行う機能ブロックの集合として動作する。すなわち、コンピュータ２０には、拍動抽出部２２と、拍動間隔算出部２４と、睡眠ステージ判定部２６と、呼吸波形抽出部３２と、呼吸間隔算出部３４と、覚醒判定部３６と、睡眠状態判定部４２とが備えられる。その結果は、例えばモニター４４に出力して視覚的に示すなどの手法によって使用者の睡眠状態として提示される。このように、圧力波形取得部と、コンピュータ２０に実現される各機能ブロックによって、本発明の実施形態の睡眠状態モニタリング装置が構成される。なお、当業者には明らかであるように、コンピュータ２０に備えられる適当な記憶手段によってデータが蓄積されたり、入出力手段によって処理の開始や休止、終了、中断などが制御されたりすることができる。【００１５】
次に、本発明の睡眠状態モニタリング装置の実施形態における処理のフローの概要について説明する。本発明の睡眠状態モニタリング装置は、睡眠か覚醒かの判定と睡眠ステージの判定とからなる睡眠状態の判定を行うことができる。具体的には、図２に示したように、睡眠状態の判定のためには、まず、体振動データを取得する（Ｓ１０）。そして、睡眠か覚醒かの判定する（Ｓ２０）。睡眠である場合には（Ｓ３０、ＹＥＳ）、睡眠ステージの判定をさらに行う（Ｓ４０）。その後、あるいは、睡眠ではない場合（覚醒である場合、Ｓ３０、ＮＯ）には、判定結果を出力する。このように、本発明の実施形態における体振動を用いた睡眠状態の判定の処理は、覚醒しているか睡眠しているかを判定する処理と、睡眠ステージを判定する処理とに分かれる。図示していないが、モニタリングの目的としては、この一連の処理を一定期間ごと（例えば１分ごと）に繰り返して、データの取得から睡眠状態の判定までを定期的に行ったり、あるいは、データを取得する処理のみを定期的に行い、適当な記憶手段にデータを格納しておいて、判定の処理を事後的に行ったりすることができる。【００１６】
以下、本発明の実施形態において実施される睡眠ステージの判定手法と睡眠状態の判定手法について説明する。図３には、睡眠か覚醒かの判定するステップ（図２、Ｓ２０）の詳細を示し、また、図４〜６には、睡眠ステージの判定のステップ（Ｓ４０）の詳細を示す。【００１７】
１．１．睡眠と覚醒との判定（Ｓ２０）
使用者の睡眠状態のうち睡眠か覚醒かを判定するために、本発明は、体振動波形から呼吸波形を抽出して利用する形態により実施することができる。呼吸波形は、人の安静時の呼吸数からおのずとその範囲が定まり、例えば体振動のうちの０．１３〜０．７０Ｈｚの周波数領域成分を呼吸波形とすることができる。このため、体振動波形をフィルタ処理し、０．１３〜０．７０Ｈｚの周波数領域成分を抽出する（図３、Ｓ２０２）。呼吸波形からは、呼吸間隔のゆらぎを抽出する。このために、呼吸波形が極大値を示す時刻を算出してその時刻を呼吸のタイミングを代表する時刻（呼吸ピーク時刻）とし（Ｓ２０４）、各呼吸の呼吸ピーク時刻データからから直前の呼吸の呼吸ピーク時刻データを減算することにより、呼吸間隔のデータを生成する（Ｓ２０６）。そして、呼吸間隔についての変動係数（ＣＶ）を求めて（Ｓ２０８）、これを呼吸間隔の一様性についての指標とする。というのは、使用者が睡眠中である場合には呼吸間隔は一様性が高いため変動係数の値が小さくなるのに対し、使用者が覚醒中である場合には呼吸間隔は一様が低くなって変動係数の値が大きくなるためである。具体的には、呼吸間隔のＣＶが０．２程度を閾値として設定し（Ｓ２１０）、それより大きい場合には覚醒状態（Ｓ２１２）、小さい場合には睡眠状態と判定することができる（Ｓ２１４）。【００１８】
１．２．睡眠ステージの判定（Ｓ４０）
睡眠ステージを判定する処理は、心拍または脈動による拍動のＲ波によって拍動のタイミングを特定する段階（Ｓ４０２）と、ある心拍とその直前の心拍との間の時間間隔（ＲＲ間隔）のデータを得る段階（Ｓ４０４）と、ＲＲ間隔のパワースペクトルを算出する段階（Ｓ４０６）と、睡眠ステージを判定する段階（Ｓ４２）とに分けて実行される。【００１９】
１．２．１．Ｒ波のタイミングの特定（Ｓ４０２）
睡眠ステージを判定する処理は、まず、体振動から、心電図（ＥＣＧ）で言うところのＲ波に対応する波の部分を特定する。なお、Ｒ波の記載は、理解の容易のためにＥＣＧにおけるＲ波との表現を用いたものであり、特にＥＣＧのＲ波に直接対応するものだけに限定されるものではない。連なっている心拍や脈動によって身体の表面に伝わる圧力波において拍動同士の時間間隔を測定するために各拍動に対してその拍動を代表するタイミングとして定めるための特徴波形であれば任意である。このためには、体振動波形から拍動に当たる波を抽出する必要がある（図１、拍動抽出部２２）。この抽出のためには、まず、体振動波形をバンドパスフィルターによって０．１Ｈｚ〜３０Ｈｚの範囲に帯域制限する。次に、フィルター後の体振動波形を時間に関して二次微分する処理を行って、体振動の加速度の各時間のデータを算出する。そして、その体振動の加速度の時間毎のデータを、さらにバンドパスフィルターによって３〜３０Ｈｚの範囲に帯域制限して、その波形データがピーク値を示す時刻を特定する。拍動抽出部２２はそのための微分ピーク処理部２２Ａを備えることができる。さらに、１秒未満の周期の体振動波形を用いて、波形データが急激に変化する区間に限定し、４秒程度の体振動波形の周期性とＲ波の周期性に注目して、取り出したピーク値の時刻がほぼ等間隔となるように補正して決定する。この時刻は、３〜３０Ｈｚの範囲で体振動の加速度が大きいと評価される時刻となるため、心臓の拍動のうち心室の興奮による振動つまり心電図におけるＲ波の時刻に相当するものと本願発明者は推測している。しかしながら、体振動から振動を抽出するという振動検出の性質上、その時刻には、心拍のみによる振動だけではなく、心臓による動脈による振動（脈動）などの影響があり得る。【００２０】
１．２．２．ＲＲ間隔データの導出（Ｓ４０４）
このようにして得られたＲ波の時刻データからＲＲ間隔（ある拍動とその直前の拍動との間の時間間隔）の時系列データを求める（図４、Ｓ４０４；図１、拍動間隔算出部２４）。ＲＲ間隔は、各拍動のＲ波の時刻データから直前の拍動におけるＲ波の時刻データを減算することにより与える。このＲＲ間隔のデータは、拍動のＲ波を検出した不均等間隔の各時刻における時間データである。このデータは、後の処理のため、一定のサンプリング間隔の時系列データにするのが好都合であり、補間処理を行ってＲＲ間隔の時系列データを得る。これには、例えば２次補間などの処理を行うことができる。以上のようにして得られたＲＲ間隔の時系列データから次に睡眠の深度を判定する。【００２１】
１．２．３．ＲＲ間隔のパワースペクトルの算出（Ｓ４０６）
ＲＲ間隔の時系列データからパワースペクトルを求めて利用する（Ｓ４０６；図１、パワースペクトル算出部２６Ａ）。なお、パワースペクトルの算出は既知のいずれの方法も用いることができ、例えば、高速フーリエ変換アルゴリズムを用いたり、Wiener-Khinchinの定理による自己相関に基づく手法などを用いることができる。パワースペクトルを算出するための窓関数（ウィンドウ関数）は任意であるが、判定の精度を高めるために２分間のデータを利用する。この点については後述する。【００２２】
１．２．４．睡眠ステージの判定
次いで、睡眠ステージの判定行う（図４、Ｓ４２；図１、睡眠ステージ判定部２６）。睡眠ステージの判定は、睡眠している状態が、ＲＥＭ睡眠およびＮｏｎ−ＲＥＭ睡眠（Ｓ１〜Ｓ４）のいずれであるかを判定する処理である。睡眠ステージの判定について、本発明は２つの実施の形態によって実施することができる。一つはＲＲ間隔の時系列データのパワースペクトルのパターンから睡眠ステージを判定する実施形態（スペクトルパターン法）であり、もう一つは、ＲＲ間隔の時系列データから自律神経機能の働きを示す指標を算出して睡眠ステージの判定を行うもの（自律神経成分比法）である。【００２３】
１．２．４．１．睡眠ステージの判定：スペクトルパターン法
従来の睡眠ステージの判定方法によって判定しておいた睡眠ステージとＲＲ間隔のパワースペクトルとを対応付けておいて、各睡眠ステージについての典型的なスペクトルのパターン（「典型パターン」）を得ることができる。そして、ある被験者に対して測定されたＲＲ間隔のパワースペクトルのパターン（以下、「未知パターン」）がいずれの典型パターンに類似するかを判定することによって、睡眠ステージを推定することができる。【００２４】
図５を参照しながらその詳細を説明する。未知パターンがいずれの典型パターンに類似するかの判定は、ＲＲ間隔のパワースペクトルのパターンのうち、ＬＦ範囲（Ｌｏｗ  Ｆｒｅｑｕｅｎｃｙ；０．０４〜０．１５Ｈｚ）およびＨＦ範囲（Ｈｉｇｈ  Ｆｒｅｑｕｅｎｃｙ；０．１５〜０．４０Ｈｚ）においてパワーの最大値を求めて（図５、Ｓ４２２）、そのパワーの値のｎ分の１の値（ｎは２以上の整数）を求めて、その値における周波数幅を評価対象として行う。具体的には、まず、ｎを２にして（Ｓ４２４）、ＲＲ間隔のパワースペクトルの最大値の１／ｎの値を求める（Ｓ４２６）。そして、パワースペクトルに極大値を与える周波数の両側において、最大値の１／ｎとなる周波数を求めてその周波数同士の差（周波数幅）を計算する（Ｓ４２８）。この処理を、ｎを逐次増加（インクリメント）して（Ｓ４３２）所定の上限Ｎまで繰り返す（Ｓ４３０）。各ｎに対してそのときの周波数幅が与えられるので、周波数幅のｎ依存性がえられる。典型パターンについて同様に求めておいた周波数幅のｎ依存性に対して、未知パターンの周波数幅を対比させて、典型パターンと未知パターンの周波数幅の差を残差として算出する（Ｓ４３４）。こうして、各睡眠ステージの典型パターンに対する未知パターンの一致度を数値として得ることができる。ここで、未知パターンのＬＦ範囲のピークは典型パターンのＬＦ範囲のパターンのピークになるべく一致するようにマッチングさせ、同様に未知パターンのＨＦ範囲のピークは典型パターンのＨＦ範囲のパターンのピークとマッチングさせるようにして、異なる周波数範囲の間でのマッチングは行わない。そして、未知パターンに対して最も高い一致度（少ない残差）を示す典型パターンとなるような睡眠ステージを、その未知パターンの睡眠ステージとする（Ｓ４３６）。なお、ここに記載したＬＦ範囲とＨＦ範囲の周波数の範囲は、後述する自律神経成分比法における記載も含めて典型的なものであり、厳密に上述の周波数範囲にのみ限定されず、例えば、ＬＦ範囲を０．０５〜０．１５Ｈｚなどとすることも本発明の開示の一部とする。【００２５】
なお、スペクトルパターン法では、睡眠ステージの判定を行うことができるが、睡眠と覚醒のどちらにあるかの判定、すなわち、ＲＥＭ睡眠と覚醒との判定を行うことはできない。【００２６】
本願発明者は、スペクトルパターン法において判定精度を高めるために、ＲＲ間隔のパワースペクトルについて、上述のＬＦ範囲とＨＦ範囲に加えて、ＶＬＦ範囲（Ｖｅｒｙ  Ｌｏｗ  Ｆｒｅｑｕｅｎｃｙ；０．０１Ｈｚ〜０．０４Ｈｚ）のスペクトルを利用することが有用であることを見出している。ＶＬＦ範囲のスペクトルとは、ＲＲ間隔のパワースペクトルにおいてＬＦ範囲よりも低周波に現れるスペクトルである。ＶＬＦ範囲のスペクトルを正しく求めるためには、ＲＲ間隔の時系列データはそれに見合った期間だけのデータが必要となり、必要なスペクトル範囲の下限周波数の逆数をウインドウ期間とすることが必要となる。例えば、０．０１Ｈｚのパワースペクトルを正しく求めるためには、このウインドウ期間を１００秒程度の適当な時間にする必要があり、実際に測定する場合には、例えば２分間以上に設定する。【００２７】
スペクトルパターン法の実施形態として典型パターンを用いた一致判定を行う手法を上述したが、本発明は任意のパターン判定手法によって実施することができる。例えば、スペクトルのパターンに対してそれが判定されるべき結果（睡眠ステージのいずれか）を対応させる学習処理をニューラルネットワークに対して行っておいて、未知パターンに対してそのニューラルネットワークを作用させて判定する、いわゆる教師付ニューラルネットワーク学習システムを用いることができる。その他の任意の教師付学習の手法を実装することもできる。【００２８】
１．２．４．２．睡眠ステージの判定：自律神経成分比法
自律神経機能の働きを示す自律神経成分比によっても、睡眠ステージの判定は可能である。このためには、まず、ＲＲ間隔のパワースペクトルにおいて、ＬＦ成分とＨＦ成分とを求める。そして、ＬＦ成分とＨＦ成分の比、すなわちＬＦ／ＨＦ比（自律神経成分比）を算出する。これは、被験者の自律神経機能の働きがＬＦ／ＨＦ比により評価できるためである（非特許文献２〜４参照）。この解析には、図６に示したように、ＬＦ範囲のパワースペクトルを算出してその範囲の積分値を得てＬＦ成分とし（Ｓ４４２）、ＨＦ範囲のパワースペクトルを算出してその範囲の積分値を得てＨＦ成分とする（Ｓ４４４）。そして、ＬＦ成分をＨＦ成分で除算して得られるＬＦ／ＨＦ比を求める（Ｓ４４６）。各睡眠ステージにおける典型的なＬＦ／ＨＦ比とこのＬＦ／ＨＦ比を比較して、睡眠ステージを判定する（Ｓ４４８）。ＬＦ／ＨＦ比は、小さいほど睡眠が深く大きいほど睡眠が浅いが、睡眠と覚醒のどちらにあるかの判定、すなわち、ＲＥＭ睡眠と覚醒との判定を行うことはできない。また、ＲＥＭ睡眠、Ｓ１〜Ｓ４の各睡眠ステージとＬＦ／ＨＦ比との関係には個人差が大きいので、例えばスペクトルパターン法による睡眠ステージの判定結果との対応関係を予め測定しておいて、使用者ごとに適応させることも有効である。【００２９】
以上の実施の形態により、使用者の睡眠状態に影響を与えにくく、使用者が在宅で利用することができ、また、睡眠か覚醒かを判定できるような睡眠状態モニタリングシステムまたは方法が実現される。【００３０】
２．睡眠時の拍動間隔のモニタリング
次に、睡眠時の拍動のモニタリングの手法について説明する。図７にこのモニタリング処理を行うための構成を示す。本発明では、コンピュータ５０と圧力波形取得部１０を集合的にモニタリング部２と呼ぶ。図１と同様に、体振動波形の取得には、圧力波形取得部１０、すなわち、空気圧センサー１２、圧力を電気信号に変換するトランスデューサ部１４、および、その電気信号をアナログ／デジタル変換するＡ／Ｄコンバータ１６を利用する。空気圧センサー１２がベッド１００に配置されていることができるのも同様である。また、図１と同様にトランスデューサ部１４およびＡ／Ｄ変換機１６が接続されて、デジタル化された圧力波信号が出力される。【００３１】
デジタル化された圧力波信号は、コンピュータ５０に入力される。コンピュータ５０は、プログラムされた命令に従って、所定の処理を行う機能ブロックの集合として動作する。すなわち、コンピュータ５０には、拍動抽出部２２と、拍動間隔算出部２４と、データ送信部４６とが備えられる。【００３２】
データ送信部４６は、例えばインターネットやイントラネットなどのコンピュータネットワーク６０を通じて記憶部７０にデータをアップロードすることができる。記憶部７０には、使用者のＩＤと対応付けたデータが複数回アップロードされており、コンピュータネットワーク６０に接続されたいずれかのコンピュータにおいて、複数回にわたって受信したデータをつなげることによって複数回のアップロードされたデータを一つのグラフに表示することができる。このとき、使用者ごとにグラフに表示するために、同じ使用者のデータを検索して連結させるデータ生成部８０が用いられる。そして、例えばモニター８６に表示することができる。【００３３】
アップロードされるデータは、拍動間隔のデータを含む任意の時系列データである。例えば、拍動間隔の時系列データそれのみであってもよいし、図１の睡眠ステージ判定部２６と同様の睡眠ステージ判定部（図示しない）が備えられている場合には、拍動間隔のデータに加えて睡眠ステージの時系列データとすることもできる。図１の呼吸波形抽出部３２と、呼吸間隔算出部３４と、覚醒判定部３６と、睡眠状態判定部４２とに相当する機能ブロック（図示しない）を備えていれば、拍動間隔のデータに加えて、覚醒か睡眠の判定結果が含まれている時系列データとすることもできる。【００３４】
データ生成部８０は、モニター８６に表示するためのデータを生成するほか、拍動間隔のデータに見られるゆらぎの強さに、通常とは異なる明らかな変化が生じるのを検知してアラーム信号を出力するゆらぎ検知部８２を備えることができる。なお、ゆらぎ検知部８２は、データ生成部８０の一部としているが、コンピュータネットワーク６０に接続されて使用者ごとのデータを記憶部７０から呼び出すことができれば、コンピュータネットワーク６０に接続された任意のコンピュータ機器に備えられていることができる。【００３５】
ゆらぎ検知部８２は、目的に応じて長さが定まるある期間（例えば１年間にわたるデータを記憶部７０に蓄積する場合には、一ヶ月などの任意の所定期間）における拍動間隔の変化を監視している。この変化を検出するために、拍動間隔は、例えば拍動間隔データのもつ標準偏差の値、変動係数の値などの数値となるように整理され、その値がある閾値を下回るかどうかという点で監視される。そして、例えば、その閾値を下回ったときにアラーム信号が送信されて、適当な手段によって使用者にフィードバックされる。このためには、例えば自動送信される電子メールや、パーソナライズされた管理用Ｗｅｂページにアラートを表示するなどの手法を用いることができる。【００３６】
これにより、在宅で測定される拍動間隔に現れる体調変化をいち早くキャッチして、睡眠時の拍動間隔として現れる使用者の体調の変化を知ることができ、体調管理を行うことができる。
【実施例１】
【００３７】
３．睡眠ステージの判定における従来の手法との相関の確認
睡眠ステージの判定については、従来から確立されている手法と本発明の実施形態とを比較するために、被験者に、従来の手法による睡眠ステージの判定と本発明の実施の形態による睡眠の判定を行った実験結果について説明する。２０歳代の男女２５名を被験者とし、被験者に睡眠（昼眠）を取らせ、そのうちの９０分間を測定に利用した。【００３８】
３．１. 睡眠ステージの判定における従来の手法
被験者から脳波、アゴ筋電図、眼球電位図、心電図を採取して従来の手法による睡眠ステージの判定を行った。各測定値を得るために、被験者の身体の各部には皿電極をつけた。具体的には、脳波のためには頭頂に、アゴ筋電図のためにアゴに、眼球電位図のために眼の左右周辺に、そして、心電図のために足首接地により両手首に、それぞれ電極を装着した。皿電極には低周波増幅器を接続して電圧信号を増幅し、ＡＤ変換器によってディジタル値に変換したのち、データログ用のコンピュータによって記録をとった。脳波によって睡眠ステージを判定し、その際、前述のＲｅｃｈｔｓｃｈａｆｆｅｎ−Ｋａｌｅｓ（Ｒ−Ｋ）の方法に従った。【００３９】
３．２．本発明の実施の形態の睡眠ステージの判定の手法
心拍動波形からＲ波を抽出し自律神経成分比を算出すると共に、呼吸波形も取得した。【００４０】
本実施の形態の手法のための測定および解析のための装置としては、図１に記載の装置である圧力波形取得部１０およびコンピュータ２０を用い、そのコンピュータでは図１のような各機能ブロックとしてコンピュータを動作させ、図４〜６のような動作を行わせるために、体振動から心拍動成分と呼吸成分を抽出するソフトウエアモジュールと、心拍動波形からＲ波を抽出してスペクトルパターン法の解析を行うソフトウエアモジュールと、自律神経成分比法を実施するためのソフトウエアモジュールとを利用した。【００４１】
３．３．Ｒ波のタイミングの相関
図８に、体振動からＲ波のタイミングを特定した例を示す。図８の上段のチャートＣ８Ａには、体振動波形Ｃ８Ａ−１０と、それから算出した呼吸波形Ｃ８Ａ−２０も示している。中心線Ｃ８Ａ−００の上部に示された波形Ｃ８Ａ−２は、体振動波形Ｃ８Ａ−１０からＲ波の候補となる時刻を検出したものである。この検出は、拍動抽出部２２によって行った。このとき、波形Ｃ８Ａ−２のピークＰ１〜Ｐ８のうち、ピークＰ２´は誤検出されたものである。このような誤検出が生じたのはノイズなどの影響であるため、誤検出のピークを除去する処理を行う。この処理は、Ｒ波のある候補が得られた時点と次の候補が得られた時点が、予め設定した時間（例えば、０．５秒）以下になった場合には、その次の候補は採用しない処理などであり、拍動によらずに得られた可能性の高い波形をＲ波として採用しないようにする処理である。この処理を行って候補から誤検出のピークを除去し、中心線Ｃ８Ａ−００の下部に示したようにＲ波の時刻を特定している。なお、この図は８秒間の現象を拡大して表示しているものである。【００４２】
そして、従来の手法によって心電図から得たＲ波と、本発明の実施の形態によって得た体振動のデータから特定したＲ波との間で発生時点を比較した。その結果、空気圧センサーにより測定された体振動波形から求めたＲ波の発生時点の検出精度は、心電図から得たＲ波発生時点と比較して、両者の一致度は平均８７％であった。また、このようにして算出したＲ波の時系列データ（０．１秒間隔の２次補間したもの）について、体振動から判定したＲＲ間隔と心電図によるＲＲ間隔との間の相関係数は、０．８０〜０．９５であり、平均相関係数は、０．８７であった。【００４３】
なお、図８の中段のチャートＣ８Ｂに現れる波形Ｃ８Ｂ−１０は呼吸波形であり、横軸は時間で１１分間の波形であり、縦軸は圧力波形（任意単位）である。また、チャートＣ８Ｂの波形Ｃ８Ｂ−２０は推定したＲＲ間隔である。【００４４】
図８の下段のチャートＣ８Ｃは、ＲＲ間隔から上記実施形態の方法によって求めたパワースペクトルＣ８Ｃ−１０である。ここに示したパワースペクトルＣ８Ｃ−１０は、深い睡眠ステージにおける結果で、グラフに表示された周波数範囲は０．０１〜０．４０Ｈｚの範囲のものであり、パワースペクトルには高周波（ＨＦ）範囲の０．２５〜０．３０Ｈｚの位置に顕著なピークＣ８Ｃ−Ｐ１が現れた。【００４５】
３．４．睡眠ステージ判定
そして、体振動から本実施の形態の手法によって睡眠ステージを判定した。睡眠ステージの判定のためには、前述のように２分間のデータを用いた。【００４６】
３．４．１.スペクトルパターン法
具体的には、図８に関連して説明したような手法によって得られたパワースペクトルを未知の睡眠状態にある被験者から得て、そこにスペクトルパターン法を適用して睡眠ステージを判定した。まず、スペクトルパターン法の判断対象になるスペクトルの観察結果について説明する。【００４７】
図９の（Ａ）〜（Ｅ）に、覚醒、Ｓ１〜Ｓ４のそれぞれの睡眠状態の場合のパワースペクトルを示す。また、図９の（Ｆ）にＲＥＭ睡眠のパワースペクトルを示す。これらのスペクトルにおいて、グラフ上は０．００〜０．０１Ｈｚおよび０．４０〜０．５０Ｈｚの値も表示されているものの、利用する周波数範囲は０．０１Ｈｚ〜０．４０Ｈｚである。これらの図に示すように、ＲＲ間隔のパワースペクトルは、図９の（Ａ）の覚醒時ではＲＲ間隔の変動が不規則なものである。その結果、ＲＲ間隔のパワースペクトルはＬＦ、ＨＦの範囲に小さな成分しか有さず、０．０４Ｈｚ以下のＶＬＦ範囲が主要なスペクトルとなった。【００４８】
また、図９の（Ｂ）および（Ｃ）に示す浅い睡眠ステージＳ１とＳ２においては、ＲＲ間隔の変動の波形（図示しない）にリズミカルな変化が見られ、パワースペクトルにもこれを反映するような明確なピークがＬＦ範囲のスペクトルとＨＦ範囲のスペクトルに現れた。このＬＦ範囲のスペクトルは、睡眠ステージがＳ１からＳ２へ移行するとより大きくなった。【００４９】
睡眠ステージＳ１は覚醒段階からＳ２との段階の移行状態にあり、覚醒時の特徴に近い分布すなわちＬＦとＨＦに亘る小さいパワースペクトル値の分布が見られた。睡眠ステージＳ１に覚醒状態のパワースペクトルが混入してくる場合もあるため、覚醒状態とＳ１の識別を実際のパワースペクトルのみから行うことは難しい。これに対し、図９の（Ｃ）に示す睡眠ステージＳ２は、ＬＦ範囲のスペクトルが存在するもののＨＦ範囲のスペクトルがＳ１よりも明確に出現してくるので、Ｓ１とＳ２の識別はＳ１と覚醒状態の識別よりも容易である。【００５０】
図９の（Ｄ）およびＥに示す深い睡眠ステージであるＳ３とＳ４においては、睡眠が深くなるにつれてＨＦ範囲のパワースペクトルがＬＦ範囲のパワースペクトルより大きくなり、Ｓ４においてＨＦ範囲のピークが主要なパワースペクトルとなる。従って、Ｓ４とＳ３の識別は、Ｓ１とＳ２との識別と同様に容易である。【００５１】
そして、図９の（Ｆ）はＲＥＭ睡眠のパワースペクトルパターンである。睡眠が覚醒に近づいたときに現れるＲＥＭ睡眠では、覚醒及びＳ１とパターンが類似しているが、ＬＦ範囲のスペクトルが、Ｓ１の場合よりもより大きく複雑になるという特徴があった。【００５２】
本発明者の検討によれば、ＲＥＭ睡眠と覚醒の区別は、ＲＲ間隔の時系列データから変動係数（ＣＶ）を求めることによって区別することができた。また、覚醒の推定には、呼吸のＣＶも考慮した。この手法は、入眠時、中途覚醒、睡眠終了後の覚醒の判定に有効であった。以上、観察者がＲＲ間隔のパワースペクトルを直接観察した場合の特徴について述べた。【００５３】
そして、以上のパターンの特徴から睡眠ステージの判定を再現性が高いものとするために、睡眠ステージの識別のために上述のスペクトルパターン法を適用した。その結果を図１０と図１１に示す。【００５４】
図１０は空気圧センサーから得た体振動波形Ｃ１０Ａ−１０および呼吸波形Ｃ１０Ａ−２０（最上段のチャートＣ１０Ａ）、検出したＲ波時系列Ｃ１０Ｂ−１０（上から２段目のチャートＣ１０Ｂ）、抽出した呼吸数Ｃ１０Ｃ−１０と心拍数Ｃ１０Ｃ−２０（上から３段目のチャートＣ１０Ｃ）、ＲＲ間隔のパワースペクトルＣ１０Ｄ−１０（最下段左のチャートＣ１０Ｄ）、および、睡眠ステージの判定例Ｃ１０Ｅ−１０およびＬＦ／ＨＦ比Ｃ１０Ｅ−２０（最下段右のチャートＣ１０Ｅ）である。また、図１１は、睡眠状態の判定結果を比較したものであり、折れ線Ｒ０は脳波を用いて従来の方法によって判定した結果であり、折れ線Ｒ１（破線）は、スペクトルパターン法によって自動判定した結果である。この自動判定結果の折れ線Ｒ１は、図１０のチャートＣ１０Ｅに示したものと同じ結果であり、縦軸は、１から−４にかけて、順に、覚醒、ＲＥＭ睡眠、Ｓ１〜Ｓ４を表わしており、７６分の体振動の測定データから得られたものである。従来の手法と実施例の手法との間での一致度は約７５％となった。【００５５】
３．４．２.自律神経成分比法（ＬＦ／ＨＦ比を用いる方法）
次に、自律神経成分比法の実施例について説明する。上述の条件により取得した体振動波形から、浅い睡眠ステージ（Ｓ１とＳ２）と深い睡眠ステージ（Ｓ３とＳ４）との区別を自律神経成分比法により行った。このとき、浅い睡眠と深い睡眠の判定、ＬＦ／ＨＦ比が０．３より大きいときに浅い睡眠とし、０．３以下のときに深い睡眠と判定した。取得したデータの範囲では、この判定結果と、脳波による睡眠ステージの判定結果との一致度は共に８０％以上となった。なお、このＬＦ／ＨＦ比については、個人差が大きいため上記の判定基準は一例であり、例えば、判定の用いる使用者ごとの値は、参照可能なデータベースに別途記録しておくことができる。【００５６】
３．５．睡眠／覚醒の判定
以上のスペクトルパターン法および自律神経成分比法では覚醒時とＲＥＭ睡眠の区別は難しいが、呼吸ピークの変動係数（ＣＶ）やＲＲ間隔の変動係数（ＣＶ）を用いることで、両睡眠状態が区別可能となった。本発明の実施の形態の睡眠と覚醒との判定について、呼吸の変動係数（ＣＶ）によって判定を行った。その結果、呼吸のＣＶの閾値を０．２として、それ以上の値を覚醒と判定することが有効であった。また、ＲＲ間隔の変動係数（ＣＶ）には個人差が生じたが、特定の個人では十分に高い再現性が実現できた。【００５７】
４．拍動モニタリング
次に、図７に示した装置を用いて拍動間隔（ＲＲ間隔）をモニタリングする睡眠状態モニタリングシステム６の作製例を説明する。モニタリング部２は、Ｒ波のデータに対するＲＲ間隔の算出処理を行う。そしてそのＲＲ間隔のデータを、記憶部７０を管理するデータサーバに一日一度アップロードする。アップロードは任意のデータ通信手段を用いることができるが、例えば、ＦＴＰ（ファイルトランスファープロトコル）によってプットすることにより行うことができる。記憶部７０におけるデータの蓄積の態様は任意であり、例えば、データサーバのファイルシステムのフォルダ（ディレクトリ）を利用して使用者のＩＤごとおよび日付ごとのＣＳＶ（Ｃｏｍｍａ  Ｓｅｐａｒａｔｅｄ  Ｖａｌｕｅｓ）形式のコンピュータファイルを複数蓄積することができる。別の例としては、データベースマネージメントシステムなどを利用して適当なデータベースに蓄積することができる。なお、睡眠開始（入眠）のタイミングを適切に判定するため、コンピュータ５０には図示しない光センサーを接続しておいて、周囲の照度情報によって測定開始のタイミングを制御することができる。さらに、測定の時間を入眠後の睡眠サイクル２回分に限定するために、測定開始後３時間のデータのみ取得することができる。【００５８】
このＣＳＶ形式のファイルに含まれるＲＲ間隔のデータから、ＲＲ間隔のゆらぎを算出する。このゆらぎは変動係数（ＣＶ）の算出によって行うことができる。複数のＣＳＶファイルに分かれている各日のデータは連結されて計算され、日付を追ってＲＲ間隔のゆらぎを表示することができる。【００５９】
また、ゆらぎ検知部８２には、例えば、過去一ヶ月間のＲＲ間隔のゆらぎの変動を監視する短期変動監視モジュール（図示しない）と、過去１年間のＲＲ間隔のゆらぎの変動を監視する長期変動監視モジュール（図示しない）を実装することができる。ここで、短期変動監視モジュールおよび長期変動監視モジュールは、ＲＲ間隔のゆらぎに現れる変動がどのように推移するかをそれぞれの期間に合わせて監視している。これらのモジュールでは、使用者の健康状態を監視するために、睡眠の時間の中から監視に適した睡眠状態にある時間において取得されたＲＲ間隔データを抽出することができる。このような監視に適した状態にある時間は、例えば、ＲＲ間隔に長周期のゆらぎが生じていない時間を選択したり、あるいは、睡眠ステージがＳ３あるいはＳ４にある時間を選択することができる。このような目的で、図７に示した睡眠状態モニタリングシステム６には、図１に示した睡眠ステージ判定部２６、呼吸波形抽出部３２、呼吸間隔算出部３４、覚醒判定部３６、睡眠状態判定部４２などの機能ブロックをさらに追加しておいて、短期変動監視モジュールおよび長期変動監視モジュールがそれらの機能ブロックからの出力によって動作が制御されることができる。これらにより、図７に示した睡眠状態モニタリングシステム６においては、ＲＲ間隔のゆらぎがどのように変動するかを、在宅の使用者に特段の負担をかけることなく、継続して監視することができる。【００６０】
５．まとめ
以上、本発明の実施例において、空気圧センサーから心拍変動を抽出して、Ｒ波を検出し、心電図との比較において、高い精度でＲ波を特定することができた。また、呼吸ピークの変動係数（ＣＶ）やＲ波の変動係数によって睡眠と覚醒とを適切に区別することができた。さらに、特定したＲ波からＲＲ間隔を求めてＲＲ間隔のパワースペクトルを算出し、スペクトルパターン法及び自律神経成分比法によって睡眠ステージの判定が精度よく実行できることを示した。【００６１】
睡眠／覚醒の判定においてＲ波の変動係数を用いたとき、あるいは、睡眠ステージの判定に自律神経成分比法を用いたときには、個人差の問題が残っているものの、特定の被験者に対しては高い再現性が得られた。【００６２】
さらに、ＲＲ間隔のゆらぎがどのように変動するかを、在宅の使用者に特段の負担をかけることなく、継続して監視することができる。【００６３】
以上、本発明の実施の形態につき述べたが、本発明は既述の実施の形態に限定されるものではなく、本発明の技術的思想に基づいて各種の変形、変更および組み合わせが可能である。たとえば、パターンの判定において、適当な人工知能システムを用いて、ＲＲ間隔のパワースペクトルパターンと、同時に従来の手法によって行った判定結果との関係を学習させて、パワースペクトルから睡眠ステージの判定結果を得ることができる。また、個人差が避けられない判定（例えばＬＦ／ＨＦ比による判定）については、被験者別の判定基準を決めるようなファイルを利用して、被験者ごとに適応させた判定を実行することもできる。【００６４】
また、図７に示したシステムによって、拍動の間隔における変化を監視して使用者や医療従事者に伝達することができる。使用者は、その情報に基づいて医療機関にて検査を受けるなどの対応を取ることができる。
【図１】
2009150744.tif 000003
【図２】
2009150744.tif 000004
【図３】
2009150744.tif 000005
【図４】
2009150744.tif 000006
【図５】
2009150744.tif 000007
【図６】
2009150744.tif 000008
【図７】
2009150744.tif 000009
【図８】
2009150744.tif 000010
【図９】
2009150744.tif 000011
【図１０】
2009150744.tif 000012
【図１１】
2009150744.tif 000013
【国際調査報告】
2009150744.tif 000014
2009150744.tif 000015
2009150744.tif 000016
2009150744.tif 000017
2009150744.tif 000018
2009150744.tif 000019
2009150744.tif 000020