831 research outputs found
Bilinearity rank of the cone of positive polynomials and related cones
For a proper cone K ⊂ Rn and its dual cone K the complementary slackness condition xT s = 0 defines an n-dimensional manifold C(K) in the space { (x, s) | x ∈ K, s ∈ K^* }. When K is a symmetric cone, this manifold can be described by a set of n bilinear equalities. When K is a symmetric cone, this fact translates to a set of n linearly independent bilinear identities (optimality conditions) satisfied by every (x, s) ∈ C(K). This proves to be very useful when
optimizing over such cones, therefore it is natural to look for similar optimality conditions for non-symmetric cones. In this paper we define the bilinearity rank of a cone, which is the number of linearly independent bilinear identities valid for the cone, and describe a linear algebraic technique to bound this quantity. We examine several well-known cones, in particular
the cone of positive polynomials P2n+1 and its dual, the closure of the moment cone M2n+1, and compute their bilinearity ranks. We show that there are exactly four linearly independent bilinear identities which hold for all (x,s) ∈ C(P2n+1), regardless of the dimension of the cones. For nonnegative polynomials over an interval or half-line there are only two linearly independent bilinear identities. These results are extended to trigonometric and exponential
polynomials
Bilinearity rank of the cone of positive polynomials and related cones
For a proper cone K ⊂ Rn and its dual cone K the complementary slackness condition xT s = 0 defines an n-dimensional manifold C(K) in the space { (x, s) | x ∈ K, s ∈ K^* }. When K is a symmetric cone, this manifold can be described by a set of n bilinear equalities. When K is a symmetric cone, this fact translates to a set of n linearly independent bilinear identities (optimality conditions) satisfied by every (x, s) ∈ C(K). This proves to be very useful when
optimizing over such cones, therefore it is natural to look for similar optimality conditions for non-symmetric cones. In this paper we define the bilinearity rank of a cone, which is the number of linearly independent bilinear identities valid for the cone, and describe a linear algebraic technique to bound this quantity. We examine several well-known cones, in particular
the cone of positive polynomials P2n+1 and its dual, the closure of the moment cone M2n+1, and compute their bilinearity ranks. We show that there are exactly four linearly independent bilinear identities which hold for all (x,s) ∈ C(P2n+1), regardless of the dimension of the cones. For nonnegative polynomials over an interval or half-line there are only two linearly independent bilinear identities. These results are extended to trigonometric and exponential
polynomials
New Formulation and Strong MISOCP Relaxations for AC Optimal Transmission Switching Problem
As the modern transmission control and relay technologies evolve,
transmission line switching has become an important option in power system
operators' toolkits to reduce operational cost and improve system reliability.
Most recent research has relied on the DC approximation of the power flow model
in the optimal transmission switching problem. However, it is known that DC
approximation may lead to inaccurate flow solutions and also overlook stability
issues. In this paper, we focus on the optimal transmission switching problem
with the full AC power flow model, abbreviated as AC OTS. We propose a new
exact formulation for AC OTS and its mixed-integer second-order conic
programming (MISOCP) relaxation. We improve this relaxation via several types
of strong valid inequalities inspired by the recent development for the closely
related AC Optimal Power Flow (AC OPF) problem. We also propose a practical
algorithm to obtain high quality feasible solutions for the AC OTS problem.
Extensive computational experiments show that the proposed formulation and
algorithms efficiently solve IEEE standard and congested instances and lead to
significant cost benefits with provably tight bounds
A Parametric Multi-Convex Splitting Technique with Application to Real-Time NMPC
A novel splitting scheme to solve parametric multiconvex programs is
presented. It consists of a fixed number of proximal alternating minimisations
and a dual update per time step, which makes it attractive in a real-time
Nonlinear Model Predictive Control (NMPC) framework and for distributed
computing environments. Assuming that the parametric program is semi-algebraic
and that its KKT points are strongly regular, a contraction estimate is derived
and it is proven that the sub-optimality error remains stable if two key
parameters are tuned properly. Efficacy of the method is demonstrated by
solving a bilinear NMPC problem to control a DC motor.Comment: To appear in Proceedings of the 53rd IEEE Conference on Decision and
Control 201
Partitioning Procedure for Polynomial Optimization: Application to Portfolio Decisions with Higher Order Moments
We consider the problem of finding the minimum of a real-valued multivariate polynomial function constrained in a compact set defined by polynomial inequalities and equalities. This problem, called polynomial optimization problem (POP), is generally nonconvex and has been of growing interest to many researchers in recent years. Our goal is to tackle POPs using decomposition. Towards this goal we introduce a partitioning procedure. The problem manipulations are in line with the pattern used in the Benders decomposition [1], namely relaxation preceded by projection. Stengle’s and Putinar’s Positivstellensatz are employed to derive the so-called feasibility and optimality constraints, respectively. We test the performance of the proposed method on a collection of benchmark problems and we present the numerical results. As an application, we consider the problem of selecting an investment portfolio optimizing the mean, variance, skewness and kurtosis of the portfolio.Polynomial optimization, Semidefinite relaxations, Positivstellensatz, Sum of squares, Benders decomposition, Portfolio optimization
A Parametric Non-Convex Decomposition Algorithm for Real-Time and Distributed NMPC
A novel decomposition scheme to solve parametric non-convex programs as they
arise in Nonlinear Model Predictive Control (NMPC) is presented. It consists of
a fixed number of alternating proximal gradient steps and a dual update per
time step. Hence, the proposed approach is attractive in a real-time
distributed context. Assuming that the Nonlinear Program (NLP) is
semi-algebraic and that its critical points are strongly regular, contraction
of the sequence of primal-dual iterates is proven, implying stability of the
sub-optimality error, under some mild assumptions. Moreover, it is shown that
the performance of the optimality-tracking scheme can be enhanced via a
continuation technique. The efficacy of the proposed decomposition method is
demonstrated by solving a centralised NMPC problem to control a DC motor and a
distributed NMPC program for collaborative tracking of unicycles, both within a
real-time framework. Furthermore, an analysis of the sub-optimality error as a
function of the sampling period is proposed given a fixed computational power.Comment: 16 pages, 9 figure
Matrix Minor Reformulation and SOCP-based Spatial Branch-and-Cut Method for the AC Optimal Power Flow Problem
Alternating current optimal power flow (AC OPF) is one of the most
fundamental optimization problems in electrical power systems. It can be
formulated as a semidefinite program (SDP) with rank constraints. Solving AC
OPF, that is, obtaining near optimal primal solutions as well as high quality
dual bounds for this non-convex program, presents a major computational
challenge to today's power industry for the real-time operation of large-scale
power grids. In this paper, we propose a new technique for reformulation of the
rank constraints using both principal and non-principal 2-by-2 minors of the
involved Hermitian matrix variable and characterize all such minors into three
types. We show the equivalence of these minor constraints to the physical
constraints of voltage angle differences summing to zero over three- and
four-cycles in the power network. We study second-order conic programming
(SOCP) relaxations of this minor reformulation and propose strong cutting
planes, convex envelopes, and bound tightening techniques to strengthen the
resulting SOCP relaxations. We then propose an SOCP-based spatial
branch-and-cut method to obtain the global optimum of AC OPF. Extensive
computational experiments show that the proposed algorithm significantly
outperforms the state-of-the-art SDP-based OPF solver and on a simple personal
computer is able to obtain on average a 0.71% optimality gap in no more than
720 seconds for the most challenging power system instances in the literature
Nonlinear Integer Programming
Research efforts of the past fifty years have led to a development of linear
integer programming as a mature discipline of mathematical optimization. Such a
level of maturity has not been reached when one considers nonlinear systems
subject to integrality requirements for the variables. This chapter is
dedicated to this topic.
The primary goal is a study of a simple version of general nonlinear integer
problems, where all constraints are still linear. Our focus is on the
computational complexity of the problem, which varies significantly with the
type of nonlinear objective function in combination with the underlying
combinatorial structure. Numerous boundary cases of complexity emerge, which
sometimes surprisingly lead even to polynomial time algorithms.
We also cover recent successful approaches for more general classes of
problems. Though no positive theoretical efficiency results are available, nor
are they likely to ever be available, these seem to be the currently most
successful and interesting approaches for solving practical problems.
It is our belief that the study of algorithms motivated by theoretical
considerations and those motivated by our desire to solve practical instances
should and do inform one another. So it is with this viewpoint that we present
the subject, and it is in this direction that we hope to spark further
research.Comment: 57 pages. To appear in: M. J\"unger, T. Liebling, D. Naddef, G.
Nemhauser, W. Pulleyblank, G. Reinelt, G. Rinaldi, and L. Wolsey (eds.), 50
Years of Integer Programming 1958--2008: The Early Years and State-of-the-Art
Surveys, Springer-Verlag, 2009, ISBN 354068274
- …