190 likes | 285 Views
Ellipsoid-type Confidential Bounds on Semi-algebraic Sets via SDP Relaxation . Makoto Yamashita Masakazu Kojima Tokyo Institute of Technology. SDP relaxation (convex region). SDP solution. Confidential Bounds in Polynomial Optimization Problem . Semi-algebraic Sets .
E N D
Ellipsoid-type Confidential Bounds on Semi-algebraic Sets via SDP Relaxation Makoto Yamashita Masakazu Kojima Tokyo Institute of Technology
SDP relaxation (convex region) SDP solution Confidential Bounds in Polynomial Optimization Problem Semi-algebraic Sets We compute this ellipsoid by SDP. Optimal solutions exist in this ellipsoid. (Polynomials) Feasible region Local adjustment for feasible region min Optimal
Ellipsoid research • . • MVEE (the minimum volume enclosing ellipsoid) • Our approach based on SDP relaxation • Solvable by SDP • Small computation cost⇒We can execute multiple times changing
Outline • Math Form of Ellipsoids • SDP relaxation • Examples of POP • Tightness of Ellipsoids
Mathematical Formulation • . • EllipsoidWe define • . By some steps, we consider SDP relaxation
Lifting • . • . • Note that • Furthermore ⇒ (convex hull) quadratic linear (easier) Still difficult
SDP relaxation • . • . relaxation
Inner minimization • . • . • Gradient • Optimal attained at • . • Cover
Relations of SDP
Example from POP • ex9_1_2 from GLOBAL library(http://www.gamsworld.org/global/global.htm) • We use SparsePOP to solve this by SDP relaxation SparsePOP http://sparsepop.sourceforge.net
Reduced POP Optimal Solutions:
Optimal Solutions: Ellipsoids for Reduced SDP Very tight bound
Results on POP • Very good objective values • ex_9_1_2 & ex_9_1_8 have multiple optimal solutions ⇒ large radius
Tightness of Ellipsoids • Target set • 6 Shape Matricies • We draw 2D picture ,
The case p=2 (2 constraints) • The ellipsoids are tight. 6 ellipsoids by SDP Target set
More constraints Ellipsoids shrink. But its speed is slower than the target set. p=2 p=32 p=128
Conclusion & Future works • An enclosing ellipsoid by SDP relaxation • Improve the SDP solution of POP • Very low computation cost • Successive ellipsoid for POP sometimes stops before bounding the region appropriately • Ellipsoids may become loose in the case of many constraints
Thank you very much for your attention. This talk is based on the following technical paper Masakazu Kojima and Makoto Yamashita, “Enclosing Ellipsoids and Elliptic Cylinders of Semialgebraic Sets and Their Application to Error Boundsin Polynomial Optimization”, Research Report B-459, Dept. of Math. and Comp. Sciences,Tokyo Institute of Technology, Oh-Okayama, Meguro, Tokyo 152-8552,January 2010.