Exercises

Optimization of ellipsoid shapes

  1. For P in mathbf{R}^{n times n}, with P symmetric and positive-definite, we define the ellipsoid

 mathbf{E}(P) := left{ x ~:~ x^TP x le 1 right} .

A measure of the ‘‘size’’ of the ellipsoid is mbox{bf Tr} P^{-1}, with mbox{bf Tr} the trace (the sum of the diagonal elements of the matrix argument).

    1. Motivate our choice of the size function. Hint: Figure out the the semi-axis lengths of the ellipsoid as a function of P.

    2. Show that mathbf{E}(P) = { R^{-1}u ::: |u|_2 le 1}, where R is a factor of P (any matrix R such that P = R^TR; in matlab, R can be obtained via the command chol.)

    3. Show that for any Q=Q^T, Q positive-definite, the set {cal E}(Q) is contained in {cal E}(P) if and only if Q-P is positive semi-definite.

    4. For given n times n symmetric matrices P_i, i=1,2, both positive-definite, show how to compute an ellipsoid, centered at the origin, that contains both {cal E}(P_1),{cal E}(P_2) and has minimum size.

    5. Implement and plot your result with the data contained in here. (The function P rightarrow mbox{bf Tr} P^{-1} for P symmetric and positive definite, is implemented in CVX using trace_inv.)