Real solving polynomial equations with semidefinite programming
|
|
- Marsha Cooper
- 5 years ago
- Views:
Transcription
1 .5cm. Real solving polynomial equations with semidefinite programming Jean Bernard Lasserre - Monique Laurent - Philipp Rostalski LAAS, Toulouse - CWI, Amsterdam - ETH, Zürich LAW 2008 Real solving polynomial equations with semidefinite programming p.1
2 The problem Given polynomials h 1,...,h m R[x] = R[x 1,...,x n ] Compute all common real roots (assuming finitely many), i.e. compute the real variety V R (I) of the ideal I := (h 1,...,h m ) Find a basis of the real radical ideal I(V R (I)) V R (I) := {v R n f(v) = 0 f I} I(V R (I)) := {f R[x] f(v) = 0 v V R (I)} }{{} = {f R[x] m N s i R[x] f 2m + i s2 i I} Real Nullstellensatz Real solving polynomial equations with semidefinite programming p.2
3 Our contribution 1. A semidefinite characterization of I(V R (I)) [as the kernel of some positive semidefinite moment matrix] 2. Assuming V R (I) <, an algorithm for finding: a generating set (border or Gröbner basis) of I(V R (I)) the real variety V R (I) Remarks about the method: real algebraic in nature: no complex roots computed works if V R (I) is finite (even if V C (I) is not) no preliminary Gröbner basis of I is needed numerical, based on semidefinite programming (SDP) Real solving polynomial equations with semidefinite programming p.3
4 Plan of the talk 1. The moment-matrix method for V R (I) 2. Adapt the moment-matrix method for V C (I) [drop PSD] 3. Relate to the prolongation-projection algorithm of Zhi and Reid for V C (I) 4. Adapt the prolongation-projection algorithm for V R (I) [add PSD] Real solving polynomial equations with semidefinite programming p.4
5 The complex case is well understood Given an ideal I R[x] with V C (I) <, find the (complex) variety V C (I) and the radical ideal I(V C (I)). Linear algebra in the finite dimensional space R[x]/I Need a linear basis of R[x]/I and a normal form algorithm V C (I) can be computed e.g. with: Linear algebra methods: Eigenvalue method [Stetter-Möller, Stickelberger, Rouillier] Homotopy methods [Verschelde]... Seidenberg [1974]: I(V C (I)) = (I {q 1,...,q n }), where q i is the square-free part of p i, the monic generator of I R[x i ]. Real solving polynomial equations with semidefinite programming p.5
6 The eigenvalue method for V C (I) <, i.e. dim R[x]/I < Stickelberger theorem: Let m f be the multiplication by f linear operator in R[x]/I. 1. The eigenvalues of m f are {f(v) v V C (I)}. 2. The eigenvectors of m T f give the points v V C(I). M T f ζ B,v = f(v)ζ B,v v V C (I) where M f is the matrix of m f in a base B of R[x]/I and ζ B,v := (b(v)) b B Moreover, when B is a set of monomials and 1 B, a border basis of I can be read directly from the multiplication matrices M x1,...,m xn. Real solving polynomial equations with semidefinite programming p.6
7 Finding a linear basis B of R[x]/I and a basis G of the ideal I Typically: G is a Gröbner basis and B is the set of standard monomials for a given monomial ordering (e.g. via Buchberger s algorithm) More generally: Assume B = {b 1 = 1,b 2,...,b N } is a set of monomials with border B := (x 1 B... x n B) \ B. Write any border monomial x i b j = N k=1 a (ij) k b k } {{ } Span(B) + g (ij) }{{} I Then: G := {g (ij) x i b j B} is a (border) basis of I and carries the same information as the multiplication matrices M x1,...,m xn Real solving polynomial equations with semidefinite programming p.7
8 Counting real roots with the Hermite quadratic form For f R[x] Hermite bilinear form: H f : R[x]/I R[x]/I R (g,h) Tr(M fgh ) Theorem: For f = 1 rank(h 1 ) = V C (I), Sign(H 1 ) = V R (I), Rad (H 1 ) = I(V C (I)) rank(h f ) = {v V C (I) f(v) 0} Sign(H f ) = {v V R (I) f(v) > 0} {v V R (I) f(v) < 0} Real solving polynomial equations with semidefinite programming p.8
9 To find V R (I) and a basis of the real radical ideal I(V R (I)) it suffices to have a linear basis B of R[x]/I(V R (I)) and the multiplication matrices in R[x]/I(V R (I))! New tool: Moment matrices y R Nn 2s Ms (y) := (y α+β ) α,β N n s N n s := {α N n α = i α i s} monomials x α of degree s Motivation: For y = (v α ) α N n 2s =: ζ 2s,v where v R n M s (y) = ζ s,v ζ T s,v 0 and KerM s (y) I(v) Real solving polynomial equations with semidefinite programming p.9
10 Real roots of I = (h 1,...,h m ) and PSD moment matrices Lemma: For v V R (I) and t D := max j deg(h j ) the vector y = ζ t,v = (v α ) α t satisfies: the linear constraints (LC): y T [v V C (I)] (h j x α ) = 0 j = 1...m α s.t. α + deg(h j ) t the PSD constraint: M t/2 (y) 0 [v R n ] Set: K t := {y R Nn t (LC), M t/2 (y) 0} Obviously: K t cone(ζ t,v v V R (I)} Theorem: t s D π s (K t )=cone(ζ s,v v V R (I)} Real solving polynomial equations with semidefinite programming p.10
11 Semidefinite characterization of I(V R (I)) Theorem 1: Let y be a generic element of K t, i.e. y lies in the relative interior of the cone K t. Then (KerM t/2 (y)) I(V R (I)) with equality for t large enough. Geometric property of SDP: y is generic rankm t/2 (y) is maximum KerM t/2 (y) KerM t/2 (z) z K t Thus: for v V R (I), KerM t/2 (y) KerM t/2 (ζ t,v ) I(v). Let {g 1,...,g L } be a basis of I(V R (I)). Real Nullstellensatz: gl 2m + i s2 i = m j=1 u jh j. This implies: g l KerM t/2 (y) for t large enough. Real solving polynomial equations with semidefinite programming p.11
12 Stopping criterion when V R (I) < Theorem 2: Let y be a generic element of K t. Assume one of the following two flatness conditions holds: (F1) rankm s (y) = rankm s 1 (y) for some D s t/2 (Fd) rankm s (y) = rankm s d (y) for some d = D/2 s t/2. Then: I(V R (I)) = (KerM s (y)) Any base B of the column space of M s 1 (y) is a base of R[x]/I(V R (I)) The multiplication matrices can be constructed from M s (y). Real solving polynomial equations with semidefinite programming p.12
13 Sketch of proof: Assume rankm s (y) = rankm s 1 (y) Thm [Curto-Fialkow 1996] π 2s (y) has a flat extension ỹ R Nn, i.e. such that rankm(ỹ) = rankm s (y). Thm [La 2005] As M(ỹ) 0, (KerM s (y))=kerm(ỹ) is a real radical 0-dimensional ideal. I }{{} (LC) (KerM s (y)) I(V }{{} R (I)) y generic Thus: (KerM s (y))=i(v R (I)) B indexes a base of M s 1 (y) = B indexes a base of M(ỹ) = B is a base of R[x]/KerM(ỹ) = R[x]/I(V R (I)) Use linear dependencies in M s (y) to construct the multiplication matrices. Real solving polynomial equations with semidefinite programming p.13
14 The moment-matrix algorithm for V R (I) Input: h 1,...,h m R[x] Output: B base of R[x]/I(V R (I)) The multiplication matrices M xi in R[x]/I(V R (I)) Algorithm: For t D Step 1: Compute a generic element y K t. Step 2: Check if (F1) or (Fd) holds. If yes, return a column basis B of M s 1 (y) and M xi = M 1 B P i, M B := principal submatrix of M s 1 (y) indexed by B P i := submatrix of M s (y) with rows in B and columns in x i B. If no, go to Step 1 with t t + 1. Theorem: The algorithm terminates. Real solving polynomial equations with semidefinite programming p.14
15 The algorithm terminates: (F1) holds for t large enough. For t t 0, KerM t/2 (y) contains a Gröbner base {g 1,...,g L } of I(V R (I)) for a total degree ordering. B := {b 1,...,b N }: set of standard monomials base of R[x]/I(V R (I)). Set: s := 1 + max b B deg(b) and assume t t 0, t/2 > s. N L For α s, write x α = λ i b i + u l g l i=1 }{{} deg s 1 Thus: x α N i=1 λ ib i KerM t/2 (y). That is: rankm s (y) = rankm s 1 (y). l=1 }{{} deg α s< t/2 Real solving polynomial equations with semidefinite programming p.15
16 Two small examples Ex. 1: I = (h := x x2 2 ) V R (I) = {0}, V C (I) =. M 1 (y) 0, 0 = y T h = y20 + y 02 = y α = 0 α 0. Any generic y K 2 is y = (y 0, 0,..., 0) with y 0 > 0. Thus: (KerM 1 (y)) = (x 1,x 2 ) = I(V R (I)). Ex. 2: I = (h i := x i (x 2 i + 1) i = 1,...,n) V R (I) = {0}, V C (I) = 3 n. M 2 (y) 0, 0 = y T (x i h i ) = y 4ei +y 2ei i = y α = 0 α 0. Any generic y K 4 is y = (y 0, 0,..., 0) with y 0 > 0. Thus: (KerM 1 (y)) = (x 1,...,x n ) = I(V R (I)). Real solving polynomial equations with semidefinite programming p.16
17 Some algorithmic issues How to find a generic y K t, i.e. with rankm t (y) max.? Solve the SDP program: min y Kt 1 with a SDP solver using the extended self-dual embedding property. Then the central path converges to a solution in the relative interior of the optimum face, i.e., to a generic point y K t. How to compute ranks of matrices? We use SVD decomposition, but this is a sensitive numerical issue... The method may work without (F1) or (Fd): If rankm B (y) = rankm B B (y) and the formal multiplication matrices commute. Real solving polynomial equations with semidefinite programming p.17
18 Extension of the moment-matrix algorithm to V C (I) Omit the PSD condition and work with the linear space: K t = {y R Nn t y T (h j x α ) = 0 j,α with α + deg(h j ) t} The same algorithm applies: For t D Pick generic y K t, i.e. rankm s (y) maximum s t/2 [choose y K t randomly] Check if the flatness condition (F1) or (Fd) holds. If yes, find a basis of R[x]/J where J := (KerM s (y)) satisfies I J I(V C (I)) and thus V C (J) = V C (I). If not, iterate with t + 1. Real solving polynomial equations with semidefinite programming p.18
19 Find the ideal (KerM s (y)) = I in the Gorenstein case The inclusion I (KerM s (y)) I(V C (I)) may be strict for any generic y. Example: For I = (x 2 1,x2 2,x 1x 2 ), V C (I) = {0}, I(V C (I)) = (x 1,x 2 ), dim R[x]/I = 3, dim R[x]/I(V C (I)) = 1, while dim R[x]/(KerM s (y)) = 2 for any generic y! Recall: The algebra A := R[x]/I is Gorenstein if there exists a non-degenerate bilinear form on A satisfying (f,gh) = (fg,h) f,g,h A, i.e. if there exists y K with I = KerM(y) Hence: y K t s.t. rankm s (y) = rankm s 1 (y) and I = (KerM s (y)) IFF A is Gorenstein. Real solving polynomial equations with semidefinite programming p.19
20 Example: the moment-matrix algorithm for real/complex roots I = (x 2 1 2x 1x 3 + 5, x 1 x x 2x 3 + 1, 3x 2 2 8x 1x 3 ), D = 3, d = 2 Ranks of M s (y) for generic y K t, K t : t = s = s = s = s = s = no PSD 8 complex roots t = s = s = s = s = 3 10 with PSD extract 2 real roots Real solving polynomial equations with semidefinite programming p.20
21 8 complex / 2 real roots: v 1 = v 2 = v 3 = v 4 = v 5 = v 6 = v 7 = v 8 = h h h h h h h h i 1.101, 2.878, i i, i, i i i, i, i i, i, i i, i, i i i, i, i i i, i, i i 0.966, 2.813, i i Real solving polynomial equations with semidefinite programming p.21
22 Extracting real roots without (F1) or (Fd) I = (5x 9 1 6x5 1 x 2 + x 1 x x 1x 3, 2x 6 1 x 2 + 2x 2 1 x x 2x 3,x x ) D = 9, d = 5, V R (I) = 8, V C (I) = 20 order rank sequence of extract. order s accuracy comm. error t M s (y) (1 s t/2 ) MON/SVD MON/SVD MON/SVD (3)/ ( ) / / (3)/3(3) e-5/ e-5/ Quotient basis: B = {1,x 1,x 2,x 3,x 2 1,x 1x 2,x 1 x 3,x 2 x 3 } border basis G of size 10 8 x 1 = ( 0.515, , ) x 2 = ( 0.502,0.119, ) >< x 3 = (0.502, 0.119,0.0124) x 4 = (0.515, , ) Real solutions: x 5 = (0.262, 0.444, ) x 6 = ( 2.07e-5,0.515, 1.27e-6) >: x 7 = ( 0.262, 0.444, ) x 8 = ( 1.05e-5, 0.515, 7.56e-7) Real solving polynomial equations with semidefinite programming p.22
23 Link with the elimination method of Zhi and Reid Theorem: If (F1) holds, i.e. for some D s t/2 rankm s (y) = rankm s 1 (y) for generic y K t, then dimπ 2s (K t ) = dim π 2s 1 (K t ) = dimπ 2s (K t+1 ) Theorem (based on [Zhi-Reid 2004]): If for some D s t (ZR) dim π s (K t ) = dim π s 1 (K t ) = dimπ s (K t+1 ) then one can construct a base of R[x]/I and the multiplication matrices in R[x]/I [and thus extract V C (I)]. Hence: The Zhi-Reid criterion (ZR) may be satisfied earlier than the flatness criterion (F1). Real solving polynomial equations with semidefinite programming p.23
24 Example: I = (x 2 1 2x 1 x 3 + 5,x 1 x x 2 x 3 + 1, 3x 2 2 8x 1 x 3 ) t = s = s = s = s = s = rankm 3 (y)=rankm 2 (y) for y K 9 t = s = s = s = s = s = s = s = s = dim π 3 (K 6 ) = dim π 2 (K 6 ) = dim π 3 (K 7 ) Real solving polynomial equations with semidefinite programming p.24
25 Extending the Zhi-Reid criterion to the real case In the complex case, K t = H t where H t := {h j x α j,α with deg(h j x α ) t}. In the real case, K t is a cone, contained in the linear space P t, with the same dimensions: dim K t = dim P t, where P t := H t {fx α f KerM t/2 (y), deg(x α ) t/2 } Theorem: If for some D s t (ZR+) dim π s (P t ) = dim π s 1(P t ) = dimπ s((p t P t ) ) then one can construct a base of J with I J I(V R (I)) and thus extract V R (I) = V C (J) R n. Real solving polynomial equations with semidefinite programming p.25
26 Link with the flatness criterion Theorem: In the PSD case, the flatness criterion (F1): rankm s (y) = rankm s 1 (y) for generic y K t is equivalent to the stronger version of the (ZR) criterion: (ZR++) dimπ s 1 (P t ) = dim π 2s(P t ) = dimπ 2s((P t P t ) ) in which case we find the real radical ideal J = I(V R (I)). Hence: the algorithm based on (ZR) may stop earlier than the moment-matrix algorithm, based on (F1). Future work: Adapt other known efficient algorithms for complex roots to real roots by incorporating SDP conditions. Real solving polynomial equations with semidefinite programming p.26
Chapter 2 The Approach of Moments for Polynomial Equations
Chapter 2 The Approach of Moments for Polynomial Equations Monique Laurent and Philipp Rostalski 2.1 Introduction Computing all points x K n (K = R or C) at which a given system of polynomials in n variables
More informationSUMS OF SQUARES, MOMENT MATRICES AND OPTIMIZATION OVER POLYNOMIALS. February 5, 2008
SUMS OF SQUARES, MOMENT MATRICES AND OPTIMIZATION OVER POLYNOMIALS MONIQUE LAURENT February 5, 2008 Abstract. We consider the problem of minimizing a polynomial over a semialgebraic set defined by polynomial
More informationOptimization over Polynomials with Sums of Squares and Moment Matrices
Optimization over Polynomials with Sums of Squares and Moment Matrices Monique Laurent Centrum Wiskunde & Informatica (CWI), Amsterdam and University of Tilburg Positivity, Valuations and Quadratic Forms
More informationA Sparse Flat Extension Theorem for Moment Matrices
A Sparse Flat Extension Theorem for Moment Matrices Monique Laurent and Bernard Mourrain Abstract. In this note we prove a generalization of the flat extension theorem of Curto and Fialkow [4] for truncated
More informationLMI MODELLING 4. CONVEX LMI MODELLING. Didier HENRION. LAAS-CNRS Toulouse, FR Czech Tech Univ Prague, CZ. Universidad de Valladolid, SP March 2009
LMI MODELLING 4. CONVEX LMI MODELLING Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ Prague, CZ Universidad de Valladolid, SP March 2009 Minors A minor of a matrix F is the determinant of a submatrix
More informationADVANCED TOPICS IN ALGEBRAIC GEOMETRY
ADVANCED TOPICS IN ALGEBRAIC GEOMETRY DAVID WHITE Outline of talk: My goal is to introduce a few more advanced topics in algebraic geometry but not to go into too much detail. This will be a survey of
More informationTensor decomposition and moment matrices
Tensor decomposition and moment matrices J. Brachat (Ph.D.) & B. Mourrain GALAAD, INRIA Méditerranée, Sophia Antipolis Bernard.Mourrain@inria.fr SIAM AAG, Raleigh, NC, Oct. -9, 211 Collaboration with A.
More informationDetecting global optimality and extracting solutions in GloptiPoly
Detecting global optimality and extracting solutions in GloptiPoly Didier Henrion 1, Jean-Bernard Lasserre 1 September 5, 5 Abstract GloptiPoly is a Matlab/SeDuMi add-on to build and solve convex linear
More informationGlobal polynomial optimization with Moments Matrices and Border Basis
Global polynomial optimization with Moments Matrices and Border Basis Marta Abril Bucero, Bernard Mourrain, Philippe Trebuchet Inria Méditerranée, Galaad team Juillet 15-19, 2013 M. Abril-Bucero, B. Mourrain,
More informationGRÖBNER BASES AND POLYNOMIAL EQUATIONS. 1. Introduction and preliminaries on Gróbner bases
GRÖBNER BASES AND POLYNOMIAL EQUATIONS J. K. VERMA 1. Introduction and preliminaries on Gróbner bases Let S = k[x 1, x 2,..., x n ] denote a polynomial ring over a field k where x 1, x 2,..., x n are indeterminates.
More informationSparse representations from moments
Sparse representations from moments Bernard Mourrain Inria Méditerranée, Sophia Antipolis BernardMourrain@inriafr Sparse representation of signals Given a function or signal f (t): decompose it as f (t)
More informationLecture 7: Positive Semidefinite Matrices
Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.
More informationHilbert s 17th Problem to Semidefinite Programming & Convex Algebraic Geometry
Hilbert s 17th Problem to Semidefinite Programming & Convex Algebraic Geometry Rekha R. Thomas University of Washington, Seattle References Monique Laurent, Sums of squares, moment matrices and optimization
More informationThe Graph Realization Problem
The Graph Realization Problem via Semi-Definite Programming A. Y. Alfakih alfakih@uwindsor.ca Mathematics and Statistics University of Windsor The Graph Realization Problem p.1/21 The Graph Realization
More informationNotes on the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre by M. Laurent, December 13, 2012
Notes on the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre by M. Laurent, December 13, 2012 We present the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre
More informationLecture 15: Algebraic Geometry II
6.859/15.083 Integer Programming and Combinatorial Optimization Fall 009 Today... Ideals in k[x] Properties of Gröbner bases Buchberger s algorithm Elimination theory The Weak Nullstellensatz 0/1-Integer
More informationConvex Optimization & Parsimony of L p-balls representation
Convex Optimization & Parsimony of L p -balls representation LAAS-CNRS and Institute of Mathematics, Toulouse, France IMA, January 2016 Motivation Unit balls associated with nonnegative homogeneous polynomials
More informationSemidefinite Relaxations Approach to Polynomial Optimization and One Application. Li Wang
Semidefinite Relaxations Approach to Polynomial Optimization and One Application Li Wang University of California, San Diego Institute for Computational and Experimental Research in Mathematics September
More informationLecture 5. 1 Goermans-Williamson Algorithm for the maxcut problem
Math 280 Geometric and Algebraic Ideas in Optimization April 26, 2010 Lecture 5 Lecturer: Jesús A De Loera Scribe: Huy-Dung Han, Fabio Lapiccirella 1 Goermans-Williamson Algorithm for the maxcut problem
More informationCOURSE ON LMI PART I.2 GEOMETRY OF LMI SETS. Didier HENRION henrion
COURSE ON LMI PART I.2 GEOMETRY OF LMI SETS Didier HENRION www.laas.fr/ henrion October 2006 Geometry of LMI sets Given symmetric matrices F i we want to characterize the shape in R n of the LMI set F
More informationHW2 - Due 01/30. Each answer must be mathematically justified. Don t forget your name.
HW2 - Due 0/30 Each answer must be mathematically justified. Don t forget your name. Problem. Use the row reduction algorithm to find the inverse of the matrix 0 0, 2 3 5 if it exists. Double check your
More informationThe moment-lp and moment-sos approaches
The moment-lp and moment-sos approaches LAAS-CNRS and Institute of Mathematics, Toulouse, France CIRM, November 2013 Semidefinite Programming Why polynomial optimization? LP- and SDP- CERTIFICATES of POSITIVITY
More informationSemidefinite Programming
Semidefinite Programming Notes by Bernd Sturmfels for the lecture on June 26, 208, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra The transition from linear algebra to nonlinear algebra has
More informationPolynomials, Ideals, and Gröbner Bases
Polynomials, Ideals, and Gröbner Bases Notes by Bernd Sturmfels for the lecture on April 10, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra We fix a field K. Some examples of fields
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 2
Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 2 Instructor: Farid Alizadeh Scribe: Xuan Li 9/17/2001 1 Overview We survey the basic notions of cones and cone-lp and give several
More informationDetecting global optimality and extracting solutions in GloptiPoly
Detecting global optimality and extracting solutions in GloptiPoly Didier HENRION 1,2 Jean-Bernard LASSERRE 1 1 LAAS-CNRS Toulouse 2 ÚTIA-AVČR Prague Part 1 Description of GloptiPoly Brief description
More informationFILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS.
FILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS. Let A be a ring, for simplicity assumed commutative. A filtering, or filtration, of an A module M means a descending sequence of submodules M = M 0
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationAlgebraic Varieties. Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra
Algebraic Varieties Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra Algebraic varieties represent solutions of a system of polynomial
More informationREPRESENTATION THEORY. WEEKS 10 11
REPRESENTATION THEORY. WEEKS 10 11 1. Representations of quivers I follow here Crawley-Boevey lectures trying to give more details concerning extensions and exact sequences. A quiver is an oriented graph.
More informationDistance Geometry-Matrix Completion
Christos Konaxis Algs in Struct BioInfo 2010 Outline Outline Tertiary structure Incomplete data Tertiary structure Measure diffrence of matched sets Def. Root Mean Square Deviation RMSD = 1 n x i y i 2,
More informationContinuous Optimisation, Chpt 9: Semidefinite Optimisation
Continuous Optimisation, Chpt 9: Semidefinite Optimisation Peter J.C. Dickinson DMMP, University of Twente p.j.c.dickinson@utwente.nl http://dickinson.website/teaching/2017co.html version: 28/11/17 Monday
More informationTopics in linear algebra
Chapter 6 Topics in linear algebra 6.1 Change of basis I want to remind you of one of the basic ideas in linear algebra: change of basis. Let F be a field, V and W be finite dimensional vector spaces over
More informationAn explicit construction of distinguished representations of polynomials nonnegative over finite sets
An explicit construction of distinguished representations of polynomials nonnegative over finite sets Pablo A. Parrilo Automatic Control Laboratory Swiss Federal Institute of Technology Physikstrasse 3
More informationA NEW SEMIDEFINITE PROGRAMMING HIERARCHY FOR CYCLES IN BINARY MATROIDS AND CUTS IN GRAPHS
A NEW SEMIDEFINITE PROGRAMMING HIERARCHY FOR CYCLES IN BINARY MATROIDS AND CUTS IN GRAPHS JOÃO GOUVEIA, MONIQUE LAURENT, PABLO A. PARRILO, AND REKHA THOMAS Abstract. The theta bodies of a polynomial ideal
More informationGlobal Optimization with Polynomials
Global Optimization with Polynomials Deren Han Singapore-MIT Alliance National University of Singapore Singapore 119260 Email: smahdr@nus.edu.sg Abstract The class of POP (Polynomial Optimization Problems)
More informationStructure of rings. Chapter Algebras
Chapter 5 Structure of rings 5.1 Algebras It is time to introduce the notion of an algebra over a commutative ring. So let R be a commutative ring. An R-algebra is a ring A (unital as always) together
More information6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC
6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC 2003 2003.09.02.10 6. The Positivstellensatz Basic semialgebraic sets Semialgebraic sets Tarski-Seidenberg and quantifier elimination Feasibility
More informationA new look at nonnegativity on closed sets
A new look at nonnegativity on closed sets LAAS-CNRS and Institute of Mathematics, Toulouse, France IPAM, UCLA September 2010 Positivstellensatze for semi-algebraic sets K R n from the knowledge of defining
More informationNon-commutative polynomial optimization
Non-commutative polynomial optimization S. Pironio 1, M. Navascuès 2, A. Acín 3 1 Laboratoire d Information Quantique (Brussels) 2 Department of Mathematics (Bristol) 3 ICFO-The Institute of Photonic Sciences
More informationLecture Note 5: Semidefinite Programming for Stability Analysis
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State
More informationHow to generate weakly infeasible semidefinite programs via Lasserre s relaxations for polynomial optimization
CS-11-01 How to generate weakly infeasible semidefinite programs via Lasserre s relaxations for polynomial optimization Hayato Waki Department of Computer Science, The University of Electro-Communications
More informationIntroduction to Semidefinite Programming I: Basic properties a
Introduction to Semidefinite Programming I: Basic properties and variations on the Goemans-Williamson approximation algorithm for max-cut MFO seminar on Semidefinite Programming May 30, 2010 Semidefinite
More informationGraph invariants in the spin model
Graph invariants in the spin model Alexander Schrijver 1 Abstract. Given a symmetric n n matrix A, we define, for any graph G, f A (G) := a φ(u),φ(v). φ:v G {1,...,n} uv EG We characterize for which graph
More informationLecture 2: Gröbner Basis and SAGBI Basis
Lecture 2: Gröbner Basis and SAGBI Basis Mohammed Tessema Suppose we have a graph. Suppose we color the graph s vertices with 3 colors so that if the vertices are adjacent they are not the same colors.
More informationAlgebraic Optimization with Applications to System Theory
VRIJE UNIVERSITEIT Algebraic Optimization with Applications to System Theory ACADEMISCH PROEFSCHRIFT ter verkrijging van de graad van doctor aan de Vrije Universiteit Amsterdam, op gezag van de rector
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 3: Positive-Definite Systems; Cholesky Factorization Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 11 Symmetric
More informationQR Decomposition. When solving an overdetermined system by projection (or a least squares solution) often the following method is used:
(In practice not Gram-Schmidt, but another process Householder Transformations are used.) QR Decomposition When solving an overdetermined system by projection (or a least squares solution) often the following
More information12. Hilbert Polynomials and Bézout s Theorem
12. Hilbert Polynomials and Bézout s Theorem 95 12. Hilbert Polynomials and Bézout s Theorem After our study of smooth cubic surfaces in the last chapter, let us now come back to the general theory of
More informationMoment Matrices, Trace Matrices and the Radical of Ideals
Moment Matrices, Trace Matrices and the Radical of Ideals Itnuit Janovitz-Freireich, Ágnes Szántó North Carolina State University Campus Box 805 Raleigh, NC, 95 USA {ijanovi, aszanto}@ncsu.edu Bernard
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More information8. Prime Factorization and Primary Decompositions
70 Andreas Gathmann 8. Prime Factorization and Primary Decompositions 13 When it comes to actual computations, Euclidean domains (or more generally principal ideal domains) are probably the nicest rings
More informationUnderstanding and Implementing F5
Understanding and Implementing F5 John Perry john.perry@usm.edu University of Southern Mississippi Understanding and Implementing F5 p.1 Overview Understanding F5 Description Criteria Proofs Implementing
More informationNo OPTIMIZATION OF UNIVARIATE FUNCTIONS ON BOUNDED INTERVALS BY INTERPOLATION AND SEMIDEFINITE PROGRAMMING
No. 2006 26 OPTIMIZATION OF UNIVARIATE FUNCTIONS ON BOUNDED INTERVALS BY INTERPOLATION AND SEMIDEFINITE PROGRAMMING By E. de Klerk, G. Elabwabi, D. den Hertog April 2006 ISSN 0924-7815 Optimization of
More informationRanks of Real Symmetric Tensors
Ranks of Real Symmetric Tensors Greg Blekherman SIAM AG 2013 Algebraic Geometry of Tensor Decompositions Real Symmetric Tensor Decompositions Let f be a form of degree d in R[x 1,..., x n ]. We would like
More informationOptimization of univariate functions on bounded intervals by interpolation and semidefinite programming
Optimization of univariate functions on bounded intervals by interpolation and semidefinite programming E. de Klerk G. Elabwabi D. den Hertog April 24, 2006 Abstract We consider the problem of minimizing
More informationReview problems for MA 54, Fall 2004.
Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on
More informationReal Symmetric Matrices and Semidefinite Programming
Real Symmetric Matrices and Semidefinite Programming Tatsiana Maskalevich Abstract Symmetric real matrices attain an important property stating that all their eigenvalues are real. This gives rise to many
More informationTropical Algebraic Geometry 3
Tropical Algebraic Geometry 3 1 Monomial Maps solutions of binomial systems an illustrative example 2 The Balancing Condition balancing a polyhedral fan the structure theorem 3 The Fundamental Theorem
More informationMath Matrix Algebra
Math 44 - Matrix Algebra Review notes - (Alberto Bressan, Spring 7) sec: Orthogonal diagonalization of symmetric matrices When we seek to diagonalize a general n n matrix A, two difficulties may arise:
More informationToric Ideals, an Introduction
The 20th National School on Algebra: DISCRETE INVARIANTS IN COMMUTATIVE ALGEBRA AND IN ALGEBRAIC GEOMETRY Mangalia, Romania, September 2-8, 2012 Hara Charalambous Department of Mathematics Aristotle University
More informationMoments and Positive Polynomials for Optimization II: LP- VERSUS SDP-relaxations
Moments and Positive Polynomials for Optimization II: LP- VERSUS SDP-relaxations LAAS-CNRS and Institute of Mathematics, Toulouse, France Tutorial, IMS, Singapore 2012 LP-relaxations LP- VERSUS SDP-relaxations
More informationAlgebra II. Paulius Drungilas and Jonas Jankauskas
Algebra II Paulius Drungilas and Jonas Jankauskas Contents 1. Quadratic forms 3 What is quadratic form? 3 Change of variables. 3 Equivalence of quadratic forms. 4 Canonical form. 4 Normal form. 7 Positive
More informationSome notes on Coxeter groups
Some notes on Coxeter groups Brooks Roberts November 28, 2017 CONTENTS 1 Contents 1 Sources 2 2 Reflections 3 3 The orthogonal group 7 4 Finite subgroups in two dimensions 9 5 Finite subgroups in three
More informationMCS 563 Spring 2014 Analytic Symbolic Computation Friday 31 January. Quotient Rings
Quotient Rings In this note we consider again ideals, but here we do not start from polynomials, but from a finite set of points. The application in statistics and the pseudo code of the Buchberger-Möller
More informationSince G is a compact Lie group, we can apply Schur orthogonality to see that G χ π (g) 2 dg =
Problem 1 Show that if π is an irreducible representation of a compact lie group G then π is also irreducible. Give an example of a G and π such that π = π, and another for which π π. Is this true for
More informationFlat extension and ideal projection
Flat extension and ideal projection Thomas Kunkle Department of Mathematics, College of Charleston, Charleston, SC 29424 0001 Abstract A generalization of the flat extension theorems of Curto and Fialkow
More informationLinGloss. A glossary of linear algebra
LinGloss A glossary of linear algebra Contents: Decompositions Types of Matrices Theorems Other objects? Quasi-triangular A matrix A is quasi-triangular iff it is a triangular matrix except its diagonal
More informationSemidefinite Programming
Semidefinite Programming Basics and SOS Fernando Mário de Oliveira Filho Campos do Jordão, 2 November 23 Available at: www.ime.usp.br/~fmario under talks Conic programming V is a real vector space h, i
More informationSolutions of exercise sheet 8
D-MATH Algebra I HS 14 Prof. Emmanuel Kowalski Solutions of exercise sheet 8 1. In this exercise, we will give a characterization for solvable groups using commutator subgroups. See last semester s (Algebra
More information2. Intersection Multiplicities
2. Intersection Multiplicities 11 2. Intersection Multiplicities Let us start our study of curves by introducing the concept of intersection multiplicity, which will be central throughout these notes.
More informationLecture 8 : Eigenvalues and Eigenvectors
CPS290: Algorithmic Foundations of Data Science February 24, 2017 Lecture 8 : Eigenvalues and Eigenvectors Lecturer: Kamesh Munagala Scribe: Kamesh Munagala Hermitian Matrices It is simpler to begin with
More information3.1. Derivations. Let A be a commutative k-algebra. Let M be a left A-module. A derivation of A in M is a linear map D : A M such that
ALGEBRAIC GROUPS 33 3. Lie algebras Now we introduce the Lie algebra of an algebraic group. First, we need to do some more algebraic geometry to understand the tangent space to an algebraic variety at
More informationRational Univariate Representation
Rational Univariate Representation 1 Stickelberger s Theorem a rational univariate representation (RUR) 2 The Elbow Manipulator a spatial robot arm with three links 3 Application of the Newton Identities
More informationRational Sums of Squares and Applications
f ΣR[X] 2 " f ΣQ[X] 2 Rational Sums of Squares and Applications Christopher Hillar (MSRI & Berkeley) A 2008 study found that adding a picture of a brain scan to a scientific argument about human nature
More informationarxiv: v1 [cs.sc] 29 Nov 2008
Moment Matrices, Trace Matrices and the Radical of Ideals arxiv:0812.0088v1 [cs.sc] 29 Nov 2008 I. Janovitz-Freireich & A. Szántó & B. Mourrain & L. Rónyai November 5, 2018 Abstract Let f 1,...,f s K[x
More informationSEMIDEFINITE PROGRAM BASICS. Contents
SEMIDEFINITE PROGRAM BASICS BRIAN AXELROD Abstract. A introduction to the basics of Semidefinite programs. Contents 1. Definitions and Preliminaries 1 1.1. Linear Algebra 1 1.2. Convex Analysis (on R n
More information5 Quiver Representations
5 Quiver Representations 5. Problems Problem 5.. Field embeddings. Recall that k(y,..., y m ) denotes the field of rational functions of y,..., y m over a field k. Let f : k[x,..., x n ] k(y,..., y m )
More informationx i e i ) + Q(x n e n ) + ( i<n c ij x i x j
Math 210A. Quadratic spaces over R 1. Algebraic preliminaries Let V be a finite free module over a nonzero commutative ring F. Recall that a quadratic form on V is a map Q : V F such that Q(cv) = c 2 Q(v)
More informationAlgebra Workshops 10 and 11
Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on
More informationCopositive Plus Matrices
Copositive Plus Matrices Willemieke van Vliet Master Thesis in Applied Mathematics October 2011 Copositive Plus Matrices Summary In this report we discuss the set of copositive plus matrices and their
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5
Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5 Instructor: Farid Alizadeh Scribe: Anton Riabov 10/08/2001 1 Overview We continue studying the maximum eigenvalue SDP, and generalize
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 16: Eigenvalue Problems; Similarity Transformations Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 18 Eigenvalue
More informationThe moment-lp and moment-sos approaches in optimization
The moment-lp and moment-sos approaches in optimization LAAS-CNRS and Institute of Mathematics, Toulouse, France Workshop Linear Matrix Inequalities, Semidefinite Programming and Quantum Information Theory
More informationEigenvectors. Prop-Defn
Eigenvectors Aim lecture: The simplest T -invariant subspaces are 1-dim & these give rise to the theory of eigenvectors. To compute these we introduce the similarity invariant, the characteristic polynomial.
More informationMa/CS 6b Class 23: Eigenvalues in Regular Graphs
Ma/CS 6b Class 3: Eigenvalues in Regular Graphs By Adam Sheffer Recall: The Spectrum of a Graph Consider a graph G = V, E and let A be the adjacency matrix of G. The eigenvalues of G are the eigenvalues
More informationA Field Extension as a Vector Space
Chapter 8 A Field Extension as a Vector Space In this chapter, we take a closer look at a finite extension from the point of view that is a vector space over. It is clear, for instance, that any is a linear
More informationExplicit Sensor Network Localization using Semidefinite Programming and Clique Reductions
Explicit Sensor Network Localization using Semidefinite Programming and Clique Reductions Nathan Krislock and Henry Wolkowicz Dept. of Combinatorics and Optimization University of Waterloo Haifa Matrix
More information2.1 Laplacian Variants
-3 MS&E 337: Spectral Graph heory and Algorithmic Applications Spring 2015 Lecturer: Prof. Amin Saberi Lecture 2-3: 4/7/2015 Scribe: Simon Anastasiadis and Nolan Skochdopole Disclaimer: hese notes have
More informationExact algorithms for linear matrix inequalities
0 / 12 Exact algorithms for linear matrix inequalities Mohab Safey El Din (UPMC/CNRS/IUF/INRIA PolSys Team) Joint work with Didier Henrion, CNRS LAAS (Toulouse) Simone Naldi, Technische Universität Dortmund
More information4. Algebra and Duality
4-1 Algebra and Duality P. Parrilo and S. Lall, CDC 2003 2003.12.07.01 4. Algebra and Duality Example: non-convex polynomial optimization Weak duality and duality gap The dual is not intrinsic The cone
More informationMinimum Polynomials of Linear Transformations
Minimum Polynomials of Linear Transformations Spencer De Chenne University of Puget Sound 30 April 2014 Table of Contents Polynomial Basics Endomorphisms Minimum Polynomial Building Linear Transformations
More informationThe Pythagoras numbers of projective varieties.
The Pythagoras numbers of projective varieties. Grigoriy Blekherman (Georgia Tech, USA) Rainer Sinn (FU Berlin, Germany) Gregory G. Smith (Queen s University, Canada) Mauricio Velasco (Universidad de los
More informationOperator-valued Herglotz kernels and functions of positive real
Operator-valued Herglotz kernels and functions of positive real part on the ball University of Florida July 24, 2008 Let D = {z C : z < 1}. Theorem (Riesz-Herglotz) Let f Hol(D). Then Rf 0 in D iff a postitive
More informationMath (P)refresher Lecture 8: Unconstrained Optimization
Math (P)refresher Lecture 8: Unconstrained Optimization September 2006 Today s Topics : Quadratic Forms Definiteness of Quadratic Forms Maxima and Minima in R n First Order Conditions Second Order Conditions
More informationFinal Report on the Newton Institute Programme Polynomial Optimisation
Final Report on the Newton Institute Programme Polynomial Optimisation July 15th to August 19th 2013 1 Committees The organising committee consisted of: Adam N. Letchford (Department of Management Science,
More informationLecture 2: Linear operators
Lecture 2: Linear operators Rajat Mittal IIT Kanpur The mathematical formulation of Quantum computing requires vector spaces and linear operators So, we need to be comfortable with linear algebra to study
More information4 Hilbert s Basis Theorem and Gröbner basis
4 Hilbert s Basis Theorem and Gröbner basis We define Gröbner bases of ideals in multivariate polynomial rings and see how they work in tandem with the division algorithm. We look again at the standard
More information(1) A frac = b : a, b A, b 0. We can define addition and multiplication of fractions as we normally would. a b + c d
The Algebraic Method 0.1. Integral Domains. Emmy Noether and others quickly realized that the classical algebraic number theory of Dedekind could be abstracted completely. In particular, rings of integers
More information10. Smooth Varieties. 82 Andreas Gathmann
82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It
More information