阶梯状黄土边坡稳定性分析的改进遗传算法及应用

阶梯状黄土边坡稳定性分析的关键是估算其稳定系数的最小值.稳定系数的求解涉及诸多因素且计算过程繁杂,传统优化算法往往不能有效地搜索到其全局最小解.为此,提出一种改进的自适应遗传算法.算法对基因变量空间进行网格状划分,采用迭代选优法建立均匀分布的初始种群,运用优质个体保留遗传策略,并按照特定的准则自适应地调整交叉概率和变异概率,提高算法的全局搜索能力和收敛速度.实例应用表明算法能够快速有效地收敛于土坡稳定系数的全局最小解,且计算结果与实际情况更加吻合.     

Evaluating Performance of Honey Bee Mating Optimization

The honey bee mating optimization (HBMO) algorithm is presented and tested with various test functions, and its performance is compared with the genetic algorithm (GA). It is shown that the HBMO algorithm can overcome the weaknesses of the GA. The HBMO converges faster than the GA. Even when the HMBO starts from a more improper initial condition than the GA, it can reach a better solution in a smaller number of function evaluations. Furthermore, in some cases, the GA was not able to reach the global minimum.     

改进GA机制融合的二进制蜂群算法研究

为了克服人工蜂群算法蜜源更新过程中的随机性并保留蜜源中个体序列合理的组合形式,通过分析基本蜂群算法更新公式的机理,提出一种改进GA(Genetic A1gorithm)机制融合的二进制蜂群算法.算法以二进制编码,首先依概率对任意两蜜源进行"去同存异"操作后随机排列,将排列结果放入到其中某个体中形成新个体.然后依概率进行二进制个体的"翻转"操作,上述两种操作从其本质上相当于GA的类交叉和类变异操作;其次利用GA机制收敛性的证明方式在理论上证明算法是收敛的.最后通过应用不同特性的多维基准函数和算法之间的比较验证改进蜂群算法具有良好的收敛能力和鲁棒性.     

On a Global Optimization Algorithm for Bivariate Smooth Functions

The problem of approximating the global minimum of a function of two variables is considered. A method is proposed rooted in the statistical approach to global optimization. The proposed algorithm partitions the feasible region using a Delaunay triangulation. Only the objective function values are required by the optimization algorithm. The asymptotic convergence rate is analyzed for a class of smooth functions. Numerical examples are provided.     

Order reduction of linear discrete systems using a genetic algorithm

Order reduction of linear discrete systems using classical methods of optimization is well understood and developed by various workers. The present effort is towards development of a method of linear discrete system order reduction using a genetic algorithm (GA) to get rid of usual difficulties of classical methods. The method developed is applied to a variety of systems and the reduced order models are obtained using the method. The results reported are encouraging and more work can be initiated in this area using a GA.     

The Development of Information Guided Evolution Algorithm for Global Optimization

Evolutionary algorithm (EA) has become popular in global optimization with applications widely used in many industrial areas. However, there exists probable premature convergence problem when rugged contour situation is encountered. As to the original genetic algorithm (GA), no matter single population or multi-population cases, the ways to prevent the problem of probable premature convergence are to implement various selection methods, penalty functions and mutation approaches. This work proposes a novel approach to perform very efficient mutation to prevent from premature convergence by introducing the concept of information theory. Information-guided mutation is implemented to several variables, which are selected based on the information entropy derived in this work. The areas of search are also determined on the basis of the information amount obtained from previous searches. Several benchmark problems are solved to show the superiority of this information-guided EA. An industrial scale problem is also presented in this work.     

A combined direction stochastic approximation algorithm

The stochastic approximation problem is to find some root or minimum of a nonlinear function in the presence of noisy measurements. The classical algorithm for stochastic approximation problem is the Robbins-Monro (RM) algorithm, which uses the noisy negative gradient direction as the iterative direction. In order to accelerate the classical RM algorithm, this paper gives a new combined direction stochastic approximation algorithm which employs a weighted combination of the current noisy negative gradient and some former noisy negative gradient as iterative direction. Both the almost sure convergence and the asymptotic rate of convergence of the new algorithm are established. Numerical experiments show that the new algorithm outperforms the classical RM algorithm.     

SUSTAINABLE FOREST MANAGEMENT FOR OPTIMIZING MULTISPECIES WILDLIFE HABITAT: A COASTAL DOUGLAS-FIR EXAMPLE

Wildlife species viability optimization models are developed to convert a given set of initial forest conditions, through a combination of natural growth and management treatments, to a forest system which addresses the joint habitat needs of multispecies populations over time. A linear model of forest cover and wildlife populations is used to form a system of forest management control variables for wildlife habitat modification. The paper examines two objective functions coupled to this system for optimizing sustainable joint species viability. The first maximizes the product of periodic joint viabilities over all time periods, focusing management resources on long-term equilibria, with less emphasis on conversion strategy. The second iteratively maximizes the minimum periodic joint viability over all time periods. This focuses management resources on the most limiting time periods, typically the conversion phase periods. Both objective functions resulted in either point or cyclic equilibria, with cycle lengths equal to minimum forest treatment ages. A third objective, based on maximizing the minimum individual species periodic viability is used to examine single species emphasis. Examples are developed through a case study of 92 vertebrate species found in coastal Douglas-fir stands of northwestern California.     

OMEGA one multi ethnic genetic approach

The genetic algorithm (GA) is a quite efficient paradigm to solve several optimization problems. It is substantially a search technique that uses an ever-changing neighborhood structure related to a population which evolves according to a number of genetic operators. In the GA framework many techniques have been devised to escape from a local optimum when the algorithm fails in locating the global one. To this aim we present a variant of the GA which we call OMEGA (One Multi Ethnic Genetic Approach). The main difference is that, starting from an initial population, \(k\) different sub-populations are produced at each iteration and they independently evolve in \(k\) different environments. The resulting sub–populations are then recombined and the process is iterated. Our basic algorithmic scheme is tested on a recent and well-studied variant of the classic problem of the minimum spanning tree: the Minimum Labeling Spanning Tree problem. We compare our algorithm with several approaches drawn from the literature. The results are encouraging in view of future application of OMEGA to other classes of problems.     

Augmented lagrangian nonlinear programming algorithm that uses SQP and trust region techniques

An augmented Lagrangian nonlinear programming algorithm has been developed. Its goals are to achieve robust global convergence and fast local convergence. Several unique strategies help the algorithm achieve these dual goals. The algorithm consists of three nested loops. The outer loop estimates the Kuhn-Tucker multipliers at a rapid linear rate of convergence. The middle loop minimizes the augmented Lagrangian functions for fixed multipliers. This loop uses the sequential quadratic programming technique with a box trust region stepsize restriction. The inner loop solves a single quadratic program. Slack variables and a constrained form of the fixed-multiplier middleloop problem work together with curved line searches in the inner-loop problem to allow large penalty wieghts for rapid outer-loop convergence. The inner-loop quadratic programs include quadratic onstraint terms, which complicate the inner loop, but speed the middle-loop progress when the constraint curvature is large.The new algorithm compares favorably with a commercial sequential quadratic programming algorithm on five low-order test problems. Its convergence is more robust, and its speed is not much slower.This research was supported in part by the National Aeronautics and Space Administration under Grant No. NAG-1-1009.     

A Finite Algorithm for Global Minimization of Separable Concave Programs

Researchers first examined the problem of separable concave programming more than thirty years ago, making it one of the earliest branches of nonlinear programming to be explored. This paper proposes a new algorithm that finds the exact global minimum of this problem in a finite number of iterations. In addition to proving that our algorithm terminates finitely, the paper extends a guarantee of finiteness to all branch-and-bound algorithms for concave programming that (1) partition exhaustively using rectangular subdivisions and (2) branch on the incumbent solution when possible. The algorithm uses domain reduction techniques to accelerate convergence; it solves problems with as many as 100 nonlinear variables, 400 linear variables and 50 constraints in about five minutes on an IBM RS/6000 Power PC. An industrial application with 152 nonlinear variables, 593 linear variables, and 417 constraints is also solved in about ten minutes.     

A Combined Global & Local Search (CGLS) Approach to Global Optimization

This paper presents a general approach that combines global search strategies with local search and attempts to find a global minimum of a real valued function of n variables. It assumes that derivative information is unreliable; consequently, it deals with derivative free algorithms, but derivative information can be easily incorporated. This paper presents a nonmonotone derivative free algorithm and shows numerically that it may converge to a better minimum starting from a local nonglobal minimum. This property is then incorporated into a random population to globalize the algorithm. Convergence to a zero order stationary point is established for nonsmooth convex functions, and convergence to a first order stationary point is established for strictly differentiable functions. Preliminary numerical results are encouraging. A Java implementation that can be run directly from the Web allows the interested reader to get a better insight of the performance of the algorithm on several standard functions. The general framework proposed here, allows the user to incorporate variants of well known global search strategies. Research done under the cooperation agreement between Universidade de Vigo and Universidad Simón Bolívar.     

Estimating the parameters of a fatigue model using Benders’ decomposition

This paper shows how Benders decomposition can be used for estimating the parameters of a fatigue model. The objective function of such model depends on five parameters of different nature. This makes the parameter estimation problem of the fatigue model suitable for the Benders decomposition, which allows us to use well-behaved and robust parameter estimation methods for the different subproblems. To build the Benders cuts, explicit formulas for the sensitivities (partial derivatives) are obtained. This permits building the classical iterative method, in which upper and lower bounds of the optimal value of the objective function are obtained until convergence. Two alternative objective functions to be optimized are the likelihood and the sum of squares error functions, which relate to the maximum likelihood and the minimum error principles, respectively. The method is illustrated by its application to a real-world problem.     

Global Solution of Optimization Problems with Parameter-Embedded Linear Dynamic Systems

This paper develops a theory for the global solution of nonconvex optimization problems with parameter-embedded linear dynamic systems. A quite general problem formulation is introduced and a solution is shown to exists. A convexity theory for integrals is then developed to construct convex relaxations for utilization in a branch-and-bound framework to calculate a global minimum. Interval analysis is employed to generate bounds on the state variables implied by the bounds on the embedded parameters. These bounds, along with basic integration theory, are used to prove convergence of the branch-and-bound algorithm to the global minimum of the optimization problem. The implementation of the algorithm is then considered and several numerical case studies are examined thoroughly     

一类新的求解P_*(κ)阵线性互补问题的多项式内点算法

利用核函数及其性质,对P_*(k)阵线性互补问题提出了一种新的宽邻域不可行内点算法.对核函数作了一些适当的改进,所以是不同于Peng等人介绍的自正则障碍函数.最后证明了算法具有近似O((1+2k)n3/4log(nμ~0)/ε)多项式复杂性,是优于传统的基于对数障碍函数求解宽邻域内点算法的复杂性.     

Annealing a Genetic Algorithm for Constrained Optimization

In this paper, we adapt a genetic algorithm for constrained optimization problems. We use a dynamic penalty approach along with some form of annealing, thus forcing the search to concentrate on feasible solutions as the algorithm progresses. We suggest two different general-purpose methods for guaranteeing convergence to a globally optimal (feasible) solution, neither of which makes any assumptions on the structure of the optimization problem. The former involves modifying the GA evolution operators to yield a Boltzmann-type distribution on populations. The latter incorporates a dynamic penalty along with a slow annealing of acceptance probabilities. We prove that, with probability one, both of these methods will converge to a globally optimal feasible state.     

Flying-V型仓储布局货位分配优化方法研究

Flying-V是一种典型的非传统布局方式,根据其布局方式的特性,针对仓储货位分配优化问题,以货物出入库效率最高和货物存放的重心最低为优化目标,建立了货位分配多目标优化模型,并采用自适应策略的遗传算法(GA),以及粒子群算法(PSO)进行求解。根据货位分配的优化特点,在GA算法的选择、交叉和变异环节均采用自适应策略, 同时采用惯性权重线性递减的方法设计了PSO算法,有效地解决了两种算法收敛速度慢和易“早熟”的问题,提高了算法的寻优性能。为了更好地表现两种优化求解算法的有效性和优越性,结合具体的货位分配实例利用MATLAB软件编程实现。通过对比分析优化结果表明,PSO算法在收敛速度和优化效果方面相比于自适应GA算法更具有优势,更加合适于解决Flying-V型仓储布局货位分配优化问题。     

Inspired grey wolf optimizer for solving large-scale function optimization problems

Grey wolf optimizer algorithm was recently presented as a new heuristic search algorithm with satisfactory results in real-valued and binary encoded optimization problems that are categorized in swarm intelligence optimization techniques. This algorithm is more effective than some conventional population-based algorithms, such as particle swarm optimization, differential evolution and gravitational search algorithm. Some grey wolf optimizer variants were developed by researchers to improve the performance of the basic grey wolf optimizer algorithm. Inspired by particle swarm optimization algorithm, this study investigates the performance of a new algorithm called Inspired grey wolf optimizer which extends the original grey wolf optimizer by adding two features, namely, a nonlinear adjustment strategy of the control parameter, and a modified position-updating equation based on the personal historical best position and the global best position. Experiments are performed on four classical high-dimensional benchmark functions, four test functions proposed in the IEEE Congress on Evolutionary Computation 2005 special session, three well-known engineering design problems, and one real-world problem. The results show that the proposed algorithm can find more accurate solutions and has higher convergence rate and less number of fitness function evaluations than the other compared techniques.     

A convergence theory for an overlapping Schwarz algorithm using discontinuous iterates

Summary. A new type of overlapping Schwarz methods, using discontinuous iterates, is constructed by modifying the classical overlapping Schwarz algorithm. This new algorithm allows for discontinuous iterates across the artificial interface. For Poissons equation, this algorithm can be considered as an overlapping version of Lions Robin iteration method for which little is known concerning the rate of convergence. Since overlap improves the performance of the classical algorithms considerably, the existence of a uniform convergence factor is the fundamental question for our algorithm. A new theory using Lagrange multipliers is developed and conditions are found for the existence of an almost uniform convergence factor for the dual variables, which implies rapid convergence of the primal variables, in the two overlapping subdomain case. Our result also shows a relation between the boundary conditions of the given problem and the artificial interface condition. Numerical results for the general case with cross points are also presented. They indicate possible extensions of our results to this more general case.Mathematics Subject Classification (2000): 65F10, 65N30, 65N55Acknowledgement I would like to thank my advisor Olof Widlund for suggesting this problem, for many helpful and interesting discussions, and for all his encouragement. I also thank the referees for their helpful corrections and suggestions.     

A memetic algorithm based on multiple learning procedures for global optimal design of composite structures

The mechanical properties of composite structures are sensitive to the choice of the material systems, the distribution of laminates on the structure, and the sizing at the laminate level. This sensitivity depends on many variables. To overcome the difficulties in attaining a global solution for the optimisation problem of composite structures, a new strategy is proposed based on the hybridisation of the memetic algorithm (MA) and the selfish gene (SG) algorithm. Instead of a local search, as performed in MAs, the selfish gene (SG) theory is applied, which follows a different learning scheme in which the conventional population of the individuals is replaced by a virtual population of alleles. The proposed approach, which is called the memetic-based selfish gene algorithm (MA + SG), is a mixed model that applies multiple learning procedures to explore the synergy of different cultural transmission rules in the evolutionary process. The principal aspects of the approach are as follows: co-evolution of multiple populations, species conservation, migration rules, self-adaptive multiple crossovers, local search in hybrid crossover with local genetic improvements, controlled mutation, individual age control, and feature-based allele statistical analysis. To discuss the capabilities of the proposed approach, numerical examples are represented to compare the results of MA + SG with those obtained using genetic algorithms (GA) and MA. The numerical results of the comparison tests showed that the GA and the MA maintained long periods without the evolution of the best-fitted individual/solutions. This behaviour during evolution is associated with the slow maturation of the elite group of populations in the GA and MA approaches. This behaviour is avoided when the MA + SG proposed approach is used with computational cost benefits.