Low-stimulus electro-acupuncture nerve electric stimulation parameter optimization method and system

Abstract

The present application relates to the technical field of electric needle nerve electric stimulation, in particular to a low-stimulation electric needle nerve electric stimulation parameter optimization method and system, which comprises the following steps: obtaining first data; inserting an electric needle into an acupuncture point of a stimulation object, setting the actual stimulation frequency of the electric needle as a reference stimulation frequency, the actual pulse width as a reference pulse width, and the actual current intensity as a reference current intensity; reading a reference blood pressure value, a reference heart rate value, and a reference electromyogram amplitude value; setting the actual stimulation frequency, the actual pulse width, and the actual current intensity of the electric needle dynamically at a preset adjustment interval time, obtaining a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and an electromyogram amplitude value three-dimensional function; and calculating electric needle nerve electric stimulation optimization parameters through a parameter optimization algorithm. The present application calculates electric needle nerve electric stimulation optimization parameters through a parameter optimization algorithm, thereby ensuring that a stimulation object maintains optimal physiological indicators during nerve electric stimulation.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of electro-acupuncture nerve electrical stimulation, and particularly to a low-stimulation electro-acupuncture nerve electrical stimulation parameter optimization method and system. BACKGROUND

[0002] Low-stimulation electro-acupuncture is a new medical technology that combines traditional acupuncture and electrical stimulation treatment, and has the advantages of simple operation and good treatment effect. However, in existing electro-acupuncture treatment, the setting of nerve electrical stimulation parameters often depends on the experience of doctors, lacks scientific optimization methods, and is prone to cause the inadaptation of stimulation intensity to physiological indicators of patients, resulting in problems such as rapid changes in blood pressure, sudden changes in heart rate, and muscle convulsions. How to realize accurate setting of electro-acupuncture nerve electrical stimulation parameters according to individual differences of patients is a technical problem to be solved in the field at present. In view of this technical problem, the existing technology mainly adopts the technical means of dividing the nerve electrical stimulation parameter value range into several sections, and selecting appropriate stimulation parameters in the divided parameter value sections according to different purposes of nerve electrical stimulation by artificial experience. However, this method still has the following defects: first, the parameter setting lacks individualization, only considers the treatment purpose, and fails to fully consider the differences in sensitivity of different patients to electrical stimulation; second, the parameter adjustment lacks real-time feedback, and cannot be dynamically optimized according to the physiological response of patients; third, the parameter selection lacks a systematic method, and it is difficult to ensure that the optimal physiological indicators are maintained during the stimulation process.

[0003] Therefore, it is necessary to adaptively optimize the nerve electrical stimulation parameters based on real-time measurement data of physiological parameters of the stimulation object, so as to ensure that the stimulation object maintains the optimal physiological indicators during the nerve electrical stimulation process. SUMMARY

[0004] (1) Technical problem to be solved

[0005] The purpose of the present application is to provide a low-stimulation electro-acupuncture nerve electrical stimulation parameter optimization method and system to realize adaptive optimization of electro-acupuncture nerve electrical stimulation parameters.

[0006] (2) Technical solution

[0007] To achieve the above-mentioned purpose, the present application provides a low-stimulation electro-acupuncture nerve electrical stimulation parameter optimization method, which comprises the following steps:

[0008] S1, obtaining first data, the first data comprising a preset nerve electric stimulation reference parameter, a nerve electric stimulation lower limit parameter, and a nerve electric stimulation upper limit parameter; the nerve electric stimulation reference parameter comprising a reference stimulation frequency, a reference pulse width, and a reference current intensity; the nerve electric stimulation lower limit parameter comprising a lower limit stimulation frequency, a lower limit pulse width, and a lower limit current intensity; and the nerve electric stimulation upper limit parameter comprising an upper limit stimulation frequency, an upper limit pulse width, and an upper limit current intensity.

[0009] S2, inserting the electric button needle into an acupoint of a stimulation subject, setting an actual stimulation frequency of the electric button needle as the reference stimulation frequency, an actual pulse width as the reference pulse width, and an actual current intensity as the reference current intensity; reading blood pressure values, heart rate values, and muscle electric signal amplitudes of the stimulation subject by using a physiological parameter acquisition device to obtain reference blood pressure values, reference heart rate values, and reference muscle electric signal amplitudes.

[0010] S3, setting the actual stimulation frequency, the actual pulse width, and the actual current intensity of the electric button needle dynamically at preset adjustment interval times as intervals, and reading blood pressure values, heart rate values, and muscle electric signal amplitudes of the stimulation subject by using the physiological parameter acquisition device to obtain a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and a muscle electric signal amplitude three-dimensional function.

[0011] S4, calculating electric button needle nerve electric stimulation optimization parameters by using a parameter optimization algorithm according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the muscle electric signal amplitude three-dimensional function, the electric button needle nerve electric stimulation optimization parameters comprising an optimal stimulation frequency, an optimal pulse width, and an optimal current intensity.

[0012] Further, the method of setting the actual stimulation frequency, the actual pulse width, and the actual current intensity of the electric button needle dynamically at preset adjustment interval times as intervals, and reading blood pressure values, heart rate values, and muscle electric signal amplitudes of the stimulation subject by using the physiological parameter acquisition device to obtain a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and a muscle electric signal amplitude three-dimensional function comprises:

[0013] setting the actual stimulation frequency, the actual pulse width, and the actual current intensity of the electric button needle as the lower limit stimulation frequency, the reference pulse width, and the reference current intensity respectively; adjusting the actual stimulation frequency by using a preset first frequency as a step at preset adjustment interval times as intervals, so that the actual stimulation frequency increases step by step until the first difference between the upper limit stimulation frequency and the first frequency is greater than the actual stimulation frequency, and reading blood pressure values, heart rate values, and muscle electric signal amplitudes of the stimulation subject to obtain a first blood pressure array, a first heart rate array, and a first muscle electric signal amplitude array.

[0014] The actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle are respectively set as the reference stimulation frequency, the lower limit pulse width and the reference current intensity; the actual pulse width is adjusted by taking the preset adjustment interval time as an interval and taking the first pulse width as a step, so that the actual pulse width is increased step by step until the first difference between the upper limit pulse width and the first pulse width, the blood pressure value, the heart rate value and the electromyographic signal amplitude of the stimulation object are read, and the second blood pressure array, the second heart rate array and the second electromyographic signal amplitude array are obtained.

[0015] The actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle are respectively set as the reference stimulation frequency, the lower limit pulse width and the reference current intensity; the actual pulse width is adjusted by taking the preset adjustment interval time as an interval and taking the first pulse width as a step, so that the actual pulse width is increased step by step until the first difference between the upper limit pulse width and the first pulse width, the blood pressure value, the heart rate value and the electromyographic signal amplitude of the stimulation object are read, and the second blood pressure array, the second heart rate array and the second electromyographic signal amplitude array are obtained.

[0016] 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 blood pressure value three-dimensional function, a heart rate value three-dimensional function and an electromyographic signal amplitude three-dimensional function are obtained through a fitting algorithm.

[0017] Further, the method of obtaining the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude three-dimensional function 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 through the fitting algorithm comprises:

[0018] According to the first blood pressure array, a first blood pressure function is obtained through a nonlinear least square fitting algorithm; the first blood pressure function is expressed as ; wherein represents the actual stimulation frequency; according to the first heart rate array, a first heart rate function is obtained through a nonlinear least square fitting algorithm; the first heart rate function is expressed as ; and according to the first electromyographic signal amplitude array, a first electromyographic signal amplitude function is obtained through a nonlinear least square fitting algorithm; the first electromyographic signal amplitude function is expressed as .

[0019] According to the second blood pressure array, a second blood pressure function is obtained through a nonlinear least square fitting algorithm; the second blood pressure function is expressed as ; wherein represents the actual pulse width; a second heart rate function is fitted according to the second heart rate array by using a nonlinear least squares fitting algorithm; the second heart rate function is expressed as

[0020] a third blood pressure function is fitted according to the third blood pressure array by using a nonlinear least squares fitting algorithm; the third blood pressure function is expressed as represents the actual current intensity; a third heart rate function is fitted according to the third heart rate array by using a nonlinear least squares fitting algorithm; the third heart rate function is expressed as

[0021] a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and an electromyogram amplitude three-dimensional function are calculated by a three-dimensional construction formula 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 electromyogram amplitude function, the second electromyogram amplitude function, and the third electromyogram amplitude function; the three-dimensional construction formula is:

[0022]

[0023] represents the blood pressure value three-dimensional function, represents the heart rate value three-dimensional function, represents the electromyogram amplitude three-dimensional function, represents the reference stimulation frequency, represents the reference pulse width, represents the reference current intensity, represents the reference blood pressure value, represents the reference heart rate value, represents the reference electromyogram amplitude.

[0024] Further, the method for calculating the electro-acupuncture nerve electrical stimulation optimization parameters including the optimal stimulation frequency, the optimal pulse width, and the optimal current intensity according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the electromyogram amplitude three-dimensional function by using a parameter optimization algorithm includes:

[0025] ​​​​​​​The first target function is constructed according to a blood pressure value three-dimensional function, a heart rate value three-dimensional function and an electromyogram signal amplitude value three-dimensional function.

[0026] The first constraint condition is constructed according to a lower limit stimulation frequency, a lower limit pulse width, a lower limit current intensity, an upper limit stimulation frequency, an upper limit pulse width and an upper limit current intensity.

[0027] ;

[0028] Wherein, represents the lower limit stimulation frequency, represents the lower limit pulse width, represents the lower limit current intensity, represents the upper limit stimulation frequency, represents the upper limit pulse width, represents the upper limit current intensity.

[0029] The particle swarm optimization algorithm is adopted to obtain the optimal stimulation frequency, the optimal pulse width and the optimal current intensity, with the first target function value being minimum as the target and the first constraint condition being the constraint condition.

[0030] Further, the method of constructing the first target function according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyogram signal amplitude value three-dimensional function comprises:

[0031] The first deviation function, the second deviation function and the third deviation function are constructed according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, the electromyogram signal amplitude value three-dimensional function and the pre-set lower limit blood pressure allowable value , the pre-set upper limit blood pressure allowable value , the pre-set lower limit heart rate allowable value , the pre-set upper limit heart rate allowable value , the pre-set lower limit electromyogram signal amplitude allowable value , the pre-set upper limit electromyogram signal amplitude allowable value ; the first deviation function, the second deviation function and the third deviation function are respectively:

[0032] ;

[0033] ;

[0034] ;

[0035] Wherein, represents the first deviation function, represents the second deviation function, represents the third deviation function.

[0036] The first target function is constructed according to the first deviation function, the second deviation function and the third deviation function The first target function is:

[0037] .

[0038] Based on the same inventive concept, in another aspect, the present application also provides a low-stimulation electric button needle nerve electric stimulation parameter optimization system,

[0039] The system comprises:

[0040] The data reading module is configured to obtain first data, wherein the first data comprises pre-set nerve electric stimulation reference parameters, nerve electric stimulation lower limit parameters and nerve electric stimulation upper limit parameters; the nerve electric stimulation reference parameters comprise reference stimulation frequency, reference pulse width and reference current intensity; the nerve electric stimulation lower limit parameters comprise lower limit stimulation frequency, lower limit pulse width and lower limit current intensity; and the nerve electric stimulation upper limit parameters comprise upper limit stimulation frequency, upper limit pulse width and upper limit current intensity.

[0041] The parameter collection module is connected with the data reading module and is configured to insert the electric button needle into an acupuncture point of a stimulation subject, set the actual stimulation frequency of the electric button 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, and read blood pressure value, heart rate value and muscle electric signal amplitude value of the stimulation subject through a physiological parameter collection device to obtain reference blood pressure value, reference heart rate value and reference muscle electric signal amplitude value.

[0042] The function fitting module is connected with the parameter collection module and is configured to set the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric button needle dynamically at a pre-set adjustment interval time, and read blood pressure value, heart rate value and muscle electric signal amplitude value of the stimulation subject through the physiological parameter collection device to obtain blood pressure value three-dimensional function, heart rate value three-dimensional function and muscle electric signal amplitude value three-dimensional function.

[0043] The optimization calculation module is connected with the function fitting module and is configured to calculate electric button needle nerve electric stimulation optimization parameters according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the muscle electric signal amplitude value three-dimensional function through a parameter optimization algorithm, wherein the electric button needle nerve electric stimulation optimization parameters comprise optimal stimulation frequency, optimal pulse width and optimal current intensity.

[0044] Further, the function fitting module comprises:

[0045] The first measuring module is configured to set the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle as the lower limit stimulation frequency, the reference pulse width and the reference current intensity respectively, adjust the actual stimulation frequency by the preset adjustment interval time as the interval and the preset first frequency as the step length, and make the actual stimulation frequency increase step by step until the first difference between the upper limit stimulation frequency and the first frequency, read the blood pressure value, the heart rate value and the muscle electrical signal amplitude of the stimulation object, and obtain the first blood pressure array, the first heart rate array and the first muscle electrical signal amplitude array.

[0046] The second measuring module is connected with the first measuring module and is configured to set the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle as the reference stimulation frequency, the lower limit pulse width and the reference current intensity respectively, adjust the actual pulse width by the preset adjustment interval time as the interval and the preset first pulse width as the step length, and make the actual pulse width increase step by step until the first difference between the upper limit pulse width and the first pulse width, read the blood pressure value, the heart rate value and the muscle electrical signal amplitude of the stimulation object, and obtain the second blood pressure array, the second heart rate array and the second muscle electrical signal amplitude array.

[0047] The third measuring module is connected with the second measuring module and is configured to set the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle as the reference stimulation frequency, the reference pulse width and the lower limit current intensity respectively, adjust the actual current intensity by the preset adjustment interval time as the interval and the preset first current intensity as the step length, and make the actual current intensity increase step by step until the first difference between the upper limit current intensity and the first current intensity, read the blood pressure value, the heart rate value and the muscle electrical signal amplitude of the stimulation object, and obtain the third blood pressure array, the third heart rate array and the third muscle electrical signal amplitude array.

[0048] The fitting calculation module is connected with the third measuring module and is configured to obtain the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the muscle electrical signal amplitude three-dimensional function by a fitting algorithm according to the first blood pressure array, the first heart rate array, the first muscle electrical signal amplitude array, the second blood pressure array, the second heart rate array, the second muscle electrical signal amplitude array, the third blood pressure array, the third heart rate array and the third muscle electrical signal amplitude array.

[0049] Further, the fitting calculation module comprises:

[0050] The first fitting module is configured to obtain the first blood pressure function by a nonlinear least square fitting algorithm according to the first blood pressure array; the first blood pressure function is expressed as ; wherein represents the actual stimulation frequency; the first heart rate function is obtained by a nonlinear least square fitting algorithm according to the first heart rate array; the first heart rate function is expressed as ; according to the first myoelectric signal amplitude array, a first myoelectric signal amplitude function is fitted by using a nonlinear least square fitting algorithm; the first myoelectric signal amplitude function is expressed as .

[0051] A second fitting module is connected with the first fitting module, and is used for fitting a second blood pressure function according to the second blood pressure array by using a nonlinear least square fitting algorithm; the second blood pressure function is expressed as ; wherein represents an actual pulse width; a second heart rate function is fitted according to the second heart rate array by using a nonlinear least square fitting algorithm; the second heart rate function is expressed as ; a second myoelectric signal amplitude function is fitted according to the second myoelectric signal amplitude array by using a nonlinear least square fitting algorithm; the second myoelectric signal amplitude function is expressed as .

[0052] A third fitting module is connected with the second fitting module, and is used for fitting a third blood pressure function according to the third blood pressure array by using a nonlinear least square fitting algorithm; the third blood pressure function is expressed as ; wherein represents an actual current intensity; a third heart rate function is fitted according to the third heart rate array by using a nonlinear least square fitting algorithm; the third heart rate function is expressed as ; a third myoelectric signal amplitude function is fitted according to the third myoelectric signal amplitude array by using a nonlinear least square fitting algorithm; the third myoelectric signal amplitude function is expressed as .

[0053] A three-dimensional construction module is connected with the third fitting module, and is used for calculating a blood pressure value three-dimensional function, a heart rate value three-dimensional function and a myoelectric signal amplitude three-dimensional function by using a three-dimensional construction formula 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 myoelectric signal amplitude function, the second myoelectric signal amplitude function and the third myoelectric signal amplitude function; the three-dimensional construction formula is:

[0054] ;

[0055] wherein, represents the blood pressure value three-dimensional function, represents the heart rate value three-dimensional function, represents the myoelectric signal amplitude three-dimensional function, represents a reference stimulation frequency, represents a reference pulse width, represents a reference current intensity, represents a reference blood pressure value, a reference heart rate value, a reference electromyogram signal amplitude value.

[0056] Further, the optimization calculation module comprises:

[0057] a target function construction module, configured to construct a first target function according to a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and an electromyogram signal amplitude value three-dimensional function.

[0058] a constraint condition construction module, connected to the target function construction module, configured to construct a first constraint condition according to a lower limit stimulation frequency, a lower limit pulse width, a lower limit current intensity, an upper limit stimulation frequency, an upper limit pulse width, and an upper limit current intensity; the first constraint condition is:

[0059] ;

[0060] wherein, the lower limit stimulation frequency is represented by fmin, the lower limit pulse width is represented by wmin, the lower limit current intensity is represented by imin, the upper limit stimulation frequency is represented by fmax, the upper limit pulse width is represented by wmax, the upper limit current intensity is represented by imax.

[0061] a particle swarm optimization module, connected to the constraint condition construction module, configured to adopt a particle swarm optimization algorithm to obtain an optimal stimulation frequency, an optimal pulse width, and an optimal current intensity, with the first target function value being minimum as a target and the first constraint condition being a constraint condition.

[0062] Further, the target function construction module comprises:

[0063] a deviation function generation module, configured to construct a first deviation function, a second deviation function, and a third deviation function according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, the electromyogram signal amplitude value three-dimensional function, and a pre-set lower limit blood pressure allowable value , a pre-set upper limit blood pressure allowable value , a pre-set lower limit heart rate allowable value , a pre-set upper limit heart rate allowable value , a pre-set lower limit electromyogram signal amplitude allowable value , and a pre-set upper limit electromyogram signal amplitude allowable value ; the first deviation function, the second deviation function, and the third deviation function are respectively:

[0064] ;

[0065] ;

[0066] ;

[0067] wherein, represents a first deviation function, represents a second deviation function, represents a third deviation function.

[0068] a target function generation module, connected with the deviation function generation module, for constructing a first target function according to the first deviation function, the second deviation function and the third deviation function ; the first target function is:

[0069] .

[0070] (3) Advantageous effects

[0071] Compared with the prior art, the advantageous effects of the present application are:

[0072] According to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude value three-dimensional function, the electric acupuncture needle nerve electric stimulation optimization parameters are calculated through a parameter optimization algorithm, so as to ensure that the stimulation object maintains the optimal physiological index in the nerve electric stimulation process. BRIEF DESCRIPTION OF DRAWINGS

[0073] Figure 1 is a flowchart of the low-stimulation electric acupuncture needle nerve electric stimulation parameter optimization method of embodiment 1 of the present application;

[0074] Figure 2 is a module composition schematic diagram of the low-stimulation electric acupuncture needle nerve electric stimulation parameter optimization system of embodiment 2 of the present application. DETAILED DESCRIPTION

[0075] The technical solutions in the embodiments of the present application will be described clearly and completely below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0076] Before examples are given, the application scenario of the concept of the present application needs to be described. The present application is applied to dynamically adjusting nerve electric stimulation parameters according to the physiological index of a stimulation object, so as to ensure that the stimulation object maintains the optimal physiological index in the nerve electric stimulation process.

[0077] Embodiment 1: As shown in the following table, the present embodiment provides a low-stimulation electric acupuncture needle nerve electric stimulation parameter optimization method, which comprises the following steps: Figure 1

[0078] ​S1, acquiring first data, the first data comprising preset nerve electric stimulation reference parameters, nerve electric stimulation lower limit parameters, nerve electric stimulation upper limit parameters; the nerve electric stimulation reference parameters comprising reference stimulation frequency, reference pulse width, reference current intensity; the nerve electric stimulation lower limit parameters comprising lower limit stimulation frequency, lower limit pulse width, lower limit current intensity; the nerve electric stimulation upper limit parameters comprising upper limit stimulation frequency, upper limit pulse width, upper limit current intensity.

[0079] Exemplarily, the lower limit parameter and the upper limit parameter of the nerve electric stimulation are obtained by querying the electric acupuncture needle specification. 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 reference parameter of the nerve electric stimulation is set differently according to different purposes of the electric acupuncture needle. The frequency segment between the lower limit stimulation frequency and the upper limit stimulation frequency is divided into a first frequency segment, a second frequency segment and a third frequency segment according to a first preset stimulation frequency segmentation node and a second preset stimulation frequency segmentation node. The first preset stimulation frequency segmentation node is 20 Hz, and the second preset stimulation frequency segmentation node is 80 Hz, so that the first frequency segment is 1 Hz to 20 Hz, the second frequency segment is 20 Hz to 80 Hz, and the third frequency segment is 80 Hz to 120 Hz. The reference stimulation frequency of the nerve electric stimulation for promoting the functional recovery of nerves and muscles is set as the frequency corresponding to the midpoint of the first frequency segment, i.e. 10.5 Hz; the reference stimulation frequency of the nerve electric stimulation for inhibiting the conduction of pain nerve impulses is set as the frequency corresponding to the midpoint of the third frequency segment, i.e. 100 Hz; and the reference stimulation frequency of the nerve electric stimulation for other purposes is set as the frequency corresponding to the midpoint of the second frequency segment, i.e. 50 Hz. The pulse width segment between the lower limit pulse width and the upper limit pulse width is divided into a first pulse width segment, a second pulse width segment and a third pulse width segment according to a first preset pulse width segmentation node and a second preset pulse width segmentation node. The first preset pulse width segmentation node is 100 microseconds, and the second preset pulse width segmentation node is 200 microseconds, so that the first pulse width segment is 20 microseconds to 100 microseconds, the second pulse width segment is 100 microseconds to 200 microseconds, and the third pulse width segment is 200 microseconds to 1000 microseconds. The reference pulse width of the nerve electric stimulation for promoting the functional recovery of nerves and muscles is set as the pulse width corresponding to the midpoint of the third pulse width segment, i.e. 600 microseconds; the reference pulse width of the nerve electric stimulation for inhibiting the conduction of pain nerve impulses is set as the pulse width corresponding to the midpoint of the first pulse width segment, i.e. 60 microseconds; and the reference pulse width of the nerve electric stimulation for other purposes is set as the pulse width corresponding to the midpoint of the second pulse width segment, i.e. 150 microseconds. The reference current intensity is set as 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 as 10.5 Hz, the reference pulse width is set as 600 microseconds, and the reference current intensity is set as 4 mA for the purpose of promoting the functional recovery of nerves and muscles.

[0080] S2, inserting the button needle into the acupoint of the stimulation object, setting the actual stimulation frequency of the button 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 muscle electrical signal amplitude of the stimulation object by the physiological parameter acquisition device to obtain the reference blood pressure value, reference heart rate value, and reference muscle electrical signal amplitude.

[0081] Exemplarily, the button needle is inserted into the acupoint of the stimulation object, and the actual stimulation frequency of the button needle is set as 10.5 Hz, the actual pulse width is set as 600 microseconds, and the actual current intensity is set as 4 milliamperes. The button needle is used to perform nerve electrical stimulation on the acupoint of the stimulation object. Within a pre-set sampling time segment, a sphygmomanometer is used to measure the blood pressure value of the stimulation object at a pre-set first interval time as a sampling interval to obtain a blood pressure data sequence, and the blood pressure data sequence is averaged to obtain a reference blood pressure value. The blood pressure value is diastolic blood pressure, and the unit is millimeters of mercury. Within the pre-set sampling time segment, a heart rate meter is used to measure the heart rate value of the stimulation object at the pre-set first interval time as a sampling interval to obtain a heart rate data sequence, and the heart rate data sequence is averaged to obtain a reference heart rate value. The unit of the heart rate value is beats per minute. Within the pre-set sampling time segment, an electromyogram measuring instrument is used to measure the maximum value of the muscle electrical signal amplitude of the stimulation object within the first interval time at the pre-set first interval time as a sampling interval to obtain a muscle electrical signal amplitude data sequence, and the muscle electrical signal amplitude data sequence is averaged to obtain a reference muscle electrical signal amplitude. The unit of the muscle electrical signal amplitude is millivolts. The pre-set sampling time segment is 300 seconds, and the first interval time is 20 seconds. Since the first interval time is less than 1 minute, when the heart rate value of the stimulation object is measured, the heart rate value of the stimulation object is calculated by an equal proportion magnification method according to the number of heartbeats of the stimulation object within 20 seconds directly measured.

[0082] S3, setting the actual stimulation frequency, actual pulse width, and actual current intensity of the button needle dynamically at a pre-set adjustment interval time as an interval, and reading the blood pressure value, heart rate value, and muscle electrical signal amplitude of the stimulation object by the physiological parameter acquisition device to obtain a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and a muscle electrical signal amplitude three-dimensional function.

[0083] S4, calculating the button needle nerve electrical stimulation optimization parameters by a parameter optimization algorithm according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the muscle electrical signal amplitude three-dimensional function to obtain the button needle nerve electrical stimulation optimization parameters, wherein the button needle nerve electrical stimulation optimization parameters include an optimal stimulation frequency, an optimal pulse width, and an optimal current intensity.

[0084] Exemplarily, according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the myoelectric signal amplitude three-dimensional function, the optimal stimulation frequency is 14.2 Hz, the optimal pulse width is 582.8 microseconds, and the optimal current intensity is 3.9 mA, which are calculated by a parameter optimization algorithm.

[0085] Further, the method for dynamically setting the actual stimulation frequency, the actual pulse width, and the actual current intensity of the acupressure needle at the preset adjustment interval time interval, and reading the blood pressure value, the heart rate value, and the myoelectric signal amplitude of the stimulation object by the physiological parameter acquisition device to obtain the blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the myoelectric signal amplitude three-dimensional function comprises:

[0086] The actual stimulation frequency, the actual pulse width, and the actual current intensity of the acupressure needle are set to the lower limit stimulation frequency, the reference pulse width, and the reference current intensity, respectively. The actual stimulation frequency is adjusted by the preset first frequency as a step at the preset adjustment interval time interval, so that the actual stimulation frequency is increased step by step until the first difference between the upper limit stimulation frequency and the first frequency, the blood pressure value, the heart rate value, and the myoelectric signal amplitude of the stimulation object are read, and the first blood pressure array, the first heart rate array, and the first myoelectric signal amplitude array are obtained.

[0087] Exemplarily, the actual stimulation frequency, the actual pulse width, and the actual current intensity of the acupressure 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, the actual stimulation frequency is 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 step by step at the time interval of 60 seconds, the blood pressure value, the heart rate value, and the myoelectric signal amplitude of the stimulation object are read, and the first blood pressure array, the first heart rate array, and the first myoelectric signal amplitude array are obtained.

[0088] The actual stimulation frequency, the actual pulse width, and the actual current intensity of the acupressure needle are set to the reference stimulation frequency, the lower limit pulse width, and the reference current intensity, respectively. The actual pulse width is adjusted by the preset first pulse width as a step at the preset adjustment interval time interval, so that the actual pulse width is increased step by step until the first difference between the upper limit pulse width and the first pulse width, the blood pressure value, the heart rate value, and the myoelectric signal amplitude of the stimulation object are read, and the second blood pressure array, the second heart rate array, and the second myoelectric signal amplitude array are obtained.

[0089] Exemplarily, the actual stimulation frequency, the actual pulse width, and the actual current intensity of the acupoint stimulator are set to 10.5 Hz, 20 microseconds, and 4 milliampere respectively. The pre-set adjustment interval is 60 seconds, and the pre-set first pulse width is 100 microseconds. The difference between the upper limit pulse width and the first pulse width is 900 microseconds. Therefore, the actual pulse width is set to 20 microseconds, 120 microseconds, 220 microseconds, 320 microseconds, 420 microseconds, 520 microseconds, 620 microseconds, 720 microseconds, 820 microseconds, and 920 microseconds in a time interval of 60 seconds, and the blood pressure value, the heart rate value, and the muscle electrical signal amplitude of the stimulation subject are read to obtain the second blood pressure array, the second heart rate array, and the second muscle electrical signal amplitude array.

[0090] The actual stimulation frequency, the actual pulse width, and the actual current intensity of the acupoint stimulator are set to the reference stimulation frequency, the reference pulse width, and the lower limit current intensity respectively. The actual current intensity is adjusted in a pre-set adjustment interval and a pre-set first current intensity as a step, and the actual current intensity is increased step by step until the first time that the actual current intensity is greater than the difference between the upper limit current intensity and the first current intensity. The blood pressure value, the heart rate value, and the muscle electrical signal amplitude of the stimulation subject are read to obtain the third blood pressure array, the third heart rate array, and the third muscle electrical signal amplitude array.

[0091] Exemplarily, the actual stimulation frequency, the actual pulse width, and the actual current intensity of the acupoint stimulator are set to 10.5 Hz, 600 microseconds, and 0 milliampere respectively. The pre-set adjustment interval is 60 seconds, and the pre-set first current intensity is 1 milliampere. The difference between the upper limit current intensity and the first current intensity is 7 milliampere. Therefore, the actual current intensity is set to 0 milliampere, 1 milliampere, 2 milliampere, 3 milliampere, 4 milliampere, 5 milliampere, 6 milliampere, 7 milliampere, and 8 milliampere in a time interval of 60 seconds, and the blood pressure value, the heart rate value, and the muscle electrical signal amplitude of the stimulation subject are read to obtain the third blood pressure array, the third heart rate array, and the third muscle electrical signal amplitude array.

[0092] The blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the muscle electrical signal amplitude three-dimensional function are obtained by a fitting algorithm according to the first blood pressure array, the first heart rate array, the first muscle electrical signal amplitude array, the second blood pressure array, the second heart rate array, the second muscle electrical signal amplitude array, the third blood pressure array, the third heart rate array, and the third muscle electrical signal amplitude array.

[0093] Further, the method of obtaining the blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the muscle electrical signal amplitude three-dimensional function according to the first blood pressure array, the first heart rate array, the first muscle electrical signal amplitude array, the second blood pressure array, the second heart rate array, the second muscle electrical signal amplitude array, the third blood pressure array, the third heart rate array, and the third muscle electrical signal amplitude array by a fitting algorithm comprises:

[0094] According to the first blood pressure array, a first blood pressure function is fitted by using a nonlinear least square fitting algorithm; the first blood pressure function is expressed as ; wherein represents an actual stimulation frequency; according to the first heart rate array, a first heart rate function is fitted by using a nonlinear least square fitting algorithm; the first heart rate function is expressed as ; according to the first myoelectric signal amplitude array, a first myoelectric signal amplitude function is fitted by using a nonlinear least square fitting algorithm; the first myoelectric signal amplitude function is expressed as .

[0095] According to the second blood pressure array, a second blood pressure function is fitted by using a nonlinear least square fitting algorithm; the second blood pressure function is expressed as ; wherein represents an actual pulse width; according to the second heart rate array, a second heart rate function is fitted by using a nonlinear least square fitting algorithm; the second heart rate function is expressed as ; according to the second myoelectric signal amplitude array, a second myoelectric signal amplitude function is fitted by using a nonlinear least square fitting algorithm; the second myoelectric signal amplitude function is expressed as .

[0096] According to the third blood pressure array, a third blood pressure function is fitted by using a nonlinear least square fitting algorithm; the third blood pressure function is expressed as ; wherein represents an actual current intensity; according to the third heart rate array, a third heart rate function is fitted by using a nonlinear least square fitting algorithm; the third heart rate function is expressed as ; according to the third myoelectric signal amplitude array, a third myoelectric signal amplitude function is fitted by using a nonlinear least square fitting algorithm; the third myoelectric signal amplitude function is expressed as .

[0097] 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 myoelectric signal amplitude function, the second myoelectric signal amplitude function, the third myoelectric signal amplitude function, a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and a myoelectric signal amplitude three-dimensional function are calculated by using a three-dimensional construction formula; the three-dimensional construction formula is:

[0098] ;

[0099] wherein, represents the blood pressure value three-dimensional function, represents the heart rate value three-dimensional function, represents the myoelectric signal amplitude three-dimensional function, Indicates the baseline stimulation frequency. Indicates the reference pulse width. Indicates the reference current intensity. This represents the baseline blood pressure value. This represents the baseline heart rate value. This indicates the amplitude of the baseline electromyographic signal.

[0100] For example, the obtained first blood pressure function, first heart rhythm function, first electromyography (EMG) signal amplitude function, second blood pressure function, second heart rhythm function, second EMG signal amplitude function, third blood pressure function, third heart rhythm function, and third EMG signal amplitude function are respectively:

[0101] ;

[0102] ;

[0103] ;

[0104] ;

[0105] ;

[0106] ;

[0107] ;

[0108] ;

[0109] .

[0110] Furthermore, the method for calculating optimized parameters for electroacupuncture nerve stimulation based on the three-dimensional functions of blood pressure, heart rate, and electromyography signal amplitude using a parameter optimization algorithm, wherein the optimized parameters include optimal stimulation frequency, optimal pulse width, and optimal current intensity, includes:

[0111] The first objective function is constructed based on the three-dimensional functions of blood pressure, heart rate, and electromyography signal amplitude.

[0112] The first constraint condition is constructed based on the lower limit stimulation frequency, lower limit pulse width, lower limit current intensity, upper limit stimulation frequency, upper limit pulse width, and upper limit current intensity; the first constraint condition is:

[0113] ;

[0114] in, This indicates the lower limit of stimulation frequency. Indicates the lower limit pulse width. Indicates the lower limit current intensity. Indicates the upper limit of stimulation frequency. Indicates the upper limit pulse width. This indicates the upper limit of the current intensity.

[0115] For example, 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.

[0116] With the objective function being minimized and the first constraint condition being applied, the optimal stimulation frequency, optimal pulse width, and optimal current intensity are obtained by using a particle swarm optimization algorithm.

[0117] Furthermore, the method for constructing the first objective function based on the three-dimensional functions of blood pressure, heart rate, and electromyography signal amplitude includes:

[0118] Based on the three-dimensional function of blood pressure, the three-dimensional function of heart rate, the three-dimensional function of electromyographic signal amplitude, and the preset lower limit of permissible blood pressure. Pre-set upper limit blood pressure allowable value Pre-set lower limit heart rate Preset upper limit of allowable heart rate Pre-set lower limit allowable value for electromyographic signal amplitude Pre-set upper limit allowable value for electromyographic signal amplitude 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:

[0119] ;

[0120] ;

[0121] ;

[0122] in, Represents the first deviation function. This represents the second deviation function. This represents the third deviation function.

[0123] 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:

[0124] .

[0125] Exemplarily, the pre-set lower limit blood pressure allowable value is 60 mmHg , the pre-set upper limit blood pressure allowable value is 90 mmHg , the pre-set lower limit heart rate allowable value is 60 beats per minute , the pre-set upper limit heart rate allowable value is 100 beats per minute , the pre-set lower limit electromyogram amplitude allowable value is 0.01 mV , and the pre-set upper limit electromyogram amplitude allowable value is 5 mV.

[0126] Embodiment 2: Based on the same inventive concept, as shown in the following table, this embodiment also provides a low-stimulation electro-acupuncture nerve electrical stimulation parameter optimization system, which comprises: Figure 2 a data reading module for acquiring first data, the first data comprising pre-set nerve electrical stimulation reference parameters, nerve electrical stimulation lower limit parameters, and nerve electrical stimulation upper limit parameters; the nerve electrical stimulation reference parameters comprising reference stimulation frequency, reference pulse width, and reference current intensity; the nerve electrical stimulation lower limit parameters comprising lower limit stimulation frequency, lower limit pulse width, and lower limit current intensity; and the nerve electrical stimulation upper limit parameters comprising upper limit stimulation frequency, upper limit pulse width, and upper limit current intensity.

[0127] a parameter collection module connected to the data reading module, for inserting an electro-acupuncture needle into an acupoint of a stimulation subject, setting the actual stimulation frequency of the electro-acupuncture 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; and reading the blood pressure value, heart rate value, and electromyogram amplitude value of the stimulation subject through a physiological parameter collection device to obtain a reference blood pressure value, a reference heart rate value, and a reference electromyogram amplitude value.

[0128] a function fitting module connected to the parameter collection module, for setting the actual stimulation frequency, actual pulse width, and actual current intensity of the electro-acupuncture needle at a pre-set adjustment interval time, and reading the blood pressure value, heart rate value, and electromyogram amplitude value of the stimulation subject through the physiological parameter collection device to obtain a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and an electromyogram amplitude value three-dimensional function.

[0129] an optimization calculation module connected to the function fitting module, for calculating electro-acupuncture nerve electrical stimulation optimization parameters comprising optimal stimulation frequency, optimal pulse width, and optimal current intensity through a parameter optimization algorithm according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, and the electromyogram amplitude value three-dimensional function.

[0130] Further, the function fitting module comprises:

[0131] a blood pressure value three-dimensional function fitting sub-module for fitting the blood pressure value three-dimensional function according to the blood pressure value of the stimulation subject at the pre-set adjustment interval time.

[0132] ​The first measuring module is configured to set the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle as the lower limit stimulation frequency, the reference pulse width and the reference current intensity respectively, adjust the actual stimulation frequency by the preset adjustment interval time as the interval and the preset first frequency as the step length, and make the actual stimulation frequency increase step by step until the first difference between the upper limit stimulation frequency and the first frequency, read the blood pressure value, the heart rate value and the muscle electrical signal amplitude of the stimulation object, and obtain the first blood pressure array, the first heart rate array and the first muscle electrical signal amplitude array.

[0133] The second measuring module is connected with the first measuring module and is configured to set the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle as the reference stimulation frequency, the lower limit pulse width and the reference current intensity respectively, adjust the actual pulse width by the preset adjustment interval time as the interval and the preset first pulse width as the step length, and make the actual pulse width increase step by step until the first difference between the upper limit pulse width and the first pulse width, read the blood pressure value, the heart rate value and the muscle electrical signal amplitude of the stimulation object, and obtain the second blood pressure array, the second heart rate array and the second muscle electrical signal amplitude array.

[0134] The third measuring module is connected with the second measuring module and is configured to set the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric acupressure needle as the reference stimulation frequency, the reference pulse width and the lower limit current intensity respectively, adjust the actual current intensity by the preset adjustment interval time as the interval and the preset first current intensity as the step length, and make the actual current intensity increase step by step until the first difference between the upper limit current intensity and the first current intensity, read the blood pressure value, the heart rate value and the muscle electrical signal amplitude of the stimulation object, and obtain the third blood pressure array, the third heart rate array and the third muscle electrical signal amplitude array.

[0135] The fitting calculation module is connected with the third measuring module and is configured to obtain the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the muscle electrical signal amplitude three-dimensional function by a fitting algorithm according to the first blood pressure array, the first heart rate array, the first muscle electrical signal amplitude array, the second blood pressure array, the second heart rate array, the second muscle electrical signal amplitude array, the third blood pressure array, the third heart rate array and the third muscle electrical signal amplitude array.

[0136] Further, the fitting calculation module comprises:

[0137] The first fitting module is configured to obtain the first blood pressure function by a nonlinear least square fitting algorithm according to the first blood pressure array; the first blood pressure function is expressed as ; wherein represents the actual stimulation frequency; the first heart rate function is obtained by a nonlinear least square fitting algorithm according to the first heart rate array; the first heart rate function is expressed as ; according to the first myoelectric signal amplitude array, a first myoelectric signal amplitude function is fitted by using a nonlinear least square fitting algorithm; the first myoelectric signal amplitude function is expressed as .

[0138] The second fitting module is connected with the first fitting module, and is used for fitting a second blood pressure function according to the second blood pressure array by using a nonlinear least square fitting algorithm; the second blood pressure function is expressed as ; wherein represents an actual pulse width; a second heart rate function is fitted according to the second heart rate array by using a nonlinear least square fitting algorithm; the second heart rate function is expressed as ; a second myoelectric signal amplitude function is fitted according to the second myoelectric signal amplitude array by using a nonlinear least square fitting algorithm; the second myoelectric signal amplitude function is expressed as .

[0139] The third fitting module is connected with the second fitting module, and is used for fitting a third blood pressure function according to the third blood pressure array by using a nonlinear least square fitting algorithm; the third blood pressure function is expressed as ; wherein represents an actual current intensity; a third heart rate function is fitted according to the third heart rate array by using a nonlinear least square fitting algorithm; the third heart rate function is expressed as ; a third myoelectric signal amplitude function is fitted according to the third myoelectric signal amplitude array by using a nonlinear least square fitting algorithm; the third myoelectric signal amplitude function is expressed as .

[0140] The three-dimensional construction module is connected with the third fitting module, and is used for calculating a blood pressure value three-dimensional function, a heart rate value three-dimensional function and a myoelectric signal amplitude three-dimensional function by using a three-dimensional construction formula 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 myoelectric signal amplitude function, the second myoelectric signal amplitude function and the third myoelectric signal amplitude function; the three-dimensional construction formula is:

[0141] ;

[0142] wherein, represents the blood pressure value three-dimensional function, represents the heart rate value three-dimensional function, represents the myoelectric signal amplitude three-dimensional function, represents a reference stimulation frequency, represents a reference pulse width, represents a reference current intensity, represents a reference blood pressure value, a reference heart rate value, a reference electromyogram signal amplitude.

[0143] Further, the optimization calculation module comprises:

[0144] a target function construction module, configured to construct a first target function according to a blood pressure value three-dimensional function, a heart rate value three-dimensional function, and an electromyogram signal amplitude three-dimensional function.

[0145] a constraint condition construction module, connected to the target function construction module, configured to construct a first constraint condition according to a lower limit stimulation frequency, a lower limit pulse width, a lower limit current intensity, an upper limit stimulation frequency, an upper limit pulse width, and an upper limit current intensity; the first constraint condition is:

[0146] ;

[0147] wherein, the lower limit stimulation frequency is represented by fmin, the lower limit pulse width is represented by wmin, the lower limit current intensity is represented by imin, the upper limit stimulation frequency is represented by fmax, the upper limit pulse width is represented by wmax, the upper limit current intensity is represented by imax.

[0148] a particle swarm optimization module, connected to the constraint condition construction module, configured to adopt a particle swarm optimization algorithm to obtain an optimal stimulation frequency, an optimal pulse width, and an optimal current intensity, with the first target function value being minimum as a target and the first constraint condition being a constraint condition.

[0149] Further, the target function construction module comprises:

[0150] a deviation function generation module, configured to construct a first deviation function, a second deviation function, and a third deviation function according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, the electromyogram signal amplitude three-dimensional function, and a pre-set lower limit blood pressure allowable value , a pre-set upper limit blood pressure allowable value , a pre-set lower limit heart rate allowable value , a pre-set upper limit heart rate allowable value , a pre-set lower limit electromyogram signal amplitude allowable value , and a pre-set upper limit electromyogram signal amplitude allowable value ; the first deviation function, the second deviation function, and the third deviation function are respectively:

[0151] ;

[0152] ;

[0153] ;

[0154] wherein, represents a first deviation function, represents a second deviation function, represents a third deviation function.

[0155] a target function generating module, connected with the deviation function generating module, configured to construct a first target function according to the first deviation function, the second deviation function and the third deviation function ; the first target function is:

[0156] .

[0157] It should be noted that, as for the system in the above-mentioned embodiments, the specific manner in which each module performs operations has been described in detail in the embodiments related to the method, and will not be described in detail here.

[0158] Finally, it should be noted that: although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacements for part of the technical features, any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for optimizing low irritation electro-acupuncture nerve electrical stimulation parameters, characterized in that, The method comprises the following steps: S1, acquiring first data, the first data comprising pre-set nerve electric stimulation reference parameters, nerve electric stimulation lower limit parameters and nerve electric stimulation upper limit parameters; the nerve electric stimulation reference parameters comprising reference stimulation frequency, reference pulse width and reference current intensity; the nerve electric stimulation lower limit parameters comprising lower limit stimulation frequency, lower limit pulse width and lower limit current intensity; the nerve electric stimulation upper limit parameters comprising upper limit stimulation frequency, upper limit pulse width and upper limit current intensity; S2, inserting an electric button needle into an acupuncture point of a stimulation subject, setting the actual stimulation frequency of the electric button 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 muscle electric signal amplitude of the stimulation subject by a physiological parameter acquisition device to obtain the reference blood pressure value, reference heart rate value and reference muscle electric signal amplitude; S3, setting the actual stimulation frequency, actual pulse width and actual current intensity of the electric button needle at intervals of a pre-set adjustment interval time, reading the blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation subject by the physiological parameter acquisition device to obtain a blood pressure value three-dimensional function, a heart rate value three-dimensional function and a muscle electric signal amplitude three-dimensional function; S4, calculating the electric button needle nerve electric stimulation optimization parameters by a parameter optimization algorithm according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the muscle electric signal amplitude three-dimensional function, the electric button needle nerve electric stimulation optimization parameters comprising optimal stimulation frequency, optimal pulse width and optimal current intensity; The method for setting the actual stimulation frequency, actual pulse width and actual current intensity of the electric button needle at intervals of a pre-set adjustment interval time, reading the blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation subject by the physiological parameter acquisition device to obtain a blood pressure value three-dimensional function, a heart rate value three-dimensional function and a muscle electric signal amplitude three-dimensional function comprises: setting the actual stimulation frequency, actual pulse width and actual current intensity of the electric button needle as the lower limit stimulation frequency, reference pulse width and reference current intensity respectively; adjusting the actual stimulation frequency at intervals of a pre-set adjustment interval time and by a pre-set first frequency as a step length, so that the actual stimulation frequency increases step by step until the first difference between the upper limit stimulation frequency and the first frequency, reading the blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation subject to obtain a first blood pressure array, a first heart rate array and a first muscle electric signal amplitude array; setting the actual stimulation frequency, actual pulse width and actual current intensity of the electric button needle as the reference stimulation frequency, lower limit pulse width and reference current intensity respectively; adjusting the actual pulse width at intervals of a pre-set adjustment interval time and by a pre-set first pulse width as a step length, so that the actual pulse width increases step by step until the first difference between the upper limit pulse width and the first pulse width, reading the blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation subject to obtain a second blood pressure array, a second heart rate array and a second muscle electric signal amplitude array; The actual stimulation frequency, the actual pulse width and the actual current intensity of the electro-acupuncture needle are respectively set as a reference stimulation frequency, a reference pulse width and a lower limit current intensity; the actual current intensity is adjusted by a preset adjustment interval time as an interval and a preset first current intensity as a step, so that the actual current intensity is increased step by step until the first difference between the upper limit current intensity and the first current intensity, and the blood pressure value, the heart rate value and the electromyographic signal amplitude of the stimulation object are read to obtain a third blood pressure array, a third heart rate array and a third electromyographic signal amplitude array; The blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude three-dimensional function are obtained through a fitting algorithm 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.

2. The low irritation electro-acupuncture nerve electro-stimulation parameter optimization method of claim 1, wherein, The method for obtaining the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude three-dimensional function 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 through the fitting algorithm comprises: According to the first blood pressure array, a first blood pressure function is fitted by using a nonlinear least square fitting algorithm; the first blood pressure function is expressed as ; wherein represents an actual stimulation frequency; according to the first heart rate array, a first heart rate function is fitted by using a nonlinear least square fitting algorithm; the first heart rate function is expressed as ; according to the first electromyogram amplitude array, a first electromyogram amplitude function is fitted by using a nonlinear least square fitting algorithm; the first electromyogram amplitude function is expressed as ; According to the second blood pressure array, a second blood pressure function is fitted by using a nonlinear least square fitting algorithm; the second blood pressure function is expressed as ; wherein represents an actual pulse width; according to the second heart rate array, a second heart rate function is fitted by using a nonlinear least square fitting algorithm; the second heart rate function is expressed as ; according to the second myoelectric signal amplitude array, a second myoelectric signal amplitude function is fitted by using a nonlinear least square fitting algorithm; the second myoelectric signal amplitude function is expressed as ; According to the third blood pressure array, a third blood pressure function is fitted by using a nonlinear least square fitting algorithm; the third blood pressure function is expressed as ; wherein represents an actual current intensity; according to the third heart rate array, a third heart rate function is fitted by using a nonlinear least square fitting algorithm; the third heart rate function is expressed as ; according to the third myoelectric signal amplitude array, a third myoelectric signal amplitude function is fitted by using a nonlinear least square fitting algorithm; the third myoelectric signal amplitude function is expressed as ; The blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude three-dimensional function are calculated through a three-dimensional construction formula 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; the three-dimensional construction formula is: ； wherein, represents a blood pressure value three-dimensional function, represents a heart rate value three-dimensional function, represents a myoelectric signal amplitude three-dimensional function, represents a reference stimulation frequency, represents a reference pulse width, represents a reference current intensity, represents a reference blood pressure value, represents a reference heart rate value, represents a reference myoelectric signal amplitude.

3. The low irritation electro-acupuncture nerve electro-stimulation parameter optimization method of claim 2, wherein, The electro-acupuncture needle nerve electric stimulation optimization parameter is calculated through a parameter optimization algorithm according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude three-dimensional function; the electro-acupuncture needle nerve electric stimulation optimization parameter comprises an optimal stimulation frequency, an optimal pulse width and an optimal current intensity; the method comprises: A first target function is constructed according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude three-dimensional function; A first constraint condition is constructed according to 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: ； wherein, represents a lower limit stimulation frequency, represents a lower limit pulse width, represents a lower limit current intensity, represents an upper limit stimulation frequency, represents an upper limit pulse width, represents an upper limit current intensity; The optimal stimulation frequency, the optimal pulse width and the optimal current intensity are obtained by solving through a particle swarm optimization algorithm with the first target function taking the minimum value as the target and the first constraint condition as the constraint condition.

4. The low irritation electro-acupuncture nerve electro-stimulation parameter optimization method of claim 3, wherein, The method for constructing the first target function according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the electromyographic signal amplitude three-dimensional function comprises: According to the blood pressure value three-dimensional function, the heart rate value three-dimensional function, the electromyogram signal amplitude value three-dimensional function and the pre-set lower limit blood pressure allowable value , the pre-set upper limit blood pressure allowable value , the pre-set lower limit heart rate allowable value , the pre-set upper limit heart rate allowable value , the pre-set lower limit electromyogram signal amplitude value allowable value , the pre-set upper limit electromyogram signal amplitude value allowable value , a first deviation function, a second deviation function and a third deviation function are constructed; the first deviation function, the second deviation function and the third deviation function are respectively: ； ； ； wherein denotes a first deviation function, denotes a second deviation function, denotes a third deviation function; constructing a first objective function according to the first deviation function, the second deviation function and the third deviation function ; the first objective function is: 。 5. A low irritation electro-acupuncture nerve electro-stimulation parameter optimization system characterized by, The system comprises: The data reading module is used for acquiring first data, and the first data includes preset neural electric stimulation reference parameters, neural electric stimulation lower limit parameters and neural electric stimulation upper limit parameters; the neural electric stimulation reference parameters include reference stimulation frequency, reference pulse width and reference current intensity; the neural electric stimulation lower limit parameters include lower limit stimulation frequency, lower limit pulse width and lower limit current intensity; and the neural electric stimulation upper limit parameters include upper limit stimulation frequency, upper limit pulse width and upper limit current intensity; The parameter acquisition module is connected with the data reading module and is used for inserting the electric thimble into an acupuncture point of a stimulation object, setting the actual stimulation frequency of the electric thimble as the reference stimulation frequency, the actual pulse width as the reference pulse width and the actual current intensity as the reference current intensity, and reading blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation object by the physiological parameter acquisition device to obtain reference blood pressure value, reference heart rate value and reference muscle electric signal amplitude; The function fitting module is connected with the parameter acquisition module and is used for setting the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric thimble at intervals of a preset adjustment interval time, and reading blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation object by the physiological parameter acquisition device to obtain blood pressure value three-dimensional function, heart rate value three-dimensional function and muscle electric signal amplitude three-dimensional function; The optimization calculation module is connected with the function fitting module and is used for calculating electric thimble neural electric stimulation optimization parameters including optimal stimulation frequency, optimal pulse width and optimal current intensity by a parameter optimization algorithm according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the muscle electric signal amplitude three-dimensional function; The function fitting module includes: The first measurement module is used for setting the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric thimble as the lower limit stimulation frequency, the reference pulse width and the reference current intensity respectively, adjusting the actual stimulation frequency at intervals of a preset adjustment interval time and at a preset first frequency as a step length, making the actual stimulation frequency increase step by step until the first difference is greater than the upper limit stimulation frequency minus the first frequency, and reading blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation object to obtain a first blood pressure array, a first heart rate array and a first muscle electric signal amplitude array; The second measurement module is connected with the first measurement module and is used for setting the actual stimulation frequency, the actual pulse width and the actual current intensity of the electric thimble as the reference stimulation frequency, the lower limit pulse width and the reference current intensity respectively, adjusting the actual pulse width at intervals of a preset adjustment interval time and at a preset first pulse width as a step length, making the actual pulse width increase step by step until the first difference is greater than the upper limit pulse width minus the first pulse width, and reading blood pressure value, heart rate value and muscle electric signal amplitude of the stimulation object to obtain a second blood pressure array, a second heart rate array and a second muscle electric signal amplitude array. The third measuring module is connected with the second measuring module, and is configured to set an actual stimulation frequency, an actual pulse width and an actual current intensity of the electric acupuncture needle as a reference stimulation frequency, a reference pulse width and a lower limit current intensity respectively; take a preset adjustment interval time as an interval, and adjust the actual current intensity by a preset first current intensity as a step to make the actual current intensity increase step by step until the actual current intensity is greater than a difference between an upper limit current intensity and the first current intensity for the first time; read a blood pressure value, a heart rate value and a muscle electrical signal amplitude value of the stimulation object to obtain a third blood pressure array, a third heart rate array and a third muscle electrical signal amplitude array; The fitting calculation module is connected with the third measuring module, and is configured to obtain a blood pressure value three-dimensional function, a heart rate value three-dimensional function and a muscle electrical signal amplitude three-dimensional function by a fitting algorithm according to the first blood pressure array, the first heart rate array, the first muscle electrical signal amplitude array, the second blood pressure array, the second heart rate array, the second muscle electrical signal amplitude array, the third blood pressure array, the third heart rate array and the third muscle electrical signal amplitude array.

6. The low irritation electro-acupressure nerve electro-stimulation parameter optimization system of claim 5, wherein, The fitting calculation module comprises: The first fitting module is configured to fit a first blood pressure function according to the first blood pressure array by using a nonlinear least square fitting algorithm, wherein the first blood pressure function is expressed as ; wherein represents an actual stimulation frequency; the second fitting module is configured to fit a first heart rate function according to the first heart rate array by using the nonlinear least square fitting algorithm, wherein the first heart rate function is expressed as ; the third fitting module is configured to fit a first electromyogram amplitude function according to the first electromyogram amplitude array by using the nonlinear least square fitting algorithm, wherein the first electromyogram amplitude function is expressed as ; The second fitting module is connected with the first fitting module, and is configured to obtain a second blood pressure function by using a nonlinear least square fitting algorithm according to the second blood pressure array; the second blood pressure function is expressed as ; wherein represents an actual pulse width; a second heart rate function is obtained by using a nonlinear least square fitting algorithm according to the second heart rate array; the second heart rate function is expressed as ; and a second electromyogram amplitude function is obtained by using a nonlinear least square fitting algorithm according to the second electromyogram amplitude array; the second electromyogram amplitude function is expressed as ; The third fitting module is connected with the second fitting module, and is configured to obtain a third blood pressure function by using a nonlinear least square fitting algorithm according to the third blood pressure array; the third blood pressure function is expressed as ; wherein represents an actual current intensity; a third heart rate function is obtained by using a nonlinear least square fitting algorithm according to the third heart rate array; the third heart rate function is expressed as ; a third electromyogram amplitude function is obtained by using a nonlinear least square fitting algorithm according to the third electromyogram amplitude array; the third electromyogram amplitude function is expressed as ; The three-dimensional construction module is connected with the third fitting module, and is configured to obtain the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the muscle electrical signal amplitude three-dimensional function by a three-dimensional construction formula 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 muscle electrical signal amplitude function, the second muscle electrical signal amplitude function and the third muscle electrical signal amplitude function; the three-dimensional construction formula is: ； wherein, represents a blood pressure value three-dimensional function, represents a heart rate value three-dimensional function, represents a myoelectric signal amplitude three-dimensional function, represents a reference stimulation frequency, represents a reference pulse width, represents a reference current intensity, represents a reference blood pressure value, represents a reference heart rate value, represents a reference myoelectric signal amplitude.

7. The low irritation electro-acupressure nerve electro-stimulation parameter optimization system of claim 6, wherein, The optimization calculation module comprises: The target function construction module is configured to construct a first target function according to the blood pressure value three-dimensional function, the heart rate value three-dimensional function and the muscle electrical signal amplitude three-dimensional function; The constraint condition construction module is connected with the target function construction module, and is configured to construct a first constraint condition according to 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: ； wherein, represents a lower limit stimulation frequency, represents a lower limit pulse width, represents a lower limit current intensity, represents an upper limit stimulation frequency, represents an upper limit pulse width, represents an upper limit current intensity; The particle swarm optimization module is connected with the constraint condition construction module, and is configured to take the first target function as a target and the first constraint condition as a constraint condition, and solve the optimal stimulation frequency, the optimal pulse width and the optimal current intensity by using a particle swarm optimization algorithm.

8. The low irritation electro-acupressure nerve electro-stimulation parameter optimization system of claim 7, wherein, The target function construction module comprises: The bias function generation module is configured to generate a first bias function, a second bias function, and a third bias function according to a blood pressure value three-dimensional function, a heart rate value three-dimensional function, a myoelectric signal amplitude three-dimensional function, a pre-set lower limit blood pressure allowable value, a pre-set upper limit blood pressure allowable value, a pre-set lower limit heart rate allowable value, a pre-set upper limit heart rate allowable value, a pre-set lower limit myoelectric signal amplitude allowable value, and a pre-set upper limit myoelectric signal amplitude allowable value. The first bias function, the second bias function, and the third bias function are respectively:​​​​​ ； ； ； wherein denotes a first deviation function, denotes a second deviation function, denotes a third deviation function; The target function generation module is connected with the bias function generation module, and is configured to construct a first target function according to the first bias function, the second bias function and the third bias function The first target function is: 。

Images