Recovery of algebraic-exponential data from moments Jean B. Lasserre LAAS-CNRS and Institute of Mathematics, Toulouse, France ICERM, Brown University, June 2014 ⋆ Part of this work is joint with M. Putinar Jean B. Lasserre Recovery of algebraic-exponential data from moments
Motivation An important property of Positively Homogeneous Functions (PHF) Some properties (convexity, polarity) Sub-level sets of minimum volume containing K Exact reconstruction from moments Recovery of the defining function of a semi-algebraic set Jean B. Lasserre Recovery of algebraic-exponential data from moments
Motivation An important property of Positively Homogeneous Functions (PHF) Some properties (convexity, polarity) Sub-level sets of minimum volume containing K Exact reconstruction from moments Recovery of the defining function of a semi-algebraic set Jean B. Lasserre Recovery of algebraic-exponential data from moments
Motivation An important property of Positively Homogeneous Functions (PHF) Some properties (convexity, polarity) Sub-level sets of minimum volume containing K Exact reconstruction from moments Recovery of the defining function of a semi-algebraic set Jean B. Lasserre Recovery of algebraic-exponential data from moments
Motivation An important property of Positively Homogeneous Functions (PHF) Some properties (convexity, polarity) Sub-level sets of minimum volume containing K Exact reconstruction from moments Recovery of the defining function of a semi-algebraic set Jean B. Lasserre Recovery of algebraic-exponential data from moments
Motivation An important property of Positively Homogeneous Functions (PHF) Some properties (convexity, polarity) Sub-level sets of minimum volume containing K Exact reconstruction from moments Recovery of the defining function of a semi-algebraic set Jean B. Lasserre Recovery of algebraic-exponential data from moments
Motivation An important property of Positively Homogeneous Functions (PHF) Some properties (convexity, polarity) Sub-level sets of minimum volume containing K Exact reconstruction from moments Recovery of the defining function of a semi-algebraic set Jean B. Lasserre Recovery of algebraic-exponential data from moments
Exact reconstruction Reconstruction of a shape K ⊂ R n (convex or not) from knowledge of finitely many moments � x α 1 1 · · · x α n α ∈ N n y α = n dx , d , K for some integer d , is a difficult and challenging problem! EXACT recovery of K from y = ( y α ) , α ∈ N n d , is even more difficult and challenging! Jean B. Lasserre Recovery of algebraic-exponential data from moments
Exact reconstruction Reconstruction of a shape K ⊂ R n (convex or not) from knowledge of finitely many moments � x α 1 1 · · · x α n α ∈ N n y α = n dx , d , K for some integer d , is a difficult and challenging problem! EXACT recovery of K from y = ( y α ) , α ∈ N n d , is even more difficult and challenging! Jean B. Lasserre Recovery of algebraic-exponential data from moments
Exact recovery (continued) Examples of exact recovery: Quadrature (planar) Domains in ( R 2 ) (Gustafsson, He, Milanfar and Putinar (Inverse Problems, 2000)) • via an exponential transform Convex Polytopes (in R n ) (Gravin, L., Pasechnik and Robins (Discrete & Comput. Geometry (2012)) • Use Brion-Barvinok-Khovanski-Lawrence-Pukhlikov � � c , x � j dx combined with moment formula for projections P a Prony-type method to recover the vertices of P . and extension to Non convex polyhedra by Pasechnik et al. • via inversion of Fantappié transform Jean B. Lasserre Recovery of algebraic-exponential data from moments
Exact recovery (continued) Examples of exact recovery: Quadrature (planar) Domains in ( R 2 ) (Gustafsson, He, Milanfar and Putinar (Inverse Problems, 2000)) • via an exponential transform Convex Polytopes (in R n ) (Gravin, L., Pasechnik and Robins (Discrete & Comput. Geometry (2012)) • Use Brion-Barvinok-Khovanski-Lawrence-Pukhlikov � � c , x � j dx combined with moment formula for projections P a Prony-type method to recover the vertices of P . and extension to Non convex polyhedra by Pasechnik et al. • via inversion of Fantappié transform Jean B. Lasserre Recovery of algebraic-exponential data from moments
Approximate recovery can de done in multi-dimensions (Cuyt, Golub, Milanfar and Verdonk, 2005) via : (multi-dimensional versions of) homogeneous Padé approximants applied to the Stieltjes transform. cubature formula at each point of grid solving a linear system of equations to retrieve the indicator function of K Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: I Exact recovery. K = { x ∈ R n : g ( x ) ≤ 1 } has finite Lebesgue volume. g is a nonnegative homogeneous polynomial Data are finitely many moments: � x α d x , α ∈ N n y α = d . K • Also works for Quasi-homogeneous polynomials, i.e., when g ( λ u 1 x 1 , . . . , λ u n x n ) = λ g ( x ) , x ∈ R n , λ > 0 for some vector u ∈ Q n . ( d -Homogeneous =u-quasi homogeneous with u i = 1 d for all i ). Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: I Exact recovery. K = { x ∈ R n : g ( x ) ≤ 1 } has finite Lebesgue volume. g is a nonnegative homogeneous polynomial Data are finitely many moments: � x α d x , α ∈ N n y α = d . K • Also works for Quasi-homogeneous polynomials, i.e., when g ( λ u 1 x 1 , . . . , λ u n x n ) = λ g ( x ) , x ∈ R n , λ > 0 for some vector u ∈ Q n . ( d -Homogeneous =u-quasi homogeneous with u i = 1 d for all i ). Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: I Exact recovery. K = { x ∈ R n : g ( x ) ≤ 1 } has finite Lebesgue volume. g is a nonnegative homogeneous polynomial Data are finitely many moments: � x α d x , α ∈ N n y α = d . K • Also works for Quasi-homogeneous polynomials, i.e., when g ( λ u 1 x 1 , . . . , λ u n x n ) = λ g ( x ) , x ∈ R n , λ > 0 for some vector u ∈ Q n . ( d -Homogeneous =u-quasi homogeneous with u i = 1 d for all i ). Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: I Exact recovery. K = { x ∈ R n : g ( x ) ≤ 1 } has finite Lebesgue volume. g is a nonnegative homogeneous polynomial Data are finitely many moments: � x α d x , α ∈ N n y α = d . K • Also works for Quasi-homogeneous polynomials, i.e., when g ( λ u 1 x 1 , . . . , λ u n x n ) = λ g ( x ) , x ∈ R n , λ > 0 for some vector u ∈ Q n . ( d -Homogeneous =u-quasi homogeneous with u i = 1 d for all i ). Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: I Exact recovery. K = { x ∈ R n : g ( x ) ≤ 1 } has finite Lebesgue volume. g is a nonnegative homogeneous polynomial Data are finitely many moments: � x α d x , α ∈ N n y α = d . K • Also works for Quasi-homogeneous polynomials, i.e., when g ( λ u 1 x 1 , . . . , λ u n x n ) = λ g ( x ) , x ∈ R n , λ > 0 for some vector u ∈ Q n . ( d -Homogeneous =u-quasi homogeneous with u i = 1 d for all i ). Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: II Exact recovery. G ⊂ R n is open with G = int G and with real algebraic boundary ∂ G . A polynomial of degree d vanishes on ∂ G . Data are finitely many moments: � x α d x , α ∈ N n y α = d . K Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: II Exact recovery. G ⊂ R n is open with G = int G and with real algebraic boundary ∂ G . A polynomial of degree d vanishes on ∂ G . Data are finitely many moments: � x α d x , α ∈ N n y α = d . K Jean B. Lasserre Recovery of algebraic-exponential data from moments
This talk: II Exact recovery. G ⊂ R n is open with G = int G and with real algebraic boundary ∂ G . A polynomial of degree d vanishes on ∂ G . Data are finitely many moments: � x α d x , α ∈ N n y α = d . K Jean B. Lasserre Recovery of algebraic-exponential data from moments
A little detour Positively Homogeneous functions (PHF) form a wide class of functions encountered in many applications. As a consequence of homogeneity, they enjoy very particular properties, and among them the celebrated and very useful Euler’s identity which allows to deduce additional properties of PHFs in various contexts. Another (apparently not well-known) property of PHFs yields surprising and unexpected results, some of them already known in particular cases. The case of homogeneous polynomials is even more interesting! Jean B. Lasserre Recovery of algebraic-exponential data from moments
A little detour Positively Homogeneous functions (PHF) form a wide class of functions encountered in many applications. As a consequence of homogeneity, they enjoy very particular properties, and among them the celebrated and very useful Euler’s identity which allows to deduce additional properties of PHFs in various contexts. Another (apparently not well-known) property of PHFs yields surprising and unexpected results, some of them already known in particular cases. The case of homogeneous polynomials is even more interesting! Jean B. Lasserre Recovery of algebraic-exponential data from moments
Recommend
More recommend