Optimal cycle for a signalized intersection using Global Optimization and Complementarity

A queueing system resulting from a signalised intersection regulated by pre-timed control in an urban traffic network is considered in this paper. Subsequently, we analyse the manner in which Global Optimisation and Complementarity may be used to determine the optimal cycle length and green split allocation for an isolated signalised intersection. The model in question has been formulated as a Mathematical Program with Equilibrium (or Complementarity) Constraints (MPEC). A?sequential complementarity algorithm for computing a global minimum for the MPEC is also subject to analysis in this paper. Furthermore, computational experience is included to demonstrate the efficiency of this method as an effective solution for the problem in question.     

On the numerical solution of the quadratic eigenvalue complementarity problem

The Quadratic Eigenvalue Complementarity Problem (QEiCP) is an extension of the Eigenvalue Complementarity Problem (EiCP) that has been introduced recently. Similar to the EiCP, the QEiCP always has a solution under reasonable hypotheses on the matrices included in its definition. This has been established in a previous paper by reducing a QEiCP of dimension n to a special 2n-order EiCP. In this paper we propose an enumerative algorithm for solving the QEiCP by exploiting this equivalence with an EiCP. The algorithm seeks a global minimum of a special Nonlinear Programming Problem (NLP) with a known global optimal value. The algorithm is shown to perform very well in practice but in some cases terminates with only an approximate optimal solution to NLP. Hence, we propose a hybrid method that combines the enumerative method with a fast and local semi-smooth method to overcome the latter drawback. This algorithm is also shown to be useful for computing a positive eigenvalue for an EiCP under similar assumptions. Computational experience is reported to demonstrate the efficacy and efficiency of the hybrid enumerative method for solving the QEiCP.     

单调线性互补问题的最小原则(II)

在文［１］的基础上，对单调线性互补问题（ＭＬＣＰ）给出了不同于文［１７］的最小原则的另一形式，并提出了一个在有限步内求出单调线性互补问题解集的新算法；给出了单调线性互补问题的三个误差界公式．这些公式推广了文［６］的有关结果，并且较文［８］中的误差界表示形式简洁和易于检验．     

Branch and bound algorithm with applications to robust stability

We discuss a branch and bound algorithm for global optimization of NP-hard problems related to robust stability. This includes computing the distance to instability of a system with uncertain parameters, computing the minimum stability degree of a system over a given set of uncertain parameters, and computing the worst case \(H_\infty \) norm over a given parameter range. The success of our method hinges (1) on the use of an efficient local optimization technique to compute lower bounds fast and reliably, (2) a method with reduced conservatism to compute upper bounds, and (3) the way these elements are favorably combined in the algorithm.     

一类求解MPEC的全局收敛性算法

概述MPEC求解上存在的困难；提出将求解普通约束优化问题的信赖域方法用于求解MPEC；在适当条件下建立算法的全局收敛性定理。     

A Benders approach for computing lower bounds for the mirrored traveling tournament problem

A Benders approach for computing lower bounds for the mirrored Traveling Tournament Problem is proposed. The method obtained improved lower bounds for a number of benchmark instances at the Challenge Traveling Tournament Problems homepage http://mat.gsia.cmu.edu/TOURN/.     

Survey on Vector Complementarity Problems

The paper aims at summarizing the main results on Vector Complementarity Problems (VCP), including the existence of a solution and the relations with Vector Variational Inequalities and Vector Optimization Problems. Particular attention will be given to a VCP with a variable domination structure, where the ordering cone depends on the unknown variable.     

Solution of a General Linear Complementarity Problem Using Smooth Optimization and Its Application to Bilinear Programming and LCP

This paper addresses a General Linear Complementarity Problem (GLCP) that has found applications in global optimization. It is shown that a solution of the GLCP can be computed by finding a stationary point of a differentiable function over a set defined by simple bounds on the variables. The application of this result to the solution of bilinear programs and LCPs is discussed. Some computational evidence of its usefulness is included in the last part of the paper. Accepted 28 June 1999. Online publication 4 December 2000.     

Multiplier convergence in trust-region methods with application to convergence of decomposition methods for MPECs

We study piecewise decomposition methods for mathematical programs with equilibrium constraints (MPECs) for which all constraint functions are linear. At each iteration of a decomposition method, one step of a nonlinear programming scheme is applied to one piece of the MPEC to obtain the next iterate. Our goal is to understand global convergence to B-stationary points of these methods when the embedded nonlinear programming solver is a trust-region scheme, and the selection of pieces is determined using multipliers generated by solving the trust-region subproblem. To this end we study global convergence of a linear trust-region scheme for linearly-constrained NLPs that we call a trust-search method. The trust-search has two features that are critical to global convergence of decomposition methods for MPECs: a robustness property with respect to switching pieces, and a multiplier convergence result that appears to be quite new for trust-region methods. These combine to clarify and strengthen global convergence of decomposition methods without resorting either to additional conditions such as eventual inactivity of the trust-region constraint, or more complex methods that require a separate subproblem for multiplier estimation.      

Computing closely matching upper and lower bounds on textile nesting problems

We consider an industrial cutting problem in textile manufacturing and report on algorithms for computing cutting images and lower bounds on waste for this problem. For the upper bounds we use greedy strategies based on hodographs and global optimization based on simulated annealing. For the lower bounds we use branch-and-bound methods for computing optimal solutions of placement subproblems that determine the performance of the overall subproblem. The upper bounds are computed in less than an hour on a common-day workstation and are competitive in quality with results obtained by human nesters. The lower bounds take a few hours to compute and are within 0.4% of the upper bound for certain types of clothing (e.g., for pants).     

Equivalent Lipschitz surrogates for zero-norm and rank optimization problems

This paper proposes a mechanism to produce equivalent Lipschitz surrogates for zero-norm and rank optimization problems by means of the global exact penalty for their equivalent mathematical programs with an equilibrium constraint (MPECs). Specifically, we reformulate these combinatorial problems as equivalent MPECs by the variational characterization of the zero-norm and rank function, show that their penalized problems, yielded by moving the equilibrium constraint into the objective, are the global exact penalization, and obtain the equivalent Lipschitz surrogates by eliminating the dual variable in the global exact penalty. These surrogates, including the popular SCAD function in statistics, are also difference of two convex functions (D.C.) if the function and constraint set involved in zero-norm and rank optimization problems are convex. We illustrate an application by designing a multi-stage convex relaxation approach to the rank plus zero-norm regularized minimization problem.     

Parallel Implementation of Successive Convex Relaxation Methods for Quadratic Optimization Problems

As computing resources continue to improve, global solutions for larger size quadrically constrained optimization problems become more achievable. In this paper, we focus on larger size problems and get accurate bounds for optimal values of such problems with the successive use of SDP relaxations on a parallel computing system called Ninf (Network-based Information Library for high performance computing).     

隐互补约束优化问题的磨光SQP算法

本文提出了一类隐互补约束优化问题的磨光SQP算法.首先,我们给出了这类优化问题的最优性和约束规范性条件.然后,在适当假设条件下,我们证明了算法具有全局收敛性.     

A DC programming approach for solving the symmetric Eigenvalue Complementarity Problem

In this paper, we investigate a DC (Difference of Convex functions) programming technique for solving large scale Eigenvalue Complementarity Problems (EiCP) with real symmetric matrices. Three equivalent formulations of EiCP are considered. We first reformulate them as DC programs and then use DCA (DC Algorithm) for their solution. Computational results show the robustness, efficiency, and high speed of the proposed algorithms.     

On the solution of NP-hard linear complementarity problems

In this paper two enumerative algorithms for the Linear Complementarity Problems (LCP) are discussed. These procedures exploit the equivalence of theLCP into a nonconvex quadratic and a bilinear programs. It is shown that these algorithms are efficient for processing NP-hardLCPs associated with reformulations of the Knapsack problem and should be recommended to solve difficultLCPs.     

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     

On the computation of all eigenvalues for the eigenvalue complementarity problem

In this paper, a parametric algorithm is introduced for computing all eigenvalues for two Eigenvalue Complementarity Problems discussed in the literature. The algorithm searches a finite number of nested intervals \([\bar{l}, \bar{u}]\) in such a way that, in each iteration, either an eigenvalue is computed in \([\bar{l}, \bar{u}]\) or a certificate of nonexistence of an eigenvalue in \([\bar{l}, \bar{u}]\) is provided. A hybrid method that combines an enumerative method [1] and a semi-smooth algorithm [2] is discussed for dealing with the Eigenvalue Complementarity Problem over an interval \([\bar{l}, \bar{u}]\) . Computational experience is presented to illustrate the efficacy and efficiency of the proposed techniques.     

An algorithm for global optimization of Lipschitz continuous functions

An algorithm is presented which locates the global minimum or maximum of a function satisfying a Lipschitz condition. The algorithm uses lower bound functions defined on a partitioned domain to generate a sequence of lower bounds for the global minimum. Convergence is proved, and some numerical results are presented.     

Outer approximation algorithms for separable nonconvex mixed-integer nonlinear programs

A rigorous decomposition approach to solve separable mixed-integer nonlinear programs where the participating functions are nonconvex is presented. The proposed algorithms consist of solving an alternating sequence of Relaxed Master Problems (mixed-integer linear program) and two nonlinear programming problems (NLPs). A sequence of valid nondecreasing lower bounds and upper bounds is generated by the algorithms which converge in a finite number of iterations. A Primal Bounding Problem is introduced, which is a convex NLP solved at each iteration to derive valid outer approximations of the nonconvex functions in the continuous space. Two decomposition algorithms are presented in this work. On finite termination, the first yields the global solution to the original nonconvex MINLP and the second finds a rigorous bound to the global solution. Convergence and optimality properties, and refinement of the algorithms for efficient implementation are presented. Finally, numerical results are compared with currently available algorithms for example problems, illuminating the potential benefits of the proposed algorithm.     

On the Guignard constraint qualification for mathematical programs with equilibrium constraints

We recapitulate the well-known fact that most of the standard constraint qualifications are violated for mathematical programs with equilibrium constraints (MPECs). We go on to show that the Abadie constraint qualification is only satisfied in fairly restrictive circumstances. In order to avoid this problem, we fall back on the Guignard constraint qualification (GCQ). We examine its general properties and clarify the position it occupies in the context of MPECs. We show that strong stationarity is a necessary optimality condition under GCQ. Also, we present several sufficient conditions for GCQ, showing that it is usually satisfied for MPECs.