A sleep quality analysis system and method based on a medical bed sensor

Abstract

This invention belongs to the technical field of sleep quality analysis, specifically relating to a sleep quality analysis system and method based on sensors in a medical bed, to address the technical problem of existing methods relying too heavily on duration statistics while neglecting the continuity and circadian rhythm of sleep. The analysis method includes: S1, acquiring vital sign signal data collected by sensors in the medical bed within a sleep cycle and dividing it into multiple continuous time segments; S2, extracting physiological indicators such as respiratory rate, heart rate, and body movement amplitude within each time segment, and calculating the variability depth of the corresponding time segment; S3, generating a variability depth baseline; calculating non-steady-state disturbance factors based on the frequency and duration of micro-awakening disturbance segments; and S4, generating a comprehensive sleep quality score. The analysis method of this invention improves the comprehensiveness and objectivity of the assessment results, more realistically reflecting sleep health status.

Description

Technical Field

[0001] This invention belongs to the technical field of sleep quality analysis, specifically relating to a sleep quality analysis system and method based on sensors in a medical bed. Background Technology

[0002] Traditional sleep quality assessment primarily relies on the gold standard—polysomnography (PSG). This technology collects various physiological signals, such as electroencephalograms (EEG), electrooculograms (EOG), electromyograms (EMG), and electrocardiograms (ECG), by attaching electrodes to various parts of the subject's body. It can divide sleep stages and diagnose various sleep disorders. However, PSG equipment is expensive, complex to operate, and the monitoring process is highly invasive. Subjects need to undergo the process in a professional sleep laboratory environment, and the cables connected throughout the body can severely interfere with normal sleep, resulting in monitoring results that may not fully reflect the natural sleep state.

[0003] In contrast, sleep monitoring solutions based on non-contact sensors, such as pressure, vibration, or piezoelectric sensors installed under medical beds or in mattresses, extract basic physiological indicators like respiratory rate, heart rate, and body movement by sensing weak vibration signals caused by changes in breathing, heart rate, and body posture during sleep. However, generating an assessment report solely by statistically analyzing the duration of each sleep stage, while intuitive, is rather crude and neglects two crucial dimensions of the sleep process: first, the continuity and stability of sleep, avoiding the identification of transient disturbances such as micro-awakenings and body movements that disrupt sleep structure; and second, the rhythmicity of the sleep cycle, failing to represent the degree of deviation between the actual sleep cycle structure and the standard healthy sleep cycle template. Therefore, existing sensor-based assessment models for medical beds are insufficient in terms of comprehensiveness and accuracy, and the scoring results often fail to fully reflect true sleep quality. Summary of the Invention

[0004] This invention provides a sleep quality analysis system and method based on sensors in a medical bed, to solve the technical problem that existing methods rely too much on duration statistics and ignore the continuity and circadian rhythm of sleep.

[0005] In a first aspect, the present invention provides a sleep quality analysis method based on sensors in a medical bed, comprising the following steps:

[0006] S1, acquire vital sign signal data collected by sensors on the medical bed during a sleep cycle, and divide it into multiple consecutive time segments;

[0007] S2, extract the physiological indicators of respiratory rate, heart rate and body movement amplitude in each time segment, and preliminarily determine the sleep stage of the corresponding time segment based on the preset physiological indicator correlation matrix of sleep stage; at the same time, calculate the normalized power spectral entropy of the physiological indicators in each time segment as the variation depth of the corresponding time segment.

[0008] S3. Based on the mean variation depth of the same sleep stage in historical data and the preset stage sensitivity coefficient, a variation depth baseline is generated; when the variation depth of a certain time segment exceeds the variation depth baseline of the corresponding sleep stage, the time segment is marked as a micro-awakening disturbance segment; the non-steady-state disturbance factor is calculated according to the frequency and duration of the micro-awakening disturbance segment.

[0009] S4. The determined sleep stage sequence is matched with the standard sleep cycle template by time regularization, and the cycle rhythm deviation is calculated. The distribution duration of deep sleep, light sleep and REM sleep stages in the whole sleep cycle is statistically analyzed to generate the basic evaluation coefficient. The basic evaluation coefficient, cycle rhythm deviation and non-steady-state disturbance factor are fused and non-linear weighted calculation is performed to generate a comprehensive sleep quality score.

[0010] Furthermore, in S1, the vital sign data is divided into non-overlapping window functions with a time segment length of 30 seconds to generate a series of continuous time segment sequences.

[0011] Furthermore, physiological indicators of respiratory rate, heart rate, and body movement amplitude were extracted for each time segment, including:

[0012] Bandpass filtering is applied to the signal in each time segment to separate heartbeat and respiratory signals; heart rate and respiratory rate are obtained by detecting the peak points of the filtered signals and calculating the peak intervals; the standard deviation of the original signal amplitude within a time segment is calculated as the body motion amplitude index for the corresponding time segment.

[0013] Further, the normalized power spectral entropy of the physiological indicators within each time segment is calculated, including:

[0014] Perform a Fast Fourier Transform on the heart rate and respiratory rate change sequence within the time segment to obtain the power spectral density distribution; then, extract the energy values ​​of each frequency component of the power spectrum. Divide by total energy The normalized probability is obtained. This forms a probability distribution; according to the Shannon entropy formula... Calculate the information entropy of the probability distribution, and use the information entropy as the mutation depth, where, This represents each individual frequency component in the power spectral density distribution. It represents the total number of all frequency components in the power spectral density distribution.

[0015] Furthermore, when generating the variation depth baseline, at least 50 historical segments in the historical database that belong to the same sleep stage as the current time segment are queried, and the arithmetic mean of the variation depth of the historical segments is calculated; the arithmetic mean is multiplied by an adjustment factor, which is equal to 1 plus the sensitivity coefficient of the corresponding sleep stage, wherein the sensitivity coefficient of the light sleep stage is 0.2, the sensitivity coefficient of the deep sleep stage is 0.3, and the sensitivity coefficient of the REM sleep stage is 0.25, thereby generating the variation depth baseline of the sleep stage.

[0016] Furthermore, based on the frequency and duration of micro-awakening perturbation segments, the non-steady-state perturbation factor is calculated, including:

[0017] The frequency of micro-awakening disturbances is calculated by dividing the total number of micro-awakening disturbances throughout the entire sleep cycle by the total number of sleep hours. ;

[0018] Calculate the average duration of all micro-awakening perturbation segments. ;

[0019] By analyzing the frequency of occurrence and average duration Perform normalized weighted summation to calculate the unsteady-state perturbation factor. The calculation formula is: ;

[0020] in, and These are the preset normalized baseline values ​​for frequency and duration, respectively. and For the corresponding weight coefficients, and .

[0021] Furthermore, when calculating the circadian rhythm deviation, the sleep stages of "deep sleep", "light sleep", "rapid eye movement" and "wakefulness" are mapped to the values ​​4, 3, 2 and 1 respectively, thus forming a numerical sequence of sleep stages;

[0022] The cumulative normalized path cost is obtained by performing time-warping matching between the numerical sequence of sleep stages and the standard sleep cycle template sequence that has also undergone numerical mapping.

[0023] The cumulative normalized path cost is divided by the length of the matching path and normalized. This value is used as the periodic rhythm deviation.

[0024] Furthermore, when generating the baseline assessment coefficient, the proportion of total sleep time spent in deep sleep, light sleep, and REM sleep stages throughout the entire sleep cycle is statistically analyzed and recorded as follows: , , ;

[0025] The basic evaluation coefficient is obtained by weighting and summing the proportions using preset weights. Basic evaluation coefficient The calculation formula is: ;

[0026] in, , , These are the weighting coefficients for each sleep stage, and their sum is 1.

[0027] Furthermore, the formula for generating a comprehensive sleep quality score is as follows: ;

[0028] in, For comprehensive scoring, Based on the evaluation coefficient, This represents the degree of deviation from the periodic rhythm. It is an unsteady-state perturbation factor; , , These are the weighting coefficients for each item. and It is a non-linear adjustment index.

[0029] Secondly, the present invention provides a sleep quality analysis system based on sensors in a medical bed, including a memory and a processor. The memory stores computer program instructions, and when the computer program instructions are executed by the processor, the above-mentioned sleep quality analysis method based on sensors in a medical bed is implemented.

[0030] The beneficial effects are as follows: This invention calculates the normalized power spectral entropy of physiological indicators of vital signs and uses it as the depth of variation to analyze subtle fluctuations during sleep. Based on the baseline of the depth of variation, it identifies micro-awakening disturbances that disrupt sleep continuity, thereby generating non-steady-state disturbance factors and assessing sleep stability. Simultaneously, it performs time-normalization matching between the user's actual sleep stage sequence and a standard sleep cycle template to calculate the circadian rhythm deviation, thus representing the regularity and rhythmicity of sleep structure. By non-linearly fusing and weighting the basic assessment coefficient reflecting sleep volume, the non-steady-state disturbance factor reflecting stability, and the circadian rhythm deviation representing rhythmicity, the comprehensive sleep quality score derived by this invention not only includes an assessment of sleep volume but also incorporates sleep quality, thereby improving the comprehensiveness and objectivity of the assessment results and more realistically reflecting sleep health status. Attached Figure Description

[0031] Figure 1 This is a flowchart of a sleep quality analysis method based on sensors in a medical bed;

[0032] Figure 2This is a schematic diagram illustrating the extraction of physiological indicators and the determination of stages.

[0033] Figure 3 This is a schematic diagram of signal acquisition and segmentation;

[0034] Figure 4 A diagram illustrating the overall sleep quality score;

[0035] Figure 5 This is a schematic diagram of a sleep quality analysis system based on sensors on a medical bed. Detailed Implementation

[0036] An embodiment of the sleep quality analysis method based on medical bed sensors provided by this invention:

[0037] like Figure 1 As shown, a sleep quality analysis based on sensors in a medical bed includes the following steps:

[0038] S1: Acquire vital sign signal data collected by sensors on the medical bed during a sleep cycle and divide it into multiple consecutive time segments.

[0039] Specifically, by using a piezoelectric film sensor deployed under the medical mattress, pressure change signals caused by breathing, heartbeat and body movement are continuously collected throughout the user's sleep at night, forming raw vital sign time series data. The raw vital sign time series data is then divided into multiple continuous time segments according to a fixed time window (e.g., every 30 seconds), with each segment serving as the basic unit for subsequent analysis.

[0040] In an optional embodiment, in S1, the vital sign data is divided into non-overlapping window functions with a time segment length of 30 seconds to generate a series of continuous time segment sequences.

[0041] Specifically, a complete sleep cycle of data is acquired from the piezoelectric thin-film sensor of the medical bed, for example, a continuous signal stream lasting 8 hours. Assuming a sampling frequency of 50Hz, a total of 1.44 million data points are obtained. A fixed time window length of 30 seconds is set, corresponding to 1500 data points. Figure 3 The entire signal stream was sliced ​​in a non-overlapping manner. The first time slice contained data points from 1 to 1500, the second time slice contained data points from 1501 to 3000, and so on, until the end of the signal stream. The slicing process transformed the continuous 8-hour vital sign data into a sequence of 960 independent time slices, each with a duration of 30 seconds. Each slice is the basic unit for subsequent physiological indicator analysis.

[0042] S2, extract the physiological indicators of respiratory rate, heart rate and body movement amplitude in each time segment, and preliminarily determine the sleep stage of the corresponding time segment based on the preset physiological indicator correlation matrix of sleep stage; at the same time, calculate the normalized power spectral entropy of the physiological indicators in each time segment as the variation depth of the corresponding time segment.

[0043] Specifically, for each 30-second time segment of data, digital signal processing techniques are first employed. Low-frequency respiratory waves are extracted using low-pass filtering, and the respiratory rate is obtained by calculating the peak interval of the respiratory waves. Higher-frequency cardiac impulse signals are separated using band-pass filtering, and the heart rate is obtained using a peak detection algorithm. The amplitude of body movement is represented by calculating the energy or variance of the signal. For example... Figure 2 The extracted respiratory rate, heart rate, and body movement amplitude are used as a feature vector and input into a classification model that has been trained with a large amount of clinical data, namely the preset sleep stage physiological index association matrix (such as decision tree or support vector machine), so as to preliminarily determine the most likely sleep stage of the time segment, such as deep sleep, light sleep, REM sleep, or wakefulness.

[0044] Meanwhile, taking heart rate as an example, a continuous heartbeat interval sequence is extracted within this time segment, and a fast Fourier transform is performed on this sequence to obtain the power spectral density. This power spectral density is then normalized so that the sum of the power of all frequency components is one, forming a probability distribution. The information entropy of this probability distribution is calculated based on the Shannon entropy formula, and the resulting entropy value is the variation depth of this time segment. This value can represent the complexity and uncertainty of heart rate changes.

[0045] In an optional embodiment, physiological parameters such as respiratory rate, heart rate, and body movement amplitude are extracted for each time segment, including:

[0046] Bandpass filtering is applied to the signal in each time segment to separate heartbeat and respiratory signals; heart rate and respiratory rate are obtained by detecting the peak points of the filtered signals and calculating the peak intervals; the standard deviation of the original signal amplitude within a time segment is calculated as the body motion amplitude index for the corresponding time segment.

[0047] The generated 30-second time segment data is then processed for signal separation and index extraction. Specifically, a bandpass filter of 0.8 to 2.5 Hz is used to extract the weak vibration signal caused by heartbeats (i.e., cardiac impact signal); simultaneously, another bandpass filter of 0.1 to 0.5 Hz is used to extract the respiratory signal caused by chest cavity fluctuations. In the filtered cardiac impact signal, a peak detection algorithm is used to identify all heartbeat peaks within 30 seconds. For example, if 35 peaks are detected, the average time interval between peaks (e.g., 0.86 s) is used to calculate a heart rate of approximately 70 beats per minute. The same method is applied to the respiratory signal; if 8 respiratory peaks are detected, the respiratory rate is calculated to be approximately 16 breaths per minute. Returning to the unfiltered original 30-second signal segment, the standard deviation of the amplitude of 1500 data points is calculated. This standard deviation is used as an index of body movement amplitude for that segment. For example, the standard deviation is 0.1 when the user's body is stable, but it may increase to 1.5 if there are larger movements such as turning over.

[0048] In an optional embodiment, calculating the normalized power spectral entropy of the physiological indicator within each time segment includes:

[0049] Perform a Fast Fourier Transform on the heart rate and respiratory rate change sequence within the time segment to obtain the power spectral density distribution; then, extract the energy values ​​of each frequency component of the power spectrum. Divide by total energy The normalized probability is obtained. This forms a probability distribution; according to the Shannon entropy formula... Calculate the information entropy of the probability distribution, and use the information entropy as the mutation depth, where, This represents each individual frequency component in the power spectral density distribution. It represents the total number of all frequency components in the power spectral density distribution.

[0050] Taking a heart rate variability sequence as an example, based on the time intervals of the detected consecutive heartbeat peaks, a heartbeat interval sequence is first generated, for example, sequences of 0.85s, 0.87s, and 0.84s within a 30s segment. A Fast Fourier Transform is then performed on the interval sequence to transform the time domain to the frequency domain, resulting in a power spectral density map. This map shows the energy distribution of heart rate variability at different frequencies. Subsequently, the energy value of each frequency component is calculated. And add up the energy of all components to get the total energy. By using the energy of each component Divide by total energy This gives the proportion of that component in the total energy, i.e., the normalized probability. ,all The sum of these probabilities is 1. Substituting these probability values ​​into the Shannon entropy formula, the resulting information entropy value is the variation depth of that time segment. For example, if the energy is mainly concentrated on a few frequencies, a lower entropy value, such as 1.2, is calculated, indicating simple variability; conversely, if the energy is evenly distributed across multiple frequencies, a higher entropy value, such as 3.5, is obtained, indicating complex variability.

[0051] S3. Based on the mean variation depth of the same sleep stage in historical data and the preset stage sensitivity coefficient, a variation depth baseline is generated. When the variation depth of a certain time segment exceeds the variation depth baseline of the corresponding sleep stage, the time segment is marked as a micro-awakening disturbance segment. The non-steady-state disturbance factor is calculated according to the occurrence frequency and duration of the micro-awakening disturbance segment.

[0052] Specifically, historical sleep data is analyzed in advance to calculate the average level of variation depth under different sleep stages such as deep sleep and light sleep. This average level is then multiplied by a set stage sensitivity coefficient (e.g., 1.2 for deep sleep and 1.5 for light sleep) to generate a specific variation depth threshold, i.e., variation depth baseline, for each sleep stage.

[0053] When analyzing the current sleep, the depth of variation calculated for each time segment is compared with the baseline of the corresponding sleep stage. If the depth of variation of the current segment is greater than the corresponding baseline, the sleep is considered to be disturbed, and the time segment is marked as a micro-wake disturbance segment.

[0054] After analyzing the entire sleep cycle, the total number of all segments marked as micro-awakening disturbances is counted and divided by the total sleep duration to obtain the average number of disturbances per hour. The total duration of the time segments is summed, and a comprehensive non-steady-state disturbance factor is calculated by using a weighting formula, such as adding the normalized average number of disturbances to the proportion of the total disturbance duration. The higher the value of the non-steady-state disturbance factor, the more unstable and fragmented the sleep process is.

[0055] In an optional embodiment, when generating the variation depth baseline, at least 50 historical segments in the historical database that belong to the same sleep stage as the current time segment are queried, and the arithmetic mean of the variation depth of the historical segments is calculated; the arithmetic mean is multiplied by an adjustment factor, which is equal to 1 plus the sensitivity coefficient of the corresponding sleep stage, wherein the sensitivity coefficient of the light sleep stage is 0.2, the sensitivity coefficient of the deep sleep stage is 0.3, and the sensitivity coefficient of the REM sleep stage is 0.25, to generate the variation depth baseline of the sleep stage.

[0056] Specifically, assuming the current time segment is classified as deep sleep with a variation depth value of 2.1, the database storing the user's historical sleep data is accessed to retrieve all past time segments marked as deep sleep. The 60 most recent segments are then selected randomly or chronologically. The variation depth values ​​of these 60 historical segments are read (e.g., 1.9, 2.2, 2.0, etc.), and their arithmetic mean is calculated, assuming a result of 2.05. Since the current stage is deep sleep, the corresponding sensitivity coefficient is 0.3, therefore the adjustment factor is calculated to be 1.3. The variation depth baseline is obtained by multiplying the arithmetic mean by the adjustment factor, resulting in 2.665. This value of 2.665 is the personalized judgment threshold for the current deep sleep segment, used for subsequent comparisons to determine whether the physiological state of the current segment is stable.

[0057] In an optional embodiment, the nonsteady-state perturbation factor is calculated based on the frequency and duration of the micro-awakening perturbation segments, including:

[0058] The frequency of micro-awakening disturbances is calculated by dividing the total number of micro-awakening disturbances throughout the entire sleep cycle by the total number of sleep hours. ;

[0059] Calculate the average duration of all micro-awakening perturbation segments. ;

[0060] By analyzing the frequency of occurrence and average duration Perform normalized weighted summation to calculate the unsteady-state perturbation factor. The calculation formula is: ;

[0061] in, and These are the preset normalized baseline values ​​for frequency and duration, respectively. and For the corresponding weight coefficients, and .

[0062] Specifically, after the entire night's sleep monitoring is completed, all 30-second segments identified as micro-awakening disturbances are summarized. Assuming a 7.5-hour sleep period, a total of 15 micro-awakening disturbance segments were detected. Frequency of occurrence... That is, 2 times per hour. Since each segment is 30 seconds long, and the 15 segments may occur sporadically, their average duration is calculated. If they are 15 independent 30-second segments, the average duration is... That is, 30 seconds. Set the baseline values ​​and weights required for standardization, such as frequency baseline values. 4 times / hour, duration baseline For 60 seconds, weight It is 0.7. The value is 0.3. Substitute this value into the formula to calculate the unsteady-state disturbance factor. The calculated result is 0.5.

[0063] S4. The determined sleep stage sequence is matched with the standard sleep cycle template by time regularization, and the cycle rhythm deviation is calculated. The distribution duration of deep sleep, light sleep and REM sleep stages in the whole sleep cycle is statistically analyzed to generate the basic evaluation coefficient. The basic evaluation coefficient, cycle rhythm deviation and non-steady-state disturbance factor are fused and non-linear weighted calculation is performed to generate a comprehensive sleep quality score.

[0064] Specifically, firstly, the circadian rhythm deviation is calculated: the sleep stage determination results of the whole night (e.g., the sequence of light sleep-deep sleep-light sleep-REM sleep) are taken as an actual sleep rhythm sequence; an ideal sleep cycle template sequence is built in, which is constructed from data of healthy people. This template usually contains several complete 90-minute sleep cycles; a time warping algorithm is used to calculate the minimum matching cost between the actual sequence and the ideal template sequence. After normalization, the cost is the circadian rhythm deviation, which represents the difference between the user's sleep cycle structure and the standard health model.

[0065] Secondly, generate basic evaluation coefficients: traverse all time segments throughout the night, count the segments that are determined to be deep sleep, light sleep, and REM sleep, and multiply each count value by the duration of a single segment (e.g., 30 seconds) to calculate the total duration of each sleep stage; use the ratio of key indicators such as deep sleep duration and REM sleep duration to the total sleep duration as core parameters to form basic evaluation coefficients.

[0066] Finally, a nonlinear weighted calculation is performed: the three dimensions of indicators obtained above—basic evaluation coefficient, circadian rhythm deviation, and non-steady-state disturbance factor—are standardized; the three standardized values ​​are then input into a preset nonlinear fusion model, such as a small neural network or a fuzzy logic system. This model can intelligently weight the indicators according to different combinations. For example, when both the circadian rhythm deviation and non-steady-state disturbance factor are poor, the final score will be lowered, even if the deep sleep duration reaches the target. The model will output a comprehensive sleep quality score from 0 to 100.

[0067] In an optional embodiment, when calculating the circadian rhythm deviation, the sleep stages “deep sleep”, “light sleep”, “REM sleep”, and “wake” are mapped to the values ​​4, 3, 2, and 1, respectively, thereby forming a numerical sequence of sleep stages.

[0068] The cumulative normalized path cost is obtained by performing time-warping matching between the numerical sequence of sleep stages and the standard sleep cycle template sequence that has also undergone numerical mapping.

[0069] The cumulative normalized path cost is divided by the length of the matching path and normalized. This value is used as the periodic rhythm deviation.

[0070] Specifically, the analysis results of sleep stages throughout the night, such as wakefulness, light sleep, deep sleep, and REM sleep, are converted into a numerical sequence of sleep stages. According to the mapping rules, this sequence may be 1, 3, 3, 4, 2, etc. Simultaneously, a template sequence representing the ideal sleep rhythm is built-in; for example, a complete 90-minute cycle might be represented as 1, 3, 3, 4, 4, 4, 3, 2, 2. A time warping algorithm is used to compare the user's actual sleep sequence with the ideal template sequence. This algorithm finds an optimal alignment path that minimizes the sum of numerical differences between corresponding points in the two sequences. This minimum sum of differences is the cumulative warping path cost, assumed to be 250. The algorithm also records the length of the optimal alignment path, assumed to be 300. Normalizing the cost by dividing by the path length yields a value of 0.83, which is the circadian rhythm deviation. The smaller the value, the closer the user's sleep cycle structure is to the ideal pattern.

[0071] In an optional embodiment, when generating the baseline evaluation coefficient, the proportion of total sleep time in deep sleep, light sleep, and REM sleep stages throughout the entire sleep cycle is calculated and denoted as follows: , , ;

[0072] The basic evaluation coefficient is obtained by weighting and summing the proportions using preset weights. Basic evaluation coefficient The calculation formula is: ;

[0073] in, , , These are the weighting coefficients for each sleep stage, and their sum is 1.

[0074] Specifically, after analyzing sleep data from the entire night, the total duration of each major sleep stage is calculated. For example, on a night with a total sleep time of 450 minutes, deep sleep lasts 90 minutes, light sleep lasts 240 minutes, and REM sleep lasts 120 minutes. The proportion of each stage to the total sleep time is then calculated. =0.2; ≈0.533; ≈0.267. The weights of different sleep stages on sleep quality are preset, for example, the weight of deep sleep. The weight of light sleep is 0.5. The weight of REM sleep is 0.2. The sum of the three is 1, and the initial value is 0.3. Substitute these values ​​into the formula to calculate the basic evaluation coefficient. The calculated result is approximately 0.2867.

[0075] In an optional embodiment, the formula for generating a comprehensive sleep quality score is as follows: ;

[0076] in, For comprehensive scoring, Based on the evaluation coefficient, This represents the degree of deviation from the periodic rhythm. It is an unsteady-state perturbation factor; , , These are the weighting coefficients for each item. and It is a non-linear adjustment index.

[0077] Specifically, the previously calculated baseline evaluation coefficients were integrated in the final fusion stage. Circadian rhythm deviation and unsteady-state disturbance factors A total score is generated using multiple dimensions of sleep indicators. For example... Figure 4 Assuming the basic evaluation coefficients have been obtained The deviation from the periodic rhythm is 0.2867. It is 0.83, and the unsteady-state perturbation factor The value is 0.5. A corresponding weight and non-linear adjustment index are set for the indicator, for example, It is 200. It is 15. For 20, and the index It is 1.2. The value is 1.5. Nonlinear adjustments are made to the deviation and disturbance factor, and the calculated value is... Approximately equal to 0.799; calculated It is approximately equal to 0.354. Substituting all the values ​​into the formula, the user's overall sleep quality score is 38.275.

[0078] An embodiment of the sleep quality analysis system based on medical bed sensors provided by this invention:

[0079] like Figure 5 As shown, the sleep quality analysis system based on medical bed sensors includes a processor and a memory. The memory stores computer program instructions, which, when executed by the processor, implement the aforementioned sleep quality analysis method based on medical bed sensors.

[0080] The sleep quality analysis system based on sensors in medical beds also includes other components well known to those skilled in the art, such as communication interfaces. Their settings and functions are known in the art and will not be described in detail here.

[0081] In addition, in the description of this specification, "multiple" means at least two, such as two, three or more, etc., unless otherwise expressly and specifically defined.

Claims

1. A sleep quality analysis method based on sensors in a medical bed, characterized in that, Includes the following steps: S1, acquire vital sign signal data collected by sensors on the medical bed during a sleep cycle, and divide it into multiple consecutive time segments; S2, extract the physiological indicators of respiratory rate, heart rate and body movement amplitude in each time segment, and preliminarily determine the sleep stage of the corresponding time segment based on the preset physiological indicator correlation matrix of sleep stage; at the same time, calculate the normalized power spectral entropy of the physiological indicators in each time segment as the variation depth of the corresponding time segment. S3, based on the mean variation depth of the same sleep stage in historical data and the preset stage sensitivity coefficient, generates a variation depth baseline; when the variation depth of a certain time segment exceeds the variation depth baseline of the corresponding sleep stage, the time segment is marked as a micro-awakening disturbance segment. Based on the frequency and duration of micro-awakening disturbances, the non-steady-state disturbance factor is calculated, including: counting the total number of micro-awakening disturbances throughout the entire sleep cycle, dividing by the total number of sleep hours to obtain the frequency. ; Calculate the average duration of all micro-awakening perturbation segments. By analyzing the frequency of occurrence and average duration Perform normalized weighted summation to calculate the unsteady-state perturbation factor. The calculation formula is: ; and These are the preset normalized baseline values ​​for frequency and duration, respectively. and For the corresponding weight coefficients, and ; S4. The determined sleep stage sequence is time-normalized and matched with a standard sleep cycle template to calculate the cycle rhythm deviation. The distribution duration of deep sleep, light sleep, and REM sleep stages throughout the entire sleep cycle is statistically analyzed to generate a baseline assessment coefficient. The baseline assessment coefficient, cycle rhythm deviation, and non-steady-state perturbation factor are fused and a non-linear weighted calculation is performed to generate a comprehensive sleep quality score. When generating the baseline assessment coefficient, the proportion of total sleep time in deep sleep, light sleep, and REM sleep stages throughout the entire sleep cycle is statistically analyzed and denoted as follows: , , The basic evaluation coefficient is obtained by weighting and summing the proportions using preset weights. Basic evaluation coefficient The calculation formula is: ; , , These are the weighting coefficients for each sleep stage, and their sum is 1.

2. The sleep quality analysis method based on medical bed sensors according to claim 1, characterized in that, In S1, the vital sign data is divided into non-overlapping window functions with a time segment length of 30 seconds to generate a series of continuous time segment sequences.

3. The sleep quality analysis method based on medical bed sensors according to claim 1, characterized in that, Physiological parameters such as respiratory rate, heart rate, and body movement amplitude were extracted for each time segment, including: Bandpass filtering is applied to the signal in each time segment to separate heartbeat and respiratory signals; heart rate and respiratory rate are obtained by detecting the peak points of the filtered signals and calculating the peak intervals; the standard deviation of the original signal amplitude within a time segment is calculated as the body motion amplitude index for the corresponding time segment.

4. The sleep quality analysis method based on medical bed sensors according to claim 3, characterized in that, Calculate the normalized power spectral entropy of the physiological indicators within each time segment, including: Perform a Fast Fourier Transform on the heart rate and respiratory rate change sequence within the time segment to obtain the power spectral density distribution; then, extract the energy values ​​of each frequency component of the power spectrum. Divide by total energy The normalized probability is obtained. This forms a probability distribution; according to the Shannon entropy formula... Calculate the information entropy of the probability distribution, and use the information entropy as the mutation depth, where, This represents each individual frequency component in the power spectral density distribution. It represents the total number of all frequency components in the power spectral density distribution.

5. The sleep quality analysis method based on medical bed sensors according to claim 1, characterized in that, When generating the variation depth baseline, at least 50 historical segments in the historical database that belong to the same sleep stage as the current time segment are queried, and the arithmetic mean of the variation depth of the historical segments is calculated. The arithmetic mean is multiplied by an adjustment factor, which is equal to 1 plus the sensitivity coefficient of the corresponding sleep stage, wherein the sensitivity coefficient of the light sleep stage is 0.2, the sensitivity coefficient of the deep sleep stage is 0.3, and the sensitivity coefficient of the REM sleep stage is 0.25, and the variation depth baseline of the sleep stage is generated.

6. The sleep quality analysis method based on medical bed sensors according to any one of claims 1-5, characterized in that, When calculating the deviation of the circadian rhythm, the sleep stages of "deep sleep", "light sleep", "rapid eye movement" and "wakefulness" are mapped to the values ​​4, 3, 2 and 1 respectively, thus forming a numerical sequence of sleep stages; The cumulative normalized path cost is obtained by performing time-warping matching between the numerical sequence of sleep stages and the standard sleep cycle template sequence that has also undergone numerical mapping. The cumulative normalized path cost is divided by the length of the matching path and normalized. This value is used as the periodic rhythm deviation.

7. The sleep quality analysis method based on medical bed sensors according to claim 6, characterized in that, The formula for generating a comprehensive sleep quality score is as follows: ; in, For comprehensive scoring, Based on the evaluation coefficient, This represents the degree of deviation from the periodic rhythm. It is an unsteady-state perturbation factor; , , These are the weighting coefficients for each item. and It is a non-linear adjustment index.

8. A sleep quality analysis system based on sensors in a medical bed, characterized in that, It includes a memory and a processor, wherein the memory stores computer program instructions, and when the computer program instructions are executed by the processor, the sleep quality analysis method based on medical bed sensors as described in any one of claims 1-7 is implemented.

Images