CN103263260B - Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter - Google Patents

Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter Download PDF

Info

Publication number
CN103263260B
CN103263260B CN201310157508.8A CN201310157508A CN103263260B CN 103263260 B CN103263260 B CN 103263260B CN 201310157508 A CN201310157508 A CN 201310157508A CN 103263260 B CN103263260 B CN 103263260B
Authority
CN
China
Prior art keywords
sleep
bed
connecting shaft
test object
comb filter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201310157508.8A
Other languages
Chinese (zh)
Other versions
CN103263260A (en
Inventor
宋军
栗原阳介
渡边嘉二郎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING BOSHI LINKAG TECHNOLOGY Co.,Ltd.
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 CN201310157508.8A priority Critical patent/CN103263260B/en
Publication of CN103263260A publication Critical patent/CN103263260A/en
Application granted granted Critical
Publication of CN103263260B publication Critical patent/CN103263260B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

The present invention provides a kind of detecting system, in order to improve certainty of measurement, have employed exclusive technology and eliminates the external interference signals passing to bed from floor.The sensor arranged gathers physiology signal, comb filter is used to carry out data process, from obtaining heart rate number, Respiration Rate and standing up the signal waveform of frequency, extract the data of sleep stage continuously, move the state of data presumption sleep according to high-precision pulse and body.Speculate the sleep stage algorithm of Depth of sleep, use sleep stage sorting technique and exclusive sleep presumption technology, improve the accuracy rate of sleep stage presumption, thus the mensuration system of a kind of sleep index of detection accurately is provided.

Description

Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter
Technical field
The present invention relates to physiological detection signal supervisory instrument and detection signal is carried out statistical disposition, and according to the life extracted Reason signal speculates dormant sleep measurement system.
Background technology
Sleep has from the mental and physical important function of fatigue recovery, by assuring that high-quality sleep can promote to be good for The maintenance of health and strengthening.
For in order to carry out sleep quality measurement and without constraint, contactless in the case of measure the technology of pulse, Paper publishing and Japanese patent application has been carried out by present inventor.Such as, it has been suggested that a kind of sleep stage determining method and sentence Determining device, it utilizes the mattress detection Human Physiology information having enclosed air, it is possible to the sleep stage of discriminating test object Change, or calculate the time of each sleep stage.
According to described method, from the viewpoint of the sleep quality improving test object, in order to monitor and analyze sleep letter Breath, can be used for detecting physiology signal during test object sleep, uses comb filter to carry out data process, make this letter During breath is got up with relevant test object, the objective and subjective sleep parameters of mensuration is associated.The present invention proposes a germline System, it measures the sleep state of test object from the viewpoint of the sleep quality improving test object accurately.
Mattress is utilized to detect biological information in existing bed main body, by the change detection of the sleep stage of test object In the sleep stage determining method of sleep stage and decision maker, sometimes because of the factor in terms of external environment condition, big amount measurement data Destroyed, thus affect presumption.Owing to biological information signal is small, use hypersensitivity sensor, therefore, also can pick up survey The signal beyond signal that examination object sends, main incoming by ground or floor.
It addition, nearest mensuration environment, because being difficult to ensure that the reasons such as environment is consistent, sometimes it is set to 2 layers or more than 2 layers Bed, the vibration rank of bed is the highest sometimes.Accordingly, there exist the hash how got rid of beyond the Human Physiology information of testee Problem.
In anti-locking system is transmitted in vibration, the most also develop the vibrating sensor being incorporated with for detecting device with strong Property processed ground controls the system of the energy dynamic formula control device of the amplitude of vibration.
But, it is not reaching to practical due to larger-scale unit and the most high reason of cost.It addition, currently with The little wheel of elastomeric material of the anti-skidding cushioning effect having elastomeric material concurrently being directly mounted on pillar, it is difficult to as can be efficient Prevent the system of transmission that bed vibrates from playing a role.Exploitation is needed to prevent floor vibration on this device with low cost and realization For the purpose of transmission, although do not suppress the effect of the vibration occurred from device itself as described determinator forcibly, but It is prevented from floor vibration to transmit to bed main body, the Anti-vibration system that the vibration making determinator itself occur decays efficiently.
Summary of the invention
The present invention sets up in view of the above problems, its object is to, it is provided that a kind of detecting system, and it is without about Bundle, measure the multiple objective parameter relevant with the sleep quality testing object going to bed on bed under noncontact, prevent from passing through The pillar external vibration incoming from floor of bed, thus, eliminates the signal beyond physiology signal, with currently used bed Compare, it is possible to detect human body signal accurately.It addition, rely on sensor of the invention structure, and together with mattress combinations Sleep measurement device, it is possible to REM sleep and the cycle of non, the corresponding occurrence rate of Depth of sleep grade, rejecting The various information of the actual length of one's sleep of midway awakening, Sleep efficiency, apnea state etc carry out high accuracy and judge.As this Invention an embodiment, for the purpose of the impact of the floor vibration alleviating set place, as shown in figure 13, it is possible to from Floor vertically configures air spring, by vertical component and horizontal component with horizontal direction in the vibration transfer path of bed Vibration transmission suppress relatively low.
Described detecting system is provided with detection bed body vibration, the body motion detection device being made up of acceleration transducer.
Described detecting system is provided with detection with the body of h.d. test object is dynamic on bed, breathe, the body of heartbeat moves, The detection device being made up of position sensor.
In the present invention, high sensitivity pressure sensor is used to measure the pulse signal in sleep, high-precision from pulse signal The method that degree ground separation pulse component and body move component, uses comb filter.Move, to table according to the pulse separated and body The index of the index showing REM sleep state and the degree of depth representing sleep defines.It addition, devise according to the age distinguish each The calculating function of sleep stage average originating rate and standard deviation.The index of motion relevant the two sleep and use function presumption The algorithm of sleep stage.10 male sex are tested the sleep at 20 nights of object (22.2 years old mean age), with this motion The sleep stage that method calculates compares with the result of the sleep stage of R-K method, sleep level be divided into 6 stages (awakening, REM, Non-REM1,2,3,4), be divided into 5 stages (awakening, REM, Non-REM1,2,3/4), be divided into 3 stages (awakening, REM, Non-REM), time, the mean value of respective concordance rate is 51.6%, 56.2%, 77.5%.It addition, the mean value of k statistic is 0.29,0.39,0.48.
Currently, the international standard as the state for holding sleep uses R-K method.
The method based on E.E.G, eye movement, jaw muscle electricity, by sleep state classification be awakening, REM sleep, 6 stages such as Non-REM sleep 1,2,3,4.But, the method needs to measure E.E.G, eye movement, the electricity of jaw muscle electricity by being used for Very viscous head, face and the jaw of being affixed on is therefore, restrictive high.It addition, device is extensive and price is high, therefore, individual exists every day Family uses difficulty.Monitor that dormant basic research has multiple report so far without restrictions at home.These reports All in accordance with dynamic etc. and estimate sleep state by the pulse measured by the sensor without restricted type being arranged on bed and body.
On the one hand, the present invention provides a kind of detecting system, on the bed of test object, arranges at least more than one nothing Constraint, contactless sensor, detect that the above-mentioned test object h.d. in bed body relevant with sleep quality be dynamic, breathing, Multiple physiology signals of heartbeat etc, then use comb filter that described physiological signal is carried out data process, extract Along with dormant pulse signal and body movement signal.
On the other hand, the present invention provides a kind of detecting system, on the bed of test object, arranges at least more than one Without constraint, contactless sensor, detect that the relevant body of the sleep quality of above-mentioned test object h.d. in bed moves, exhales Multiple physiology signals of suction, heartbeat etc., then use comb filter that described physiological signal is carried out data process, extract Move with dormant heartbeat and body.
On the other hand, the present invention provides a kind of Depth of sleep monitoring system, and it includes aforesaid detecting system, also includes According to the heartbeat extracted by described signal extracting device, breathe and stand up each frequency signal waveform that body is dynamic, estimating continuously also Exporting the estimating device of the sleep stage of described test object, the variation figure of the sleep stage exported according to this estimating device is sentenced Determine Depth of sleep.
Accompanying drawing explanation
Fig. 1 air pressure mode measures the principle of pulse;
The human body signal that Fig. 2 measures in air pressure mode;
Frequency-the gain characteristic of Fig. 3 comb filter;
Fig. 4 uses the heart rate of comb filter group to estimate;
The average originating rate of each sleep stage at Fig. 5 each age and standard deviation;
The distribution in Fig. 6 REM sleep stage;
Fig. 7 awakening stage, the distribution of Non-REM sleep stage;
Fig. 8 measures system;
The heart rate number that Fig. 9 comb filter/FFT calculates respectively with the comparison of the heart rate number of ECG;
Figure 10 REM sleep index and δ ripple, spindle wave, the most oculomotor comparison;
Figure 11 Depth of sleep index and the comparison of δ ripple;
The comparison of Figure 12 sleep stage;
The example of the anti-vibration sensor signal acceptance agencies in Figure 13 footpost portion.
Detailed description of the invention
2. the embodiment of pneumatic mode
2.1 pneumatic mode
Fig. 1 represents the principle of the pneumatic mode of inventor's motion.Below the mattress of home-use bed (or below cotton-padded mattress), paving If the airbed of the ethene of thickness about 5mm.The internal pressure of airbed is the pressure identical with atmospheric pressure, and people lies in mattress Time upper, the dynamic air being broadcast in airbed by mattress of body that pulse and standing up etc. is formed.
The pressure change of this air is measured with hypersensitivity pressure sensor (P Co., Ltd. RIMO S11-M2).
Pressure sensor can measure the pressure oscillation of 0.2Pa-2Pa, has flat frequency characteristic with 0.1-3kHz [13], pressure sensor the signal measured is transfused to human body signal split circuit.
After passing through the bandpass filter of 5-10Hz in human body signal split circuit, implement full-wave rectification, envelope processing also Output.The heart rate number of the people in sleep is about 0.8Hz~1.5Hz, but there is harmonic wave in the waveform measured by pressure sensor Component, is filtered than of a relatively high 5~10Hz with S/N.
2.2 problem log
In the present invention, the motion using the pneumatic mode measurement signal with Fig. 1 is as follows.
(P1) according to output signal, the method that separation pulse component and body move component accurately.
(P2) definition represents the index of REM sleep.
(P3) definition represents the index of Depth of sleep.
(P4) letter calculating the average originating rate of each sleep stage, standard deviation defined according to the age bracket of healthy population Number.
(P5) algorithm of presumption R-K method sequentially sleep stage is proposed.
3. comb filter
Studying problem (P1), inventors, to sensor output data, implement FFT process, according to the base of pulse This ripple, harmonic component and other component, move pulse component component and separate [9-12] with body.But, with standing up etc. How containing the low frequency components near heart rate number during body is dynamic, therefore, the heart rate number asked for sometimes according to the peaks spectrum of FFT contains wrong Difference.
In the present invention, proposing there is a kind of method, it is in order to effectively utilize the harmonic component of pulse component, accurately Separate pulse and body moves, employ comb filter.
In body signal separates, the signal Δ t=0.01sec of output is sampled, and saves as the time shown in Fig. 2 Sequence data.The discrete time that 1min is spaced is set to k, k=1,2 ..., Tib。TibTest object couches and measures in bed Total time.The discrete time of the individual Δ t in the 1min of discrete time k is set to l, l=1,2 ..., N (=60/ Δ t).To Discrete time k, the pulse component of l are set to xk (l), body is moved component are set to nkTime (l), the output yk of body signal split circuit L () is moved the linear of component with pulse component and body as (1) formula and represents.
yk(l)=xk(l)+nk(l) (1)
As long as moving component according to (1) formula separable pulse component and body, it becomes possible to presumption heart rate number accurately.By pectination When the Dead Time of wave filter is set to T, τ=T/ Δ t in discrete time.It addition, the feedback oscillator of comb filter is set to g (0≤g≤1)。
At this moment, discrete time k, the pulse component of l and body move component as (2) formula, the output of the comb filter of (3) formula Give.
x ^ k ( l ) = y k ( l ) + y k ( l - τ ) + g · x ^ k ( l - τ ) - - - ( 2 )
n ^ k ( l ) = y k ( l ) - y k ( l - τ ) + g · n ^ k ( l - τ ) - - - ( 3 )
Fig. 3 (a), (b) represent the frequency-gain characteristic of the comb filter of (2) formula, (3) formula.
(2) under the frequency of DC component and the integral multiple of 1/T, there is peak value shown in the wave filter of formula such as Fig. 3 (a).Anti-mistake Coming, shown in the wave filter such as Fig. 3 (b) of (3) formula, under the frequency identical with the crest frequency of Fig. 3 (a), gain is 0.Fig. 3 (a), B the acuity of the peak value of the frequency characteristic of () changes according to feedback oscillator g, become sharp when close to 1.In (2) formula, (3) In formula, T and genuine heart rate cycle close in the case of, the wave filter of (2) formula makes the amplitude amplification of pulse component, in turn, (3) wave filter of formula makes pulse component decay, and makes body move aliquot amplification.Here, as shown in Figure 4, arranged in parallel has 78 extremely District time T is respectively 0.66s, 0.67s, 0.68s ..., the comb filter of (2) formula of 1.43s.1st comb filter dead District time T=0.66s, crest frequency is the integral multiple of 1.51Hz, therefore, with the resonance pulse of 90 times/min.Similarly, the 78th The Dead Time T=1.43s of individual comb filter, crest frequency is the integral multiple of 0.7Hz, therefore, with the pulse of 42 times/min Resonance.In terms of the heart rate number of every 1 point of kind, the resolution of the Dead Time of these comb filter is 60 Δ t/T2, 90 times/ It is 1.35 times/min during min, is 0.29 time/min when 42 times/min.Therefore, these 78 comb filter groups are at 42 times/min ~the pulse of average 0.82 time/min of resolution in the range of 90 times/min, can be covered.By observation signal ykL () input in parallel should Comb filter group, compares the output from each comb filterAmplitude.Presumption wherein output amplitude maximumThe Dead Time of comb filter, as the heart rate cycle.Then, to the Dead Time T estimated as this heart rate cycle, (3) output of formulaMake pulse component decay, on the contrary, make body move aliquot amplification.Thereby, it is possible to the size that putative aspect is dynamic.
4. sleep index
4.1 indexs being characterized with REM sleep
Problem (P2) is studied.
The comb filter definition of (2) formula of use is about the index of REM sleep.
REM sleep has a following feature [2,14]:
(1) represent and Non-REM1 and similar E.E.G of awakening.
(2) δ ripple, the occurrence frequency of spindle wave reduce.
(3) anxiety of antigravity muscle is wholly absent.
(4) eye movement rapidly occurs.
(5) Pulse Rate, Respiration Rate increase, and rhythm is irregular.It addition, blood pressure rises.
(6) in the case of adult, with the cycle of average 90~100min, REM sleep occurs.
(7) before and after REM sleep, body is dynamic concentrates.
In R-K method, shown in the judgement in REM sleep stage such as feature (1)~(4), emphatically for E.E.G, jaw muscle electricity, eyeball Motion.In the present invention, according to the feature (5) of pneumatic mode, it is conceived to that " rhythm of heart rate number becomes not advise when REM sleep Then ".
Data y by first half 30s of the 1min data of discrete time kk(1),yk(2),yk(3)…,yk(N/2) input Fig. 4 Comb filter group.At this moment outputAmplitude be maximum, the Dead Time being estimated as the heart rate cycle is set to Similarly, data y to the later half 30s of discrete time kk(N/2+1),yk(N/2+2),yk(N/2+3)…,yk(N), by amplitudeThe Dead Time becoming maximum is set toThus, the first half in the 1min of this discrete time k and later half heart rate number Change is usedRepresent.With k=1,2 ..., TibCarry out same process, move flat by q data of front and back All, ask for representing the index of the state of REM sleep.
R S I ( k ) = 1 2 q + 1 Σ i = - q q 60 | 1 / T k + i 1 - 1 / T k + i 2 | - - - ( 4 )
Knowable to the feature (5) of REM sleep, this index keeps big value when REM sleep.
4.2 indexs being characterized with awakening, Non-REM sleep
Problem (P3) is studied.Use (2) formula, the comb filter of (3) formula, define the finger about the degree of depth slept Mark.
Non-REM sleeps
(1) proportional to the degree of depth of sleep, the occurrence frequency of δ ripple uprises.
(2) in sleep stage Non-REM2, there is spindle wave.
(3) gradually deepening from waking state along with sleep, the size that body moves diminishes, and occurrence frequency tails off.
(4) Pulse Rate deepens and step-down along with sleep.
(5) Non-REM1 is sometimes after Non-REM3,4 or REM sleep occur, and moves along with big body and occurs.
Speciality according to pneumatic mode in the same manner as the situation of REM sleep, by feature (3) feature " along with sleep deepens, The size that body moves diminishes, and occurrence frequency tails off ", the depth index of definition sleep.
The output of the comb filter of (3) formula of the Dead Time T close with the heart rate cyclePulse component is made to decline Subtract, make body move aliquot amplification.The mean amplitude of tide of this output is set to by discrete time kDynamic the sending out of body in shallow sleep Raw frequency is many,Become bigger value, deep sleep becomes less value.
Sensor level output signal because of bedding and test object or sleeping position and different.In order to by this difference benchmark, incite somebody to action The sensor output y of discrete time kkL the mean amplitude of tide of () is set toWithRemoveAnd obtainIt addition,Having Significantly change in the case of body is dynamic and in the case of not having body dynamic, even if also having small change when moving almost without body Change.The big change caused to control this body to move, and small change is amplified, obtainLogarithm, as table Show the index of Depth of sleep.
S D I ( k ) = 1 2 log 2 ( P k n P k y ) - - - ( 5 )
This index becomes little value when deep sleep, becomes big value along with shoaling.
5. the presumption of sleep stage
The occurrence rate of the sleep stage that 5.1 ages are different
Problem (P4) is studied, determines each sleep stage at its age during the actual age of input test object Average occurrence rate and the function of outputting standard deviation.This function to 3 years old~92 years old of bibliography [15-18] 306 The occurrence rate of each sleep stage of people, determines with 4 curve approximations.
The occurrence rate of each sleep stage (6) formula represents, standard deviation (7) formula represents.(6) formula, (7) formula are to each sleep Stage represents respectively.Table 1 represents coefficient c0~c4And coefficient of determination R2
fx=c4a4+c3a3+c2a2+c1a+c0 (6)
δfx=c4a4+c3a3+c2a2+c1a+c0 (7)
Table 1 (6) formula, the coefficient of (7) formula and the coefficient of determination
X represents each sleep stage.The definition territory of a is set as becoming the age of 306 people on the data basis making this function Amplitude 3~92.(6) formula, (7) are even if number of times is set as more than 5 times by formula, and the coefficient of determination is the most almost unchanged, therefore, is set as 4 Secondary.Fig. 5 (a) represents the average occurrence rate of each sleep stage represented by (6) formula.
It addition, Fig. 5 (b) represents the standard deviation of each sleep stage average originating rate represented by (7) formula.
As shown in Fig. 5 (a), fWake+fREM+fNR1+fNR2+fNR3+fNR4=100 [%] sets up.
The occurrence rate that the occurrence rate of sleep increases along with the age and awakens rises, and the occurrence rate of contrary deep sleep reduces.
The presumption of 5.2 sleep stages
Problem (P5) is studied.
By REM sleep index RSI (k), (6) formula, the time of (7) formula presumption REM sleep.
The actual age of test object is substituted into (6) formula, a of (7) formula, seeks fREM,δfREM
At this moment the occurrence rate of REM sleep is set in fREM±δfREMIn the range of.
As shown in Figure 6, by REM sleep index RSI (k) by the order sequence declined, at fREM±δfREMScope in, will The time-bands of the RSI (k) that the body that maximum occurs keeps with the above when moving was allocated as the REM sleep time.This is to make use of REM The feature (7) of sleep.
Secondly, according to SDI (k), (6) formula, the distribution awakening of (7) formula, the time of Non-REM sleep 1,2,3,4.
The order of Depth of sleep index S DI (k) of the time-bands except REM sleep decline as shown in Figure 7 is sorted.
In the same manner as the situation of REM sleep, the actual age of test object is substituted into (6) formula, (7) formula,
Obtain awakening, the occurrence rate of each sleep stage in Non-REM1,2,3,4.
In the range of the occurrence rate of the most each sleep stage is set in as follows respectively:
Awakening: [fWake-δfNR1,fWake+δfWake]
Non-REM1:[fNR1-δfNR2,fNR1+δfNR1]
Non-REM2:[fNR2-δfNR3,fNR2+δfNR2]
Non-REM3:[fNR3-δfNR4,fNR3+δfNR3]
Non-REM4:[fNR4-δfNR3,fNR4+δfNR4]
Body in respective scope, as REM sleep, is moved big position and sets by the occurrence rate of awakening and Non-REM1 For its occurrence rate.The occurrence rate of Non-REM2,3,4, by SDI (k) by the order sequence declined, tilts big place and is set as Its occurrence rate, so distributes awakening, the sleep stage of Non-REM1,2,3,4.
After being assigned with each sleep stage, return to original time series, as R-K method sequentially sleep stage.
6. test checking
6.1 test object and test environments
Test object is with the health adult male sex 10 (A~J), 22.2 years old mean age as object, after obtaining informed consent, Measure 20 evenings.
Fig. 8 represents Measurement situation.
Use pneumatic mode pulse and the dynamic measurement of body.As pulse, the comparison use of sleep stage, traced by multiple fluctuation Device (SANYOFIT2500NEC tri-is flourish), measures ECG, E.E.G and eye movement simultaneously.
The checking of the validity of 6.2 comb filter
About the validity of problem (P1), to the heart rate number using comb filter presumption and the heart rate number estimated by FFT Compare checking.
The data of the 2nd (being later set to J-2) of test object J are inputted the comb filter group of Fig. 4 and estimates heart rate Number.The feedback oscillator of this comb filter group is g=0.95.
It addition, seek the heart rate number of 1 minute of J-2 with FFT.As heart rate number with reference to from the digital independent R-R of ECG Interval, seeks the average heart rate number of 1min.Fig. 9 (a) represents the heart rate number estimated by comb filter and the heart rate asked for by ECG Number, Fig. 9 (b) represents the heart rate number estimated by FFT and the heart rate number asked for by ECG.In 10min~100min, pectination filter Heart rate number presumption this side of result that ripple device is carried out is compared with the presumption result of FFT, little with the error of the heart rate number of ECG.Due to This time being mixed into of low-frequency noise such as moves with respiratory components and body, it is taken as that with having the amplification pectination to the characteristic of harmonic wave The presumption result of wave filter is correct.Data to all test objects, have estimated heart rate number with comb filter, FFT, ECG.Table The 2 square mean errors representing the heart rate number estimated by the FFT of comb filter and the heart rate number asked for by ECG.
In the presumption result of comb filter, there are 17 kinds of situation errors little in 20 kinds of situations.Pushing away in comb filter Determining in result, the mean value of error is also little than the error of FFT result.
The difference of the heart rate number of table 2 pectination ripple device and the ratio of the heart rate number difference of FFT
The checking of the appropriate property of 6.3REM sleep index
The validity of validation problem (P2).
Figure 10 represents that the δ ripple of E.E.G that the RSI (k) with C-3 and polygraph measure, spindle involve eye rapidly The incidence of ball motion.
(4) the rolling average number of times of the RSI (k) of formula is q=10.When RSI (k) keeps big value, the variation of heart rate number Greatly, therefore, compared with the feature of REM sleep (5), the possibility of REM sleep is high.
The time-bands of the grey of Figure 10 represent δ ripple, spindle wave incidence little, the most oculomotor occurrence number become Height, the feature of REM sleep is (2), (4) feature.
In identical time-bands, REM sleep index RSI (k) represents big value.
It addition, the big position of the value of RSI (k) occurred the feature of (6) of REM sleep occur with about 100 minute cycle.This Being to represent in REM sleep, the state of its feature (2), (4) and (5) occurs simultaneously, therefore, index RSI (k) can estimate REM The generation of sleep.
The checking of the appropriate property of 6.4 Depth of sleep indexs
Validity is verified for problem (P3).
(3) formula is set to the mean amplitude of tide of (5) formula of g=0.8 Obtained by the rolling average of 10 data.
Figure 11 represents SDI (k) and the δ ripple of C-3.
In the time-bands of the grey of Figure 11, the occurrence rate of δ ripple is high, knowable to Non-REM sleep characteristics (1), is positioned at deep Sleep state.In this time-bands, SDI (k) represents little value.
This be represent due to Non-REM sleep feature (1), the state of (3) occur simultaneously, therefore, according to sleep deep Degree index S DI (k) can estimate the degree of depth of sleep.
The presumption of 6.5 sleep stages
For the validity of (P5), use the algorithm of narration in the 5th chapter to judge sleep stage, sentence with by R-K method Fixed sleep stage compares.Figure 12 represents the RSI (k) of the data of C-3, SDI (k), the sleep stage of presumption and uses R-K The sleep stage that method determines.
The sleep stage that relatively estimated by, SDI (k) and during with the sleep stage of R-K method presumption, as being in REM sleep spy Levy the happening part of the rhythm in about 100 minute cycle as (6) or towards the appearance etc. that dawn, direction sleep stage shoaled Overall shape is similar.But, 100 minutes, near 280 minutes, be judged to the feelings of Non-REM4 relative to the method proposed Condition, is judged to Non-REM2 in R-K method.
It addition, near 470min, the method for proposition is judged to Wake, but is judged to REM in R-K method.To with carrying Sleep stage that the method that goes out judges, the sleep stage judged by R-K method respectively with 6 stages (awakening, REM, Non-REM1, 2,3,4), by Non-REM3, (awakening, REM, Non-REM1,2,3/4) [19], 3 stages in 45 stages united two into one (awakening, REM, Non-REM) in compare.
Comparative approach uses TibThe ratio (consistent number/T of the consistent number of the result of determination of two kinds of methodsib)*100 The concordance rate of [%], k statistic.Table 3 represents a when sleep stage entirely testing object is divided into 6 stages, 5 stages, 3 stage Cause rate k statistic, respective mean value and standard deviation.
The mode of table 3 motion and the concordance rate of the sleep section of R-K method and κ are measured
The sleep stage of the sleep=537min of the one of Figure 12 be concordance rate during 6 stage be 52.1%, k statistic is 0.25, sleep stage be concordance rate during 5 stage be 58.7%, k statistic is 0.38, consistent when sleep stage was 3 stage Rate is 79.5%, and k statistic is 0.46.
7. as the test object of C-3, the concordance rate of each sleep stage number of the test object of table 3 and k statistic with Number of stages to tail off and rise.
In the case of sleep stage was 6 stages, relative to average the 36.4% of the concordance rate of document [12], concordance rate Mean value rises to 51.6%.
The mean value of k statistic at this moment is 0.29, also includes the accidental error unanimously caused.
It addition, when sleep stage was 5 stage, rise 4.6% when the mean value of concordance rate is 56.2%, in contrast, 3 stages were 77.5%.
This is because inconsistent rate between Non-REM1, Non-REM2, Non-REM3/4 than Non-REM3 and The inconsistent rate of the judgement between Non-REM4 is many.
K statistic is 0.39 when 5 stage, is 0.48 when 3 stage, it is possible to hold figure and the REM of the overall situation of sleep Sleep, the rhythm of Non-REM sleep.
According to method as above, can be used for, during test object constraint, nothing invasion and attack sleep, detecting and testing object The relevant body of sleep quality dynamic, breathe, the health bio signal of heartbeat etc, use comb filter to carry out at data Reason, the parameter of the objective and subjective sleep measured during making this information and getting up about test object is associated.
From the viewpoint of the sleep quality improving test object, need to measure accurately the sleep of test object.But It is especially to implement can carrying out in the measuring chamber that above-mentioned data obtain of soundproof, vibrationproof, but in general society It is economically relatively difficult.Therefore, when the bedroom that setting place is general ward, family, that be utilized at present, straight Connect in the vibration proof transmission such as rubber rumble of the anti-skidding buffer of the elastomeric material being installed on pillar, mobile wheelchair, for coming From the vibration on floor, prevent vibration transmission and measure with carrying out high-quality the most difficult.
Formation, in view of above-mentioned problem, is prevented the external vibration of the pillar of incoming bed from floor from carrying out the knot transmitted by the present invention The vibration proof transfer device of structure, is set as constituting as follows.
Such as, motion has a kind of vibration proof transfer device, and it is on the bed using the test object being arranged in nothing sleep Multiple sensors, detect the physiology signal of h.d. in bed, use comb filter to carry out data process, extract companion Along with in the detecting system that dormant pulse and body are dynamic, from floor in the vibration transfer path of bed vertically and Horizontal direction configuration vibration proof transferring element, cuts off the vibration transmission of vertical component and horizontal component.
Figure 13 represent the 1st scheme of vibration proof transfer device examples of implementation.In actual device, use multiple (generally 4) buffer cell, but this figure is one side view.
In this configuration example, in order to prevent in the transmission of the floor vibration arranging position, it is generally configured with: for elastic bearing The air spring 304 of the Z-direction (above-below direction) of device 1, for supporting the vibration proof mechanism of this vertical direction and being fixed on bed and prop up Connecting shaft 302 on post, fixing air spring 304 and in buffer moving range, keep the urceolus of movement of horizontal direction 303 and in order to improve the damping characteristics of the Z-direction of air spring 304 and and the buffer 305 that sets.This buffer 305 is consolidated Due to the bottom of device, and when the high vibrationproof for the buffering carrying out Z-direction, urceolus 303 is fixed on floor 300.
During test object bunk bed, owing to one end applies the load of more than behaviour in service, therefore, become connecting shaft 302 One end 302a abuts with urceolus 303, and air spring 304 will not apply the load of more than necessity.On the other hand, do not connect at bed pillar In the case of ground, the other end 302b of the groove of connecting shaft 302 and urceolus 303 still abut and will not expand to required scope with Outward, this just makes air spring 304 keep intact.
The pneumatics P0kg/cm**2 of air spring 304 is body weight W1kg and the weight W2kg phase of bed entirety of test object Add, for W0kg, if the bed of 4 supportings, even if load partially, be the most also 1/4, setting W0/4kg air spring 304 During contact area S cm**2, available pneumatics P=W0/4S kg/cm**2 calculates.Pass in vibration as near this calculating Pass flexible buffer flexible in path, be set as vertically configuring with horizontal direction, it is possible to by vertical component and level The vibration transmission suppression of component is relatively low flexible buffer, is generally effective.
2nd scheme, it is characterised in that implement vibration proof transferring element, because spring maintains the capable of expansion and contraction of a pillar, Via being installed on the pulley of a pillar and hanging from bed main body with rope, and can vertically up-down vibration, the vibration of vertical component The stretching spring being installed on rope flexible in absorbed, the vibration of horizontal component by the rope of hanging as put Movement in absorbed.By lengthening the parts for pendency, it is also possible to suppress to be suitable by the intrinsic vibration number of horizontal direction Low.
The vibration proof transfer device of the third aspect, it is characterised in that vibration proof transferring element uses magnet, and it is from floor The most vertically contrary magnet or electromagnet is configured with horizontal direction in the vibration transfer path of bed, anti- The only vibration transmission of vertical component and horizontal component.
The vibration proof transfer device of fourth aspect, it is characterised in that vibration proof transferring element makes to be provided with floating material Bed main body floats over a pillar or is set directly in the tank that there was added liquid on floor, with tank and described bed main body not phase The mode connect installs holding member, makes liquid absorb vibration between bed main body and tank.
The vibration proof transfer device of the 5th aspect, it is characterised in that from the vibration on ground through footpost conduction to bed body, Under footpost, installation one has the axis of guide of sliding eye as damper, and when being subject to the pressure from bed body, bed body can guide Sliding in axle, equipped with elastomer (air packing ring) in the axis of guide, the position horizontally and vertically of bed body is moved and is all added Being pressed on elastomer (air packing ring), elastomer (air packing ring) is equivalent to the compression spring of passive vertical telescopic, Added with lubriation material between the axis of guide and elastomer.

Claims (8)

1. a detecting system, uses on the bed being arranged on test object, and at least one is without constraint, contactless sensor, inspection The body relevant with sleep quality measuring the described test object of h.d. in bed moves, breathes and multiple human bodies of heartbeat Physiological signal, it is characterised in that use the comb filter group being made up of multiple comb filter that physiological signal is carried out data Process, it is determined that the resonance pulse signal corresponding with described comb filter Dead Time, thus extract accurately with sleeping The pulse signal of dormancy state and body movement signal;Described detecting system includes vibration proof transfer device, described vibration proof transfer device Including buffer, urceolus, air spring and connecting shaft, described buffer has inwall, and described inwall forms concave station;Described urceolus Being arranged on concave station, the outer surface of described urceolus matches with the inner surface of described concave station;Described air spring is arranged on urceolus In, with the contact internal walls of urceolus;Described connecting shaft includes connecting shaft, base and lower connecting shaft, and described upper connecting shaft is arranged on On the upper surface of base, described lower connecting shaft is arranged on the lower surface of base, and described upper connecting shaft contacts with bed pillar, institute The upper surface of the lower surface and described air spring of stating base contacts, described lower connecting shaft and the inner surface of described air spring Contact.
2. detecting system as claimed in claim 1, it is characterised in that when using mattress in bed, described detecting system sets There are the internal pressure detecting this mattress, the body motion detection device being made up of pressure sensor.
3. detecting system as claimed in claim 1, it is characterised in that described detecting system is provided with detection bed body vibration, by The body motion detection device that acceleration transducer is constituted.
4. detecting system as claimed in claim 1, it is characterised in that described detecting system is provided with detection with h.d. on bed The body of test object is dynamic, breathe, the body of heartbeat moves, the detection device being made up of position sensor.
5. a detecting system, uses on the bed being arranged on test object, and at least one is without constraint, contactless sensor, inspection Survey the body relevant with the sleep quality of described test object h.d. in bed to move, breathe and multiple Human Physiologies of heartbeat Signal, it is characterised in that use the comb filter group being made up of multiple comb filter that physiological signal is carried out data process, Judge the resonance pulse signal corresponding with described comb filter Dead Time, thus extract accurately with sleep state Pulse signal and body movement signal, and estimate the heart rate cycle;Described detecting system includes vibration proof transfer device, described vibration proof Transfer device includes that buffer, urceolus, air spring and connecting shaft, described buffer have inwall, and described inwall forms concave station; Described urceolus is arranged on concave station, and the outer surface of described urceolus matches with the inner surface of described concave station;Described air spring sets Put in urceolus, with the contact internal walls of urceolus;Described connecting shaft includes connecting shaft, base and lower connecting shaft, described upper connection Axle is arranged on the upper surface of base, and described lower connecting shaft is arranged on the lower surface of base, described upper connecting shaft and bed pillar Contacting, the lower surface of described base contacts with the upper surface of described air spring, described lower connecting shaft and described air bullet The inner surface of spring contacts;Carry out the setting not on the sensor that the sleep generation of test object affects, be used for monitoring life Reason parameter and the multiple described sensor of ambient parameter.
6. a Depth of sleep monitoring system, it includes the detecting system described in claim 1 or 5, also includes according to by signal Heartbeat that extraction element extracts, breathe and stand up each frequency signal waveform that body is dynamic, estimate continuously and export described survey The estimating device of the sleep stage of examination object, judges sleep according to the variation figure of the sleep stage exported by described estimating device The degree of depth.
7. Depth of sleep monitoring system as claimed in claim 6, it is characterised in that join between health and the bed of test object Put the elastomer mattress built with sensor, from the output signal extraction heart of the body motion detection device that the body of detection test object moves The dirty signal beating, breathe and standing up body dynamic frequency, judges Depth of sleep according to the evolutionary mode of sleep stage.
8. Depth of sleep monitoring system as claimed in claim 6, it is characterised in that go to bed by mattress structure at test object Time on the elastomer bed become, Depth of sleep monitoring system is provided with the pressure sensor of the internal pressure detecting this mattress and is constituted Body motion detection device.
CN201310157508.8A 2013-05-02 2013-05-02 Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter Active CN103263260B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310157508.8A CN103263260B (en) 2013-05-02 2013-05-02 Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310157508.8A CN103263260B (en) 2013-05-02 2013-05-02 Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter

Publications (2)

Publication Number Publication Date
CN103263260A CN103263260A (en) 2013-08-28
CN103263260B true CN103263260B (en) 2016-09-07

Family

ID=49006960

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310157508.8A Active CN103263260B (en) 2013-05-02 2013-05-02 Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter

Country Status (1)

Country Link
CN (1) CN103263260B (en)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103263261B (en) * 2013-05-02 2016-09-07 宋军 Without constraint physiological parameter measurement method presumption sleep index and the system of sleep stage
CN104107037A (en) * 2014-07-15 2014-10-22 北京博实联创科技有限公司 Physiological information acquiring and processing system
CN104382574B (en) * 2014-11-06 2017-09-29 深圳市维亿魄科技有限公司 A kind of method and system based on pulse wave data monitoring sleep quality
CN104545870B (en) * 2015-01-23 2017-01-11 苏州本草芙源医疗设备有限公司 Heart rate detection mattress and method
CN105559751B (en) 2015-12-14 2018-06-12 安徽华米信息科技有限公司 Monitor the method, apparatus and wearable device of micro- active state
CN105595973A (en) * 2016-01-08 2016-05-25 四川东鼎里智信息技术有限责任公司 Sleeping abnormity warning device
CN105615836A (en) * 2016-01-08 2016-06-01 四川东鼎里智信息技术有限责任公司 Parasomnia monitoring system
CN107307846A (en) * 2016-04-27 2017-11-03 南京理工大学 Contactless sleep stage method
CN106562762B (en) * 2016-08-10 2019-10-25 深圳市格兰莫尔科技有限公司 A kind of sleep monitor and improve device and method
CN106562761B (en) * 2016-08-10 2020-04-10 深圳市格兰莫尔科技有限公司 Physiological signal acquisition device and method and sleep monitoring and improving device
WO2019061081A1 (en) * 2017-09-27 2019-04-04 深圳和而泰智能控制股份有限公司 Method, apparatus, and device for monitoring physiological information and smart pad
WO2019075757A1 (en) * 2017-10-20 2019-04-25 深圳和而泰智能控制股份有限公司 Sleep body motion monitoring method, apparatus and sleep body motion monitor
CN108042108B (en) * 2017-12-06 2020-12-08 中国科学院苏州生物医学工程技术研究所 Sleep quality monitoring method and system based on body vibration signals
CN109199355B (en) * 2018-09-18 2021-09-28 深圳和而泰数据资源与云技术有限公司 Heart rate information detection method and device and detection equipment
JP2020054782A (en) * 2018-09-26 2020-04-09 日本電信電話株式会社 Biological information analysis device, biological information analysis method, and biological information analysis system
CN109924961A (en) * 2019-04-01 2019-06-25 中国医学科学院生物医学工程研究所 Attitude adjusting method, device, terminal and the storage medium of nursing bed
TWI726352B (en) * 2019-07-09 2021-05-01 國立臺灣科技大學 Non-direct contact method for monitoring physiological and activity signal

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3616611B1 (en) * 2006-06-01 2020-12-30 ResMed Sensor Technologies Limited Apparatus, system, and method for monitoring physiological signs
JP2008188079A (en) * 2007-02-01 2008-08-21 Matsushita Electric Ind Co Ltd Biological signal detecting device and sleeping apparatus using the same
CN201090656Y (en) * 2007-06-01 2008-07-23 成都飞机工业(集团)有限责任公司 Air spring for automobile
JP5403161B2 (en) * 2010-06-21 2014-01-29 アイシン精機株式会社 Biological information detection system
CN102415879A (en) * 2011-09-09 2012-04-18 北京大学深圳研究生院 Sleep monitoring device based on piezoelectric film sensor
CN102697488B (en) * 2012-06-28 2013-11-20 黑龙江大学 Unrestrained sleep physiological information detection device
CN104107037A (en) * 2014-07-15 2014-10-22 北京博实联创科技有限公司 Physiological information acquiring and processing system

Also Published As

Publication number Publication date
CN103263260A (en) 2013-08-28

Similar Documents

Publication Publication Date Title
CN103263260B (en) Use physio-parameter detection system and the Depth of sleep monitoring system of comb filter
CN103263261B (en) Without constraint physiological parameter measurement method presumption sleep index and the system of sleep stage
EP3078948B1 (en) Acoustic and vibration information accumulation mechanism, acoustic and vibration sensing system, and computer program
Paalasmaa et al. Adaptive heartbeat modeling for beat-to-beat heart rate measurement in ballistocardiograms
Chee et al. Air mattress sensor system with balancing tube for unconstrained measurement of respiration and heart beat movements
JP4611206B2 (en) Fatigue degree measuring device, fatigue detecting device, and computer program
EP2783631B1 (en) Biological status estimation device and computer program
CN101657156B (en) Ballistocardiographic sensor system with a sensor arrangement and method of ballistocardiographic detection of body movements
CN105769122B (en) A kind of measuring apparatus of apnea hypopnea indexes
CN102088904B (en) Ballistocardiogram analysis method and device
EP2862507B1 (en) Physiological state-analyzing device and computer program
CN107082034A (en) A kind of intelligent automobile seat cushion and its measuring method
KR101628262B1 (en) Method and apparatus for heart rate calculation using frequency analysis
CN108289638A (en) Biological state determination device and organism state determination method
JP5352814B2 (en) Autonomic nerve component index estimating apparatus and autonomic nerve component index estimating method
JP3877615B2 (en) Sleep depth estimation device
Chow et al. Respiratory monitoring using an air-mattress system
CN108354610A (en) A kind of Falls Among Old People detection method and detecting system based on three-axis sensor and EGC sensor
WO2004107978A1 (en) Sleep stage judgment method and judgment device
JP2000325315A (en) Sleep stage determinig method and sleep stage determining device
TWI629049B (en) A method for analyzing a heart shock signal for calculating a short-term heart rate value
Sloboda et al. A simple sleep stage identification technique for incorporation in inexpensive electronic sleep screening devices
CN209966365U (en) Portable sleep monitoring system
CN109805909A (en) Human body physical sign measuring system and method when a kind of seat
CN103349547B (en) For the vibration proof transfer device of physio-parameter detection system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: SONG JUN

Free format text: FORMER OWNER: BEIJING BOSHI LINKAGE TECHNOLOGY CO., LTD.

Effective date: 20141114

Free format text: FORMER OWNER: SONG JUN KURIHARA YOSUKE WATANABE KAJIRO

Effective date: 20141114

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100193 HAIDIAN, BEIJING TO: 100085 HAIDIAN, BEIJING

TA01 Transfer of patent application right

Effective date of registration: 20141114

Address after: 100085 Beijing city Haidian District Qinghe woolen Road No. 26 Building No. 1, 1B1401 language Qing Yuan

Applicant after: Song Jun

Address before: 100193, room 2, building 2344, incubator of Beijing Zhongguancun Software Park, Haidian District, Beijing, Northeast China

Applicant before: Beijing Boshi Linkage Technology Co., Ltd.

Applicant before: Song Jun

Applicant before: Kurihara Yosuke

Applicant before: Watanabe Kajiro

C14 Grant of patent or utility model
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Physiological parameter detecting system with comb filter and monitoring system for depth of sleep

Effective date of registration: 20180814

Granted publication date: 20160907

Pledgee: Zhongguancun Beijing technology financing Company limited by guarantee

Pledgor: Song Jun

Registration number: 2018990000687

PC01 Cancellation of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20210323

Granted publication date: 20160907

Pledgee: Zhongguancun Beijing technology financing Company limited by guarantee

Pledgor: Song Jun

Registration number: 2018990000687

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210408

Address after: 100193203, 2 / F, building 10, Zhongguancun Software Park, 8 Dongbeiwang West Road, Haidian District, Beijing

Patentee after: BEIJING BOSHI LINKAG TECHNOLOGY Co.,Ltd.

Address before: 1b1401, building 1, Yanyu Qingyuan, No.26, Qinghe Maofang Road, Haidian District, Beijing 100085

Patentee before: Song Jun