Portable autonomic nerve function evaluation device, method, equipment and program product

Abstract

The invention discloses a portable autonomic nerve function evaluation device, method, equipment and program product, and is specifically characterized in that the breathing of a user is individually adjusted to a specific frequency through an intelligent breathing monitoring and guiding program, and the heart rate fluctuation condition of the user in the process is recorded by a heart rate acquisition device; the heart and lung coupling strength of a user is quantitatively analyzed based on a feature extraction algorithm, and a mapping model between physiological features and clinical indexes is constructed to realize quantitative evaluation and portable monitoring of autonomic nerve functions. Based on accurate decoding of the cardiopulmonary coupling process, the portable autonomic nerve function evaluation device, method, equipment and program product which are rapid, low in manufacturing cost and easy and convenient to operate and have popularization significance are achieved, and the needed measuring device is simple and convenient and can be used for individualized home health evaluation.

Description

Technical Field

[0001] The invention relates to a portable evaluation technology for autonomic nervous function, and belongs to the technical field of signal processing. Background Art

[0002] The autonomic nervous system is closely related to the regulation of heart rate, blood pressure, digestion, urination, and body temperature, and is essential for maintaining the balance of the body's internal environment. In a fast-paced society, as the sub-health status of individuals gradually deepens, the negative impact of autonomic nervous function imbalance or disorder on the quality of life becomes more and more significant. Therefore, autonomic nervous function monitoring has guiding significance in assessing health levels and adjusting lifestyles.

[0003] The measurement principle of heart rate variability (HRV) is simple and is often used to monitor autonomic nervous function, but this analysis method is greatly affected by individual differences. HRV analysis often relies on a longer signal monitoring time in the evaluation of autonomic nervous system function, especially in the evaluation of sympathetic nervous function. In general, in order to achieve more accurate quantitative analysis of autonomic nervous function, researchers will conduct a comprehensive analysis based on indicators such as blood pressure variability (BPV) or pre-ejection period (PEP). However, BPV analysis requires the use of a continuous dynamic blood pressure monitor, and PEP analysis requires the use of an impedance electrocardiogram. Both of these are complex in principle and expensive, making them difficult to use for individualized daily health management. Therefore, a sensitive, accurate, and economical method for evaluating autonomic nervous function is needed. Summary of the invention

[0004] Purpose of the invention: In view of the above-mentioned deficiencies in the prior art, the purpose of the present invention is to provide a portable assessment device, method, equipment and program product for autonomic nervous function, so as to realize quantitative assessment and portable monitoring of autonomic nervous function.

[0005] Technical solution: To achieve the above-mentioned invention object, the present invention adopts the following technical solution:

[0006] A portable device for evaluating autonomic nervous function, comprising:

[0007] The heart rate extraction and preprocessing module is used to collect physiological signals and obtain the user's heart rate fluctuation through preprocessing operations;

[0008] Intelligent breathing test and guidance module, used to adaptively measure the user's breathing ability, determine the personalized breathing rhythm, and guide the user to follow the breathing rhythm;

[0009] The cardiopulmonary coupling feature extraction module is used to perform spectrum analysis and causal analysis on the user's heart rate signal during breathing training, and calculate the cardiopulmonary coupling index, including the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during long inhalation and short exhalation, and the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during short inhalation and long exhalation;

[0010] The autonomic nervous function analysis module is used to integrate physiological characteristics and determine the user's autonomic nervous system state through an autonomic nervous balance analysis model and an autonomic nervous activity analysis model; the autonomic nervous balance is characterized by the ratio of the cardiopulmonary follow-up response index in different processes, and the autonomic nervous activity is characterized by the cardiopulmonary coupling strength.

[0011] Furthermore, in the heart rate extraction and preprocessing module, the preprocessing operation steps specifically include:

[0012] Find the peak point in the physiological signal and obtain the peak point sequence {r(n), n=1, 2, ..., N} derived from the physiological signal, where N is the total number of points in the sequence;

[0013] Performing a first-order difference operation on the peak point sequence to calculate the peak point interval sequence {rr(n), n=1, 2, ..., N-1};

[0014] The user's beat-to-beat heart rate is calculated based on the peak point interval sequence:

[0015] Furthermore, in the intelligent breathing test and guidance module, the steps of measuring breathing capacity and formulating and guiding breathing rhythm specifically include:

[0016] Introducing the designated breathing method and rhythm through voice and / or video, guiding the user to continue to inhale evenly until no more gas can be inhaled, and then continue to exhale evenly until no more gas can be exhaled;

[0017] The first-order difference of the user's heart rate during the calculation process is: hr′(n)=hr(n)-hr(n-1), n=2,...,N-1, where N is the total number of peak points in the physiological signal sequence, hr i The user's beat-to-beat heart rate;

[0018] The first point in the sequence where hr′(n) ≥ 2bmp is defined as the starting inhalation key frame n. 1 ; The first point after the inhalation key frame in time sequence and with hr′(n)≤0.5bmp is defined as the transition key frame n 2 ; The first point after the transition key frame in time sequence and hr′(n)≤-2bmp is defined as the starting exhalation key frame n 2 ; The first point after the start exhalation key frame in time sequence and with hr′(n)≥-0.5bmp is defined as the end key frame n 4 ;

[0019] Calculate the user's maximum constant inspiratory time: And maximum constant exhalation time: in to is the time corresponding to the corresponding key frame;

[0020] The user's breathing rhythm is set in a personalized and adaptive manner based on the measured maximum uniform inhalation time and maximum uniform exhalation time of the user, and the user is guided to breathe according to the set rhythm through voice, vibration and / or video, including long inhalation and short exhalation training and short inhalation and long exhalation training.

[0021] Furthermore, in the cardiopulmonary coupling feature extraction module, the steps of calculating the cardiopulmonary follow-up response index and the cardiopulmonary coupling strength specifically include:

[0022] The heart rate fluctuation is converted from the time domain to the frequency domain by discrete Fourier transform:

[0023]

[0024] Where i is an imaginary unit, N is the number of sampling points, hr(n) is the heart rate signal, and HR[k] is the frequency domain distribution of the heart rate signal;

[0025] Calculate the power of each point in the discrete frequency domain to obtain the power spectrum density of the heart rate signal: where f s is the sampling frequency of the original physiological signal.

[0026] Calculate the cardiopulmonary response index during long inhalation and short exhalation: where f 1 It is the fixed respiratory rate during long inhalation and short exhalation, peak s1 Indicates the maximum frequency of the peak point of the heart rate signal spectrum, C 1 yes The integer part of ;

[0027] Calculate the cardiopulmonary response index during short inhalation and long exhalation: where f 2 It is the fixed breathing rate during short inhalation and long exhalation, peak s2 Indicates the maximum frequency of the peak point of the heart rate signal spectrum, C 2 yes The integer part of ;

[0028] The heart rate signal hr(n) is decomposed into multiple physically meaningful components HR through EEMD j (n), j = 1, 2, ..., m;

[0029] The respiratory timing signal is transmitted through EEMD Disassemble into multiple components with physical meaning RESP j (n), j = 1, 2, ..., m;

[0030] HR is obtained by Hilbert transform j The phase φ of (n) x (n), and RESP j The phase φ of (n) y (n);

[0031] Calculate HR j (n) and RESP j (n) coupling strength: Where i is the imaginary unit;

[0032] The cardiopulmonary coupling strength of the user is calculated as: CCS = max{Coh 1 , Coh 2 , ..., Coh m}.

[0033] Furthermore, in the autonomic nervous function analysis module, the step of determining the state of the user's autonomic nervous system specifically includes:

[0034] Calculates the user's autonomic balance: CRI 1 CRI is the cardiopulmonary response index during long inhalation and short exhalation. 2 It is the cardiopulmonary response index during short inhalation and long exhalation;

[0035] Based on the statistical analysis results, determine whether the user's autonomic nervous system state deviates from the normal range;

[0036] If the user's autonomic nervous system state deviates from the normal range, the autonomic nervous system activity is calculated: CCS 1 CCS is the cardiopulmonary coupling strength during long inhalation and short exhalation. 2 It is the cardiopulmonary coupling strength during short inhalation and long exhalation; SNA represents the activity of sympathetic nerves, and PNA represents the activity of parasympathetic nerves;

[0037] A feature vector [SNA, SNA] is constructed, and based on the feature vector, the distance between the user's autonomic nervous activity and the respective autonomic nervous dysfunction types obtained by clustering is calculated, so as to determine the specific type of the user's autonomic nervous dysfunction.

[0038] Furthermore, the autonomic dysfunction type is obtained by performing K-means clustering on the collected sample data, and the number of clusters is set to 6, corresponding to 6 types of autonomic imbalance:

[0039] a. Biased towards the sympathetic side: sympathetic nerves are strong and parasympathetic nerves are normal;

[0040] b. Biased towards the sympathetic side: the sympathetic nerves are normal, and the parasympathetic nerves are weak;

[0041] c. Biased towards the sympathetic side: sympathetic nerves are strong and parasympathetic nerves are weak;

[0042] d. Biased towards the parasympathetic nerve side: the sympathetic nerves are normal and the parasympathetic nerves are strong;

[0043] e. Biased towards the parasympathetic side: sympathetic nerves are weak and parasympathetic nerves are normal;

[0044] f. Biased towards the parasympathetic nerve side: the sympathetic nerves are weak and the parasympathetic nerves are strong.

[0045] A portable method for evaluating autonomic nervous function comprises the following steps:

[0046] Collect physiological signals and obtain the user's heart rate fluctuation through preprocessing operations;

[0047] Adaptively measure the user's breathing ability, determine a personalized breathing rhythm, and guide the user to follow that breathing rhythm;

[0048] Perform spectrum analysis and causal analysis on the user's heart rate signal during breathing training, and calculate the cardiopulmonary coupling index, including the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during long inhalation and short exhalation, and the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during short inhalation and long exhalation;

[0049] Physiological characteristics are integrated and the user's autonomic nervous system state is determined through an autonomic nervous balance analysis model and an autonomic nervous activity analysis model; the autonomic nervous balance is characterized by the ratio of the cardiopulmonary follow-up response index in different processes, and the autonomic nervous activity is characterized by the cardiopulmonary coupling strength.

[0050] An electronic device comprises a memory, a processor and a computer program / instruction stored in the memory and executable on the processor, wherein the steps of the portable autonomic nervous function assessment method are implemented when the computer program / instruction is executed by the processor.

[0051] Furthermore, the electronic device also includes a physiological signal acquisition device, or acquires the collected physiological signals or the user's heart rate signal by communicating with the physiological signal acquisition device.

[0052] A computer program product comprises a computer program / instruction, wherein when the computer program / instruction is executed by a processor, the steps of the portable assessment method of autonomic nervous function are implemented.

[0053] Beneficial effects: Compared with the prior art, the present invention has the following beneficial effects:

[0054] As an improvement and innovation, the present invention (1) proposes a portable evaluation device for autonomic nervous function based on the principle of cardiopulmonary coupling, which only needs to extract and analyze the heart rate signal of the user in the dynamic process, has the characteristics of simplicity and lightness, and can be transplanted to various types of wearable health devices or other terminals with heart rate monitoring functions that are currently widely used; (2) proposes an intelligent breathing test and guidance method, which adaptively determines the user's breathing ability through the heart rate signal, sets a personalized breathing frequency, and guides the user to grasp the breathing rhythm through multiple sensory media such as vision, hearing, and touch; (3) proposes a cardiopulmonary coupling analysis framework, which combines spectrum and causal analysis to accurately and quantitatively extract the user's cardiopulmonary coupling characteristics, and realize objective and detailed evaluation of the user's autonomic nervous system function.

[0055] The beneficial effects of the present invention are: based on the accurate decoding of the cardiopulmonary coupling process, a fast, low-cost, easy-to-operate, and popularizable portable autonomic nervous function assessment device, method, equipment, and program product are realized. The required measurement device is simple and can be used for individualized home health assessment. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] Figure 1 It is a schematic diagram of the technical solution of the present invention.

[0057] Figure 2 It is a schematic diagram of the system composition of the present invention.

[0058] Figure 3 The figure is a schematic diagram of a specific implementation process of a certain embodiment of the present invention. DETAILED DESCRIPTION

[0059] The technical details of the technical solution of the present invention are further described below in conjunction with the accompanying drawings and specific embodiments.

[0060] Rhythmic breathing can induce obvious cardiopulmonary coupling, that is, a high degree of follow-up response between breathing and heart rate, which is closely related to the autonomic nervous system. When the autonomic nervous system loses balance or is damaged, the cardiopulmonary coupling phenomenon related to the autonomic nervous function will also weaken. Therefore, the present invention proposes a new idea, that is, measuring the user's heart rate under a specific breathing rhythm set adaptively, combining a feature extraction algorithm to quantitatively construct a correlation between the cardiopulmonary coupling intensity and the autonomic nervous function, in the hope of realizing a fast, accurate, low-cost, easy to operate, and popularizable autonomic nervous function assessment method. Related embodiments are as follows:

[0061] Example 1: Figure 1As shown, an autonomic nervous function portable evaluation device disclosed in an embodiment of the present invention mainly includes: a heart rate extraction and preprocessing module, an intelligent breathing test and guidance module, a cardiopulmonary coupling feature extraction module, an autonomic nervous function analysis module, etc. The heart rate extraction and preprocessing module is used to collect physiological signals and obtain the user's heart rate fluctuation through preprocessing operations; the intelligent breathing test and guidance module is used to adaptively measure the user's breathing ability, determine the personalized breathing rhythm, and guide the user to follow the breathing rhythm; the cardiopulmonary coupling feature extraction module is used to perform spectrum analysis and causal analysis on the user's heart rate signal during breathing training, and calculate the cardiopulmonary coupling index, including the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during long inhalation and short exhalation, and the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during short inhalation and long exhalation; the autonomic nervous function analysis module is used to integrate physiological characteristics and determine the user's autonomic nervous system state through the autonomic nervous balance analysis model and the autonomic nervous activity analysis model; wherein the autonomic nervous balance is characterized by the ratio of the cardiopulmonary follow-up response index in different processes, and the autonomic nervous activity is characterized by the cardiopulmonary coupling strength. The central controller is used to control each module.

[0062] The specific processing idea of ​​the heart rate extraction and preprocessing module can be: collect original physiological signals containing heart rate information such as PPG or ECG through the device, use filtering processing to remove baseline and high-frequency noise, and calculate the heart rate information in the detection process based on the peak feature points in the original physiological signal.

[0063] Specifically, in some embodiments, the preprocessing operation steps in the heart rate extraction and preprocessing module specifically include:

[0064] Find the peak point in the physiological signal and obtain the peak point sequence {r(n), n=1,2,…,N} derived from the physiological signal, where N is the total number of points in the sequence;

[0065] Perform a first-order difference operation on the peak point sequence to calculate the peak point interval sequence {rr(n), i=1, 2, ..., N-1};

[0066] The user's beat-to-beat heart rate is calculated based on the peak point interval sequence:

[0067] The specific processing idea of ​​the intelligent breathing test and guidance module can be: analyze the heart rate signal of the user during the uniform deep breathing process, calculate the comprehensive index for evaluating the user's breathing ability, set the appropriate breathing rhythm based on the index, and guide the user to follow the breathing rhythm based on sensory media such as vision, hearing, and touch.

[0068] Specifically, in some embodiments, the steps of measuring breathing capacity and formulating and guiding breathing rhythm in the intelligent breathing test and guidance module specifically include:

[0069] Introducing the designated breathing method and rhythm through voice and / or video, guiding the user to continue to inhale evenly until no more gas can be inhaled, and then continue to exhale evenly until no more gas can be exhaled;

[0070] The first-order difference of the user's heart rate during the calculation process is: hr′(n)=hr(n)-hr(n-1), n=2,...,N-1, where N is the total number of peak points in the physiological signal sequence, hr i The user's beat-to-beat heart rate;

[0071] The first point in the sequence where hr′(n) ≥ 2bmp is defined as the starting inhalation key frame n. 1 ; The first point after the inhalation key frame in time sequence and with hr′(n)≤0.5bmp is defined as the transition key frame n 2 ; The first point after the transition key frame in time sequence and hr′(n)≤-2bmp is defined as the starting exhalation key frame n 3 ; The first point after the start exhalation key frame in time sequence and with hr′(n)≥-0.5bmp is defined as the end key frame n 4 ;

[0072] Calculate the user's maximum constant inspiratory time: And maximum constant exhalation time: in to is the time corresponding to the corresponding key frame;

[0073] The user's breathing rhythm is set in a personalized and adaptive manner based on the measured maximum uniform inhalation time and maximum uniform exhalation time of the user, and the user is guided to breathe according to the set rhythm through voice, vibration and / or video, including long inhalation and short exhalation training and short inhalation and long exhalation training.

[0074] The specific processing idea of ​​the cardiopulmonary coupling feature extraction module can be: based on methods such as spectrum analysis and causal analysis, the cardiopulmonary coupling features related to the autonomic nervous system are extracted.

[0075] Specifically, in some embodiments, the calculation steps of the cardiopulmonary following response index and the cardiopulmonary coupling strength in the cardiopulmonary coupling feature extraction module specifically include:

[0076] The heart rate fluctuation is converted from the time domain to the frequency domain by discrete Fourier transform:

[0077]

[0078] Where i is an imaginary unit, N is the number of sampling points, hr(n) is the heart rate signal, and HR[k] is the frequency domain distribution of the heart rate signal;

[0079] Calculate the power of each point in the discrete frequency domain to obtain the power spectrum density of the heart rate signal: where f s is the sampling frequency of the original physiological signal.

[0080] Calculate the cardiopulmonary response index during long inhalation and short exhalation: where f 1 It is the fixed respiratory rate during long inhalation and short exhalation, peak s1 Indicates the maximum frequency of the peak point of the heart rate signal spectrum, C 1 yes The integer part of ;

[0081] Calculate the cardiopulmonary response index during short inhalation and long exhalation: where f 2 It is the fixed breathing rate during short inhalation and long exhalation, peak s2 Indicates the maximum frequency of the peak point of the heart rate signal spectrum, C 2 yes The integer part of ;

[0082] The heart rate signal hr(n) is decomposed into multiple physically meaningful components HR through EEMD j (n), j = 1, 2, ..., m;

[0083] The respiratory timing signal is transmitted through EEMD Disassemble into multiple components with physical meaning RESP j (n), j = 1, 2, ..., m;

[0084] HR is obtained by Hilbert transform j The phase φ of (n) x (n), and RESP j The phase φ of (n) y (n);

[0085] Calculate HR j (n) and RESP j (n) coupling strength: Where i is the imaginary unit;

[0086] The cardiopulmonary coupling strength of the user is calculated as: CCS = max{Coh 1 , Coh 2 , ..., Coh m}.

[0087] The specific processing idea of ​​the autonomic nervous function analysis module can be: based on the statistical analysis results, a multidimensional comprehensive model is constructed to clarify the mapping relationship between the cardiopulmonary coupling characteristics and the autonomic nervous state, and to achieve quantitative analysis of physiological indicators including but not limited to autonomic nervous balance, sympathetic nervous activity, parasympathetic nervous activity, etc.

[0088] Specifically, in some embodiments, the step of determining the state of the user's autonomic nervous system in the autonomic nervous function analysis module specifically includes:

[0089] Calculates the user's autonomic balance: CRI 1 CRI is the cardiopulmonary response index during long inhalation and short exhalation. 2 It is the cardiopulmonary response index during short inhalation and long exhalation;

[0090] Based on the statistical analysis results, determine whether the user's autonomic nervous system state deviates from the normal range;

[0091] If the user's autonomic nervous system state deviates from the normal range, the autonomic nervous system activity is calculated: CCS 1 CCS is the cardiopulmonary coupling strength during long inhalation and short exhalation. 2 It is the cardiopulmonary coupling strength during short inhalation and long exhalation; SNA represents the activity of sympathetic nerves, and PNA represents the activity of parasympathetic nerves;

[0092] A feature vector [SNA, SNA] is constructed, and based on the feature vector, the distance between the user's autonomic nervous activity and the respective autonomic nervous dysfunction types obtained by clustering is calculated, so as to determine the specific type of the user's autonomic nervous dysfunction.

[0093] Example 2: Figure 2 As shown, a portable method for evaluating autonomic nervous function disclosed in an embodiment of the present invention comprises the following steps:

[0094] (1) Heart rate signal extraction: Collect physiological signals and obtain the user's heart rate fluctuation through preprocessing operations;

[0095] (2) Intelligent breathing test and guidance: Adaptively measure the user's breathing ability, determine the personalized breathing rhythm, and guide the user to follow the breathing rhythm;

[0096] (3) Extraction of cardiopulmonary coupling features: Spectral analysis and causal analysis are performed on the user's heart rate signal during breathing training to calculate the cardiopulmonary coupling index, including the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during long inhalation and short exhalation, and the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during short inhalation and long exhalation;

[0097] (4) Construction of an autonomic nervous function analysis model: Integrating physiological characteristics, and determining the user's autonomic nervous system state through an autonomic nervous balance analysis model and an autonomic nervous activity analysis model; the autonomic nervous balance is represented by the ratio of the cardiopulmonary follow-up response index in different processes, and the autonomic nervous activity is represented by the cardiopulmonary coupling strength.

[0098] The executor of the heart rate signal extraction method in step (1) may be a terminal such as a wearable device, including but not limited to any terminal or system capable of collecting physiological signals such as PPG signals, ECG signals, millimeter wave radar, etc. that can be used to calculate heart rate information.

[0099] This embodiment can record the original physiological signal containing heart rate information through a physiological signal acquisition device integrated with a PPG photoelectric sensor or an integrated ECG sensor. The PPG photoelectric sensor can include visible light such as green light and red light, and can also include invisible light of multiple wavelengths such as infrared light. The ECG sensor can include a single lead or multiple lead electrodes.

[0100] The influence of baseline drift, power frequency interference and other noises in the original physiological signal on heart rate extraction can be reduced by filtering processing. The filtering processing can be bandpass filtering. The filtering can use a zero-phase shift digital filter to bandpass filter the signal. The signal can be processed by bandpass filtering to obtain a physiological signal with baseline drift and high-frequency noise removed.

[0101] The findpeaks function in MATLAB software can be used to find the peak points in the physiological signal and obtain the peak point sequence {r(n), n=1,2,...,N} derived from the physiological signal.

[0102] The peak point interval sequence {rr(n), i=1, 2, ..., N-1} can be calculated by performing a first-order difference operation on the peak point sequence.

[0103] The user's beat-to-beat heart rate can be calculated based on the peak point interval sequence:

[0104] The execution subject of the implementation method of the intelligent breathing test and guidance program in step (2) can be a mobile phone, tablet computer, wearable device and other devices, or a laptop computer, desktop computer and other devices with interaction, storage and data processing capabilities, or other intelligent terminal devices.

[0105] In this embodiment, the method in step (2) can be used to execute a physical key or a virtual key to trigger the intelligent respiratory monitoring and guidance program. The program can include the following modules: a test guidance module, a collection module, a test analysis module, and a training guidance module.

[0106] The test guidance module can introduce the specified breathing method and rhythm through voice (or video) and other means. The breathing method and rhythm specified in this module can be: continue to inhale evenly until no more gas can be inhaled, and then continue to exhale evenly until no more gas can be exhaled.

[0107] The function of the acquisition module can be to capture the user's heart rate fluctuations during the breathing monitoring process. The module can be implemented based on the heart rate extraction method in step (1).

[0108] The function of the test analysis module can be: based on the principle of cardiopulmonary coupling, convert the heart rate fluctuations captured by the acquisition module into breathing-related features. Breathing-related features can include maximum uniform inhalation time and maximum uniform exhalation time. The analysis module can perform a first-order difference operation on the user's heart rate fluctuations during the test to obtain the heart rate change rate: hr′(n)=hr(n)-hr(n-1), n=2, ..., N-1.

[0109] The first point in the timing where hr′(n)≥2bmp can be defined as the starting inhalation key frame n 1 ; The first point after the inhalation key frame in time sequence and with hr′(n)≤0.5bmp is defined as the transition key frame n 2 ; The first point after the transition key frame in time sequence and hr′(n)≤-2bmp is defined as the starting exhalation key frame n 3 ; The first point after the start exhalation key frame in time sequence and with hr′(n)≥-0.5bmp is defined as the end key frame n 4 .

[0110] The measured maximum uniform inspiratory time can be expressed as:

[0111] The maximum uniform exhalation time measured can be expressed as:

[0112] The function of the training guidance module can be: according to the maximum uniform inhalation time and maximum uniform exhalation time of the user measured by the test analysis module, the appropriate breathing rhythm is set in a personalized and adaptive manner, and the user is guided to breathe according to the set rhythm through voice, vibration, video, etc. In order to avoid the influence of individual breathing habit differences on the subsequent analysis results, the breathing training of the training guidance module can be: 3 minutes of long inhalation and short exhalation training and 3 minutes of short inhalation and long exhalation training.

[0113] Among them, the inhalation time in long inhalation and short exhalation training can be set to (0.75×Δt 1 ) seconds, the exhalation time can be set to (0.55×Δt 2 ) seconds; the inhalation time in short inhalation and long exhalation training can be set to (0.55×Δt1 ) seconds, the exhalation time can be set to (0.75×Δt 2 )Second.

[0114] The execution subject of the cardiopulmonary coupling feature extraction method in step (3) can be a device such as a mobile phone, a tablet computer, a wearable device, or a device with data processing capabilities such as a laptop computer and a desktop computer.

[0115] This embodiment can extract the user's cardiopulmonary coupling characteristics through power spectrum analysis and causal analysis.

[0116] Power spectrum analysis can convert heart rate fluctuations from the time domain to the frequency domain through discrete Fourier transform:

[0117] The power spectrum density of the heart rate signal can be obtained by calculating the power of each point in the discrete frequency domain:

[0118] Among them, the power spectrum density corresponding to the i-th peak point of the heart rate signal power spectrum can be expressed as: PSD (peak i ).

[0119] For the long inhalation and short exhalation process, the fixed respiratory frequency can be expressed as:

[0120] From the heart rate-respiration follow-up response, we can know that the cardiopulmonary follow-up response index during long inhalation and short exhalation can be expressed as:

[0121]

[0122] For the short inhalation and long exhalation process, the fixed respiratory frequency can be expressed as:

[0123] From the heart rate-breathing follow-up response, we can know that the cardiopulmonary follow-up response index during short inhalation and long exhalation can be expressed as:

[0124]

[0125] Causal analysis can decompose the heart rate signal into multiple physically meaningful components HR through EEMD j (n), j = 1, 2, ..., m.

[0126] According to the set breathing rhythm, the state value of the inhalation process can be set to "1", and the state value of the exhalation process can be set to "0" to obtain the timing signal of the specified breathing rhythm Similarly, the respiratory timing signal can be decomposed into multiple physically meaningful components RESP by EEMD. j(n), j = 1, 2, ..., m.

[0127] HR can be solved by Hilbert transform j The phase φ of (t) x (n), and RESP j The phase φ of (t) y (n).

[0128] HR j (t) and RESP j The coupling strength between (t) can be expressed as: where i is the imaginary unit.

[0129] From the heart rate-breathing follow-up response, it can be seen that the cardiopulmonary coupling strength during breathing training can be expressed as: CCS = max{Coh 1 , Coh 2 , ..., Coh m}.

[0130] Among them, the cardiopulmonary coupling strength during long inhalation and short exhalation can be recorded as CCS 1 The cardiopulmonary coupling strength during short inhalation and long exhalation can be recorded as CCS 2 .

[0131] The execution subject of the method for constructing the autonomic nervous function analysis model in step (4) can be a device such as a mobile phone, a tablet computer, a wearable device, or a device with data processing capabilities such as a laptop computer and a desktop computer.

[0132] The autonomic nervous balance analysis model can be constructed by integrating the cardiopulmonary response index in different processes:

[0133] The normal range of ANB can be selected based on the statistical distribution of large sample data, where the high value can be selected as ANB ≥ 1.25 and the low value can be selected as ANB ≤ 0.8. A higher ANB indicates that the autonomic balance is biased towards the sympathetic side, while a lower ANB indicates that the autonomic balance is biased towards the parasympathetic side.

[0134] After preliminary judgment through power spectrum analysis, the autonomic nervous activity analysis model can be constructed by integrating the cardiopulmonary follow-up response index of individuals with autonomic nervous imbalance in different processes:

[0135] Among them, SNA represents the activity of sympathetic nerves, and PNA represents the activity of parasympathetic nerves.

[0136] In order to achieve a more comprehensive representation of the autonomic nervous system and provide targeted personalized suggestions, it is necessary to classify individuals with different autonomic nervous system imbalances. i, its autonomic nervous function characteristics can be represented by the feature vector ANS(U i )=[SNA i , PNA i ] indicates that the set {ANS(U 1 ), ANS(U 2 ), ..., ANS(U N )} for clustering, the number of clusters k can be selected as 6. Cluster analysis obtains the set {K n (x n ,y n ), n = 1, 2, ..., 6}, corresponding to the following 6 types of autonomic imbalance:

[0137] a. Biased towards the sympathetic side: sympathetic nerves are strong and parasympathetic nerves are normal;

[0138] b. Biased towards the sympathetic side: the sympathetic nerves are normal, and the parasympathetic nerves are weak;

[0139] c. Biased towards the sympathetic side: sympathetic nerves are strong and parasympathetic nerves are weak;

[0140] d. Biased towards the parasympathetic nerve side: the sympathetic nerves are normal and the parasympathetic nerves are strong;

[0141] e. Biased towards the parasympathetic side: sympathetic nerves are weak and parasympathetic nerves are normal;

[0142] f. Biased towards the parasympathetic side: sympathetic nerves are weak and parasympathetic nerves are strong;

[0143] For U i , we can calculate the Euclidean distance from each set Then D k ={D n} min , and its corresponding index k is the user's autonomic nervous system imbalance category. Based on the classification result, a personalized adjustment plan can be automatically selected from the constructed knowledge graph and pushed to the specified user. The specific analysis result is shown in the following example. Figure 3 .

[0144] As shown in the above implementation process, this method can first determine the autonomic balance based on spectrum analysis, and then solve the autonomic activity of individuals with autonomic dysfunction / disorder through causal analysis, and determine the type of disorder. Among them, the power spectrum analysis principle is simple, the operation is fast, and the processor and storage resources occupied are relatively small; while the causal analysis is more complicated to calculate due to the introduction of EEMD and Hilbert transform, and occupies more processor and storage resources. Therefore, the diversion step processing strategy introduced in this embodiment, that is, firstly make a rough judgment through power spectrum analysis with less computing power consumption. If the user's autonomic imbalance is found, more computing power is further called to solve the type of imbalance, which can save onboard computing power to a certain extent.

[0145] Embodiment 3: An electronic device disclosed in an embodiment of the present invention includes a memory, a processor, and a computer program / instruction stored in the memory and executable on the processor, wherein the steps of the portable evaluation method for autonomic nervous function are implemented when the computer program / instruction is executed by the processor. The electronic device may be a mobile phone, a tablet computer, a wearable device, etc., including a physiological signal acquisition device; or a mobile phone, a tablet computer, a laptop computer, a desktop computer, etc., which acquires collected physiological signals or processed heart rate signals by communicating with the physiological signal acquisition device.

[0146] Embodiment 4: A computer program product disclosed in an embodiment of the present invention comprises a computer program / instruction, which implements the steps of the portable method for assessing autonomic nervous function when executed by a processor. The program / instruction code for implementing the method of the present invention can be written in any combination of one or more programming languages. These program / instruction codes can be provided to a processor or controller of a general-purpose computer, a special-purpose computer or other programmable data processing device, so that the steps of the method of the present invention are implemented when the program / instruction code is executed by the processor or controller. The program / instruction code can be executed entirely on the machine, partially on the machine, partially on the machine as an independent software package and partially on a remote machine, or entirely on a remote machine or server. Anything not described in detail in the present invention is a well-known technology for those skilled in the art.

Claims

1. A portable device for evaluating autonomic nervous function, characterized in that: The device comprises: The heart rate extraction and preprocessing module is used to collect physiological signals and obtain the user's heart rate fluctuation through preprocessing operations; Intelligent breathing test and guidance module, used to adaptively measure the user's breathing ability, determine the personalized breathing rhythm, and guide the user to follow the breathing rhythm; The cardiopulmonary coupling feature extraction module is used to perform spectrum analysis and causal analysis on the user's heart rate signal during breathing training, and calculate the cardiopulmonary coupling index, including the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during long inhalation and short exhalation, and the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during short inhalation and long exhalation; The autonomic nervous function analysis module is used to integrate physiological characteristics and determine the user's autonomic nervous system state through an autonomic nervous balance analysis model and an autonomic nervous activity analysis model; the autonomic nervous balance is characterized by the ratio of the cardiopulmonary follow-up response index in different processes, and the autonomic nervous activity is characterized by the cardiopulmonary coupling strength.

2. The portable device for evaluating autonomic nervous function according to claim 1, characterized in that: In the heart rate extraction and preprocessing module, the preprocessing operation steps specifically include: Find the peak point in the physiological signal and obtain the peak point sequence {r(n), n=1, 2, ..., N} derived from the physiological signal, where N is the total number of points in the sequence; Perform a first-order difference operation on the peak point sequence to calculate the peak point interval sequence {rr(n), i=1, 2, ..., N-1}; The user's beat-to-beat heart rate is calculated based on the peak point interval sequence:

3. The portable device for evaluating autonomic nervous function according to claim 1, characterized in that: In the intelligent breathing test and guidance module, the steps of measuring breathing capacity and formulating and guiding breathing rhythm specifically include: Introducing the designated breathing method and rhythm through voice and / or video, guiding the user to continue to inhale evenly until no more gas can be inhaled, and then continue to exhale evenly until no more gas can be exhaled; The first-order difference of the user's heart rate during the calculation process: hr′(n)=hr(n)-hr(n-1), n=2,...,N-1, where N is the total number of peak points in the physiological signal sequence, and hr(n) is the user's beat-to-beat heart rate; The first point in the sequence where hr′(n)≥2bmp is defined as the starting inhalation key frame n1; the first point in the sequence where hr′(n)≤0.5bmp is located after the inhalation key frame is defined as the transition key frame n2; the first point in the sequence where hr′(n)≤-2bmp is located after the transition key frame is defined as the starting exhalation key frame n3; the first point in the sequence where hr′(n)≥-0.5bmp is located after the starting exhalation key frame is defined as the ending key frame n4; Calculate the user's maximum constant inspiratory time: And maximum constant exhalation time: in to is the time corresponding to the corresponding key frame; The user's breathing rhythm is set in a personalized and adaptive manner based on the measured maximum uniform inhalation time and maximum uniform exhalation time of the user, and the user is guided to breathe according to the set rhythm through voice, vibration and / or video, including long inhalation and short exhalation training and short inhalation and long exhalation training.

4. The portable device for evaluating autonomic nervous function according to claim 1, characterized in that: In the cardiopulmonary coupling feature extraction module, the steps of calculating the cardiopulmonary follow-up response index and the cardiopulmonary coupling strength specifically include: The heart rate fluctuation is converted from the time domain to the frequency domain by discrete Fourier transform: Where i is an imaginary unit, N is the number of sampling points, hr(n) is the heart rate signal, and HR[k] is the frequency domain distribution of the heart rate signal; Calculate the power of each point in the discrete frequency domain to obtain the power spectrum density of the heart rate signal: where f s is the sampling frequency of the original physiological signal; Calculate the cardiopulmonary response index during long inhalation and short exhalation: Where f1 is the fixed respiratory frequency during long inhalation and short exhalation, peak s1 Indicates the maximum frequency of the peak point of the heart rate signal spectrum, C1 is The integer part of ; Calculate the cardiopulmonary response index during short inhalation and long exhalation: Where f2 is the fixed breathing frequency during short inhalation and long exhalation, peak s2 Indicates the maximum frequency of the peak point of the heart rate signal spectrum, C2 is The integer part of ; The heart rate signal hr(n) is decomposed into multiple physically meaningful components HR through EEMD j (n),j=1,2,…,m; The respiratory timing signal is transmitted through EEMD Disassemble into multiple components with physical meaning RESP j (n),j=1,2,…,m; HR is obtained by Hilbert transform j The phase φ of (n) x (n), and RESP j The phase φ of (n) y (n); Calculate HR j (n) and RESP j (n) coupling strength: Where i is the imaginary unit; The cardiopulmonary coupling strength of the user is calculated as: CCS = max{Coh1,Coh2,…,Coh m }.

5. The portable device for evaluating autonomic nervous function according to claim 1, characterized in that: In the autonomic nervous function analysis module, the step of determining the state of the user's autonomic nervous system specifically includes: Calculates the user's autonomic balance: Among them, CRI1 is the cardiopulmonary response index during long inhalation and short exhalation, and CRI2 is the cardiopulmonary response index during short inhalation and long exhalation; Based on the statistical analysis results, determine whether the user's autonomic nervous system state deviates from the normal range; If the user's autonomic nervous system state deviates from the normal range, the autonomic nervous system activity is calculated: Among them, CCS1 is the cardiopulmonary coupling strength during long inhalation and short exhalation, and CCS2 is the cardiopulmonary coupling strength during short inhalation and long exhalation; SNA represents the activity of sympathetic nerves, and PNA represents the activity of parasympathetic nerves; A feature vector [SNA, SNA] is constructed, and based on the feature vector, the distance between the user's autonomic nervous activity and the respective autonomic nervous dysfunction types obtained by clustering is calculated, so as to determine the specific type of the user's autonomic nervous dysfunction.

6. The portable device for evaluating autonomic nervous function according to claim 5, characterized in that: The autonomic dysfunction type is obtained by performing K-means clustering on the collected sample data. The number of clusters is set to 6, corresponding to 6 types of autonomic imbalance: a. Biased towards the sympathetic side: sympathetic nerves are strong and parasympathetic nerves are normal; b. Biased towards the sympathetic side: the sympathetic nerves are normal, and the parasympathetic nerves are weak; c. Biased towards the sympathetic side: sympathetic nerves are strong and parasympathetic nerves are weak; d. Biased towards the parasympathetic nerve side: the sympathetic nerves are normal and the parasympathetic nerves are strong; e. Biased towards the parasympathetic side: sympathetic nerves are weak and parasympathetic nerves are normal; f. Biased towards the parasympathetic nerve side: the sympathetic nerves are weak and the parasympathetic nerves are strong.

7. A portable method for evaluating autonomic nervous function, characterized in that: The steps include: Collect physiological signals and obtain the user's heart rate fluctuation through preprocessing operations; Adaptively measure the user's breathing ability, determine a personalized breathing rhythm, and guide the user to follow that breathing rhythm; Perform spectrum analysis and causal analysis on the user's heart rate signal during breathing training, and calculate the cardiopulmonary coupling index, including the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during long inhalation and short exhalation, and the cardiopulmonary follow-up response index and cardiopulmonary coupling strength during short inhalation and long exhalation; Physiological characteristics are integrated and the user's autonomic nervous system state is determined through an autonomic nervous balance analysis model and an autonomic nervous activity analysis model; the autonomic nervous balance is characterized by the ratio of the cardiopulmonary follow-up response index in different processes, and the autonomic nervous activity is characterized by the cardiopulmonary coupling strength.

8. An electronic device comprising a memory, a processor, and a computer program / instruction stored in the memory and executable on the processor, characterized in that: When the computer program / instructions are executed by a processor, the steps of the method according to claim 7 are implemented.

9. The electronic device according to claim 8, characterized in that: It also includes a physiological signal acquisition device, or acquires the collected physiological signals or the user's heart rate signal by communicating with the physiological signal acquisition device.

10. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instructions are executed by a processor, the steps of the method of claim 7 are implemented.

Images