WO2009150765A1 - Sleeping condition monitoring apparatus, monitoring system, and computer program - Google Patents

Sleeping condition monitoring apparatus, monitoring system, and computer program Download PDF

Info

Publication number
WO2009150765A1
WO2009150765A1 PCT/JP2008/072697 JP2008072697W WO2009150765A1 WO 2009150765 A1 WO2009150765 A1 WO 2009150765A1 JP 2008072697 W JP2008072697 W JP 2008072697W WO 2009150765 A1 WO2009150765 A1 WO 2009150765A1
Authority
WO
WIPO (PCT)
Prior art keywords
interval
sleep
beat
unit
waveform
Prior art date
Application number
PCT/JP2008/072697
Other languages
French (fr)
Japanese (ja)
Inventor
博明 鈴木
和義 坂本
Original Assignee
ハートメトリクス株式会社
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 ハートメトリクス株式会社 filed Critical ハートメトリクス株式会社
Priority to JP2009141394A priority Critical patent/JP2011115188A/en
Publication of WO2009150765A1 publication Critical patent/WO2009150765A1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/113Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0816Measuring devices for examining respiratory frequency
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4806Sleep evaluation
    • A61B5/4812Detecting sleep stages or cycles
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6887Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
    • A61B5/6892Mats
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/40Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4029Detecting, measuring or recording for evaluating the nervous system for evaluating the peripheral nervous systems
    • A61B5/4035Evaluating the autonomic nervous system

Definitions

  • the present invention relates to a sleep state monitoring device, a monitoring system, and a computer program.
  • sleep is classified into six categories, namely, awakening, REM (Rapid Eye Movements) sleep, and stage 1 (S1) to stage 4 (S4).
  • REM Rapid Eye Movements
  • S1 to S4 stage 4
  • electrical measurement is usually used.
  • awakening, REM sleep, and S1 to S4 increase the depth of sleep in this order.
  • REM sleep does not occur immediately after the start of sleep, and deep sleep (S3 or When the sleep becomes shallower toward awakening after S4), the REM sleep state is entered.
  • electrodes are attached to several parts of a person's (user or subject) body, for example, an electroencephalogram from the scalp, a myoelectric potential from the jaw, an ocular potential from the periphery of the eye, And you need to use ECG.
  • a technique for determining a sleep state from an electroencephalogram is known as the Rechtschaffen-Kales (RK) method.
  • a technique using an electrocardiogram is disclosed in, for example, Patent Document 1 (Japanese Patent Laid-Open No. 2003-225211).
  • Non-Patent Document 1 a method of defining the body vibration magnitude as an index
  • Non-Patent Document 3 discloses that an autonomic component ratio is obtained based on an electrocardiogram waveform and a determination is made between shallow sleep (S1 or S2) and deep sleep (S3 or S4).
  • Non-Patent Document 5 discloses changes in the spectrum of heartbeats in heart disease
  • Non-Patent Document 6 applies the power spectrum of heartbeat variability to provide an autonomic nervous function at an early stage due to diabetes. It is disclosed that a failure can be detected.
  • a device or method that allows the user to easily check the quality of sleep while staying at home without the assistance of a specialist or health care professional makes it possible to grasp his / her health condition and manage his / her health. Becomes easier. Moreover, if a decrease in the quality of sleep can be detected easily, accurate information for diagnosis by a specialist or a healthcare professional can be given.
  • An object of the present invention is to solve the above several problems.
  • a pressure waveform acquisition unit that obtains a pressure waveform as an electric signal by a pressure detection unit that is in direct contact with any part of the user's body or through clothing, and a predetermined process is performed on the pressure waveform.
  • a pulsation extraction unit that extracts a pulsation waveform of the user and calculates a pulsation interval that is a time interval from the previous pulsation for each pulsation in the pulsation waveform;
  • a sleep state monitoring device including a sleep stage determination unit that determines a sleep stage that is a stage of the user's sleep depth from the beat interval.
  • the sleep state monitoring of the present invention can also be implemented as a method or as a computer program in a computer connected to the pressure detection unit.
  • the sleep state is used to indicate all states of awakening, REM sleep, awakening and sleep of stage 1 (S1) to stage 4 (S4), and the sleep stage refers to awakening from a sleeping state.
  • the exclusion ie, REM sleep, S1-S4.
  • the pulsation includes the pulsation of the heart and vibrations generated by the pulsation, for example, vibrations of the body surface due to arterial pulsations.
  • a pressure waveform acquisition unit that obtains a pressure waveform as an electric signal by a pressure detection unit that is in direct contact with any part of the user's body or through clothing, and a predetermined process is performed on the pressure waveform.
  • a pulsation extraction unit that extracts a pulsation waveform of a user, a pulsation interval calculation unit that calculates a pulsation interval that is a time interval from the previous pulsation for each pulsation in the pulsation waveform, and a computer
  • a monitoring unit connected to a network including a data transmission unit for transmitting data through the network; and receiving time series data including at least a pulsation interval from the device through the computer network a plurality of times.
  • a sleep state monitoring system having a storage unit for recording the time series data corresponding to a user, the computer network comprising: Any computer connected to the work displays the time-series data received multiple times in association with the user as a graph with respect to time, and changes in fluctuation of the user's pulsation interval in a past period
  • a sleep state monitoring system including a data display unit that presents This sleep state monitoring system can also be implemented as a program in a computer that can be connected to the monitoring unit through a computer network.
  • any embodiment of the present invention which sleep is awake or sleeping, or which user is sleeping, without restraining the user and affecting the sleep state It becomes possible to monitor a sleep state such as whether it is a stage. In addition, according to any embodiment of the present invention, it is possible to monitor the health condition of the user based on fluctuations in the pulsation interval.
  • the flowchart which shows embodiment of the sleep state monitoring method of this invention. The flowchart which shows the process which determines whether it is sleep or awakening in embodiment of the sleep state monitoring method of this invention.
  • FIG. 9 is an example of a computer screen showing a result of obtaining a power spectrum at an RR interval from body vibration obtained from a subject in an example of the present invention, and is an arousal state (FIG. 9A) and S1 to S4 (FIG. 9B).
  • FIG. 9 (E)) and the power spectrum in the case of REM sleep FIG. 9 (F)).
  • the example of the computer screen which shows the result of having determined the sleep state by the body vibration obtained from the test subject in the Example of this invention.
  • the example of the computer screen which compares and shows the result of having determined the sleep state automatically by the body vibration obtained from the subject in the embodiment of the present invention and the result determined by the conventional method.
  • the graph which shows the fluctuation (SDNN) of RR interval measured from the test subject in the Example of this invention.
  • the sleep state is determined from the body vibration waveform.
  • the body vibration waveform is acquired by the sleep state monitoring apparatus 1 including the pressure waveform acquisition unit 10 and the computer 20.
  • the pressure waveform acquisition unit 10 includes an air pressure sensor 12, a transducer unit 14 that converts pressure into an electrical signal, and an A / D converter 16 that performs analog / digital conversion of the electrical signal. These are collectively referred to as a pressure waveform acquisition unit 10.
  • the air pressure sensor 12 is disposed on the bed 100, for example, and is made so that a user can lie on the air pressure sensor 12.
  • a sponge-like stuffing is enclosed in a highly air-tight bag-like envelope, and the overall thickness is about 1 cm to 3 cm, and the cushion shape is a square shape of 30 cm square. It is designed to be flexible and soft enough that you can sleep even when you lie down.
  • a tube that can communicate with internal air and propagate pressure is connected to one peripheral portion of the envelope, and the tube extends to the transducer unit 14.
  • the transducer section 14 has an airtight chamber to which the tube is connected.
  • An opening is provided on one side of the airtight chamber, and a ceramic piezoelectric element is disposed so as to close the opening. .
  • a pressure wave inside the air pressure sensor 12 is transmitted to the piezoelectric element to generate an analog value piezoelectric signal, which becomes a pressure wave signal and is amplified by an appropriate amplifier circuit.
  • the amplified pressure wave signal is input to the A / D converter 16, and a digitized pressure wave signal is output.
  • the digitized pressure wave signal is input to the computer 20.
  • the computer 20 operates as a set of functional blocks that perform predetermined processing in accordance with programmed instructions. That is, the computer 20 includes a pulsation extraction unit 22, a pulsation interval calculation unit 24, a sleep stage determination unit 26, a respiration waveform extraction unit 32, a respiration interval calculation unit 34, an awakening determination unit 36, and sleep.
  • the result is presented as the user's sleep state by a technique such as outputting to the monitor 44 and visually indicating the result.
  • the sleep state monitoring apparatus according to the embodiment of the present invention is configured by the pressure waveform acquisition unit and each functional block realized by the computer 20.
  • data can be accumulated by appropriate storage means provided in the computer 20, and the start, stop, end, interruption, etc. of the process can be controlled by the input / output means. .
  • the sleep state monitoring apparatus of the present invention can determine a sleep state consisting of determination of sleep or awakening and determination of a sleep stage. Specifically, as shown in FIG. 2, in order to determine the sleep state, body vibration data is first acquired (S10). And it is judged whether it is sleep or awakening (S20). If it is sleep (S30, YES), the sleep stage is further determined (S40). After that, or when it is not sleep (when it is awakening, S30, NO), a determination result is output.
  • the process for determining the sleep state using the body vibration in the embodiment of the present invention is divided into a process for determining whether the person is awake or sleeping and a process for determining the sleep stage.
  • this series of processing is repeated at regular intervals (for example, every minute) to periodically perform the process from data acquisition to sleep state determination, or Only the processing to be acquired can be performed periodically, the data can be stored in an appropriate storage means, and the determination processing can be performed afterwards.
  • FIG. 3 shows details of the step of determining whether it is sleep or awakening (FIG. 2, S20), and FIGS. 4 to 6 show details of the step of determining the sleep stage (S40).
  • the present invention can be implemented by a form in which a respiratory waveform is extracted from a body vibration waveform and used.
  • the range of the respiration waveform is naturally determined from the respiration rate of the person at rest. For example, a frequency domain component of 0.13 to 0.70 Hz in body vibration can be used as the respiration waveform. For this reason, the body vibration waveform is filtered and a frequency domain component of 0.13 to 0.70 Hz is extracted (FIG. 3, S202). From the respiration waveform, the fluctuation of the respiration interval is extracted.
  • the time when the respiration waveform shows the maximum value is calculated, and that time is set as the time representative of the respiration timing (respiration peak time) (S204).
  • respiration interval data is generated (S206).
  • a coefficient of variation (CV) for the breathing interval is obtained (S208), and this is used as an index for the uniformity of the breathing interval.
  • CV coefficient of variation
  • the CV of the breathing interval is set as a threshold value of about 0.2 (S210), and when it is larger than that, it can be determined as an awake state (S212), and when smaller, it can be determined as a sleep state (S214). .
  • the process of determining the sleep stage includes the step of identifying the timing of the pulsation by the R wave of the pulsation due to the heartbeat or pulsation (S402), and the data of the time interval (RR interval) between a certain heartbeat and the immediately preceding heartbeat The process is divided into the step (S404) of obtaining the power spectrum, the step of calculating the power spectrum of the RR interval (S406), and the step of determining the sleep stage (S42).
  • R wave timing S402
  • ECG electrocardiogram
  • the description of the R wave uses an expression of the R wave in the ECG for easy understanding, and is not limited to the one that directly corresponds to the EC G R wave.
  • any characteristic waveform can be used as long as it is determined as a timing representative of each beat. is there. For this purpose, it is necessary to extract a wave corresponding to a pulsation from the body vibration waveform (FIG.
  • the body vibration waveform is band-limited to a range of 0.1 Hz to 30 Hz by a band pass filter.
  • the body vibration waveform after the filter is subjected to a second-order differentiation with respect to time, and data of acceleration of body vibration at each time is calculated.
  • the time-dependent data of the acceleration of the body vibration is further band-limited to a range of 3 to 30 Hz by a bandpass filter, and the time when the waveform data shows the peak value is specified.
  • the pulsation extraction unit 22 may include a differential peak processing unit 22A for that purpose.
  • the present invention corresponds to the time of the vibration due to the ventricular excitement of the heart beat, that is, the time of the R wave in the electrocardiogram Guesses.
  • the time may be affected not only by vibration due to the heartbeat but also by vibration (pulsation) caused by the artery due to the heart.
  • RR interval data Time series data of the RR interval (time interval between a certain beat and the immediately preceding beat) is obtained from the R wave time data thus obtained (FIG. 4, S404; FIG. 1, beat interval). Calculation unit 24).
  • the RR interval is given by subtracting the R wave time data of the immediately preceding beat from the R wave time data of each beat.
  • the data of the RR interval is time data at each time of the unequal interval at which the pulsating R wave is detected. This data is conveniently converted into time-series data with a constant sampling interval for later processing, and interpolation processing is performed to obtain time-series data with an RR interval. For this, for example, processing such as secondary interpolation can be performed.
  • the sleep depth is determined from the time-series data of the RR interval obtained as described above.
  • the sleep stage is determined (FIG. 4, S42; FIG. 1, sleep stage determination unit 26).
  • the determination of the sleep stage is a process for determining whether the sleeping state is REM sleep or Non-REM sleep (S1 to S4).
  • the present invention can be implemented by two embodiments.
  • One is an embodiment (spectrum pattern method) in which the sleep stage is determined from the power spectrum pattern of the time series data of the RR interval, and the other is an index indicating the function of the autonomic nervous function from the time series data of the RR interval. Is calculated to determine the sleep stage (autonomic nerve component ratio method).
  • the sleep stage determined by the conventional sleep stage determination method is associated with the power spectrum of the RR interval, and a typical spectrum pattern (“ A typical pattern ”) can be obtained.
  • the sleep stage can be determined by determining which typical pattern is similar to the pattern of the power spectrum of the RR interval measured for a subject (hereinafter, “unknown pattern”).
  • n is set to 2 (S424), and a value of 1 / n of the maximum value of the power spectrum in the RR interval is obtained (S426).
  • a frequency that is 1 / n of the maximum value is obtained, and a difference (frequency width) between the frequencies is calculated (S428).
  • This process is repeated by incrementing n (S432) and repeating until a predetermined upper limit N (S430). Since the frequency width at that time is given to each n, the n dependence of the frequency width can be obtained.
  • the frequency width of the unknown pattern is compared with the n dependence of the frequency width obtained in the same manner for the typical pattern, and the difference between the frequency width of the typical pattern and the unknown pattern is calculated as a residual (S434).
  • the coincidence degree of the unknown pattern with respect to the typical pattern of each sleep stage can be obtained as a numerical value.
  • the peaks of the LF range of the unknown pattern are matched so as to match the peaks of the pattern of the LF range of the typical pattern as much as possible.
  • the peaks of the HF range of the unknown pattern are matched with the peaks of the pattern of the HF range of the typical pattern.
  • no matching is performed between different frequency ranges.
  • a sleep stage that becomes a typical pattern showing the highest degree of coincidence (small residual) with respect to the unknown pattern is set as the sleep stage of the unknown pattern (S436).
  • the frequency ranges of the LF range and the HF range described here are typical, including the description in the autonomic nerve component ratio method described later, and are not limited to the above-described frequency range strictly. For example, Setting the LF range to 0.05 to 0.15 Hz is also part of the disclosure of the present invention.
  • the spectral pattern method can determine the sleep stage, but cannot determine whether the sleep stage is awake or wakefulness, that is, the REM sleep or the awakening.
  • the inventor of the present application has a VLF range (Very Low Frequency: 0.01 Hz to 0.04 Hz) in addition to the above LF range and HF range for the power spectrum of the RR interval.
  • VLF range Very Low Frequency: 0.01 Hz to 0.04 Hz
  • the spectrum in the VLF range is a spectrum that appears at a lower frequency than the LF range in the power spectrum of the RR interval.
  • the time series data of the RR interval needs to have data corresponding to the corresponding period, and the reciprocal of the lower limit frequency of the necessary spectrum range needs to be the window period.
  • this window period in order to correctly obtain a power spectrum of 0.01 Hz, it is effective to set this window period to an appropriate time of about 100 seconds. When actually measuring, for example, it is set to 2 minutes or more. If the spectrum in the VLF range is not used, a shorter time, for example, 60 seconds can be used.
  • the present invention can be implemented by any pattern determination method.
  • a learning process is performed on a neural network corresponding to a result (one of sleep stages) to be determined for a spectrum pattern, and the neural network is applied to an unknown pattern.
  • a so-called supervised neural network learning system for determination can be used. Any other supervised learning technique can be implemented.
  • the sleep stage can also be determined by the autonomic nerve component ratio indicating the function of the autonomic nerve function.
  • an LF component and an HF component are obtained in the power spectrum of the RR interval.
  • the ratio between the LF component and the HF component that is, the LF / HF ratio (autonomic nerve component ratio) is calculated. This is because the function of the subject's autonomic nerve function can be evaluated by the LF / HF ratio (see Non-Patent Documents 2 to 4). In this analysis, as shown in FIG.
  • the power spectrum of the LF range is calculated to obtain the integral value of the range to obtain the LF component (S442), and the power spectrum of the HF range is calculated to integrate the range.
  • a value is obtained and set as an HF component (S444).
  • an LF / HF ratio obtained by dividing the LF component by the HF component is obtained (S446).
  • the sleep stage is determined by comparing the typical LF / HF ratio in each sleep stage with this LF / HF ratio (S448). The smaller the LF / HF ratio is, the deeper the sleep is, the deeper the sleep is, and the lighter is the sleep. However, it is not possible to determine whether the sleep is wake or wake.
  • each sleep stage of S1 to S4 and the LF / HF ratio for example, the correspondence relationship with the determination result of the sleep stage by the spectral pattern method is measured in advance. Adapting to each user is also effective.
  • the inventors of the present application have found that the data acquired at a time when the HF component does not reach a certain value has few fluctuations of high frequency in the HF component, that is, the RR interval, and an error is likely to occur in the sleep stage determination. Yes. An error hardly occurs in the determination of the sleep stage when the value of the HF component is, for example, 300 msec 2 or more. For this reason, it is possible to estimate the superiority or inferiority of the accuracy of the sleep stage determination by determining whether the HF component is larger than a predetermined value before determining the sleep stage. In particular, setting a lower limit for the value of such an HF component and using it in advance is particularly effective when it is desired to accurately determine whether or not it is in the deep sleep stage S3 or S4.
  • a sleep state monitoring system or method is realized that hardly affects the sleep state of the user, can be used at home by the user, and can determine whether it is sleep or awakening. .
  • FIG. 7 shows a configuration for performing this monitoring process.
  • the computer 50 and the pressure waveform acquisition unit 10 are collectively referred to as the monitoring unit 2.
  • the body vibration waveform is acquired by acquiring the pressure waveform acquisition unit 10, that is, the air pressure sensor 12, the transducer unit 14 that converts pressure into an electrical signal, and the A / A that converts the electrical signal into analog / digital.
  • a D converter 16 is used.
  • the air pressure sensor 12 can be disposed on the bed 100.
  • the transducer unit 14 and the A / D converter 16 are connected to output a digitized pressure wave signal.
  • the digitized pressure wave signal is input to the computer 50.
  • the computer 50 operates as a set of functional blocks that perform predetermined processing in accordance with programmed instructions. That is, the computer 50 includes a pulsation extraction unit 22, a pulsation interval calculation unit 24, and a data transmission unit 46.
  • the data transmission unit 46 can upload data to the storage unit 70 through a computer network 60 such as the Internet or an intranet.
  • Data associated with the user's ID is uploaded to the storage unit 70 a plurality of times, and the upload is performed a plurality of times by connecting the data received a plurality of times in any computer connected to the computer network 60.
  • Displayed data can be displayed in one graph.
  • a data generation unit 80 for searching and connecting data of the same user is used. For example, it can be displayed on the monitor 86.
  • Uploaded data is arbitrary time series data including beat interval data.
  • time series data of the pulsation interval may be used alone, or when a sleep stage determination unit (not shown) similar to the sleep stage determination unit 26 of FIG. In addition to data, it can be time series data of the sleep stage. If functional blocks (not shown) corresponding to the respiratory waveform extraction unit 32, the respiratory interval calculation unit 34, the awakening determination unit 36, and the sleep state determination unit 42 of FIG. In addition, time-series data including a determination result of arousal or sleep can be used.
  • the data generation unit 80 detects an obvious change that is different from the normal in the intensity of fluctuations found in the data of the beat interval, and generates an alarm signal.
  • a fluctuation detecting unit 82 for outputting can be provided.
  • the fluctuation detection unit 82 is a part of the data generation unit 80. However, if the fluctuation detection unit 82 is connected to the computer network 60 and can call up data for each user from the storage unit 70, the fluctuation detection unit 82 can be arbitrarily connected to the computer network 60. It can be provided in computer equipment.
  • the fluctuation detection unit 82 monitors a change in the pulsation interval in a certain period of time determined according to the purpose (for example, an arbitrary predetermined period such as one month when data for one year is accumulated in the storage unit 70). is doing.
  • the pulsation interval can be, for example, “SDNN” (Standard Deviation of Normal-to-Normal Intervals) data, which is the standard deviation value of the pulsation interval data, and the coefficient of variation. It is arranged so that it becomes a numerical value such as the value of and is monitored in terms of whether the value falls below a certain threshold. Then, for example, an alarm signal is transmitted when the value falls below the threshold value, and is fed back to the user by appropriate means.
  • SDNN Standard Deviation of Normal-to-Normal Intervals
  • an alarm signal is transmitted when the value falls below the threshold value, and is fed back to the user by appropriate means.
  • a method of displaying an alert on an automatically transmitted e-mail or a personalized management Web page
  • the present inventors when the health state is monitored by the value of SDNN, the present inventors have found that it is effective to combine the determination result of the sleep stage. That is, when using the pulsation interval data acquired at the time when the sleep stage of the subject is S1 or S2, the SDNN calculated therefrom may not be a good reflection of the health condition of the subject.
  • the SDNN calculated from the beat interval data acquired at the time when the stage is S3 or S4 is a good reflection of the health condition of the subject.
  • the inventor of the present application has a sleep state of shallow S1 and S2, and the sleep state of the subject is not stable.
  • the brain stage, electromyogram, electrooculogram, and electrocardiogram were collected from subjects and the sleep stage was determined by a conventional method.
  • a plate electrode was attached to each part of the subject's body. Specifically, electrodes are applied to the top of the head for electroencephalograms, to the jaws for electromyograms, to the left and right sides of the eye for electrocardiograms, and to both wrists by ankle contact for electrocardiograms. Attached.
  • a low-frequency amplifier was connected to the plate electrode to amplify the voltage signal, converted to a digital value by an AD converter, and recorded by a data logging computer.
  • the sleep stage was determined by an electroencephalogram, and the method of Rechtschaffen-Kales (RK) described above was followed.
  • the pressure waveform acquisition unit 10 and the computer 20 which are the devices described in FIG. 1 are used, and in the computer, each functional block as shown in FIG.
  • a software module for extracting a heartbeat component and a respiratory component from body vibration, and an R wave from the heartbeat waveform to extract a spectrum wave A software module for analysis and a software module for performing the autonomic nerve component ratio method were used.
  • FIG. 8 shows an example in which R-wave timing is specified from body vibration.
  • the body vibration waveform C8A-10 and the respiration waveform C8A-20 calculated therefrom are also shown.
  • a waveform C8A-2 shown in the upper part of the center line C8A-00 is obtained by detecting a candidate time of the R wave from the body vibration waveform C8A-10. This detection was performed by the pulsation extraction unit 22. At this time, among the peaks P1 to P8 of the waveform C8A-2, the peak P2 ′ is a false detection.
  • processing for removing the erroneous detection peak is performed. This process is performed when the time point when a candidate with an R wave is obtained and the time point when the next candidate is obtained are less than or equal to a preset time (for example, 0.5 seconds), the next candidate is This is processing that is not employed, such as processing that does not employ a waveform that is highly likely to be obtained regardless of pulsation as an R wave. This processing is performed to remove erroneous detection peaks from the candidates, and the time of the R wave is specified as shown below the center line C8A-00. In this figure, the phenomenon of 8 seconds is enlarged and displayed.
  • the generation time was compared between the R wave obtained from the electrocardiogram by the conventional method and the R wave specified from the body vibration data obtained by the embodiment of the present invention.
  • the detection accuracy of the R wave generation time point obtained from the body vibration waveform measured by the air pressure sensor was 87% on average compared to the R wave generation time point obtained from the electrocardiogram.
  • the correlation coefficient between the RR interval determined from the body vibration and the RR interval according to the electrocardiogram for the R-wave time-series data calculated in this way (secondarily interpolated at 0.1 second intervals) is:
  • the average correlation coefficient was 0.87 between 0.80 and 0.95.
  • the waveform C8B-10 appearing in the chart C8B in the middle of FIG. 8 is a respiratory waveform
  • the horizontal axis is a waveform of 11 minutes in time
  • the vertical axis is a pressure waveform (arbitrary unit).
  • the waveform C8B-20 in the chart C8B is the estimated RR interval.
  • the power spectrum C8C-10 shown here is a result in a deep sleep stage, and the frequency range displayed in the graph is in the range of 0.01 to 0.40 Hz, and the power spectrum has a high frequency (HF) range.
  • HF high frequency
  • Spectral pattern method Specifically, the power spectrum obtained by the method described in relation to FIG. 8 is obtained from a subject in an unknown sleep state, and the sleep pattern is determined by applying the spectral pattern method there. did. First, the observation result of the spectrum to be determined by the spectrum pattern method will be described.
  • FIG. 9A to 9E show power spectra in the case of awakening and sleep states S1 to S4.
  • FIG. 9F shows a power spectrum of REM sleep.
  • the frequency range to be used is 0.01 Hz to 0.40 Hz.
  • the power spectrum of the RR interval has irregular fluctuations in the RR interval at the time of awakening in FIG.
  • the power spectrum of the RR interval has only a small component in the LF and HF ranges, and the VLF range of 0.04 Hz or less is the main spectrum.
  • the sleep stage S1 is in a state of transition from the awakening stage to the S2 stage, and a distribution close to the characteristics at the time of awakening, that is, a distribution of small power spectrum values over LF and HF was observed. Since the power spectrum of the awake state may be mixed into the sleep stage S1, it is difficult to distinguish the awake state from S1 only from the actual power spectrum. On the other hand, in the sleep stage S2 shown in FIG. 9C, the spectrum of the HF range appears more clearly than the spectrum of S1 although the spectrum of the LF range exists. Easier than state identification.
  • FIG. 9 is a power spectrum pattern of REM sleep.
  • REM sleep that appears when sleep approaches wakefulness, the pattern is similar to wakefulness and S1, but the spectrum of the LF range is larger and more complex than in S1.
  • REM sleep and awakening can be distinguished by obtaining a coefficient of variation (CV) from time series data of RR intervals.
  • CV of respiration was also considered in the estimation of arousal. This method was effective in determining sleep awakening, mid-wake awakening, and awakening after the end of sleep. The characteristics when the observer directly observes the power spectrum at the RR interval have been described.
  • FIG. 10 shows body vibration waveform C10A-10 and respiratory waveform C10A-20 (top chart C10A) obtained from the air pressure sensor, and detected R wave time series C10B-10 (second chart C10B from the top). Respiration rate C10C-10 and heart rate C10C-20 (third chart C10C from the top), RR interval power spectrum C10D-10 (bottom left chart C10D), and sleep stage determination example C10E-10 and The LF / HF ratio is C10E-20 (the bottom right chart C10E).
  • FIG. 11 compares the determination result of the sleep state, the broken line R0 is a result determined by a conventional method using an electroencephalogram, and the broken line R1 (broken line) is a result of automatic determination by a spectrum pattern method.
  • the line R1 of the automatic determination result is the same as that shown in the chart C10E of FIG. 10, and the vertical axis represents awakening, REM sleep, and S1 to S4 in order from 1 to ⁇ 4. It was obtained from the measurement data of the minute body vibration. The degree of agreement between the conventional method and the method of the example was about 75%.
  • the autonomic component ratio method was used to distinguish between the shallow sleep stages (S1 and S2) and the deep sleep stages (S3 and S4) from the body vibration waveform acquired under the above conditions. At this time, shallow sleep and deep sleep were judged as shallow sleep when the LF / HF ratio was greater than 0.3, and deep sleep was judged when 0.3 or less. In the range of the acquired data, the degree of coincidence between the determination result and the determination result of the sleep stage by the electroencephalogram was 80% or more.
  • the above-mentioned determination standard is an example.
  • the value for each user used for the determination can be separately recorded in a referable database.
  • the HF component value itself reaches a predetermined value as described above, that is, when the sleep stage determination accuracy reaches a certain level, control is performed so that the sleep stage is determined. It is also possible to do.
  • the HF component does not reach a predetermined value for example, 300 msec 2 , that is, when the sleep stage determination accuracy does not reach a certain level, the result of the sleep stage determination is obtained.
  • a flag indicating that the accuracy is not good when outputting, it is possible to attach a flag indicating that the accuracy is not good.
  • the monitoring unit 2 performs an RR interval calculation process for R wave data.
  • the data of the RR interval is uploaded once a day to a data server that manages the storage unit 70.
  • any data communication means can be used for uploading, for example, it can be performed by putting using FTP (File Transfer Protocol).
  • the form of data storage in the storage unit 70 is arbitrary. For example, a computer file in CSV (Comma Separated Values) format for each user ID and each date using a file system folder (directory) of the data server. Multiple data can be stored. As another example, data can be stored in an appropriate database using a database management system or the like.
  • an optical sensor (not shown) is connected to the computer 50, and the timing of measurement start can be controlled by ambient illuminance information. Furthermore, in order to limit the measurement time to two sleep cycles after falling asleep, only data for 3 hours after the start of measurement can be acquired. In addition, in order to be able to collect beat interval data only when sleeping at a certain depth or more, and to perform data processing using the determination result of the sleep stage later, the RR interval at each time The determination result of the sleep stage at each time can be acquired together with the data. Thereby, for example, using the pulsation data when the sleep stages are S3 and S4, it is possible to acquire SDNN data that favorably reflects the health condition of the subject.
  • the fluctuation of the RR interval is calculated from the data of the RR interval included in the CSV file. This fluctuation can be calculated by calculating a standard deviation from the data of the RR interval acquired during a certain period (for example) and calculating it as SDNN, or by calculating a coefficient of variation (CV).
  • the data of each day divided into a plurality of CSV files are concatenated and calculated, and the fluctuation of the RR interval can be displayed following the date.
  • the fluctuation detection unit 82 includes, for example, a short-term fluctuation monitoring module (not shown) that monitors fluctuations in the RR interval during the past month, and a long-term fluctuation that monitors fluctuations in the RR interval over the past year.
  • a monitoring module (not shown) can be implemented.
  • the short-term fluctuation monitoring module and the long-term fluctuation monitoring module monitor how the fluctuation appearing in the fluctuation of the RR interval changes according to each period.
  • RR interval data acquired in the time in a sleep state suitable for monitoring can be extracted from the sleep time.
  • the sleep state monitoring system 6 shown in FIG. 7 includes the sleep stage determination unit 26, the respiratory waveform extraction unit 32, the respiratory interval calculation unit 34, the awakening determination unit 36, and the sleep state determination shown in FIG.
  • the operations of the short-term fluctuation monitoring module and the long-term fluctuation monitoring module can be controlled by outputs from these functional blocks.
  • SDNN was measured with a pulsation monitoring device prepared in this example for one subject with a high blood glucose level. The measurement is performed twice on the day corresponding to the first day when the subject starts exercise therapy and the day corresponding to the 105th day after starting exercise therapy, and the subject takes about 100 minutes of sleep during each day. RR interval data was acquired during the sleep. Then, RR interval data was accumulated for 60 seconds, a standard deviation was calculated from the RR interval data for 60 seconds, and SDNN was obtained every 60 seconds. At this time, in order to determine the sleep stage in accordance with the SDNN data, sleep stage data was also acquired.
  • the sleep stage was determined by the above-described spectrum pattern method using the power spectrum of the RR interval obtained by the window period of 60 seconds in order to meet the above conditions. Since the window period is 60 seconds, the power spectrum in the VLF range is not used.
  • the said test subject continued for 3 days or more as a exercise therapy between the 1st day and the 105th day, and the walking where the walk distance of each week becomes about 40 km.
  • FIGS. 12A and 12B are both SDNN data on the first day and the 105th day, the horizontal axis is the elapsed time after sleep (unit: minutes), and the vertical axis is SDNN (unit: milliseconds). ).
  • FIG. 12A is a graph in which SDNN measured every minute is plotted as it is, and FIG. 12B is obtained at the time when the subject is in the sleep stage S3 or S4 among the measured SDNN. It is the graph which plotted data. Looking at the trend from FIG. 12A alone, it can be said that the SDNN on the first day has a smaller value than the SDNN on the 105th day. On the other hand, for example, there is no significant difference until about 30 minutes after falling asleep. In addition, this relationship is reversed after 90 minutes after falling asleep.
  • FIG. 12 (B) in which only the cases where the sleep stage is S3 and S4 are plotted, the data greatly deviated in each data series is excluded from the data of any day, Also, data immediately after falling asleep and in the latter half are also excluded.
  • the SDNN plotted in FIG. 12B is distributed over 20 to 40 msec for the first day, whereas it is distributed over about 40 msec for the 105th day.
  • the deep time zone of sleep on day 1 and day 105 is relatively early at 25 to 80 minutes after falling asleep on the first day.
  • the 105th day is also supported by the lateness of 40-100 minutes after falling asleep. That is, the above experimental result means that it is difficult to extract the SDNN of only deep sleep depending on the elapsed time after falling asleep without determining the sleep stage based on the measurement result.
  • heart rate variability was extracted from the air pressure sensor, the R wave was detected, and the R wave could be specified with high accuracy in comparison with the electrocardiogram.
  • sleep and wakefulness could be appropriately distinguished by the coefficient of variation (CV) of the respiratory peak and the coefficient of variation of the R wave.
  • the power spectrum of the RR interval was calculated from the specified R wave, and it was shown that the sleep stage can be accurately determined by the spectrum pattern method and the autonomic component ratio method.
  • a suitable artificial intelligence system is used to learn the relationship between the power spectrum pattern of the RR interval and the determination result made by the conventional method at the same time, and the sleep stage determination result is obtained from the power spectrum. Obtainable.
  • a determination adapted to each subject can be performed using a file that determines a determination criterion for each subject.
  • the system shown in FIG. 7 can monitor and transmit changes in the interval between pulsations to the user and medical staff.
  • the user can take measures such as receiving an examination at a medical institution based on the information.

Abstract

Healthcare is executed by using a sleeping condition monitoring apparatus including a pressure waveform obtaining portion (10) for obtaining a pressure waveform from a pressure detecting portion being in contact with a user body, a pulsation extracting portion (22) for extracting a pulsating waveform by performing predetermined processing on the pressure waveform, a pulsating interval calculating portion (24) for calculating a pulsating interval, and a sleeping level determining portion (26) for determining a sleeping level, that is, a level of sleeping depth of the user from the pulsating interval, and collecting data of the pulsating interval during the sleep through a computer network, in order to measure sleeping conditions at home.

Description

睡眠状態モニタリング装置、モニタリングシステムおよびコンピュータプログラムSleep state monitoring device, monitoring system, and computer program
 本発明は、睡眠状態モニタリング装置、モニタリングシステムおよびコンピュータプログラムに関する。 The present invention relates to a sleep state monitoring device, a monitoring system, and a computer program.
 現在、人の健康管理を行ったり心身へのストレスを適切に管理するための手法のひとつとして、睡眠の質を管理することが求められている。睡眠の質を具体的に評価するには、睡眠を6つのカテゴリ、すなわち、覚醒、REM(Rapid Eye Movements)睡眠、ステージ1(S1)~ステージ4(S4)に分類して評価することが行われており、このためには、通常、電気的計測が利用されている。なお、覚醒、REM睡眠、S1~S4、は、この順に睡眠の深度が大きくなるが、通常の睡眠の状態の推移では、睡眠開始直後にREM睡眠になることはなく、一旦深い眠り(S3またはS4)になってから覚醒に向かって眠りが浅くなる際にREM睡眠の状態になる。この測定では、人(使用者または被験者)の身体のいくつかの部位に電極(皿電極)を装着して、例えば、頭皮からの脳波、アゴからの筋電位、眼の周辺からの眼球電位、および心電図を利用する必要がある。脳波から睡眠状態を判定する技術が、Rechtschaffen-Kales(R-K)の方法として知られている。また、心電図を利用する技術は、例えば特許文献1(特開2003-225211号公報)に開示されている。これらの測定手法による睡眠状態の評価は精度が高いという特徴を有するものの、正しい評価結果を得るためには、専門的知識を有する専門家または医療従事者がデータの解析を行う必要がある。また、電気的測定のために電極を測定時に装着するため、使用者に対する拘束性が高く、睡眠中の使用者に対して精神的な負担を与えてしまう。このため、これらの従来の技術は、睡眠の状態に影響すると共に、長期間にわたる継続的な測定を行うことも難しい。 Currently, it is required to manage the quality of sleep as one of the methods for managing human health and appropriately managing mental and physical stress. To specifically evaluate the quality of sleep, sleep is classified into six categories, namely, awakening, REM (Rapid Eye Movements) sleep, and stage 1 (S1) to stage 4 (S4). For this purpose, electrical measurement is usually used. In addition, awakening, REM sleep, and S1 to S4 increase the depth of sleep in this order. However, in the transition of the normal sleep state, REM sleep does not occur immediately after the start of sleep, and deep sleep (S3 or When the sleep becomes shallower toward awakening after S4), the REM sleep state is entered. In this measurement, electrodes (dish electrodes) are attached to several parts of a person's (user or subject) body, for example, an electroencephalogram from the scalp, a myoelectric potential from the jaw, an ocular potential from the periphery of the eye, And you need to use ECG. A technique for determining a sleep state from an electroencephalogram is known as the Rechtschaffen-Kales (RK) method. A technique using an electrocardiogram is disclosed in, for example, Patent Document 1 (Japanese Patent Laid-Open No. 2003-225211). Although the evaluation of the sleep state by these measurement methods has a feature that the accuracy is high, in order to obtain a correct evaluation result, it is necessary for a specialist or a medical worker having specialized knowledge to analyze the data. In addition, since the electrodes are attached at the time of measurement for electrical measurement, the user's restraint is high, and a mental burden is imposed on the sleeping user. For this reason, these conventional techniques affect the sleep state and are difficult to perform continuous measurement over a long period of time.
 また、睡眠状態の評価を上記以外の生体情報によって実現する手法について、多くの報告がなされている。この手法は、体振動の情報の採り方によりいくつかのバリエーションがあるが、その一つとして、体振動の大きさを指標として定義する方法(非特許文献1)が知られている。 In addition, many reports have been made on techniques for realizing sleep state evaluation using biological information other than the above. This method has some variations depending on how the body vibration information is taken. One of them is a method of defining the body vibration magnitude as an index (Non-Patent Document 1).
 さらに、自律神経成分比の研究によって睡眠ステージをごく大きく分類することが報告されている。非特許文献3には、心電図波形に基づいて自律神経成分比を求めて、浅い睡眠(S1またはS2)と、深い睡眠(S3またはS4)との間で判定を行うことが開示されている。 Furthermore, it has been reported that sleep stages are classified into very large groups by studying the autonomic component ratio. Non-Patent Document 3 discloses that an autonomic component ratio is obtained based on an electrocardiogram waveform and a determination is made between shallow sleep (S1 or S2) and deep sleep (S3 or S4).
 そして、 非特許文献5には、心疾患における心拍のスペクトルの変化について開示されており、また、非特許文献6には、心拍変動のパワースペクトルを応用して、糖尿病による初期段階の自律神経機能不全を検出しうることについて開示されている。
特開2003-225211号公報 渡辺崇士、渡辺嘉二郎: 就寝時無拘束計測生体データによる睡眠段階の推定、計測自動制御学会論文集、38(7), 581-589, (2002) 早野順一郎他: 心拍変動と自律神経機能、生物物理、28, 198―202, (1988) Baharav A, et al. : Fluctuations in autonomic nervous activity during sleep displayed by power spectrum analysis of heart rate variability, Neurology, 45, 1183―1187, (1995) Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology: Heart Rate Variability(HRV):Standards of Measurement, Physiological Interpretation, and Clinical Use, Circulation, 93(5), 1043-1065, (1996) Saul JP, Arai Y, Berger RD. et al.: Assessment of Autonomic regulation in chronic congestive heart failure by heart rate spectral analysis, American Journal of Cardiology, 61, 1292-1299, (1988) Yamasaki, Y, Ueda, N., Kishimoto, M et al.: Assessment of early stage autonomic nerve dysfunction in diabetic subjects; application of power spectral analysis of heart rate variability, Diabetes Research, 17(2), 73-80, (1991)
Non-Patent Document 5 discloses changes in the spectrum of heartbeats in heart disease, and Non-Patent Document 6 applies the power spectrum of heartbeat variability to provide an autonomic nervous function at an early stage due to diabetes. It is disclosed that a failure can be detected.
JP 2003-225211 A Takashi Watanabe, Yoshijiro Watanabe: Sleep stage estimation based on unconstrained measurement biological data, Proceedings of the Society of Instrument and Control Engineers, 38 (7), 581-589, (2002) Junichiro Hayano et al .: Heart rate variability and autonomic function, Biophysics, 28, 198-202, (1988) Baharav A, et al. : Fluctuations in autonomic novel activity draining sleep displayed by power spectrum of heart rate variability, 118, 118 Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology: Heart Rate Variability (HRV): Standards of Measurement, Physiological Interpretation, and Clinical Use, Circulation, 93 (5), 1043-1065, (1996) Saul JP, Arai Y, Berger RD. et al .: Assessment of Autonomous regulation in chronic congruent heart failure by heart analysis, American Journal of Cardiology, 1988, American Journal of Cardio. Yamazaki, Y, Ueda, N .; Kishimoto, M et al. : Application of early stage of humanization, diabetic sub-objects, application of power-specific idiosyncratic, 17
 専門家または医療従事者の助けを必要とせずに使用者が在宅のまま睡眠の質を簡単にチェックできるような装置や方法があれば、自己の健康状態を把握することが可能となり、健康管理が容易になる。また、睡眠の質の低下が手軽に検出できれば、専門家または医療従事者が診断するための的確な情報を与えることができる。また、非特許文献3の手法では、心電図を測定する必要があるばかりでなく、睡眠ステージを細分するような判定は難しく、また、睡眠か覚醒かを判定することもできない。さらに、在宅のままの使用者の拍動間隔のゆらぎのデータを収集して使用者の健康状態を監視することも容易に行うことはできない。本発明は、上記いくつかの課題を解決することを課題とする。 A device or method that allows the user to easily check the quality of sleep while staying at home without the assistance of a specialist or health care professional makes it possible to grasp his / her health condition and manage his / her health. Becomes easier. Moreover, if a decrease in the quality of sleep can be detected easily, accurate information for diagnosis by a specialist or a healthcare professional can be given. In the method of Non-Patent Document 3, it is not only necessary to measure an electrocardiogram, but it is difficult to determine whether the sleep stage is subdivided, and it is not possible to determine whether it is sleep or wakefulness. Furthermore, it is not easy to monitor the health condition of a user by collecting data on fluctuations in the pulsation interval of the user who remains at home. An object of the present invention is to solve the above several problems.
 本発明においては、使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、該拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定部とを備える睡眠状態モニタリング装置が提供される。また本発明の睡眠状態モニタリングは、方法として、あるいは、圧力検出部に接続されたコンピュータにおいてコンピュータプログラムとして実施することもできる。 In the present invention, a pressure waveform acquisition unit that obtains a pressure waveform as an electric signal by a pressure detection unit that is in direct contact with any part of the user's body or through clothing, and a predetermined process is performed on the pressure waveform. A pulsation extraction unit that extracts a pulsation waveform of the user and calculates a pulsation interval that is a time interval from the previous pulsation for each pulsation in the pulsation waveform; There is provided a sleep state monitoring device including a sleep stage determination unit that determines a sleep stage that is a stage of the user's sleep depth from the beat interval. The sleep state monitoring of the present invention can also be implemented as a method or as a computer program in a computer connected to the pressure detection unit.
 本明細書において、睡眠状態とは、覚醒、REM睡眠、ステージ1(S1)~ステージ4(S4)という覚醒と睡眠の全ての状態を指すために用い、睡眠ステージとは、睡眠状態から覚醒を除いたもの、すなわち、REM睡眠、S1~S4を指すために用いる。また、拍動とは、心臓の拍動およびそれに影響されて生じる振動、例えば、動脈の脈動による身体表面の振動を含む。 In this specification, the sleep state is used to indicate all states of awakening, REM sleep, awakening and sleep of stage 1 (S1) to stage 4 (S4), and the sleep stage refers to awakening from a sleeping state. Used to refer to the exclusion, ie, REM sleep, S1-S4. The pulsation includes the pulsation of the heart and vibrations generated by the pulsation, for example, vibrations of the body surface due to arterial pulsations.
 本発明においては、使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、コンピュータネットワークを通じてデータを送信するデータ送信部とを備えたネットワークに接続されたモニタリング部と、該コンピュータネットワークを通じて該装置から少なくとも拍動間隔を含む時系列データを複数回受信し、受信した時系列データを使用者に対応させて該時系列データを記録するための記憶部とを有する睡眠状態モニタリングシステムであって、前記コンピュータネットワークに接続されたいずれかのコンピュータが、複数回にわたって受信した前記時系列データを使用者に対応付けて時間に対するグラフとして表示し、過去のある期間における該使用者の拍動間隔のゆらぎの変化を提示するデータ表示部を備えている、睡眠状態モニタリングシステムが提供される。この睡眠状態モニタリングシステムは、モニタリング部にコンピュータネットワークを通じて接続可能にされているコンピュータにおいてプログラムとして実施することもできる。 In the present invention, a pressure waveform acquisition unit that obtains a pressure waveform as an electric signal by a pressure detection unit that is in direct contact with any part of the user's body or through clothing, and a predetermined process is performed on the pressure waveform. A pulsation extraction unit that extracts a pulsation waveform of a user, a pulsation interval calculation unit that calculates a pulsation interval that is a time interval from the previous pulsation for each pulsation in the pulsation waveform, and a computer A monitoring unit connected to a network including a data transmission unit for transmitting data through the network; and receiving time series data including at least a pulsation interval from the device through the computer network a plurality of times. A sleep state monitoring system having a storage unit for recording the time series data corresponding to a user, the computer network comprising: Any computer connected to the work displays the time-series data received multiple times in association with the user as a graph with respect to time, and changes in fluctuation of the user's pulsation interval in a past period A sleep state monitoring system including a data display unit that presents This sleep state monitoring system can also be implemented as a program in a computer that can be connected to the monitoring unit through a computer network.
 本発明のいずれかの実施形態によれば、使用者を拘束せず睡眠状態にも影響を与えることなく、覚醒しているか睡眠しているか、あるいは、使用者が睡眠しているのがどの睡眠ステージであるかなどの睡眠状態を監視することが可能になる。また、本発明のいずれかの実施形態によれば、拍動間隔のゆらぎによって使用者の健康状態を監視することが可能になる。 According to any embodiment of the present invention, which sleep is awake or sleeping, or which user is sleeping, without restraining the user and affecting the sleep state It becomes possible to monitor a sleep state such as whether it is a stage. In addition, according to any embodiment of the present invention, it is possible to monitor the health condition of the user based on fluctuations in the pulsation interval.
本発明の睡眠状態モニタリング装置の実施形態を示すブロック図。The block diagram which shows embodiment of the sleep state monitoring apparatus of this invention. 本発明の睡眠状態モニタリング方法の実施形態を示すフローチャート。The flowchart which shows embodiment of the sleep state monitoring method of this invention. 本発明の睡眠状態モニタリング方法の実施形態において、睡眠か覚醒かを判定する処理を示すフローチャート。The flowchart which shows the process which determines whether it is sleep or awakening in embodiment of the sleep state monitoring method of this invention. 本発明の睡眠状態モニタリング方法の実施形態において、睡眠ステージを判定する処理を示すフローチャート。The flowchart which shows the process which determines a sleep stage in embodiment of the sleep state monitoring method of this invention. 本発明の睡眠状態モニタリング方法の実施形態において、スペクトルパターン法によって睡眠ステージを判定する処理を示すフローチャート。The flowchart which shows the process which determines a sleep stage by the spectrum pattern method in embodiment of the sleep state monitoring method of this invention. 本発明の睡眠状態モニタリング方法の実施形態において、自律神経成分比法によって睡眠ステージを判定する処理を示すフローチャート。The flowchart which shows the process which determines a sleep stage by the autonomic nerve component ratio method in embodiment of the sleep state monitoring method of this invention. 本発明の睡眠状態モニタリングシステムの実施形態を示すブロック図。The block diagram which shows embodiment of the sleep state monitoring system of this invention. 本発明の実施例においてR波を特定する処理を行った結果を示すコンピュータ画面の例。The example of the computer screen which shows the result of having performed the process which specifies R wave in the Example of this invention. 本発明の実施例において被験者から得た体振動からRR間隔のパワースペクトルを得た結果を示すコンピュータ画面の例であり、覚醒状態(図9(A))、S1~S4(図9(B)~図9(E))およびREM睡眠(図9(F))の場合のパワースペクトルを示す。FIG. 9 is an example of a computer screen showing a result of obtaining a power spectrum at an RR interval from body vibration obtained from a subject in an example of the present invention, and is an arousal state (FIG. 9A) and S1 to S4 (FIG. 9B). FIG. 9 (E)) and the power spectrum in the case of REM sleep (FIG. 9 (F)). 本発明の実施例において被験者から得た体振動によって睡眠状態を判定した結果を示すコンピュータ画面の例。The example of the computer screen which shows the result of having determined the sleep state by the body vibration obtained from the test subject in the Example of this invention. 本発明の実施例において被験者から得た体振動によって睡眠状態を自動判定した結果と、従来の手法によって判定した結果を対比して示すコンピュータ画面の例。The example of the computer screen which compares and shows the result of having determined the sleep state automatically by the body vibration obtained from the subject in the embodiment of the present invention and the result determined by the conventional method. 本発明の実施例において被験者から測定したRR間隔のゆらぎ(SDNN)を示すグラフ。The graph which shows the fluctuation (SDNN) of RR interval measured from the test subject in the Example of this invention.
符号の説明Explanation of symbols
 1 睡眠状態モニタリング装置
 2 モニタリング部
 10 圧力波形取得部
 12 空気圧センサー
 14 トランスデューサ部
 16 A/D変換機
 20、50 コンピュータ
 24 拍動間隔算出部
 26 睡眠ステージ判定部
 32 呼吸波形抽出部
 34 呼吸間隔算出部
 36 覚醒判定部
 42 睡眠状態判定部
 46 データ送信部
 6 睡眠状態モニタリングシステム
 60 コンピュータネットワーク
 70 記憶部
 80 データ生成部
DESCRIPTION OF SYMBOLS 1 Sleep state monitoring apparatus 2 Monitoring part 10 Pressure waveform acquisition part 12 Air pressure sensor 14 Transducer part 16 A / D converter 20, 50 Computer 24 Beating interval calculation part 26 Sleep stage determination part 32 Respiration waveform extraction part 34 Respiration interval calculation part 36 Awakening determination unit 42 Sleep state determination unit 46 Data transmission unit 6 Sleep state monitoring system 60 Computer network 70 Storage unit 80 Data generation unit
1.睡眠状態の判定
 以下図面を参照して、本発明の実施の形態について説明する。本発明の実施の形態においては、体振動波形から睡眠状態を判定する。体振動波形の取得は、図1に示したように、圧力波形取得部10とコンピュータ20とを備える睡眠状態モニタリング装置1によって行う。圧力波形取得部10は、空気圧センサー12、圧力を電気信号に変換するトランスデューサ部14、および、その電気信号をアナログ/デジタル変換するA/Dコンバータ16を含んでおり、本発明においては、これらを集合的に圧力波形取得部10と呼ぶ。空気圧センサー12は、例えばベッド100に配置されていて、その上に使用者が横たわることができるように作られている。例えば、気密性の高い袋状のエンベロープにスポンジ状の詰め物が封入されていて、全体として1cm~3cm程度の厚み、30cm四方の正方形形状となるクッション状になっていて、使用者がその上に身体を横たえても睡眠が可能なような柔軟性や肌触りとなるように作られている。そして、そのエンベロープの周縁部の一箇所には、内部の空気と連通して圧力を伝播させることができるチューブが接続されていて、そのチューブがトランスデューサ部14に延びている。トランスデューサ部14には、そのチューブが接続される気密室があり、その気密室の一つの側面には開口部が設けられて、その開口部を塞ぐようにセラミックス製の圧電素子が配置されている。その圧電素子には、空気圧センサー12の内部の圧力波が伝達してアナログ値の圧電信号が生成され、その圧電信号は圧力波信号となって適当な増幅回路によって増幅される。増幅された圧力波信号は、A/D変換機16に入力され、デジタル化された圧力波信号が出力される。
1. Determination of Sleep State An embodiment of the present invention will be described below with reference to the drawings. In the embodiment of the present invention, the sleep state is determined from the body vibration waveform. As shown in FIG. 1, the body vibration waveform is acquired by the sleep state monitoring apparatus 1 including the pressure waveform acquisition unit 10 and the computer 20. The pressure waveform acquisition unit 10 includes an air pressure sensor 12, a transducer unit 14 that converts pressure into an electrical signal, and an A / D converter 16 that performs analog / digital conversion of the electrical signal. These are collectively referred to as a pressure waveform acquisition unit 10. The air pressure sensor 12 is disposed on the bed 100, for example, and is made so that a user can lie on the air pressure sensor 12. For example, a sponge-like stuffing is enclosed in a highly air-tight bag-like envelope, and the overall thickness is about 1 cm to 3 cm, and the cushion shape is a square shape of 30 cm square. It is designed to be flexible and soft enough that you can sleep even when you lie down. A tube that can communicate with internal air and propagate pressure is connected to one peripheral portion of the envelope, and the tube extends to the transducer unit 14. The transducer section 14 has an airtight chamber to which the tube is connected. An opening is provided on one side of the airtight chamber, and a ceramic piezoelectric element is disposed so as to close the opening. . A pressure wave inside the air pressure sensor 12 is transmitted to the piezoelectric element to generate an analog value piezoelectric signal, which becomes a pressure wave signal and is amplified by an appropriate amplifier circuit. The amplified pressure wave signal is input to the A / D converter 16, and a digitized pressure wave signal is output.
 デジタル化された圧力波信号は、コンピュータ20に入力される。コンピュータ20は、プログラムされた命令に従って、所定の処理を行う機能ブロックの集合として動作する。すなわち、コンピュータ20には、拍動抽出部22と、拍動間隔算出部24と、睡眠ステージ判定部26と、呼吸波形抽出部32と、呼吸間隔算出部34と、覚醒判定部36と、睡眠状態判定部42とが備えられる。その結果は、例えばモニター44に出力して視覚的に示すなどの手法によって使用者の睡眠状態として提示される。このように、圧力波形取得部と、コンピュータ20に実現される各機能ブロックによって、本発明の実施形態の睡眠状態モニタリング装置が構成される。なお、当業者には明らかであるように、コンピュータ20に備えられる適当な記憶手段によってデータが蓄積されたり、入出力手段によって処理の開始や休止、終了、中断などが制御されたりすることができる。 The digitized pressure wave signal is input to the computer 20. The computer 20 operates as a set of functional blocks that perform predetermined processing in accordance with programmed instructions. That is, the computer 20 includes a pulsation extraction unit 22, a pulsation interval calculation unit 24, a sleep stage determination unit 26, a respiration waveform extraction unit 32, a respiration interval calculation unit 34, an awakening determination unit 36, and sleep. A state determination unit 42. The result is presented as the user's sleep state by a technique such as outputting to the monitor 44 and visually indicating the result. As described above, the sleep state monitoring apparatus according to the embodiment of the present invention is configured by the pressure waveform acquisition unit and each functional block realized by the computer 20. As will be apparent to those skilled in the art, data can be accumulated by appropriate storage means provided in the computer 20, and the start, stop, end, interruption, etc. of the process can be controlled by the input / output means. .
 次に、本発明の睡眠状態モニタリング装置の実施形態における処理のフローの概要について説明する。本発明の睡眠状態モニタリング装置は、睡眠か覚醒かの判定と睡眠ステージの判定とからなる睡眠状態の判定を行うことができる。具体的には、図2に示したように、睡眠状態の判定のためには、まず、体振動データを取得する(S10)。そして、睡眠か覚醒かの判定する(S20)。睡眠である場合には(S30、YES)、睡眠ステージの判定をさらに行う(S40)。その後、あるいは、睡眠ではない場合(覚醒である場合、S30、NO)には、判定結果を出力する。このように、本発明の実施形態における体振動を用いた睡眠状態の判定の処理は、覚醒しているか睡眠しているかを判定する処理と、睡眠ステージを判定する処理とに分かれる。図示していないが、モニタリングの目的としては、この一連の処理を一定期間ごと(例えば1分ごと)に繰り返して、データの取得から睡眠状態の判定までを定期的に行ったり、あるいは、データを取得する処理のみを定期的に行い、適当な記憶手段にデータを格納しておいて、判定の処理を事後的に行ったりすることができる。 Next, an outline of a processing flow in the embodiment of the sleep state monitoring apparatus of the present invention will be described. The sleep state monitoring apparatus of the present invention can determine a sleep state consisting of determination of sleep or awakening and determination of a sleep stage. Specifically, as shown in FIG. 2, in order to determine the sleep state, body vibration data is first acquired (S10). And it is judged whether it is sleep or awakening (S20). If it is sleep (S30, YES), the sleep stage is further determined (S40). After that, or when it is not sleep (when it is awakening, S30, NO), a determination result is output. As described above, the process for determining the sleep state using the body vibration in the embodiment of the present invention is divided into a process for determining whether the person is awake or sleeping and a process for determining the sleep stage. Although not shown, for the purpose of monitoring, this series of processing is repeated at regular intervals (for example, every minute) to periodically perform the process from data acquisition to sleep state determination, or Only the processing to be acquired can be performed periodically, the data can be stored in an appropriate storage means, and the determination processing can be performed afterwards.
 以下、本発明の実施形態において実施される睡眠ステージの判定手法と睡眠状態の判定手法について説明する。図3には、睡眠か覚醒かの判定するステップ(図2、S20)の詳細を示し、また、図4~6には、睡眠ステージの判定のステップ(S40)の詳細を示す。 Hereinafter, a sleep stage determination method and a sleep state determination method implemented in the embodiment of the present invention will be described. FIG. 3 shows details of the step of determining whether it is sleep or awakening (FIG. 2, S20), and FIGS. 4 to 6 show details of the step of determining the sleep stage (S40).
1.1.睡眠と覚醒との判定(S20)
 使用者の睡眠状態のうち睡眠か覚醒かを判定するために、本発明は、体振動波形から呼吸波形を抽出して利用する形態により実施することができる。呼吸波形は、人の安静時の呼吸数からおのずとその範囲が定まり、例えば体振動のうちの0.13~0.70Hzの周波数領域成分を呼吸波形とすることができる。このため、体振動波形をフィルタ処理し、0.13~0.70Hzの周波数領域成分を抽出する(図3、S202)。呼吸波形からは、呼吸間隔のゆらぎを抽出する。このために、呼吸波形が極大値を示す時刻を算出してその時刻を呼吸のタイミングを代表する時刻(呼吸ピーク時刻)とし(S204)、各呼吸の呼吸ピーク時刻データからから直前の呼吸の呼吸ピーク時刻データを減算することにより、呼吸間隔のデータを生成する(S206)。そして、呼吸間隔についての変動係数(CV)を求めて(S208)、これを呼吸間隔の一様性についての指標とする。というのは、使用者が睡眠中である場合には呼吸間隔は一様性が高いため変動係数の値が小さくなるのに対し、使用者が覚醒中である場合には呼吸間隔は一様性が低くなって変動係数の値が大きくなるためである。具体的には、呼吸間隔のCVが0.2程度を閾値として設定し(S210)、それより大きい場合には覚醒状態(S212)、小さい場合には睡眠状態と判定することができる(S214)。
1.1. Determination of sleep and awakening (S20)
In order to determine whether the user's sleep state is sleep or wakefulness, the present invention can be implemented by a form in which a respiratory waveform is extracted from a body vibration waveform and used. The range of the respiration waveform is naturally determined from the respiration rate of the person at rest. For example, a frequency domain component of 0.13 to 0.70 Hz in body vibration can be used as the respiration waveform. For this reason, the body vibration waveform is filtered and a frequency domain component of 0.13 to 0.70 Hz is extracted (FIG. 3, S202). From the respiration waveform, the fluctuation of the respiration interval is extracted. For this purpose, the time when the respiration waveform shows the maximum value is calculated, and that time is set as the time representative of the respiration timing (respiration peak time) (S204). By subtracting the peak time data, respiration interval data is generated (S206). Then, a coefficient of variation (CV) for the breathing interval is obtained (S208), and this is used as an index for the uniformity of the breathing interval. This is because when the user is sleeping, the breathing interval is highly uniform, so the coefficient of variation is small, whereas when the user is awake, the breathing interval is uniform. This is because the value of the coefficient of variation increases as the value decreases. Specifically, the CV of the breathing interval is set as a threshold value of about 0.2 (S210), and when it is larger than that, it can be determined as an awake state (S212), and when smaller, it can be determined as a sleep state (S214). .
1.2.睡眠ステージの判定(S40)
 睡眠ステージを判定する処理は、心拍または脈動による拍動のR波によって拍動のタイミングを特定する段階(S402)と、ある心拍とその直前の心拍との間の時間間隔(RR間隔)のデータを得る段階(S404)と、RR間隔のパワースペクトルを算出する段階(S406)と、睡眠ステージを判定する段階(S42)とに分けて実行される。
1.2. Sleep stage determination (S40)
The process of determining the sleep stage includes the step of identifying the timing of the pulsation by the R wave of the pulsation due to the heartbeat or pulsation (S402), and the data of the time interval (RR interval) between a certain heartbeat and the immediately preceding heartbeat The process is divided into the step (S404) of obtaining the power spectrum, the step of calculating the power spectrum of the RR interval (S406), and the step of determining the sleep stage (S42).
1.2.1.R波のタイミングの特定(S402)
 睡眠ステージを判定する処理は、まず、体振動から、心電図(ECG)で言うところのR波に対応する波の部分を特定する。なお、R波の記載は、理解の容易のためにECGにおけるR波との表現を用いたものであり、特にECGのR波に直接対応するものだけに限定されるものではない。連なっている心拍や脈動によって身体の表面に伝わる圧力波において拍動同士の時間間隔を測定するために各拍動に対してその拍動を代表するタイミングとして定めるための特徴波形であれば任意である。このためには、体振動波形から拍動に当たる波を抽出する必要がある(図1、拍動抽出部22)。この抽出のためには、まず、体振動波形をバンドパスフィルターによって0.1Hz~30Hzの範囲に帯域制限する。次に、フィルター後の体振動波形を時間に関して二次微分する処理を行って、体振動の加速度の各時間のデータを算出する。そして、その体振動の加速度の時間毎のデータを、さらにバンドパスフィルターによって3~30Hzの範囲に帯域制限して、その波形データがピーク値を示す時刻を特定する。拍動抽出部22はそのための微分ピーク処理部22Aを備えることができる。さらに、1秒未満の周期の体振動波形を用いて、波形データが急激に変化する区間に限定し、4秒程度の体振動波形の周期性とR波の周期性に注目して、取り出したピーク値の時刻がほぼ等間隔となるように補正して決定する。この時刻は、3~30Hzの範囲で体振動の加速度が大きいと評価される時刻となるため、心臓の拍動のうち心室の興奮による振動つまり心電図におけるR波の時刻に相当するものと本願発明者は推測している。しかしながら、体振動から振動を抽出するという振動検出の性質上、その時刻には、心拍のみによる振動だけではなく、心臓による動脈による振動(脈動)などの影響があり得る。
1.2.1. Specification of R wave timing (S402)
In the process of determining the sleep stage, first, the part of the wave corresponding to the R wave referred to in the electrocardiogram (ECG) is specified from the body vibration. The description of the R wave uses an expression of the R wave in the ECG for easy understanding, and is not limited to the one that directly corresponds to the EC G R wave. In order to measure the time interval between beats in a pressure wave transmitted to the surface of the body by a continuous heartbeat or pulsation, any characteristic waveform can be used as long as it is determined as a timing representative of each beat. is there. For this purpose, it is necessary to extract a wave corresponding to a pulsation from the body vibration waveform (FIG. 1, pulsation extracting unit 22). For this extraction, first, the body vibration waveform is band-limited to a range of 0.1 Hz to 30 Hz by a band pass filter. Next, the body vibration waveform after the filter is subjected to a second-order differentiation with respect to time, and data of acceleration of body vibration at each time is calculated. Then, the time-dependent data of the acceleration of the body vibration is further band-limited to a range of 3 to 30 Hz by a bandpass filter, and the time when the waveform data shows the peak value is specified. The pulsation extraction unit 22 may include a differential peak processing unit 22A for that purpose. Furthermore, using a body vibration waveform with a period of less than 1 second, it was limited to the section in which the waveform data changed rapidly, and was extracted focusing on the periodicity of the body vibration waveform of about 4 seconds and the periodicity of the R wave. The peak time is corrected and determined so as to be substantially equally spaced. Since this time is a time when it is evaluated that the acceleration of the body vibration is large in the range of 3 to 30 Hz, the present invention corresponds to the time of the vibration due to the ventricular excitement of the heart beat, that is, the time of the R wave in the electrocardiogram Guesses. However, due to the nature of vibration detection in which vibration is extracted from body vibration, the time may be affected not only by vibration due to the heartbeat but also by vibration (pulsation) caused by the artery due to the heart.
1.2.2.RR間隔データの導出(S404)
 このようにして得られたR波の時刻データからRR間隔(ある拍動とその直前の拍動との間の時間間隔)の時系列データを求める(図4、S404;図1、拍動間隔算出部24)。RR間隔は、各拍動のR波の時刻データから直前の拍動におけるR波の時刻データを減算することにより与える。このRR間隔のデータは、拍動のR波を検出した不均等間隔の各時刻における時間データである。このデータは、後の処理のため、一定のサンプリング間隔の時系列データにするのが好都合であり、補間処理を行ってRR間隔の時系列データを得る。これには、例えば2次補間などの処理を行うことができる。以上のようにして得られたRR間隔の時系列データから次に睡眠の深度を判定する。
1.2.2. Derivation of RR interval data (S404)
Time series data of the RR interval (time interval between a certain beat and the immediately preceding beat) is obtained from the R wave time data thus obtained (FIG. 4, S404; FIG. 1, beat interval). Calculation unit 24). The RR interval is given by subtracting the R wave time data of the immediately preceding beat from the R wave time data of each beat. The data of the RR interval is time data at each time of the unequal interval at which the pulsating R wave is detected. This data is conveniently converted into time-series data with a constant sampling interval for later processing, and interpolation processing is performed to obtain time-series data with an RR interval. For this, for example, processing such as secondary interpolation can be performed. Next, the sleep depth is determined from the time-series data of the RR interval obtained as described above.
1.2.3.RR間隔のパワースペクトルの算出(S406)
 RR間隔の時系列データからパワースペクトルを求めて利用する(S406;図1、パワースペクトル算出部26A)。なお、パワースペクトルの算出は既知のいずれの方法も用いることができ、例えば、高速フーリエ変換アルゴリズムを用いたり、Wiener-Khinchinの定理による自己相関に基づく手法などを用いることができる。パワースペクトルを算出するための窓関数(ウィンドウ関数)は任意であるが、判定の精度を高めるために2分間のデータを利用する。この点については後述する。
1.2.3. Calculation of power spectrum of RR interval (S406)
The power spectrum is obtained from the time series data of the RR interval and used (S406; FIG. 1, power spectrum calculation unit 26A). Note that any known method can be used to calculate the power spectrum. For example, a fast Fourier transform algorithm or a method based on autocorrelation based on the Wiener-Khinchin theorem can be used. The window function (window function) for calculating the power spectrum is arbitrary, but data for 2 minutes is used to increase the accuracy of the determination. This point will be described later.
1.2.4.睡眠ステージの判定
 次いで、睡眠ステージの判定行う(図4、S42;図1、睡眠ステージ判定部26)。睡眠ステージの判定は、睡眠している状態が、REM睡眠およびNon-REM睡眠(S1~S4)のいずれであるかを判定する処理である。睡眠ステージの判定について、本発明は2つの実施の形態によって実施することができる。一つはRR間隔の時系列データのパワースペクトルのパターンから睡眠ステージを判定する実施形態(スペクトルパターン法)であり、もう一つは、RR間隔の時系列データから自律神経機能の働きを示す指標を算出して睡眠ステージの判定を行うもの(自律神経成分比法)である。
1.2.4. Next, the sleep stage is determined (FIG. 4, S42; FIG. 1, sleep stage determination unit 26). The determination of the sleep stage is a process for determining whether the sleeping state is REM sleep or Non-REM sleep (S1 to S4). Regarding the determination of the sleep stage, the present invention can be implemented by two embodiments. One is an embodiment (spectrum pattern method) in which the sleep stage is determined from the power spectrum pattern of the time series data of the RR interval, and the other is an index indicating the function of the autonomic nervous function from the time series data of the RR interval. Is calculated to determine the sleep stage (autonomic nerve component ratio method).
1.2.4.1.睡眠ステージの判定:スペクトルパターン法
 従来の睡眠ステージの判定方法によって判定しておいた睡眠ステージとRR間隔のパワースペクトルとを対応付けておいて、各睡眠ステージについての典型的なスペクトルのパターン(「典型パターン」)を得ることができる。そして、ある被験者に対して測定されたRR間隔のパワースペクトルのパターン(以下、「未知パターン」)がいずれの典型パターンに類似するかを判定することによって、睡眠ステージを判定することができる。
1.2.4.1. Sleep Stage Determination: Spectral Pattern Method The sleep stage determined by the conventional sleep stage determination method is associated with the power spectrum of the RR interval, and a typical spectrum pattern (“ A typical pattern ") can be obtained. The sleep stage can be determined by determining which typical pattern is similar to the pattern of the power spectrum of the RR interval measured for a subject (hereinafter, “unknown pattern”).
 図5を参照しながらその詳細を説明する。未知パターンがいずれの典型パターンに類似するかの判定は、RR間隔のパワースペクトルのパターンのうち、LF範囲(Low Frequency;0.04~0.15Hz)およびHF範囲(High Frequency;0.15~0.40Hz)においてパワーの最大値を求めて(図5、S422)、そのパワーの値のn分の1の値(nは2以上の整数)を求めて、その値における周波数幅を評価対象として行う。具体的には、まず、nを2にして(S424)、RR間隔のパワースペクトルの最大値の1/nの値を求める(S426)。そして、パワースペクトルに極大値を与える周波数の両側において、最大値の1/nとなる周波数を求めてその周波数同士の差(周波数幅)を計算する(S428)。この処理を、nを逐次増加(インクリメント)して(S432)所定の上限Nまで繰り返す(S430)。各nに対してそのときの周波数幅が与えられるので、周波数幅のn依存性がえられる。典型パターンについて同様に求めておいた周波数幅のn依存性に対して、未知パターンの周波数幅を対比させて、典型パターンと未知パターンの周波数幅の差を残差として算出する(S434)。こうして、各睡眠ステージの典型パターンに対する未知パターンの一致度を数値として得ることができる。ここで、未知パターンのLF範囲のピークは典型パターンのLF範囲のパターンのピークになるべく一致するようにマッチングさせ、同様に未知パターンのHF範囲のピークは典型パターンのHF範囲のパターンのピークとマッチングさせるようにして、異なる周波数範囲の間でのマッチングは行わない。そして、未知パターンに対して最も高い一致度(少ない残差)を示す典型パターンとなるような睡眠ステージを、その未知パターンの睡眠ステージとする(S436)。なお、ここに記載したLF範囲とHF範囲の周波数の範囲は、後述する自律神経成分比法における記載も含めて典型的なものであり、厳密に上述の周波数範囲にのみ限定されず、例えば、LF範囲を0.05~0.15Hzなどとすることも本発明の開示の一部とする。 Details will be described with reference to FIG. The determination of which typical pattern is similar to the unknown pattern is made by determining the LF range (Low Frequency; 0.04 to 0.15 Hz) and the HF range (High Frequency; 0.15 to 0.15) among the power spectrum patterns of the RR interval. 0.40 Hz) to obtain the maximum value of power (FIG. 5, S422), to obtain a value of 1 / n of the power value (n is an integer of 2 or more), and evaluate the frequency width at that value Do as. Specifically, first, n is set to 2 (S424), and a value of 1 / n of the maximum value of the power spectrum in the RR interval is obtained (S426). Then, on both sides of the frequency that gives the maximum value to the power spectrum, a frequency that is 1 / n of the maximum value is obtained, and a difference (frequency width) between the frequencies is calculated (S428). This process is repeated by incrementing n (S432) and repeating until a predetermined upper limit N (S430). Since the frequency width at that time is given to each n, the n dependence of the frequency width can be obtained. The frequency width of the unknown pattern is compared with the n dependence of the frequency width obtained in the same manner for the typical pattern, and the difference between the frequency width of the typical pattern and the unknown pattern is calculated as a residual (S434). Thus, the coincidence degree of the unknown pattern with respect to the typical pattern of each sleep stage can be obtained as a numerical value. Here, the peaks of the LF range of the unknown pattern are matched so as to match the peaks of the pattern of the LF range of the typical pattern as much as possible. Similarly, the peaks of the HF range of the unknown pattern are matched with the peaks of the pattern of the HF range of the typical pattern. Thus, no matching is performed between different frequency ranges. Then, a sleep stage that becomes a typical pattern showing the highest degree of coincidence (small residual) with respect to the unknown pattern is set as the sleep stage of the unknown pattern (S436). The frequency ranges of the LF range and the HF range described here are typical, including the description in the autonomic nerve component ratio method described later, and are not limited to the above-described frequency range strictly. For example, Setting the LF range to 0.05 to 0.15 Hz is also part of the disclosure of the present invention.
 なお、スペクトルパターン法では、睡眠ステージの判定を行うことができるが、睡眠と覚醒のどちらにあるかの判定、すなわち、REM睡眠と覚醒との判定を行うことはできない。 Note that the spectral pattern method can determine the sleep stage, but cannot determine whether the sleep stage is awake or wakefulness, that is, the REM sleep or the awakening.
 本願発明者は、スペクトルパターン法において判定精度を高めるために、RR間隔のパワースペクトルについて、上述のLF範囲とHF範囲に加えて、VLF範囲(Very Low Frequency;0.01Hz~0.04Hz)のスペクトルを利用することが有用であることを見出している。VLF範囲のスペクトルとは、RR間隔のパワースペクトルにおいてLF範囲よりも低周波に現れるスペクトルである。VLF範囲のスペクトルを正しく求めるためには、RR間隔の時系列データはそれに見合った期間だけのデータが必要となり、必要なスペクトル範囲の下限周波数の逆数をウインドウ期間とすることが必要となる。例えば、0.01Hzのパワースペクトルを正しく求めるためには、このウインドウ期間を100秒程度の適当な時間にすることが有効であり、実際に測定する場合には、例えば2分間以上に設定する。なお、VLF範囲のスペクトルを利用しない場合にはこれよりも短い時間、例えば60秒とすることができる。 In order to increase the determination accuracy in the spectrum pattern method, the inventor of the present application has a VLF range (Very Low Frequency: 0.01 Hz to 0.04 Hz) in addition to the above LF range and HF range for the power spectrum of the RR interval. We have found it useful to use spectra. The spectrum in the VLF range is a spectrum that appears at a lower frequency than the LF range in the power spectrum of the RR interval. In order to correctly obtain the spectrum in the VLF range, the time series data of the RR interval needs to have data corresponding to the corresponding period, and the reciprocal of the lower limit frequency of the necessary spectrum range needs to be the window period. For example, in order to correctly obtain a power spectrum of 0.01 Hz, it is effective to set this window period to an appropriate time of about 100 seconds. When actually measuring, for example, it is set to 2 minutes or more. If the spectrum in the VLF range is not used, a shorter time, for example, 60 seconds can be used.
 スペクトルパターン法の実施形態として典型パターンを用いた一致判定を行う手法を上述したが、本発明は任意のパターン判定手法によって実施することができる。例えば、スペクトルのパターンに対してそれが判定されるべき結果(睡眠ステージのいずれか)を対応させる学習処理をニューラルネットワークに対して行っておいて、未知パターンに対してそのニューラルネットワークを作用させて判定する、いわゆる教師付ニューラルネットワーク学習システムを用いることができる。その他の任意の教師付学習の手法を実装することもできる。 Although the method of performing coincidence determination using a typical pattern has been described above as an embodiment of the spectral pattern method, the present invention can be implemented by any pattern determination method. For example, a learning process is performed on a neural network corresponding to a result (one of sleep stages) to be determined for a spectrum pattern, and the neural network is applied to an unknown pattern. A so-called supervised neural network learning system for determination can be used. Any other supervised learning technique can be implemented.
1.2.4.2.睡眠ステージの判定:自律神経成分比法
 自律神経機能の働きを示す自律神経成分比によっても、睡眠ステージの判定は可能である。このためには、まず、RR間隔のパワースペクトルにおいて、LF成分とHF成分とを求める。そして、LF成分とHF成分の比、すなわちLF/HF比(自律神経成分比)を算出する。これは、被験者の自律神経機能の働きがLF/HF比により評価できるためである(非特許文献2~4参照)。この解析には、図6に示したように、LF範囲のパワースペクトルを算出してその範囲の積分値を得てLF成分とし(S442)、HF範囲のパワースペクトルを算出してその範囲の積分値を得てHF成分とする(S444)。そして、LF成分をHF成分で除算して得られるLF/HF比を求める(S446)。各睡眠ステージにおける典型的なLF/HF比とこのLF/HF比を比較して、睡眠ステージを判定する(S448)。LF/HF比は、小さいほど睡眠が深く大きいほど睡眠が浅いが、睡眠と覚醒のどちらにあるかの判定、すなわち、REM睡眠と覚醒との判定を行うことはできない。また、REM睡眠、S1~S4の各睡眠ステージとLF/HF比との関係には個人差が大きいので、例えばスペクトルパターン法による睡眠ステージの判定結果との対応関係を予め測定しておいて、使用者ごとに適応させることも有効である。
1.2.4.2. Sleep Stage Determination: Autonomic Nerve Component Ratio Method The sleep stage can also be determined by the autonomic nerve component ratio indicating the function of the autonomic nerve function. For this purpose, first, an LF component and an HF component are obtained in the power spectrum of the RR interval. Then, the ratio between the LF component and the HF component, that is, the LF / HF ratio (autonomic nerve component ratio) is calculated. This is because the function of the subject's autonomic nerve function can be evaluated by the LF / HF ratio (see Non-Patent Documents 2 to 4). In this analysis, as shown in FIG. 6, the power spectrum of the LF range is calculated to obtain the integral value of the range to obtain the LF component (S442), and the power spectrum of the HF range is calculated to integrate the range. A value is obtained and set as an HF component (S444). Then, an LF / HF ratio obtained by dividing the LF component by the HF component is obtained (S446). The sleep stage is determined by comparing the typical LF / HF ratio in each sleep stage with this LF / HF ratio (S448). The smaller the LF / HF ratio is, the deeper the sleep is, the deeper the sleep is, and the lighter is the sleep. However, it is not possible to determine whether the sleep is wake or wake. In addition, since there is a large individual difference in the relationship between the REM sleep, each sleep stage of S1 to S4 and the LF / HF ratio, for example, the correspondence relationship with the determination result of the sleep stage by the spectral pattern method is measured in advance. Adapting to each user is also effective.
 なお、本願発明者らは、HF成分がある程度の値に満たない時刻に取得したデータでは、HF成分つまりRR間隔に高い周波数のゆらぎが少なく、睡眠ステージの判定に誤りが生じやすいことを見出している。睡眠ステージの判定に誤りが生じにくいのは、HF成分の値が例えば300msec以上の値になっている場合である。このため、睡眠ステージの判定の前に、HF成分が所定の値より大きいかどうかを判定することによって、睡眠ステージの判定の精度の優劣を推測することができる。特に、このようなHF成分の値に下限を設定して事前に用いることは、深い睡眠ステージであるS3あるいはS4にあるかどうかを正確に判定したい場合に特に有効である。 The inventors of the present application have found that the data acquired at a time when the HF component does not reach a certain value has few fluctuations of high frequency in the HF component, that is, the RR interval, and an error is likely to occur in the sleep stage determination. Yes. An error hardly occurs in the determination of the sleep stage when the value of the HF component is, for example, 300 msec 2 or more. For this reason, it is possible to estimate the superiority or inferiority of the accuracy of the sleep stage determination by determining whether the HF component is larger than a predetermined value before determining the sleep stage. In particular, setting a lower limit for the value of such an HF component and using it in advance is particularly effective when it is desired to accurately determine whether or not it is in the deep sleep stage S3 or S4.
 以上の実施の形態により、使用者の睡眠状態に影響を与えにくく、使用者が在宅で利用することができ、また、睡眠か覚醒かを判定できるような睡眠状態モニタリングシステムまたは方法が実現される。 According to the above embodiment, a sleep state monitoring system or method is realized that hardly affects the sleep state of the user, can be used at home by the user, and can determine whether it is sleep or awakening. .
2.睡眠時の拍動間隔のモニタリング
 次に、睡眠時の拍動のモニタリングの手法について説明する。図7にこのモニタリング処理を行うための構成を示す。本発明では、コンピュータ50と圧力波形取得部10を集合的にモニタリング部2と呼ぶ。図1と同様に、体振動波形の取得には、圧力波形取得部10、すなわち、空気圧センサー12、圧力を電気信号に変換するトランスデューサ部14、および、その電気信号をアナログ/デジタル変換するA/Dコンバータ16を利用する。空気圧センサー12がベッド100に配置されていることができるのも同様である。また、図1と同様にトランスデューサ部14およびA/D変換機16が接続されて、デジタル化された圧力波信号が出力される。
2. Monitoring of beat interval during sleep Next, a technique for monitoring beat during sleep will be described. FIG. 7 shows a configuration for performing this monitoring process. In the present invention, the computer 50 and the pressure waveform acquisition unit 10 are collectively referred to as the monitoring unit 2. As in FIG. 1, the body vibration waveform is acquired by acquiring the pressure waveform acquisition unit 10, that is, the air pressure sensor 12, the transducer unit 14 that converts pressure into an electrical signal, and the A / A that converts the electrical signal into analog / digital. A D converter 16 is used. Similarly, the air pressure sensor 12 can be disposed on the bed 100. Similarly to FIG. 1, the transducer unit 14 and the A / D converter 16 are connected to output a digitized pressure wave signal.
 デジタル化された圧力波信号は、コンピュータ50に入力される。コンピュータ50は、プログラムされた命令に従って、所定の処理を行う機能ブロックの集合として動作する。すなわち、コンピュータ50には、拍動抽出部22と、拍動間隔算出部24と、データ送信部46とが備えられる。 The digitized pressure wave signal is input to the computer 50. The computer 50 operates as a set of functional blocks that perform predetermined processing in accordance with programmed instructions. That is, the computer 50 includes a pulsation extraction unit 22, a pulsation interval calculation unit 24, and a data transmission unit 46.
 データ送信部46は、例えばインターネットやイントラネットなどのコンピュータネットワーク60を通じて記憶部70にデータをアップロードすることができる。記憶部70には、使用者のIDと対応付けたデータが複数回アップロードされており、コンピュータネットワーク60に接続されたいずれかのコンピュータにおいて、複数回にわたって受信したデータをつなげることによって複数回のアップロードされたデータを一つのグラフに表示することができる。このとき、使用者ごとにグラフに表示するために、同じ使用者のデータを検索して連結させるデータ生成部80が用いられる。そして、例えばモニター86に表示することができる。 The data transmission unit 46 can upload data to the storage unit 70 through a computer network 60 such as the Internet or an intranet. Data associated with the user's ID is uploaded to the storage unit 70 a plurality of times, and the upload is performed a plurality of times by connecting the data received a plurality of times in any computer connected to the computer network 60. Displayed data can be displayed in one graph. At this time, in order to display on a graph for each user, a data generation unit 80 for searching and connecting data of the same user is used. For example, it can be displayed on the monitor 86.
 アップロードされるデータは、拍動間隔のデータを含む任意の時系列データである。例えば、拍動間隔の時系列データそれのみであってもよいし、図1の睡眠ステージ判定部26と同様の睡眠ステージ判定部(図示しない)が備えられている場合には、拍動間隔のデータに加えて睡眠ステージの時系列データとすることもできる。図1の呼吸波形抽出部32と、呼吸間隔算出部34と、覚醒判定部36と、睡眠状態判定部42とに相当する機能ブロック(図示しない)を備えていれば、拍動間隔のデータに加えて、覚醒か睡眠の判定結果が含まれている時系列データとすることもできる。 * Uploaded data is arbitrary time series data including beat interval data. For example, only the time series data of the pulsation interval may be used alone, or when a sleep stage determination unit (not shown) similar to the sleep stage determination unit 26 of FIG. In addition to data, it can be time series data of the sleep stage. If functional blocks (not shown) corresponding to the respiratory waveform extraction unit 32, the respiratory interval calculation unit 34, the awakening determination unit 36, and the sleep state determination unit 42 of FIG. In addition, time-series data including a determination result of arousal or sleep can be used.
 データ生成部80は、モニター86に表示するためのデータを生成するほか、拍動間隔のデータに見られるゆらぎの強さに、通常とは異なる明らかな変化が生じるのを検知してアラーム信号を出力するゆらぎ検知部82を備えることができる。なお、ゆらぎ検知部82は、データ生成部80の一部としているが、コンピュータネットワーク60に接続されて使用者ごとのデータを記憶部70から呼び出すことができれば、コンピュータネットワーク60に接続された任意のコンピュータ機器に備えられていることができる。 In addition to generating data to be displayed on the monitor 86, the data generation unit 80 detects an obvious change that is different from the normal in the intensity of fluctuations found in the data of the beat interval, and generates an alarm signal. A fluctuation detecting unit 82 for outputting can be provided. The fluctuation detection unit 82 is a part of the data generation unit 80. However, if the fluctuation detection unit 82 is connected to the computer network 60 and can call up data for each user from the storage unit 70, the fluctuation detection unit 82 can be arbitrarily connected to the computer network 60. It can be provided in computer equipment.
 ゆらぎ検知部82は、目的に応じて長さが定まるある期間(例えば1年間にわたるデータを記憶部70に蓄積する場合には、一ヶ月などの任意の所定期間)における拍動間隔の変化を監視している。この変化を検出するために、拍動間隔は、例えば拍動間隔データのもつ標準偏差の値である「SDNN」(Standard Deviation of Normal-to-Normal Intervals)のデータとすることもでき、変動係数の値などの数値となるように整理され、その値がある閾値を下回るかどうかという点で監視される。そして、例えば、その閾値を下回ったときにアラーム信号が送信されて、適当な手段によって使用者にフィードバックされる。このためには、例えば自動送信される電子メールや、パーソナライズされた管理用Webページにアラートを表示するなどの手法を用いることができる。 The fluctuation detection unit 82 monitors a change in the pulsation interval in a certain period of time determined according to the purpose (for example, an arbitrary predetermined period such as one month when data for one year is accumulated in the storage unit 70). is doing. In order to detect this change, the pulsation interval can be, for example, “SDNN” (Standard Deviation of Normal-to-Normal Intervals) data, which is the standard deviation value of the pulsation interval data, and the coefficient of variation. It is arranged so that it becomes a numerical value such as the value of and is monitored in terms of whether the value falls below a certain threshold. Then, for example, an alarm signal is transmitted when the value falls below the threshold value, and is fed back to the user by appropriate means. For this purpose, for example, a method of displaying an alert on an automatically transmitted e-mail or a personalized management Web page can be used.
 特にSDNNの値によって健康状態を監視する場合には、本願発明者らは、睡眠ステージの判定結果を組み合わせることが有効であることを見出した。すなわち、被験者の睡眠ステージがS1やS2である時刻に取得した拍動間隔データを用いると、そこから算出したSDNNが被験者の健康状態を良好に反映したものとならない場合があるのに対し、睡眠ステージがS3またはS4となる時刻において取得した拍動間隔データから算出したSDNNは、被験者の健康状態を良好に反映したものとなる。この原因は必ずしも明らかではないが、本願発明者は、睡眠ステージが浅いS1、S2のときには、被験者の睡眠状態が安定しておらず、その状態で取得されるSDNNには被験者の健康状態以外にも睡眠が浅いことを原因とする拍動間隔の変動が影響するのに対し、睡眠ステージが深いときにはそのような健康状態以外の原因が減少するためではないかと推測している。また、本発明の場合には、拍動の監視を圧力波形取得部によって行うことそれ自体が深い睡眠ステージS3あるいはS4の判定に適したものでもあるため、その点からも、拍動間隔のゆらぎであるSDNNを深い睡眠ステージにおいて取得することは測定システムの特質に適うものである。 In particular, when the health state is monitored by the value of SDNN, the present inventors have found that it is effective to combine the determination result of the sleep stage. That is, when using the pulsation interval data acquired at the time when the sleep stage of the subject is S1 or S2, the SDNN calculated therefrom may not be a good reflection of the health condition of the subject. The SDNN calculated from the beat interval data acquired at the time when the stage is S3 or S4 is a good reflection of the health condition of the subject. Although this cause is not necessarily clear, the inventor of the present application has a sleep state of shallow S1 and S2, and the sleep state of the subject is not stable. However, it is speculated that this may be caused by a decrease in causes other than the health condition when the sleep stage is deep, while fluctuations in the pulsation interval due to shallow sleep are affected. In the case of the present invention, since the monitoring of pulsation by the pressure waveform acquisition unit itself is suitable for the determination of the deep sleep stage S3 or S4, the fluctuation of the pulsation interval is also from that point. Obtaining SDNN in the deep sleep stage is suitable for the characteristics of the measurement system.
 これにより、在宅で測定される拍動間隔に現れる体調変化をいち早くキャッチして、睡眠時の拍動間隔として現れる使用者の体調の変化を知ることができ、体調管理を行うことができる。 This makes it possible to quickly catch a change in physical condition that appears in the beat interval measured at home, know the change in the user's physical condition that appears as the beat interval during sleep, and manage the physical condition.
3.睡眠ステージの判定における従来の手法との相関の確認
 睡眠ステージの判定については、従来から確立されている手法と本発明の実施形態とを比較するために、被験者に、従来の手法による睡眠ステージの判定と本発明の実施の形態による睡眠の判定を行った実験結果について説明する。20歳代の男女25名を被験者とし、被験者に睡眠(昼眠)を取らせ、そのうちの90分間を測定に利用した。
3. Confirmation of correlation with conventional method in determination of sleep stage Regarding determination of sleep stage, in order to compare the method established in the past with the embodiment of the present invention, the subject is informed of the sleep stage according to the conventional method. An experimental result of determination and sleep determination according to the embodiment of the present invention will be described. Twenty-five men and women in their 20s were subjects, and the subjects were allowed to sleep (nap), of which 90 minutes were used for measurement.
3.1. 睡眠ステージの判定における従来の手法
 被験者から脳波、アゴ筋電図、眼球電位図、心電図を採取して従来の手法による睡眠ステージの判定を行った。各測定値を得るために、被験者の身体の各部には皿電極をつけた。具体的には、脳波のためには頭頂に、アゴ筋電図のためにアゴに、眼球電位図のために眼の左右周辺に、そして、心電図のために足首接地により両手首に、それぞれ電極を装着した。皿電極には低周波増幅器を接続して電圧信号を増幅し、AD変換器によってディジタル値に変換したのち、データログ用のコンピュータによって記録をとった。脳波によって睡眠ステージを判定し、その際、前述のRechtschaffen-Kales(R-K)の方法に従った。
3.1. Conventional method for determination of sleep stage The brain stage, electromyogram, electrooculogram, and electrocardiogram were collected from subjects and the sleep stage was determined by a conventional method. In order to obtain each measurement value, a plate electrode was attached to each part of the subject's body. Specifically, electrodes are applied to the top of the head for electroencephalograms, to the jaws for electromyograms, to the left and right sides of the eye for electrocardiograms, and to both wrists by ankle contact for electrocardiograms. Attached. A low-frequency amplifier was connected to the plate electrode to amplify the voltage signal, converted to a digital value by an AD converter, and recorded by a data logging computer. The sleep stage was determined by an electroencephalogram, and the method of Rechtschaffen-Kales (RK) described above was followed.
3.2.本発明の実施の形態の睡眠ステージの判定の手法
 心拍動波形からR波を抽出し自律神経成分比を算出すると共に、呼吸波形も取得した。
3.2. Method for Determining Sleep Stage According to Embodiment of the Present Invention An R wave was extracted from a heartbeat waveform to calculate an autonomic component ratio, and a respiratory waveform was also acquired.
 本実施の形態の手法のための測定および解析のための装置としては、図1に記載の装置である圧力波形取得部10およびコンピュータ20を用い、そのコンピュータでは図1のような各機能ブロックとしてコンピュータを動作させ、図4~6のような動作を行わせるために、体振動から心拍動成分と呼吸成分を抽出するソフトウエアモジュールと、心拍動波形からR波を抽出してスペクトルパターン法の解析を行うソフトウエアモジュールと、自律神経成分比法を実施するためのソフトウエアモジュールとを利用した。 As a device for measurement and analysis for the method of the present embodiment, the pressure waveform acquisition unit 10 and the computer 20 which are the devices described in FIG. 1 are used, and in the computer, each functional block as shown in FIG. In order to operate the computer and perform the operations shown in FIGS. 4 to 6, a software module for extracting a heartbeat component and a respiratory component from body vibration, and an R wave from the heartbeat waveform to extract a spectrum wave A software module for analysis and a software module for performing the autonomic nerve component ratio method were used.
3.3.R波のタイミングの相関
 図8に、体振動からR波のタイミングを特定した例を示す。図8の上段のチャートC8Aには、体振動波形C8A-10と、それから算出した呼吸波形C8A-20も示している。中心線C8A-00の上部に示された波形C8A-2は、体振動波形C8A-10からR波の候補となる時刻を検出したものである。この検出は、拍動抽出部22によって行った。このとき、波形C8A-2のピークP1~P8のうち、ピークP2´は誤検出されたものである。このような誤検出が生じたのはノイズなどの影響であるため、誤検出のピークを除去する処理を行う。この処理は、R波のある候補が得られた時点と次の候補が得られた時点が、予め設定した時間(例えば、0.5秒)以下になった場合には、その次の候補は採用しない処理などであり、拍動によらずに得られた可能性の高い波形をR波として採用しないようにする処理である。この処理を行って候補から誤検出のピークを除去し、中心線C8A-00の下部に示したようにR波の時刻を特定している。なお、この図は8秒間の現象を拡大して表示しているものである。
3.3. Correlation of R-wave timing FIG. 8 shows an example in which R-wave timing is specified from body vibration. In the upper chart C8A of FIG. 8, the body vibration waveform C8A-10 and the respiration waveform C8A-20 calculated therefrom are also shown. A waveform C8A-2 shown in the upper part of the center line C8A-00 is obtained by detecting a candidate time of the R wave from the body vibration waveform C8A-10. This detection was performed by the pulsation extraction unit 22. At this time, among the peaks P1 to P8 of the waveform C8A-2, the peak P2 ′ is a false detection. Since such erroneous detection occurs due to the influence of noise or the like, processing for removing the erroneous detection peak is performed. This process is performed when the time point when a candidate with an R wave is obtained and the time point when the next candidate is obtained are less than or equal to a preset time (for example, 0.5 seconds), the next candidate is This is processing that is not employed, such as processing that does not employ a waveform that is highly likely to be obtained regardless of pulsation as an R wave. This processing is performed to remove erroneous detection peaks from the candidates, and the time of the R wave is specified as shown below the center line C8A-00. In this figure, the phenomenon of 8 seconds is enlarged and displayed.
 そして、従来の手法によって心電図から得たR波と、本発明の実施の形態によって得た体振動のデータから特定したR波との間で発生時点を比較した。その結果、空気圧センサーにより測定された体振動波形から求めたR波の発生時点の検出精度は、心電図から得たR波発生時点と比較して、両者の一致度は平均87%であった。また、このようにして算出したR波の時系列データ(0.1秒間隔の2次補間したもの)について、体振動から判定したRR間隔と心電図によるRR間隔との間の相関係数は、0.80~0.95であり、平均相関係数は、0.87であった。 Then, the generation time was compared between the R wave obtained from the electrocardiogram by the conventional method and the R wave specified from the body vibration data obtained by the embodiment of the present invention. As a result, the detection accuracy of the R wave generation time point obtained from the body vibration waveform measured by the air pressure sensor was 87% on average compared to the R wave generation time point obtained from the electrocardiogram. The correlation coefficient between the RR interval determined from the body vibration and the RR interval according to the electrocardiogram for the R-wave time-series data calculated in this way (secondarily interpolated at 0.1 second intervals) is: The average correlation coefficient was 0.87 between 0.80 and 0.95.
 なお、図8の中段のチャートC8Bに現れる波形C8B-10は呼吸波形であり、横軸は時間で11分間の波形であり、縦軸は圧力波形(任意単位)である。また、チャートC8Bの波形C8B-20は推定したRR間隔である。 The waveform C8B-10 appearing in the chart C8B in the middle of FIG. 8 is a respiratory waveform, the horizontal axis is a waveform of 11 minutes in time, and the vertical axis is a pressure waveform (arbitrary unit). The waveform C8B-20 in the chart C8B is the estimated RR interval.
 図8の下段のチャートC8Cは、RR間隔から上記実施形態の方法によって求めたパワースペクトルC8C-10である。ここに示したパワースペクトルC8C-10は、深い睡眠ステージにおける結果で、グラフに表示された周波数範囲は0.01~0.40Hzの範囲のものであり、パワースペクトルには高周波(HF)範囲の0.25~0.30Hzの位置に顕著なピークC8C-P1が現れた。 8 is a power spectrum C8C-10 obtained from the RR interval by the method of the above embodiment. The power spectrum C8C-10 shown here is a result in a deep sleep stage, and the frequency range displayed in the graph is in the range of 0.01 to 0.40 Hz, and the power spectrum has a high frequency (HF) range. A remarkable peak C8C-P1 appeared at a position of 0.25 to 0.30 Hz.
3.4.睡眠ステージ判定
 そして、体振動から本実施の形態の手法によって睡眠ステージを判定した。睡眠ステージの判定のためには、前述のように2分間のデータを用いた。
3.4. Sleep stage determination And the sleep stage was determined by the method of this Embodiment from the body vibration. For determination of the sleep stage, data for 2 minutes was used as described above.
3.4.1.スペクトルパターン法
 具体的には、図8に関連して説明したような手法によって得られたパワースペクトルを未知の睡眠状態にある被験者から得て、そこにスペクトルパターン法を適用して睡眠ステージを判定した。まず、スペクトルパターン法の判断対象になるスペクトルの観察結果について説明する。
3.4.1. Spectral pattern method Specifically, the power spectrum obtained by the method described in relation to FIG. 8 is obtained from a subject in an unknown sleep state, and the sleep pattern is determined by applying the spectral pattern method there. did. First, the observation result of the spectrum to be determined by the spectrum pattern method will be described.
 図9の(A)~(E)に、覚醒、S1~S4のそれぞれの睡眠状態の場合のパワースペクトルを示す。また、図9の(F)にREM睡眠のパワースペクトルを示す。これらのスペクトルにおいて、グラフ上は0.00~0.01Hzおよび0.40~0.50Hzの値も表示されているものの、利用する周波数範囲は0.01Hz~0.40Hzである。これらの図に示すように、RR間隔のパワースペクトルは、図9の(A)の覚醒時ではRR間隔の変動が不規則なものである。その結果、RR間隔のパワースペクトルはLF、HFの範囲に小さな成分しか有さず、0.04Hz以下のVLF範囲が主要なスペクトルとなった。 9A to 9E show power spectra in the case of awakening and sleep states S1 to S4. FIG. 9F shows a power spectrum of REM sleep. In these spectra, although values of 0.00 to 0.01 Hz and 0.40 to 0.50 Hz are also displayed on the graph, the frequency range to be used is 0.01 Hz to 0.40 Hz. As shown in these figures, the power spectrum of the RR interval has irregular fluctuations in the RR interval at the time of awakening in FIG. As a result, the power spectrum of the RR interval has only a small component in the LF and HF ranges, and the VLF range of 0.04 Hz or less is the main spectrum.
 また、図9の(B)および(C)に示す浅い睡眠ステージS1とS2においては、RR間隔の変動の波形(図示しない)にリズミカルな変化が見られ、パワースペクトルにもこれを反映するような明確なピークがLF範囲のスペクトルとHF範囲のスペクトルに現れた。このLF範囲のスペクトルは、睡眠ステージがS1からS2へ移行するとより大きくなった。 In addition, in the shallow sleep stages S1 and S2 shown in FIGS. 9B and 9C, a rhythmic change is observed in the waveform of RR interval fluctuation (not shown), and this is reflected in the power spectrum. Clear peaks appeared in the LF and HF range spectra. The spectrum of this LF range became larger when the sleep stage shifted from S1 to S2.
 睡眠ステージS1は覚醒段階からS2との段階の移行状態にあり、覚醒時の特徴に近い分布すなわちLFとHFに亘る小さいパワースペクトル値の分布が見られた。睡眠ステージS1に覚醒状態のパワースペクトルが混入してくる場合もあるため、覚醒状態とS1の識別を実際のパワースペクトルのみから行うことは難しい。これに対し、図9の(C)に示す睡眠ステージS2は、LF範囲のスペクトルが存在するもののHF範囲のスペクトルがS1よりも明確に出現してくるので、S1とS2の識別はS1と覚醒状態の識別よりも容易である。 The sleep stage S1 is in a state of transition from the awakening stage to the S2 stage, and a distribution close to the characteristics at the time of awakening, that is, a distribution of small power spectrum values over LF and HF was observed. Since the power spectrum of the awake state may be mixed into the sleep stage S1, it is difficult to distinguish the awake state from S1 only from the actual power spectrum. On the other hand, in the sleep stage S2 shown in FIG. 9C, the spectrum of the HF range appears more clearly than the spectrum of S1 although the spectrum of the LF range exists. Easier than state identification.
 図9の(D)およびEに示す深い睡眠ステージであるS3とS4においては、睡眠が深くなるにつれてHF範囲のパワースペクトルがLF範囲のパワースペクトルより大きくなり、S4においてHF範囲のピークが主要なパワースペクトルとなる。従って、S4とS3の識別は、S1とS2との識別と同様に容易である。 In S3 and S4, which are deep sleep stages shown in FIGS. 9D and 9E, the power spectrum of the HF range becomes larger than the power spectrum of the LF range as sleep becomes deeper, and the peak of the HF range is the main in S4. It becomes a power spectrum. Therefore, identification of S4 and S3 is as easy as identification of S1 and S2.
 そして、図9の(F)はREM睡眠のパワースペクトルパターンである。睡眠が覚醒に近づいたときに現れるREM睡眠では、覚醒及びS1とパターンが類似しているが、LF範囲のスペクトルが、S1の場合よりもより大きく複雑になるという特徴があった。 And (F) of FIG. 9 is a power spectrum pattern of REM sleep. In REM sleep that appears when sleep approaches wakefulness, the pattern is similar to wakefulness and S1, but the spectrum of the LF range is larger and more complex than in S1.
 本発明者の検討によれば、REM睡眠と覚醒の区別は、RR間隔の時系列データから変動係数(CV)を求めることによって区別することができた。また、覚醒の推定には、呼吸のCVも考慮した。この手法は、入眠時、中途覚醒、睡眠終了後の覚醒の判定に有効であった。以上、観察者がRR間隔のパワースペクトルを直接観察した場合の特徴について述べた。 According to the study of the present inventor, REM sleep and awakening can be distinguished by obtaining a coefficient of variation (CV) from time series data of RR intervals. Moreover, CV of respiration was also considered in the estimation of arousal. This method was effective in determining sleep awakening, mid-wake awakening, and awakening after the end of sleep. The characteristics when the observer directly observes the power spectrum at the RR interval have been described.
 そして、以上のパターンの特徴から睡眠ステージの判定を再現性が高いものとするために、睡眠ステージの識別のために上述のスペクトルパターン法を適用した。その結果を図10と図11に示す。 Then, in order to make the sleep stage determination highly reproducible from the characteristics of the above pattern, the above-described spectral pattern method was applied to identify the sleep stage. The results are shown in FIGS.
 図10は空気圧センサーから得た体振動波形C10A-10および呼吸波形C10A-20(最上段のチャートC10A)、検出したR波時系列C10B-10(上から2段目のチャートC10B)、抽出した呼吸数C10C-10と心拍数C10C-20(上から3段目のチャートC10C)、RR間隔のパワースペクトルC10D-10(最下段左のチャートC10D)、および、睡眠ステージの判定例C10E-10およびLF/HF比C10E-20(最下段右のチャートC10E)である。また、図11は、睡眠状態の判定結果を比較したものであり、折れ線R0は脳波を用いて従来の方法によって判定した結果であり、折れ線R1(破線)は、スペクトルパターン法によって自動判定した結果である。この自動判定結果の折れ線R1は、図10のチャートC10Eに示したものと同じ結果であり、縦軸は、1から-4にかけて、順に、覚醒、REM睡眠、S1~S4を表わしており、76分の体振動の測定データから得られたものである。従来の手法と実施例の手法との間での一致度は約75%となった。 FIG. 10 shows body vibration waveform C10A-10 and respiratory waveform C10A-20 (top chart C10A) obtained from the air pressure sensor, and detected R wave time series C10B-10 (second chart C10B from the top). Respiration rate C10C-10 and heart rate C10C-20 (third chart C10C from the top), RR interval power spectrum C10D-10 (bottom left chart C10D), and sleep stage determination example C10E-10 and The LF / HF ratio is C10E-20 (the bottom right chart C10E). Moreover, FIG. 11 compares the determination result of the sleep state, the broken line R0 is a result determined by a conventional method using an electroencephalogram, and the broken line R1 (broken line) is a result of automatic determination by a spectrum pattern method. It is. The line R1 of the automatic determination result is the same as that shown in the chart C10E of FIG. 10, and the vertical axis represents awakening, REM sleep, and S1 to S4 in order from 1 to −4. It was obtained from the measurement data of the minute body vibration. The degree of agreement between the conventional method and the method of the example was about 75%.
3.4.2.自律神経成分比法(LF/HF比を用いる方法)
 次に、自律神経成分比法の実施例について説明する。上述の条件により取得した体振動波形から、浅い睡眠ステージ(S1とS2)と深い睡眠ステージ(S3とS4)との区別を自律神経成分比法により行った。このとき、浅い睡眠と深い睡眠の判定、LF/HF比が0.3より大きいときに浅い睡眠とし、0.3以下のときに深い睡眠と判定した。取得したデータの範囲では、この判定結果と、脳波による睡眠ステージの判定結果との一致度は共に80%以上となった。なお、このLF/HF比については、個人差が大きいため上記の判定基準は一例であり、例えば、判定の用いる使用者ごとの値は、参照可能なデータベースに別途記録しておくことができる。また、好適な例として、上述のようにHF成分の値そのものが所定の値に達した場合、すなわち、睡眠ステージの判定精度が一定程度に達した場合に、睡眠ステージの判定を行うように制御することも可能である。また、別の好適な例として、HF成分が所定の値、例えば300msecに達しない場合、すなわち、睡眠ステージの判定精度が一定程度に達しない場合には、睡眠ステージの判定を行った結果を出力する際、精度が良好でないことを示すフラグを付すことが可能である。
3.4.2. Autonomic component ratio method (method using LF / HF ratio)
Next, an example of the autonomic nerve component ratio method will be described. The autonomic component ratio method was used to distinguish between the shallow sleep stages (S1 and S2) and the deep sleep stages (S3 and S4) from the body vibration waveform acquired under the above conditions. At this time, shallow sleep and deep sleep were judged as shallow sleep when the LF / HF ratio was greater than 0.3, and deep sleep was judged when 0.3 or less. In the range of the acquired data, the degree of coincidence between the determination result and the determination result of the sleep stage by the electroencephalogram was 80% or more. In addition, about this LF / HF ratio, since the individual difference is large, the above-mentioned determination standard is an example. For example, the value for each user used for the determination can be separately recorded in a referable database. As a preferred example, when the HF component value itself reaches a predetermined value as described above, that is, when the sleep stage determination accuracy reaches a certain level, control is performed so that the sleep stage is determined. It is also possible to do. As another preferred example, when the HF component does not reach a predetermined value, for example, 300 msec 2 , that is, when the sleep stage determination accuracy does not reach a certain level, the result of the sleep stage determination is obtained. When outputting, it is possible to attach a flag indicating that the accuracy is not good.
3.5.睡眠/覚醒の判定
 以上のスペクトルパターン法および自律神経成分比法では覚醒時とREM睡眠の区別は難しいが、呼吸ピークの変動係数(CV)やRR間隔の変動係数(CV)を用いることで、両睡眠状態が区別可能となった。本発明の実施の形態の睡眠と覚醒との判定について、呼吸の変動係数(CV)によって判定を行った。その結果、呼吸のCVの閾値を0.2として、それ以上の値を覚醒と判定することが有効であった。また、RR間隔の変動係数(CV)には個人差が生じたが、特定の個人では十分に高い再現性が実現できた。
3.5. Sleep / wake determination The spectral pattern method and autonomic component ratio method described above makes it difficult to distinguish between wakefulness and REM sleep, but by using the coefficient of variation of the respiratory peak (CV) and the coefficient of variation of the RR interval (CV), Both sleep states became distinguishable. The determination of sleep and wakefulness according to the embodiment of the present invention was performed based on the coefficient of variation (CV) of respiration. As a result, it was effective to set the respiration CV threshold to 0.2 and determine a value higher than that as arousal. In addition, although individual differences occurred in the coefficient of variation (CV) of the RR interval, sufficiently high reproducibility could be realized for specific individuals.
4.拍動モニタリング
 次に、図7に示した装置を用いて拍動間隔(RR間隔)をモニタリングする睡眠状態モニタリングシステム6の作製例を説明する。
4). Beat monitoring Next, a production example of the sleep state monitoring system 6 that monitors the beat interval (RR interval) using the apparatus shown in FIG. 7 will be described.
4.1.拍動モニタリング装置の動作
 モニタリング部2は、R波のデータに対するRR間隔の算出処理を行う。そしてそのRR間隔のデータを、記憶部70を管理するデータサーバに一日一度アップロードする。アップロードは任意のデータ通信手段を用いることができるが、例えば、FTP(ファイルトランスファープロトコル)によってプットすることにより行うことができる。記憶部70におけるデータの蓄積の態様は任意であり、例えば、データサーバのファイルシステムのフォルダ(ディレクトリ)を利用して使用者のIDごとおよび日付ごとのCSV(Comma Separated Values)形式のコンピュータファイルを複数蓄積することができる。別の例としては、データベースマネージメントシステムなどを利用して適当なデータベースに蓄積することができる。なお、睡眠開始(入眠)のタイミングを適切に判定するため、コンピュータ50には図示しない光センサーを接続しておいて、周囲の照度情報によって測定開始のタイミングを制御することができる。さらに、測定の時間を入眠後の睡眠サイクル2回分に限定するために、測定開始後3時間のデータのみ取得することができる。また、一定以上の深さの睡眠にある場合のみ拍動間隔のデータを収集したり、後に睡眠ステージの判定結果を用いてデータ処理を行うことを可能とするために、各時刻でのRR間隔データに、それぞれの時刻での睡眠ステージの判定結果を合わせて取得することができる。これにより、例えば睡眠ステージがS3とS4の場合の拍動データを利用して、被験者の健康状態を良好に反映するようなSDNNのデータを取得することができる。
4.1. Operation of the pulsation monitoring device The monitoring unit 2 performs an RR interval calculation process for R wave data. The data of the RR interval is uploaded once a day to a data server that manages the storage unit 70. Although any data communication means can be used for uploading, for example, it can be performed by putting using FTP (File Transfer Protocol). The form of data storage in the storage unit 70 is arbitrary. For example, a computer file in CSV (Comma Separated Values) format for each user ID and each date using a file system folder (directory) of the data server. Multiple data can be stored. As another example, data can be stored in an appropriate database using a database management system or the like. In order to appropriately determine the timing of sleep start (sleeping), an optical sensor (not shown) is connected to the computer 50, and the timing of measurement start can be controlled by ambient illuminance information. Furthermore, in order to limit the measurement time to two sleep cycles after falling asleep, only data for 3 hours after the start of measurement can be acquired. In addition, in order to be able to collect beat interval data only when sleeping at a certain depth or more, and to perform data processing using the determination result of the sleep stage later, the RR interval at each time The determination result of the sleep stage at each time can be acquired together with the data. Thereby, for example, using the pulsation data when the sleep stages are S3 and S4, it is possible to acquire SDNN data that favorably reflects the health condition of the subject.
 このCSV形式のファイルに含まれるRR間隔のデータから、RR間隔のゆらぎを算出する。このゆらぎの算出は、一定期間(例えば)に取得されたRR間隔のデータから標準偏差を算出してSDNNとして算出すること、あるいは、変動係数(CV)を算出することによって行うことができる。複数のCSVファイルに分かれている各日のデータは連結されて計算され、日付を追ってRR間隔のゆらぎを表示することができる。 * The fluctuation of the RR interval is calculated from the data of the RR interval included in the CSV file. This fluctuation can be calculated by calculating a standard deviation from the data of the RR interval acquired during a certain period (for example) and calculating it as SDNN, or by calculating a coefficient of variation (CV). The data of each day divided into a plurality of CSV files are concatenated and calculated, and the fluctuation of the RR interval can be displayed following the date.
 また、ゆらぎ検知部82には、例えば、過去一ヶ月間のRR間隔のゆらぎの変動を監視する短期変動監視モジュール(図示しない)と、過去1年間のRR間隔のゆらぎの変動を監視する長期変動監視モジュール(図示しない)を実装することができる。ここで、短期変動監視モジュールおよび長期変動監視モジュールは、RR間隔のゆらぎに現れる変動がどのように推移するかをそれぞれの期間に合わせて監視している。これらのモジュールでは、使用者の健康状態を監視するために、睡眠の時間の中から監視に適した睡眠状態にある時間において取得されたRR間隔データを抽出することができる。このような監視に適した状態にある時間としては、例えば、RR間隔に長周期のゆらぎが生じていない時間を選択したり、あるいは、睡眠ステージがS3あるいはS4にある時間を選択することができる。このような目的で、図7に示した睡眠状態モニタリングシステム6には、図1に示した睡眠ステージ判定部26、呼吸波形抽出部32、呼吸間隔算出部34、覚醒判定部36、睡眠状態判定部42などの機能ブロックをさらに追加しておいて、短期変動監視モジュールおよび長期変動監視モジュールがそれらの機能ブロックからの出力によって動作が制御されることができる。これらにより、図7に示した睡眠状態モニタリングシステム6においては、RR間隔のゆらぎがどのように変動するかを、在宅の使用者に特段の負担をかけることなく、継続して監視することができる。 The fluctuation detection unit 82 includes, for example, a short-term fluctuation monitoring module (not shown) that monitors fluctuations in the RR interval during the past month, and a long-term fluctuation that monitors fluctuations in the RR interval over the past year. A monitoring module (not shown) can be implemented. Here, the short-term fluctuation monitoring module and the long-term fluctuation monitoring module monitor how the fluctuation appearing in the fluctuation of the RR interval changes according to each period. In these modules, in order to monitor a user's health state, RR interval data acquired in the time in a sleep state suitable for monitoring can be extracted from the sleep time. As the time suitable for such monitoring, for example, the time when the long-period fluctuation is not generated in the RR interval can be selected, or the time when the sleep stage is in S3 or S4 can be selected. . For such a purpose, the sleep state monitoring system 6 shown in FIG. 7 includes the sleep stage determination unit 26, the respiratory waveform extraction unit 32, the respiratory interval calculation unit 34, the awakening determination unit 36, and the sleep state determination shown in FIG. By further adding functional blocks such as the unit 42, the operations of the short-term fluctuation monitoring module and the long-term fluctuation monitoring module can be controlled by outputs from these functional blocks. Thus, in the sleep state monitoring system 6 shown in FIG. 7, it is possible to continuously monitor how the fluctuation of the RR interval fluctuates without placing a special burden on the user at home. .
4.2.拍動モニタリングの測定例
 以上のようにして測定したSDNNの実際のデータについて説明する。血糖値の高い1名の被験者を対象として、本実施例にて作製した拍動モニタリング装置によってSDNNを測定した。測定は、その被験者が運動療法を開始する第1日にあたる日と、運動療法を開始して第105日にあたる日の2回行い、それぞれの日の日中に約100分の睡眠を被験者に取らせて、その睡眠中にRR間隔のデータを取得した。そして、RR間隔のデータを60秒分だけ蓄積し、その60秒分のRR間隔のデータから標準偏差を算出して、60秒ごとにSDNNを得た。このとき、SDNNのデータに合わせて睡眠ステージの判定を行うために、睡眠ステージのデータも取得した。ここでの睡眠ステージの判定は、上記の条件に合わせるために、60秒のウインドウ期間によって求めたRR間隔のパワースペクトルを用い、上述のスペクトルパターン法による判定を行った。ウインドウ期間が60秒であるため、VLF範囲のパワースペクトルは利用していない。なお、上記被験者は、第1日と第105日の間に、運動療法として、週に3日以上実施し、各週の歩行距離が40km程度となるウォーキングを続けていた。
4.2. Measurement example of pulsation monitoring Actual data of SDNN measured as described above will be described. SDNN was measured with a pulsation monitoring device prepared in this example for one subject with a high blood glucose level. The measurement is performed twice on the day corresponding to the first day when the subject starts exercise therapy and the day corresponding to the 105th day after starting exercise therapy, and the subject takes about 100 minutes of sleep during each day. RR interval data was acquired during the sleep. Then, RR interval data was accumulated for 60 seconds, a standard deviation was calculated from the RR interval data for 60 seconds, and SDNN was obtained every 60 seconds. At this time, in order to determine the sleep stage in accordance with the SDNN data, sleep stage data was also acquired. Here, the sleep stage was determined by the above-described spectrum pattern method using the power spectrum of the RR interval obtained by the window period of 60 seconds in order to meet the above conditions. Since the window period is 60 seconds, the power spectrum in the VLF range is not used. In addition, the said test subject continued for 3 days or more as a exercise therapy between the 1st day and the 105th day, and the walking where the walk distance of each week becomes about 40 km.
 図12(A)および(B)は、共に、第1日および第105日のSDNNのデータであり、横軸は入眠後の経過時間(単位:分)、縦軸はSDNN(単位:ミリ秒)とした。図12(A)は、1分ごとに測定したSDNNをそのままプロットしたグラフであり、図12(B)は、測定されたSDNNのうち、被験者が睡眠ステージS3またはS4にある時刻に得られたデータをプロットしたグラフである。図12(A)のみから傾向を観ると、第1日のSDNNは第105日のSDNNより小さな値が多いといえる。その一方で、例えば、入眠後30分程度までは大きな違いが見られない。また、入眠後90分以降はこの関係はむしろ逆転している。 FIGS. 12A and 12B are both SDNN data on the first day and the 105th day, the horizontal axis is the elapsed time after sleep (unit: minutes), and the vertical axis is SDNN (unit: milliseconds). ). FIG. 12A is a graph in which SDNN measured every minute is plotted as it is, and FIG. 12B is obtained at the time when the subject is in the sleep stage S3 or S4 among the measured SDNN. It is the graph which plotted data. Looking at the trend from FIG. 12A alone, it can be said that the SDNN on the first day has a smaller value than the SDNN on the 105th day. On the other hand, for example, there is no significant difference until about 30 minutes after falling asleep. In addition, this relationship is reversed after 90 minutes after falling asleep.
 これに対し、睡眠ステージがS3とS4の場合のみをプロットした図12(B)では、いずれの日のデータからも、それぞれのデータ系列において大きく偏位しているデータが除かれており、また、入眠直後や後半のデータも除かれている。そして、図12(B)においてプロットされたSDNNは、第1日のものは20~40msecに分布しているのに対し、第105日のものは概ね40msec以上に分布している。 On the other hand, in FIG. 12 (B) in which only the cases where the sleep stage is S3 and S4 are plotted, the data greatly deviated in each data series is excluded from the data of any day, Also, data immediately after falling asleep and in the latter half are also excluded. The SDNN plotted in FIG. 12B is distributed over 20 to 40 msec for the first day, whereas it is distributed over about 40 msec for the 105th day.
 さらに、図12(A)と(B)とを比較すると、図12(A)において、例えば第1日のデータのうちSDNNが大きな値を示しているものは、睡眠が浅いためにRR間隔の変動が大きくなり、SDNNが見かけ上大きな値を示していたことも分かる。
 これに加えて、入眠後の経過時間によって深い睡眠のSDNNが抽出できるか見てみても、図12(A)と(B)との比較から、睡眠中は睡眠の深度が時間と共に複雑に変化していて図12(A)にはデータがあるのに(B)にはデータが残らない時刻が散発的に生じている。このため、入眠後の経過時間によって深い睡眠の時間を定めることはできない。この点は、主に図12(B)から、第1日と第105日の睡眠の深い時間帯が第1日は概ね入眠後25~80分であり比較的早期であるのに対し、第105日は、入眠後40~100分と遅くなっていることによっても裏づけられる。すなわち、上記実験結果は、測定結果によって睡眠ステージの判定を行うことなく入眠後の経過時間に依拠して深い睡眠のみのSDNNを抽出することが困難であることを意味している。   
Further, when comparing FIG. 12A and FIG. 12B, in FIG. 12A, for example, among the data for the first day, the SDNN showing a large value indicates that the RR interval is low because sleep is shallow. It can also be seen that the fluctuations increased and that SDNN showed an apparently large value.
In addition to this, even if we look at whether deep sleep SDNN can be extracted according to the elapsed time after falling asleep, from the comparison between FIG. 12 (A) and FIG. In FIG. 12 (A), there are sporadic times when there is data but no data remains in (B). For this reason, the deep sleep time cannot be determined by the elapsed time after falling asleep. This point is mainly from FIG. 12 (B), in which the deep time zone of sleep on day 1 and day 105 is relatively early at 25 to 80 minutes after falling asleep on the first day. The 105th day is also supported by the lateness of 40-100 minutes after falling asleep. That is, the above experimental result means that it is difficult to extract the SDNN of only deep sleep depending on the elapsed time after falling asleep without determining the sleep stage based on the measurement result.
5.まとめ
 以上、本発明の実施例において、空気圧センサーから心拍変動を抽出して、R波を検出し、心電図との比較において、高い精度でR波を特定することができた。また、呼吸ピークの変動係数(CV)やR波の変動係数によって睡眠と覚醒とを適切に区別することができた。さらに、特定したR波からRR間隔を求めてRR間隔のパワースペクトルを算出し、スペクトルパターン法及び自律神経成分比法によって睡眠ステージの判定が精度よく実行できることを示した。
5). Summary As described above, in the embodiment of the present invention, heart rate variability was extracted from the air pressure sensor, the R wave was detected, and the R wave could be specified with high accuracy in comparison with the electrocardiogram. Moreover, sleep and wakefulness could be appropriately distinguished by the coefficient of variation (CV) of the respiratory peak and the coefficient of variation of the R wave. Furthermore, the power spectrum of the RR interval was calculated from the specified R wave, and it was shown that the sleep stage can be accurately determined by the spectrum pattern method and the autonomic component ratio method.
 睡眠/覚醒の判定においてR波の変動係数を用いたとき、あるいは、睡眠ステージの判定に自律神経成分比法を用いたときには、個人差の問題が残っているものの、特定の被験者に対しては高い再現性が得られた。 When using the R-wave coefficient of variation in sleep / wake determination or when using the autonomic component ratio method for sleep stage determination, although individual differences remain, for certain subjects High reproducibility was obtained.
 さらに、RR間隔のゆらぎがどのように変動するかを、在宅の使用者に特段の負担をかけることなく、継続して監視することができる。 Furthermore, it is possible to continuously monitor how the fluctuation of the RR interval fluctuates without placing a special burden on the user at home.
 以上、本発明の実施の形態につき述べたが、本発明は既述の実施の形態に限定されるものではなく、本発明の技術的思想に基づいて各種の変形、変更および組み合わせが可能である。たとえば、パターンの判定において、適当な人工知能システムを用いて、RR間隔のパワースペクトルパターンと、同時に従来の手法によって行った判定結果との関係を学習させて、パワースペクトルから睡眠ステージの判定結果を得ることができる。また、個人差が避けられない判定(例えばLF/HF比による判定)については、被験者別の判定基準を決めるようなファイルを利用して、被験者ごとに適応させた判定を実行することもできる。 Although the embodiments of the present invention have been described above, the present invention is not limited to the above-described embodiments, and various modifications, changes and combinations can be made based on the technical idea of the present invention. . For example, in determining a pattern, a suitable artificial intelligence system is used to learn the relationship between the power spectrum pattern of the RR interval and the determination result made by the conventional method at the same time, and the sleep stage determination result is obtained from the power spectrum. Obtainable. In addition, for a determination that an individual difference cannot be avoided (for example, a determination based on the LF / HF ratio), a determination adapted to each subject can be performed using a file that determines a determination criterion for each subject.
 また、図7に示したシステムによって、拍動の間隔における変化を監視して使用者や医療従事者に伝達することができる。使用者は、その情報に基づいて医療機関にて検査を受けるなどの対応を取ることができる。 In addition, the system shown in FIG. 7 can monitor and transmit changes in the interval between pulsations to the user and medical staff. The user can take measures such as receiving an examination at a medical institution based on the information.

Claims (28)

  1.  使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、
     該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、
     該拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、
     前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定部と
     を備える睡眠状態モニタリング装置。
    A pressure waveform acquisition unit that obtains a pressure waveform as an electric signal by a pressure detection unit in contact with any part of the user's body directly or through clothing;
    A pulsation extraction unit that performs a predetermined process on the pressure waveform to extract a pulsation waveform of the user;
    For each beat in the beat waveform, a beat interval calculating unit that calculates a beat interval that is a time interval from the immediately preceding beat;
    A sleep state monitoring device comprising: a sleep stage determination unit that determines a sleep stage that is a stage of the user's sleep depth from the pulsation interval.
  2.  前記圧力波形に所定の処理を行って呼吸波形を抽出する呼吸波形抽出部と、
     該呼吸波形における各回の呼吸について、直前の呼吸からの時間間隔である呼吸間隔を算出する呼吸間隔算出部と、
     使用者が睡眠しているか覚醒しているかを該呼吸間隔から判定する覚醒判定部と
     をさらに備える請求項1に記載の睡眠状態モニタリング装置。
    A respiration waveform extraction unit that performs a predetermined process on the pressure waveform to extract a respiration waveform;
    For each breath in the breathing waveform, a breath interval calculation unit that calculates a breath interval that is a time interval from the immediately preceding breath;
    The sleep state monitoring apparatus according to claim 1, further comprising: an awakening determination unit that determines whether the user is sleeping or awake from the breathing interval.
  3.  前記覚醒判定部が、前記拍動間隔から拍動間隔の変動係数を算出し、該拍動間隔の変動係数によって睡眠か覚醒かを判定するものである、請求項2に記載の睡眠状態モニタリング装置。 The sleep state monitoring device according to claim 2, wherein the awakening determination unit calculates a fluctuation coefficient of a beat interval from the beat interval, and determines whether it is sleep or awakening based on the fluctuation coefficient of the beat interval. .
  4.  前記覚醒判定部が、前記呼吸間隔から呼吸間隔の変動係数を算出し、該呼吸間隔の変動係数によって睡眠か覚醒かを判定する、請求項2に記載の睡眠状態モニタリング装置。 The sleep state monitoring device according to claim 2, wherein the awakening determination unit calculates a variation coefficient of a breathing interval from the breathing interval, and determines whether it is sleep or awakening based on the variation coefficient of the breathing interval.
  5.  前記拍動間隔算出部は、微分処理とピーク判定とを組み合わせて前記拍動波形から各拍動のタイミングを示す時刻を特定する微分ピーク処理部を有する、請求項1に記載の睡眠状態モニタリング装置。 The sleep state monitoring apparatus according to claim 1, wherein the pulsation interval calculation unit includes a differential peak processing unit that identifies a time indicating the timing of each pulsation from the pulsation waveform by combining differential processing and peak determination. .
  6.  前記睡眠ステージ判定部は、
     拍動間隔の時系列データから拍動間隔のゆらぎについてのパワースペクトルを算出するパワースペクトル算出部と、
     該パワースペクトルの示すパターンを判定対象として睡眠ステージを判定するスペクトルパターン解析部と
     を備えている、請求項1に記載の睡眠状態モニタリング装置。
    The sleep stage determination unit
    A power spectrum calculation unit for calculating a power spectrum about fluctuation of the beat interval from time series data of the beat interval;
    The sleep state monitoring apparatus according to claim 1, further comprising: a spectrum pattern analysis unit that determines a sleep stage using a pattern indicated by the power spectrum as a determination target.
  7.  前記スペクトルパターン解析部が、ニューラルネットワークの学習によるパターン判定部を備える、請求項6に記載の睡眠状態モニタリング装置。 The sleep state monitoring apparatus according to claim 6, wherein the spectrum pattern analysis unit includes a pattern determination unit based on learning of a neural network.
  8.  前記睡眠ステージ判定部は、拍動間隔の時系列データに含まれる拍動間隔のゆらぎについての低周波成分と高周波成分とを算出して該低周波成分と該高周波成分との比を計算するLF/HF比算出部を備え、該比の値によって睡眠ステージを判定するものである、請求項1に記載の睡眠状態モニタリング装置。 The sleep stage determination unit calculates a low frequency component and a high frequency component for fluctuation of the beat interval included in the time series data of the beat interval, and calculates a ratio between the low frequency component and the high frequency component. The sleep state monitoring apparatus according to claim 1, further comprising a / HF ratio calculation unit, wherein the sleep stage is determined based on a value of the ratio.
  9.  前記睡眠ステージ判定部が、前記高周波成分の値が所定の値を超える場合の前記比の値によって睡眠ステージを判定するものである、請求項8に記載の睡眠状態モニタリング装置。   The sleep state monitoring device according to claim 8, wherein the sleep stage determination unit determines a sleep stage based on the value of the ratio when the value of the high frequency component exceeds a predetermined value. *
  10.  前記所定の値が300msecである、請求項9に記載の睡眠状態モニタリング装置。 The sleep state monitoring apparatus according to claim 9, wherein the predetermined value is 300 msec 2 .
  11.  前記圧力波形取得部が、使用者の身体を支える寝具と身体との間に配置される空気圧センサーを含む、請求項1に記載の睡眠状態モニタリング装置。 The sleep state monitoring apparatus according to claim 1, wherein the pressure waveform acquisition unit includes an air pressure sensor disposed between a bedding supporting a user's body and the body.
  12.  使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得ステップと、
     該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出ステップと、
     該拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出ステップと、
     前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定ステップと
     を含む睡眠状態モニタリング方法
    A pressure waveform acquisition step of obtaining a pressure waveform as an electrical signal by a pressure detection unit in contact with any part of the user's body directly or through clothing;
    A pulsation extraction step of extracting a pulsation waveform of the user by performing a predetermined process on the pressure waveform;
    For each beat in the beat waveform, a beat interval calculating step for calculating a beat interval that is a time interval from the previous beat;
    A sleep stage monitoring method comprising: a sleep stage determination step of determining a sleep stage, which is a stage of the depth of sleep of the user, from the beat interval
  13.  前記圧力波形に所定の処理を行って呼吸波形を抽出する呼吸波形抽出ステップと、
     該呼吸波形における各回の呼吸について、直前の呼吸からの時間間隔である呼吸間隔を算出する呼吸間隔算出ステップと、
     使用者が睡眠しているか覚醒しているかを呼吸間隔から判定する覚醒判定ステップと
     をさらに含む請求項12に記載の睡眠状態モニタリング方法。
    A respiration waveform extracting step for extracting a respiration waveform by performing a predetermined process on the pressure waveform;
    For each breath in the breathing waveform, a breath interval calculating step for calculating a breath interval that is a time interval from the immediately preceding breath;
    The sleep state monitoring method according to claim 12, further comprising: an arousal determination step of determining whether the user is sleeping or awake from a breathing interval.
  14.  前記覚醒判定ステップが、前記拍動間隔から拍動間隔の変動係数を算出し、該拍動間隔の変動係数によって睡眠か覚醒かを判定するものである、請求項13に記載の睡眠状態モニタリング方法。 The sleep state monitoring method according to claim 13, wherein the awakening determination step calculates a fluctuation coefficient of a beat interval from the beat interval, and determines whether it is sleep or awakening based on the fluctuation coefficient of the beat interval. .
  15.  前記覚醒判定ステップが、前記呼吸間隔から呼吸間隔の変動係数を算出し、該呼吸間隔の変動係数によって睡眠か覚醒かを判定するものである、請求項13に記載の睡眠状態モニタリング方法。 The sleep state monitoring method according to claim 13, wherein the awakening determination step calculates a fluctuation coefficient of a breathing interval from the breathing interval and determines whether it is sleep or awakening based on the fluctuation coefficient of the breathing interval.
  16.  前記拍動間隔算出ステップは、微分処理とピーク判定とを組み合わせて前記拍動波形から各拍動のタイミングを示す時刻を特定する微分ピーク処理ステップを含む、請求項12に記載の睡眠状態モニタリング方法。 The sleep state monitoring method according to claim 12, wherein the pulsation interval calculating step includes a differential peak processing step of identifying a time indicating the timing of each pulsation from the pulsation waveform by combining differentiation processing and peak determination. .
  17.  前記睡眠ステージ判定ステップは、
     拍動間隔の時系列データから拍動間隔のゆらぎについてのパワースペクトルを算出するパワースペクトル算出ステップと、
     該パワースペクトルの示すパターンを判定対象として睡眠ステージを判定するスペクトルパターン解析ステップと
     を含む、請求項12に記載の睡眠状態モニタリング方法。
    The sleep stage determination step includes
    A power spectrum calculating step for calculating a power spectrum about fluctuation of the beat interval from time series data of the beat interval;
    The sleep pattern monitoring method according to claim 12, further comprising: a spectrum pattern analyzing step of determining a sleep stage using a pattern indicated by the power spectrum as a determination target.
  18.  前記スペクトルパターン解析ステップが、学習を行わせたニューラルネットワークにパターンを判定させるステップを含む、請求項17に記載の睡眠状態モニタリング方法。 The sleep state monitoring method according to claim 17, wherein the spectral pattern analysis step includes a step of causing a learned neural network to determine a pattern.
  19.  前記睡眠ステージ判定ステップは、拍動間隔の時系列データに含まれる拍動間隔のゆらぎについての低周波成分と高周波成分とを算出して該低周波成分と該高周波成分との比を計算するLF/HF比算出ステップを含み、該比の値によって睡眠ステージを判定するステップである、請求項12に記載の睡眠状態モニタリング方法。 In the sleep stage determination step, an LF that calculates a low frequency component and a high frequency component with respect to fluctuations of the beat interval included in the time series data of the beat interval, and calculates a ratio between the low frequency component and the high frequency component. The sleep state monitoring method according to claim 12, comprising a step of calculating a / HF ratio, wherein the sleep stage is determined based on a value of the ratio.
  20.  前記睡眠ステージ判定ステップが前記高周波成分の値が所定の値を超えるかどうかを判定するステップをさらに含む、請求項19に記載の睡眠状態モニタリング方法。 The sleep state monitoring method according to claim 19, wherein the sleep stage determination step further includes a step of determining whether or not a value of the high frequency component exceeds a predetermined value.
  21.  前記所定の値が300msecである、請求項20に記載の睡眠状態モニタリング方法。 The sleep state monitoring method according to claim 20, wherein the predetermined value is 300 msec 2 .
  22.  前記圧力波形取得ステップが、使用者の身体を支える寝具と身体との間に配置される空気圧センサーによって圧力波形を取得するステップを含む、請求項12に記載の睡眠状態モニタリング方法。 The sleep state monitoring method according to claim 12, wherein the pressure waveform acquisition step includes a step of acquiring a pressure waveform by an air pressure sensor disposed between a bedding supporting a user's body and the body.
  23.  使用者の身体のいずれかの部分に直接または着衣を介して接しており電気信号として圧力波形を得る圧力検出部に接続されたコンピュータに、請求項12~22に記載のいずれかの方法を実行させるためのコンピュータプログラム。 The method according to any one of claims 12 to 22 is executed on a computer connected to a pressure detection unit that obtains a pressure waveform as an electric signal by contacting any part of a user's body directly or through clothing. Computer program to let you.
  24.  使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、コンピュータネットワークを通じてデータを送信するデータ送信部とを備えたコンピュータネットワークに接続されたモニタリング部と、
     該コンピュータネットワークを通じて該モニタリング部から少なくとも拍動間隔を含む時系列データを複数回受信し、受信した時系列データを使用者に対応させて該時系列データを記録するための記憶部と
     を有する睡眠状態モニタリングシステムであって、前記コンピュータネットワークに接続されたいずれかのコンピュータが、複数回にわたって受信した前記時系列データを使用者に対応付けて時間に対するグラフとして表示し、過去のある期間における該使用者の拍動間隔のゆらぎの変化を提示するデータ表示部を備えている、睡眠状態モニタリングシステム。
    A pressure waveform acquisition unit that obtains a pressure waveform as an electric signal by a pressure detection unit that is in contact with any part of the user's body directly or through clothing, and performs a predetermined process on the pressure waveform to perform Data is transmitted via a computer network, a beat extraction unit that extracts a beat waveform, a beat interval calculation unit that calculates a beat interval that is a time interval from the previous beat for each beat in the beat waveform, and a computer network A monitoring unit connected to a computer network including a data transmission unit to perform,
    A sleep unit comprising: a storage unit configured to receive time series data including at least a pulsation interval from the monitoring unit a plurality of times through the computer network, and to record the time series data corresponding to the received time series data; A state monitoring system, wherein any one of the computers connected to the computer network displays the time-series data received over a plurality of times as a graph against time in association with a user, and the use in a past period A sleep state monitoring system comprising a data display unit that presents a change in fluctuation of a person's beat interval.
  25.  前記コンピュータネットワークに接続されたいずれかのコンピュータが、過去のある期間における該使用者の拍動間隔のデータにおいて、拍動間隔のゆらぎを示す数値についての前記期間より短い短期間における変動を検知してアラーム信号を出力するゆらぎ変化探知部を有している請求項24に記載の睡眠状態モニタリングシステム。 Any one of the computers connected to the computer network detects a fluctuation in a short period shorter than the period with respect to a numerical value indicating fluctuation of the pulsation interval in the data of the pulsation interval of the user in a past period. The sleep state monitoring system according to claim 24, further comprising a fluctuation change detection unit that outputs an alarm signal.
  26.  前記コンピュータネットワークに接続されたいずれかのコンピュータまたは前記モニタリング部が、前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定部をさらに有しており、
     前記記憶部が、前記拍動間隔のゆらぎを示す数値を睡眠ステージ判定部が判定した睡眠ステージと関連付けて記憶する、請求項24に記載の睡眠状態モニタリングシステム。
    Any one of the computers connected to the computer network or the monitoring unit further includes a sleep stage determination unit that determines a sleep stage that is a stage of the depth of sleep of the user from the beat interval,
    The sleep state monitoring system according to claim 24, wherein the storage unit stores a numerical value indicating fluctuation of the beat interval in association with the sleep stage determined by the sleep stage determination unit.
  27.  前記コンピュータネットワークに接続されたいずれかのコンピュータまたは前記モニタリング部が、前記使用者の睡眠の深度の段階である睡眠ステージを該拍動間隔から判定する睡眠ステージ判定部をさらに有しており、
     前記記憶部が、睡眠ステージ判定部が判定した睡眠ステージがステージ3またはステージ4の場合にのみ前記拍動間隔のゆらぎを示す数値を記憶する、請求項24に記載の睡眠状態モニタリングシステム。
    Any one of the computers connected to the computer network or the monitoring unit further includes a sleep stage determination unit that determines a sleep stage that is a stage of the depth of sleep of the user from the beat interval,
    The sleep state monitoring system according to claim 24, wherein the storage unit stores a numerical value indicating fluctuation of the pulsation interval only when the sleep stage determined by the sleep stage determination unit is stage 3 or stage 4.
  28.  使用者の身体のいずれかの部分に直接または着衣を介して接している圧力検出部により、電気信号として圧力波形を得る圧力波形取得部と、該圧力波形に所定の処理を行って使用者の拍動波形を抽出する拍動抽出部と、拍動波形における各拍動について、直前の拍動からの時間間隔である拍動間隔を算出する拍動間隔算出部と、コンピュータネットワークを通じてデータを送信するデータ送信部とを備え、コンピュータネットワークに接続されているモニタリング部に該コンピュータネットワークを通じて接続可能にされたコンピュータを、請求項24または25のいずれかに記載のコンピュータとして動作させるためのコンピュータプログラム。 A pressure waveform acquisition unit that obtains a pressure waveform as an electric signal by a pressure detection unit that is in contact with any part of the user's body directly or through clothing, and performs a predetermined process on the pressure waveform to perform Data is transmitted via a computer network, a beat extraction unit that extracts a beat waveform, a beat interval calculation unit that calculates a beat interval that is a time interval from the previous beat for each beat in the beat waveform, and a computer network 26. A computer program for operating a computer that is connected to a monitoring unit connected to a computer network through the computer network as a computer according to claim 24 or 25.
PCT/JP2008/072697 2008-06-13 2008-12-12 Sleeping condition monitoring apparatus, monitoring system, and computer program WO2009150765A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009141394A JP2011115188A (en) 2008-06-13 2009-06-12 Sleeping condition monitoring apparatus, monitoring system, and computer program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
PCT/JP2008/060857 WO2009150744A1 (en) 2008-06-13 2008-06-13 Sleeping state monitoring device, monitoring system, and computer program
JPPCT/JP2008/060857 2008-06-13

Publications (1)

Publication Number Publication Date
WO2009150765A1 true WO2009150765A1 (en) 2009-12-17

Family

ID=41416465

Family Applications (2)

Application Number Title Priority Date Filing Date
PCT/JP2008/060857 WO2009150744A1 (en) 2008-06-13 2008-06-13 Sleeping state monitoring device, monitoring system, and computer program
PCT/JP2008/072697 WO2009150765A1 (en) 2008-06-13 2008-12-12 Sleeping condition monitoring apparatus, monitoring system, and computer program

Family Applications Before (1)

Application Number Title Priority Date Filing Date
PCT/JP2008/060857 WO2009150744A1 (en) 2008-06-13 2008-06-13 Sleeping state monitoring device, monitoring system, and computer program

Country Status (2)

Country Link
JP (1) JPWO2009150744A1 (en)
WO (2) WO2009150744A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011148793A1 (en) * 2010-05-28 2011-12-01 Sharp Kabushiki Kaisha Method and system for reliable respiration parameter estimation from acoustic physiological signal
CN105496362A (en) * 2015-12-15 2016-04-20 成都麦杰康科技有限公司 Sleep state monitoring system and method
US20220125376A1 (en) * 2019-02-13 2022-04-28 Kyoto University Sleep apnea syndrome determination apparatus, sleep apnea syndrome determination method, and sleep apnea syndrome determination program

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2012043820A1 (en) * 2010-09-30 2014-02-24 国立大学法人 東京医科歯科大学 Method and apparatus for myocardial toxicity testing and cardiomyocyte evaluation
US20130231579A1 (en) * 2010-11-08 2013-09-05 Toyota Jidosha Kabushiki Kaisha Sleep state estimation device
CN103284727B (en) * 2012-02-23 2016-01-27 广东乐心医疗电子股份有限公司 Sleep monitor method and apparatus
JP5874489B2 (en) * 2012-03-27 2016-03-02 富士通株式会社 Sleep state determination device and sleep state determination method
JP6622455B2 (en) * 2014-11-19 2019-12-18 シャープ株式会社 Sleep state determination system
JP6466729B2 (en) * 2015-02-10 2019-02-06 アイシン精機株式会社 Activity determination system
JP6439729B2 (en) 2016-03-24 2018-12-19 トヨタ自動車株式会社 Sleep state estimation device
JP6875981B2 (en) * 2016-11-22 2021-05-26 パラマウントベッド株式会社 Terminal device, output method and computer program
JP7301681B2 (en) * 2019-09-05 2023-07-03 住友理工株式会社 Sleep state estimation method and sleep state estimation device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03277346A (en) * 1990-03-27 1991-12-09 Matsushita Electric Works Ltd Sleep deciding device
JPH0698863A (en) * 1992-09-21 1994-04-12 Matsushita Electric Ind Co Ltd Evaluating device for sleep
JP2000325315A (en) * 1999-05-19 2000-11-28 Arata Nemoto Sleep stage determinig method and sleep stage determining device
JP2004089267A (en) * 2002-08-29 2004-03-25 Ritsumeikan Sleeping depth estimation device and bedding equipped with the same
WO2004107978A1 (en) * 2003-06-03 2004-12-16 Cb System Co. Sleep stage judgment method and judgment device
JP2005118151A (en) * 2003-10-14 2005-05-12 Sanyo Electric Co Ltd Sleeping state estimating device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03277346A (en) * 1990-03-27 1991-12-09 Matsushita Electric Works Ltd Sleep deciding device
JPH0698863A (en) * 1992-09-21 1994-04-12 Matsushita Electric Ind Co Ltd Evaluating device for sleep
JP2000325315A (en) * 1999-05-19 2000-11-28 Arata Nemoto Sleep stage determinig method and sleep stage determining device
JP2004089267A (en) * 2002-08-29 2004-03-25 Ritsumeikan Sleeping depth estimation device and bedding equipped with the same
WO2004107978A1 (en) * 2003-06-03 2004-12-16 Cb System Co. Sleep stage judgment method and judgment device
JP2005118151A (en) * 2003-10-14 2005-05-12 Sanyo Electric Co Ltd Sleeping state estimating device

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011148793A1 (en) * 2010-05-28 2011-12-01 Sharp Kabushiki Kaisha Method and system for reliable respiration parameter estimation from acoustic physiological signal
CN105496362A (en) * 2015-12-15 2016-04-20 成都麦杰康科技有限公司 Sleep state monitoring system and method
US20220125376A1 (en) * 2019-02-13 2022-04-28 Kyoto University Sleep apnea syndrome determination apparatus, sleep apnea syndrome determination method, and sleep apnea syndrome determination program

Also Published As

Publication number Publication date
WO2009150744A1 (en) 2009-12-17
JPWO2009150744A1 (en) 2011-11-10

Similar Documents

Publication Publication Date Title
WO2009150765A1 (en) Sleeping condition monitoring apparatus, monitoring system, and computer program
US20220296176A1 (en) Processing biological data
JP6687263B2 (en) Biological information analysis device, system, program, and biological information analysis method
JP2011115188A (en) Sleeping condition monitoring apparatus, monitoring system, and computer program
KR101656611B1 (en) Method for obtaining oxygen desaturation index using unconstrained measurement of bio-signals
US7909771B2 (en) Diagnosis of sleep apnea
EP3334337B1 (en) Monitoring of sleep phenomena
CN106999065B (en) Wearable pain monitor using accelerometry
JP5632986B2 (en) Sleep stage determination device and sleep stage determination method
JP6513005B2 (en) Fatigue meter
US20200359909A1 (en) Monitoring device including vital signals to identify an infection and/or candidates for autonomic neuromodulation therapy
Ghosh et al. Detection of essential hypertension with physiological signals from wearable devices
JP7018621B2 (en) Biosignal analysis device, biosignal analysis method and biosignal analysis system
WO2004107978A1 (en) Sleep stage judgment method and judgment device
US20180055373A1 (en) Monitoring device to identify candidates for autonomic neuromodulation therapy
US10327648B2 (en) Blood vessel mechanical signal analysis
GB2469547A (en) Measurement of heart rate variability
JP2021535817A (en) Providing time information of the target person
EP4327743A1 (en) Sleeping state estimation system
JP6557489B2 (en) Biological state estimation device and computer program
JP6865438B2 (en) Method and device for creating indicators to determine neuropsychiatric status
Bassiouni et al. Combination of ECG and PPG signals for smart healthcare systems: Techniques, applications, and challenges
EP3708071A1 (en) Device, system, method and computer program for detecting atrial fibrillation
JP6927492B2 (en) Method and device for creating indicators to determine sleep disorders
US20230346234A1 (en) Wearable photoplethysmography device for detecting clinical decompensation based on heart rate variability

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 08874618

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 08874618

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP