Nonnegative Polynomials and Sums of Squares

Grigoriy Blekherman May 17, 2011 SIAM Optimization

Nonnegative Polynomials and Sums of Squares Accentuate the positive; delete the negative.

Grigoriy Blekherman May 17, 2011 SIAM Optimization

Dramatis Personae • A multivariate polynomial p is called nonnegative (psd) if

p(x) ≥ 0 for all x ∈ Rn . • A polynomial p is a sum of squares (sos) if we can write

p(x) =

P

qi2 for some polynomials qi .

• Sums of squares are clearly nonnegative from their

presentation. Abstract Question: Can we always write a nonnegative polynomial in a way that makes its nonnegativity apparent? Practical Version: Can we efficiently compute such representations?

The First Date

• If a polynomial p is nonnegative then we can make it

homogeneous and it will remain nonnegative. • n number of variables • 2d degree • R[x]2d vector space of homogeneous polynomials (forms) in n

variables, of degree 2d. Nonnegative polynomials and sums of squares form full dimensional closed convex cones Pn,2d and Σn,2d in R[x]2d .

Hilbert’s Theorem

Hilbert’s Theorem Pn,2d = Σn,2d in the following three cases: n = 2 (univariate non-homogeneous case), 2d = 2 (quadratic forms), and n = 3, 2d = 4 (ternary quartics). In all other cases there exist nonnegative forms that are not sums of squares. The standard proof of equality cases provides a separate argument for each case. Today we will see a unified argument and applications to the 17th problem in 3 variables based on studying the dual cone Σ∗n,2d and Gorenstein ideals.

Hilbert’s 17th Problem Is it true that we can write every nonnegative polynomial as a sum of squares of rational functions? YES (Artin,Schreier 1920’s). Equivalently, there exists h ∈ Σn,2m such that ph is a sum of squares. However, the multiplier h may have very large degree.

Hilbert’s 17th Problem Is it true that we can write every nonnegative polynomial as a sum of squares of rational functions? YES (Artin,Schreier 1920’s). Equivalently, there exists h ∈ Σn,2m such that ph is a sum of squares. However, the multiplier h may have very large degree.

Hilbert’s Ternary Forms Theorem Suppose that p ∈ P3,2d . Then there exists q ∈ P3,2d−4 such that pq is a sum of squares. To obtain a sum of squares multiplier, apply the Theorem to q and continue knocking down the degree.

The Dual Cone to Sums of Squares If Σn,2d 6= Pn,2d then there exists a linear functional ` : R[x]2d → R such that `(q 2 ) ≥ 0,

so that ` is supporting to Σn,2d

but ` cuts through the interior of Pn,2d . How to construct linear functionals that are nonnegative on squares? Point evaluations work: for v ∈ Rn , let `v (f ) = f (v ). But point evaluations are also supporting to Pn,2d .

Explicit Example Take 9 points in R3 : {v1 , . . . , v9 } = {(x, y , 1)} with x, y ∈ {−1, 0, 1}. Set v9 = (0, 0, 1). Then for any sum of squares p ∈ Σ3,6 we have `(p) = p(v1 ) + . . . + p(v8 ) −

4p(v9 ) ≥ 0. 5

The points vi come from intersection of two real cubics: q1 = x(x − z)(x + z) and q2 = y (y − z)(y + z). Let p = z(3x 2 + 3y 2 − 5z 2 ). Then `(q12 + q22 + p 2 ) = `(q12 ) + `(q22 ) + `(p 2 ) = 0, but q12 + q22 + p 2 is a strictly positive form.

Gorenstein Ideals To a linear functional ` in the dual space R[x]∗2d we can associate the Gorenstein Ideal I (`) with socle `. I (`) consists of all forms of degree greater than 2d and forms q of degree at most 2d, such that `(rq) = 0 for all forms r of degree 2d − deg q.

Gorenstein Ideals To a linear functional ` in the dual space R[x]∗2d we can associate the Gorenstein Ideal I (`) with socle `. I (`) consists of all forms of degree greater than 2d and forms q of degree at most 2d, such that `(rq) = 0 for all forms r of degree 2d − deg q. Easy Example: If we take `v (f ) = f (v ) then the Gorenstein Ideal I (`v ) of `v consists of all forms of degree at most 2d that have a zero at v . We are interested in positive socles: `(p 2 ) ≥ 0

for all

p ∈ R[x]d .

No Point Evaluations We only want the supporting hyperplanes that intersect Σn,2d maximally. This translates to positive Gorenstein ideals in which Id generates I2d . We are also not interested in functionals ` that come from point evaluations. This translates to The forms in I (`)d have no common zeroes over C.

No Point Evaluations We only want the supporting hyperplanes that intersect Σn,2d maximally. This translates to positive Gorenstein ideals in which Id generates I2d . We are also not interested in functionals ` that come from point evaluations. This translates to The forms in I (`)d have no common zeroes over C. Ternary quartics: the forms in I (`)2 have no common zeroes. There exist q1 , q2 and q3 in I (`)d with no common zeroes. Take the ideal J =< q1 , q2 , q3 > generated by qi . The dimension of I (`)4 is already 3 ∗ 6 − 32 = 15. This is the dimension of R[x]3,4 , so we have generated too much.

Border Cases For the cases n = 3, 2d = 6 and n = 4, 2d = 4 we get the following dimension counts: 3 cubics in R[x]3,3 with no common zeroes generate in degree 6  3 ideal of dimension 3 ∗ 10 − 2 = 27, which is exactly 1 less than dimension of R[x]3,6 . Then can construct an appropriate linear functional. 4 quadratics in R[x]4,2 with no common zeroes generate in degree 4 ideal of dimension 4 ∗ 10 − 42 = 34, which is exactly 1 less than dimension of R[x]4,4 . Then can construct an appropriate linear functional.

”Positive” Results Equivalent Reformulation of Hilbert’s Ternary Forms Theorem: ∗ Suppose that a linear functional ` ∈ H3,4d defines a positive Gorenstein ideal I (`). Then I (`) does not contain any strictly positive forms of degree 2d + 2.

Main idea of proof: Assume not. We can define the linear functional ` by evaluation on zeroes of p. Since p has no real zeroes, we have to evaluate on complex zeroes. It is possible to limit the number of points at which we need to evaluate. A dimension count yields that a linear functional ` defined only on a few complex points cannot be positive. Is the bound tight? Can show that there exist a nonnegative form p of degree 4k + 6 so that pq is not a sum of squares for all q of degree 2k.

THANK YOU!

Nonnegative Polynomials and Sums of Squares Why ...

May 17, 2011 - variables, of degree 2d. Nonnegative polynomials and sums of squares form full dimensional closed convex cones Pn,2d and Σn,2d in R[x]2d .

214KB Sizes 1 Downloads 217 Views

Recommend Documents

Nonnegative Polynomials and Sums of Squares Why ...
Jan 24, 2011 - Why do bad things happen to good polynomials? Grigoriy .... There exist constants c1(d) and c2(d), dependent on the degree d only, such that.

Nonnegative Polynomials and Sums of Squares Why ...
Mar 17, 2011 - We can compute instead the best sos lower bound: γ. ∗. = max f −γ is sos ... For small n and 2d find good bounds on degree of q and explicit.

NONNEGATIVE MATRIX FACTORIZATION AND SPATIAL ...
ABSTRACT. We address the problem of blind audio source separation in the under-determined and convolutive case. The contribution of each source to the mixture channels in the time-frequency domain is modeled by a zero-mean Gaussian random vector with

SUMS OF KLOOSTERMAN SUMS OVER ARITHMETIC ...
Sep 24, 2010 - 10. SATADAL GANGULY AND JYOTI SENGUPTA. 2.3. Special functions. The importance of Bessel functions in the theory of automorphic forms can be gauged from the result of Sears and Titchmarsh. (see [ST09] or Chapter 16, [IK04]) which says

A Generalization of Riemann Sums
For every continuous function f on the interval [0,1], lim n→∞. 1 nα n. ∑ k=1 f .... 4. ∫ 1. 0 dt. 1 + t. − I. Hence, I = π. 8 log 2. Replacing back in (5) we obtain (1).

Proof Without Words: Perfect Numbers and Sums of ...
Theorem. Every even perfect number Np = 2p−1(2p − 1) with p ≥ 3 prime is the sum of the first n odd cubes for n = 2(p−1)/2, i.e., Np = 13 + 33 +···+ (2n − 1)3 [1].

Toward Faster Nonnegative Matrix Factorization: A New Algorithm and ...
College of Computing, Georgia Institute of Technology. Atlanta, GA ..... Otherwise, a complementary ba- ...... In Advances in Neural Information Pro- cessing ...

Sums of distances in normed spaces
where the sum is taken over all integers i,j, satisfying 1 < i < j < r. Let. S = { x: I[x[[ = 1} be the unit sphere of X. Martelli and Busenberg [8] use inequalities in ...

Difference of two squares - presentations.pdf
Page 1 of 9. Page 1 of 9. Page 2 of 9. Page 2 of 9. Page 3 of 9. Page 3 of 9. Page 4 of 9. Page 4 of 9. Difference of two squares - presentations.pdf. Difference of two squares - presentations.pdf. Open. Extract. Open with. Sign In. Main menu. Displa

COLLECTION SHT_25 squares
Fill up this sheet with 25 BOX TOPS. (no campbells labels or “bonus box tops”, please!) & you'll be donating this much cash without even opening your wallet! +.

Difference of two squares - presentations.pdf
There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. Difference of two ...

COLLECTION SHT_25 squares
Thank You! TEACHER________________________ STUDENT'S NAME: Fill up this sheet with 25 BOX TOPS. (no campbells labels or “bonus box tops”, please!) & you'll be donating this much cash without even opening your wallet! +. +. +. +. = $.50. +. +. +.

Multiplying Polynomials "Why is a stick of gum like a sneeze".pdf ...
Multiplying Polynomials "Why is a stick of gum like a sneeze".pdf. Multiplying Polynomials "Why is a stick of gum like a sneeze".pdf. Open. Extract. Open with.

10.1 Adding and Subtracting Polynomials
Apr 6, 2014 - Vocabulary: polynomial degree of polynomial monomial standard form leading coefficient binomial degree trinomial. Example 1 Identifying Polynomial Coefficients. What are coefficients? Name ALL the coefficients. Rewrite the polynomials i

12.11 Applications of Taylor Polynomials
Goal: Approximate a Function with a Taylor Polynomial ó Why? ◦ Polynomials are the easiest functions to work with ó How? ó How? ◦ Physics & Engineering.

Dividing Polynomials
A2. I AM. ID: 1. Dividing Polynomials. Date. Block. Divide. 1) (5n3 + 3n? + 2n) + 6n? 2) (4x + 32x+ + 2x3) + 8x. 2. 3) (2k + 12kº + 5k) + 4k? 4) (2x + 4x + 16x4) + 4x3. 5) (k® +k? – 20k + 22) + (k - 3). 6) (2x + 5x2 + 8x + 10) + (x + 2). 7) (a3 -

Equimodular Polynomials and the Tritangency ...
puts the natural topology of Cn on the set of monic polynomials of degree n. When- ..... ALEX RYBA received his B.A. and Ph.D. from the University of Cambridge. ... Department of Computer Science, Queens College, Flushing NY 11367.

The Magic of Magic Squares
third order magic square is similar to the nine points of inflection on an ... Some mathematicians call this formula “the most remarkable formula in mathemat-.

Representations of Orthogonal Polynomials
classical orthogonal polynomials if it is the solution of a di erence equation of the ... in the discrete case, in terms of the coe cients a; b; c; d and e of the given di ...

ASYMPTOTICS OF CHEBYSHEV POLYNOMIALS, I ...
Mar 25, 2015 - argument shows that for any γ ∈ (−Tn e,Tn e) all n solutions of Tn(x) = γ are ... there is a probability measure whose Coulomb energy is R(e). Since ... The Green's function, Ge(z), of a compact subset, e ⊂ C, is defined by.

Inscribed Squares
Dec 6, 2004 - on a sideline of ABC can be constructed easily from a homothety of a square erected on the side of ABC. Consider the two squares erected on the side BC. Their centers are the points with homogeneous barycentric coordinates (−a2 : SC +

Estimates for sums of eigenvalues of the Laplacian
Now suppose that that there exists a bi-Lischitz map f which maps Ω onto an open ball B in Rn. Let CΩ be a .... (Our definition of F is motivated by the definition of.