John ellipsoid

Outer Löwner–John ellipsoid containing a set of a points in ⁠$\mathbb {R} ^{2}$⁠

In mathematics, the John ellipsoid or Löwner–John ellipsoid E(K) associated to a convex body K in n-dimensional Euclidean space ⁠$\mathbb {R} ^{n}$⁠ can refer to the n-dimensional ellipsoid of maximal volume contained within K or the ellipsoid of minimal volume that contains K.

Often, the minimal volume ellipsoid is called the Löwner ellipsoid, and the maximal volume ellipsoid is called the John ellipsoid (although John worked with the minimal volume ellipsoid in his original paper). One can also refer to the minimal volume circumscribed ellipsoid as the outer Löwner–John ellipsoid, and the maximum volume inscribed ellipsoid as the inner Löwner–John ellipsoid.

The German-American mathematician Fritz John proved in 1948 that each convex body in ⁠$\mathbb {R} ^{n}$⁠ is circumscribed by a unique ellipsoid of minimal volume, and that the dilation of this ellipsoid by factor 1/n is contained inside the convex body. That is, the outer Lowner-John ellipsoid is larger than the inner one by a factor of at most n. For a balanced body, this factor can be reduced to ${\sqrt {n}}.$

Properties

The inner Löwner–John ellipsoid E(K) of a convex body $K\subset \mathbb {R} ^{n}$ is a closed unit ball B in ⁠$\mathbb {R} ^{n}$⁠ if and only if B ⊆ K and there exists an integer m ≥ n and, for i = 1, ..., m, real numbers c_i > 0 and unit vectors $u_{i}\in S^{n-1}\cap \partial K$ (where S is the unit n-sphere) such that

$$\sum _{i=1}^{m}c_{i}u_{i}=0$$

and, for all $x\in \mathbb {R} ^{n}:$

$$x=\sum _{i=1}^{m}c_{i}(x\cdot u_{i})u_{i}.$$

Computation

In general, computing the John ellipsoid of a given convex body is a hard problem. However, for some specific cases, explicit formulas are known. Some cases are particularly important for the ellipsoid method.

Let E(A, a) be an ellipsoid in ⁠$\mathbb {R} ^{n},$⁠ defined by a matrix A and center a. Let c be a nonzero vector in ⁠$\mathbb {R} ^{n}.$⁠ Let E'(A, a, c) be the half-ellipsoid derived by cutting E(A, a) at its center using the hyperplane defined by c. Then, the Lowner-John ellipsoid of E'(A, a, c) is an ellipsoid E(A', a') defined by: $${\begin{aligned}\mathbf {a'} &=\mathbf {a} -{\frac {1}{n+1}}\mathbf {b} \\\mathbf {A'} &={\frac {n^{2}}{n^{2}-1}}\left(\mathbf {A} -{\frac {2}{n+1}}\mathbf {bb} ^{\mathrm {T} }\right)\end{aligned}}$$ where b is a vector defined by: $$\mathbf {b} ={\frac {1}{\sqrt {\mathbf {c} ^{\mathrm {T} }\mathbf {Ac} }}}\mathbf {Ac}$$ Similarly, there are formulas for other sections of ellipsoids, not necessarily through its center.

Applications

The computation of Löwner–John ellipsoids (and in more general, the computation of minimal-volume polynomial level sets enclosing a set) has found many applications in control and robotics. In particular, computing Löwner–John ellipsoids has applications in obstacle collision detection for robotic systems, where the distance between a robot and its surrounding environment is estimated using a best ellipsoid fit.

Löwner–John ellipsoids has also been used to approximate the optimal policy in portfolio optimization problems with transaction costs.

See also

• Banach–Mazur compactum – Concept in functional analysis
• Ellipsoid method
• Steiner inellipse, the special case of the inner Löwner–John ellipsoid for a triangle.
• Fat object, related to radius of largest contained ball.

References

• Gardner, Richard J. (2002). "The Brunn-Minkowski inequality". Bull. Amer. Math. Soc. (N.S.). 39 (3): 355–405 (electronic). doi:10.1090/S0273-0979-02-00941-2. ISSN 0273-0979.