1、电气工程及其自动化毕业设计英语翻译郑州航空工业管理学院英 文 翻 译 2011 届 电气工程及其自动化 专业 班级题 目 遗传算法在非线性模型中的应用 姓 名 学号 指导教师 黄文力 职称 副教授 二一 一 年 三 月 三十 日英语原文:Application of Genetic Programming to Nonlinear ModelingIntroductionIdentification of nonlinear models which are based in part at least on the underlying physics of the real system
2、presents many problems since both the structure and parameters of the model may need to be determined. Many methods exist for the estimation of parameters from measures response data but structural identification is more difficult. Often a trial and error approach involving a combination of expert k
3、nowledge and experimental investigation is adopted to choose between a number of candidate models. Possible structures are deduced from engineering knowledge of the system and the parameters of these models are estimated from available experimental data. This procedure is time consuming and sub-opti
4、mal. Automation of this process would mean that a much larger range of potential model structure could be investigated more quickly.Genetic programming (GP) is an optimization method which can be used to optimize the nonlinear structure of a dynamic system by automatically selecting model structure
5、elements from a database and combining them optimally to form a complete mathematical model. Genetic programming works by emulating natural evolution to generate a model structure that maximizes (or minimizes) some objective function involving an appropriate measure of the level of agreement between
6、 the model and system response. A population of model structures evolves through many generations towards a solution using certain evolutionary operators and a “survival-of-the-fittest” selection scheme. The parameters of these models may be estimated in a separate and more conventional phase of the
7、 complete identification process.ApplicationGenetic programming is an established technique which has been applied to several nonlinear modeling tasks including the development of signal processing algorithms and the identification of chemical processes. In the identification of continuous time syst
8、em models, the application of a block diagram oriented simulation approach to GP optimization is discussed by Marenbach, Bettenhausen and Gray, and the issues involved in the application of GP to nonlinear system identification are discussed in Grays another paper. In this paper, Genetic programming
9、 is applied to the identification of model structures from experimental data. The systems under investigation are to be represented as nonlinear time domain continuous dynamic models.The model structure evolves as the GP algorithm minimizes some objective function involving an appropriate measure of
10、 the level of agreement between the model and system responses. One examples is (1) Where is the error between model output and experimental data for each of N data points. The GP algorithm constructs and reconstructs model structures from the function library. Simplex and simulated annealing method
11、 and the fitness of that model is evaluated using a fitness function such as that in Eq.(1). The general fitness of the population improves until the GP eventually converges to a model description of the system.The Genetic programming algorithm For this research, a steady-state Genetic-programming a
12、lgorithm was used. At each generation, two parents are selected from the population and the offspring resulting from their crossover operation replace an existing member of the same population. The number of crossover operations is equal to the size of the population i.e. the crossover rate is 100.
13、The crossover algorithm used was a subtree crossover with a limit on the depth of the resulting tree. Genetic programming parameters such as mutation rate and population size varied according to the application. More difficult problems where the expected model structure is complex or where the data
14、are noisy generally require larger population sizes. Mutation rate did not appear to have a significant effect for the systems investigated during this research. Typically, a value of about 2 was chosen. The function library varied according to application rate and what type of nonlinearity might be
15、 expected in the system being identified. A core of linear blocks was always available. It was found that specific nonlinearity such as look-up tables which represented a physical phenomenon would only be selected by the Genetic Programming algorithm if that nonlinearity actually existed in the dyna
16、mic system. This allows the system to be tested for specific nonlinearities.Programming model structure identification Each member of the Genetic Programming population represents a candidate model for the system. It is necessary to evaluate each model and assign to it some fitness value. Each candi
17、date is integrated using a numerical integration routine to produce a time response. This simulation time response is compared with experimental data to give a fitness value for that model. A sum of squared error function (Eq.(1) is used in all the work described in this paper, although many other f
18、itness functions could be used. The simulation routine must be robust. Inevitably, some of the candidate models will be unstable and therefore, the simulation program must protect against overflow error. Also, all system must return a fitness value if the GP algorithm is to work properly even if tho
19、se systems are unstable.Parameter estimation Many of the nodes of the GP trees contain numerical parameters. These could be the coefficients of the transfer functions, a gain value or in the case of a time delay, the delay itself. It is necessary to identify the numerical parameters of each nonlinea
20、r model before evaluating its fitness. The models are randomly generated and can therefore contain linearly dependent parameters and parameters which have no effect on the output. Because of this, gradient based methods cannot be used. Genetic Programming can be used to identify numerical parameters
21、 but it is less efficient than other methods. The approach chosen involves a combination of the Nelder-Simplex and simulated annealing methods. Simulated annealing optimizes by a method which is analogous to the cooling process of a metal. As a metal cools, the atoms organize themselves into an orde
22、red minimum energy structure. The amount of vibration or movement in the atoms is dependent on temperature. As the temperature decreases, the movement, though still random, become smaller in amplitude and as long as the temperature decreases slowly enough, the atoms order themselves slowly enough, t
23、he atoms order themselves into the minimum energy structure. In simulated annealing, the parameters start off at some random value and they are allowed to change their values within the search space by an amount related to a quantity defined as system temperature. If a parameter change improves over
24、all fitness, it is accepted, if it reduces fitness it is accepted with a certain probability. The temperature decreases according to some predetermined cooling schedule and the parameter values should converge to some solution as the temperature drops. Simulated annealing has proved particularly eff
25、ective when combines with other numerical optimization techniques. One such combination is simulated annealing with Nelder-simplex is an (n+1) dimensional shape where n is the number of parameters. This simples explores the search space slowly by changing its shape around the optimum solution .The s
26、imulated annealing adds a random component and the temperature scheduling to the simplex algorithm thus improving the robustness of the method . This has been found to be a robust and reasonably efficient numerical optimization algorithm. The parameter estimation phase can also be used to identify o
27、ther numerical parameters in part of the model where the structure is known but where there are uncertainties about parameter values.Representation of a GP candidate modelNonlinear time domain continuous dynamic models can take a number of different forms. Two common representations involve sets of
28、differential equations or block diagrams. Both these forms of model are well known and relatively easy to simulate .Each has advantages and disadvantages for simulation, visualization and implementation in a Genetic Programming algorithm. Block diagram and equation based representations are consider
29、ed in this paper along with a third hybrid representation incorporating integral and differential operators into an equation based representation.Choice of experimental data setexperimental designThe identification of nonlinear systems presents particular problems regarding experimental design. The
30、system must be excited across the frequency range of interest as with a linear system, but it must also cover the range of any nonlinearities in the system. This could mean ensuring that the input shape is sufficiently varied to excite different modes of the system and that the data covers the opera
31、tional range of the system state space.A large training data set will be required to identify an accurate model. However the simulation time will be proportional to the number of data points, so optimization time must be balanced against quantity of data. A recommendation on how to select efficient
32、step and PRBS signals to cover the entire frequency rage of interest may be found in Godfrey and Ljungs texts.Model validation An important part of any modeling procedure is model validation. The new model structure must be validated with a different data set from that used for the optimization. There are many techniques for validation of nonlinear models, the simplest of which is analogue matching where the time response of the model is compared with available response data from the real system. The model validation results can be used to refine the Genetic Programming algorithm as
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1