We consider the problem of approximating a $p\times m$ rational transfer function $H(s)$ of high degree by another $p\times m$ rational transfer function $\hat{H}(s)$ of much smaller degree. We derive the gradients of the $H_2$-norm of the approximation error and show how stationary points can be described via tangential interpolation. CitationTechnical report UCL-INMA-2007.034, Department of … Read more
Given a separable strongly self-concordant function f:Rn -> R, we show the associated spectral function F(X)= (foL)(X) is also strongly self-concordant function. In addition, there is a universal constant O such that, if f(x) is separable self-concordant barrier then O^2F(X) is a self-concordant barrier. We estimate that for the universal constant we have O
A convex body K has associated with it a unique circumscribed ellipsoid CE(K) with minimum volume, and a unique inscribed ellipsoid IE(K) with maximum volume. We first give a unified, modern exposition of the basic theory of these extremal ellipsoids using the semi-infinite programming approach pioneered by Fritz John in his seminal 1948 paper. We … Read more
Elementary symmetric polynomials can be thought of as derivative polynomials of $E_n(x)=\prod_{i=1,\ldots,n} x_i$. Their associated hyperbolicity cones give a natural sequence of relaxations for $\Re^n_+$. We establish a recursive structure for these cones, namely, that the coordinate projections of these cones are themselves hyperbolicity cones associated with elementary symmetric polynomials. As a consequence of this … Read more
In this paper we present the convergence analysis of the inexact infeasible path-following(IIPF) interior point algorithm. In this algorithm the preconditioned conjugate gradient method is used to solve the reduced KKT system (the augmented system). The augmented system is preconditioned by using a block triangular matrix. The KKT system is solved approximately. Therefore, it becomes … Read more
In this paper we present a cutting plane algorithm for the Set Covering problem. Cutting planes are generated by a “general” (i.e. not based on the “template paradigm”) separation algorithm based on the following idea: i) identify a suitably small subproblem defined by a subset of the constraints of the formulation; ii) run an exact … Read more
This paper describes an active-set algorithm for nonlinear programming that solves a parametric linear programming subproblem at each iteration to generate an estimate of the active set. A step is then computed by solving an equality constrained quadratic program based on this active-set estimate. This approach respresents an extension of the standard sequential linear-quadratic programming … Read more
A regularized Newton-like method for solving nonnegative least-squares problems is proposed and analysed in this paper. A preconditioner for KKT systems arising in the method is introduced and spectral properties of the preconditioned matrix are analysed. A bound on the condition number of the preconditioned matrix is provided. The bound does not depend on the … Read more
The idea of a finite collection of closed sets having “strongly regular intersection” at a given point is crucial in variational analysis. We show that this central theoretical tool also has striking algorithmic consequences. Specifically, we consider the case of two sets, one of which we assume to be suitably “regular” (special cases being convex … Read more
The transportation paradox is related to the classical transportation problem. For certain instances of this problem an increase in the amount of goods to be transported may lead to a decrease in the optimal total transportation cost. Even though the paradox has been known since the early days of linear programming, it has got very … Read more