Signal processing method and related device

Abstract

The invention provides a signal processing method and a related device. The signal processing method comprises the following steps: acquiring a respiration signal and a heartbeat signal of a target object; determining an instantaneous heart rate signal of the target object according to the heartbeat signal; based on the respiratory signal and the instantaneous heart rate signal, a cardiopulmonary coupling degree index of the target object is determined, and the cardiopulmonary coupling degree index is used for identifying the degree of respiratory sinus arrhythmia of the target object. By means of the method, the respiratory signal and the heartbeat signal of the target object can be synthesized, the cardiopulmonary coupling degree index used for identifying the degree of the respiratory sinus arrhythmia of the target object can be obtained from the angle of cardiopulmonary coupling, and therefore the respiratory sinus arrhythmia of the target object can be analyzed more accurately.

Description

Technical Field

[0001] The present application relates to the field of data processing, and in particular, to a signal processing method and related devices. Background Art

[0002] Respiratory Sinus Arrhythmia (RSA) refers to a normal physiological phenomenon in which the heart rate increases during inspiration and decreases during expiration. RSA can match intermittent expiratory ventilation with relatively continuous pulmonary blood perfusion, namely ventilation / perfusion match (V / QMatch), which means that RSA can improve the efficiency of pulmonary gas exchange by adjusting the heart rate and changing the distribution of cardiac output in the respiratory cycle.

[0003] From the perspective of the interaction between the heart and lungs, respiratory sinus arrhythmia reflects the driving effect of respiration on heart rate and the compensatory effect of heart rate on respiration. This good coupling can effectively ensure that the heart and lungs can efficiently adapt to the wide range of fluctuations in the body's demand for tissue oxygen supply and metabolism under different internal and external environments, and maintain the stability of body functions. When affected by intrinsic pathological factors in the heart or respiratory aspects, the appropriate coupling between the heart and lungs reflected by respiratory sinus arrhythmia can also play a key benign or malignant compensatory role. It can be seen that the appropriate degree of respiratory sinus arrhythmia is a sign of the human body's adaptive maintenance of body homeostasis. In practical applications, quantitative evaluation of the degree of respiratory sinus arrhythmia is of great significance.

[0004] In related technologies, the degree of respiratory sinus arrhythmia is usually evaluated in the analysis of heart rate variability (HRV). Since respiratory sinus arrhythmia is the main source of short-term heart rate variability, standard deviation indicators such as the Standard Deviation of All NN Intervals (SDNN) and the root mean square of the successive differences (RMSSD) of adjacent RR intervals show a certain correlation with the degree of respiratory sinus arrhythmia and can be used to measure the degree of respiratory sinus arrhythmia. However, the above standard deviation indicators such as SDNN and RMSSD are not specific indicators of respiratory sinus arrhythmia. In this regard, although the maximum heart rate and minimum heart rate difference (HRmax-HRmin) indicator in the heart rate variability analysis imposes respiratory cycle constraints on the basis of describing the discrete degree of the instantaneous heart rate sequence, and can more specifically describe the degree of respiratory sinus arrhythmia, the HRmax-HRmin indicator is still a simple statistical indicator. In many cases, it is still unable to accurately assess the degree of respiratory sinus arrhythmia and has a high degree of uncertainty. At the same time, since it is generally believed that the power of the high-frequency band has a strong correlation with the degree of respiratory sinus arrhythmia in the frequency domain heart rate variability analysis, the high-frequency band power index can also be used to measure the degree of respiratory sinus arrhythmia, but the high-frequency band power index is not a specific indicator of respiratory sinus arrhythmia.

[0005] That is, in the relevant technology, the time domain standard deviation indicators and high-frequency power indicators that can reflect the degree of respiratory sinus arrhythmia in heart rate variability analysis are not specific indicators for respiratory sinus arrhythmia, but related indicators in heart rate variability analysis. The correlation between these indicators and the degree of respiratory sinus arrhythmia is usually attributed to the greater contribution of respiratory sinus arrhythmia to short-term heart rate variability. When the degree of respiratory sinus arrhythmia is weak, the above indicators cannot be used to accurately analyze the degree of respiratory sinus arrhythmia. Summary of the invention

[0006] The embodiments of the present application at least provide a signal processing method and related devices, which can integrate the respiratory signal and heartbeat signal of the target object, and obtain a cardiopulmonary coupling index for identifying the degree of respiratory sinus arrhythmia of the target object from the perspective of cardiopulmonary coupling, thereby performing a more accurate analysis of the respiratory sinus arrhythmia of the target object.

[0007] In a first aspect, the present application provides a signal processing method, comprising: Acquire the breathing signal and heartbeat signal of the target object; Determine the instantaneous heart rate signal of the target object according to the heartbeat signal; Based on the breathing signal and the instantaneous heart rate signal, a cardiopulmonary coupling index of the target object is determined, and the cardiopulmonary coupling index is used to identify the degree of respiratory sinus arrhythmia of the target object.

[0008] In a possible implementation, determining the instantaneous heart rate signal of the target object according to the heartbeat signal includes: The instantaneous heart rate signal of the target object is determined according to the RR interval information of the heartbeat signal.

[0009] In a possible implementation, determining the instantaneous heart rate signal of the target object according to the RR interval information of the heartbeat signal includes: The instantaneous heart rate signal of the target object is determined by taking the inverse of the RR interval information of the heartbeat signal and performing spline interpolation.

[0010] In a possible implementation, determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal includes: determining the strength of the time-frequency correlation between the respiration signal and the instantaneous heart rate signal; Based on the strength of the time-frequency correlation, the cardiopulmonary coupling index of the target subject is determined.

[0011] In a possible implementation, determining the time-frequency correlation strength between the respiratory signal and the instantaneous heart rate signal includes: The breathing signal and the instantaneous heart rate signal are transformed into a time-frequency space to obtain a breathing time-frequency signal corresponding to the breathing signal and an instantaneous heart rate time-frequency signal corresponding to the instantaneous heart rate signal; The time-frequency correlation strength is determined based on the same-frequency power strength or the same-phase power strength between the respiratory time-frequency signal and the instantaneous heart rate time-frequency signal.

[0012] In a possible implementation, determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal includes: Determine the causal strength between the breathing signal and the instantaneous heart rate signal, the causal strength is used to identify the unidirectional driving response strength between the breathing and the heart rate of the target object; Based on the causal strength, a cardiopulmonary coupling index of the target subject is determined.

[0013] In a possible implementation, the method further includes: Determine the heartbeat response sensitivity index of the target object based on the breathing signal and the heartbeat signal; The respiratory driving force index of the target subject is determined according to the cardiopulmonary coupling index and the heartbeat response sensitivity index.

[0014] In a possible implementation, determining the heartbeat response sensitivity index of the target object based on the breathing signal and the heartbeat signal includes: Based on the breathing signal, determine the average breathing frequency of the target object; based on the heartbeat signal, determine the average heart rate and average diastolic duration of the target object; The heart rate response sensitivity index is determined by the following formula: RSI=T D ×(f e ) 2 / f r

[0015] Among them, RSI represents the heart rate response sensitivity index, T D represents the average diastolic duration, f e represents the average heart rate, f r Represents the average respiratory rate.

[0016] In a possible implementation, determining the respiratory driving force index of the target object according to the cardiopulmonary coupling index and the heartbeat response sensitivity index includes: The respiratory drive index is determined by the following formula: RDI = (q / RSI) × K

[0017] Among them, RDI represents the respiratory driving force index, q represents the cardiopulmonary coupling index, RSI represents the heartbeat response sensitivity index, and K represents the preset normalization constant.

[0018] In a possible implementation, obtaining a breathing signal and a heartbeat signal of a target object includes: Synchronously collect breathing signals and heartbeat signals.

[0019] In a possible implementation, obtaining a breathing signal and a heartbeat signal of a target object includes: Collect heartbeat signals; Determine the breathing signal based on the heartbeat signal.

[0020] In a second aspect, the present application further provides a signal processing device, comprising: An acquisition module, used to acquire a breathing signal and a heartbeat signal of a target object; A first determination module is used to determine the instantaneous heart rate signal of the target object according to the heartbeat signal; The second determination module is used to determine the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal, where the cardiopulmonary coupling index is used to identify the degree of respiratory sinus arrhythmia of the target object.

[0021] In a third aspect, the present application also provides an electronic device, comprising: a processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor and the memory communicate through the bus, and when the machine-readable instructions are executed by the processor, the signal processing method provided in the present application is executed.

[0022] In a fourth aspect, the present application further provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the signal processing method provided by the present application is executed.

[0023] In a fifth aspect, the present application also provides a computer program product, including a computer program, which executes the signal processing method provided in the present application when executed by a processor.

[0024] In summary, the present application provides a signal processing method and related devices, including: obtaining a breathing signal and a heartbeat signal of a target object; determining the instantaneous heart rate signal of the target object according to the heartbeat signal; determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal, the cardiopulmonary coupling index being used to identify the degree of respiratory sinus arrhythmia of the target object. Through the above method, the breathing signal and the heartbeat signal of the target object can be combined to obtain a cardiopulmonary coupling index for identifying the degree of respiratory sinus arrhythmia of the target object from the perspective of cardiopulmonary coupling, thereby performing a more accurate analysis of the respiratory sinus arrhythmia of the target object.

[0025] Other advantages of the present application will be explained in more detail in conjunction with the following description and drawings.

[0026] It should be understood that the above description is only an overview of the technical solution of the present application, so that the technical means of the present application can be generally understood and then implemented according to the contents of the specification. In order to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the specific implementation methods of the present application are specifically described below by way of example. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following is a brief introduction to the drawings required for use in the embodiments. The drawings herein are incorporated into the specification and constitute a part of the specification. These drawings illustrate embodiments that conform to the present application and are used together with the specification to illustrate the technical solutions of the present application. It should be understood that the drawings only illustrate certain embodiments of the present application and should not be regarded as limiting the scope of protection. For ordinary technicians in this field, other relevant drawings can be obtained based on these drawings without creative work. Moreover, the same reference numerals are used to represent the same components throughout the drawings. In the drawings: Figure 1 A method flow chart of a signal processing method provided in an embodiment of the present application; Figure 2 A method flow chart of another signal processing method provided in an embodiment of the present application; Figure 3 A device diagram of a signal processing device provided in an embodiment of the present application. Specific embodiments

[0028] The exemplary embodiments of the present application will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present application are shown in the accompanying drawings, it should be understood that the present application can be implemented in various forms and should not be limited by the embodiments described herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present application and to fully convey the scope of the present application to those skilled in the art.

[0029] In the description of the embodiments of the present application, it should be understood that terms such as "including" or "having" are intended to indicate the presence of disclosed features, numbers, steps, behaviors, components, parts, or a combination thereof in the present specification, and do not exclude the possibility of the presence of one or more other features, numbers, steps, behaviors, components, parts, or a combination thereof.

[0030] Unless otherwise specified, “ / ” means or. For example, A / B can mean A or B. The “and / or” in this article is merely a way to describe the association relationship of associated objects, indicating that three relationships can exist. For example, A and / or B can mean: A exists alone, A and B exist at the same time, and B exists alone.

[0031] The terms "first", "second", etc. are used only to distinguish the same or similar technical features for the convenience of description, and should not be understood as indicating or implying the relative importance or quantity of these technical features. Thus, the features defined by "first", "second", etc. may explicitly or implicitly include one or more of such features. In the description of the embodiments of the present application, unless otherwise specified, the term "plurality" means two or more than two.

[0032] In the related art, although the degree of respiratory sinus arrhythmia is evaluated in the heart rate variability analysis, the time domain standard deviation indicators and high-frequency power indicators that can reflect the degree of respiratory sinus arrhythmia in the heart rate variability analysis are not specific indicators for respiratory sinus arrhythmia, but related indicators in the heart rate variability analysis. The correlation between these indicators and the degree of respiratory sinus arrhythmia is usually attributed to the greater contribution of respiratory sinus arrhythmia to short-term heart rate variability. When the degree of respiratory sinus arrhythmia is weak, the above indicators cannot be used to accurately analyze the degree of respiratory sinus arrhythmia.

[0033] In view of this, the embodiments of the present application at least provide a signal processing method and related devices, which can integrate the respiratory signal and heartbeat signal of the target object, and obtain a cardiopulmonary coupling index for identifying the degree of respiratory sinus arrhythmia of the target object from the perspective of cardiopulmonary coupling, thereby performing a more accurate analysis of the respiratory sinus arrhythmia of the target object.

[0034] The signal processing method provided in the embodiment of the present application can be implemented by a computer device, which can be a terminal device or a server, wherein the server can be an independent physical server, or a server cluster or distributed system composed of multiple physical servers, or a cloud server providing cloud computing services. Terminal devices include but are not limited to mobile phones, computers, intelligent voice interaction devices, smart home appliances, vehicle terminals, aircraft, etc. The terminal device and the server can be directly or indirectly connected via wired or wireless communication, and this application does not limit this.

[0035] The signal processing method provided by the present application is described below through a method embodiment. Figure 1 As shown, Figure 1 A method flow chart of a signal processing method provided in an embodiment of the present application, wherein the aforementioned computer device may be a server, and the method includes: S101. Acquire a breathing signal and a heartbeat signal of a target object.

[0036] Specifically, the target object is a subject who needs to undergo respiratory sinus arrhythmia detection.

[0037] The respiratory signal refers to a physiological signal caused by respiratory movement. The characteristics of the respiratory signal can be collected through various means, such as using a chest wall respiratory movement sensor, a nasopharyngeal pressure sensor, etc. Preferably, the respiratory signal can be a respiratory impedance signal, or a respiratory movement signal collected by a stress strain sensor, which can reflect the fluctuation of the respiratory cycle. In this application, the respiratory signal can be X r (t) indicates.

[0038] A heartbeat signal refers to a physiological signal that reflects the heartbeat. The heartbeat signal can be recorded and analyzed by medical equipment such as an electrocardiogram (ECG). Preferably, the heartbeat signal can be an electrocardiogram signal so as to reflect the cardiac cycle and the systolic and diastolic periods. In this application, the heartbeat signal can be represented by X e (t) indicates.

[0039] In order for the server to comprehensively analyze the cardiopulmonary coupling degree of the target object in subsequent steps, the server may obtain the breathing signal and the heartbeat signal of the target object.

[0040] In a possible implementation, obtaining a breathing signal and a heartbeat signal of the target object in S101 includes: Synchronously collect breathing signals and heartbeat signals.

[0041] Specifically, the server can directly and synchronously collect the respiratory signal and the heartbeat signal through a collection device including two analog-to-digital converters. In practical applications, the collection time of the collection device can usually be 60s. In order to ensure the synchronous collection of the respiratory signal and the heartbeat signal, the collection front end corresponding to the heartbeat signal and the collection front end corresponding to the respiratory signal need to be controlled by a synchronous clock.

[0042] In a possible implementation, obtaining a breathing signal and a heartbeat signal of the target object in S101 includes: Collect heartbeat signals; Determine the breathing signal based on the heartbeat signal.

[0043] Specifically, the server can first collect the heartbeat signal through an acquisition device including an analog-to-digital converter. Since the rise and fall of the chest during breathing will cause changes in the electrical activity of the heart, the server can determine the respiration signal based on the heartbeat signal. For example, the respiration signal can be determined based on the heartbeat signal by extracting respiration information through an electrocardiogram (ECG-derived Respiration, EDR).

[0044] S102: Determine the instantaneous heart rate signal of the target object according to the heartbeat signal.

[0045] Specifically, after obtaining the breathing signal and heartbeat signal of the target object in S101, in order to perform better cardiopulmonary coupling analysis on the target object later, the server can first determine the instantaneous heart rate signal of the target object according to the heartbeat signal. The instantaneous heart rate signal refers to the rate of each heartbeat of the object. In this application, the instantaneous heart rate signal can be expressed as X h (t) indicates.

[0046] In a possible implementation, determining the instantaneous heart rate signal of the target object according to the heartbeat signal in S102 includes: The instantaneous heart rate signal of the target object is determined according to the RR interval information of the heartbeat signal.

[0047] Specifically, the RR interval is a parameter on the electrocardiogram, which refers to the time interval between the R wave peaks in two consecutive ventricular depolarization complex (QRS) wave groups. The ventricular depolarization complex represents the depolarization process of the ventricles, and the R wave is the most significant peak in the ventricular depolarization complex.

[0048] The server can determine the instantaneous heart rate signal of the target object based on the RR interval information of the heartbeat signal.

[0049] In a possible implementation, determining the instantaneous heart rate signal of the target object according to the RR interval information of the heartbeat signal includes: The instantaneous heart rate signal of the target object is determined by taking the inverse of the RR interval information of the heartbeat signal and performing spline interpolation.

[0050] Specifically, since spline interpolation can be used to estimate and fill in missing heart rate data to estimate the instantaneous heart rate value at each time point, the server can take the inverse of the RR interval information of the heartbeat signal and determine the instantaneous heart rate signal of the target object through spline interpolation, so that the determined instantaneous heart rate signal is of the same length as the breathing signal. For example, the server can use the Pan-Tompkins algorithm to identify the position sequence of the R peaks in the heartbeat signal, and then obtain the initial instantaneous heart rate signal, and then use cubic spline interpolation to interpolate the initial instantaneous heart rate signal into an instantaneous heart rate signal of the same length as the breathing signal to ensure the synchronization between the instantaneous heart rate signal and the breathing signal.

[0051] S103: Determine the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal.

[0052] Specifically, based on the respiratory signal and the instantaneous heart rate signal, the target object's heart and lung are subjected to an overall coupling analysis, and the server can determine the target object's heart and lung coupling index. Since respiratory sinus arrhythmia is mainly manifested as the heart rate accelerating during inhalation and decelerating during exhalation, the heart and lung coupling index, as an index that can identify the degree of coupling between the respiratory signal and the instantaneous heart rate signal, can be used to identify the degree of respiratory sinus arrhythmia of the target object.

[0053] It should be noted that, unlike the heart rate variability analysis in the related art which only relies on the heartbeat signal to conduct an empirical analysis of the degree of respiratory sinus arrhythmia, the heartbeat signal is processed in the present application to obtain the corresponding instantaneous heart rate signal, and an overall cardiopulmonary coupling analysis is performed through the instantaneous heart rate signal and the respiratory signal to obtain the cardiopulmonary coupling index, so as to conduct a more accurate analysis of the target object's respiratory sinus arrhythmia through the cardiopulmonary coupling index.

[0054] In a possible implementation, determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal includes: determining the strength of the time-frequency correlation between the respiration signal and the instantaneous heart rate signal; Based on the strength of the time-frequency correlation, the cardiopulmonary coupling index of the target subject is determined.

[0055] Specifically, the server may first determine the time-frequency correlation strength between the respiratory signal and the instantaneous heart rate signal. The time-frequency correlation strength is used to identify the correlation between the respiratory signal and the instantaneous heart rate signal from two dimensions: time and frequency.

[0056] After determining the time-frequency correlation, the server may determine the cardiopulmonary coupling index of the target object based on the time-frequency correlation strength.

[0057] In a possible implementation, determining the time-frequency correlation strength between the respiratory signal and the instantaneous heart rate signal includes: The breathing signal and the instantaneous heart rate signal are transformed into a time-frequency space to obtain a breathing time-frequency signal corresponding to the breathing signal and an instantaneous heart rate time-frequency signal corresponding to the instantaneous heart rate signal; The time-frequency correlation strength is determined based on the same-frequency power strength or the same-phase power strength between the respiratory time-frequency signal and the instantaneous heart rate time-frequency signal.

[0058] Specifically, in order to analyze the breathing signal and the instantaneous heart rate signal from the two dimensions of time and frequency, the server can transform the breathing signal and the instantaneous heart rate signal into the time-frequency space to obtain the breathing time-frequency signal corresponding to the breathing signal and the instantaneous heart rate time-frequency signal corresponding to the instantaneous heart rate signal. In this application, the breathing time-frequency signal in the time-frequency space can be expressed as S r (t,w), the instantaneous heart rate time-frequency signal can be expressed as S h (t,w).

[0059] After obtaining the respiratory time-frequency signal and the instantaneous heart rate time-frequency signal, the server can determine the time-frequency correlation strength based on the same-frequency power strength or same-phase power strength between the respiratory time-frequency signal and the instantaneous heart rate time-frequency signal, that is, analyze from two dimensions of time and frequency.

[0060] In a possible implementation, the time-frequency correlation intensity quantification method used by the server to calculate the cardiopulmonary coupling degree may be a cross wavelet transform.

[0061] Specifically, the server can perform cross wavelet transform on the instantaneous heart rate signal and the respiratory signal, that is, calculate the average power density of the instantaneous heart rate signal and the respiratory signal in the cross wavelet time spectrum as the cardiopulmonary coupling degree, wherein the time window width used for calculating the cross wavelet time spectrum can be 60 s, and the selected frequency domain scale range can be 0.2 Hz~0.5 Hz, that is, the general respiratory frequency range.

[0062] In a possible implementation, determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal includes: Determine the causal strength between the breathing signal and the instantaneous heart rate signal, the causal strength is used to identify the unidirectional driving response strength between the breathing and the heart rate of the target object; Based on the causal strength, a cardiopulmonary coupling index of the target subject is determined.

[0063] Specifically, since respiratory sinus arrhythmia is mainly manifested as the heart rate accelerating during inhalation and decelerating during exhalation, the respiratory signal and the instantaneous heart rate signal can be processed from the dimension of the one-way driving response between respiration and heart rate, that is, the server can determine the causal strength between the respiratory signal and the instantaneous heart rate signal, and the causal strength is used to identify the one-way driving response strength between the respiration and heart rate of the target object.

[0064] After determining the causal strength, the server may determine a cardiopulmonary coupling index of the target object based on the causal strength.

[0065] In this embodiment, the server can determine the causal strength between the respiratory signal and the instantaneous heart rate signal through the following four methods: (1) Granger causality analysis based on time series prediction technology; (2) transfer entropy method based on conditional probability distribution estimation; (3) convergent cross mapping method based on phase space reconstruction; (4) coupling term estimation method based on multi-oscillator dynamics coupling model (such as Winfree coupling model).

[0066] In a possible implementation, the server may determine the causal strength between the respiratory signal and the instantaneous heart rate signal by a convergent cross-mapping method of phase space reconstruction.

[0067] The convergent cross-mapping method refers to the directional manifold structure comparison of two signals in the phase space to analyze the causal strength between the two. Simply put, compared with the time-frequency correlation evaluation method, this method further introduces directional information in the process of similarity comparison of the phase space structure, thereby being able to reflect the causality between the two signals.

[0068] Although both the respiratory signal and the instantaneous heart rate signal are nonlinear signals, they both have obvious periodic rhythmicity, so they can be reconstructed into a clearer phase space manifold structure diagram through a simple delayed coordinate method; compared with methods based on linear model assumptions such as Granger causality analysis, the convergent cross mapping method, which operates in phase space, does not require linear model assumptions, is particularly suitable for nonlinear signal analysis, and has a small amount of calculation, so it is the preferred method.

[0069] After accurately and comprehensively analyzing the respiratory sinus arrhythmia of the target object from the perspective of cardiopulmonary coupling, in a possible implementation, the cardiopulmonary coupling index may also be decoupled and quantified, that is, the signal processing method may also include: S201, determining a heartbeat response sensitivity index of a target object based on a breathing signal and a heartbeat signal; S202: Determine the respiratory driving force index of the target object according to the cardiopulmonary coupling index and the heartbeat response sensitivity index.

[0070] Specifically, through the aforementioned steps S101-S103, the cardiopulmonary coupling index of the target object is obtained from the perspective of the cardiopulmonary coupling of the target object, thereby performing an overall analysis of the respiratory sinus arrhythmia of the target object. On this basis, since respiratory sinus arrhythmia is mainly manifested as the heart rate accelerating during inhalation and decelerating during exhalation, that is, breathing is the driver and the heart rate is the follower, the cardiopulmonary coupling index can be expressed as the product of the respiratory drive index (Respiratory Drive Index, RDI) and the heartbeat response sensitivity index (Response Sensitivity Index, RSI), wherein the respiratory drive index is an indicator used to measure the respiratory drive strength, which can be used to reflect the contribution of the respiratory aspect to the respiratory sinus arrhythmia phenomenon, and the heartbeat response sensitivity index can be used to describe the sensitivity of the heartbeat to certain stimuli or condition changes, which can be used to reflect the contribution of the heart to the respiratory sinus arrhythmia phenomenon.

[0071] Based on this, the server can first determine the target object's heartbeat response sensitivity index based on the respiratory signal and the heartbeat signal, and then determine the target object's respiratory drive index based on the cardiopulmonary coupling index and the heartbeat response sensitivity index. That is to say, in this embodiment, the cardiopulmonary coupling index is decoupled and quantified from the multiplicative coupling perspective of "driver-slave", so as to obtain a respiratory drive index that can independently evaluate the contribution of the respiratory aspect and a heartbeat response sensitivity index that independently evaluates the contribution of the heartbeat aspect, so that the target object's respiratory sinus arrhythmia can be more accurately analyzed through the cardiopulmonary coupling index, the respiratory drive index and the heartbeat response sensitivity index.

[0072] In a possible implementation, determining the heartbeat response sensitivity index of the target object based on the breathing signal and the heartbeat signal in S201 includes: Based on the breathing signal, determine the average breathing frequency of the target object; based on the heartbeat signal, determine the average heart rate and average diastolic duration of the target object; The heart rate response sensitivity index is determined by the following formula: RSI=T D ×(f e ) 2 / f r

[0073] Among them, RSI represents the heart rate response sensitivity index, T Drepresents the average diastolic duration, f e represents the average heart rate, f r Represents the average respiratory rate.

[0074] Specifically, the heartbeat response sensitivity index can be calculated through the breathing signal and heartbeat signal of the target object, that is, the server can determine the average breathing frequency of the target object based on the breathing signal, and determine the average heart rate and average diastolic duration of the target object based on the heartbeat signal. For example, the server can determine the average breathing frequency of the target object based on the breathing signal with a window width of 60s, and determine the average heart rate and average diastolic duration of the target object based on the heartbeat signal with a window width of 60s; the server can then determine the heartbeat response sensitivity index through the above formula.

[0075] In a possible implementation, determining the respiratory driving force index of the target object according to the cardiopulmonary coupling index and the heartbeat response sensitivity index in S202 includes: The respiratory drive index is determined by the following formula: RDI = (q / RSI) × K

[0076] Among them, RDI represents the respiratory driving force index, q represents the cardiopulmonary coupling index, RSI represents the heartbeat response sensitivity index, and K represents the preset normalization constant.

[0077] Specifically, according to the multiplicative coupling hypothesis of "cardiopulmonary coupling index = respiratory drive index × heartbeat response sensitivity index" in the respiratory sinus arrhythmia phenomenon, the server can determine the respiratory drive index through the above formula, wherein, in practical applications, the preset normalization constant can be 3-5, that is, the ratio of human heart rate to respiratory frequency f e / f r Experience value range.

[0078] It can be seen that the present application provides a signal processing method, including: obtaining the breathing signal and heartbeat signal of the target object; determining the instantaneous heart rate signal of the target object according to the heartbeat signal; determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal, and the cardiopulmonary coupling index is used to identify the degree of respiratory sinus arrhythmia of the target object. Through the above method, the breathing signal and heartbeat signal of the target object can be combined to obtain the cardiopulmonary coupling index used to identify the degree of respiratory sinus arrhythmia of the target object from the perspective of cardiopulmonary coupling, thereby performing a more accurate analysis of the respiratory sinus arrhythmia of the target object.

[0079] In addition, in a further preferred embodiment, the cardiopulmonary coupling index is decoupled and quantified from the perspective of multiplicative coupling of "driver-slave", thereby obtaining a respiratory drive index that can independently evaluate the contribution of the respiratory aspect and a heartbeat response sensitivity index that can independently evaluate the contribution of the heartbeat aspect, so that a more comprehensive and accurate analysis of the target object's respiratory sinus arrhythmia can be performed through the cardiopulmonary coupling index, the respiratory drive index and the heartbeat response sensitivity index.

[0080] In the description of this specification, the description with reference to the terms "some possible embodiments", "some embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application, and the above terms do not necessarily represent the same embodiment or example. Moreover, the described specific features, structures, materials or characteristics may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are contradictory.

[0081] About the method flow chart of the present application embodiment, some operations are described as different steps performed in a certain order. Such flow chart belongs to illustrative and non-restrictive. Some steps described in this article can be grouped together and performed in a single operation, or some steps can be divided into multiple sub-steps and can be performed in an order different from that shown in this article. Each step shown in the flow chart can be realized in any way by any circuit structure and / or tangible mechanism (for example, by software, hardware (for example, the logical function realized by processor or chip) etc. running on computer equipment and / or any combination thereof).

[0082] Those skilled in the art will appreciate that, in the method described in the above specific implementation, the writing order of each step does not mean a strict execution order, and the specific execution order of each step should be determined by its function and possible internal logic.

[0083] Based on the aforementioned Figure 1 and Figure 2 The present application is described below through a device embodiment. Figure 3 As shown, the signal processing device 300 includes: An acquisition module 301 is used to acquire a breathing signal and a heartbeat signal of a target object; A first determination module 302 is used to determine the instantaneous heart rate signal of the target object according to the heartbeat signal; The second determination module 303 is used to determine the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal. The cardiopulmonary coupling index is used to identify the degree of respiratory sinus arrhythmia of the target object.

[0084] In a possible implementation, the first determining module 302 is configured to: The instantaneous heart rate signal of the target object is determined according to the RR interval information of the heartbeat signal.

[0085] In a possible implementation, the first determining module 302 is configured to: The instantaneous heart rate signal of the target object is determined by taking the inverse of the RR interval information of the heartbeat signal and performing spline interpolation.

[0086] In a possible implementation, the second determining module 303 is configured to: determining the strength of the time-frequency correlation between the respiration signal and the instantaneous heart rate signal; Based on the strength of the time-frequency correlation, the cardiopulmonary coupling index of the target subject is determined.

[0087] In a possible implementation, the second determining module 303 is configured to: The breathing signal and the instantaneous heart rate signal are transformed into a time-frequency space to obtain a breathing time-frequency signal corresponding to the breathing signal and an instantaneous heart rate time-frequency signal corresponding to the instantaneous heart rate signal; The time-frequency correlation strength is determined based on the same-frequency power strength or the same-phase power strength between the respiratory time-frequency signal and the instantaneous heart rate time-frequency signal.

[0088] In a possible implementation, the second determining module 303 is configured to: Determine the causal strength between the breathing signal and the instantaneous heart rate signal, the causal strength is used to identify the unidirectional driving response strength between the breathing and the heart rate of the target object; Based on the causal strength, a cardiopulmonary coupling index of the target subject is determined.

[0089] In a possible implementation, the second determining module 303 is further configured to: Determine the heartbeat response sensitivity index of the target object based on the breathing signal and the heartbeat signal; The respiratory driving force index of the target subject is determined according to the cardiopulmonary coupling index and the heartbeat response sensitivity index.

[0090] In a possible implementation, the second determining module 303 is further configured to: Based on the breathing signal, determine the average breathing frequency of the target object; based on the heartbeat signal, determine the average heart rate and average diastolic duration of the target object; The heart rate response sensitivity index is determined by the following formula: RSI=T D ×(f e ) 2 / f r

[0091] Among them, RSI represents the heart rate response sensitivity index, T D represents the average diastolic duration, f e represents the average heart rate, f r Represents the average respiratory rate.

[0092] In a possible implementation, the second determining module 303 is further configured to: The respiratory drive index is determined by the following formula: RDI = (q / RSI) × K

[0093] Among them, RDI represents the respiratory driving force index, q represents the cardiopulmonary coupling index, RSI represents the heartbeat response sensitivity index, and K represents the preset normalization constant.

[0094] In a possible implementation, the acquisition module 301 is used to: Synchronously collect breathing signals and heartbeat signals.

[0095] In a possible implementation, the acquisition module 301 is used to: Collect heartbeat signals; Determine the breathing signal based on the heartbeat signal.

[0096] It should be noted that the device in the implementation mode of the present application can implement each process of the implementation mode of the aforementioned method and achieve the same effects and functions, which will not be repeated here.

[0097] The embodiment of the present application also provides an electronic device, including: a processor, a memory and a bus. The memory stores machine-readable instructions executable by the processor. When the electronic device is running, the processor and the memory communicate through the bus. When the machine-readable instructions are executed by the processor, the following processing is performed: Acquire the breathing signal and heartbeat signal of the target object; Determine the instantaneous heart rate signal of the target object according to the heartbeat signal; Based on the breathing signal and the instantaneous heart rate signal, a cardiopulmonary coupling index of the target object is determined, and the cardiopulmonary coupling index is used to identify the degree of respiratory sinus arrhythmia of the target object.

[0098] The present application also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the signal processing method described in the above method embodiment are executed. The storage medium can be a volatile or non-volatile computer-readable storage medium.

[0099] An embodiment of the present application also provides a computer program product, including a computer program. The computer program product carries a program code. The instructions included in the program code can be used to execute the steps of the signal processing method described in the above method embodiment. For details, please refer to the above method embodiment, which will not be repeated here.

[0100] The computer program product may be implemented in hardware, software or a combination thereof. In one optional embodiment, the computer program product is implemented as a computer storage medium. In another optional embodiment, the computer program product is implemented as a software product, such as a software development kit (SDK).

[0101] Each embodiment in this application is described in a progressive manner, and the same or similar parts between the various embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the device, equipment, and computer-readable storage medium embodiments, since they are basically similar to the method embodiments, their descriptions are simplified, and the relevant parts can be referred to the partial description of the method embodiments.

[0102] The apparatus, equipment and computer-readable storage medium provided in the embodiments of the present application correspond one-to-one to the method. Therefore, the apparatus, equipment and computer-readable storage medium also have similar beneficial technical effects as the corresponding methods. Since the beneficial technical effects of the methods have been described in detail above, the beneficial technical effects of the apparatus, equipment and computer-readable storage medium will not be repeated here.

[0103] Although the spirit and principle of the present application have been described above with reference to several specific embodiments, it should be understood that the present application is not limited to the disclosed specific embodiments, and the division of various aspects does not mean that the features in these aspects cannot be combined. The present application is intended to cover various modifications and equivalent arrangements included in the spirit and scope of the attached claims.

Claims

1. A signal processing method, characterized in that: include: Acquire the breathing signal and heartbeat signal of the target object; Determining the instantaneous heart rate signal of the target object according to the heartbeat signal; A cardiopulmonary coupling index of the target object is determined based on the breathing signal and the instantaneous heart rate signal, where the cardiopulmonary coupling index is used to identify the degree of respiratory sinus arrhythmia of the target object.

2. The method according to claim 1, characterized in that Determining the instantaneous heart rate signal of the target object according to the heartbeat signal includes: The instantaneous heart rate signal of the target object is determined according to the RR interval information of the heartbeat signal.

3. The method according to claim 2, characterized in that The step of determining the instantaneous heart rate signal of the target object according to the RR interval information of the heartbeat signal comprises: The instantaneous heart rate signal of the target object is determined by taking the inverse of the RR interval information of the heartbeat signal and performing spline interpolation.

4. The method according to claim 1, characterized in that: Determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal includes: Determining the time-frequency correlation strength between the respiration signal and the instantaneous heart rate signal; A cardiopulmonary coupling index of the target object is determined based on the time-frequency correlation strength.

5. The method according to claim 4, characterized in that The determining the time-frequency correlation strength between the breathing signal and the instantaneous heart rate signal comprises: Transforming the breathing signal and the instantaneous heart rate signal into a time-frequency space to obtain a breathing time-frequency signal corresponding to the breathing signal and an instantaneous heart rate time-frequency signal corresponding to the instantaneous heart rate signal; The time-frequency correlation strength is determined based on the same-frequency power strength or the same-phase power strength between the respiratory time-frequency signal and the instantaneous heart rate time-frequency signal.

6. The method according to claim 1, characterized in that Determining the cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal includes: Determining a causal strength between the respiration signal and the instantaneous heart rate signal, the causal strength being used to identify a one-way driving response strength between the respiration and the heart rate of the target object; Based on the causality strength, a cardiopulmonary coupling index of the target subject is determined.

7. The method according to claim 1, characterized in that The method further comprises: Determining a heartbeat response sensitivity index of the target object based on the breathing signal and the heartbeat signal; A respiratory driving force index of the target object is determined according to the cardiopulmonary coupling index and the heartbeat response sensitivity index.

8. The method according to claim 7, characterized in that The step of determining the heartbeat response sensitivity index of the target object based on the breathing signal and the heartbeat signal includes: Based on the breathing signal, determine the average breathing frequency of the target object; based on the heartbeat signal, determine the average heart rate and average diastolic duration of the target object; The heartbeat response sensitivity index is determined by the following formula: RSI=T D ×(f e ) 2 / f r Wherein, RSI represents the heartbeat response sensitivity index, T D represents the average diastolic duration, f e represents the average heart rate, f r represents the average respiratory rate.

9. The method according to claim 7, characterized in that: Determining the respiratory driving force index of the target object according to the cardiopulmonary coupling index and the heartbeat response sensitivity index includes: The respiratory drive index is determined by the following formula: RDI=(q / RSI)×K Wherein, RDI represents the respiratory driving force index, q represents the cardiopulmonary coupling index, RSI represents the heartbeat response sensitivity index, and K represents a preset normalization constant.

10. The method according to claim 1, characterized in that The step of obtaining the breathing signal and the heartbeat signal of the target object includes: The breathing signal and the heartbeat signal are collected synchronously.

11. The method according to claim 1, characterized in that The step of obtaining the breathing signal and the heartbeat signal of the target object includes: collecting the heartbeat signal; The breathing signal is determined according to the heartbeat signal.

12. A signal processing device, characterized in that: include: An acquisition module, used to acquire a breathing signal and a heartbeat signal of a target object; A first determining module, configured to determine an instantaneous heart rate signal of the target object according to the heartbeat signal; The second determination module is used to determine a cardiopulmonary coupling index of the target object based on the breathing signal and the instantaneous heart rate signal, wherein the cardiopulmonary coupling index is used to identify the degree of respiratory sinus arrhythmia of the target object.

13. An electronic device, characterized in that: include: A processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor and the memory communicate via the bus, and when the machine-readable instructions are executed by the processor, the signal processing method as described in any one of claims 1 to 11 is performed.

14. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the signal processing method according to any one of claims 1 to 11 is executed.

15. A computer program product, characterized in that The invention comprises a computer program, which, when executed by a processor, performs the signal processing method according to any one of claims 1 to 11.

Images