CN113974568B - Slope interference-free method for calculating metabolic rate of sleep process - Google Patents
Slope interference-free method for calculating metabolic rate of sleep process Download PDFInfo
- Publication number
- CN113974568B CN113974568B CN202111322781.2A CN202111322781A CN113974568B CN 113974568 B CN113974568 B CN 113974568B CN 202111322781 A CN202111322781 A CN 202111322781A CN 113974568 B CN113974568 B CN 113974568B
- Authority
- CN
- China
- Prior art keywords
- sleep
- steady state
- state
- weight
- calculating
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 230000037323 metabolic rate Effects 0.000 title claims abstract description 45
- 230000008569 process Effects 0.000 title claims abstract description 25
- 238000005259 measurement Methods 0.000 claims abstract description 44
- 239000013585 weight reducing agent Substances 0.000 claims abstract description 22
- 230000037396 body weight Effects 0.000 claims abstract description 21
- 238000012216 screening Methods 0.000 claims abstract description 9
- 230000008030 elimination Effects 0.000 claims abstract description 5
- 238000003379 elimination reaction Methods 0.000 claims abstract description 5
- 230000004580 weight loss Effects 0.000 claims description 9
- 230000004060 metabolic process Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 abstract description 11
- 230000002503 metabolic effect Effects 0.000 abstract description 6
- 230000013632 homeostatic process Effects 0.000 abstract description 3
- 230000009467 reduction Effects 0.000 abstract description 3
- 230000008859 change Effects 0.000 description 10
- 238000005303 weighing Methods 0.000 description 7
- 208000022249 Sleep-Wake Transition disease Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002688 persistence Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000031648 Body Weight Changes Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000004579 body weight change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000035622 drinking Effects 0.000 description 1
- 239000003651 drinking water Substances 0.000 description 1
- 235000020188 drinking water Nutrition 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000004622 sleep time Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4866—Evaluating metabolism
Abstract
The invention relates to a calculation method of metabolic rate of sleep process based on slope interference elimination, comprising measuring weight of a bed once every first preset time interval, identifying steady state according to the measured value, and screening effective sleep steady state according to slope and duration of each steady state in sleep period; calculating the sum of weight reduction values during each effective sleep homeostasis period as an absolute weight reduction value during sleep; sleep metabolic rate was calculated using the absolute decrease in body weight during sleep. In the invention, the metabolic quantity is calculated by measuring the absolute reduction of the weight of the person in the period of sleeping in the bed and the weight reduction in the sleeping process, and the metabolic condition of the person can be reflected more truly because the weight measurement is more direct than the heat measurement; by adopting the method for calculating the weight reduction value of each stable state duration, each interference factor affecting the weight in the sleeping process can be eliminated, and the stable state containing interference can be removed by setting a slope threshold to remove the stable state with larger slope deviation.
Description
Technical Field
The invention belongs to the technical field of sleep process metabolic rate monitoring, and relates to a slope interference-free calculation method of a metabolic rate of a sleep process.
Background
The metabolic rate represents the energy consumed by the human body per unit time. In the prior art, the metabolic rate is generally calculated by using the metabolism of heat, but the monitoring of the heat is difficult, and the heat can be converted only by other indexes, so that a large error exists in a calculation result.
Because people consume energy at all times, under the condition of not drinking water and not eating food, the weight is gradually reduced, so the metabolic rate can be calculated through the change of the weight, but the change value of the weight is very small, the traditional weight measurement means cannot accurately monitor the reduced value, and the metabolic rate can also generate larger fluctuation due to the movement quantity of each period of time in the process of activities. However, the metabolism rate of a person is basically stable during sleep, and if the body weight during sleep can be monitored with high accuracy, various interference factors affecting the body weight during sleep are eliminated, and the absolute reduction of the body weight during the period of sleeping of the person in a bed is measured, the metabolism rate can be calculated through the reduction of the body weight during sleep.
Disclosure of Invention
Accordingly, the present invention is directed to a method for calculating a metabolic rate of a sleep process based on slope disturbance elimination.
In order to achieve the above purpose, the present invention provides the following technical solutions:
a method for calculating metabolic rate of sleep process based on slope disturbance removal, comprising the following steps:
step S1, a first preset time T is set every interval 0 Measuring the weight of the primary bed;
step S2, defining the length as uT 0 The time period of the time frame is a long time window, whether the long time window taking the current measurement time as the end time is in a stable state is judged according to the measured value of the weight in the long time window, and if the long time window is in the stable state, the current state is judged; recognizing the starting time and the ending time of sleep according to the stable states, and calculating the duration time and the slope of each stable state in the sleep period;
step S3, screening effective sleep stable states according to the slope and duration time of each stable state in the sleep period, and calculating the sum of weight reduction values in each effective sleep stable state period as an absolute weight reduction value D in the sleep period;
and S4, calculating the sleep metabolism rate by using the value of the absolute weight reduction during the sleep period.
Further, in the step S1, a plurality of pressure sensors are arranged under the bed, and each pressure sensor is arranged at intervals of a first preset time T 0 Measuring the weight of the primary bed; defining total weight S as the sum of the measured values of the pressure sensors, and defining empty bed weight B as the empty bed weight of the pressure sensorsThe sum of the measured values defines the instantaneous body weight w=s-B;
in the step S2, the method for identifying the sleep start time is as follows: judging whether a loading event occurs according to the measured total weight S; and judging whether to enter a steady state after the occurrence of the loading event, if so, determining to enter a sleep state, and taking the starting time of the steady state as the starting time of sleep.
Further, the step S2 includes the following substeps:
step S201, judging whether a loading event occurs according to the measured value of the pressure sensor; executing step S202 if a get-in event occurs;
step S202, detecting whether a steady state is entered and whether a getting-out event occurs, and judging that the sleeping state is entered if the steady state is entered; step S203 is performed; returning to execute step S201 if a get-off event occurs;
step S203, detecting whether the steady state is ended, if so, executing step S204;
step S204, calculating the duration and the slope of the steady state;
step S205, whether a steady state is re-entered or whether a getting-out event occurs is detected, and if the steady state is re-entered, the step S204 is executed again; executing step S206 if a get-off event occurs;
step S206, detecting whether a loading event occurs within the threshold of the loading and unloading time and entering a steady state, and if the loading event occurs within the threshold of the loading and unloading time and entering the steady state, returning to execute step S204; otherwise, step S207 is performed;
step S207, determining that sleep is completed, and executing step S3.
Further, the method for judging whether the long time window is in a stable state is as follows: calculating the current measurement time as the end time and the length of 2uT 0 Average S of all total weights S over a period of time TW By means of average number S TW Calculating a standard deviation of all total weight S values recorded over the long time window for the total weight mean, the standard deviation being defined as σ TWC Setting a long steady-state standard deviation threshold delta 1 When a long time window ends, if sigma of the long time window TWC ≤δ 1 Judging that the long time window is in a stable state; if sigma TWC ＞δ 1 The long time window is determined to be in an unstable state.
Further, the method for judging whether the event of getting in or out of bed occurs is as follows: if a certain measuring moment of the sensor is changed from an empty bed state to a bed state, judging that a loading event occurs at the measuring moment; if a certain measuring moment of the sensor is changed from the in-bed state to the empty-bed state, the occurrence of an out-of-bed event at the measuring moment is judged.
Further, the slope of the steady state is calculated by linearly fitting the values of all the instantaneous weights W in the steady state by using a least square method.
Further, the weight loss value DeltaW during steady state was calculated using the following formula a ：
ΔW a ＝-K×ΔT W
Wherein DeltaT W Indicating the duration of the steady state.
Further, the method for screening the effective sleep stable state comprises the following steps:
defining a slope threshold and a steady state duration threshold, firstly calculating the average value of slopes of all steady states during sleep, and then removing the steady state with the absolute value of the difference between the slope and the average value being larger than the slope threshold and the steady state with the positive slope; and removing the stable state with the duration less than the threshold value of the stable state duration, and taking the rest stable state as the effective sleep stable state.
Further, the method of calculating the sleep metabolic rate using the absolute weight loss value D during sleep includes:
the sleep metabolic rate m is calculated according to the following formula:
wherein W is s1 Indicating the body weight at the beginning of sleep.
Further, the method of calculating the sleep metabolic rate using the absolute weight loss value during sleep further includes:
the 8-hour daily sleep metabolic rate m (8 h) was calculated according to the following:
where Δt represents the duration of sleep.
In the present invention, by measuring the absolute decrease in weight of a person during a period of sleeping in a bed, the sleep metabolic rate can be calculated from the decrease in weight during sleep, and the metabolic level of the person can be reflected by long-term statistics, as with the caloric metabolic rate. Since the measurement of body weight is more direct than the measurement of calories, the metabolic situation of a person can be reflected more realistically. By introducing the steady state and calculating the weight reduction value of each steady state duration, each interference factor affecting the weight in the sleeping process can be eliminated, and the weight measurement is more accurate. By setting the slope threshold to remove steady state where the slope deviates greatly, steady state containing disturbances can be removed.
Drawings
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in the following preferred detail with reference to the accompanying drawings, in which:
FIG. 1 is a flow chart of a preferred embodiment of the method for calculating metabolic rate of sleep process based on slope disturbance elimination according to the present invention.
FIG. 2 is a flow chart for identifying the start and end times of sleep from the steady states and calculating the duration and slope of each steady state during sleep.
Fig. 3 is a graph showing the change in weighing readings for 8 hours.
Fig. 4 is a schematic representation of the details of the change in weighing readings for 8 hours.
Figure 5 is a graphical representation of steady state body weight change over 8 hours.
Detailed Description
Other advantages and effects of the present invention will become apparent to those skilled in the art from the following disclosure, which describes the embodiments of the present invention with reference to specific examples. The invention may be practiced or carried out in other embodiments that depart from the specific details, and the details of the present description may be modified or varied from the spirit and scope of the present invention. It should be noted that the illustrations provided in the following embodiments merely illustrate the basic idea of the present invention by way of illustration, and the following embodiments and features in the embodiments may be combined with each other without conflict.
As shown in fig. 1, a preferred embodiment of the method for calculating the metabolic rate of the sleep process based on slope disturbance removal according to the present invention comprises the steps of:
step S1, arranging a plurality of pressure sensors under the bed, wherein each pressure sensor is arranged at intervals of a first preset time T 0 The weight of the primary bed was measured.
Because the weight change generated by metabolism is very weak in the sleeping process, the measurement is difficult, and uncertain interference factors exist in the sleeping process; in the embodiment, n (n is more than or equal to 4) pressure sensors are arranged under the bed, preferably four pressure sensors are respectively padded under four bed legs, of course, the pressure sensors can also be padded at the contact positions of the bed board and the bed frame, and the number of the pressure sensors can also be more than four; the measurement accuracy of the pressure sensor should be not less than 10g, preferably 1g. Every first preset time T of each pressure sensor 0 (T 0 Typically set to 0.1 to 2s, preferably 1 s) to measure the weight once. In the whole measurement process, the original absolute output values of the pressure sensors are set as S1, S2, … … and Sn, and the sum of the S and the Sn is S, which is called total weight; the weight output values measured by the pressure sensors when the bed is empty are set as B1, B2, … … and Bn, and the sum of the weight output values is B, and the weight output values are called empty bed weight; (the value of the empty bed weight can be updated at regular intervals, so that the removed bed product is not usually taken or calculated). Instantaneous body weight w=s-B is defined.
When the accuracy of the sensor is high, the measured value of the sensor is not changed even in a static state, and the measured result is not remained at a constant value, so that a stable reading cannot be obtained; thus, the present embodiment introduces the concept of steady state to determine whether the measured value of the weight of the bed is stable, and the concept of steady state and steady state weight to measure the weight value.
And S2, recognizing a steady state according to the measured value of the pressure sensor, recognizing the starting time and the ending time of sleep according to the steady state, and calculating the duration and the slope of each steady state in the sleep period. The steady state can be identified by directly judging according to the value of the total weight S, or calculating the instant weight W according to the total weight S, and then judging according to the value of the instant weight W.
It should be noted that the steady state and the steady state in the present embodiment are completely different, and the steady state refers to a state in which fluctuation of the measured value is small, and is a state; and steady state refers to a period of time in steady state, which is a period of time.
The method for judging the steady state comprises the following steps:
define a length of uT 0 Is a long time window, where u is a natural number, and the value of u is preferably u=30 (i.e., uT 0 For half a minute, u may take other values such as 60, etc.), and it is determined whether or not the long time window with the current measurement time as the end time is in a steady state, and if the long time window is in a steady state, it is determined that the current measurement time is in a steady state.
The following method can be used to determine whether the long window is in a steady state: calculating the current measurement time as the end time and the length of 2uT 0 The average S of all total weights S over a period of time (i.e. twice as long window, in this example one minute) TW By means of average number S TW Calculating a standard deviation of all total weight S values recorded over the long time window for the total weight mean, the standard deviation being defined as σ TWC Setting a long steady-state standard deviation threshold delta 1 (δ 1 The value of (a) is generally in the range of 5 to 1000, preferably 200 in this embodiment), when a long window ends, if sigma of the long window TWC ≤δ 1 Judging that the long time window is in a stable state; if sigma TWC ＞δ 1 The long time window is determined to be in an unstable state. The standard deviation of this embodiment uses a window of twice as long to calculate the mean,the situation that the reading value changes greatly but the standard deviation changes little when the long time window is alternated, so that inaccurate identification can be avoided.
In this embodiment, in order to avoid a situation where the change in the reading value is large but the change in the standard deviation is not large just when the time windows are alternated, a time window of 2 times, that is, 60s is used instead of 30s when calculating the average value. And the present embodiment appropriately increases the standard deviation threshold value delta 1 The method can avoid the interruption of the steady state by tiny disturbance, and increase the persistence of the steady state, thereby improving the measurement precision and better identifying the increase and decrease events of the object; for example: delta 1 When the number is 200, the number of mobile phones on the bed can be accurately identified. Of course, the mean value can also be calculated by adopting a time window, and the long steady-state standard deviation threshold delta needs to be reduced 1 And thus reduce the ability to identify increasing or decreasing events of the object.
In a specific combination example, the length of the long time window is set to 30 seconds, and each measurement interval is set to 1 second. The current time is 8 hours and 01 minutes and 00 seconds (i.e. 08:01:00), and a long time window taking the time as the end time is 08:00:31-8:01:00; since the long time window contains 30 measurement instants, there are 30 total weight S values; then calculate the average number of 60 total weights S in a 2-time window (i.e. 08:00:01-8:01:00), and calculate the standard deviation of the values of 30 total weights S in 08:00:31-8:01:00 with the average number of 60 total weights S in 08:00:01-8:01:00 as the average value, if not greater than 200, it is indicated that the long time window of 08:00:31-08:01:00 is in a steady state, and the moment of 8.01 minutes 00 seconds is in a steady state (i.e. the moment in steady state). Otherwise, the long time window of 08:00:31 to 08:01:00 is described as being in an unsteady state, and the time of 8.01 minutes 00 seconds is not in a steady state (i.e., is not in a steady state).
Each time of measurement, whether the current time of measurement is steady state or not needs to be judged according to the above examples, and the examples are specifically combined:
assume that the first measurement time after getting in bed: 8:01:00, the end time of the long time window is 8:01:00, the forward pushing is carried out for 30 seconds (the length of the long time window) and the average of the previous 60 seconds is taken as the average value to judge whether 8:00:00 is steady state or not, namely, whether the total weight standard deviation of 30 times obtained in the 30 seconds exceeds 200 or not is judged from 08:00:31-8:01:00, if yes, the 8:01:00 is in unsteady state.
Second measurement: 8:01:01, the end time of the long time window is 8:01:01, the time is pushed forward for 30 seconds to judge whether 8:01:01 is in a steady state or not, namely, whether the total weight standard deviation of 30 times obtained in the 30 seconds is more than 200 is judged from 08:00:32-8:01:01, and if the total weight standard deviation is not more than 200, the time is 8:01:01 in the steady state.
Third measurement: 8:01:02, the end time of the long time window is 8:01:02, the time is pushed forward for 30 seconds to judge whether 8:01:02 is in a steady state or not, namely, whether the total weight standard deviation of 30 times obtained in the 30 seconds is more than 200 is judged from 08:00:33-8:01:02, and if the total weight standard deviation is not more than 200, the time is 8:01:02 in a steady state.
Fourth measurement: 8:01:03, and pushing forward for 30 seconds to determine whether 8:01:03 is in a steady state or not, namely, whether the total 30 total weight standard deviation obtained from 08:00:34-8:01:03 in 30 seconds exceeds 200, if not, the total weight standard deviation is in the steady state 8:01:03.
The fifth measurement, 8:01:04, was at steady state.
Sixth measurement, 8:01:05, is at steady state.
… … (in steady state)
8:31:59, in steady state.
8:32:00, in unsteady state.
… … (in unsteady state and no occurrence of a bedevent is detected).
8:35:00, in an unsteady state, and no occurrence of a bedevent is detected.
8:35:01, in steady state.
… … (in steady state)
11:00:00, in steady state.
11:00:01, in unsteady state.
… … (in unsteady state)
11:01:00, and detecting that an out-of-bed event occurs.
… … (in empty bed)
11:03:00, and detecting that a loading event occurs.
… … (in unsteady state)
11:03:10, in unsteady state.
11:03:11, in steady state.
… … (in steady state)
13:09:00, in steady state.
13:09:01, in unsteady state.
… … (in unsteady state)
13:10:00, detecting that an out-of-bed event occurs.
… … (in empty bed)
13:40:00, and detecting that no event of getting in bed occurs.
Steady state is defined as: if a steady state long time window is included between two adjacent non-steady state long time windows, the duration of the steady state long time window between the two adjacent non-steady state long time windows is defined as a steady state.
In the present embodiment, the time of the ith measurement is defined as t i I is a natural number; defining the starting and ending time as t i Is TW i Defining the ending time as t i+1 Is TW i+1 If the window TW is long i Is in an unstable state, and has a long time window TW i+1 Is in a steady state, then it is considered to be from the long time window TW i+1 End time t of (2) i+1 Start to enter steady state, let t i+1 Defining as the starting time of a steady state; if long time window TW i ～TW i+k Are all in a stable state, k is a natural number, and the long time window TW i+k+1 Is not stable, then it is considered to be in the long time window TW i+k End time t of (2) i+k Ending steady state, let t i+k Defined as the end time of the steady state.
For example: in the above example, the long time window with the end time 08:01:00 is in an unstable state, and the long time window with the end time 08:01:01 is in a stable state, so the end time 08:01:01 is the start time of the first steady state after getting on bed. Since the long time windows with the end time of 08:01:01 to 08:31:59 are all in a stable state and the long time window with the end time of 08:32:00 is in an unstable state, the time period of the end time of 08:01:01 to 08:31:59 is considered to be a steady state, and the end time of the steady state is considered to be the end time of the moment 08:31:59. The long time window with the end time of 8:32:00 to 8:35:00 is in an unstable state, and the long time window with the end time of 8:35:01 is in a stable state again, so that the end time of 8:35:01 is the start time of the next steady state.
As shown in fig. 2, in this embodiment, the step S2 is split into the following sub-steps:
step S201, judging whether a loading event occurs according to the measured value of the pressure sensor; executing step S202 if a get-in event occurs; in this embodiment, the instant weight W is preferably used to determine the event of getting on and off, but it is also possible to determine the event directly from the total weight S.
To eliminate errors in the sensors due to interference factors such as temperature changes during empty bed conditions, the pressure sensors may be calibrated at intervals during steady state of the empty bed conditions, for example, at intervals of 30 minutes. The calibration method is to set the value of the instantaneous weight W at the current measurement time to 0 by taking the value of the total weight S currently measured (i.e., the average or median of the total weight at each measurement time in a short time window having the current measurement time as the end time) as the value of the empty bed weight B. The influence of temperature change on the pressure sensor can be reduced through calibration, the measurement accuracy is improved, and the interference generated after adding or reducing articles on the bed can be eliminated.
The method for judging whether the event of getting in or out of bed occurs comprises the following steps: if a certain measuring moment of the sensor is changed from an empty bed state to a bed state, judging that a loading event occurs at the measuring moment; if a certain measuring moment of the sensor is changed from the in-bed state to the empty-bed state, the occurrence of an out-of-bed event at the measuring moment is judged.
The following method can be adopted for judging whether the bed is empty: define a length of vT 0 Wherein V is a natural number (V is generally not greater than 10, e.g., V may be set to 3, 5, etc.),judging whether a short time window taking the current measurement time as the end time is in a stable state or not, if the short time window is in a stable state, comparing the average value or the median of the weight W at each moment in the short time window with an empty bed threshold value, and if the average value or the median is smaller than the empty bed threshold value, judging that the measurement time is in an empty bed state.
The method for judging whether the short time window is in a stable state comprises the following steps: defining standard deviation of all instantaneous weight W values recorded in a short time window as sigma TWD Setting a short steady-state standard deviation threshold delta 0 When a short time window ends, if sigma of the short time window TWD ≤δ 0 Judging that the short time window is in a stable state; if sigma TWD ＞δ 0 It is determined that the short time window is in an unstable state. The method for judging whether the short time window is in a stable state is similar to the method for judging the long time window, except that the duration and standard deviation threshold of the time window are different.
The following method can be adopted for judging whether the bed is in a state: setting a reference weight W of the user during the measurement period r The method comprises the steps of carrying out a first treatment on the surface of the And sets a weight difference threshold value, and sets (W r -W) is compared with a weight difference threshold value if (W) r -W) is less than or equal to the weight difference threshold value. At the time of initial use, the weight of the user should be set as the reference weight W r The reference weight W can be compared with the measured value in the measuring process r Is adaptively adjusted.
Step S202, detecting whether a steady state is entered and whether a getting-out event occurs, if so, indicating that the user enters a sleep state, and executing step S203; if a get-off event occurs, it is indicated that the user is not sleeping, and the process returns to step S201. For example: in the above example, after getting on bed, the long time window with the end time 08:01:00 is in an unstable state, and the long time window with the end time 08:01:01 is in a stable state, so the start time 08:01:01 is the first stable state after getting on bed, and is also the sleep start time. Step S202 is preferably divided into step S2021 and step S2022:
step S2021, judging whether a launch event occurs, if so, returning to execute step S201; otherwise, step S2022 is performed.
Step S2022, detecting whether a steady state is entered, if so, executing step S203; otherwise, step S2021 is executed back.
Of course, step S2021 may be performed first or step S2022 may be performed first in the implementation.
Step S203, whether the steady state is ended is detected, if the steady state is ended, step S204 is executed, and statistics and calculation are carried out on the measured value during the steady state. Continuing with the description of the above example: and the measurement result of the time 08:01:01 is a long time window for measuring the stable state for the first time after the patient gets in bed, and the patient enters a sleep state. Since the long time windows with the end time of 08:01:01 to 08:31:59 are all in a stable state and the long time window with the end time of 08:32:00 is in an unstable state, the time period of the end time of 08:01:01 to 08:31:59 is considered to be a steady state, and the end time of the steady state is the end time of the time 08:31:59.
Step S204, calculating duration, slope K and weight reduction ΔW of the steady state a . Of course, the weight loss value ΔW may not be calculated in this step a After screening the effective sleep steady state in step S3, only calculating the steady state weight reduction value DeltaW of the effective sleep a 。
In this embodiment, the least square method is preferably adopted to linearly fit the values of all the instant weight W in the steady state and calculate the slope; the calculation formula is as follows:
wherein N represents the number of values of the instantaneous body weight involved in the steady state, x i A time value corresponding to the value of the body weight at each moment in the steady state;representing an average of time values in steady state; y is i Transient body representing steady state at each momentA heavy value; />Mean of values representing instantaneous body weight in steady state.
Of course, other linear fitting methods may be used to calculate the slope of the steady state, as long as the slope can be calculated from the values of all instantaneous body weights W in the steady state.
The method for calculating the weight loss value during the steady state is as follows: the weight difference was calculated from the duration and slope of the steady state, namely:
ΔW a ＝-K×ΔT W
wherein DeltaT W Indicating the duration of steady state; this eliminates the disturbance of the fluctuation at the start or end of the steady state, and the calculation result is more reliable.
Of course, the weight reduction ΔW during steady state can also be calculated directly from the steady state weight a The calculation formula is as follows:
ΔW a ＝W a+ -W a-
wherein W is a+ A steady state body weight representing the onset of effective sleep steady state; w (W) a- Indicating steady state body weight at the end of effective sleep steady state. Since the value of the empty bed weight B is not recalculated during the steady state duration, the empty bed weight B value must be equal during the same steady state, and thus the steady state reading can also be directly used to calculate the weight reduction value aw during steady state a ：
ΔW a ＝S aw+ -S aw-
Wherein S is aw+ A steady state reading representing an effective sleep steady state onset time; s is S aw- A steady state reading representing the end of active sleep steady state.
As shown in fig. 3 and 4, when the instantaneous weight is taken as the weighing reading, the change of the weighing reading for 8 hours is shown schematically, and it can be seen from fig. 4 that the weighing reading curve looks flat in the ordinate range of 0-100 kg, and has very large jump all the time after amplification, so that the reading has huge jump in unsteady state, and even in steady state, no stable reading exists, and the instantaneous weight W at a single moment cannot be directly used for calculation. Thus, two things are done before computation: 1. removing readings at non-steady state; 2. and (5) calculating the average value of the readings at the steady state.
Newly defining the value of the total weight S during the unsteady state as a null value; definition of steady state reading S a Calculating the current measurement time as the end time and the length rT 0 The average of all total weights S over a period of time of (2) as steady state reading S a Where r is a natural number and r.gtoreq.u, preferably r=2u, i.e. the steady state reading S is taken as the one minute average of the total weight S a . It should be noted that a null value is no value, and is not "0", because "0" is a valued value, and participates in the calculation at the time of averaging; for example, for the first moment of steady state, since the previous 59 moments are unsteady, the total weight is null, so there is only one total weight value, and the average S at that moment is calculated TW 59 null values are removed, and only the average value of one value is calculated; if the total weight during the unsteady state is defined as "0", the average of the total weight at that time and 59 "0" s is calculated, so that an error may occur. Definition of steady state body weight W a ＝S a -B; a steady-state body weight (one minute mean) is output every 30s as a weighing reading, and a change of the weighing reading at this time is shown in fig. 5, at this time, we consider the part outputting the null value as unsteady and equating it to an abnormal event, consider the part outputting the non-null value as steady-state and equating it to effective sleep, and the disturbance of the abnormal event can be effectively eliminated by calculating only the body weight reduction at steady-state.
Step S205, whether a steady state is re-entered or whether a getting-out event occurs is detected, and if the steady state is re-entered, the step S204 is executed again; step S206 is performed if a get-off event occurs.
Continuing with the description of the above example: after the first steady state of time 08:31:59 ends, the long time window with time 8:32:00 to 8:35:00 as the end time is in an unsteady state, and the long time window with time 8:35:01 as the end time is in a steady state again, so time 8:35:01 is the starting time of the next steady state. Since the first steady state is ended and the getting-out event does not occur, the sleep is judged not to be ended, and the step S206 is not performed to judge the sleep end and the getting-out event. After the steady state is completed at time 11:00:00, the occurrence of the getting-out event is detected at 11:01:00, so that step S206 is executed to determine the sleep completion and the getting-in/out event. Similarly, after the steady state at the time 13:09:00 is ended, the occurrence of the getting-out event is detected at the time 13:10:00, and step S206 is also executed to determine that the sleep is ended and the getting-in/out event is performed.
Step S206, detecting whether a loading event occurs within the threshold of the loading and unloading time and entering a steady state, and if the loading event occurs within the threshold of the loading and unloading time and entering the steady state, returning to execute step S204; otherwise, step S207 is performed.
Continuing with the description of the above example: after detecting the occurrence of a bedevent at 11:01:00, detecting the occurrence of a bedevent at 11:03:00, and again being in steady state at 11:03:11. Since the time interval from the time of getting out of bed (11:01:00) to the time of getting in the steady state again (11:03:11) is smaller than the threshold value of getting in and out of bed (30 minutes), it is judged that the event of getting in and out of bed occurs, and sleep is not ended. After detecting the occurrence of the getting-out event at the time 13:10:00, until the getting-in event is not detected at the time 11:40:00, the time interval exceeds the getting-in and getting-out threshold value, so that the sleep is judged to be ended; time 13:09:00 is the end time of steady state before the occurrence of the beddown event, and is also taken as the end time of sleep.
Step S207, determining that sleep is completed, and executing step S3.
And S3, screening out effective sleep stable states according to the slope and duration time of each stable state in the sleep period, and calculating the sum of weight reduction values in each effective sleep stable state period as an absolute weight reduction value D in the sleep period.
The method for screening the effective sleep stable state comprises the following steps:
defining a slope threshold and a steady state duration threshold, firstly calculating the average value of slopes of all steady states during sleep, and then removing the steady state with the absolute value of the difference between the slope and the average value being larger than the slope threshold and the steady state with the positive slope; and removing the stable state with the duration less than the threshold value of the stable state duration, and taking the rest stable state as the effective sleep stable state. The slope threshold may be adjusted according to the actual measurement situation, so as to screen out a steady state with obvious abnormality, and the steady state duration threshold is generally selected to be 5min to 15min, and in this embodiment, preferably 10min.
As shown in table 1, 23 steady states were identified for the slope and duration statistics of the steady states in fig. 5.
TABLE 1
| Steady state sequence number | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
| Slope of | -1.01143 | 0.00414 | -0.01466 | -0.00938 | -0.01547 | -0.02469 | -0.12304 | -0.01671 |
| Time hours | 0.04167 | 0.31667 | 0.90000 | 0.27500 | 0.10833 | 0.05833 | 0.23333 | 0.08333 |
| Steady state sequence number | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
| Slope of | -0.02102 | 0.04059 | -0.02770 | 0.03900 | -0.04728 | -0.04366 | -0.01842 | -0.07339 |
| Time hours | 1.0000 | 0.11667 | 0.48333 | 0.06667 | 0.25000 | 0.64167 | 0.10000 | 0.15000 |
| Steady state sequence number | 17 | 18 | 19 | 20 | 21 | 22 | 23 | |
| Slope of | -0.01483 | -0.00897 | -0.05706 | -0.14544 | -0.03605 | -1.12640 | -0.39966 | |
| Time hours | 0.09166 | 0.04167 | 0.12500 | 0.15833 | 0.47500 | 0.02500 | 0.12500 |
Calculating the average value of slopes of 23 stable states in the table, and then removing the stable states with larger slope deviation average value, positive slope and stable state duration less than 10min, and screening out the stable states 3, 4, 7, 9, 11, 13, 14 and 21 as effective sleep stable states. The slope threshold is generally set to 0.1-0.2, and in this embodiment, the slope threshold is 0.15.
The formula for calculating the sum of the weight loss values during each effective sleep homeostasis period as the absolute weight loss value D during sleep is as follows:
D＝∑ΔW
s4, calculating a sleep metabolic rate by using the absolute weight reduction value D in the sleep period; preferably, the sleep metabolic rate m and the 8-hour daily sleep metabolic rate m (8 h) are defined to calculate the sleep metabolic rate; of course, other indicators (e.g., 12 hour daily sleep metabolic rate, etc.) may be used to calculate the sleep metabolic rate.
The calculation formula of the sleep metabolic rate m is as follows:
wherein W is s1 Indicating the body weight at the beginning of sleep.
The calculation formula of the 8-hour daily sleep metabolic rate m (8 h) is as follows:
where Δt represents the duration of sleep, and in this embodiment represents the total duration of each steady state during sleep. Sleep metabolic rate m and 8 hours daily sleep metabolic rate m (8 hours) are preferably expressed in thousands of points; this index is used to normalize the sleep metabolic rate at different sleep times.
As shown in table 2, the effective sleep homeostasis selected from table 1, and the calculated absolute weight loss value D, sleep metabolic rate m, and 8-hour daily sleep metabolic rate m (8 h) during sleep were calculated.
TABLE 2
In this embodiment, by measuring the absolute decrease in weight of a person during a period of sleeping in a bed, the sleep metabolic rate can be calculated from the decrease in weight during sleep, and thus the metabolic level of the person can be reflected by long-term statistics, as with the caloric metabolic rate. Since the measurement of body weight is more direct than the measurement of calories, the metabolic situation of a person can be reflected more realistically. By introducing the steady state and calculating the weight reduction value of each steady state duration, various interference factors (such as toilet, water drinking, clothes adding and the like) affecting the weight in the sleeping process can be eliminated, and the weight measurement is more accurate. By setting a larger standard deviation threshold delta 1 The method can avoid the interruption of the steady state by tiny disturbance, and increase the persistence of the steady state, thereby improving the measurement precision and better identifying the increase and decrease events of the object. By setting the slope threshold to remove the steady state with larger deviation of slope, the steady state containing interference can be removed, and the increase of the standard deviation threshold delta can be eliminated 1 Resulting in unrecognized interference.
Finally, it is noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the spirit and scope of the present invention, which is intended to be covered by the claims of the present invention.
Claims (6)
1. A method for calculating a metabolic rate of a sleep process based on slope disturbance removal, comprising the steps of:
step S1, a first preset time T is set every interval 0 Measuring the weight of the primary bed; the method comprises the following steps:
a plurality of pressure sensors are arranged under the bed, and each pressure sensor is arranged at intervals of a first preset time T 0 Measuring the weight of the primary bed; defining total weight S as the sum of the measurements of the pressure sensors, defining empty bed weight B as the sum of the measurements of the pressure sensors when the bed is empty, defining instantaneous weight w=s-B;
step S2, defining the length as uT 0 The time period of the time frame is a long time window, whether the long time window taking the current measurement time as the end time is in a stable state is judged according to the measured value of the weight in the long time window, and if the long time window is in the stable state, the current state is judged; recognizing the starting time and the ending time of sleep according to the stable states, and calculating the duration time and the slope of each stable state in the sleep period; the method for identifying the starting time of sleeping comprises the following steps:
judging whether a loading event occurs according to the measured total weight S; and judging whether to enter a steady state after the occurrence of the loading event, if so, judging to enter a sleep state, and taking the starting time of the steady state as the starting time of sleep;
step S3, screening effective sleep stable states according to the slope and duration time of each stable state in the sleep period, and calculating the sum of weight reduction values in each effective sleep stable state period as an absolute weight reduction value D in the sleep period;
the calculating method of the slope of the steady state comprises the steps of adopting a least square method to linearly fit the values of all the instant weight W in the steady state and calculating the slope;
the method for screening the effective sleep stable state comprises the following steps:
defining a slope threshold and a steady state duration threshold, firstly calculating the average value of slopes of all steady states during sleep, and then removing the steady state with the absolute value of the difference between the slope and the average value being larger than the slope threshold and the steady state with the positive slope; removing the stable state with the duration less than the threshold value of the stable state duration, and taking the rest stable state as an effective sleep stable state;
the weight loss value DeltaW during steady state was calculated using the following a ：
ΔW a ＝-K×ΔT W
Wherein DeltaT W Indicating the duration of steady state;
and S4, calculating the sleep metabolism rate by using the value of the absolute weight reduction during the sleep period.
2. The method for calculating the metabolic rate of a sleep process based on slope disturbance elimination according to claim 1, wherein said step S2 includes the sub-steps of:
step S201, judging whether a loading event occurs according to the measured value of the pressure sensor; executing step S202 if a get-in event occurs;
step S202, detecting whether a steady state is entered and whether a getting-out event occurs, and judging that the sleeping state is entered if the steady state is entered; step S203 is performed; returning to execute step S201 if a get-off event occurs;
step S203, detecting whether the steady state is ended, if so, executing step S204;
step S204, calculating the duration and the slope of the steady state;
step S205, whether a steady state is re-entered or whether a getting-out event occurs is detected, and if the steady state is re-entered, the step S204 is executed again; executing step S206 if a get-off event occurs;
step S206, detecting whether a loading event occurs within the threshold of the loading and unloading time and entering a steady state, and if the loading event occurs within the threshold of the loading and unloading time and entering the steady state, returning to execute step S204; otherwise, step S207 is performed;
step S207, determining that sleep is completed, and executing step S3.
3. The method for calculating the metabolic rate of a sleep process based on slope interference elimination according to claim 1, wherein the method for judging whether the long time window is in a steady state is as follows: calculating the current measurement time as the end time and the length of 2uT 0 Average S of all total weights S over a period of time TW By means of average number S TW Standard for calculating the value of all total weights S recorded in the long time window for the total weight meanThe standard deviation is defined as sigma TWC Setting a long steady-state standard deviation threshold delta 1 When a long time window ends, if sigma of the long time window TWC ≤δ 1 Judging that the long time window is in a stable state; if sigma TWC ＞δ 1 The long time window is determined to be in an unstable state.
4. The method for calculating the metabolic rate of a sleep process based on slope disturbance removal according to claim 2, wherein the method for determining whether a get-in event and a get-out event occur is: if a certain measuring moment of the sensor is changed from an empty bed state to a bed state, judging that a loading event occurs at the measuring moment; if a certain measuring moment of the sensor is changed from the in-bed state to the empty-bed state, the occurrence of an out-of-bed event at the measuring moment is judged.
5. The method for calculating a metabolic rate of a sleep process based on slope interference according to claim 1, wherein the method for calculating a sleep metabolic rate using an absolute weight reduction value D during sleep comprises:
the sleep metabolic rate m is calculated according to the following formula:
wherein W is s1 Indicating the body weight at the beginning of sleep.
6. The method for calculating a metabolic rate for a sleep process based on slope interference as set forth in claim 5, wherein the method for calculating a sleep metabolic rate using an absolute weight reduction value during sleep further comprises:
the 8-hour daily sleep metabolic rate m (8 h) was calculated according to the following:
where Δt represents the duration of sleep.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111322781.2A CN113974568B (en) | 2021-11-09 | 2021-11-09 | Slope interference-free method for calculating metabolic rate of sleep process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111322781.2A CN113974568B (en) | 2021-11-09 | 2021-11-09 | Slope interference-free method for calculating metabolic rate of sleep process |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113974568A CN113974568A (en) | 2022-01-28 |
| CN113974568B true CN113974568B (en) | 2024-03-26 |
Family
ID=79747454
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111322781.2A Active CN113974568B (en) | 2021-11-09 | 2021-11-09 | Slope interference-free method for calculating metabolic rate of sleep process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113974568B (en) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006106516A2 (en) * | 2005-04-05 | 2006-10-12 | Andante Medical Devices Ltd. | Rehabilitation system |
| CN1883527A (en) * | 2006-06-06 | 2006-12-27 | 赵一弘 | A health product of lucid ganoderma and preparation method thereof |
| JP3967365B1 (en) * | 2006-12-12 | 2007-08-29 | 株式会社クロスウェル | Biological information management system, control method thereof, and program |
| CN101568925A (en) * | 2006-12-22 | 2009-10-28 | Icw公司 | Automatic priorization of instances of medical care |
| JP2013198574A (en) * | 2012-03-23 | 2013-10-03 | Tanita Corp | Biometric apparatus and biometric system |
| WO2013157307A1 (en) * | 2012-04-20 | 2013-10-24 | Emura Tokuo | Pedometer |
| CN103748465A (en) * | 2011-08-05 | 2014-04-23 | 美艾利尔圣地亚哥公司 | Methods and compositions for monitoring heart failure |
| CN107153764A (en) * | 2017-05-08 | 2017-09-12 | 杭州变啦网络科技有限公司 | A kind of data correcting method of intelligent human-body composition scale |
| CN107368768A (en) * | 2016-05-11 | 2017-11-21 | 欧姆龙健康医疗事业株式会社 | Suitable for the user identification method and scale of scale |
| CN107669687A (en) * | 2010-09-01 | 2018-02-09 | 艾尼纳制药公司 | 5 HT for body weight control2CThe modified release dosage form of activator |
| CN108430314A (en) * | 2015-12-18 | 2018-08-21 | 苏州大学 | Biological rhythm adjusts system and biological rhythm method of adjustment |
| CN108523899A (en) * | 2018-04-23 | 2018-09-14 | 深圳市友宏科技有限公司 | Monitoring method, storage medium in sleep procedure from bed state and device |
| CN112274901A (en) * | 2019-07-22 | 2021-01-29 | 财团法人工业技术研究院 | Weight training method and system |
| CN112367909A (en) * | 2018-07-02 | 2021-02-12 | 三星电子株式会社 | Electronic device and method for providing personalized biometric information based on biometric signals using the same |
| CN112914884A (en) * | 2021-01-19 | 2021-06-08 | 重庆火后草科技有限公司 | Method for measuring weight value in sleep state through steady-state duration confidence |
| CN112924007A (en) * | 2021-01-19 | 2021-06-08 | 重庆火后草科技有限公司 | Weight measurement method based on target sleep |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE60335899D1 (en) * | 2002-09-19 | 2011-03-10 | Univ Ramot | METHOD, DEVICE AND SYSTEM FOR CHARACTERIZING SLEEP |
| US8202224B2 (en) * | 2006-11-13 | 2012-06-19 | Pacesetter, Inc. | System and method for calibrating cardiac pressure measurements derived from signals detected by an implantable medical device |
| US20090076347A1 (en) * | 2007-09-17 | 2009-03-19 | Shape Medical Systems, Inc. | Pattern Recognition System for Classifying the Functional Status of Patients with Chronic Disease |
| EP2489307B1 (en) * | 2009-10-14 | 2018-10-03 | Delta Tooling Co., Ltd. | Biological state estimation device and computer program |
| JP5440434B2 (en) * | 2010-07-22 | 2014-03-12 | オムロンヘルスケア株式会社 | Weight management device |
| CN107106029A (en) * | 2014-12-18 | 2017-08-29 | 皇家飞利浦有限公司 | System from Health ＆ Fitness Tip to person under inspection, equipment, method and computer program for providing |
-
2021
- 2021-11-09 CN CN202111322781.2A patent/CN113974568B/en active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006106516A2 (en) * | 2005-04-05 | 2006-10-12 | Andante Medical Devices Ltd. | Rehabilitation system |
| CN1883527A (en) * | 2006-06-06 | 2006-12-27 | 赵一弘 | A health product of lucid ganoderma and preparation method thereof |
| JP3967365B1 (en) * | 2006-12-12 | 2007-08-29 | 株式会社クロスウェル | Biological information management system, control method thereof, and program |
| CN101568925A (en) * | 2006-12-22 | 2009-10-28 | Icw公司 | Automatic priorization of instances of medical care |
| CN107669687A (en) * | 2010-09-01 | 2018-02-09 | 艾尼纳制药公司 | 5 HT for body weight control2CThe modified release dosage form of activator |
| CN103748465A (en) * | 2011-08-05 | 2014-04-23 | 美艾利尔圣地亚哥公司 | Methods and compositions for monitoring heart failure |
| JP2013198574A (en) * | 2012-03-23 | 2013-10-03 | Tanita Corp | Biometric apparatus and biometric system |
| WO2013157307A1 (en) * | 2012-04-20 | 2013-10-24 | Emura Tokuo | Pedometer |
| CN108430314A (en) * | 2015-12-18 | 2018-08-21 | 苏州大学 | Biological rhythm adjusts system and biological rhythm method of adjustment |
| CN107368768A (en) * | 2016-05-11 | 2017-11-21 | 欧姆龙健康医疗事业株式会社 | Suitable for the user identification method and scale of scale |
| CN107153764A (en) * | 2017-05-08 | 2017-09-12 | 杭州变啦网络科技有限公司 | A kind of data correcting method of intelligent human-body composition scale |
| CN108523899A (en) * | 2018-04-23 | 2018-09-14 | 深圳市友宏科技有限公司 | Monitoring method, storage medium in sleep procedure from bed state and device |
| CN112367909A (en) * | 2018-07-02 | 2021-02-12 | 三星电子株式会社 | Electronic device and method for providing personalized biometric information based on biometric signals using the same |
| CN112274901A (en) * | 2019-07-22 | 2021-01-29 | 财团法人工业技术研究院 | Weight training method and system |
| CN112914884A (en) * | 2021-01-19 | 2021-06-08 | 重庆火后草科技有限公司 | Method for measuring weight value in sleep state through steady-state duration confidence |
| CN112924007A (en) * | 2021-01-19 | 2021-06-08 | 重庆火后草科技有限公司 | Weight measurement method based on target sleep |
Non-Patent Citations (2)
| Title |
|---|
| Zhang, K ; Sun, M ; Werner, P ; 等.Sleeping metabolic rate in relation to body mass index and body composition.《INTERNATIONAL JOURNAL OF OBESITY》.2002,全文. * |
| 睡眠过程心率变异性分析及睡眠呼吸暂停综合征的判别研究;江丽仪;《中国优秀硕士学位论文全文数据库医药卫生科技辑》;20111215;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113974568A (en) | 2022-01-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104169680B (en) | Use the Level Change of baroceptor monitoring device | |
| US10448928B2 (en) | Method and device for detecting physiological index | |
| CN112924007B (en) | Weight measurement method based on target sleep | |
| CN114027792B (en) | Metabolic rate detection method for sleep process based on linear correlation coefficient interference elimination | |
| CN112284556B (en) | Calibration system and method with temperature compensation | |
| CN110849459B (en) | Creep deformation correction method for weighing sensor | |
| CN104921736B (en) | A kind of continuous blood sugar monitoring device for including parameter Estimation function filtration module | |
| US20070238938A1 (en) | Activity-induced energy expenditure estimating instrument | |
| US20070168230A1 (en) | Method for conditioning body weight and device for implementing this method | |
| JP2012026846A5 (en) | ||
| CN105258779A (en) | Automatic calibration method of electronic scale | |
| CN103901942B (en) | For the calibration method and device of the clock accuracy of terminal | |
| CN113974568B (en) | Slope interference-free method for calculating metabolic rate of sleep process | |
| CN107992991A (en) | Annual electricity sales amount Forecasting Methodology based on external environmental factor and Co-integration Theory | |
| CN112891098B (en) | Body weight measuring method for health monitor | |
| CN103961066A (en) | Temperature measurement method and temperature measurement device | |
| CN112462846A (en) | Clock correction method, clock correction device, computer equipment and storage medium | |
| CN112914884B (en) | Method for measuring weight value in sleep state through steady-state duration confidence | |
| CN111238624A (en) | Object weight measuring method, weighing device and readable storage medium | |
| US10062453B1 (en) | Calibrating I/O impedances using estimation of memory die temperature | |
| CN112914883B (en) | Method for measuring weight value in sleep state through eccentricity confidence | |
| CN103344742A (en) | Method for predicting criticality of entomophthora in grains stored in granary | |
| CN113974567B (en) | Method for calculating metabolic rate in sleeping process | |
| CN107576710B (en) | Electrochemical gas sensor and its calibration method, air conditioner | |
| CN112284469B (en) | Zero drift processing method of ultrasonic water meter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |