CN104330968A - Inverse model PID compound control method based on improved support vector regression - Google Patents
Inverse model PID compound control method based on improved support vector regression Download PDFInfo
- Publication number
- CN104330968A CN104330968A CN201410668815.7A CN201410668815A CN104330968A CN 104330968 A CN104330968 A CN 104330968A CN 201410668815 A CN201410668815 A CN 201410668815A CN 104330968 A CN104330968 A CN 104330968A
- Authority
- CN
- China
- Prior art keywords
- support vector
- vector regression
- inversion model
- particle
- svr
- 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.)
- Pending
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
The invention requests for protecting an inverse model PID compound control method based on improved support vector regression. The method utilizes a multi-agent particle swarm optimization algorithm to optimize related parameters of a support vector regression to effectively improve modeling precision and generalization ability. On the above basis, an inverse model of MAPSO-SVR (Multi-Agent Particle Swarm Optimization-Support Vector Regression) is established, so that precision of the inverse model is improved to a certain extent; PID control is introduced, the MAPSO-SVR inverse model/PID compound control method is provided, own continuous amendment of a nonlinear system is achieved, and tracking and control ability of the system can be effectively improved.
Description
Technical field
The invention belongs to a kind of inverse model control method based on support vector regression, particularly relate to a kind of method optimizing the inversion model/PID complex controll nonlinear system of SVR correlation parameter based on multi-agent particle swarm (MAPSO).
Background technology
Method of inverse is set up than more complete design theory in the nonlinear system of general type in recent years, can nonlinear system in general form as research object, and special requirement is not had to equation form, do not need to introduce differential, several how abstract mathematical theory, thus there is general Research Significance yet.Have at present and much utilize the method for neural network inversion model extensively to be studied, also achieve theoretical and real achievement preferably.But in modeling, neural network exists Local Minimum, network number of plies difficulty such as to determine at the problem.And utilize support vector regression (SVR) to set up inversion model, because SVR has the advantages such as global optimum, generalization ability be strong, the antijamming capability of model effectively can be improved.The tracking performance of the direct influential system of degree of accuracy of model, the optimum choice of parameter affects significantly the precision of prediction of SVR and generalization ability, and common optimized algorithm has a lot, as genetic algorithm (GA), particle cluster algorithm (PSO) etc.Genetic algorithm is easily absorbed in local minimum, and this algorithm has very strong dependence to choosing of inherent parameters, particle cluster algorithm is only shared the information of optimal particle in colony and have ignored the information of other particle, particle may be caused just to attracted in early days near certain point at algorithm, thus reduce the diversity of population.
The present invention proposes the method utilizing multi-Agent particle cluster algorithm (MAPSO) to optimize SVR correlation parameter.MAPSO algorithm combines multiparticle information sharing and adaptable inertia weigh strategy, multiparticle information sharing adopts multiparticle information to revise each particle action strategy next time, to reduce the possibility that particle is absorbed in local optimum, inertia weight self-adaptative adjustment is according to colony's prematurity convergence degree, by individual concrete condition adjustment inertia weight, particle is made to jump out local optimum, by realizing global optimum with the competition of its neighbor particle, cooperation.Because method of inverse is " pseudo-linear system " structure in theory, exist and control deficiency, the problem that stability is not high, introduce PID FEEDBACK CONTROL, improve the tracking of system, control ability.
Summary of the invention
For above deficiency of the prior art, the object of the present invention is to provide a kind of improve system tracking, control ability based on the inversion model PID composite control method improving support vector regression, technical scheme of the present invention is as follows:
Based on the inversion model PID composite control method improving support vector regression, it comprises the following steps:
101, Initialize installation groupy phase related parameter, generating training data collection, test data set, and training dataset and test data set are normalized;
102, utilize multi-Agent particle cluster algorithm/MAPSO to carry out optimizing to support vector regression SVR correlation parameter, by fitness function and particle position, speed more new formula to support vector regression SVR parameter optimization;
103, judge whether to meet end condition, if meet, optimizing terminates thus determines the parameter optimization result of optimal particle and MAPSO-SVR, jumps to step 104; Otherwise return step 102;
104, determine the inversion model matching factor, set up MAPSO-SVR inversion model;
105, reference input function y is set
routthe u value in (k) and before moment, and according to y
routk (), inversion model and object model, calculate the value uctr (k) predicting out through inversion model;
106, by known y
rout(k), y
outk () value, according to definition and increment type PID principle △ u (the k)=k of error
p* x1+k
d* x2+k
i* x3, wherein k
prepresent proportional control factor, k
drepresent derivative control coefficient, k
irepresent integral control coefficient, x1=e (k)-e (k-1), x2=e (k)-2*e (k-1)+e (k-2), x3=e (k), and according to the uctr (k) in step 105, try to achieve current time u (k) value according to formula u (k)=uctr (k)+△ u (k);
107, subsequent time real output value y is calculated according to step 106 and object model
out(k+1);
If the value of 108 k does not also exceed the sampling time, then proceed to step 105, terminate until k value reaches sampling time setting value.
Further, the initialization groupy phase related parameter in step 101 comprises definition solution space, environment scale, maximum permission iterations, Inertia Weight scope, Studying factors.
Further, the particle fitness function F in step 102 is:
Wherein y
i, y '
irepresent SVR respectively and train output valve and desired output, N is number of training,
Further, the position of the information updating particle according to optimal particle in step 102, the formula of speed are:
In formula, subscript t is iterations, x
tparticles spatial position when being the t time iteration, v
tbe particle rapidity during the t time iteration, wherein C
1, C
2for Studying factors, r
1, r
2for the random number between (0,1), w is inertia constant.
Further, the inversion model matching factor in step 104 is
Advantage of the present invention and beneficial effect as follows:
The present invention utilizes MAPSO algorithm optimization SVR correlation parameter, introduces PID simultaneously and controls to improve whole system static and dynamic performance.Set up the not high and control system of inversion model precision for SVR and control the problems such as not enough, utilize MAPSO algorithm optimization SVR parameter, avoid the randomness of artificial selection parameter, improve the modeling accuracy that SVR sets up inversion model; Introduce PID to control, compensate for the problem of " pseudo-linear " Systematical control deficiency, the system that realizes constantly revises real output value makes input and output error little as much as possible, thus improves tracking, the control ability of system.The present invention carries out optimization, the improvement of control method mainly for the parameter optimization of SVR and the deficiency of Inverted control system, and its effect of optimization is obvious, improves the tracking of system, control ability.
Accompanying drawing explanation
Fig. 1 is according to multiplex control system program flow diagram of the present invention.
Embodiment
The invention will be further elaborated to provide an infinite embodiment below in conjunction with accompanying drawing.But should be appreciated that, these describe just example, and do not really want to limit the scope of the invention.In addition, in the following description, the description to known features and technology is eliminated, to avoid unnecessarily obscuring concept of the present invention.
As shown in Figure 1, the present invention uses inverse system principle, on the basis of inversion model, introduces MAPSO algorithm and PID control method, and improve modeling accuracy and improve system control performance, making the tracking of system, control ability gets a promotion, specific implementation step is as follows:
Step 1: Initialize installation correlation parameter, generating training data, test data data are normalized.
Step 2:MAPSO initialization of population, optimum configurations also produces initial position and the speed of population.
Step 3: according to step 2 and fitness function:
calculate current particle adaptive value, wherein y
i, y '
irepresent SVR respectively and train output valve and desired output, N is number of training.
Step 4: the personal best particle according to neighbor information more new particle: each particle predefine neighbours' environment, hypothetical particle L
i,j, and L
i,j=(l
1, l
2... l
d) be its position in optimization problem solution space, M
i,jl
i,jm neighbours in have the particle of optimal adaptation value, and M
i,j=(m
1, m
2..., m
d).If L
i,jmeet f (L
i,j)≤f (M
i,j), then it is a winner, otherwise it is a loser.If L
i,jbe a winner, it remains unchanged in the position of solution space; If it is loser, then a L
i,jin the position of solution space then according to L
i,j=M
i,j+ rand (-1,1) (L
i,j-M
i,j) upgrade.
Step 5: according to step 4 (i.e. P
best, g
best), and formula:
The more position of new particle, speed, in formula, subscript t is iterations, x
tparticles spatial position when being the t time iteration, v
tbe particle rapidity during the t time iteration, wherein C
1, C
2for Studying factors, r
1, r
2for the random number between (0,1), w is inertia constant.
Step 6: judge whether to meet end condition.If meet, optimizing terminates thus determines the parameter optimization result of optimal particle and MAPSO-SVR; Otherwise go to step 3.
Step 7: according to the inversion model modeling principle determination matching factor:
By object model
{ X, Y} data set, wherein the input control signal of object model is the random signal between 0 to 1, is trained, set up inversion model by SVR to it in the inputoutput data foundation generated.
Step 8: setting expectation function y
rout(k), and according to y
routk (), inversion model and object model draw uctr (k) value.
Step 9: according to definition e (the k)=y of error
rout(k)-y
out(k) and increment type PID principle △ u (k)=k
p* x1+k
d* x2+k
i* x3, wherein k
pexpression ratio controls, k
drepresent that differential controls, k
irepresent integration control, x1=e (k)-e (k-1), x2=e (k)-2*e (k-1)+e (k-2), x3=e (k).And integrating step 8, try to achieve current time u (k) value by u (k)=uctr (k)+△ u (k).
Step 10: according to step 9, the y in object model and before moment
outvalue, obtains subsequent time y
out(k+1) value.
Step 11: if the value of k does not also exceed the sampling time, then proceed to step 8, terminates until k value reaches sampling time setting value.
The control research of nonlinear system, in the field such as industrial or agricultural, manufacturing industry, has a wide range of applications.This method proposes to optimize SVR parameter with MAPSO and promotes inversion model modeling accuracy, introduces this conventional control of PID as feedback fraction, improves the tracking of nonlinear system, control ability, the performance of nonlinear control system is further optimized.
These embodiments are interpreted as only being not used in for illustration of the present invention limiting the scope of the invention above.After the content of reading record of the present invention, technician can make various changes or modifications the present invention, and these equivalence changes and modification fall into the scope of the claims in the present invention equally.
Claims (5)
1., based on the inversion model PID composite control method improving support vector regression, it is characterized in that: comprise the following steps:
101, Initialize installation groupy phase related parameter, generating training data collection, test data set, and training dataset and test data set are normalized;
102, utilize multi-Agent particle cluster algorithm/MAPSO to carry out optimizing to support vector regression SVR correlation parameter, by fitness function and particle position, speed more new formula to support vector regression SVR parameter optimization;
103, judge whether to meet end condition, if meet, optimizing terminates thus determines the parameter optimization result of optimal particle and MAPSO-SVR, jumps to step 104; Otherwise return step 102;
104, determine the inversion model matching factor, set up MAPSO-SVR inversion model;
105, reference input function y is set
routthe u value in (k) and before moment, and according to y
routk (), inversion model and object model, calculate the value uctr (k) predicting out through inversion model;
106, by known y
rout(k), y
outk () value, according to definition and increment type PID principle △ u (the k)=k of error
p* x1+k
d* x2+k
i* x3, wherein k
prepresent proportional control factor, k
drepresent derivative control coefficient, k
irepresent integral control coefficient, x1=e (k)-e (k-1), x2=e (k)-2*e (k-1)+e (k-2), x3=e (k), and according to the uctr (k) in step 105, try to achieve current time u (k) value according to formula u (k)=uctr (k)+△ u (k);
107, subsequent time real output value y is calculated according to step 106 and object model
out(k+1);
If the value of 108 k does not also exceed the sampling time, then proceed to step 105, terminate until k value reaches sampling time setting value.
2. the inversion model PID composite control method based on improving support vector regression according to claim 1, is characterized in that: the initialization groupy phase related parameter in step 101 comprises definition solution space, environment scale, maximum permission iterations, Inertia Weight scope, Studying factors.
3. the inversion model PID composite control method based on improving support vector regression according to claim 1, is characterized in that: the particle fitness function F in step 102 is:
Wherein y
i, y,
irepresent SVR respectively and train output valve and desired output, N is number of training.
4. the inversion model PID composite control method based on improving support vector regression according to claim 1, is characterized in that: the position of the information updating particle according to optimal particle in step 102, the formula of speed are:
In formula, subscript t is iterations, x
tparticles spatial position when being the t time iteration, v
tbe particle rapidity during the t time iteration, wherein C
1, C
2for Studying factors, r
1, r
2for the random number between (0,1), w is inertia constant.
5. the inversion model PID composite control method based on improving support vector regression according to claim 1, is characterized in that: the inversion model matching factor in step 104 is
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410668815.7A CN104330968A (en) | 2014-11-18 | 2014-11-18 | Inverse model PID compound control method based on improved support vector regression |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410668815.7A CN104330968A (en) | 2014-11-18 | 2014-11-18 | Inverse model PID compound control method based on improved support vector regression |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104330968A true CN104330968A (en) | 2015-02-04 |
Family
ID=52405715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410668815.7A Pending CN104330968A (en) | 2014-11-18 | 2014-11-18 | Inverse model PID compound control method based on improved support vector regression |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104330968A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105955032A (en) * | 2016-06-23 | 2016-09-21 | 上海电机学院 | Inverter control method for optimization of extreme learning machine on the basis of bat algorithm |
CN106527148A (en) * | 2016-12-15 | 2017-03-22 | 南京航空航天大学 | Sliding mode robust control method of discrete time delay uncertain system |
CN108233778A (en) * | 2018-01-19 | 2018-06-29 | 中国科学院合肥物质科学研究院 | A kind of more intelligent cell connection in series-parallel Zhan Shou mechanisms Collaborative Control optimization method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100252714B1 (en) * | 1997-10-16 | 2000-04-15 | 전주범 | Operating control device for virtual reality system |
CN103926930A (en) * | 2014-05-07 | 2014-07-16 | 重庆邮电大学 | Multi-robot cooperation map building method based on Hilbert curve detection |
-
2014
- 2014-11-18 CN CN201410668815.7A patent/CN104330968A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100252714B1 (en) * | 1997-10-16 | 2000-04-15 | 전주범 | Operating control device for virtual reality system |
CN103926930A (en) * | 2014-05-07 | 2014-07-16 | 重庆邮电大学 | Multi-robot cooperation map building method based on Hilbert curve detection |
Non-Patent Citations (8)
Title |
---|
唐贤伦等: "基于多Agent粒子群优化的多步SVR模型预测控制", 《***工程与电子技术》 * |
唐贤伦等: "多智能体粒子群优化的SVR模型预测控制", 《控制与决策》 * |
宋夫华等: "支持向量机α阶逆***控制-离散非线性***", 《浙江大学学报(工学版)》 * |
柯玉波等: "基于SVM的非线性***逆控制的设计与仿真", 《2007中国控制与决策学术年会论文集》 * |
穆朝絮等: "基于粒子群优化的非线性***最小二乘支持向量机预测控制方法", 《控制理论与应用》 * |
肖健梅等: "基于支持向量机的非线性***逆控制", 《中南大学学报(自然科学版)》 * |
袁小芳等: "基于支持向量机的非线性逆控制及仿真研究", 《湖南大学学报(自然科学版)》 * |
陈涛: "基于差分进化算法的支持向量回归机参数优化", 《计算机仿真》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105955032A (en) * | 2016-06-23 | 2016-09-21 | 上海电机学院 | Inverter control method for optimization of extreme learning machine on the basis of bat algorithm |
CN106527148A (en) * | 2016-12-15 | 2017-03-22 | 南京航空航天大学 | Sliding mode robust control method of discrete time delay uncertain system |
CN106527148B (en) * | 2016-12-15 | 2019-04-05 | 南京航空航天大学 | A kind of sliding Mode Robust Control method of Discrete-Delay uncertain system |
CN108233778A (en) * | 2018-01-19 | 2018-06-29 | 中国科学院合肥物质科学研究院 | A kind of more intelligent cell connection in series-parallel Zhan Shou mechanisms Collaborative Control optimization method |
CN108233778B (en) * | 2018-01-19 | 2020-05-12 | 中国科学院合肥物质科学研究院 | Multi-intelligent-unit series-parallel unfolding and folding mechanism cooperative control optimization method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Panda et al. | Hybrid BFOA–PSO algorithm for automatic generation control of linear and nonlinear interconnected power systems | |
Baillieul et al. | Encyclopedia of systems and control | |
Mohanty et al. | Differential evolution algorithm based automatic generation control for interconnected power systems with non-linearity | |
Lu et al. | A kernel extreme learning machine algorithm based on improved particle swam optimization | |
Sayed et al. | Gain tuning PI controllers for boiler turbine unit using a new hybrid jump PSO | |
CN101315544A (en) | Greenhouse intelligent control method | |
Kalyan et al. | Higher order degree of freedom controller for load frequency control of multi area interconnected power system with time delays | |
Gao et al. | Forecasting annual natural gas consumption via the application of a novel hybrid model | |
CN114217524A (en) | Power grid real-time self-adaptive decision-making method based on deep reinforcement learning | |
Huang et al. | Optimal design of load frequency active disturbance rejection control via double-chains quantum genetic algorithm | |
CN109254530A (en) | MFA control method based on grinding process basis circuit | |
Li et al. | A free search krill herd algorithm for functions optimization | |
CN104050505A (en) | Multilayer-perceptron training method based on bee colony algorithm with learning factor | |
CN104330968A (en) | Inverse model PID compound control method based on improved support vector regression | |
CN104616072A (en) | Method for improving concentration of glutamic acid fermented product based on interval optimization | |
Liu et al. | Extreme learning machine based on particle swarm optimization for estimation of reference evapotranspiration | |
Chen et al. | Sliding mode control based on disturbance observer for greenhouse climate systems | |
CN103559541A (en) | Back propagation method for out-of-order data stream in big data | |
Wang et al. | A hybrid model with combined feature selection based on optimized VMD and improved multi-objective coati optimization algorithm for short-term wind power prediction | |
CN107168066A (en) | A kind of greenhouse self-adaptation control method | |
Minqiang et al. | Temperature prediction of hydrogen producing reactor using SVM regression with PSO | |
CN103475608A (en) | Simulation annealing and fruit fly hybrid optimization wavelet generalized discrete multi-modulus blind equalization method | |
Gurban et al. | Greenhouse climate control enhancement by using genetic algorithms | |
Castañón | Dynamic noncooperative games | |
Kawan | Data rate of nonlinear control systems and feedback entropy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150204 |
|
RJ01 | Rejection of invention patent application after publication |