CN101816821B - Walking aid functional electrical stimulation precision control method based on ant colony fuzzy controller - Google Patents

Walking aid functional electrical stimulation precision control method based on ant colony fuzzy controller Download PDF

Info

Publication number
CN101816821B
CN101816821B CN2010101822151A CN201010182215A CN101816821B CN 101816821 B CN101816821 B CN 101816821B CN 2010101822151 A CN2010101822151 A CN 2010101822151A CN 201010182215 A CN201010182215 A CN 201010182215A CN 101816821 B CN101816821 B CN 101816821B
Authority
CN
China
Prior art keywords
fuzzy controller
kfuzzi
ant
domain
fuzzy
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN2010101822151A
Other languages
Chinese (zh)
Other versions
CN101816821A (en
Inventor
明东
张广举
邱爽
刘秀云
徐瑞
万柏坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
DATIAN MEDICAL SCIENCE ENGINEERING (TIANJIN) Co.,Ltd.
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN2010101822151A priority Critical patent/CN101816821B/en
Publication of CN101816821A publication Critical patent/CN101816821A/en
Application granted granted Critical
Publication of CN101816821B publication Critical patent/CN101816821B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to the field of rehabilitation devices and discloses a walking aid functional electrical stimulation precision control method based on an ant colony fuzzy controller, aiming at effectively improving the accuracy and the stability of an FES (Functional Electrical Stimulation) system. In the technical scheme of the invention, the walking aid FES precision control method based on the ant colony fuzzy controller comprises the steps of: firstly, converting a quantitative factor and a proportional factor of the fuzzy controller and the selection of 12 decision factors of a membership function parameter into a combination optimization problem applicable to an ant colony algorithm and carrying out encoding on the combination optimization problem and generating n initial urban agglomerations formed by individuals randomly; then, establishing a reasonable corresponding relationship target function of an actual joint angle and a muscle model output joint angle and determining the parameter configuration of the ant colony algorithm; entering an optimizing process; and regulating ant colony information quantity according to deviation, entering a next optimizing process, repeating the process, finally realizing the self-adaption on-line setting of the parameters of the fuzzy controller and applying to the FES system. The invention is mainly used for improving the accuracy and the stability of the FES system.

Description

Walk-aiding functional electric stimulation precision control method based on ant crowd fuzzy controller
Technical field
The present invention relates to the rehabilitation appliances field, especially based on the walk-aiding functional electric stimulation precision control method of ant crowd fuzzy controller.
Background technology
(Functional Electrical Stimulation is to stimulate limb motion muscle group and peripheral nerve thereof through current pulse sequence FES) to functional electrostimulation, recovers or rebuild the technology of the componental movement function of paralytic patient effectively.At present, because the spinal cord regeneration ability is faint, to the spinal cord injury paralysed patient, Shang Weiyou can directly repair effective treatment method of damage, and implementing function rehabilitation training is effective measures.According to statistics, spinal cord injury paralysed patient number increases year by year, and function rehabilitation training is a technology of demanding demand urgently.The sixties in 20th century, Liberson successfully utilizes the electro photoluminescence nervus peronaeus to correct the gait of hemiplegic patient's drop foot first, has started the new way that functional electrostimulation is used to move and Sensory rehabilitation is treated.At present, FES has become the componental movement function of recovering or rebuilding paralytic patient, is important rehabilitation means.Yet how accurate triggering sequential and the pulse current intensity of controlling FES can accurately be accomplished the key problem in technology that the intended function action is still FES with assurance electro photoluminescence action effect.According to statistics; The mode of the triggering of FES control is at present studied still few; And according to action effect and predetermined action deviation; Come adjustment FES stimulus intensity and time sequence parameter automatically with closed-loop control, thereby improved real-time, accuracy and the stability of FES system greatly, but present effectively control method is still among exploring.
Fuzzy controller is a kind of method through fuzzy logic and approximate resoning; People's experience formalization, modelling; Become computing machine acceptable controlling models; Let computer generation replace the people to control the high-level policy and novel technical method of controlled device in real time; Can improve controllability, adaptability and the rationality of control algolithm effectively, especially be difficult to modeling and have the problem of enriching manual experience to have peculiar advantage with math equation, and human muscle's complicacy and time variation operating environment make it set up its mathematical model to complicacy; Cause traditional control method to be difficult to adapt to the strict demand in FES field, fuzzy controller is that the precision control of FES provides new departure.The fuzzy controller core technology is exactly to confirm the factors such as method of structure, the fuzzy rule that is adopted, compositional rule of inference algorithm and the fuzzy decision of fuzzy controller, and the key that fuzzy control will obtain optimum control effect is to adjust to fuzzy controller quantizing factor, scale factor and fuzzy control rule are isoparametric.In the FES field, system stability is required very strictness, so select also particularly important to Fuzzy Controller Parameters.
Ant in the ant group algorithm simulation biological world is having no prompting to seek down by the foraging behavior proposition of the shortest path of ant cave to the food source simulated evolutionary algorithm based on population; Has stronger adaptability; Distributed parallel calculates, and is easy to the integrated advantage of other algorithms.At present, also not on top of, experimental knowledge relatively is short of whole mechanism of muscle.Utilize ant group algorithm adjusting of Fuzzy Controller Parameters to be helped the control of the precision of functional electrostimulation.
Summary of the invention
For overcoming the deficiency of prior art, purport of the present invention is the precision control method that proposes a kind of new FES, the accurately stable current-mode of controlling the FES system in real time.The present invention can improve FES system accuracy and stability effectively, and obtains considerable social benefit and economic benefit.For achieving the above object, the technical scheme that the present invention adopts is: the walk-aiding functional electric stimulation precision control method based on ant crowd fuzzy controller comprises the following steps:
At first the selection with 12 decision variables of quantizing factor, scale factor and the subordinate function parameter of fuzzy controller is converted into the combinatorial optimization problem that ant group algorithm is suitable for, and it is encoded and produces the initial group of cities that the n individuals is formed at random;
Next sets up rational actual joint angles and the corresponding relation objective function of muscle model output joint angles and the parameter setting of definite ant group algorithm;
Searching process: utilize the ant random search; Optimize subordinate function and the quantizing factor and the scale factor of fuzzy controller; And call the fuzzy controller of having adjusted, whether checking reaches goal-selling, does not repeat above operation repeatedly if having; Parameter convergence up to ant group algorithm perhaps reaches predetermined index, and final output promptly gets the decision variable of fuzzy controller and the number of times of ant crowd operation;
The decision variable of the fuzzy controller that the aforementioned final output of foundation promptly gets; Calculate the deviation of output and this output and muscle model output by fuzzy controller; According to deviation ant crowd quantity of information is adjusted, and got into next searching process, repeatedly this process; The final self-adaptation on-line tuning of realizing Fuzzy Controller Parameters, and be used for the FES system.
The said suitable combinatorial optimization problem of ant group algorithm that is converted into is to be converted into to find the solution shortest route problem, and concrete grammar is that the quantizing factor of ant crowd fuzzy controller and scale factor are that the basic quantization factor of fuzzy controller and basic scale factor multiply by factor kfuzzi (the i) (i=1 that is optimized by ant group algorithm respectively; 2,3), the adjustment of ant crowd fuzzy controller subordinate function promptly is adjusting to basic fuzzy controller subordinate function domain; On the basic domain basis of fuzzy controller, add or deduct factor kfuzzi (i) (i=4,5,6 that ant group algorithm is optimized; 7,8,9; 10,11,12); And to optimization factor kfuzzi (i) (i=1,2 of ant group algorithm ... The initial group of cities of 12) encoding and producing n individuals composition at random, the said binary coding that is encoded to.
Said binary coding, corresponding decoding formula does
kfuzzi ( i ) = x i min + ( Σ j = 0 l b j × 2 j - 1 ) × x i max - x i min 2 l - 1
Wherein kfuzzi (i) is the variable quantity of the fuzzy controller adjusted, and l is the length of coding, b ∈ [0,1], x ImaxAnd x IminBe respectively the maximal value and the minimum value of decision content.
Adopt following method to confirm described fuzzy controller:
Input fuzzy controller initialization module variable is respectively the error e (k) and the error change rate ec (k) of actual output joint angles and expectation joint degree, and its domain is FE=[E, E]; FEC=[EC, EC], the stimulating current intensity u (k) of output; Its domain is FU=[U, U]
The quantification domain of error does
X={-n,-n+1,…0,…,n-1,n} (1)
The quantification domain of error rate does
X 1={-m,-m+1,…0,…,m-1,m}; (2)
The quantification domain of controlled quentity controlled variable does
Y={-k,-k+1,…0,…,k-1,k} (3)
Quantizing factor is respectively
K e=n/X e (4)
K ec=m/X ec (5)
Scale factor does
K u=k/Y u (6)
The opinion domain of error: { 3-2-1 012 3}; The domain of error rate be the domain of 3-2-1 012 3} output valves 3-2-1 012 3}, the control law of fuzzy control is: if E1 and EC1 then U1, if E2 and EC2 then U2 ... Ep and ECp be Up then;
Its total fuzzy control rule is:
Figure GDA0000156290990000031
R=(E i×CE i) T1οC i(8)
E wherein 1=(a 1iA Ni), EC 1=(b 1iB Mi), U 1=(c 1iC Ti) (i=1 ... P)
The reverse gelatinization method that adopts is a method of weighted mean:
u c = ( Σ i = - s s ik i ) / ( Σ i = - s s u i ) - - - ( 9 )
For each concrete observed value deviation E *With its error rate EC *, use quantizing factor formula separately to become the element that quantizes in the domain more respectively, again its fuzzy E that turns to *And EC *,
Figure GDA0000156290990000033
Figure GDA0000156290990000034
E wherein *=(e 1E n), EC *=(f 1F m)
Can be by formula 8 in the hope of the accurate amount of output;
The method of said definite fuzzy controller adopts angular error and error rate thereof on [90 90], and domain is [3 3], then can use formula
X = ( 2 n b - a ( x ′ - a + b 2 ) ) .
Said subordinate function and quantizing factor and the scale factor that utilizes the ant random search to make its variable optimize fuzzy controller be,
kfuzzi ( i ) = x i min + ( Σ j = 0 l b j × 2 j - 1 ) × x i max - x i min 2 l - 1 - - - ( 17 )
Wherein kfuzzi (i) is the variable quantity of the fuzzy controller adjusted, and l is the length of coding, b ∈ [0,1], x ImaxAnd x IminBe respectively the maximal value and the minimum value of decision content.
Ant crowd fuzzy controller is to utilize ant group algorithm to the optimization of adjusting of the relevant parameter of fuzzy controller, and concrete steps are following:
Ant group algorithm is to the control of Fuzzy Controller Parameters, like formula (18), (19), shown in (20):
K e1=kfuzzi(1)*K e (18)
K c1=kfuzzi(2)*K c (19)
K u1=kfuzzi(3)*K u (20)
The error domain of fuzzy controller does
{-3-kfuzzi(4),-2-kfuzzi(5),-1-kfuzzi(6),0,1+kfuzzi(6),2+kfuzzi(5),3+kfuzzi(4)}
The error rate domain of fuzzy controller does
{-3-kfuzzi(7),-2-kfuzzi(8),-1-kfuzzi(9),0,1+kfuzzi(9),2+kfuzzi(8),3+kfuzzi(7)}
The domain of the output valve of fuzzy controller
{-3-kfuzzi(10),-2-kfuzzi(11),-1-kfuzzi(12),0,1+kfuzzi(12),2+kfuzzi(11),3+kfuzzi(10)}。
Said subordinate function and quantizing factor and the scale factor course of work of utilizing the ant random search to make its variable optimize fuzzy controller is:
Step1: parameter initialization.Make time t=0 and cycle index N Max=0, maximum cycle N is set Cmax, m ant placed starting point.
Step2: ant number and cycle index are set
Step3: the ant random search, after the end of once creeping, determine the actual input variable of the selected conduct of which characteristic variable, revise the taboo list index, after promptly choosing ant is moved to new element, and move this element in the individual taboo table of this ant
Step4: calculate subordinate function, the probability that ant individual state transition probability formula calculates is selected element
Step5: utilize the condition of the information generating fuzzy rule that training sample provides, the accuracy of check fuzzy model
Step6:, change Step3, otherwise be Step7 if the ant element has not traveled through
Step7: the pheromone concentration that the plain concentration of lastest imformation is divided the high characteristic variable of accuracy is enhanced, and next time can be selected with bigger probability when searching for
Step8: satisfy and finish to regulate the end of adjusting.
Characteristics of the present invention are: utilize the ant group algorithm Fuzzy Controller Parameters of adjusting, and the accurately stable then strength of current of controlling the FES system in real time effectively, and can improve FES system real time, accuracy and stability effectively.
Description of drawings
Fig. 1 ant group algorithm structured flowchart of Fuzzy Controller Parameters of adjusting.
Fig. 2 ant group algorithm structural map.
Fig. 3 experiment scene figure.
The fuzzy controller that Fig. 4 ant group algorithm is adjusted is followed the trail of the result.
The adjust relative error of the preset down input joint angles of Fuzzy Controller Parameters control and actual output of Fig. 5 ant group algorithm.
Embodiment
Proposed to adjust fuzzy controller (Fuzzy Controller) parameter with accurate control function property electro photoluminescence (Functional Electrical Stimulation, FES) new method of current-mode by ant group algorithm (Ant Algorithm) self-adaptation.Its techniqueflow is: optimize subordinate function and the quantizing factor and the scale factor of fuzzy controller through the ant group algorithm self-adjusting, control the current-mode of FES system then.This method is a kind of brand-new functional electrostimulation accurate control technique.
Structure based on the application of the walk-aiding functional electric stimulation precision control method of ant crowd fuzzy controller is as shown in Figure 1.Its workflow is: at first the selection with 12 decision variables of quantizing factor, scale factor and the subordinate function parameter of fuzzy controller is converted into the combinatorial optimization problem that ant group algorithm is suitable for; And to its initial group of cities of encoding and producing n individuals composition at random; Next sets up rational actual joint angles and the corresponding relation objective function of muscle model output joint angles and the parameter setting of definite ant group algorithm; Utilize the ant random search to make its variable optimize subordinate function and the quantizing factor and the scale factor of fuzzy controller; And call the fuzzy controller of having adjusted; Whether checking reaches goal-selling, does not repeat above operation repeatedly if having, and perhaps reaches predetermined index up to parameter convergence; Final output promptly gets the decision variable of fuzzy controller and the number of times of ant crowd operation.Computing system output under the new fuzzy controller and with the deviation of muscle model after to the adjustment of ant crowd quantity of information, make it get into next searching process.This process finally realizes the self-adaptation on-line tuning of Fuzzy Controller Parameters repeatedly, and is used for the FES system.
1 design of fuzzy control
Because people's singularity, the FES field is strict to controller stability, robustness, real-time, designs fuzzy controller equalization stable property and real-time and has selected two-dimensional fuzzy controller; Promptly two input variable difference reality are exported joint angles and the error e (k) and the error change rate ec (k) that expect the joint degree, and its domain is FE=[E, E]; FEC=[EC, EC], the stimulating current intensity u (k) of output; Its domain is FU=[U, U].
The quantification domain of error does
X={-n,-n+1,…0,…,n-1,n} (1)
The quantification domain of error rate does
X 1={-m,-m+1,…0,…,m-1,m}; (2)
The quantification domain of controlled quentity controlled variable does
Y={-k,-k+1,…0,…,k-1,k} (3)
Quantizing factor is respectively
K e=n/X e (4)
K ec=m/X ec (5)
Scale factor does
K u=k/Y u (6)
The present invention adopts the opinion domain of error: { 3-2-1 012 3}; The domain of error rate is the { domain of 3-2-1 01 23} output valves { 3-2-1 012 3}.Control rule tables is: if E1 and EC1 then U1, if E2 and EC2 then U2 ... Ep and ECp be Up then;
Its total fuzzy control rule does
Figure GDA0000156290990000051
R=(E i×CE i) T1οC i(8)
E wherein 1=(a 1iA Ni), EC 1=(b 1iB Mi), U 1=(c 1iC Ti) (i=1 ... P)
The reverse gelatinization method that adopts is a method of weighted mean:
u c = ( Σ i = - s s ik i ) / ( Σ i = - s s u i ) - - - ( 9 )
For each concrete observed value deviation E *With its error rate EC *, use quantizing factor formula separately to become the element that quantizes in the domain more respectively, again its fuzzy E that turns to *And EC *,
Figure GDA0000156290990000062
Figure GDA0000156290990000063
E wherein *=(e 1E n), EC *=(f 1F m)
Have formula 8 can in the hope of output accurate amount.
Angular error of the present invention and error rate thereof are on [90 90], and domain is [3 3], then can use formula
X = ( 2 n b - a ( x ′ - a + b 2 ) ) - - - ( 11 )
The 2 ant group algorithms Fuzzy Controller Parameters of adjusting
Ant group algorithm is a kind of novel bionic Algorithm that comes from the Nature biological world, when finding the solution optimization problem with ant group algorithm, at first optimization problem is transformed in order to find the solution shortest route problem.Every ant is from initial contact N 00Or N 01Set out, N in proper order passes by 1, N 2, a wherein child node, up to destination node N K0, N K1Form path (N 0tN 1tN Kt), t ∈ [0,1].A binary feasible solution can be represented in its path.Following characteristic is arranged during each ant visit city:
The state transformation rule: the state transformation rule that ant group algorithm uses is the rule of ratio at random that proposes based on the TSP problem, and it provides the probability that the ant k that is positioned at city i selects to move to city j,
Figure GDA0000156290990000065
τ wherein Ij(i j) is (i, fitness j), η Ij(i j) is the inverse of distance.α is the relative significance level of residual risk, the relative significance level that β is expectation value.
In ant group algorithm, selection mode does
Figure GDA0000156290990000066
Wherein, q is for being evenly distributed on a random number on [0,1], q 0Be the parameter on [0,1].
Overall situation update rule: ant algorithm has different update algorithm; The overall situation that ant crowd system adopts is upgraded principle; Only allowing the ant release pheromone of globally optimal solution, is the neighborhood that mainly concentrates on the best path of being found out till the current circulation for the search that makes ant like this.
τ ij(i,j)←(1-ρ)□τ ij(i,j)+ρ·Δτ ij(i,j) (14)
Figure GDA0000156290990000071
Wherein ρ is that information is counted volatility coefficient, L GbBe the global optimum's path local updating information that finds so far: every ant is set up the renewal that the plain mark of the information of carrying out number is also arranged in the process of separating
τ ij(i,j)←(1-γ)□τ ij(i,j)+γ·Δτ ij(i,j) (16)
γ ∈ [0,1] wherein.
Basic fuzzy controller steady-state behaviour can not reach the requirement in FES field; It promptly is that correlation parameter with fuzzy controller utilizes binary coding to be converted into to find the solution shortest route problem that the present invention will adjust to quantizing factor and scale factor and fuzzy control rule, and the structure of ant group algorithm is as shown in Figure 2.When basic domain and word set were constant, the quantizing factor variation can cause deviation and the pairing language value of rate of change thereof to change, and the controlled quentity controlled variable that the variation of scale factor can directly cause acting on controlled device changes.Be specially: K eBig more, system's rise time is short more, otherwise long more; K EcBig more, the reaction of system is sensitiveer, otherwise blunt more, K uBig more, system's rise time is short more, but causes vibration easily, and K uToo smallly be prone to make the dynamic process of system elongated.
What the present invention adopted is binary coding, and then corresponding decoding formula does
kfuzzi ( i ) = x i min + ( Σ j = 0 l b j × 2 j - 1 ) × x i max - x i min 2 l - 1 - - - ( 17 )
Wherein kfuzzi (i) is the variable quantity of the fuzzy controller adjusted, and l is the length of coding, b ∈ [0,1], x ImaxAnd x IminBe respectively the maximal value and the minimum value of decision content.
Ant group algorithm is to the control such as the following formula of Fuzzy Controller Parameters:
K e1=kfuzzi(1)*K e (18)
K c1=kfuzzi(2)*K c (19)
K u1=kfuzzi(3)*K u (20)
The error domain does
{-3-kfuzzi(4),-2-kfuzzi(5),-1-kfuzzi(6),0,1+kfuzzi(6),2+kfuzzi(5),3+kfuzzi(4)}
The error rate domain does
{-3-kfuzzi(7),-2-kfuzzi(8),-1-kfuzzi(9),0,1+kfuzzi(9),2+kfuzzi(8),3+kfuzzi(7)}
The domain of output valve
{-3-kfuzzi(10),-2-kfuzzi(11),-1-kfuzzi(12),0,1+kfuzzi(12),2+kfuzzi(11),3+kfuzzi(10)}
The adjust concrete workflow of Fuzzy Controller Parameters of ant group algorithm is:
Step1: parameter initialization.Make time t=0 and cycle index N Max=0, maximum cycle N is set Cmax, m ant placed starting point.
Step2: ant number and cycle index are set
Step3: the ant random search, after the end of once creeping, determine the actual input variable of the selected conduct of which characteristic variable, revise the taboo list index, after promptly choosing ant is moved to new element, and move this element in the individual taboo table of this ant
Step4: calculate subordinate function, the probability that ant individual state transition probability formula calculates is selected element
Step5: utilize the condition of the information generating fuzzy rule that training sample provides, the accuracy of check fuzzy model
Step6:, change Step3, otherwise be Step7 if the ant element has not traveled through
Step7: the pheromone concentration that the plain concentration of lastest imformation is divided the high characteristic variable of accuracy is enhanced, and next time can be selected with bigger probability when searching for
Step8: satisfy and finish to regulate the end of adjusting.
3 experimental programs
Parastep functional electrostimulation walk help system and the PS-2137 of PASCO company protractor and Data Studio software that experimental provision adopts U.S. SIGMEDICS company to produce.The Parastep system comprises microprocessor and boost pulse generation circuit, contains six stimulation channels, powered battery.Experiment content is: utilize the FES system that the relevant muscle group of lower limb is stimulated, utilize the PS-2137 of PASCO company protractor to gather knee joint angle and the measured knee joint angle of Data Studio software records.Require the experimenter healthy, no lower limb muscles, bone illness, impassivity illness and severe cardiac pulmonary disease.The experimenter sits idly on test board during experiment, and stimulating electrode is fixed in the end positions of quadriceps muscle of thigh, and protractor is fixed on thigh and the shank, makes the joint motion point press close to knee joint moving point position.Shank does not loosen, keeps vertical vacant state when applying electro photoluminescence, and the FES experiment scene is as shown in Figure 3.The electric stimulation pulse sequence adopts classical Lilly waveform, and pulsed frequency is 25Hz, pulsewidth 150 μ s, and pulse current is adjustable in 0~120m scope.Can adjust stimulus intensity to change the knee joint angle that produces by stimulating through changing the pulse current size in the experiment.Before the experiment, set the knee joint angle movement locus of expectation, utilize the measurement of angle meter to detect the knee joint subtended angle in real time in the experiment and change.The experimental data sampling rate is 128Hz, and the data recording duration is 60s.
The Fuzzy Controller Parameters new algorithm that the ant crowd adjusts is calculated the FES pulse current amplitude and is adjusted, the knee joint angle that the FES effect is produced move the movement locus of expection.Fig. 4 follows the trail of the result for the fuzzy controller that the ant group algorithm adaptive optimization is adjusted.Red line representes that desired movement track, blue line are actual output joint angles among the figure.The X axle is the time, and the Y axle is the motion of knee joint angle.For more clearly observing the departure that ant group algorithm is adjusted fuzzy controller; Shown in the relative error of preset input knee joint angle and actual knee joint angle under Fig. 5 ant group algorithm Tuning PID Controller; Can find out that then error all within 3%, can reach accurate control.
Purport of the present invention is the precision control method that proposes a kind of new FES, through the ant group algorithm self-adaptation Fuzzy Controller Parameters of adjusting, the accurately stable then strength of current of controlling the FES system in real time effectively.This invention can improve FES system real time, accuracy and stability effectively, and obtains considerable social benefit and economic benefit.Optimum implementation is intended and is adopted patent transfer, technological cooperation or product development.

Claims (4)

1. the walk-aiding functional electric stimulation precision control method based on ant crowd fuzzy controller is characterized in that, comprises the following steps:
At first the selection with 12 decision variables of quantizing factor, scale factor and the subordinate function parameter of fuzzy controller is converted into the combinatorial optimization problem that ant group algorithm is suitable for; And to its initial group of cities of encoding and producing n individuals composition at random; What said coding adopted is binary coding, and corresponding decoding formula is:
Figure FDA0000156290980000011
Wherein kfuzzi (i) is the variable quantity of the fuzzy controller adjusted, and l is the length of coding, b ∈ [0,1], x ImaxAnd x IminBe respectively the maximal value and the minimum value of decision content,
Ant crowd fuzzy controller promptly is to utilize ant group algorithm to the optimization of adjusting of the relevant parameter of fuzzy controller, and concrete steps are following:
Ant group algorithm is to the control of Fuzzy Controller Parameters, like formula (18), (19), shown in (20):
K e1=kfuzzi(1)*K e (18)
K c1=kfuzzi(2)*K c (19)
K u1=kfuzzi(3)*K u (20)
The error domain of fuzzy controller is:
{-3-kfuzzi(4),-2-kfuzzi(5),-1-kfuzzi(6),0,1+kfuzzi(6),2+kfuzzi(5),3+kfuzzi(4)}
The error rate domain of fuzzy controller is:
{-3-kfuzzi(7),-2-kfuzzi(8),-1-kfuzzi(9),0,1+kfuzzi(9),2+kfuzzi(8),3+kfuzzi(7)}
The domain of the output valve of fuzzy controller:
{-3-kfuzzi(10),-2-kfuzzi(11),-1-kfuzzi(12),0,1+kfuzzi(12),2+kfuzzi(11),3+kfuzzi(10)};
Next sets up rational actual joint angles and the corresponding relation objective function of muscle model output joint angles and the parameter setting of definite ant group algorithm;
Searching process: utilize the ant random search; Optimize subordinate function and the quantizing factor and the scale factor of fuzzy controller; And call the fuzzy controller of having adjusted, whether checking reaches goal-selling, does not repeat above operation repeatedly if having; Parameter convergence up to ant group algorithm perhaps reaches predetermined index, and final output promptly gets the decision variable of fuzzy controller and the number of times of ant crowd operation;
The decision variable of the fuzzy controller that the aforementioned final output of foundation promptly gets; Calculate the deviation of output and this output and muscle model output by fuzzy controller; According to deviation ant crowd quantity of information is adjusted, and got into next searching process, repeatedly this process; The final self-adaptation on-line tuning of realizing Fuzzy Controller Parameters, and be used for the FES system.
2. a kind of walk-aiding functional electric stimulation precision control method based on ant crowd fuzzy controller according to claim 1 is characterized in that, the said suitable combinatorial optimization problem of ant group algorithm that is converted into is to be converted into to find the solution shortest route problem; Concrete grammar is that the quantizing factor of ant crowd fuzzy controller and scale factor are that the basic quantization factor of fuzzy controller and basic scale factor multiply by the factor kfuzzi (i) that is optimized by ant group algorithm, i=1,2 respectively; 3, the adjustment of ant crowd fuzzy controller subordinate function promptly is the adjustment to the subordinate function domain of basic fuzzy controller, on the basic domain basis of fuzzy controller, adds or deduct the factor kfuzzi (i) that ant group algorithm is optimized, i=4; 5,6,7,8; 9,10,11; 12, and to the optimization factor kfuzzi (i) of ant group algorithm, i=1; 2 ... 12, the initial group of cities of encoding and producing n individuals composition at random, the said binary coding that is encoded to.
3. a kind of walk-aiding functional electric stimulation precision control method based on ant crowd fuzzy controller according to claim 1 is characterized in that, adopts following method to confirm described fuzzy controller:
Input fuzzy controller initialization module variable is respectively the error e (k) and the error change rate ec (k) of actual output joint angles and expectation joint degree, and its domain is FE=[E, E]; FEC=[EC, EC], the stimulating current intensity u (k) of output; Its domain is FU=[U, U]
The quantification domain of error does
X={-n,-n+1,…0,…,n-1,n} (1)
The quantification domain of error rate does
X 1={-m,-m+1,…0,…,m-1,m}; (2)
The quantification domain of controlled quentity controlled variable does
Y={-k,-k+1,…0,…,k-1,k} (3)
Quantizing factor is respectively
K e=n/X e (4)
K ec=m/X ec (5)
Scale factor does
K u=k/Y u (6)
Adopt the domain of error: { 3-2-1 012 3}; The domain of error rate be the domain of 3-2-1 012 3} output valves 3-2-1 012 3}, control rule tables is: if E1 and EC1 then U1, if E2 and EC2 then U2 ... Ep and ECp be Up then;
Its total fuzzy control rule is:
Figure FDA0000156290980000021
R=(E i×EC i) T1οC i?(8)
E wherein 1=(a 1iA Ni), EC 1=(b 1iB Mi), U 1=(c 1iC Ti) (i=1 ... P)
The reverse gelatinization method that adopts is a method of weighted mean:
Figure FDA0000156290980000022
For each concrete observed value deviation E *With its error rate EC *, use quantizing factor formula separately to become the element that quantizes in the domain more respectively, again its fuzzy E that turns to *And EC *,
Figure FDA0000156290980000023
Figure FDA0000156290980000024
E wherein *=(e 1E n), EC *=(f 1F m)
Ask the accurate amount of output by formula (9);
The method of said definite fuzzy controller adopts angular error and error rate thereof on [90 90], and domain is [3 3], then uses formula
Figure FDA0000156290980000031
transforms, and angular error and error rate commentaries on classics thereof are calculated its domain space.
4. a kind of walk-aiding functional electric stimulation precision control method according to claim 1 based on ant crowd fuzzy controller; It is characterized in that said subordinate function and quantizing factor and the scale factor course of work of utilizing the ant random search to make its variable optimize fuzzy controller is:
Step1: parameter initialization makes time t=0 and cycle index N Max=0, maximum cycle N is set Cmax, m ant placed starting point;
Step2: ant number and cycle index are set;
Step3: ant random search; After the end of once creeping, determine the actual input variable of the selected conduct of which characteristic variable, revise the taboo list index; After promptly choosing ant is moved to new element, and move to this element in the individual taboo table of ant;
Step4: calculate subordinate function, the probability that ant individual state transition probability formula calculates is selected element;
Step5: utilize the condition of the information generating fuzzy rule that training sample provides, the accuracy of check fuzzy model;
Step6:, change Step3, otherwise be Step7 if the ant element has not traveled through;
Step7: the pheromone concentration that the plain concentration of lastest imformation is divided the high characteristic variable of accuracy is enhanced, and next time can be selected with bigger probability when searching for;
Step8: satisfy termination condition, the end of adjusting.
CN2010101822151A 2010-05-26 2010-05-26 Walking aid functional electrical stimulation precision control method based on ant colony fuzzy controller Active CN101816821B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101822151A CN101816821B (en) 2010-05-26 2010-05-26 Walking aid functional electrical stimulation precision control method based on ant colony fuzzy controller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101822151A CN101816821B (en) 2010-05-26 2010-05-26 Walking aid functional electrical stimulation precision control method based on ant colony fuzzy controller

Publications (2)

Publication Number Publication Date
CN101816821A CN101816821A (en) 2010-09-01
CN101816821B true CN101816821B (en) 2012-08-15

Family

ID=42652189

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101822151A Active CN101816821B (en) 2010-05-26 2010-05-26 Walking aid functional electrical stimulation precision control method based on ant colony fuzzy controller

Country Status (1)

Country Link
CN (1) CN101816821B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102274581B (en) * 2011-05-18 2013-07-10 天津大学 Precise control method for functional electric stimulation
CN102284137A (en) * 2011-05-20 2011-12-21 天津大学 Functional electrical stimulation multisource information fusion control method
CN102662317B (en) * 2012-03-27 2014-05-21 中国人民解放军国防科学技术大学 PID controller based on prokaryotic bionic array
CN108363292B (en) * 2018-02-11 2020-11-03 广东电网有限责任公司电力科学研究院 Thermal power generating unit AGC control method and device based on fuzzy control strategy
CN109709958B (en) * 2018-12-26 2020-10-27 南京航空航天大学 Extensible control method for AGV electromagnetic navigation control system
CN113398467B (en) * 2021-06-02 2024-03-29 天津市第一中心医院 Fuzzy automatic adjustment method for dynamic range parameters of artificial cochlea

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6171239B1 (en) * 1998-08-17 2001-01-09 Emory University Systems, methods, and devices for controlling external devices by signals derived directly from the nervous system
CN1598719A (en) * 2004-09-27 2005-03-23 北京交通大学 Nerve network optimization controller and control method
CN101596338A (en) * 2009-04-29 2009-12-09 天津大学 Functional electric stimulation precision control method based on BP neural network tuned proportion integration differentiation PID
CN201431704Y (en) * 2009-06-10 2010-03-31 燕铁斌 Functional low-frequency electrical stimulation therapy device based on mode of human walking

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6171239B1 (en) * 1998-08-17 2001-01-09 Emory University Systems, methods, and devices for controlling external devices by signals derived directly from the nervous system
CN1598719A (en) * 2004-09-27 2005-03-23 北京交通大学 Nerve network optimization controller and control method
CN101596338A (en) * 2009-04-29 2009-12-09 天津大学 Functional electric stimulation precision control method based on BP neural network tuned proportion integration differentiation PID
CN201431704Y (en) * 2009-06-10 2010-03-31 燕铁斌 Functional low-frequency electrical stimulation therapy device based on mode of human walking

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
宋申民等.基于蚁群算法的模糊控制规则的过滤简化.《计算机仿真》.2006,第23卷(第3期),157-163. *

Also Published As

Publication number Publication date
CN101816821A (en) 2010-09-01

Similar Documents

Publication Publication Date Title
CN101816821B (en) Walking aid functional electrical stimulation precision control method based on ant colony fuzzy controller
CN101816822B (en) Setting method of functional electrical stimulation PID (Proportion Integration Differentiation) parameter double source characteristic fusion particle swarm
CN101596338A (en) Functional electric stimulation precision control method based on BP neural network tuned proportion integration differentiation PID
CN101794114B (en) Method for tuning control parameter in walk-aiding functional electric stimulation system by utilizing genetic algorithm
CN102274581B (en) Precise control method for functional electric stimulation
CN101837165B (en) Walking aid electrostimulation fine control method based on genetic-ant colony fusion fuzzy controller
CN101846977B (en) Genetic fuzzy control method of joint angles by functional electrical stimulation
CN105615890A (en) Angle and myoelectricity continuous decoding method for human body lower limb walking joint
CN105727442A (en) Closed-loop brain controlled functional electrical stimulation system
CN104680025B (en) Oil pumper parameter optimization method based on genetic algorithm extreme learning machine
CN102764167A (en) Myoelectric prosthesis control source lead optimization method based on correlation coefficients
Li et al. Neural network based modeling and control of elbow joint motion under functional electrical stimulation
Di Febbo et al. Does Reinforcement Learning outperform PID in the control of FES-induced elbow flex-extension?
CN102488964A (en) Functional electro stimulation closed loop fuzzy proportional integral derivative (PID) control method
Qiu et al. Intelligent algorithm tuning PID method of function electrical stimulation using knee joint angle
Zheng et al. A computational model for ratbot locomotion based on cyborg intelligence
CN101837164B (en) Double source feature fusion ant colony tuning method for PID (Proportion Integration Differention) parameter in functional electro-stimulation
CN102521508B (en) Adaptive neural fuzzy muscle modeling method under functional electrical stimulation
CN104680023B (en) Oil pumper parameter optimization method based on multiobjective decision-making
Girard et al. Robust discrete-time linear control of heart rate during treadmill exercise
Koelewijn et al. Predictive simulations to replicate human gait adaptations and energetics with exoskeletons
KR101593637B1 (en) Modeling system for muscle cell activation and method thereof
Barbosa et al. Control techniques for neuromuscular electrical stimulation: A brief survey
CN106156842A (en) A kind of integration provides the method for optimally controlling of neuron models
Hmed et al. Model free control for muscular force by functional electrical stimulation using pulse width modulation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210225

Address after: Room 101, building C22, entrepreneurship headquarters base, North Fuyuan Road, Wuqing Development Zone, Wuqing District, Tianjin

Patentee after: DATIAN MEDICAL SCIENCE ENGINEERING (TIANJIN) Co.,Ltd.

Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92

Patentee before: Tianjin University