A method and device for processing sleep information

Abstract

The present invention discloses a method and device for processing sleep information, which relates to the field of sleep monitoring and signal processing technology. The method comprises calculating the instantaneous phase of an EEG signal, triggering bone conduction deviation acoustic stimulation when the instantaneous phase is within a target phase window, and dynamically optimizing the stimulation interval through reinforcement learning to obtain a high signal-to-noise ratio EEG signal; constructing a time-varying causal intensity matrix based on the high signal-to-noise ratio EEG signal and time-frequency characteristics, and performing quantum state mapping on the high signal-to-noise ratio EEG signal to extract quantum entanglement characteristics and identify sleep spindles. Simultaneously, reinforced sleep data is generated through an adversarial generative network; and analyzing the persistent homologous features of the signal-to-noise ratio EEG signal using algebraic topological analysis, and combining sleep spindles with reinforced sleep data to assess the user's sleep state through a support vector machine. The present invention effectively enhances data diversity and representativeness, providing more reliable data support for the identification and intervention of sleep disorders.

Description

Technical Field

[0001] The present invention relates to the technical field of sleep monitoring and signal processing, and in particular to a method and device for processing sleep information. Background Art

[0002] In recent years, advances in signal processing technology have introduced methods such as time-frequency analysis and causal strength analysis into sleep monitoring, further improving the accuracy of sleep staging. Furthermore, the application of artificial intelligence technologies, such as support vector machines (SVMs) and deep learning models, has opened up new possibilities for the automatic identification of sleep states. However, existing technologies still have shortcomings in terms of temporal synchronization of multimodal signals, signal noise suppression, and accurate identification of sleep spindles.

[0003] There are two major problems with existing technologies in multimodal sleep monitoring. First, the time synchronization of multimodal signals is insufficient. Although PSG can collect a variety of physiological signals, due to the inconsistency of sampling frequencies and time bases of different sensors, the time error between signals may lead to deviations in the analysis results. Second, the existing technology still has room for improvement in signal noise suppression and feature extraction. Traditional filters have limited effect in processing high-frequency noise and low-frequency baseline drift, and it is difficult to completely eliminate environmental electromagnetic interference, which will affect the signal-to-noise ratio of EEG signals and thus reduce the recognition accuracy of characteristic EEG activities such as sleep spindles. Summary of the Invention

[0004] In view of the above existing problems, the present invention is proposed.

[0005] Therefore, the present invention provides a method for processing sleep information to solve the problems of insufficient signal time synchronization and limited noise suppression effect in the prior art.

[0006] In order to solve the above technical problems, the present invention provides the following technical solutions:

[0007] In a first aspect, the present invention provides a method for processing sleep information, which includes collecting and preprocessing multimodal physiological signals and extracting time-frequency features at the same time, wherein the multimodal physiological signals include EEG signals, eye movement signals, electromyographic signals, electrocardiographic signals and respiratory signals; calculating the instantaneous phase of the EEG signal, triggering bone conduction deviation acoustic stimulation when the instantaneous phase is in the target phase window, and dynamically optimizing the stimulation interval through reinforcement learning to obtain a high signal-to-noise ratio EEG signal; constructing a time-varying causal strength matrix based on the high signal-to-noise ratio EEG signal and time-frequency features, and performing quantum state mapping on the high signal-to-noise ratio EEG signal, extracting quantum entanglement features, identifying sleep spindles, and generating reinforced sleep data through an adversarial generative network; using algebraic topological analysis to analyze the persistent homologous features of the signal-to-noise ratio EEG signal, and combining the sleep spindles and reinforced sleep data to evaluate the user's sleep state through a support vector machine.

[0008] As a preferred embodiment of the sleep information processing method of the present invention, the method further comprises: performing Hilbert transform on the pre-processed EEG signal to calculate the instantaneous phase of the EEG signal;

[0009] The target phase window is set based on the correlation between the rising edge phase of the slow wave in the historical multimodal physiological signal and the maximum response of the EEG signal. As a preferred embodiment of the sleep information processing method of the present invention, the method of dynamically optimizing the stimulation interval through reinforcement learning to obtain a high signal-to-noise ratio EEG signal includes the following steps:

[0010] The Q-learning algorithm is used to dynamically optimize the stimulus interval and initialize the state space and action space;

[0011] Based on the state space and action space, the reward value in the current state is calculated by real-time monitoring of the signal-to-noise ratio of the EEG signal, and the reward function is defined as the increment of the signal-to-noise ratio. At the same time, according to the update rule of Q-learning, Value table;

[0012] Based on the updated The value table selects the optimal action, dynamically adjusts the stimulation interval, and obtains EEG signals with high signal-to-noise ratio.

[0013] As a preferred solution of the sleep information processing method of the present invention, wherein: the extraction of time-frequency features includes the following steps:

[0014] Through Fourier transform, the power characteristics of the δ frequency band and the frequency ratio characteristics of the θ frequency band and the β frequency band are extracted;

[0015] By analyzing the number of saccades in eye movement signals, the saccade density characteristics per unit time are extracted;

[0016] The RMS amplitude feature is extracted by calculating the root mean square value of the electromyographic signal.

[0017] As a preferred embodiment of the sleep information processing method of the present invention, the method comprises the following steps: constructing a time-varying causal strength matrix based on high signal-to-noise ratio EEG signals and time-frequency features, performing quantum state mapping on the high signal-to-noise ratio EEG signals, extracting quantum entanglement features, identifying sleep spindles, and generating enhanced sleep data through a generative adversarial network.

[0018] Calculate the Granger causal strength of eye movement signals on occipital alpha waves, generate a time-varying causal strength matrix based on the Granger causal strength, and identify the dynamic causal relationship between multimodal physiological signals based on the time-varying causal strength matrix;

[0019] Calculate the fidelity of quantum states in the left and right prefrontal cortex, and identify the presence of sleep spindles by combining the dynamic causal relationships between multimodal physiological signals;

[0020] According to the medical record generator and the normal generator, combined with the Lyapunov index, sleep data is generated.

[0021] As a preferred embodiment of the sleep information processing method of the present invention, the analysis of the continuous homologous characteristics of the signal-to-noise ratio EEG signal includes the following steps:

[0022] Perform topological analysis on multi-channel data of high signal-to-noise ratio EEG signals and construct the Vietoris-Rips complex;

[0023] Set the maximum distance parameter of the Vietoris-Rips complex ;

[0024] Maximum distance parameter based on Vietoris-Rips complex , computes the one-dimensional persistent homology features of the Vietoris-Rips complex.

[0025] As a preferred solution of the sleep information processing method of the present invention, wherein: the user's sleep state is evaluated by a support vector machine, comprising the following steps:

[0026] Define the light sleep threshold Z1 and deep sleep threshold Z2;

[0027] The one-dimensional persistent homology features of sleep spindles, enhanced sleep data and Vietoris-Rips complex are integrated into a comprehensive physiological data set and input into the support vector machine to predict the sleep state index and evaluate the user's sleep state.

[0028] In a second aspect, the present invention provides a sleep information processing device, comprising a feature extraction module, an EEG signal optimization module, a sleep data generation module, and a sleep state assessment module; the feature extraction module is used to collect and preprocess multimodal physiological signals and extract time-frequency features;

[0029] The EEG signal optimization module is used to calculate the instantaneous phase of the EEG signal, trigger bone conduction deviation acoustic stimulation when the instantaneous phase is in the target phase window, and dynamically optimize the stimulation interval through reinforcement learning to obtain a high signal-to-noise ratio EEG signal; the sleep data generation module is used to construct a time-varying causal intensity matrix based on the high signal-to-noise ratio EEG signal and time-frequency characteristics, and perform quantum state mapping on the high signal-to-noise ratio EEG signal, extract quantum entanglement features, identify sleep spindles, and generate enhanced sleep data through an adversarial generative network; the sleep state evaluation module is used to use algebraic topological analysis to analyze the persistent homologous features of the signal-to-noise ratio EEG signal, and combine sleep spindles and enhanced sleep data to evaluate the user's sleep state through a support vector machine.

[0030] In a third aspect, the present invention provides a computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein when the computer program is executed by the processor, any step of the sleep information processing method according to the first aspect of the present invention is implemented.

[0031] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein: when the computer program is executed by a processor, any step of the method for processing sleep information as described in the first aspect of the present invention is implemented.

[0032] The beneficial effects of the present invention are as follows: the instantaneous phase of the EEG signal is calculated through Hilbert transform, and bone conduction deviation acoustic stimulation is triggered within the target phase window. The stimulation interval is dynamically optimized in combination with the Q-learning algorithm to maximize the signal-to-noise ratio of the EEG signal. This step significantly improves the accuracy and response speed of sleep monitoring. The high signal-to-noise ratio EEG signal is mapped to quantum entangled features through the quantum state linear combination method, and combined with the adversarial generative network to generate enhanced sleep data containing pathological and healthy features, which effectively enhances the data diversity and representativeness, and provides more reliable data support for the identification and intervention of sleep disorders. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0034] Figure 1 Flowchart of the sleep information processing method in Example 1.

[0035] Figure 2 Schematic diagram of the sleep information processing device in Example 1.

[0036] Figure 3 This is a flow chart of the dynamic optimization of stimulation intervals through reinforcement learning in Example 1.

[0037] Figure 4 This is a flowchart of the support vector machine sleep state evaluation in Example 1. DETAILED DESCRIPTION

[0038] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the specific embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0039] In the following description, many specific details are set forth to facilitate a full understanding of the present invention. However, the present invention may also be implemented in other ways different from those described herein. Those skilled in the art may make similar generalizations without violating the connotation of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.

[0040] Secondly, the term "one embodiment" or "embodiment" herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The phrase "in one embodiment" appearing in various places throughout this specification does not necessarily refer to the same embodiment, nor does it refer to a separate or selective embodiment that is mutually exclusive of other embodiments.

[0041] Example 1, reference Figures 1 to 4 , which is the first embodiment of the present invention, provides a method for processing sleep information, comprising the following steps:

[0042] S1. Collect multimodal physiological signals and preprocess them, and extract time-frequency features.

[0043] Multimodal physiological signals include EEG signals, eye movement signals, EMG signals, ECG signals and respiratory signals.

[0044] Furthermore, EEG signals are collected through an electrode array deployed on the scalp, eye movement signals are recorded through electrodes placed around the eyes, EMG signals are detected through electrodes attached to the muscle surface, ECG signals are measured through electrodes on the chest to measure heart electrical activity, and respiratory signals are monitored through a chest strap or airflow sensor to monitor breathing rate and depth.

[0045] It should be noted that the above data is obtained with the user's consent and is used for legitimate purposes;

[0046] Preprocessing includes signal noise reduction, filtering, and time alignment. The specific steps are as follows:

[0047] First, the collected original multimodal physiological input is fed into a quantum spin Hall effect filter, which eliminates environmental electromagnetic interference (such as 50Hz power frequency noise) to retain the effective multimodal physiological signal components; for example, the high-frequency noise in the EEG signal is significantly suppressed, thereby improving the signal-to-noise ratio of the EEG signal.

[0048] A band-pass filter was applied to the eye movement signal with a frequency band of 0.1–10 Hz to retain effective eye movement features. A high-pass filter was applied to the electromyography signal with a cutoff frequency of 20 Hz to remove low-frequency baseline drift.

[0049] By providing synchronized timestamps for all devices using a unified time reference (such as GPS or network time protocol), the time error of multimodal physiological signals is controlled within ±5ms. For example, EEG signals, eye movement signals, and EMG signals are precisely aligned on the time axis to ensure consistency in multimodal analysis.

[0050] The time-frequency characteristics include delta wave power and theta / beta frequency ratio characteristics, scan density characteristics and RMS amplitude characteristics.

[0051] Furthermore, the EEG signals were Fourier transformed to convert them into frequency domain signals to obtain the power spectral density distribution; the power spectral density of the delta wave band (0.5-4 Hz) was extracted in the frequency domain, and its integral value was calculated as the delta wave power; at the same time, by calculating the power ratio of the theta wave to the beta wave, the power spectral density of the theta wave band (4-8 Hz) and the beta wave band (13-30 Hz) was extracted.

[0052] After preprocessing the eye movement signals, an event detection algorithm is used to identify rapid eye movement (saccade) events. The start and end points of the saccade are determined by setting amplitude and speed thresholds. The number of saccades per unit time is counted and the saccade density is calculated.

[0053] The EMG signal is processed in segments, and the length of each segment is usually set to a fixed time window (such as 1 second). The root mean square value of each EMG signal is calculated as an indicator of the intensity of muscle activity. For example, the root mean square value of the chin EMG signal can be used to monitor muscle tension and assist in identifying micro-arousal events during sleep.

[0054] S2. By calculating the instantaneous phase of the EEG signal, bone conduction deviation acoustic stimulation is triggered when the instantaneous phase is in the target phase window. Through reinforcement learning, the stimulation interval is dynamically optimized to obtain a high signal-to-noise ratio EEG signal.

[0055] Perform Hilbert transform on the preprocessed EEG signal and calculate its instantaneous phase. The expression is:

[0056] ;

[0057] in, Indicates the EEG signal at time The instantaneous phase value at time , is the Cauchy principal value integral, Indicates time The EEG signal, is the integration variable, Indicates the integral variable The EEG signal value at .

[0058] The target phase window is set based on the correlation between the rising edge phase (approximately 0 to π / 2) of the slow wave (0.5-2 Hz) in the historical multimodal physiological signal and the maximum response of the EEG signal. .

[0059] Correlation refers to the positive relationship between the rising edge phase (approximately 0 to π / 2) of the slow wave (0.5-2 Hz) and the response intensity of the EEG signal; for example, when the stimulus is applied at the rising edge phase of the slow wave, the amplitude of the event-related potential (ERP) is significantly higher than that at other phases, indicating that this phase window is the critical time to induce the optimal response of neural activity.

[0060] The instantaneous phase refers to the angular position of the EEG signal at a certain moment, usually expressed in radians (0 to 2π).

[0061] When the instantaneous phase of the EEG signal is detected In target window When the bone conduction deviation sound stimulation is automatically triggered, the stimulation time Depend on OK, among them Indicates the instantaneous phase.

[0062] The Q-learning algorithm is used to dynamically optimize the stimulation interval. The objective function is the maximum signal-to-noise ratio of the EEG signal. The specific steps include:

[0063] Based on the instantaneous phase of EEG signals and the current stimulus interval, initialize the state space;

[0064] Initialize the action space based on the increment or decrement of the inter-stimulus interval;

[0065] Increment and decrement refer to the time of increase and decrease in the stimulus interval;

[0066] By monitoring the signal-to-noise ratio of the EEG signal in real time, the reward value in the current state is calculated, and the reward function is defined as the increment of the signal-to-noise ratio. At the same time, according to the update rule of Q-learning, Value table, the expression is:

[0067] ;

[0068] in, Is in state Next action of value, represents the learning rate, Indicates the current state Next action The reward value obtained after (i.e., the increase in the signal-to-noise ratio of the EEG signal), is the discount factor (range 0 to 1), is the next state, Is the next state The action performed, In the next state All possible moves The largest value.

[0069] Should be explained, update The value table refers to the adjustment of state-action pairs according to the update rules of the Q-learning algorithm. value to reflect the change of expected cumulative reward under the current strategy; specifically, after each action is executed, according to the reward value (i.e., the increment of the EEG signal-to-noise ratio) and the maximum value of the next state Value, update the current state-action pair This process is directly related to the dynamic optimization of the interstimulus interval, because After the value table is updated, the Q-learning algorithm selects the one with the highest The action of optimizing the stimulation interval time (i.e., the optimal stimulation interval time adjustment) gradually approaches the goal of maximizing the SNR of the EEG signal; for example, when a stimulation interval time adjustment leads to a significant improvement in the SNR, the action of The value will increase, and the Q-learning algorithm will select this action to adjust and ultimately maximize the signal-to-noise ratio.

[0070] It should also be noted that action refers to the specific operation of adjusting the stimulation interval, such as increasing or decreasing it by 100 milliseconds, and state refers to the combination of the phase information of the current EEG signal and the stimulation interval.

[0071] Based on the updated The value table selects the optimal action, dynamically adjusts the stimulation interval, and obtains EEG signals with high signal-to-noise ratio.

[0072] S3. Based on high signal-to-noise ratio EEG signals and time-frequency characteristics, a time-varying causal strength matrix is constructed, and quantum state mapping is performed on the high signal-to-noise ratio EEG signals to extract quantum entanglement features and identify sleep spindles. At the same time, enhanced sleep data is generated through a generative adversarial network.

[0073] Time series analysis was performed on high signal-to-noise ratio EEG signals and time-frequency features, and the Granger causal analysis method was used to calculate the Granger causal strength of eye movement signals on occipital alpha waves (8-12 Hz). The expression is:

[0074] ;

[0075] in, Indicates eye movement signals occipital alpha waves The Granger causal strength of Represents historical eye movement signals Historical occipital alpha waves The predicted variance of occipital alpha waves under different conditions, Representation based on historical occipital alpha waves Occipital alpha waves under the conditions The prediction variance of Indicates the current occipital alpha wave, It is the historical occipital alpha wave, is the historical eye movement signal, is the current eye movement signal;

[0076] Furthermore, when calculating the Granger causal strength of the eye movement signal E on the occipital alpha wave g, a fixed time window is used to segment the signal, and the Granger causal strength is calculated separately in each time window.

[0077] The Granger causal strength within each time window is filled into the corresponding position of the causal strength matrix. The rows and columns of the matrix represent different multimodal physiological signals, respectively. As the time window slides, the causal strength matrix is dynamically updated to generate a time-varying causal strength matrix.

[0078] Furthermore, first, multimodal physiological signals (such as EEG signals, eye movement signals and EMG signals) are time-segmented, with each segment length being a fixed time window (such as 1 second); then, the Granger causal strength between different physiological signals is calculated within each time window, such as the causal strength of eye movement signals on occipital alpha waves. , the causal strength of electromyographic signals on electroencephalographic signals, etc.; then, the Granger causal strength values in each time window are filled into the corresponding positions of the causal strength matrix. The rows and columns of the causal strength matrix represent different multimodal physiological signals, respectively. For example, the first row of the time-varying causal strength matrix represents the causal influence of eye movement signals on other multimodal physiological signals, and the second row represents the causal influence of electromyographic signals on other signals; finally, as the time window slides, the causal strength matrix is dynamically updated to generate a time-varying causal strength matrix.

[0079] Based on the time-varying causal strength matrix, the dynamic causal relationship between multimodal physiological signals is identified. For example, in the rapid eye movement (REM) sleep stage, the Granger causal strength of eye movement signals on occipital alpha waves is Significantly enhanced, indicating that eye movement signals have a stronger influence on alpha waves; while in the non-rapid eye movement sleep stage, The value of decreases, indicating that the influence of eye movement signals on alpha waves is weakened. At the same time, the Granger causal strength of myoelectric signals on EEG signals increases significantly during micro-arousal events, indicating that the influence of muscle activity on EEG signals is enhanced.

[0080] Through the quantum state linear combination method, the high signal-to-noise ratio EEG signal is quantum-mapped and converted into quantum entanglement characteristics. The expression is:

[0081] ;

[0082] in, It is the quantum entanglement feature of EEG signals. is the total number of frequency components, is the frequency component In time The phase change when Indicates the The quantum state basis vector corresponding to the frequency component is is the index variable of the frequency component.

[0083] Furthermore, the high signal-to-noise ratio EEG signal is first Fourier transformed, converting it from the time domain to the frequency domain and extracting all frequency components. Each frequency component is then mapped to a corresponding quantum state basis vector, representing the ground state of that frequency component in the quantum state space. The phase change of each frequency component is then weighted and combined with the corresponding quantum state basis vector to generate the quantum entanglement feature of the EEG signal. Finally, normalization is performed to ensure that the total probability of the quantum entanglement feature is 1, generating the final quantum state representation. For example, when analyzing EEG signals from the left and right prefrontal cortex regions, the quantum entanglement features generated by the quantum state linear combination method can be used to assess the similarity between brain regions and thus identify sleep spindles.

[0084] Based on quantum entanglement characteristics , calculate the fidelity of the quantum state of the left and right brain regions of the prefrontal lobe , the expression is:

[0085] ;

[0086] in, represents the quantum state of the left prefrontal brain area, Represents the quantum state of the right prefrontal brain area.

[0087] Based on historical multimodal physiological signal data, we define the fidelity benchmark F1 and the causal strength benchmark G1.

[0088] Fidelity of quantum states based on the left and right prefrontal cortex , combining the dynamic causal relationship between multimodal physiological signals to analyze the interaction between brain regions. >F1 and When <G1, sleep spindles are judged to be present.

[0089] Furthermore, first, the fidelity of the quantum state of the left and right prefrontal cortex is calculated. , evaluate the similarity of the quantum states of the two brain regions; then, calculate the Granger causal strength of the eye movement signal on the occipital alpha wave, and analyze the influence of the eye movement signal on the occipital alpha wave; then, the fidelity of the quantum state Compared with the fidelity benchmark F1, the Granger causal strength is compared with the causal strength benchmark G1; finally, when the fidelity When the fidelity is greater than the reference value F1 and the Granger causal strength is less than the causal strength reference G1, it is determined that sleep spindles exist. When the F1 is 0.9 (greater than the reference value F1=0.85) and the Granger causal strength is 0.6 (less than the reference value G1=0.65), it is determined that sleep spindles exist in the current EEG signal.

[0090] It should be noted that sleep spindles are a characteristic EEG activity, manifested as spindle-shaped oscillation waves of 11-16 Hz, such as regular wave groups lasting 0.5-2 seconds; the value range of the fidelity benchmark F1 is usually 0.8-0.9, and the value range of the causal strength benchmark G1 is usually 0.6-0.7.

[0091] A dual-generator architecture consisting of a pathological generator and a normal generator is used as the adversarial generative network, and the Lyapunov exponent is calculated to constrain the dynamic stability of sleep data. The dual-generator architecture includes a pathological generator and a normal generator.

[0092] The pathological generator generates sleep data containing pathological characteristics, and the normal generator generates sleep data containing healthy characteristics. The sleep data with pathological characteristics and the sleep data with healthy characteristics are combined to generate enhanced sleep data.

[0093] Furthermore, the dual-generator architecture includes a pathological generator and a normal generator. The pathological generator is used to generate sleep data containing pathological features (such as apnea), while the normal generator is used to generate sleep data containing healthy sleep features (such as normal breathing patterns). The Lyapunov exponents of the generated data are then calculated using nonlinear dynamic analysis methods to ensure that the behavior of the generated data on long-term time scales conforms to physiological laws. The pathological and normal generators are then jointly trained so that the generated sleep data contains both diverse pathological and healthy features and meets dynamic stability conditions. Finally, pathological data containing wake events or apnea events is generated to enhance the diversity and representativeness of the sleep dataset. For example, the pathological generator generates sleep data containing apnea features, while the normal generator generates sleep data that conforms to healthy sleep characteristics. Together, the two constitute an enhanced sleep dataset.

[0094] The Lyapunov exponent is calculated by nonlinear dynamic analysis method, and the expression is:

[0095] ;

[0096] in, is the Lyapunov exponent, is the small perturbation vector at the initial moment, It's the time that has passed The perturbation vector after is the norm of the perturbation vector, which indicates the magnitude of the perturbation.

[0097] It should be noted that sleep data containing pathological characteristics is obtained by collecting physiological signals from patients with sleep disorders (such as sleep apnea, insomnia, or periodic limb movement disorder). These signals include EEG signals, eye movement signals, myoelectric signals, electrocardiographic signals, and respiratory signals. These signals are recorded through polysomnography (PSG) and annotated with pathological events (such as the number of respiratory arrests and limb movement frequency). For example, in data from patients with sleep apnea, EEG signals show frequent arousal waves, and respiratory signals show periodic interruptions.

[0098] Sleep data containing health characteristics is obtained by collecting physiological signals from healthy individuals, including EEG signals, eye movement signals, electromyographic signals, electrocardiographic signals, and respiratory signals, and recording them through polysomnography (PSG). For example, in data from healthy individuals, EEG signals show typical sleep spindles and K-complexes, and respiratory signals show a stable rhythm.

[0099] S4. Use algebraic topology analysis to analyze the persistent homologous features of the signal-to-noise ratio EEG signals, and combine them with sleep spindles and enhanced sleep data to evaluate the user's sleep state through support vector machine.

[0100] Using algebraic topological analysis, we conduct topological analysis on multi-channel data of high signal-to-noise ratio EEG signals, construct the Vietoris-Rips complex, and set the maximum distance parameter of the complex based on the data distribution of historical EEG signals. , ensuring that the complex can effectively capture the topology of the signal.

[0101] Maximum distance parameter based on Vietoris-Rips complex , calculate the one-dimensional persistent homology features of the Vietoris-Rips complex , the expression is:

[0102] ;

[0103] in, represents the number of one-dimensional persistent homology features in the Vietoris-Rips complex, is the length of the duration window, is the brain area.

[0104] Furthermore, based on the data distribution of historical EEG signals, the maximum distance between signal points is determined, and the maximum distance parameter of the Vietoris-Rips complex is set accordingly. , ensuring that the Vietoris-Rips complex can effectively capture the topological structure of the signal. Then, according to the maximum distance parameter of the Vietoris-Rips complex The multi-channel data points of the high signal-to-noise ratio EEG signal are connected into a Vietoris-Rips complex to generate a topological structure. Then, the one-dimensional persistent homology group of the complex is calculated to extract the one-dimensional persistent homology features of the Vietoris-Rips complex.

[0105] A brain region is an area of the brain with a specific function or anatomical location, such as the frontal lobe, occipital lobe, or temporal lobe.

[0106] The one-dimensional continuous homology features of sleep spindles, enhanced sleep data and Vietoris-Rips complex are integrated into a comprehensive physiological data set and input into the support vector machine to predict the sleep state index. The expression is:

[0107] ;

[0108] in, is the sleep state index, is the number of support vectors, is the index variable of the support vector, is the bias term of the support vector machine, It is The weight coefficients of the support vectors, Indicates the The labels of the support vectors, Indicates the elements of a comprehensive physiological dataset, Elements of the comprehensive physiological dataset, represents the kernel function, express and The similarity between them.

[0109] Furthermore, based on historical multimodal physiological data, the support vector machine was trained using the backpropagation algorithm and gradient descent method. The specific process is as follows: first, the historical multimodal physiological data are input into the support vector machine, and then the similarity between samples in the training set is calculated by the kernel function to construct a similarity matrix; then, based on the similarity matrix and sample labels, the gradient of the loss function with respect to the weight coefficient and bias term is calculated using the backpropagation algorithm, and the weight coefficient and bias term are updated by the gradient descent method so that the support vector machine can maximize the classification interval; finally, the classification performance of the support vector machine is verified using the test set, and the hyperparameters (such as kernel function parameters and penalty coefficient) are adjusted to optimize the support vector machine.

[0110] express and The similarity between them is expressed as follows:

[0111] ;

[0112] in, Represents the similarity decay factor.

[0113] Based on the distribution of historical sleep data, define the light sleep threshold Z1 and deep sleep threshold Z2;

[0114] when When <Z1, the user is considered to be in the awake period;

[0115] When Z1≤ When <Z2, the user is considered to be in light sleep;

[0116] when When ≥Z2, the user is considered to be in deep sleep.

[0117] It should be noted that the similarity decay factor It is determined based on the performance optimization of the kernel function, with a value range of 0.01 to 1.0, and the specific value is selected through cross-validation; the light sleep threshold Z1 is determined based on the dividing point between the awake period and the light sleep period in the historical sleep data distribution, with a value range of 0.3 to 0.5; the deep sleep threshold Z2 is determined based on the dividing point between the light sleep period and the deep sleep period in the historical sleep data distribution, with a value range of 0.6 to 0.8.

[0118] This embodiment also provides a sleep information processing device, comprising: a feature extraction module, an EEG signal optimization module, a sleep data generation module, and a sleep state assessment module;

[0119] Feature extraction module, used to collect multimodal physiological signals and perform preprocessing, while extracting time-frequency features;

[0120] The EEG signal optimization module is used to calculate the instantaneous phase of the EEG signal. When the instantaneous phase is within the target phase window, bone conduction deviation acoustic stimulation is triggered. Through reinforcement learning, the stimulation interval is dynamically optimized to obtain EEG signals with a high signal-to-noise ratio.

[0121] The sleep data generation module is used to construct a time-varying causal strength matrix based on high-SNR EEG signals and time-frequency features, perform quantum state mapping on high-SNR EEG signals, extract quantum entanglement features, identify sleep spindles, and generate enhanced sleep data through a generative adversarial network.

[0122] The sleep state assessment module is used to analyze the persistent homologous features of the signal-to-noise ratio EEG signals using algebraic topology analysis, and combine sleep spindles and enhanced sleep data to assess the user's sleep state through a support vector machine.

[0123] This embodiment further provides a computer device suitable for the sleep information processing method, comprising: a memory and a processor; the memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions to implement the sleep information processing method proposed in the above embodiment.

[0124] The computer device may be a terminal, comprising a processor, memory, a communication interface, a display, and an input device connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores an operating system and computer programs. The internal memory provides an environment for the operating system and computer programs stored in the non-volatile storage media. The communication interface of the computer device is used to communicate with external terminals via wired or wireless communication. Wireless communication may be achieved via Wi-Fi, a carrier network, NFC (near-field communication), or other technologies. The display of the computer device may be a liquid crystal display or an electronic ink display. The input device may be a touchscreen overlay on the display, buttons, a trackball, or a touchpad on the computer device housing, or an external keyboard, touchpad, or mouse.

[0125] This embodiment further provides a storage medium having a computer program stored thereon. When the program is executed by a processor, the method for processing sleep information as proposed in the above embodiment is implemented. The storage medium can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk.

[0126] In summary, the present invention achieves the maximization of the signal-to-noise ratio of the EEG signal by: calculating the instantaneous phase of the EEG signal through Hilbert transform, triggering bone conduction deviation acoustic stimulation within the target phase window, and dynamically optimizing the stimulation interval time in combination with the Q-learning algorithm. This step significantly improves the accuracy and response speed of sleep monitoring. The high signal-to-noise ratio EEG signal is mapped into quantum entangled features through the quantum state linear combination method, and combined with the adversarial generative network to generate enhanced sleep data containing pathological and healthy features, which effectively enhances the data diversity and representativeness, and provides more reliable data support for the identification and intervention of sleep disorders.

[0127] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit the present invention. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

Claims

1. A method for processing sleep information, characterized by: include, Collecting and preprocessing multimodal physiological signals, while extracting time-frequency features, the multimodal physiological signals including EEG signals, eye movement signals, EMG signals, ECG signals, and respiratory signals; Calculate the instantaneous phase of the EEG signal and trigger bone conduction deviation acoustic stimulation when the instantaneous phase is within the target phase window. Through reinforcement learning, dynamically optimize the stimulation interval to obtain EEG signals with a high signal-to-noise ratio. Based on the time-frequency characteristics of high-SNR EEG signals and multimodal physiological signals, a time-varying causal strength matrix is constructed. The high-SNR EEG signals are then quantum-mapped to extract quantum entanglement features, and sleep spindles are identified. At the same time, enhanced sleep data is generated through a generative adversarial network. The algebraic topology analysis method is used to analyze the persistent homology characteristics of the signal-to-noise ratio EEG signals, and combined with sleep spindles and enhanced sleep data, the user's sleep state is evaluated through support vector machines; The steps of identifying sleep spindles are as follows: First, calculate the fidelity of the quantum states of the left and right prefrontal cortex , evaluate the similarity of the quantum states of the two brain regions; then, calculate the Granger causal strength of the eye movement signal on the occipital alpha wave, and analyze the influence of the eye movement signal on the occipital alpha wave; then, the fidelity of the quantum state Compared with the fidelity benchmark F1, the Granger causal strength is compared with the causal strength benchmark G1; finally, when the fidelity When the value is greater than the reference value F1 and the Granger causal strength is less than the causal strength reference G1, it is determined that sleep spindles are present.

2. The sleep information processing method according to claim 1, wherein: Perform Hilbert transform on the preprocessed EEG signal to calculate the instantaneous phase of the EEG signal; The target phase window is set based on the correlation between the rising edge phase of the slow wave in the historical multimodal physiological signals and the maximum response of the EEG signal.

3. The sleep information processing method according to claim 1, wherein: The method of dynamically optimizing the stimulation interval through reinforcement learning and obtaining a high signal-to-noise ratio EEG signal includes the following steps: The Q-learning algorithm is used to dynamically optimize the stimulus interval and initialize the state space and action space; Based on the state space and action space, the reward value in the current state is calculated by real-time monitoring of the signal-to-noise ratio of the EEG signal, and the reward function is defined as the increment of the signal-to-noise ratio. At the same time, according to the update rule of Q-learning, Value table; Based on the updated The value table selects the optimal action, dynamically adjusts the stimulation interval, and obtains EEG signals with high signal-to-noise ratio.

4. The sleep information processing method according to claim 1, wherein: The time-frequency feature extraction, The following steps are included: Through Fourier transform, the power characteristics of the δ frequency band and the frequency ratio characteristics of the θ frequency band and the β frequency band are extracted; By analyzing the number of saccades in eye movement signals, the saccade density characteristics per unit time are extracted; The RMS amplitude feature is extracted by calculating the root mean square value of the electromyographic signal.

5. The sleep information processing method according to claim 3, wherein: The enhanced sleep data is generated based on a medical record generator and a normal generator in combination with the Lyapunov index.

6. The sleep information processing method according to claim 5, wherein: The analysis of the signal-to-noise ratio of the EEG signal is based on the continuous homologous characteristics of the signal-to-noise ratio. The following steps are included: Perform topological analysis on multi-channel data of high signal-to-noise ratio EEG signals and construct the Vietoris-Rips complex; Set the maximum distance parameter of the Vietoris-Rips complex ; Maximum distance parameter based on Vietoris-Rips complex , computes the one-dimensional persistent homology features of the Vietoris-Rips complex.

7. The sleep information processing method according to claim 6, wherein: The method of evaluating the user's sleep state by using a support vector machine includes the following steps: Define the light sleep threshold Z1 and deep sleep threshold Z2; The one-dimensional persistent homology features of sleep spindles, enhanced sleep data and Vietoris-Rips complex are integrated into a comprehensive physiological data set and input into the support vector machine to predict the sleep state index and evaluate the user's sleep state.

8. A sleep information processing device based on the sleep information processing method according to any one of claims 1 to 7, characterized in that: Including feature extraction module, EEG signal optimization module, sleep data generation module and sleep state assessment module; Feature extraction module, used to collect multimodal physiological signals and perform preprocessing, while extracting time-frequency features; The EEG signal optimization module is used to calculate the instantaneous phase of the EEG signal. When the instantaneous phase is within the target phase window, bone conduction deviation acoustic stimulation is triggered. Through reinforcement learning, the stimulation interval is dynamically optimized to obtain EEG signals with a high signal-to-noise ratio. The sleep data generation module is used to construct a time-varying causal strength matrix based on high-SNR EEG signals and time-frequency features, perform quantum state mapping on high-SNR EEG signals, extract quantum entanglement features, identify sleep spindles, and generate enhanced sleep data through a generative adversarial network. The sleep state assessment module is used to analyze the persistent homologous features of the signal-to-noise ratio EEG signals using algebraic topology analysis, and combine sleep spindles and enhanced sleep data to assess the user's sleep state through a support vector machine.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the sleep information processing method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the sleep information processing method according to any one of claims 1 to 7 are implemented.

Images