JP2001061820A - Sleeping state judging method and sleeping state judging device - Google Patents

Sleeping state judging method and sleeping state judging device

Info

Publication number
JP2001061820A
JP2001061820A JP24053299A JP24053299A JP2001061820A JP 2001061820 A JP2001061820 A JP 2001061820A JP 24053299 A JP24053299 A JP 24053299A JP 24053299 A JP24053299 A JP 24053299A JP 2001061820 A JP2001061820 A JP 2001061820A
Authority
JP
Japan
Prior art keywords
sleep state
time
trend curve
biological information
unit
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
JP24053299A
Other languages
Japanese (ja)
Other versions
JP3788126B2 (en
Inventor
Hiroyuki Ibe
浩行 井邊
Kaoru Fukuyori
薫 福頼
Norio Nakano
紀夫 中野
Akihiro Michimori
章弘 道盛
Hiroshi Hagiwara
啓 萩原
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.)
Panasonic Electric Works Co Ltd
Original Assignee
Matsushita Electric Works Ltd
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 Matsushita Electric Works Ltd filed Critical Matsushita Electric Works Ltd
Priority to JP24053299A priority Critical patent/JP3788126B2/en
Publication of JP2001061820A publication Critical patent/JP2001061820A/en
Application granted granted Critical
Publication of JP3788126B2 publication Critical patent/JP3788126B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

PROBLEM TO BE SOLVED: To accurately detect a sleeping state by calculating a trend curve for showing an increase/decrease tendency of a time change in time series of living body information such as the heartbeat number, and calculating the trend curve by using a window function of preset time when judging a sleeping state on the basis of the curve. SOLUTION: A measuring part 1 is arranged to measure a living body information value such as the heartbeat number per unit time as a living body information measured value, and a trend curve for showing a tendency of a time change in time series of the living body information value is calculated by a trend curve calculating part 3. An increase in the living body information value in each measured time is calculated by an increase calculating part 4 from the living body information value and the trend curve to calculate a fluctuation degree of the living body information value in a specific time before each measured time by a fluctuation degree calculating part 5. An appearance probability of respective sleeping states is calculated by a calculating part from a combination of the increase and the fluctuation degree of the living body information value, and a state of adopting a maximum numeric value among the appearance probability of the respective sleeping states is estimated as a sleeping state of the corresponding measured time by an estimating part 8.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、心拍数や脈拍数の
ような比較的容易に得られる生体の活動情報に基づい
て、生体の睡眠状態の変化を検出するようにした睡眠状
態判定方法及び睡眠状態判定装置に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a sleep state judging method for detecting a change in the sleep state of a living body on the basis of relatively easily obtained information on the activity of the living body such as a heart rate and a pulse rate. The present invention relates to a sleep state determination device.

【0002】[0002]

【従来の技術】一般に、心拍信号はその人間の状態、例
えば、睡眠・覚醒、安静・運動、姿勢の変化、摂食、さ
らには情動や精神的なリラックスやストレス等によって
絶えず変化している。また、睡眠に関しても人間の睡眠
状態は一晩を通じて一様ではなく、ノンレム睡眠期とレ
ム睡眠期とのサイクルが周期的に数回出現し、その周期
が100分程度(80〜120分)であることが従来よ
り知られている。各サイクルでは、ノンレム睡眠期にお
いて浅い睡眠から深い睡眠へと睡眠状態が次第に移行
し、深い睡眠状態が持続した後、再び浅い睡眠状態とな
り、その後、レム睡眠期が出現するという変化のしかた
が一般的である。また、各サイクルにおけるノンレム睡
眠期の睡眠の深さの変化は相対的であり、入眠から覚醒
にかけてサイクル毎に睡眠深度が浅くなる。従って、上
述のような睡眠の状態変化を検出すれば、その状態に応
じて入眠を促進したり、より深い睡眠状態に誘うよう
に、また心地よく目覚めさせたりするように適宜の刺激
を与えることができる。
2. Description of the Related Art In general, a heartbeat signal is constantly changing due to the state of the human being, for example, sleep / wake, rest / exercise, change in posture, eating, emotion, mental relaxation and stress. Regarding sleep, the human sleep state is not uniform throughout the night, and a cycle between the non-REM sleep phase and the REM sleep phase appears periodically several times, and the cycle is about 100 minutes (80 to 120 minutes). It has been known for some time. In each cycle, the sleep state gradually shifts from light sleep to deep sleep in the non-REM sleep period, after the deep sleep state is maintained, the light sleep state is returned again, and then the change in the appearance of the REM sleep period is general It is a target. Further, the change in the sleep depth during the non-REM sleep period in each cycle is relative, and the sleep depth becomes shallower for each cycle from falling asleep to awakening. Therefore, if a change in the state of sleep as described above is detected, appropriate stimuli can be given to promote sleep onset according to the state, to induce a deeper sleep state, and to wake up comfortably. it can.

【0003】そこで睡眠の状態変化を検出することが従
来から種々試みられており、例えば脳波、眼球運動、筋
電などを含む睡眠ポリグラフを用いれば良い。しかしこ
のものは、装置が大がかりであり、研究室や病院などの
計測設備を備えた場所でしか利用できず、健康機器のよ
うに日常的に使用する用途には不向きである。そこで、
睡眠ポリグラフに代わる手段によって睡眠の状態変化を
精度良く検出することが望まれている。このような要求
に応えるために、心拍数や脈拍数に着目し、これらの変
化から人問の状態変化を検出する事が考えられている。
Therefore, various attempts have been made to detect changes in the state of sleep, and for example, a polysomnogram including brain waves, eye movements, and myoelectricity may be used. However, this device is large in scale and can be used only in a place equipped with measurement equipment such as a laboratory or a hospital, and is not suitable for daily use such as health equipment. Therefore,
It is desired to accurately detect a change in the state of sleep by means that is an alternative to polysomnography. In order to respond to such a demand, it has been considered that attention is paid to a heart rate and a pulse rate, and a change in the state of a human being is detected from these changes.

【0004】すなわち、夜間睡眠においては、単位時間
あたりの心拍数は入眠とともに減少し、朝方に最低値を
示し、覚醒時期が近付くにつれて増加することが知られ
ている。また、ノンレム睡眠期では単位時間当たりの心
拍数は比較的安定しているがレム睡眠期には自律神経の
活動状態に乱れが生じるから心拍数が著しく変動し、多
くの場合に増加傾向が見られるということが知られてい
る。さらに、ノンレム睡眠期の中でも睡眠深度が浅い場
合は単位時間当たりの心拍数は比較的多めで睡眠深度が
深い場合は単位時間当たりの心拍数は比較的少なめで非
常に安定している。図17(a)に覚醒状態、レム睡眠
期、ノンレム睡眠期の浅い睡眠状態、ノンレム睡眠期の
深い睡眠状態の、4通りの睡眠状態の一晩での推移の一
例を示し、また図17(b)にこのときの心拍数の推移
の一例を示す。ここで図17において、「WAKE」は
覚醒状態、「REM」はレム睡眠期、「non−RE
M」はノンレム睡眠期の浅い睡眠状態、「SWS」はノ
ンレム睡眠期の深い睡眠状態を示している。
[0004] That is, in night sleep, it is known that the heart rate per unit time decreases with falling asleep, shows a minimum value in the morning, and increases as the awakening time approaches. In the non-REM sleep period, the heart rate per unit time is relatively stable, but in the REM sleep period, the autonomic nervous system activity is disturbed, and the heart rate remarkably fluctuates. It is known that Further, even during the non-REM sleep period, when the sleep depth is shallow, the heart rate per unit time is relatively large, and when the sleep depth is deep, the heart rate per unit time is relatively small and very stable. FIG. 17 (a) shows an example of transitions of four sleep states overnight, including an awake state, a REM sleep period, a light sleep state during a non-REM sleep period, and a deep sleep state during a non-REM sleep period, and FIG. b) shows an example of the change of the heart rate at this time. In FIG. 17, “WAKE” is awake, “REM” is REM sleep, and “non-RE”.
“M” indicates a light sleep state during the non-REM sleep period, and “SWS” indicates a deep sleep state during the non-REM sleep period.

【0005】そしてこのような知見に基づいて、睡眠状
態を検出する構成が従来から種々提供されており、例え
ば特開昭63−205592号公報に開示されているよ
うに、脈拍数の増減を指標としてレム睡眠期を検出する
ものがある。また、特開平3−41926号公報に開示
されているように、脈拍数や呼吸数の増減傾向と時間的
変動とを織り込んだ睡眠指数を用いることにより睡眠の
状態変化を精度良く検出できるようにしたものがある。
Various configurations for detecting a sleep state based on such knowledge have been conventionally provided. For example, as disclosed in Japanese Patent Application Laid-Open No. 63-205592, an index of change in pulse rate is used as an index. There is one that detects the REM sleep period. Also, as disclosed in Japanese Patent Application Laid-Open No. 3-41926, it is possible to accurately detect a change in sleep state by using a sleep index incorporating the increasing and decreasing tendency of the pulse rate and respiratory rate and the temporal variation. There is something.

【0006】[0006]

【発明が解決しようとする課題】しかしながら、上記の
従来のものでは、脈拍数や呼吸数の増減傾向と時間的変
動とを織り込んだ睡眠指数を用いることにより、入眠を
判別し、レム睡眠期を検出していたものであるから、実
際の睡眠の状態変化との一致率に個人差がかなりあり、
睡眠ポリグラフに比較して精度がかなり悪く、深い睡眠
状態と浅い睡眠状態を判別できないという問題を有して
いた。
However, in the above-mentioned conventional apparatus, falling asleep is discriminated by using a sleep index incorporating the increasing and decreasing tendency of the pulse rate and respiratory rate and temporal fluctuation, and the REM sleep period is determined. Because it was detected, there is considerable individual difference in the matching rate with the actual change in sleep state,
The accuracy was considerably lower than that of polysomnography, and there was a problem that a deep sleep state and a light sleep state could not be distinguished.

【0007】本発明は上記問題点の解決を目的とするも
のであり、睡眠状態を精度良く検出することができる睡
眠状態判定方法および睡眠状態判定装置を提供しようと
するものである。
An object of the present invention is to solve the above problems, and an object of the present invention is to provide a sleep state judging method and a sleep state judging device capable of accurately detecting a sleep state.

【0008】[0008]

【課題を解決するための手段】本発明の請求項1に係る
睡眠状態判定方法は、心拍数または脈拍数などの生体情
報を計測して、生体情報の時系列の時間変化の増減傾向
を表すトレンド曲線を算出し、このトレンド曲線をもと
に睡眠状態を判定する睡眠状態判定方法において、トレ
ンド曲線を予め設定された時間の窓関数を用いて算出す
ることを特徴とするものである。
According to a first aspect of the present invention, there is provided a sleep state determining method which measures biological information such as a heart rate or a pulse rate and indicates a tendency of a time series of the biological information to increase or decrease over time. In a sleep state determination method for calculating a trend curve and determining a sleep state based on the trend curve, the trend curve is calculated using a window function of a preset time.

【0009】また請求項2の発明は、請求項1におい
て、窓関数を用いてトレンド曲線を算出するにあたっ
て、窓関数の時間の長さを40分〜70分とすることを
特徴とするものである。
According to a second aspect of the present invention, in the first aspect, when calculating the trend curve using the window function, the length of time of the window function is set to 40 to 70 minutes. is there.

【0010】本発明の請求項3に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部1
と、測定時刻を設定する測定時刻設定部2と、生体情報
値の時系列の時間変化の傾向を表すように設定したトレ
ンド曲線を算出するトレンド曲線算出部3と、上記生体
情報値とトレンド曲線から各測定時刻での生体情報値の
増加度合いを算出する増分算出部4と、上記各測定時刻
の前の一定時間内での上記生体情報値のばらつき度合い
を算出するばらつき度算出部5と、この算出された生体
情報値の増分とばらつき度の組み合わせから、予め設定
された睡眠状態出現確率密度分布表6をもとに各睡眠状
態の出現確率を算出する睡眠状態出現確率算出部7と、
睡眠状態出現確率算出部7で算出された各睡眠状態の出
現確率のなかで最も大きな数値を採った状態をその測定
時刻の睡眠状態と推定する睡眠状態推定部8とを具備し
て成ることを特徴とするものである。
According to a third aspect of the present invention, there is provided a measuring unit for measuring a biological information value such as a heart rate or a pulse rate per unit time to obtain a biological information measurement value at each time.
A measurement time setting unit 2 for setting a measurement time; a trend curve calculation unit 3 for calculating a trend curve set so as to represent a tendency of a time series change of the biological information value; An increment calculator 4 that calculates the degree of increase of the biological information value at each measurement time, a variation degree calculator 5 that calculates the degree of variation of the biological information value within a fixed time before each of the measurement times, A sleep state appearance probability calculation unit 7 that calculates an appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table 6 from a combination of the calculated increment of the biological information value and the degree of dispersion;
A sleep state estimating unit 8 for estimating the state that takes the largest numerical value among the appearance probabilities of each sleep state calculated by the sleep state appearance probability calculation unit 7 as the sleep state at the measurement time. It is a feature.

【0011】本発明の請求項4に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部1
と、測定時刻を設定する測定時刻設定部2と、生体情報
値の時系列の時間変化の傾向を表すように設定したトレ
ンド曲線を算出するトレンド曲線算出部3と、上記生体
情報値とトレンド曲線から各測定時刻での生体情報値の
増加度合いを算出する増分算出部4と、上記各測定時刻
の前の一定時間内での上記生体情報値のばらつき度合い
を算出するばらつき度算出部5と、この算出された生体
情報値の増分とばらつき度の組み合わせから、予め設定
された睡眠状態出現確率密度分布表6をもとに各睡眠状
態の出現確率を算出する睡眠状態出現確率算出部7と、
動作開始時からの経過時間から、予め設定された各睡眠
状態の出現割合表9をもとに各睡眠状態の出現割合を算
出する睡眠状態出現割合算出部10と、上記各睡眠状態
出現確率と上記各睡眠状態出現割合を掛け合わせて各睡
眠状態出現比率を算出する睡眠状態出現比率算出部11
と、上記各睡眠状態出現比率のなかで最も大きな数値を
採った状態をその測定時刻の睡眠状態と推定する睡眠状
態推定部8とを具備して成ることを特徴とするものであ
る。
According to a fourth aspect of the present invention, there is provided a sleep state judging device which measures a biological information value such as a heart rate or a pulse rate per unit time and obtains a biological information measurement value at each time.
A measurement time setting unit 2 for setting a measurement time; a trend curve calculation unit 3 for calculating a trend curve set so as to represent a tendency of a time series change of the biological information value; An increment calculator 4 that calculates the degree of increase of the biological information value at each measurement time, a variation degree calculator 5 that calculates the degree of variation of the biological information value within a fixed time before each of the measurement times, A sleep state appearance probability calculation unit 7 that calculates an appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table 6 from a combination of the calculated increment of the biological information value and the degree of dispersion;
A sleep state appearance ratio calculation unit 10 that calculates an appearance ratio of each sleep state based on a preset appearance ratio table of each sleep state from an elapsed time from the start of the operation; Sleep state appearance ratio calculation unit 11 that calculates each sleep state appearance ratio by multiplying each sleep state appearance ratio
And a sleep state estimating unit 8 for estimating the state having the largest numerical value among the sleep state appearance ratios as the sleep state at the measurement time.

【0012】本発明の請求項5に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部1
と、測定時刻を設定する測定時刻設定部2と、上記生体
情報値の時系列の進行方向の時間変化の傾向を表すよう
に設定した前トレンド曲線を算出する前トレンド曲線算
出部12と、上記生体情報値の時系列の逆行方向の時間
変化の傾向を表すように設定した後トレンド曲線を算出
する後トレンド曲線算出部13と、上記前トレンド曲線
と上記後トレンド曲線のうち各測定時刻において大きい
ほうの値を採用してトレンド曲線を算出するトレンド曲
線算出部14と、上記生体情報値とトレンド曲線から各
測定時刻での生体情報値の増加度合いを算出する増分算
出部4と、上記各測定時刻の前の一定時間内での上記生
体情報値のばらつき度合いを算出するばらつき度算出部
5と、この算出された生体情報値の増分とばらつき度の
組み合わせから、予め設定された睡眠状態出現確率密度
分布表6をもとに各睡眠状態の出現確率を算出する睡眠
状態出現確率算出部7と、上記各睡眠状態の出現確率の
なかで最も大きな数値を採った状態をその測定時刻の睡
眠状態と推定する睡眠状態推定部8とを具備して成るこ
とを特徴とするものである。
According to a fifth aspect of the present invention, there is provided a sleep state judging device which measures a biological information value such as a heart rate or a pulse rate per unit time and obtains a biological information measurement value at each time.
A measurement time setting unit 2 for setting a measurement time; a pre-trend curve calculation unit 12 for calculating a pre-trend curve set so as to represent a tendency of a temporal change of the biological information value in a time series traveling direction; A post-trend curve calculating unit 13 that calculates a trend curve after setting so as to represent the tendency of the time-series change of the biological information value in the time series in the reverse direction, and which is larger at each measurement time among the pre-trend curve and the post-trend curve. A trend curve calculating unit 14 for calculating a trend curve by using the two values, an increment calculating unit 4 for calculating the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve, and A variation degree calculation unit 5 that calculates the variation degree of the biological information value within a certain time period before the time, and a combination of the increment of the calculated biological information value and the variation degree, The sleep state appearance probability calculation unit 7 that calculates the appearance probability of each sleep state based on the set sleep state appearance probability density distribution table 6 and the largest numerical value among the appearance probabilities of the respective sleep states are taken. And a sleep state estimating unit 8 for estimating the state as the sleep state at the measurement time.

【0013】本発明の請求項6に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部1
と、測定時刻を設定する測定時刻設定部2と、上記生体
情報値の時系列の進行方向の時間変化の傾向を表すよう
に設定した前トレンド曲線を算出する前トレンド曲線算
出部12と、上記生体情報値の時系列の逆行方向の時間
変化の傾向を表すように設定した後トレンド曲線を算出
する後トレンド曲線算出部13と、上記前トレンド曲線
と上記後トレンド曲線のうち各測定時刻において大きな
ほうの値を採用してトレンド曲線を算出するトレンド曲
線算出部14と、上記生体情報値とトレンド曲線から各
測定時刻での生体情報値の増加度合いを算出する増分算
出部4と、上記各測定時刻の前の一定時間内での上記生
体情報値のばらつき度合いを算出するばらつき度算出部
5と、この算出された生体情報値の増分とばらつき度の
組み合わせから、予め設定された睡眠状態出現確率密度
分布表6をもとに各睡眠状態の出現確率を算出する睡眠
状態出現確率算出部7と、動作開始時からの経過時間か
ら、予め設定された睡眠状態の出現割合表9をもとに各
睡眠状態の出現割合を算出する睡眠状態出現割合算出部
10と、上記各睡眠状態出現確率と上記各睡眠状態出現
割合を掛け合わせて各睡眠状態出現比率を算出する睡眠
状態出現比率算出部11と、上記各睡眠状態出現比率の
なかで最も大きな数値を採った状態をその測定時刻の睡
眠状態と推定する睡眠状態推定部8とを具備して成るこ
とを特徴とするものである。
According to a sixth aspect of the present invention, there is provided a sleep state judging device which measures a biological information value such as a heart rate or a pulse rate per unit time and obtains a biological information measurement value at each time.
A measurement time setting unit 2 for setting a measurement time; a pre-trend curve calculation unit 12 for calculating a pre-trend curve set so as to represent a tendency of a temporal change of the biological information value in a time series traveling direction; A post-trend curve calculating unit 13 that sets a trend curve of the time series of the biological information value in the reverse direction and then calculates a trend curve, and a large value at each measurement time among the pre-trend curve and the post-trend curve. A trend curve calculating unit 14 for calculating a trend curve by using the two values, an increment calculating unit 4 for calculating the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve, and A variation degree calculation unit 5 that calculates the variation degree of the biological information value within a certain time period before the time, and a combination of the increment of the calculated biological information value and the variation degree, Sleep state appearance probability calculation unit 7 that calculates the appearance probability of each sleep state based on the sleep state appearance probability density distribution table 6 set in advance, and a sleep state that is set in advance from the elapsed time from the start of operation. A sleep state appearance ratio calculation unit 10 that calculates the appearance ratio of each sleep state based on the appearance ratio table 9, and calculates each sleep state appearance ratio by multiplying each sleep state appearance probability and each sleep state appearance ratio. Sleep state appearance ratio calculation unit 11 and a sleep state estimation unit 8 for estimating the state with the largest numerical value among the sleep state appearance ratios as the sleep state at the measurement time. It is assumed that.

【0014】また請求項7の発明は、請求項3〜6にお
いて、上記睡眠状態推定部8より出力された推定結果に
対して、入眠と推定した後、予め設定された時間の間は
深い睡眠状態にならないように推定結果を訂正する睡眠
判定調整部15を具備して成ることを特徴とするもので
ある。
According to a seventh aspect of the present invention, in the third to sixth aspects, after the sleep state estimating section 8 estimates that the sleep state is estimated to fall asleep, a deep sleep is performed for a preset time. It is characterized by including a sleep determination adjusting unit 15 for correcting the estimation result so as not to be in a state.

【0015】また請求項8の発明は、請求項3〜7にお
いて、上記睡眠状態推定部8より出力された推定結果に
対して、入眠と推定した後、ばらつき度が予め設定され
た閾値未満である状態が連続して予め設定された時間の
間持続した場合に、深い睡眠状態となるように推定結果
を訂正する睡眠判定調整部16を具備して成ることを特
徴とするものである。
According to the invention of claim 8, in the claims 3 to 7, the estimated result output from the sleep state estimating unit 8 is estimated to fall asleep and then the degree of variation is less than a preset threshold value. When a certain state is continuously maintained for a preset time, a sleep determination adjusting unit 16 that corrects an estimation result so as to be in a deep sleep state is provided.

【0016】[0016]

【発明の実施の形態】以下、本発明の実施の形態を説明
する。尚、以下の各実施の形態では、生体の活動状況の
指標となる生体情報値として、心拍数を用いているが、
脈拍数でも同様の処理を行なうことができるのはもちろ
んのことである。
Embodiments of the present invention will be described below. In the following embodiments, the heart rate is used as a biological information value that is an index of the activity state of a living body.
It goes without saying that the same processing can be performed for the pulse rate.

【0017】図1は請求項1,2の発明を実施するため
の装置の制御部の基本構成の一例を示すものである。心
拍の計測は、電極を用いたり電波を用いたりして心臓の
拍動を検出する心拍センサを使用して行なうことができ
るものであり、心拍センサからの出力は無線伝送あるい
は有線伝送により本装置に伝送される。本装置では、心
拍センサから伝送された心拍信号を受信して計測部1に
入力するようになっている。計測部1では、波形整形を
行うことによりパルス状の信号を得た後、単位時間毎の
パルス数を計数し、計数値を単位時間あたりの心拍数H
(t)として出力するようになっている。パルス数を計
数する単位時間は、たとえば1分である。計測部1から
出力される単位時間あたりの心拍数H(t)は、トレン
ド曲線算出部3と睡眠状態推定部8に入力される。
FIG. 1 shows an example of a basic configuration of a control unit of an apparatus for carrying out the first and second aspects of the present invention. The measurement of the heart rate can be performed using a heart rate sensor that detects the beat of the heart using electrodes or radio waves, and the output from the heart rate sensor is transmitted by wireless transmission or wired transmission to this device. Is transmitted to In this device, the heart rate signal transmitted from the heart rate sensor is received and input to the measuring unit 1. After obtaining a pulse-like signal by performing waveform shaping, the measuring unit 1 counts the number of pulses per unit time, and counts the counted value as the heart rate H per unit time.
(T). The unit time for counting the number of pulses is, for example, one minute. The heart rate per unit time H (t) output from the measurement unit 1 is input to the trend curve calculation unit 3 and the sleep state estimation unit 8.

【0018】トレンド曲線算出部3では、生体情報値の
時系列の時間変化での増減傾向を表すように設定したト
レンド曲線を、予め設定された時間(h)の窓関数を用
いて算出するようにしてある。窓関数は、ある連続した
データに対して解析を行なう場合に、解析を行なう区間
を制限するために用いられる関数であり、本発明におい
て窓関数の時間は、心拍数が増加傾向を示すREM睡眠
期の特徴を抽出するため、REM睡眠期の持続時間より
長いほうが望ましい。また逆に、心拍数のトレンドをよ
り良く反映させるためには、窓関数の時間は短いほうが
良い。この2つの条件を合わせて考慮すると、窓関数の
時間は40分〜70分の間、なかでも60分程度が最適
であると考えられる。
The trend curve calculation unit 3 calculates a trend curve set so as to represent an increasing / decreasing tendency in a time series of a biological information value by using a window function of a preset time (h). It is. The window function is a function used to limit the section to be analyzed when performing analysis on certain continuous data. In the present invention, the time of the window function is determined by the REM sleep in which the heart rate tends to increase. In order to extract the characteristics of the period, it is preferable that the period is longer than the duration of the REM sleep period. Conversely, in order to better reflect the trend of the heart rate, the time of the window function should be shorter. Considering these two conditions together, it is considered that the window function time is optimally between 40 minutes and 70 minutes, especially about 60 minutes.

【0019】続いて、トレンド曲線算出部3は図2のフ
ローチャートに示すように動作し、入力された単位時間
あたりの心拍数H(t)より次式のようにトレンド曲線
FT(t)を定義する。
Subsequently, the trend curve calculator 3 operates as shown in the flowchart of FIG. 2, and defines the trend curve FT (t) from the input heart rate H (t) per unit time as follows. I do.

【0020】 t=1の時 FT(1)=H(1) 1<t<hの時 H(t)≧FT(t−1)の時 FT(t)=FT(t−1) H(t)<FT(t−1)の時 FT(t)=H(t) h≦tの時 FT(t)=min(H(t−h+1)..H(t)) このようにして求められるトレンド曲線FT(t)は、
例えば、図3のようになる。この図3において、心拍数
は黒細線で、トレンド曲線は黒太線で表している。
When t = 1 FT (1) = H (1) When 1 <t <h When H (t) ≧ FT (t−1) FT (t) = FT (t−1) H ( t) <FT (t-1) FT (t) = H (t) When h ≦ t FT (t) = min (H (t−h + 1)... H (t)) The trend curve FT (t) is
For example, as shown in FIG. In FIG. 3, the heart rate is represented by a thin black line, and the trend curve is represented by a thick black line.

【0021】このようにしてトレンド曲線算出部3で算
出されたトレンド曲線FT(t)は、トレンド曲線算出
部3から出力されて睡眠状態推定部8に入力される。そ
して睡眠状態推定部8では、トレンド曲線FT(t)や
心拍数H(t)を用いて、覚醒状態、レム睡眠期、ノン
レム睡眠期の浅い睡眠状態、ノンレム睡眠期の深い睡眠
状態の、4種類のいずれの睡眠状態であるかを推定して
判定することができるものである。トレンド曲線FT
(t)や心拍数H(t)を用いて4種類のいずれの睡眠
状態であるかを推定する方法は、特に限定されるもので
はないが、後述の図4や図10などの場合と同様に、ト
レンド曲線FT(t)と心拍数H(t)の差から心拍数
の増分を算出すると共に、ある1分の心拍数とその次の
1分の心拍数の差を順次計算してばらつき度を算出し、
心拍数の増分とばらつき度の組み合わせにより4種類の
睡眠状態の出現確率を求め、さらに入床からの経過時間
によって4種類の睡眠状態の出現割合を求め、出現確率
と出現割合を掛け合わせた数字(出現比率)が一番大き
な睡眠状態を、その時点の睡眠状態と推定して判定する
ことができるものである。
The trend curve FT (t) calculated by the trend curve calculating section 3 is output from the trend curve calculating section 3 and input to the sleep state estimating section 8. Then, the sleep state estimating unit 8 uses the trend curve FT (t) and the heart rate H (t) to calculate the 4th state of the awake state, the REM sleep period, the light sleep state of the non-REM sleep period, and the deep sleep state of the non-REM sleep period. It can be determined by estimating which kind of sleep state it is. Trend curve FT
(T) and the method of estimating which of the four sleep states it is using heart rate H (t) are not particularly limited, but are the same as in the case of FIGS. In addition, the increase in the heart rate is calculated from the difference between the trend curve FT (t) and the heart rate H (t), and the difference between the one-minute heart rate and the next one-minute heart rate is sequentially calculated to obtain a variation. Calculate the degree,
The probability of appearance of four types of sleep states is obtained from the combination of the increase in heart rate and the degree of variability, and the appearance ratio of four types of sleep states is obtained based on the elapsed time since entering the bed. The sleep state with the largest (appearance ratio) can be determined by estimating the sleep state at that time.

【0022】図4は請求項3の発明に係る装置の制御部
の基本構成の一例を示すものである。心拍センサからの
出力は無線伝送あるいは有線伝送により本装置に伝送さ
れ、心拍センサから伝送された心拍信号は計測部1に入
力されるようになっている。計測部1では、波形整形を
行うことによりパルス状の信号を得た後、単位時間毎の
パルス数を計数し、計数値を単位時間あたりの心拍数H
(t)として出力するようになっている。パルス数を計
数する単位時間は、たとえば1分である。計測部1から
出力される単位時間あたりの心拍数H(t)は、トレン
ド曲線算出部3と増分算出部4とばらつき度算出部5に
入力される。
FIG. 4 shows an example of the basic configuration of the control unit of the device according to the third aspect of the present invention. The output from the heart rate sensor is transmitted to the apparatus by wireless transmission or wired transmission, and the heart rate signal transmitted from the heart rate sensor is input to the measuring unit 1. After obtaining a pulse-like signal by performing waveform shaping, the measuring unit 1 counts the number of pulses per unit time, and counts the counted value as the heart rate H per unit time.
(T). The unit time for counting the number of pulses is, for example, one minute. The heart rate H (t) per unit time output from the measurement unit 1 is input to the trend curve calculation unit 3, the increase calculation unit 4, and the variation calculation unit 5.

【0023】トレンド曲線算出部3では、生体情報値の
時系列の時間変化での増減傾向を表すように設定したト
レンド曲線を算出するようにしてある。すなわち、トレ
ンド曲線算出部3は図5のフローチャートに示すように
動作し、入力された単位時間あたりの心拍数H(t)よ
り次式のようにトレンド曲線FT(t)を定義する。
The trend curve calculating section 3 calculates a trend curve set so as to represent an increasing / decreasing tendency of a biological information value with time in a time series. That is, the trend curve calculation unit 3 operates as shown in the flowchart of FIG. 5, and defines the trend curve FT (t) from the input heart rate H (t) per unit time as in the following equation.

【0024】 t=1の時 FT(1)=H(1) t>1の時 H(t)≧FT(t−1)の時 FT(t)=FT(t−1) H(t)<FT(t−1)の時 FT(t)=H(t) このようにして求められるトレンド曲線FT(t)は、
たとえば、図6のようになる。この図6において、心拍
数は黒細線で、トレンド曲線は黒太線で表している。
When t = 1 FT (1) = H (1) When t> 1 When H (t) ≧ FT (t−1) FT (t) = FT (t−1) H (t) <FT (t-1) FT (t) = H (t) The trend curve FT (t) thus obtained is
For example, as shown in FIG. In FIG. 6, the heart rate is indicated by a thin black line, and the trend curve is indicated by a thick black line.

【0025】このようにしてトレンド曲線算出部3で算
出されて出力されたトレンド曲線FT(t)は各測定時
刻での生体情報値の増加度合いである増分を算出する増
分算出部4に入力される。増分算出部4では、トレンド
曲線FT(t)からの心拍数H(t)の増加度合いを表
す指標である増分ZOU(t)が算出される。すなわ
ち、増分ZOU(t)を心拍数H(t)とトレンド曲線
FT(t)を用いて次式のように定義する。
The trend curve FT (t) calculated and output by the trend curve calculating section 3 is input to an increment calculating section 4 for calculating an increment which is the degree of increase of the biological information value at each measurement time. You. The increment calculation unit 4 calculates an increment ZOU (t) which is an index indicating the degree of increase of the heart rate H (t) from the trend curve FT (t). That is, the increment ZOU (t) is defined as follows using the heart rate H (t) and the trend curve FT (t).

【0026】ZOU(t)=H(t)−FT(t) このようにして求められる増分ZOU(t)は、たとえ
ば、図7のようになる。図7は、図6の心拍数H(t)
とトレンド曲線FT(t)から得られる一晩の増分の推
移を示すものである。
ZOU (t) = H (t) -FT (t) The thus obtained increment ZOU (t) is, for example, as shown in FIG. FIG. 7 shows the heart rate H (t) of FIG.
And the trend of the overnight increment obtained from the trend curve FT (t).

【0027】またばらつき度算出部5では、計測部1よ
り入力された心拍数H(t)のばらつき度合いを算出す
るようになっている。まず、計測された心拍数H(t)
とその1分前に計測された心拍数H(t−1)の差の絶
対値(差分SABUN(t)とする)を算出する。差分
SABUN(t)の算出式は次式のように定義する。
The variability calculating section 5 calculates the variability of the heart rate H (t) input from the measuring section 1. First, the measured heart rate H (t)
And the absolute value of the difference between the heart rate H (t-1) and the heart rate H (t-1) measured one minute before that (referred to as difference SABUN (t)). The formula for calculating the difference SABUN (t) is defined as follows.

【0028】 abs(H(t)−H(t−1))<6の場合 SABUN(t)=abs(H(t)−H(t−1)) abs(H(t)−H(t−1))≧6の場合 SABUN(t)=6 さらに、差分SABUN(t)を用いてばらつき度BA
RA(t)を次式のように算出する。
If abs (H (t) -H (t-1)) <6, SABUN (t) = abs (H (t) -H (t-1)) abs (H (t) -H (t) -1)) ≧ 6 SABUN (t) = 6 Further, the degree of variation BA is calculated using the difference SABUN (t).
RA (t) is calculated as in the following equation.

【0029】t=1の場合 BARA(1)=SABUN(1) t=2の場合 BARA(1)=SABUN(1)+SABUN(2) t=3の場合 BARA(1)=SABUN(1)+SABUN(2)
+SABUN(3) t=4の場合 BARA(1)=SABUN(1)+SABUN(2)
+SABUN(3)+SABUN(4) t≧5の場合 BARA(t)=SABUN(t−4)+SABUN
(t−3)+SABUN(t−2)+SABUN(t−
1)+SABUN(t) このようにして求めたばらつき度BARA(t)は、例
えば図8のようになる。
When t = 1 BARA (1) = SABUN (1) When t = 2 BARA (1) = SABUN (1) + SABUN (2) When t = 3 BARA (1) = SABUN (1) + SABUN (2)
+ SABUN (3) When t = 4 BARA (1) = SABUN (1) + SABUN (2)
+ SABUN (3) + SABUN (4) When t ≧ 5 BARA (t) = SABUN (t−4) + SABUN
(T-3) + SABUN (t-2) + SABUN (t-
1) + SABUN (t) The variation degree BARA (t) obtained as described above is, for example, as shown in FIG.

【0030】そして上記のようにして、増分算出部4で
算出されて出力された増分ZOU(t)と、ばらつき度
算出部5で算出されて出力されたばらつき度BARA
(t)は、睡眠状態出現確率密度分布表6をもとに各睡
眠状態の出現確率を算出する睡眠状態出現確率算出部7
に入力される。睡眠状態出現確率算出部7では、増分算
出部4から出力された増分ZOU(t)とばらつき度算
出部5から出力されたばらつき度BARA(t)の組み
合わせから、各睡眠状態ごとの睡眠状態出現確率密度分
布表6をもとにして、覚醒状態、レム睡眠期、ノンレム
睡眠期の浅い睡眠状態、ノンレム睡眠期の深い睡眠状態
の、4種類の各睡眠状態の出現確率を算出するようにな
っている。
As described above, the increment ZOU (t) calculated and output by the increment calculating section 4 and the variation BARA calculated and output by the variation calculating section 5
(T) is a sleep state appearance probability calculation unit 7 that calculates the appearance probability of each sleep state based on the sleep state appearance probability density distribution table 6.
Is input to The sleep state appearance probability calculation unit 7 calculates a sleep state appearance for each sleep state from a combination of the increment ZOU (t) output from the increment calculation unit 4 and the variation BARA (t) output from the variation degree calculation unit 5. Based on the probability density distribution table 6, the appearance probabilities of the four types of sleep states of the awake state, the REM sleep period, the light sleep state of the non-REM sleep period, and the deep sleep state of the non-REM sleep period are calculated. ing.

【0031】睡眠状態出現確率密度分布表6は、たとえ
ば89夜分(時間では約4万2000分)の夜間の睡眠
状態のデータを基に作成されているものである。すなわ
ち、この約4万2000分のデータの1分1分ごとに心
拍数の増分とばらつき度が計算されており、その1分1
分ごとに睡眠状態が脳波のデータから判定されいてる。
そしてこれらのデータから、1分の心拍数の増分とばら
つき度の組み合わせと、睡眠状態の種類との関係を解析
することによって、その1分の睡眠状態が、ノンレム睡
眠期の浅い睡眠状態である確率、レム睡眠期である確
率、ノンレム睡眠期の深い睡眠状態である確率、覚醒状
態である確率をそれぞれ算出することができる。睡眠状
態出現確率密度分布表6はこれらの確率をまとめたもの
であって、図9にその一例を示すように、ノンレム睡眠
期の浅い睡眠状態(non−REM)、レム睡眠期(R
EM)、ノンレム睡眠期の深い睡眠状態(SWS)、覚
醒状態(WAKE)の各睡眠状態ごとの出現確率密度分
布表6a,6b,6c,6dからなるものであり、各睡
眠状態ごとの睡眠状態出現確率密度分布表6a,6b,
6c,6dには増分とばらつき度の組み合わせからその
睡眠状態の出現確率が求められるように2次元の表にな
っている。例えば図9にみられるように、増分が
「3」、ばらつきが「5」であると、ノンレム睡眠期の
浅い睡眠状態(non−REM)の出現確率は68.5
%になる。そしてこのように睡眠状態出現確率算出部7
で算出されて出力された各睡眠状態ごとの出現確率は、
睡眠状態推定部8に入力される。
The sleep state appearance probability density distribution table 6 is created based on night sleep state data for 89 nights (about 42,000 minutes in time), for example. That is, the increment and the degree of variation of the heart rate are calculated every minute and minute of the data of about 42,000 minutes.
The sleep state is determined from the brain wave data every minute.
Then, by analyzing the relationship between the combination of the one-minute heart rate increment and the degree of variation and the type of sleep state from these data, the one-minute sleep state is a light sleep state in the non-REM sleep period. It is possible to calculate the probability, the probability of the REM sleep period, the probability of the deep sleep state in the non-REM sleep period, and the probability of the awake state, respectively. The sleep state appearance probability density distribution table 6 summarizes these probabilities, and as shown in FIG. 9 as an example, the light sleep state (non-REM) in the non-REM sleep period and the REM sleep period (R
EM), a non-REM sleep deep sleep state (SWS), and an awake state (WAKE) for each sleep state. Appearance probability density distribution tables 6a, 6b,
6c and 6d are two-dimensional tables so that the appearance probability of the sleep state can be obtained from the combination of the increment and the degree of variation. For example, as shown in FIG. 9, when the increment is “3” and the variation is “5”, the appearance probability of the light sleep state (non-REM) in the non-REM sleep period is 68.5.
%become. And the sleep state appearance probability calculation unit 7
The appearance probability for each sleep state calculated and output in
It is input to the sleep state estimation unit 8.

【0032】睡眠状態推定部8では、睡眠状態出現確率
算出部7から出力された4つの各睡眠状態ごとの出現確
率を比べ、最も大きな値を示した睡眠状態をその時刻の
睡眠状態と推定するようにしてあり、覚醒状態、レム睡
眠期、ノンレム睡眠期の浅い睡眠状態、ノンレム睡眠期
の深い睡眠状態の、4種類のいずれの睡眠状態であるか
を推定することができるものである。
The sleep state estimating unit 8 compares the appearance probabilities of the four sleep states output from the sleep state appearance probability calculation unit 7 and estimates the sleep state having the largest value as the sleep state at that time. As described above, it is possible to estimate any one of the four sleep states of the awake state, the REM sleep period, the light sleep state of the non-REM sleep period, and the deep sleep state of the non-REM sleep period.

【0033】次に、図10は請求項4の発明に係る装置
の制御部の基本構成の一例を示すものである。計測部
1、トレンド曲線算出部3、増分算出部4、ばらつき度
算出部5、睡眠状態出現確率算出部7は上記の図4の場
合と同じである。そして睡眠状態出現確率算出部7で算
出されて出力された各睡眠状態ごとの出現確率は、睡眠
状態出現比率算出部11に入力される。
Next, FIG. 10 shows an example of the basic configuration of the control unit of the apparatus according to the fourth aspect of the present invention. The measuring unit 1, the trend curve calculating unit 3, the increment calculating unit 4, the variation calculating unit 5, and the sleep state appearance probability calculating unit 7 are the same as those in FIG. The appearance probability for each sleep state calculated and output by the sleep state appearance probability calculation unit 7 is input to the sleep state appearance ratio calculation unit 11.

【0034】睡眠状態出現割合算出部10では、睡眠の
動作開始時からの経過時間から、各睡眠状態の出現割合
表9をもとに、覚醒状態、レム睡眠期、ノンレム睡眠期
の浅い睡眠状態、ノンレム睡眠期の深い睡眠状態の、4
種類の各睡眠状態の出現割合を算出するようになってい
る。この睡眠状態出現割合表9は、前記と同様な夜間の
睡眠状態のデータを基にして入床からの経過時間と睡眠
状態の種類との関係を解析することによって、4種類の
各睡眠状態の出現割合を算出して作成されたものであっ
て、睡眠の開始時から所定時間経過したときの睡眠状態
が、覚醒状態、レム睡眠期、ノンレム睡眠期の浅い睡眠
状態、ノンレム睡眠期の深い睡眠状態のいずれであるか
の割合を示すものであり、図11に覚醒状態(WAK
E)、レム睡眠期(REM)、ノンレム睡眠期の浅い睡
眠状態(non−REM)、ノンレム睡眠期の深い睡眠
状態(SWS)の、4通りの各睡眠状態の出現割合表9
をグラフ化したものを示す。睡眠状態出現割合算出部6
出算出されて出力された各睡眠状態の出現割合は、睡眠
状態出現比率算出部11に入力される。
The sleep state appearance ratio calculation unit 10 calculates the sleep state in the awake state, the REM sleep period, and the non-REM sleep period based on the appearance ratio table 9 of each sleep state based on the elapsed time from the start of the sleep operation. , Non-REM, deep sleep, 4
The appearance ratio of each type of sleep state is calculated. The sleep state appearance ratio table 9 analyzes the relationship between the elapsed time from the entrance and the type of sleep state based on the data of the sleep state at night, as described above, to thereby obtain four types of sleep states. It is created by calculating the appearance ratio, and the sleep state when a predetermined time has elapsed from the start of sleep is awake state, REM sleep period, light sleep state in non-REM sleep period, deep sleep in non-REM sleep period FIG. 11 shows the ratio of the awake state (WAK state).
E), REM sleep period (REM), non-REM light sleep state (non-REM), non-REM deep sleep state (SWS), appearance ratio of each of the four sleep states Table 9
Is shown as a graph. Sleep state appearance ratio calculation unit 6
The sleep state appearance ratio calculated and output is input to the sleep state appearance ratio calculation unit 11.

【0035】睡眠状態出現比率算出部11では、睡眠状
態出現確率算出部7から出力された各睡眠状態ごとの出
現確率と、睡眠状態出現割合算出部10から出力された
各睡眠状態の出現割合とを掛け合わせて、各睡眠状態の
出現比率を算出するようになっている。そして睡眠状態
出現比率算出部11で算出されて出力された各睡眠状態
の出現比率は、睡眠状態推定部8に入力される。
The sleep state appearance ratio calculation unit 11 calculates the appearance probability of each sleep state output from the sleep state appearance probability calculation unit 7 and the sleep state appearance ratio output from the sleep state appearance ratio calculation unit 10. To calculate the appearance ratio of each sleep state. The sleep state appearance ratio calculated and output by the sleep state appearance ratio calculation unit 11 is input to the sleep state estimation unit 8.

【0036】睡眠状態推定部8では、睡眠状態出現比率
算出部11から出力された各睡眠状態ごとの出現比率を
比べ、最も大きな値を示した状態をその時刻の睡眠状態
と推定するようにしてあり、覚醒状態、レム睡眠期、ノ
ンレム睡眠期の浅い睡眠状態、ノンレム睡眠期の深い睡
眠状態の、4種類のいずれの睡眠状態であるかを推定す
ることができるものである。
The sleep state estimating unit 8 compares the appearance ratios for each sleep state output from the sleep state appearance ratio calculation unit 11 and estimates the state showing the largest value as the sleep state at that time. Yes, it is possible to estimate which of the four sleep states, awake state, REM sleep period, light sleep state during non-REM sleep period, and deep sleep state during non-REM sleep period.

【0037】次に、図12は請求項5の発明に係る装置
の制御部の基本構成の一例を示すものである。計測部1
は上記の図4の場合と同じである。そして計測部1から
出力された単位時間あたりの心拍数H(t)は、前トレ
ンド曲線算出部12、後トレンド曲線算出部13、増分
算出部4、ばらつき度算出部5に出力される。
Next, FIG. 12 shows an example of the basic configuration of the control unit of the device according to the fifth aspect of the present invention. Measurement unit 1
Is the same as in the case of FIG. Then, the heart rate H (t) per unit time output from the measurement unit 1 is output to the front trend curve calculation unit 12, the rear trend curve calculation unit 13, the increment calculation unit 4, and the variation calculation unit 5.

【0038】前トレンド曲線算出部12では、生体情報
値の時系列の進行方向の時間変化での増減傾向を表すよ
うに設定した前トレンド曲線FT(t)を算出するよう
にしてある。この前トレンド曲線FT(t)を算出する
動作は図4でのトレンド曲線算出部3の動作と同じであ
る。そして前トレンド曲線算出部12から出力された前
トレンド曲線FT(t)は、トレンド曲線算出部14に
入力される。
The pre-trend curve calculating section 12 calculates a pre-trend curve FT (t) which is set so as to represent an increasing / decreasing tendency of the biological information value in the time series in the time direction. The operation of calculating the previous trend curve FT (t) is the same as the operation of the trend curve calculation unit 3 in FIG. Then, the front trend curve FT (t) output from the front trend curve calculation unit 12 is input to the trend curve calculation unit 14.

【0039】また後トレンド曲線算出部13では、生体
情報値の時系列の逆行方向の時間変化での増減傾向を表
すように設定した後トレンド曲線RT(t)を算出する
ようにしてある。すなわち、後トレンド曲線算出部13
は図13のフローチャートに示すように動作し、入力さ
れた単位時間あたりの心拍数H(t)より次式のように
後トレンド曲線RT(t)を定義する。
The post-trend curve calculating section 13 calculates a post-trend curve RT (t) that is set so as to represent an increasing / decreasing tendency in a time series in the reverse direction of the time series of the biological information value. That is, the post-trend curve calculation unit 13
Operates as shown in the flowchart of FIG. 13, and defines the post-trend curve RT (t) from the input heart rate H (t) per unit time as in the following equation.

【0040】 ds:最終の計測時間 t=dsの時 RT(ds)=H(ds) t<dsの時 H(t)≧RT(t+1)の時 RT(t)=RT(t+1) H(t)<RT(t+1)の時 RT(t)=H(t) そして後トレンド曲線算出部13で算出されて出力され
た後トレンド曲線RT(t)は、トレンド曲線算出部1
4に入力される。
Ds: Final measurement time When t = ds RT (ds) = H (ds) When t <ds H (t) ≧ RT (t + 1) RT (t) = RT (t + 1) H ( t) <RT (t + 1) RT (t) = H (t) Then, the post-trend curve RT (t) calculated and output by the post-trend curve calculation unit 13 is calculated by the trend curve calculation unit 1
4 is input.

【0041】トレンド曲線算出部14では、前トレンド
曲線算出部12から出力された前トレンド曲線FT
(t)と、後トレンド曲線算出部13から出力された後
トレンド曲線RT(t)のうち、各測定時刻においてよ
り大きな値を採るほうを採用し、トレンド曲線TRE
(t)を算出するようにしてある。このようにしてトレ
ンド曲線算出部14で算出されて出力されたトレンド曲
線TRE(t)は、各測定時刻での生体情報値の増加度
合いを算出する増分算出部4に入力される。
The trend curve calculating section 14 outputs the previous trend curve FT output from the previous trend curve calculating section 12.
(T) and the later trend curve RT (t) output from the later trend curve calculation unit 13 adopting the one that takes a larger value at each measurement time.
(T) is calculated. The trend curve TRE (t) calculated and output by the trend curve calculation unit 14 in this manner is input to the increment calculation unit 4 which calculates the degree of increase of the biological information value at each measurement time.

【0042】増分算出部4では、トレンド曲線TRE
(t)からの心拍数H(t)の増加度合いを表す指標で
ある増分ZOU(t)を算出する。すなわち、増分ZO
U(t)を心拍数H(t)とトレンド曲線TRE(t)
を用いて次式のように定義する。
In the increment calculating section 4, the trend curve TRE
An increment ZOU (t), which is an index indicating the degree of increase of the heart rate H (t) from (t), is calculated. That is, the incremental ZO
U (t) is expressed as heart rate H (t) and trend curve TRE (t).
Is defined as follows.

【0043】ZOU(t)=H(t)−TRE(t) ばらつき度算出部5については、図4のものと同じであ
る。
ZOU (t) = H (t) -TRE (t) The variation calculator 5 is the same as that in FIG.

【0044】そして、増分算出部4で算出されて出力さ
れた増分ZOU(t)と、ばらつき度算出部5で算出さ
れて出力されたばらつき度BARA(t)は、睡眠状態
出現確率密度分布表6をもとに各睡眠状態の出現確率を
算出する睡眠状態出現確率算出部7に入力され、各睡眠
状態ごとの出現確率が算出される。この各睡眠状態ごと
の出現確率は睡眠状態推定部8に入力され、覚醒状態、
レム睡眠期、ノンレム睡眠期の浅い睡眠状態、ノンレム
睡眠期の深い睡眠状態の、4種類のいずれの睡眠状態で
あるかを推定する。睡眠状態出現確率算出部7と睡眠状
態推定部8については図4のものと同じである。
The increment ZOU (t) calculated and output by the increment calculator 4 and the variance BARA (t) calculated and output by the variability calculator 5 are represented by a sleep state appearance probability density distribution table. 6 is input to a sleep state appearance probability calculation unit 7 that calculates the appearance probability of each sleep state, and the appearance probability for each sleep state is calculated. The appearance probability for each sleep state is input to the sleep state estimation unit 8, and the awake state,
It is estimated which of the four types of sleep states, the REM sleep phase, the light sleep state during the non-REM sleep phase, and the deep sleep state during the non-REM sleep phase. The sleep state appearance probability calculating section 7 and the sleep state estimating section 8 are the same as those in FIG.

【0045】次に、図14は請求項6の発明に係る装置
の制御部の基本構成の一例を示すものである。計測部
1、前トレンド曲線算出部12、後トレンド曲線算出部
13、トレンド曲線算出部14、増分算出部4、ばらつ
き度算出部5、睡眠状態出現確率算出部7については、
図12のものと同じである。また睡眠状態出現割合算出
部10、睡眠状態出現比率算出部11、睡眠状態推定部
8については図10のものと同じである。
FIG. 14 shows an example of the basic configuration of the control unit of the device according to the sixth aspect of the present invention. The measurement unit 1, the front trend curve calculation unit 12, the rear trend curve calculation unit 13, the trend curve calculation unit 14, the increment calculation unit 4, the variation degree calculation unit 5, and the sleep state appearance probability calculation unit 7
It is the same as that of FIG. The sleep state appearance ratio calculation unit 10, sleep state appearance ratio calculation unit 11, and sleep state estimation unit 8 are the same as those in FIG.

【0046】そして睡眠状態出現確率算出部7で算出さ
れて出力された各睡眠状態の出現確率と、睡眠状態出現
割合算出部10で算出されて出力された各睡眠状態の出
現割合は、それぞれ睡眠状態出現比率算出部11に入力
され、各睡眠状態の出現比率が算出される。この各睡眠
状態の出現比率は睡眠状態推定部8に入力され、覚醒状
態、レム睡眠期、ノンレム睡眠期の浅い睡眠状態、ノン
レム睡眠期の深い睡眠状態の、4種類のいずれの睡眠状
態であるかを推定する。
The appearance probability of each sleep state calculated and output by the sleep state appearance probability calculation unit 7 and the appearance ratio of each sleep state calculated and output by the sleep state appearance ratio calculation unit 10 are respectively The appearance ratio is input to the state appearance ratio calculation unit 11, and the appearance ratio of each sleep state is calculated. The appearance ratio of each sleep state is input to the sleep state estimation unit 8, and is any one of four sleep states: awake state, REM sleep period, light sleep state during non-REM sleep period, and deep sleep state during non-REM sleep period. Estimate.

【0047】次に図15は請求項7の発明に係る装置の
制御部の基本構成の一例を示すものである。計測部1か
ら睡眠状態推定部8までについては図4のものと同じで
ある。
FIG. 15 shows an example of the basic configuration of the control unit of the device according to the seventh aspect of the present invention. The sections from the measuring section 1 to the sleep state estimating section 8 are the same as those in FIG.

【0048】そして睡眠状態推定部8で推定された睡眠
状態の推定結果は、睡眠状態推定部8から出力されて睡
眠判定調整部15に入力される。睡眠判定調整部15で
は、入力された睡眠状態の推定結果をもとに、入眠と推
定した後、予め設定してあった時間(例えば10分)の
間は深い睡眠状態と判定しないように推定結果を訂正す
るようになっている。
The sleep state estimation result estimated by the sleep state estimation unit 8 is output from the sleep state estimation unit 8 and input to the sleep determination adjustment unit 15. The sleep determination adjustment unit 15 estimates that the user has fallen asleep based on the input estimation result of the sleep state, and then estimates that the user is not determined to be in the deep sleep state for a preset time (for example, 10 minutes). The result is corrected.

【0049】次に図16は請求項8の発明に係る装置の
制御部の基本構成の一例を示すものである。計測部1か
ら睡眠状態推定部8までについては図4のものと同じで
ある。
Next, FIG. 16 shows an example of the basic configuration of the control unit of the device according to the eighth aspect of the present invention. The sections from the measuring section 1 to the sleep state estimating section 8 are the same as those in FIG.

【0050】そして睡眠状態推定部8で推定された睡眠
状態の推定結果は、睡眠状態推定部8から出力されて睡
眠判定調整部16に入力される。睡眠判定調整部16で
は、入力された睡眠状態の推定結果をもとに、入眠と推
定した後、ばらつき度が予め設定してある閾値未満であ
る状態が、連続して予め設定してある時間の間持続した
場合には、深い睡眠状態と判定するように推定結果を訂
正するようになっている。
The sleep state estimation result estimated by the sleep state estimation unit 8 is output from the sleep state estimation unit 8 and input to the sleep determination adjustment unit 16. The sleep determination adjustment unit 16 estimates that the user falls asleep based on the input estimation result of the sleep state, and then continuously changes the state in which the degree of variation is less than the preset threshold for a preset period of time. If the duration is maintained for a while, the estimation result is corrected so as to determine that the subject is in a deep sleep state.

【0051】[0051]

【発明の効果】上記のように本発明の請求項1に係る睡
眠状態判定方法は、心拍数または脈拍数などの生体情報
を計測して、生体情報の時系列の時間変化の増減傾向を
表すトレンド曲線を算出し、このトレンド曲線をもとに
睡眠状態を判定する睡眠状態判定方法において、トレン
ド曲線を予め設定された時間の窓関数を用いて算出する
ようにしたので、睡眠状態を精度良く、しかもリアルタ
イムで検出することができるものである。
As described above, in the sleep state determination method according to the first aspect of the present invention, biological information such as a heart rate or a pulse rate is measured, and a time series of the biological information is represented by an increasing or decreasing tendency. In the sleep state determination method of calculating a trend curve and determining a sleep state based on the trend curve, the trend curve is calculated using a window function of a preset time, so that the sleep state can be accurately calculated. In addition, it can be detected in real time.

【0052】また請求項2の発明は、窓関数を用いてト
レンド曲線を算出するにあたって、窓関数の時間の長さ
を40分〜70分としたので、睡眠状態を精度良く検出
することができるものである。
According to the second aspect of the present invention, when calculating the trend curve using the window function, the length of time of the window function is set to 40 to 70 minutes, so that the sleep state can be detected with high accuracy. Things.

【0053】本発明の請求項3に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部
と、測定時刻を設定する測定時刻設定部と、生体情報値
の時系列の時間変化の傾向を表すように設定したトレン
ド曲線を算出するトレンド曲線算出部と、上記生体情報
値とトレンド曲線から各測定時刻での生体情報値の増加
度合いを算出する増分算出部と、上記各測定時刻の前の
一定時間内での上記生体情報値のばらつき度合いを算出
するばらつき度算出部と、この算出された生体情報値の
増分とばらつき度の組み合わせから、予め設定された睡
眠状態出現確率密度分布表をもとに各睡眠状態の出現確
率を算出する睡眠状態出現確率算出部と、睡眠状態出現
確率算出部で算出された各睡眠状態の出現確率のなかで
最も大きな数値を採った状態をその測定時刻の睡眠状態
と推定する睡眠状態推定部とを具備するので、生体情報
値の増加度合いとばらつき度合いとから各睡眠状態の出
現確率を算出することができ、睡眠状態を精度良く、し
かもリアルタイムで検出することができるものである。
According to a third aspect of the present invention, there is provided a sleep state determining apparatus which measures a biological information value such as a heart rate or a pulse rate per unit time and obtains a biological information measured value at each time; A measurement time setting unit that sets a trend curve calculation unit that calculates a trend curve set to represent the tendency of the time change of the biological information value over time, and a measurement curve at each measurement time from the biological information value and the trend curve. An increment calculation unit that calculates the degree of increase of the biological information value, and a variation degree calculation unit that calculates the degree of variation of the biological information value within a certain time period before each of the measurement times, A sleep state appearance probability calculation unit that calculates an appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table from a combination of the increment and the degree of variation, and a sleep state appearance probability calculation unit. A sleep state estimating unit that estimates the state that takes the largest numerical value among the appearance probabilities of the sleep states as the sleep state at the measurement time. The state appearance probability can be calculated, and the sleep state can be detected accurately and in real time.

【0054】本発明の請求項4に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部
と、測定時刻を設定する測定時刻設定部と、生体情報値
の時系列の時間変化の傾向を表すように設定したトレン
ド曲線を算出するトレンド曲線算出部と、上記生体情報
値とトレンド曲線から各測定時刻での生体情報値の増加
度合いを算出する増分算出部と、上記各測定時刻の前の
一定時間内での上記生体情報値のばらつき度合いを算出
するばらつき度算出部と、この算出された生体情報値の
増分とばらつき度の組み合わせから、予め設定された睡
眠状態出現確率密度分布表をもとに各睡眠状態の出現確
率を算出する睡眠状態出現確率算出部と、動作開始時か
らの経過時間から、予め設定された各睡眠状態の出現割
合表をもとに各睡眠状態の出現割合を算出する睡眠状態
出現割合算出部と、上記各睡眠状態出現確率と上記各睡
眠状態出現割合を掛け合わせて各睡眠状態出現比率を算
出する睡眠状態出現比率算出部と、上記各睡眠状態出現
比率のなかで最も大きな数値を採った状態をその測定時
刻の睡眠状態と推定する睡眠状態推定部とを具備するの
で、生体情報値の増加度合いとばらつき度合いに加え
て、睡眠状態出現割合表から各睡眠状態の出現割合を算
出することができ、時間経過の概念を加えることができ
て、睡眠状態を精度良く、しかもリアルタイムで検出す
ることができるものである。
According to a fourth aspect of the present invention, there is provided a sleep state judging device for measuring a biological information value such as a heart rate or a pulse rate per unit time to obtain a biological information measured value at each time, A measurement time setting unit that sets a trend curve calculation unit that calculates a trend curve set to represent the tendency of the time change of the biological information value over time, and a measurement curve at each measurement time from the biological information value and the trend curve. An increment calculation unit that calculates the degree of increase of the biological information value, and a variation degree calculation unit that calculates the degree of variation of the biological information value within a certain time period before each of the measurement times, From the combination of the increment and the degree of variation, from a sleep state appearance probability calculation unit that calculates the appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table, and from the elapsed time from the start of operation, The sleep state appearance ratio calculation unit that calculates the appearance ratio of each sleep state based on the appearance ratio table of each sleep state that has been set, and multiplies each sleep state appearance probability by each sleep state appearance ratio Since it has a sleep state appearance ratio calculation unit that calculates the sleep state appearance ratio, and a sleep state estimation unit that estimates the state that takes the largest numerical value among the sleep state appearance ratios as the sleep state at the measurement time. In addition to the increase degree and the variation degree of the biological information value, the appearance ratio of each sleep state can be calculated from the sleep state appearance ratio table, and the concept of time lapse can be added, and the sleep state can be accurately calculated. Moreover, it can be detected in real time.

【0055】本発明の請求項5に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部
と、測定時刻を設定する測定時刻設定部と、上記生体情
報値の時系列の進行方向の時間変化の傾向を表すように
設定した前トレンド曲線を算出する前トレンド曲線算出
部と、上記生体情報値の時系列の逆行方向の時間変化の
傾向を表すように設定した後トレンド曲線を算出する後
トレンド曲線算出部と、上記前トレンド曲線と上記後ト
レンド曲線のうち各測定時刻において大きいほうの値を
採用してトレンド曲線を算出するトレンド曲線算出部
と、上記生体情報値とトレンド曲線から各測定時刻での
生体情報値の増加度合いを算出する増分算出部と、上記
各測定時刻の前の一定時間内での上記生体情報値のばら
つき度合いを算出するばらつき度算出部と、この算出さ
れた生体情報値の増分とばらつき度の組み合わせから、
予め設定された睡眠状態出現確率密度分布表をもとに各
睡眠状態の出現確率を算出する睡眠状態出現確率算出部
と、上記各睡眠状態の出現確率のなかで最も大きな数値
を採った状態をその測定時刻の睡眠状態と推定する睡眠
状態推定部とを具備するので、前トレンド曲線と後トレ
ンド曲線から算出されるトレンド曲線をもとにして生体
情報値の増加度合いをより生体の状態を正確に反映した
値として算出することができ、この生体情報値の増加度
合いとばらつき度合いとから各睡眠状態の出現確率を算
出することができるものであり、睡眠状態を精度良く、
しかもリアルタイムで検出することができるものであ
る。
According to a fifth aspect of the present invention, there is provided a sleep state judging device for measuring a biological information value such as a heart rate or a pulse rate per unit time to obtain a biological information measured value at each time, A measurement time setting unit that sets a pre-trend curve calculation unit that calculates a pre-trend curve that is set to represent a tendency of a time change in the traveling direction of the time series of the biometric information value, The post-trend curve calculation unit that calculates the trend curve after setting to represent the tendency of the time change in the reverse direction, and adopts the larger value at each measurement time among the front trend curve and the rear trend curve. A trend curve calculation unit for calculating a trend curve, an increment calculation unit for calculating the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve, and a constant before the measurement time. A variation degree calculating unit for calculating the degree of dispersion of the biological information value within between, a combination of incremental and variation degree of the calculated biological information value,
A sleep state appearance probability calculation unit that calculates the appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table, and the state that takes the largest numerical value among the appearance probabilities of each sleep state. Since it has a sleep state estimation unit that estimates the sleep state at the measurement time, the degree of increase of the biological information value can be more accurately determined based on the trend curve calculated from the front trend curve and the rear trend curve. It is possible to calculate the appearance probability of each sleep state from the increase degree and the variation degree of the biological information value, and to accurately calculate the sleep state,
Moreover, it can be detected in real time.

【0056】本発明の請求項6に係る睡眠状態判定装置
は、単位時間あたりの心拍数または脈拍数などの生体情
報値を測定して各時刻の生体情報計測値とする計測部
と、測定時刻を設定する測定時刻設定部と、上記生体情
報値の時系列の進行方向の時間変化の傾向を表すように
設定した前トレンド曲線を算出する前トレンド曲線算出
部と、上記生体情報値の時系列の逆行方向の時間変化の
傾向を表すように設定した後トレンド曲線を算出する後
トレンド曲線算出部と、上記前トレンド曲線と上記後ト
レンド曲線のうち各測定時刻において大きなほうの値を
採用してトレンド曲線を算出するトレンド曲線算出部
と、上記生体情報値とトレンド曲線から各測定時刻での
生体情報値の増加度合いを算出する増分算出部と、上記
各測定時刻の前の一定時間内での上記生体情報値のばら
つき度合いを算出するばらつき度算出部と、この算出さ
れた生体情報値の増分とばらつき度の組み合わせから、
予め設定された睡眠状態出現確率密度分布表をもとに各
睡眠状態の出現確率を算出する睡眠状態出現確率算出部
と、動作開始時からの経過時間から、予め設定された睡
眠状態の出現割合表をもとに各睡眠状態の出現割合を算
出する睡眠状態出現割合算出部と、上記各睡眠状態出現
確率と上記各睡眠状態出現割合を掛け合わせて各睡眠状
態出現比率を算出する睡眠状態出現比率算出部と、上記
各睡眠状態出現比率のなかで最も大きな数値を採った状
態をその測定時刻の睡眠状態と推定する睡眠状態推定部
とを具備するので、前トレンド曲線と後トレンド曲線か
ら算出されるトレンド曲線をもとにして生体情報値の増
加度合いをより生体の状態を正確に反映した値として算
出することができ、しかもこの生体情報値の増加度合い
とばらつき度合いに加えて、睡眠状態出現割合表から各
睡眠状態の出現割合を算出することができ、時間経過の
概念を加えることができて、睡眠状態を精度良く、しか
もリアルタイムで検出することができるものである。
According to a sixth aspect of the present invention, there is provided a sleep state judging device for measuring a biological information value such as a heart rate or a pulse rate per unit time to obtain a biological information measured value at each time; A measurement time setting unit that sets a pre-trend curve calculation unit that calculates a pre-trend curve set so as to represent a tendency of a temporal change in the traveling direction of the time series of the biological information value, and a time series of the biological information value A post-trend curve calculation unit that calculates a trend curve after setting it to represent the tendency of the time change in the reverse direction, and adopts the larger value at each measurement time among the above-mentioned front trend curve and the above-mentioned rear trend curve. A trend curve calculation unit for calculating a trend curve, an increment calculation unit for calculating the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve, and a constant before the measurement time. A variation degree calculating unit for calculating the degree of dispersion of the biological information value within between, a combination of incremental and variation degree of the calculated biological information value,
A sleep state appearance probability calculation unit that calculates the appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table, and a predetermined sleep state appearance ratio based on the elapsed time from the start of operation. A sleep state appearance ratio calculating unit that calculates an appearance ratio of each sleep state based on a table, and a sleep state appearance ratio that calculates each sleep state appearance ratio by multiplying each sleep state appearance probability by each sleep state appearance ratio. Since it has a ratio calculation unit and a sleep state estimation unit that estimates the state with the largest numerical value among the sleep state appearance ratios as the sleep state at the measurement time, it is calculated from the front trend curve and the rear trend curve The degree of increase in the biological information value can be calculated as a value that more accurately reflects the state of the living body based on the trend curve that is obtained. In addition, the appearance ratio of each sleep state can be calculated from the sleep state appearance ratio table, and the concept of elapse of time can be added, so that the sleep state can be detected accurately and in real time. .

【0057】また請求項7の発明は、上記睡眠状態推定
部より出力された推定結果に対して、入眠と推定した
後、予め設定された時間の間は深い睡眠状態にならない
ように推定結果を訂正する睡眠判定調整部を具備するの
で、睡眠状態を精度良く検出することができるものであ
る。
According to a seventh aspect of the present invention, the estimated result output from the sleep state estimating unit is estimated such that the user falls asleep, and then the estimated result is set so as not to be in a deep sleep state for a preset time. Since the sleep state adjusting unit for correcting the sleep state is provided, the sleep state can be accurately detected.

【0058】また請求項8の発明は、上記睡眠状態推定
部より出力された推定結果に対して、入眠と推定した
後、ばらつき度が予め設定された閾値未満である状態が
連続して予め設定された時間の間持続した場合に、深い
睡眠状態となるように推定結果を訂正する睡眠判定調整
部を具備するので、睡眠状態を精度良く検出することが
できるものである。
According to the eighth aspect of the present invention, the state in which the degree of variation is less than a preset threshold is continuously determined after the sleep state estimating unit is estimated to fall asleep. Since the sleep determination adjustment unit that corrects the estimation result so as to be in the deep sleep state when the sleep state is continued for the set time is provided, the sleep state can be accurately detected.

【図面の簡単な説明】[Brief description of the drawings]

【図1】請求項1,2の発明の実施の形態の一例におけ
る制御部のブロック図である。
FIG. 1 is a block diagram of a control unit according to an embodiment of the present invention.

【図2】同上のトレンド曲線算出部の動作を示すフロー
チャートである。
FIG. 2 is a flowchart showing the operation of the above trend curve calculation unit.

【図3】同上の一晩の心拍数の推移とそれに対応するト
レンド曲線を示すグラフである。
FIG. 3 is a graph showing a change in heart rate over night and a trend curve corresponding thereto.

【図4】請求項3の発明の実施の形態の一例における制
御部のブロック図である。
FIG. 4 is a block diagram of a control unit according to an embodiment of the third invention.

【図5】同上のトレンド曲線算出部の動作を示すフロー
チャートである。
FIG. 5 is a flowchart showing the operation of the above trend curve calculation unit.

【図6】同上の一晩の心拍数の推移とそれに対応するト
レンド曲線を示すグラフである。
FIG. 6 is a graph showing a change in heart rate over night and a trend curve corresponding thereto.

【図7】同上の一晩の心拍数の増分の推移を示すグラフ
である。
FIG. 7 is a graph showing a transition of an increase in a heart rate in the same night as above.

【図8】同上の一晩の心拍数のばらつき度を示すグラフ
である。
FIG. 8 is a graph showing the degree of variation in the heart rate over the same night.

【図9】同上の睡眠状態出現確率密度分布表を示す図で
ある。
FIG. 9 is a diagram showing a sleep state appearance probability density distribution table according to the first embodiment;

【図10】請求項4の発明の実施の形態の一例における
制御部のブロック図である。
FIG. 10 is a block diagram of a control unit according to an embodiment of the present invention.

【図11】同上の睡眠状態の出現割合表のグラフであ
る。
FIG. 11 is a graph of a sleep state appearance ratio table according to the embodiment.

【図12】請求項5の発明の実施の形態の一例における
制御部のブロック図である。
FIG. 12 is a block diagram of a control unit according to an example of the embodiment of the present invention;

【図13】同上の後トレンド曲線算出部の動作を示すフ
ローチャートである。
FIG. 13 is a flowchart showing the operation of the above post-trend curve calculation unit.

【図14】請求項6の発明の実施の形態の一例における
制御部のブロック図である。
FIG. 14 is a block diagram of a control unit according to an example of the embodiment of the present invention;

【図15】請求項7の発明の実施の形態の一例における
制御部のブロック図である。
FIG. 15 is a block diagram of a control unit according to an embodiment of the present invention;

【図16】請求項8の発明の実施の形態の一例における
制御部のブロック図である。
FIG. 16 is a block diagram of a control unit according to an embodiment of the invention;

【図17】(a)は一晩の睡眠状態の推移の一例を示す
グラフ、(b)は一晩の心拍数の推移の一例を示すグラ
フである。
17A is a graph showing an example of a transition of a sleep state overnight, and FIG. 17B is a graph showing an example of a transition of a heart rate overnight.

【符号の説明】[Explanation of symbols]

1 計測部 2 測定時刻設定部 3 トレンド曲線算出部 4 増分算出部 5 ばらつき度算出部 6 睡眠状態出現確率密度分布表 7 睡眠状態出現確率算出部 8 睡眠状態推定部 9 睡眠状態出現割合表 10 睡眠状態出現割合算出部 11 睡眠状態出現比率算出部 12 前トレンド曲線算出部 13 後トレンド曲線算出部 14 トレンド曲線算出部 15 睡眠判定調整部 16 睡眠判定調整部 DESCRIPTION OF SYMBOLS 1 Measurement part 2 Measurement time setting part 3 Trend curve calculation part 4 Increment calculation part 5 Variation degree calculation part 6 Sleep state appearance probability density distribution table 7 Sleep state appearance probability calculation part 8 Sleep state estimation part 9 Sleep state appearance ratio table 10 Sleep State appearance ratio calculation unit 11 Sleep state appearance ratio calculation unit 12 Front trend curve calculation unit 13 Back trend curve calculation unit 14 Trend curve calculation unit 15 Sleep judgment adjustment unit 16 Sleep judgment adjustment unit

───────────────────────────────────────────────────── フロントページの続き (72)発明者 中野 紀夫 大阪府門真市大字門真1048番地松下電工株 式会社内 (72)発明者 道盛 章弘 大阪府門真市大字門真1048番地松下電工株 式会社内 (72)発明者 萩原 啓 大阪府門真市大字門真1048番地松下電工株 式会社内 Fターム(参考) 4C017 AA02 AA10 AC16 BC11 4C038 PP05 PQ00 PS00  ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Norio Nakano 1048 Kadoma Kadoma, Kadoma City, Osaka Pref. Inside Matsushita Electric Works, Ltd. (72) Inventor Kei Hagiwara 1048 Kazuma Kadoma, Osaka Prefecture F-term in Matsushita Electric Works Co., Ltd. (reference) 4C017 AA02 AA10 AC16 BC11 4C038 PP05 PQ00 PS00

Claims (8)

【特許請求の範囲】[Claims] 【請求項1】 心拍数または脈拍数などの生体情報を計
測して、生体情報の時系列の時間変化の増減傾向を表す
トレンド曲線を算出し、このトレンド曲線をもとに睡眠
状態を判定する睡眠状態判定方法において、トレンド曲
線を予め設定された時間の窓関数を用いて算出すること
を特徴とする睡眠状態判定方法。
1. A method of measuring biological information such as a heart rate or a pulse rate, calculating a trend curve representing an increasing / decreasing tendency of a time series change of the biological information, and determining a sleep state based on the trend curve. A sleep state determination method, wherein a trend curve is calculated using a window function for a preset time.
【請求項2】 窓関数を用いてトレンド曲線を算出する
にあたって、窓関数の時間の長さを40分〜70分とす
ることを特徴とする請求項1に記載の睡眠状態判定方
法。
2. The sleep state determination method according to claim 1, wherein, when calculating the trend curve using the window function, the length of time of the window function is set to 40 minutes to 70 minutes.
【請求項3】 単位時間あたりの心拍数または脈拍数な
どの生体情報値を測定して各時刻の生体情報計測値とす
る計測部と、測定時刻を設定する測定時刻設定部と、生
体情報値の時系列の時間変化の傾向を表すように設定し
たトレンド曲線を算出するトレンド曲線算出部と、上記
生体情報値とトレンド曲線から各測定時刻での生体情報
値の増加度合いを算出する増分算出部と、上記各測定時
刻の前の一定時間内での上記生体情報値のばらつき度合
いを算出するばらつき度算出部と、この算出された生体
情報値の増分とばらつき度の組み合わせから、予め設定
された睡眠状態出現確率密度分布表をもとに各睡眠状態
の出現確率を算出する睡眠状態出現確率算出部と、睡眠
状態出現確率算出部で算出された各睡眠状態の出現確率
のなかで最も大きな数値を採った状態をその測定時刻の
睡眠状態と推定する睡眠状態推定部とを具備して成るこ
とを特徴とする睡眠状態判定装置。
3. A measuring unit that measures a biological information value such as a heart rate or a pulse rate per unit time and obtains a biological information measurement value at each time, a measurement time setting unit that sets a measurement time, and a biological information value. A trend curve calculating unit that calculates a trend curve set to represent the tendency of the time change of the time series, and an increment calculating unit that calculates the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve. And a variation calculation unit that calculates the variation of the biological information value within a certain time period before each of the measurement times, and a combination of the increment and the variation of the calculated biological information value, which is set in advance. The sleep state appearance probability calculation unit that calculates the appearance probability of each sleep state based on the sleep state appearance probability density distribution table, and the largest among the sleep state appearance probabilities calculated by the sleep state appearance probability calculation unit A sleep state estimating device comprising: a sleep state estimating unit for estimating a state in which a numerical value is taken as a sleep state at the measurement time.
【請求項4】 単位時間あたりの心拍数または脈拍数な
どの生体情報値を測定して各時刻の生体情報計測値とす
る計測部と、測定時刻を設定する測定時刻設定部と、生
体情報値の時系列の時間変化の傾向を表すように設定し
たトレンド曲線を算出するトレンド曲線算出部と、上記
生体情報値とトレンド曲線から各測定時刻での生体情報
値の増加度合いを算出する増分算出部と、上記各測定時
刻の前の一定時間内での上記生体情報値のばらつき度合
いを算出するばらつき度算出部と、この算出された生体
情報値の増分とばらつき度の組み合わせから、予め設定
された睡眠状態出現確率密度分布表をもとに各睡眠状態
の出現確率を算出する睡眠状態出現確率算出部と、動作
開始時からの経過時間から、予め設定された各睡眠状態
の出現割合表をもとに各睡眠状態の出現割合を算出する
睡眠状態出現割合算出部と、上記各睡眠状態出現確率と
上記各睡眠状態出現割合を掛け合わせて各睡眠状態出現
比率を算出する睡眠状態出現比率算出部と、上記各睡眠
状態出現比率のなかで最も大きな数値を採った状態をそ
の測定時刻の睡眠状態と推定する睡眠状態推定部とを具
備して成ることを特徴とする睡眠状態判定装置。
4. A measuring unit that measures a biological information value such as a heart rate or a pulse rate per unit time to obtain a biological information measurement value at each time, a measurement time setting unit that sets a measurement time, and a biological information value. A trend curve calculating unit that calculates a trend curve set to represent the tendency of the time change of the time series, and an increment calculating unit that calculates the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve. And a variation calculation unit that calculates the variation of the biological information value within a certain time period before each of the measurement times, and a combination of the increment and the variation of the calculated biological information value, which is set in advance. A sleep state appearance probability calculation unit that calculates the appearance probability of each sleep state based on the sleep state appearance probability density distribution table, and an appearance ratio table of each sleep state that is set in advance from the elapsed time from the start of the operation. When A sleep state appearance ratio calculation unit that calculates the appearance ratio of each sleep state, and a sleep state appearance ratio calculation unit that calculates each sleep state appearance ratio by multiplying each sleep state appearance probability and each sleep state appearance ratio. A sleep state estimating unit for estimating a state having the largest numerical value among the sleep state appearance ratios as a sleep state at the measurement time.
【請求項5】 単位時間あたりの心拍数または脈拍数な
どの生体情報値を測定して各時刻の生体情報計測値とす
る計測部と、測定時刻を設定する測定時刻設定部と、上
記生体情報値の時系列の進行方向の時間変化の傾向を表
すように設定した前トレンド曲線を算出する前トレンド
曲線算出部と、上記生体情報値の時系列の逆行方向の時
間変化の傾向を表すように設定した後トレンド曲線を算
出する後トレンド曲線算出部と、上記前トレンド曲線と
上記後トレンド曲線のうち各測定時刻において大きいほ
うの値を採用してトレンド曲線を算出するトレンド曲線
算出部と、上記生体情報値とトレンド曲線から各測定時
刻での生体情報値の増加度合いを算出する増分算出部
と、上記各測定時刻の前の一定時間内での上記生体情報
値のばらつき度合いを算出するばらつき度算出部と、こ
の算出された生体情報値の増分とばらつき度の組み合わ
せから、予め設定された睡眠状態出現確率密度分布表を
もとに各睡眠状態の出現確率を算出する睡眠状態出現確
率算出部と、上記各睡眠状態の出現確率のなかで最も大
きな数値を採った状態をその測定時刻の睡眠状態と推定
する睡眠状態推定部とを具備して成ることを特徴とする
睡眠状態判定装置。
5. A measuring unit for measuring a biological information value such as a heart rate or a pulse rate per unit time to obtain a biological information measurement value at each time, a measurement time setting unit for setting a measurement time, and the biological information A pre-trend curve calculation unit that calculates a pre-trend curve set to represent a time-varying trend in a time-series traveling direction of a value, and a time-series trend of the biometric information value in a retrograde direction. A post-trend curve calculating unit that calculates a post-trend curve after setting, a trend curve calculating unit that calculates a trend curve by adopting a larger value at each measurement time among the pre-trend curve and the post-trend curve, An increment calculation unit that calculates the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve, and calculates the degree of variation of the biological information value within a certain time before each of the measurement times. A sleep state calculating unit that calculates an appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table from a combination of the calculated increment of the biological information value and the degree of dispersion, A sleep state comprising: an appearance probability calculating unit; and a sleep state estimating unit for estimating a state having the largest numerical value among the appearance probabilities of the respective sleep states as a sleep state at the measurement time. Judgment device.
【請求項6】 単位時間あたりの心拍数または脈拍数な
どの生体情報値を測定して各時刻の生体情報計測値とす
る計測部と、測定時刻を設定する測定時刻設定部と、上
記生体情報値の時系列の進行方向の時間変化の傾向を表
すように設定した前トレンド曲線を算出する前トレンド
曲線算出部と、上記生体情報値の時系列の逆行方向の時
間変化の傾向を表すように設定した後トレンド曲線を算
出する後トレンド曲線算出部と、上記前トレンド曲線と
上記後トレンド曲線のうち各測定時刻において大きなほ
うの値を採用してトレンド曲線を算出するトレンド曲線
算出部と、上記生体情報値とトレンド曲線から各測定時
刻での生体情報値の増加度合いを算出する増分算出部
と、上記各測定時刻の前の一定時間内での上記生体情報
値のばらつき度合いを算出するばらつき度算出部と、こ
の算出された生体情報値の増分とばらつき度の組み合わ
せから、予め設定された睡眠状態出現確率密度分布表を
もとに各睡眠状態の出現確率を算出する睡眠状態出現確
率算出部と、動作開始時からの経過時間から、予め設定
された睡眠状態の出現割合表をもとに各睡眠状態の出現
割合を算出する睡眠状態出現割合算出部と、上記各睡眠
状態出現確率と上記各睡眠状態出現割合を掛け合わせて
各睡眠状態出現比率を算出する睡眠状態出現比率算出部
と、上記各睡眠状態出現比率のなかで最も大きな数値を
採った状態をその測定時刻の睡眠状態と推定する睡眠状
態推定部とを具備して成ることを特徴とする睡眠状態判
定装置。
6. A measuring unit for measuring a biological information value such as a heart rate or a pulse rate per unit time to obtain a biological information measurement value at each time, a measurement time setting unit for setting a measurement time, and the biological information A pre-trend curve calculation unit that calculates a pre-trend curve set to represent a time-varying trend in a time-series traveling direction of a value, and a time-series trend of the biometric information value in a retrograde direction. A post-trend curve calculation unit that calculates a post-trend curve after setting, a trend curve calculation unit that calculates a trend curve by adopting a larger value at each measurement time among the pre-trend curve and the post-trend curve, An increment calculation unit that calculates the degree of increase of the biological information value at each measurement time from the biological information value and the trend curve, and calculates the degree of variation of the biological information value within a certain time before each of the measurement times. A sleep state calculating unit that calculates an appearance probability of each sleep state based on a preset sleep state appearance probability density distribution table from a combination of the calculated increment of the biological information value and the degree of dispersion, An appearance probability calculation unit, a sleep state appearance ratio calculation unit that calculates an appearance ratio of each sleep state based on a preset sleep state appearance ratio table from an elapsed time from the start of the operation, and the sleep state A sleep state appearance ratio calculation unit that calculates each sleep state appearance ratio by multiplying the appearance probability and each of the sleep state appearance ratios, and the state that takes the largest numerical value among the sleep state appearance ratios as the measurement time A sleep state determination device, comprising: a sleep state estimation unit that estimates a sleep state.
【請求項7】 上記睡眠状態推定部より出力された推定
結果に対して、入眠と推定した後、予め設定された時間
の間は深い睡眠状態にならないように推定結果を訂正す
る睡眠判定調整部を具備して成ることを特徴とする請求
項3〜6のいずれかに記載の睡眠状態判定装置。
7. A sleep determination adjustment unit that corrects the estimation result output from the sleep state estimation unit after estimating that the user has fallen asleep, so as not to be in a deep sleep state for a preset time. The sleep state judging device according to any one of claims 3 to 6, further comprising:
【請求項8】 上記睡眠状態推定部より出力された推定
結果に対して、入眠と推定した後、ばらつき度が予め設
定された閾値未満である状態が連続して予め設定された
時間の間持続した場合に、深い睡眠状態となるように推
定結果を訂正する睡眠判定調整部を具備して成ることを
特徴とする請求項3〜7のいずれかに記載の睡眠状態判
定装置。
8. A state in which the degree of variability is less than a preset threshold for the estimation result output from the sleep state estimating unit, and the state where the degree of variation is less than a preset threshold continuously continues for a preset time. The sleep state determination device according to any one of claims 3 to 7, further comprising a sleep determination adjustment unit that corrects an estimation result so as to be in a deep sleep state.
JP24053299A 1999-08-26 1999-08-26 Sleep state determination device Expired - Fee Related JP3788126B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP24053299A JP3788126B2 (en) 1999-08-26 1999-08-26 Sleep state determination device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP24053299A JP3788126B2 (en) 1999-08-26 1999-08-26 Sleep state determination device

Publications (2)

Publication Number Publication Date
JP2001061820A true JP2001061820A (en) 2001-03-13
JP3788126B2 JP3788126B2 (en) 2006-06-21

Family

ID=17060941

Family Applications (1)

Application Number Title Priority Date Filing Date
JP24053299A Expired - Fee Related JP3788126B2 (en) 1999-08-26 1999-08-26 Sleep state determination device

Country Status (1)

Country Link
JP (1) JP3788126B2 (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001061797A (en) * 1999-08-26 2001-03-13 Matsushita Electric Works Ltd Fall-asleep judgment device
JP2003052649A (en) * 2001-08-10 2003-02-25 Kawasaki Heavy Ind Ltd Care-receiver monitor
JP2003260040A (en) * 2002-03-11 2003-09-16 Sanyo Electric Co Ltd Sleeping depth estimation device
WO2004107978A1 (en) * 2003-06-03 2004-12-16 Cb System Co. Sleep stage judgment method and judgment device
JP2005198829A (en) * 2004-01-15 2005-07-28 Seiko Epson Corp Biological data analyzer, biological data analyzing method, control program and recording medium
WO2005082252A1 (en) * 2004-03-01 2005-09-09 Cb System Co. Method for judging stage of sleep
JP2006508745A (en) * 2002-12-04 2006-03-16 カーディアック ペースメーカーズ,インコーポレイテッド Sleep detection with adjustable threshold
JP2006192152A (en) * 2005-01-14 2006-07-27 Toshiba Corp Sleep state determination device, method and program
JP2012065853A (en) * 2010-09-24 2012-04-05 Sleep System Kenkyusho:Kk Sleep level determining device and sleep level determining method
JP2013198720A (en) * 2012-02-23 2013-10-03 Tanita Corp Apparatus and method for measuring blood pressure
JP2014516681A (en) * 2011-05-18 2014-07-17 ヴィー−ウォッチ ソシエテ・アノニム System and apparatus for determining a person's sleep and sleep stage
JP2017169884A (en) * 2016-03-24 2017-09-28 トヨタ自動車株式会社 Sleep state estimation apparatus
JP2019187681A (en) * 2018-04-23 2019-10-31 セイコーエプソン株式会社 Pulse analysis device, and pulse analysis method
CN111248857A (en) * 2019-09-19 2020-06-09 听枫智能科技(深圳)有限公司 Sleep quality monitoring method and device, computer equipment and storage medium
JPWO2019151195A1 (en) * 2018-01-30 2021-01-28 京セラ株式会社 Electronics, estimation systems, control methods and control programs
CN116035536A (en) * 2023-03-14 2023-05-02 安徽星辰智跃科技有限责任公司 Method, system and device for detecting and quantifying sleep activity level
US11653876B2 (en) 2019-09-19 2023-05-23 Casio Computer Co., Ltd. Cyclic alternative pattern (CAP) detection device, cyclic alternative pattern (CAP) detection method, and recording medium
CN116392085A (en) * 2023-06-06 2023-07-07 安徽星辰智跃科技有限责任公司 Sleep stability quantification and adjustment method, system and device based on trend analysis
CN116392087A (en) * 2023-06-06 2023-07-07 安徽星辰智跃科技有限责任公司 Sleep stability quantification and adjustment method, system and device based on modal decomposition
CN116509336A (en) * 2023-06-27 2023-08-01 安徽星辰智跃科技有限责任公司 Sleep periodicity detection and adjustment method, system and device based on waveform analysis
CN116509336B (en) * 2023-06-27 2024-05-03 安徽星辰智跃科技有限责任公司 Sleep periodicity detection and adjustment method, system and device based on waveform analysis

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0341926A (en) * 1989-07-07 1991-02-22 Matsushita Electric Works Ltd Detector for change in sleeping state and sleeping state controller
JPH05159173A (en) * 1991-12-04 1993-06-25 Matsushita Electric Ind Co Ltd Sleeping device
JPH06217946A (en) * 1993-01-26 1994-08-09 Matsushita Electric Works Ltd Biorhythm curve measuring apparatus
JPH06285031A (en) * 1993-03-31 1994-10-11 Nec San-Ei Instr Co Ltd Sleep analyzing system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0341926A (en) * 1989-07-07 1991-02-22 Matsushita Electric Works Ltd Detector for change in sleeping state and sleeping state controller
JPH05159173A (en) * 1991-12-04 1993-06-25 Matsushita Electric Ind Co Ltd Sleeping device
JPH06217946A (en) * 1993-01-26 1994-08-09 Matsushita Electric Works Ltd Biorhythm curve measuring apparatus
JPH06285031A (en) * 1993-03-31 1994-10-11 Nec San-Ei Instr Co Ltd Sleep analyzing system

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001061797A (en) * 1999-08-26 2001-03-13 Matsushita Electric Works Ltd Fall-asleep judgment device
JP2003052649A (en) * 2001-08-10 2003-02-25 Kawasaki Heavy Ind Ltd Care-receiver monitor
JP2003260040A (en) * 2002-03-11 2003-09-16 Sanyo Electric Co Ltd Sleeping depth estimation device
JP2006508745A (en) * 2002-12-04 2006-03-16 カーディアック ペースメーカーズ,インコーポレイテッド Sleep detection with adjustable threshold
WO2004107978A1 (en) * 2003-06-03 2004-12-16 Cb System Co. Sleep stage judgment method and judgment device
JP2005198829A (en) * 2004-01-15 2005-07-28 Seiko Epson Corp Biological data analyzer, biological data analyzing method, control program and recording medium
WO2005082252A1 (en) * 2004-03-01 2005-09-09 Cb System Co. Method for judging stage of sleep
US7699785B2 (en) 2004-03-01 2010-04-20 Sleep System Laboratory Inc. Method for determining sleep stages
JP2006192152A (en) * 2005-01-14 2006-07-27 Toshiba Corp Sleep state determination device, method and program
JP2012065853A (en) * 2010-09-24 2012-04-05 Sleep System Kenkyusho:Kk Sleep level determining device and sleep level determining method
JP2014516681A (en) * 2011-05-18 2014-07-17 ヴィー−ウォッチ ソシエテ・アノニム System and apparatus for determining a person's sleep and sleep stage
JP2013198720A (en) * 2012-02-23 2013-10-03 Tanita Corp Apparatus and method for measuring blood pressure
JP2017169884A (en) * 2016-03-24 2017-09-28 トヨタ自動車株式会社 Sleep state estimation apparatus
JPWO2019151195A1 (en) * 2018-01-30 2021-01-28 京セラ株式会社 Electronics, estimation systems, control methods and control programs
JP2019187681A (en) * 2018-04-23 2019-10-31 セイコーエプソン株式会社 Pulse analysis device, and pulse analysis method
JP7110701B2 (en) 2018-04-23 2022-08-02 セイコーエプソン株式会社 Pulse analysis device and pulse analysis method
CN111248857A (en) * 2019-09-19 2020-06-09 听枫智能科技(深圳)有限公司 Sleep quality monitoring method and device, computer equipment and storage medium
US11653876B2 (en) 2019-09-19 2023-05-23 Casio Computer Co., Ltd. Cyclic alternative pattern (CAP) detection device, cyclic alternative pattern (CAP) detection method, and recording medium
CN116035536A (en) * 2023-03-14 2023-05-02 安徽星辰智跃科技有限责任公司 Method, system and device for detecting and quantifying sleep activity level
CN116392085A (en) * 2023-06-06 2023-07-07 安徽星辰智跃科技有限责任公司 Sleep stability quantification and adjustment method, system and device based on trend analysis
CN116392087A (en) * 2023-06-06 2023-07-07 安徽星辰智跃科技有限责任公司 Sleep stability quantification and adjustment method, system and device based on modal decomposition
CN116392087B (en) * 2023-06-06 2023-09-01 安徽星辰智跃科技有限责任公司 Sleep stability quantification and adjustment method, system and device based on modal decomposition
CN116392085B (en) * 2023-06-06 2023-09-12 安徽星辰智跃科技有限责任公司 Sleep stability quantification and adjustment method, system and device based on trend analysis
CN116509336A (en) * 2023-06-27 2023-08-01 安徽星辰智跃科技有限责任公司 Sleep periodicity detection and adjustment method, system and device based on waveform analysis
CN116509336B (en) * 2023-06-27 2024-05-03 安徽星辰智跃科技有限责任公司 Sleep periodicity detection and adjustment method, system and device based on waveform analysis

Also Published As

Publication number Publication date
JP3788126B2 (en) 2006-06-21

Similar Documents

Publication Publication Date Title
JP2001061820A (en) Sleeping state judging method and sleeping state judging device
JP4582642B2 (en) Sleep stage determination device
JP3733133B2 (en) Sleep state estimation device
US7429247B2 (en) Sleep state estimating device and program product
EP2914171B1 (en) Measuring psychological stress from cardiovascular and activity signals
JPH0556902B2 (en)
JP5740006B2 (en) Respiration measurement system and REM sleep determination system
US20080234785A1 (en) Sleep controlling apparatus and method, and computer program product thereof
JP5632986B2 (en) Sleep stage determination device and sleep stage determination method
WO2005082252A1 (en) Method for judging stage of sleep
WO2013161072A1 (en) Device for calculating amount of retained physical activity, method for calculating amount of retained physical activity and system for calculating amount of retained physical activity
JP6079824B2 (en) Sleep evaluation apparatus and program
CN101095612A (en) Apparatus and method for monitoring biological information
JP5776939B2 (en) Sleep state determination device and sleep state determination method
JP5929020B2 (en) Consciousness state estimation device and program
JP6040874B2 (en) Sleep stage estimation device
JP5857609B2 (en) Sleep state discrimination device
JP2012055464A (en) Sleep evaluating device, sleep evaluating system, and program
JP2016016144A (en) Biological information processing system and method of controlling biological information processing system
JP2004089267A (en) Sleeping depth estimation device and bedding equipped with the same
JP2006271474A (en) Apparatus for estimating sleeping state
US20050137470A1 (en) Method and apparatus for low blood glucose level detection
JP6845404B2 (en) Sleep stage determination method, sleep stage determination device, and sleep stage determination program
JPH07143972A (en) Method and apparatus for judging sleeping condition
US9968293B1 (en) Detecting and estimating sleep stages

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050809

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051011

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20051122

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060123

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: 20060307

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060320

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090407

Year of fee payment: 3

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090407

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100407

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100407

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110407

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130407

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees