Low-stimulation electric thumbtack needle nerve electrical stimulation parameter optimization method and system

Abstract

The invention relates to the technical field of electrical thumbtack needle nerve electrical stimulation, in particular to a low-stimulation electrical thumbtack needle nerve electrical stimulation parameter optimization method and system, and the method comprises the steps: obtaining first data; the electric thumbtack needle is inserted into the acupuncture point of the stimulation object, the actual stimulation frequency of the electric thumbtack needle is set as the reference stimulation frequency, the actual pulse width is set as the reference pulse width, and the actual current intensity is set as the reference current intensity; reading a reference blood pressure value, a reference heart rate value and a reference electromyographic signal amplitude; the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric thumbtack needle are dynamically set with the preset adjustment interval time as the interval, and a blood pressure value three-dimensional function, a heart rate value three-dimensional function and an electromyographic signal amplitude three-dimensional function are obtained; and electric thumbtack needle nerve electrical stimulation optimization parameters are calculated through a parameter optimization algorithm. The electric thumbtack needle nerve electrical stimulation optimization parameters are calculated through a parameter optimization algorithm, and therefore it is ensured that the stimulated object maintains the optimal physiological indexes in the nerve electrical stimulation process.

Description

Technical Field

[0001] The present invention relates to the technical field of electric acupuncture nerve electrical stimulation, and in particular to a method and system for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation. Background Art

[0002] Low-stimulation electroacupuncture is a novel medical technique that combines traditional acupuncture with electrical stimulation therapy, offering advantages such as ease of use and effective treatment. However, in existing electroacupuncture treatments, the setting of nerve stimulation parameters often relies on the physician's experience and lacks scientific optimization methods. This can easily lead to a mismatch between stimulation intensity and the patient's physiological parameters, resulting in rapid changes in blood pressure, heart rate, and muscle spasms. Accurately setting electroacupuncture nerve stimulation parameters based on individual patient differences remains a pressing technical challenge in this field. To address this issue, existing techniques primarily divide the range of nerve stimulation parameter values ​​into several segments and then manually select appropriate stimulation parameters within these segments based on the intended purpose of the stimulation. However, this approach still suffers from the following drawbacks: First, parameter setting lacks personalization, considering only the therapeutic objective while failing to fully account for patient-to-patient differences in sensitivity to electrical stimulation; second, parameter adjustment lacks real-time feedback, preventing dynamic optimization based on the patient's physiological responses; and third, parameter selection lacks a systematic approach, making it difficult to ensure optimal physiological indicators during stimulation.

[0003] Therefore, it is necessary to adaptively optimize the neural electrical stimulation parameters based on the real-time measurement data of the stimulated subject's physiological parameters, so as to ensure that the stimulated subject maintains the optimal physiological indicators during the neural electrical stimulation process. Summary of the Invention

[0004] (1) Technical problems to be solved The purpose of the present invention is to provide a method and system for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation, so as to realize adaptive optimization of parameters of low-stimulation electric acupuncture nerve electrical stimulation.

[0005] (2) Technical solution To achieve the above objectives, the present invention provides a method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation, the method comprising the following steps: S1. Obtain first data, which includes pre-set neural electrical stimulation baseline parameters, neural electrical stimulation lower limit parameters, and neural electrical stimulation upper limit parameters; the neural electrical stimulation baseline parameters include baseline stimulation frequency, baseline pulse width, and baseline current intensity; the neural electrical stimulation lower limit parameters include lower limit stimulation frequency, lower limit pulse width, and lower limit current intensity; the neural electrical stimulation upper limit parameters include upper limit stimulation frequency, upper limit pulse width, and upper limit current intensity.

[0006] S2, insert the electric needle into the acupuncture point of the stimulation object, set the actual stimulation frequency of the electric needle to the reference stimulation frequency, the actual pulse width to the reference pulse width, and the actual current intensity to the reference current intensity; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation object through the physiological parameter acquisition device to obtain the reference blood pressure value, reference heart rate value, and reference electromyographic signal amplitude.

[0007] S3, dynamically setting the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle with the preset adjustment interval time as the interval, and reading the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object through the physiological parameter acquisition device to obtain a three-dimensional function of the blood pressure value, a three-dimensional function of the heart rate value, and a three-dimensional function of the electromyographic signal amplitude.

[0008] S4. Based on the three-dimensional function of blood pressure value, heart rate value and electromyographic signal amplitude, the optimized parameters of electric acupuncture nerve stimulation are calculated by parameter optimization algorithm. The optimized parameters of electric acupuncture nerve stimulation include optimal stimulation frequency, optimal pulse width and optimal current intensity.

[0009] Furthermore, the method of dynamically setting the actual stimulation frequency, actual pulse width, and actual current intensity of the electric acupuncture needle at a preset adjustment interval time, and obtaining a three-dimensional function of blood pressure value, heart rate value, and electromyographic signal amplitude by reading the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject through a physiological parameter acquisition device includes: The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the lower limit stimulation frequency, reference pulse width, and reference current intensity, respectively; the actual stimulation frequency is adjusted with a preset adjustment interval time as the interval and a preset first frequency as the step size, so that the actual stimulation frequency increases stepwise until it is greater than the difference between the upper limit stimulation frequency and the first frequency for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object are read to obtain a first blood pressure array, a first heart rate array, and a first electromyographic signal amplitude array.

[0010] The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the reference stimulation frequency, lower limit pulse width, and reference current intensity, respectively; the actual pulse width is adjusted with a preset adjustment interval time as the interval and a preset first pulse width as the step size, so that the actual pulse width increases stepwise until it is greater than the difference between the upper limit pulse width and the first pulse width for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object are read to obtain a second blood pressure array, a second heart rate array, and a second electromyographic signal amplitude array.

[0011] The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the reference stimulation frequency, reference pulse width, and lower limit current intensity, respectively; the actual current intensity is adjusted with a preset adjustment interval time as the interval and a preset first current intensity as the step size, so that the actual current intensity increases stepwise until it is greater than the difference between the upper limit current intensity and the first current intensity for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object are read to obtain a third blood pressure array, a third heart rate array, and a third electromyographic signal amplitude array.

[0012] According to the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array, a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude are obtained through a fitting algorithm.

[0013] Furthermore, the method of obtaining a three-dimensional function of blood pressure values, a three-dimensional function of heart rate values, and a three-dimensional function of electromyographic signal amplitudes by a fitting algorithm based on the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array includes: According to the first blood pressure array, a nonlinear least squares fitting algorithm is used to fit a first blood pressure function; the first blood pressure function is expressed as ;in represents the actual stimulation frequency; according to the first heart rate array, a nonlinear least squares fitting algorithm is used to fit the first heart rate function; the first heart rate function is expressed as ; According to the first electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the first electromyographic signal amplitude function; the first electromyographic signal amplitude function is expressed as .

[0014] According to the second blood pressure array, a nonlinear least squares fitting algorithm is used to fit the second blood pressure function; the second blood pressure function is expressed as ;in represents the actual pulse width; according to the second heart rate array, a nonlinear least squares fitting algorithm is used to fit the second heart rate function; the second heart rate function is expressed as ; According to the second electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the second electromyographic signal amplitude function; the second electromyographic signal amplitude function is expressed as .

[0015] According to the third blood pressure array, a nonlinear least squares fitting algorithm is used to fit a third blood pressure function; the third blood pressure function is expressed as ;in Represents the actual current intensity; According to the third heart rate array, a nonlinear least squares fitting algorithm is used to fit the third heart rate function; The third heart rate function is expressed as ; According to the third myoelectric signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the third myoelectric signal amplitude function; the third myoelectric signal amplitude function is expressed as .

[0016] According to the first blood pressure function, the second blood pressure function, the third blood pressure function, the first heart rate function, the second heart rate function, the third heart rate function, the first electromyographic signal amplitude function, the second electromyographic signal amplitude function, and the third electromyographic signal amplitude function, a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude are calculated by a three-dimensional construction formula; the three-dimensional construction formula is: ; in, represents the three-dimensional function of blood pressure value, represents the three-dimensional function of heart rate value, represents the three-dimensional function of the EMG signal amplitude, represents the base stimulation frequency, Indicates the reference pulse width, Indicates the reference current intensity, Indicates the baseline blood pressure value, Indicates the baseline heart rate value. Indicates the baseline EMG signal amplitude.

[0017] Furthermore, the method of calculating the optimized parameters of the electroacupuncture nerve electrical stimulation by a parameter optimization algorithm based on the three-dimensional function of the blood pressure value, the three-dimensional function of the heart rate value, and the three-dimensional function of the electromyographic signal amplitude, wherein the optimized parameters of the electroacupuncture nerve electrical stimulation include the optimal stimulation frequency, the optimal pulse width, and the optimal current intensity, includes: The first objective function is constructed according to the three-dimensional function of blood pressure value, the three-dimensional function of heart rate value and the three-dimensional function of electromyographic signal amplitude.

[0018] A first constraint condition is constructed based on the lower limit stimulation frequency, the lower limit pulse width, the lower limit current intensity, the upper limit stimulation frequency, the upper limit pulse width, and the upper limit current intensity; the first constraint condition is: ; in, represents the lower limit stimulation frequency, Indicates the lower limit pulse width, Indicates the lower limit current intensity, represents the upper limit of stimulation frequency, Indicates the upper limit pulse width, Indicates the upper limit of current intensity.

[0019] Taking the minimum value of the first objective function as the goal and the first constraint condition as the constraint condition, the particle swarm optimization algorithm is used to solve and obtain the optimal stimulation frequency, optimal pulse width, and optimal current intensity.

[0020] Furthermore, the method of constructing the first objective function according to the three-dimensional function of blood pressure value, the three-dimensional function of heart rate value, and the three-dimensional function of electromyographic signal amplitude includes: According to the three-dimensional function of blood pressure value, heart rate value, myoelectric signal amplitude and the pre-set lower limit blood pressure allowable value , pre-set upper limit of blood pressure , Pre-set lower limit of heart rate , Pre-set upper limit of heart rate , the preset lower limit of the electromyographic signal amplitude allowable value , the pre-set upper limit of the electromyographic signal amplitude allowable value , construct a first deviation function, a second deviation function, and a third deviation function; the first deviation function, the second deviation function, and the third deviation function are respectively: ; ; ; in, represents the first deviation function, represents the second deviation function, represents the third deviation function.

[0021] Construct a first objective function based on the first deviation function, the second deviation function, and the third deviation function ; The first objective function is: .

[0022] Based on the same inventive concept, on the other hand, the present invention also provides a low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system, The system comprises: The data reading module is used to obtain the first data, which includes pre-set neural electrical stimulation baseline parameters, neural electrical stimulation lower limit parameters, and neural electrical stimulation upper limit parameters; the neural electrical stimulation baseline parameters include baseline stimulation frequency, baseline pulse width, and baseline current intensity; the neural electrical stimulation lower limit parameters include lower limit stimulation frequency, lower limit pulse width, and lower limit current intensity; the neural electrical stimulation upper limit parameters include upper limit stimulation frequency, upper limit pulse width, and upper limit current intensity.

[0023] The parameter acquisition module is connected to the data reading module and is used to insert the electric needle into the acupuncture point of the stimulation object, set the actual stimulation frequency of the electric needle as the benchmark stimulation frequency, the actual pulse width as the benchmark pulse width, and the actual current intensity as the benchmark current intensity; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation object through the physiological parameter acquisition device to obtain the benchmark blood pressure value, benchmark heart rate value, and benchmark electromyographic signal amplitude.

[0024] The function fitting module is connected to the parameter acquisition module and is used to dynamically set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle at a preset adjustment interval time, and read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation object through the physiological parameter acquisition device to obtain a three-dimensional function of the blood pressure value, a three-dimensional function of the heart rate value, and a three-dimensional function of the electromyographic signal amplitude.

[0025] The optimization calculation module is connected to the function fitting module and is used to calculate the optimization parameters of the electric acupuncture nerve electrical stimulation based on the three-dimensional function of the blood pressure value, the three-dimensional function of the heart rate value, and the three-dimensional function of the electromyography signal amplitude through a parameter optimization algorithm. The optimization parameters of the electric acupuncture nerve electrical stimulation include the optimal stimulation frequency, the optimal pulse width, and the optimal current intensity.

[0026] Furthermore, the function fitting module includes: The first measurement module is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle as the lower limit stimulation frequency, reference pulse width, and reference current intensity, respectively; adjust the actual stimulation frequency with a preset adjustment interval time as the interval and a preset first frequency as the step size, so that the actual stimulation frequency increases step by step until it is greater than the difference between the upper limit stimulation frequency and the first frequency for the first time, read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object, and obtain a first blood pressure array, a first heart rate array, and a first electromyographic signal amplitude array.

[0027] The second measurement module is connected to the first measurement module, and is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle as the reference stimulation frequency, lower limit pulse width, and reference current intensity, respectively; adjust the actual pulse width with a preset adjustment interval time as the interval and a preset first pulse width as the step size, so that the actual pulse width increases stepwise until it is greater than the difference between the upper limit pulse width and the first pulse width for the first time, read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object, and obtain a second blood pressure array, a second heart rate array, and a second electromyographic signal amplitude array.

[0028] The third measurement module is connected to the second measurement module, and is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric press needle as the reference stimulation frequency, reference pulse width, and lower limit current intensity, respectively; adjust the actual current intensity with a preset adjustment interval time as the interval and a preset first current intensity as the step size, so that the actual current intensity increases step by step until it is greater than the difference between the upper limit current intensity and the first current intensity for the first time, read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object, and obtain a third blood pressure array, a third heart rate array, and a third electromyographic signal amplitude array.

[0029] The fitting calculation module is connected to the third measurement module and is used to obtain a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude through a fitting algorithm based on the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array.

[0030] Furthermore, the fitting calculation module includes: The first fitting module is used to fit the first blood pressure function using a nonlinear least squares fitting algorithm according to the first blood pressure array; the first blood pressure function is expressed as ;in represents the actual stimulation frequency; according to the first heart rate array, a nonlinear least squares fitting algorithm is used to fit the first heart rate function; the first heart rate function is expressed as ; According to the first electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the first electromyographic signal amplitude function; the first electromyographic signal amplitude function is expressed as .

[0031] The second fitting module is connected to the first fitting module and is used to fit the second blood pressure function using a nonlinear least squares fitting algorithm according to the second blood pressure array; the second blood pressure function is expressed as ;in represents the actual pulse width; according to the second heart rate array, a nonlinear least squares fitting algorithm is used to fit the second heart rate function; the second heart rate function is expressed as ; According to the second electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the second electromyographic signal amplitude function; the second electromyographic signal amplitude function is expressed as .

[0032] The third fitting module is connected to the second fitting module and is used to fit the third blood pressure function using a nonlinear least squares fitting algorithm according to the third blood pressure array; the third blood pressure function is expressed as ;in Represents the actual current intensity; According to the third heart rate array, a nonlinear least squares fitting algorithm is used to fit the third heart rate function; The third heart rate function is expressed as ; According to the third myoelectric signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the third myoelectric signal amplitude function; the third myoelectric signal amplitude function is expressed as .

[0033] The three-dimensional construction module is connected to the third fitting module and is used to calculate a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude based on the first blood pressure function, the second blood pressure function, the third blood pressure function, the first heart rate function, the second heart rate function, the third heart rate function, the first electromyographic signal amplitude function, the second electromyographic signal amplitude function, and the third electromyographic signal amplitude function through a three-dimensional construction formula; the three-dimensional construction formula is: ; in, represents the three-dimensional function of blood pressure value, represents the three-dimensional function of heart rate value, represents the three-dimensional function of the EMG signal amplitude, represents the base stimulation frequency, Indicates the reference pulse width, Indicates the reference current intensity, Indicates the baseline blood pressure value, Indicates the baseline heart rate value. Indicates the baseline EMG signal amplitude.

[0034] Furthermore, the optimization calculation module includes: The objective function construction module is used to construct a first objective function based on a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude.

[0035] The constraint condition construction module is connected to the objective function construction module and is used to construct a first constraint condition based on the lower limit stimulation frequency, the lower limit pulse width, the lower limit current intensity, the upper limit stimulation frequency, the upper limit pulse width, and the upper limit current intensity; the first constraint condition is: ; in, represents the lower limit stimulation frequency, Indicates the lower limit pulse width, Indicates the lower limit current intensity, represents the upper limit of stimulation frequency, Indicates the upper limit pulse width, Indicates the upper limit of current intensity.

[0036] The particle swarm optimization module is connected to the constraint condition construction module and is used to use the particle swarm optimization algorithm to solve the optimal stimulation frequency, optimal pulse width, and optimal current intensity with the goal of minimizing the value of the first objective function and the first constraint condition as the constraint condition.

[0037] Furthermore, the objective function construction module includes: Deviation function generation module, used to generate the deviation function according to the three-dimensional function of blood pressure value, heart rate value, myoelectric signal amplitude and the preset lower limit blood pressure allowable value. , pre-set upper limit of blood pressure , Pre-set lower limit of heart rate , Pre-set upper limit of heart rate , the preset lower limit of the electromyographic signal amplitude allowable value , the pre-set upper limit of the electromyographic signal amplitude allowable value , construct a first deviation function, a second deviation function, and a third deviation function; the first deviation function, the second deviation function, and the third deviation function are respectively: ; ; ; in, represents the first deviation function, represents the second deviation function, represents the third deviation function.

[0038] The objective function generation module is connected to the deviation function generation module and is used to construct the first objective function based on the first deviation function, the second deviation function, and the third deviation function. ; The first objective function is: .

[0039] (3) Beneficial effects Compared with the prior art, the present invention has the following beneficial effects: Based on the three-dimensional function of blood pressure value, heart rate value and electromyographic signal amplitude, the optimized parameters of electric acupuncture nerve stimulation are calculated by parameter optimization algorithm, thereby ensuring that the stimulated object maintains the optimal physiological indicators during the nerve stimulation process. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 This is a flowchart of the method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation according to Example 1 of the present invention; Figure 2 This is a schematic diagram of the module composition of the low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system of Example 2 of the present invention. DETAILED DESCRIPTION

[0041] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0042] Before giving examples, it is necessary to explain the application scenarios of the present invention. The present invention is used to dynamically adjust the parameters of neural electrical stimulation according to the physiological indicators of the stimulated object, so as to ensure that the stimulated object maintains the optimal physiological indicators during the neural electrical stimulation process.

[0043] Example 1: Figure 1 As shown, this embodiment provides a method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation, the method comprising the following steps: S1. Obtain first data, which includes pre-set neural electrical stimulation baseline parameters, neural electrical stimulation lower limit parameters, and neural electrical stimulation upper limit parameters; the neural electrical stimulation baseline parameters include baseline stimulation frequency, baseline pulse width, and baseline current intensity; the neural electrical stimulation lower limit parameters include lower limit stimulation frequency, lower limit pulse width, and lower limit current intensity; the neural electrical stimulation upper limit parameters include upper limit stimulation frequency, upper limit pulse width, and upper limit current intensity.

[0044] Exemplarily, the lower limit parameters and upper limit parameters of neural electrical stimulation are obtained by querying the instructions for the electro-acupuncture. Among them, the lower limit stimulation frequency is 1 Hz, the upper limit stimulation frequency is 120 Hz, the lower limit pulse width is 20 microseconds, the upper limit pulse width is 1000 microseconds, the lower limit current intensity is 0 mA, and the upper limit current intensity is 8 mA. The neural electrical stimulation reference parameters are set differently according to the different purposes of use of the electro-acupuncture. According to the pre-set first stimulation frequency segmentation node and the second stimulation frequency segmentation node, the frequency segment between the lower limit stimulation frequency and the upper limit stimulation frequency is divided into a first frequency band, a second frequency band and a third frequency band. The pre-set first stimulation frequency segmentation node is 20 Hz, and the second stimulation frequency segmentation node is 80 Hz, so the first frequency band is 1 Hz to 20 Hz, the second frequency band is 20 Hz to 80 Hz, and the third frequency band is 80 Hz to 120 Hz. For electrical nerve stimulation aimed at promoting nerve and muscle regeneration and recovery, the baseline stimulation frequency is set to the frequency corresponding to the midpoint of the first frequency band, i.e., 10.5 Hz. For electrical nerve stimulation aimed at inhibiting pain nerve impulse conduction, the baseline stimulation frequency is set to the frequency corresponding to the midpoint of the third frequency band, i.e., 100 Hz. For electrical nerve stimulation for other purposes, the baseline stimulation frequency is set to the frequency corresponding to the midpoint of the second frequency band, i.e., 50 Hz. The pulse width segment between the lower and upper pulse width limits is divided into a first pulse width segment, a second pulse width segment, and a third pulse width segment based on pre-set first and second pulse width segment points. The pre-set first pulse width segment point is 100 microseconds, and the pre-set second pulse width segment point is 200 microseconds. As a result, the first pulse width segment is 20 to 100 microseconds, the second pulse width segment is 100 to 200 microseconds, and the third pulse width segment is 200 to 1000 microseconds. For electrical nerve stimulation aimed at promoting the recovery of nerve and muscle regeneration function, the reference pulse width is set to the pulse width corresponding to the midpoint of the third pulse width segment, i.e., 600 microseconds; for electrical nerve stimulation aimed at inhibiting the conduction of painful nerve impulses, the reference pulse width is set to the pulse width corresponding to the midpoint of the first pulse width segment, i.e., 60 microseconds; for electrical nerve stimulation for other purposes, the reference pulse width is set to the pulse width corresponding to the midpoint of the second pulse width segment, i.e., 150 microseconds. The reference current intensity is set to the average of the lower limit current intensity and the upper limit current intensity, i.e., 4 mA. In this embodiment, the reference stimulation frequency is set to 10.5 Hz, the reference pulse width is set to 600 microseconds, and the reference current intensity is set to 4 mA.

[0045] S2, insert the electric needle into the acupuncture point of the stimulation object, set the actual stimulation frequency of the electric needle to the reference stimulation frequency, the actual pulse width to the reference pulse width, and the actual current intensity to the reference current intensity; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation object through the physiological parameter acquisition device to obtain the reference blood pressure value, reference heart rate value, and reference electromyographic signal amplitude.

[0046] Exemplarily, an electric needle is inserted into the acupuncture point of the stimulation object, and the actual stimulation frequency of the electric needle is set to 10.5 Hz, the actual pulse width is 600 microseconds, and the actual current intensity is 4 mA. The electric needle is used to perform electrical nerve stimulation on the acupuncture point of the stimulation object. Within a pre-set sampling time segment, the blood pressure value of the stimulation object is measured using a blood pressure meter with a pre-set first interval time as the sampling interval to obtain a blood pressure data sequence, and the blood pressure data sequence is averaged to obtain a baseline blood pressure value. The blood pressure value is diastolic pressure, and the unit is mmHg. Within a pre-set sampling time segment, the heart rate value of the stimulation object is measured using a heart rate meter with a pre-set first interval time as the sampling interval to obtain a heart rate data sequence, and the heart rate data sequence is averaged to obtain a baseline heart rate value. The unit of the heart rate value is beats per minute. Within a pre-set sampling time segment, an electromyography measuring instrument is used to measure the maximum value of the electromyographic signal amplitude of the stimulated subject within the first interval with a pre-set first interval as the sampling interval, thereby obtaining an electromyographic signal amplitude data sequence. The electromyographic signal amplitude data sequence is averaged to obtain a baseline electromyographic signal amplitude. The electromyographic signal amplitude is measured in millivolts. The pre-set sampling time segment is 300 seconds, and the first interval is 20 seconds. Since the first interval is less than 1 minute, when measuring the heart rate of the stimulated subject, the heart rate of the stimulated subject is calculated by proportionally amplifying the number of heartbeats of the stimulated subject obtained by directly measuring the number of heartbeats within 20 seconds.

[0047] S3, dynamically setting the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle with the preset adjustment interval time as the interval, and reading the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object through the physiological parameter acquisition device to obtain a three-dimensional function of the blood pressure value, a three-dimensional function of the heart rate value, and a three-dimensional function of the electromyographic signal amplitude.

[0048] S4. Based on the three-dimensional function of blood pressure value, heart rate value and electromyographic signal amplitude, the optimized parameters of electric acupuncture nerve stimulation are calculated by parameter optimization algorithm. The optimized parameters of electric acupuncture nerve stimulation include optimal stimulation frequency, optimal pulse width and optimal current intensity.

[0049] For example, based on the three-dimensional function of blood pressure value, heart rate value and electromyography signal amplitude, the optimal stimulation frequency is calculated to be 14.2 Hz, the optimal pulse width is 582.8 microseconds and the optimal current intensity is 3.9 mA through the parameter optimization algorithm.

[0050] Furthermore, the method of dynamically setting the actual stimulation frequency, actual pulse width, and actual current intensity of the electric acupuncture needle at a preset adjustment interval time, and obtaining a three-dimensional function of blood pressure value, heart rate value, and electromyographic signal amplitude by reading the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject through a physiological parameter acquisition device includes:

[0051] The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the lower limit stimulation frequency, reference pulse width, and reference current intensity, respectively; the actual stimulation frequency is adjusted with a preset adjustment interval time as the interval and a preset first frequency as the step size, so that the actual stimulation frequency increases stepwise until it is greater than the difference between the upper limit stimulation frequency and the first frequency for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object are read to obtain a first blood pressure array, a first heart rate array, and a first electromyographic signal amplitude array.

[0052] Exemplarily, the actual stimulation frequency, actual pulse width, and actual current intensity of the electric acupuncture needle are set to 1 Hz, 600 microseconds, and 4 mA, respectively. The preset adjustment time interval is 60 seconds, and the preset first frequency is 10 Hz. The difference between the upper limit stimulation frequency and the first frequency is 110 Hz. Therefore, at a time interval of 60 seconds, the actual stimulation frequency is stepwise set to 1 Hz, 11 Hz, 21 Hz, 31 Hz, 41 Hz, 51 Hz, 61 Hz, 71 Hz, 81 Hz, 91 Hz, 101 Hz, and 111 Hz, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject are read to obtain a first blood pressure array, a first heart rate array, and a first electromyographic signal amplitude array.

[0053] The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the reference stimulation frequency, lower limit pulse width, and reference current intensity, respectively; the actual pulse width is adjusted with a preset adjustment interval time as the interval and a preset first pulse width as the step size, so that the actual pulse width increases stepwise until it is greater than the difference between the upper limit pulse width and the first pulse width for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object are read to obtain a second blood pressure array, a second heart rate array, and a second electromyographic signal amplitude array.

[0054] Exemplarily, the actual stimulation frequency, actual pulse width, and actual current intensity of the electric acupuncture needle are set to 10.5 Hz, 20 microseconds, and 4 mA, respectively. The preset adjustment time interval is 60 seconds, and the preset first pulse width is 100 microseconds. The difference between the upper limit pulse width and the first pulse width is 900 microseconds. Therefore, at a time interval of 60 seconds, the actual pulse width is stepwise set to 20 microseconds, 120 microseconds, 220 microseconds, 320 microseconds, 420 microseconds, 520 microseconds, 620 microseconds, 720 microseconds, 820 microseconds, and 920 microseconds. The blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject are read to obtain a second blood pressure array, a second heart rate array, and a second electromyographic signal amplitude array.

[0055] The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the reference stimulation frequency, reference pulse width, and lower limit current intensity, respectively; the actual current intensity is adjusted with a preset adjustment interval time as the interval and a preset first current intensity as the step size, so that the actual current intensity increases stepwise until it is greater than the difference between the upper limit current intensity and the first current intensity for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object are read to obtain a third blood pressure array, a third heart rate array, and a third electromyographic signal amplitude array.

[0056] Exemplarily, the actual stimulation frequency, actual pulse width, and actual current intensity of the electric acupuncture needle are set to 10.5 Hz, 600 microseconds, and 0 mA, respectively. The preset adjustment time interval is 60 seconds, and the preset first current intensity is 1 mA. The difference between the upper limit current intensity and the first current intensity is 7 mA. Therefore, at a time interval of 60 seconds, the actual current intensity is stepwise set to 0 mA, 1 mA, 2 mA, 3 mA, 4 mA, 5 mA, 6 mA, 7 mA, and 8 mA. The blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject are read to obtain a third blood pressure array, a third heart rate array, and a third electromyographic signal amplitude array.

[0057] According to the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array, a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude are obtained through a fitting algorithm.

[0058] Furthermore, the method of obtaining a three-dimensional function of blood pressure values, a three-dimensional function of heart rate values, and a three-dimensional function of electromyographic signal amplitudes by a fitting algorithm based on the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array includes: According to the first blood pressure array, a nonlinear least squares fitting algorithm is used to fit a first blood pressure function; the first blood pressure function is expressed as ;in represents the actual stimulation frequency; according to the first heart rate array, a nonlinear least squares fitting algorithm is used to fit the first heart rate function; the first heart rate function is expressed as ; According to the first electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the first electromyographic signal amplitude function; the first electromyographic signal amplitude function is expressed as .

[0059] According to the second blood pressure array, a nonlinear least squares fitting algorithm is used to fit the second blood pressure function; the second blood pressure function is expressed as ;in represents the actual pulse width; according to the second heart rate array, a nonlinear least squares fitting algorithm is used to fit the second heart rate function; the second heart rate function is expressed as ; According to the second electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the second electromyographic signal amplitude function; the second electromyographic signal amplitude function is expressed as .

[0060] According to the third blood pressure array, a nonlinear least squares fitting algorithm is used to fit a third blood pressure function; the third blood pressure function is expressed as ;in Represents the actual current intensity; According to the third heart rate array, a nonlinear least squares fitting algorithm is used to fit the third heart rate function; The third heart rate function is expressed as ; According to the third myoelectric signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the third myoelectric signal amplitude function; the third myoelectric signal amplitude function is expressed as .

[0061] According to the first blood pressure function, the second blood pressure function, the third blood pressure function, the first heart rate function, the second heart rate function, the third heart rate function, the first electromyographic signal amplitude function, the second electromyographic signal amplitude function, and the third electromyographic signal amplitude function, a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude are calculated by a three-dimensional construction formula; the three-dimensional construction formula is: ; in, represents the three-dimensional function of blood pressure value, represents the three-dimensional function of heart rate value, represents the three-dimensional function of the EMG signal amplitude, represents the base stimulation frequency, Indicates the reference pulse width, Indicates the reference current intensity, Indicates the baseline blood pressure value, Indicates the baseline heart rate value. Indicates the baseline EMG signal amplitude.

[0062] Exemplarily, the obtained first blood pressure function, first heart rhythm function, first electromyographic signal amplitude function, second blood pressure function, second heart rhythm function, second electromyographic signal amplitude function, third blood pressure function, third heart rhythm function, and third electromyographic signal amplitude function are respectively: ; ; ; ; ; ; ; ; .

[0063] Furthermore, the method of calculating the optimized parameters of the electroacupuncture nerve electrical stimulation by a parameter optimization algorithm based on the three-dimensional function of the blood pressure value, the three-dimensional function of the heart rate value, and the three-dimensional function of the electromyographic signal amplitude, wherein the optimized parameters of the electroacupuncture nerve electrical stimulation include the optimal stimulation frequency, the optimal pulse width, and the optimal current intensity, includes: The first objective function is constructed according to the three-dimensional function of blood pressure value, the three-dimensional function of heart rate value and the three-dimensional function of electromyographic signal amplitude.

[0064] A first constraint condition is constructed based on the lower limit stimulation frequency, the lower limit pulse width, the lower limit current intensity, the upper limit stimulation frequency, the upper limit pulse width, and the upper limit current intensity; the first constraint condition is: ; in, represents the lower limit stimulation frequency, Indicates the lower limit pulse width, Indicates the lower limit current intensity, represents the upper limit of stimulation frequency, Indicates the upper limit pulse width, Indicates the upper limit of current intensity.

[0065] Exemplarily, the lower limit stimulation frequency is 1 Hz, the upper limit stimulation frequency is 120 Hz, the lower limit pulse width is 20 microseconds, the upper limit pulse width is 1000 microseconds, the lower limit current intensity is 0 mA, and the upper limit current intensity is 8 mA.

[0066] Taking the minimum value of the first objective function as the goal and the first constraint condition as the constraint condition, the particle swarm optimization algorithm is used to solve and obtain the optimal stimulation frequency, optimal pulse width, and optimal current intensity.

[0067] Furthermore, the method of constructing the first objective function according to the three-dimensional function of blood pressure value, the three-dimensional function of heart rate value, and the three-dimensional function of electromyographic signal amplitude includes: According to the three-dimensional function of blood pressure value, heart rate value, myoelectric signal amplitude and the pre-set lower limit blood pressure allowable value , pre-set upper limit of blood pressure , Pre-set lower limit of heart rate , Pre-set upper limit of heart rate , the preset lower limit of the electromyographic signal amplitude allowable value , the pre-set upper limit of the electromyographic signal amplitude allowable value , construct a first deviation function, a second deviation function, and a third deviation function; the first deviation function, the second deviation function, and the third deviation function are respectively: ; ; ; in, represents the first deviation function, represents the second deviation function, represents the third deviation function.

[0068] Construct a first objective function based on the first deviation function, the second deviation function, and the third deviation function ; The first objective function is: .

[0069] For example, the preset lower limit blood pressure allowable value The upper limit of blood pressure is 60 mmHg. 90 mmHg, the pre-set lower limit of heart rate allowed The upper limit of the heart rate is 60 beats per minute. The lower limit of the myoelectric signal amplitude is 100 times per minute. 0.01 mV, the upper limit of the EMG signal amplitude allowed by the preset is 5 mV.

[0070] Example 2: Based on the same inventive concept, Figure 2 As shown, this embodiment also provides a low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system, the system comprising: The data reading module is used to obtain the first data, which includes pre-set neural electrical stimulation baseline parameters, neural electrical stimulation lower limit parameters, and neural electrical stimulation upper limit parameters; the neural electrical stimulation baseline parameters include baseline stimulation frequency, baseline pulse width, and baseline current intensity; the neural electrical stimulation lower limit parameters include lower limit stimulation frequency, lower limit pulse width, and lower limit current intensity; the neural electrical stimulation upper limit parameters include upper limit stimulation frequency, upper limit pulse width, and upper limit current intensity.

[0071] The parameter acquisition module is connected to the data reading module and is used to insert the electric needle into the acupuncture point of the stimulation object, set the actual stimulation frequency of the electric needle as the benchmark stimulation frequency, the actual pulse width as the benchmark pulse width, and the actual current intensity as the benchmark current intensity; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation object through the physiological parameter acquisition device to obtain the benchmark blood pressure value, benchmark heart rate value, and benchmark electromyographic signal amplitude.

[0072] The function fitting module is connected to the parameter acquisition module and is used to dynamically set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle at a preset adjustment interval time, and read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation object through the physiological parameter acquisition device to obtain a three-dimensional function of the blood pressure value, a three-dimensional function of the heart rate value, and a three-dimensional function of the electromyographic signal amplitude.

[0073] The optimization calculation module is connected to the function fitting module and is used to calculate the optimization parameters of the electric acupuncture nerve electrical stimulation based on the three-dimensional function of the blood pressure value, the three-dimensional function of the heart rate value, and the three-dimensional function of the electromyography signal amplitude through a parameter optimization algorithm. The optimization parameters of the electric acupuncture nerve electrical stimulation include the optimal stimulation frequency, the optimal pulse width, and the optimal current intensity.

[0074] Furthermore, the function fitting module includes: The first measurement module is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle as the lower limit stimulation frequency, reference pulse width, and reference current intensity, respectively; adjust the actual stimulation frequency with a preset adjustment interval time as the interval and a preset first frequency as the step size, so that the actual stimulation frequency increases step by step until it is greater than the difference between the upper limit stimulation frequency and the first frequency for the first time, read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object, and obtain a first blood pressure array, a first heart rate array, and a first electromyographic signal amplitude array.

[0075] The second measurement module is connected to the first measurement module, and is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle as the reference stimulation frequency, lower limit pulse width, and reference current intensity, respectively; adjust the actual pulse width with a preset adjustment interval time as the interval and a preset first pulse width as the step size, so that the actual pulse width increases stepwise until it is greater than the difference between the upper limit pulse width and the first pulse width for the first time, read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object, and obtain a second blood pressure array, a second heart rate array, and a second electromyographic signal amplitude array.

[0076] The third measurement module is connected to the second measurement module, and is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric press needle as the reference stimulation frequency, reference pulse width, and lower limit current intensity, respectively; adjust the actual current intensity with a preset adjustment interval time as the interval and a preset first current intensity as the step size, so that the actual current intensity increases step by step until it is greater than the difference between the upper limit current intensity and the first current intensity for the first time, read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated object, and obtain a third blood pressure array, a third heart rate array, and a third electromyographic signal amplitude array.

[0077] The fitting calculation module is connected to the third measurement module and is used to obtain a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude through a fitting algorithm based on the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array.

[0078] Furthermore, the fitting calculation module includes: The first fitting module is used to fit the first blood pressure function using a nonlinear least squares fitting algorithm according to the first blood pressure array; the first blood pressure function is expressed as ;in represents the actual stimulation frequency; according to the first heart rate array, a nonlinear least squares fitting algorithm is used to fit the first heart rate function; the first heart rate function is expressed as ; According to the first electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the first electromyographic signal amplitude function; the first electromyographic signal amplitude function is expressed as .

[0079] The second fitting module is connected to the first fitting module and is used to fit the second blood pressure function using a nonlinear least squares fitting algorithm according to the second blood pressure array; the second blood pressure function is expressed as ;in represents the actual pulse width; according to the second heart rate array, a nonlinear least squares fitting algorithm is used to fit the second heart rate function; the second heart rate function is expressed as ; According to the second electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the second electromyographic signal amplitude function; the second electromyographic signal amplitude function is expressed as .

[0080] The third fitting module is connected to the second fitting module and is used to fit the third blood pressure function using a nonlinear least squares fitting algorithm according to the third blood pressure array; the third blood pressure function is expressed as ;in Represents the actual current intensity; According to the third heart rate array, a nonlinear least squares fitting algorithm is used to fit the third heart rate function; The third heart rate function is expressed as ; According to the third myoelectric signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the third myoelectric signal amplitude function; the third myoelectric signal amplitude function is expressed as .

[0081] The three-dimensional construction module is connected to the third fitting module and is used to calculate a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude based on the first blood pressure function, the second blood pressure function, the third blood pressure function, the first heart rate function, the second heart rate function, the third heart rate function, the first electromyographic signal amplitude function, the second electromyographic signal amplitude function, and the third electromyographic signal amplitude function through a three-dimensional construction formula; the three-dimensional construction formula is: ; in, represents the three-dimensional function of blood pressure value, represents the three-dimensional function of heart rate value, represents the three-dimensional function of the EMG signal amplitude, represents the base stimulation frequency, Indicates the reference pulse width, Indicates the reference current intensity, Indicates the baseline blood pressure value, Indicates the baseline heart rate value. Indicates the baseline EMG signal amplitude.

[0082] Furthermore, the optimization calculation module includes: The objective function construction module is used to construct a first objective function based on a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude.

[0083] The constraint condition construction module is connected to the objective function construction module and is used to construct a first constraint condition based on the lower limit stimulation frequency, the lower limit pulse width, the lower limit current intensity, the upper limit stimulation frequency, the upper limit pulse width, and the upper limit current intensity; the first constraint condition is: ; in, represents the lower limit stimulation frequency, Indicates the lower limit pulse width, Indicates the lower limit current intensity, represents the upper limit of stimulation frequency, Indicates the upper limit pulse width, Indicates the upper limit of current intensity.

[0084] The particle swarm optimization module is connected to the constraint condition construction module and is used to use the particle swarm optimization algorithm to solve the optimal stimulation frequency, optimal pulse width, and optimal current intensity with the goal of minimizing the value of the first objective function and the first constraint condition as the constraint condition.

[0085] Furthermore, the objective function construction module includes: Deviation function generation module, used to generate the deviation function according to the three-dimensional function of blood pressure value, heart rate value, myoelectric signal amplitude and the preset lower limit blood pressure allowable value. , pre-set upper limit of blood pressure , Pre-set lower limit of heart rate , Pre-set upper limit of heart rate , the preset lower limit of the electromyographic signal amplitude allowable value , the pre-set upper limit of the electromyographic signal amplitude allowable value , construct a first deviation function, a second deviation function, and a third deviation function; the first deviation function, the second deviation function, and the third deviation function are respectively: ; ; ; in, represents the first deviation function, represents the second deviation function, represents the third deviation function.

[0086] The objective function generation module is connected to the deviation function generation module and is used to construct the first objective function based on the first deviation function, the second deviation function, and the third deviation function. ; The first objective function is: .

[0087] It should be noted that, regarding the system in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated on here.

[0088] Finally, it should be noted that although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art can still modify the technical solutions described in the aforementioned embodiments, or make equivalent substitutions for some of the technical features therein. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation, characterized in that: The method comprises the following steps: S1, obtaining first data, the first data including pre-set neural electrical stimulation baseline parameters, neural electrical stimulation lower limit parameters, and neural electrical stimulation upper limit parameters; the neural electrical stimulation baseline parameters include a baseline stimulation frequency, a baseline pulse width, and a baseline current intensity; the neural electrical stimulation lower limit parameters include a lower limit stimulation frequency, a lower limit pulse width, and a lower limit current intensity; the neural electrical stimulation upper limit parameters include an upper limit stimulation frequency, an upper limit pulse width, and an upper limit current intensity; S2, inserting the electro-pressure needle into the acupuncture point of the stimulation subject, setting the actual stimulation frequency of the electro-pressure needle as the reference stimulation frequency, the actual pulse width as the reference pulse width, and the actual current intensity as the reference current intensity; reading the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation subject through the physiological parameter acquisition device to obtain the reference blood pressure value, reference heart rate value, and reference electromyographic signal amplitude; S3, dynamically setting the actual stimulation frequency, actual pulse width, and actual current intensity of the electro-stimulation needle at a preset adjustment interval, and reading the blood pressure, heart rate, and electromyographic signal amplitude of the stimulated subject through a physiological parameter acquisition device to obtain a three-dimensional function of the blood pressure, heart rate, and electromyographic signal amplitude; S4. Based on the three-dimensional function of blood pressure value, heart rate value and electromyographic signal amplitude, the optimized parameters of electric acupuncture nerve stimulation are calculated by parameter optimization algorithm. The optimized parameters of electric acupuncture nerve stimulation include optimal stimulation frequency, optimal pulse width and optimal current intensity.

2. The method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation according to claim 1, wherein: The method of dynamically setting the actual stimulation frequency, actual pulse width, and actual current intensity of the electro-acupuncture at a preset adjustment interval, and obtaining a three-dimensional function of blood pressure, heart rate, and electromyographic signal amplitude by reading the blood pressure, heart rate, and electromyographic signal amplitude of the stimulated subject through a physiological parameter acquisition device includes: The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the lower limit stimulation frequency, the reference pulse width, and the reference current intensity, respectively; the actual stimulation frequency is adjusted with a preset adjustment interval time as the interval and a preset first frequency as the step size, so that the actual stimulation frequency increases stepwise until it is greater than the difference between the upper limit stimulation frequency and the first frequency for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject are read to obtain a first blood pressure array, a first heart rate array, and a first electromyographic signal amplitude array; The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the reference stimulation frequency, the lower limit pulse width, and the reference current intensity, respectively; the actual pulse width is adjusted with a preset adjustment interval time as the interval and a preset first pulse width as the step size, so that the actual pulse width increases stepwise until it is greater than the difference between the upper limit pulse width and the first pulse width for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject are read to obtain a second blood pressure array, a second heart rate array, and a second electromyographic signal amplitude array; The actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle are set as the reference stimulation frequency, reference pulse width, and lower limit current intensity, respectively; the actual current intensity is adjusted with a preset adjustment interval time as an interval and a preset first current intensity as a step size, so that the actual current intensity increases stepwise until it is greater than the difference between the upper limit current intensity and the first current intensity for the first time, and the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject are read to obtain a third blood pressure array, a third heart rate array, and a third electromyographic signal amplitude array; According to the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array, a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude are obtained through a fitting algorithm.

3. The method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation according to claim 2, wherein: The method of obtaining a three-dimensional function of blood pressure values, a three-dimensional function of heart rate values, and a three-dimensional function of electromyographic signal amplitudes by a fitting algorithm based on the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array includes: According to the first blood pressure array, a nonlinear least squares fitting algorithm is used to fit a first blood pressure function; the first blood pressure function is expressed as ;in represents the actual stimulation frequency; according to the first heart rate array, a nonlinear least squares fitting algorithm is used to fit the first heart rate function; the first heart rate function is expressed as ; According to the first electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the first electromyographic signal amplitude function; the first electromyographic signal amplitude function is expressed as ; According to the second blood pressure array, a nonlinear least squares fitting algorithm is used to fit the second blood pressure function; the second blood pressure function is expressed as ;in represents the actual pulse width; according to the second heart rate array, a nonlinear least squares fitting algorithm is used to fit the second heart rate function; the second heart rate function is expressed as ; According to the second electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the second electromyographic signal amplitude function; the second electromyographic signal amplitude function is expressed as ; According to the third blood pressure array, a nonlinear least squares fitting algorithm is used to fit a third blood pressure function; the third blood pressure function is expressed as ;in Represents the actual current intensity; According to the third heart rate array, a nonlinear least squares fitting algorithm is used to fit the third heart rate function; The third heart rate function is expressed as ; According to the third myoelectric signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the third myoelectric signal amplitude function; the third myoelectric signal amplitude function is expressed as ; According to the first blood pressure function, the second blood pressure function, the third blood pressure function, the first heart rate function, the second heart rate function, the third heart rate function, the first electromyographic signal amplitude function, the second electromyographic signal amplitude function, and the third electromyographic signal amplitude function, a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude are calculated by a three-dimensional construction formula; the three-dimensional construction formula is: ； in, represents the three-dimensional function of blood pressure value, represents the three-dimensional function of heart rate value, represents the three-dimensional function of the EMG signal amplitude, represents the base stimulation frequency, Indicates the reference pulse width, Indicates the reference current intensity, Indicates the baseline blood pressure value, Indicates the baseline heart rate value. Indicates the baseline EMG signal amplitude.

4. The method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation according to claim 3, wherein: The method of calculating the electroacupuncture nerve electrical stimulation optimization parameters by a parameter optimization algorithm based on the three-dimensional function of blood pressure value, heart rate value, and electromyographic signal amplitude, wherein the electroacupuncture nerve electrical stimulation optimization parameters include the optimal stimulation frequency, the optimal pulse width, and the optimal current intensity, includes: Constructing a first objective function based on a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude; A first constraint condition is constructed based on the lower limit stimulation frequency, the lower limit pulse width, the lower limit current intensity, the upper limit stimulation frequency, the upper limit pulse width, and the upper limit current intensity; the first constraint condition is: ； in, represents the lower limit stimulation frequency, Indicates the lower limit pulse width, Indicates the lower limit current intensity, represents the upper limit of stimulation frequency, Indicates the upper limit pulse width, Indicates the upper limit current intensity; Taking the minimum value of the first objective function as the goal and the first constraint condition as the constraint condition, the particle swarm optimization algorithm is used to solve and obtain the optimal stimulation frequency, optimal pulse width, and optimal current intensity.

5. The method for optimizing parameters of low-stimulation electric acupuncture nerve electrical stimulation according to claim 4, wherein: The method for constructing the first objective function according to the three-dimensional function of blood pressure value, the three-dimensional function of heart rate value, and the three-dimensional function of electromyographic signal amplitude includes: According to the three-dimensional function of blood pressure value, heart rate value, myoelectric signal amplitude and the pre-set lower limit blood pressure allowable value , pre-set upper limit of blood pressure , Pre-set lower limit of heart rate , Pre-set upper limit of heart rate , the preset lower limit of the electromyographic signal amplitude allowable value , the pre-set upper limit of the electromyographic signal amplitude allowable value , construct a first deviation function, a second deviation function, and a third deviation function; the first deviation function, the second deviation function, and the third deviation function are respectively: ； ； ； in, represents the first deviation function, represents the second deviation function, represents the third deviation function; Construct a first objective function based on the first deviation function, the second deviation function, and the third deviation function ; The first objective function is: 。 6. Low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system, characterized by: The system comprises: A data reading module is used to obtain first data, wherein the first data includes a preset neural electrical stimulation baseline parameter, a neural electrical stimulation lower limit parameter, and a neural electrical stimulation upper limit parameter; the neural electrical stimulation baseline parameter includes a baseline stimulation frequency, a baseline pulse width, and a baseline current intensity; the neural electrical stimulation lower limit parameter includes a lower limit stimulation frequency, a lower limit pulse width, and a lower limit current intensity; the neural electrical stimulation upper limit parameter includes an upper limit stimulation frequency, an upper limit pulse width, and an upper limit current intensity; The parameter acquisition module is connected to the data reading module and is used to insert the electro-pressure needle into the acupuncture point of the stimulation subject, set the actual stimulation frequency of the electro-pressure needle as the reference stimulation frequency, the actual pulse width as the reference pulse width, and the actual current intensity as the reference current intensity; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulation subject through the physiological parameter acquisition device to obtain the reference blood pressure value, reference heart rate value, and reference electromyographic signal amplitude; The function fitting module is connected to the parameter acquisition module and is used to dynamically set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle at a preset adjustment interval time, and read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject through the physiological parameter acquisition device to obtain a three-dimensional function of the blood pressure value, the heart rate value, and the electromyographic signal amplitude; The optimization calculation module is connected to the function fitting module and is used to calculate the optimization parameters of the electric acupuncture nerve electrical stimulation based on the three-dimensional function of the blood pressure value, the three-dimensional function of the heart rate value, and the three-dimensional function of the electromyography signal amplitude through a parameter optimization algorithm. The optimization parameters of the electric acupuncture nerve electrical stimulation include the optimal stimulation frequency, the optimal pulse width, and the optimal current intensity.

7. The low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system according to claim 6, characterized in that: The function fitting module includes: The first measurement module is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle as the lower limit stimulation frequency, reference pulse width, and reference current intensity, respectively; adjust the actual stimulation frequency at a preset adjustment interval time and a preset first frequency as a step size, so that the actual stimulation frequency increases stepwise until it exceeds the difference between the upper limit stimulation frequency and the first frequency for the first time; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject, and obtain a first blood pressure array, a first heart rate array, and a first electromyographic signal amplitude array; The second measurement module is connected to the first measurement module and is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle as the reference stimulation frequency, lower limit pulse width, and reference current intensity, respectively; adjust the actual pulse width at a preset adjustment interval time and a preset first pulse width as a step size, so that the actual pulse width increases stepwise until it is greater than the difference between the upper limit pulse width and the first pulse width for the first time; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject, and obtain a second blood pressure array, a second heart rate array, and a second electromyographic signal amplitude array; The third measurement module is connected to the second measurement module and is used to set the actual stimulation frequency, actual pulse width, and actual current intensity of the electric needle as the reference stimulation frequency, reference pulse width, and lower limit current intensity, respectively; adjust the actual current intensity at a preset adjustment interval time and a preset first current intensity as a step size, so that the actual current intensity increases stepwise until it is greater than the difference between the upper limit current intensity and the first current intensity for the first time; read the blood pressure value, heart rate value, and electromyographic signal amplitude of the stimulated subject, and obtain a third blood pressure array, a third heart rate array, and a third electromyographic signal amplitude array; The fitting calculation module is connected to the third measurement module and is used to obtain a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude through a fitting algorithm based on the first blood pressure array, the first heart rate array, the first electromyographic signal amplitude array, the second blood pressure array, the second heart rate array, the second electromyographic signal amplitude array, the third blood pressure array, the third heart rate array, and the third electromyographic signal amplitude array.

8. The low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system according to claim 7, characterized in that: The fitting calculation module includes: The first fitting module is used to fit the first blood pressure function using a nonlinear least squares fitting algorithm according to the first blood pressure array; the first blood pressure function is expressed as ;in represents the actual stimulation frequency; according to the first heart rate array, a nonlinear least squares fitting algorithm is used to fit the first heart rate function; the first heart rate function is expressed as ; According to the first electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the first electromyographic signal amplitude function; the first electromyographic signal amplitude function is expressed as ; The second fitting module is connected to the first fitting module and is used to fit the second blood pressure function using a nonlinear least squares fitting algorithm according to the second blood pressure array; the second blood pressure function is expressed as ;in represents the actual pulse width; according to the second heart rate array, a nonlinear least squares fitting algorithm is used to fit the second heart rate function; the second heart rate function is expressed as ; According to the second electromyographic signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the second electromyographic signal amplitude function; the second electromyographic signal amplitude function is expressed as ; The third fitting module is connected to the second fitting module and is used to fit the third blood pressure function using a nonlinear least squares fitting algorithm according to the third blood pressure array; the third blood pressure function is expressed as ;in Represents the actual current intensity; According to the third heart rate array, a nonlinear least squares fitting algorithm is used to fit the third heart rate function; The third heart rate function is expressed as ; According to the third myoelectric signal amplitude array, a nonlinear least squares fitting algorithm is used to fit the third myoelectric signal amplitude function; the third myoelectric signal amplitude function is expressed as ; The three-dimensional construction module is connected to the third fitting module and is used to calculate a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude based on the first blood pressure function, the second blood pressure function, the third blood pressure function, the first heart rate function, the second heart rate function, the third heart rate function, the first electromyographic signal amplitude function, the second electromyographic signal amplitude function, and the third electromyographic signal amplitude function through a three-dimensional construction formula; the three-dimensional construction formula is: ； in, represents the three-dimensional function of blood pressure value, represents the three-dimensional function of heart rate value, represents the three-dimensional function of the EMG signal amplitude, represents the base stimulation frequency, Indicates the reference pulse width, Indicates the reference current intensity, Indicates the baseline blood pressure value, Indicates the baseline heart rate value. Indicates the baseline EMG signal amplitude.

9. The low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system according to claim 8, characterized in that: The optimization calculation module includes: An objective function construction module is used to construct a first objective function based on a three-dimensional function of blood pressure value, a three-dimensional function of heart rate value, and a three-dimensional function of electromyographic signal amplitude; The constraint condition construction module is connected to the objective function construction module and is used to construct a first constraint condition based on the lower limit stimulation frequency, the lower limit pulse width, the lower limit current intensity, the upper limit stimulation frequency, the upper limit pulse width, and the upper limit current intensity; the first constraint condition is: ； in, represents the lower limit stimulation frequency, Indicates the lower limit pulse width, Indicates the lower limit current intensity, represents the upper limit of stimulation frequency, Indicates the upper limit pulse width, Indicates the upper limit current intensity; The particle swarm optimization module is connected to the constraint condition construction module and is used to use the particle swarm optimization algorithm to solve the optimal stimulation frequency, optimal pulse width, and optimal current intensity with the goal of minimizing the value of the first objective function and the first constraint condition as the constraint condition.

10. The low-stimulation electric acupuncture nerve electrical stimulation parameter optimization system according to claim 9, characterized in that: The objective function construction module includes: Deviation function generation module, used to generate the deviation function according to the three-dimensional function of blood pressure value, heart rate value, myoelectric signal amplitude and the preset lower limit blood pressure allowable value. , pre-set upper limit of blood pressure , Pre-set lower limit of heart rate , Pre-set upper limit of heart rate , the preset lower limit of the electromyographic signal amplitude allowable value , the pre-set upper limit of the electromyographic signal amplitude allowable value , construct a first deviation function, a second deviation function, and a third deviation function; the first deviation function, the second deviation function, and the third deviation function are respectively: ； ； ； in, represents the first deviation function, represents the second deviation function, represents the third deviation function; The objective function generation module is connected to the deviation function generation module and is used to construct the first objective function based on the first deviation function, the second deviation function, and the third deviation function. ; The first objective function is: 。

Images