4
$\begingroup$

This sounds like a simple problem, but I can't get it done. Given the general equation $ax^2 + by^2 + cz^2 + dxy + exz + fyz + gx + hy + iz + j = 0,$ what are the requirements on the coefficients so that this equation represents an ellipsoid? All coefficients are real.

Some more context. To make the representation unique, we establish another constraint: $a + c = 1$. I'm not sure how important that might be. The question is in homework where we have to fit this equation to a set of points, using QR-decomposition. That isn't a problem, the small followup question "how can you see the coefficients determine an ellipsoid" is the problem.

I tried working on the equation to get something usable, from basic term regrouping to putting stuff in matrices and decomposing these, which didn't really help. One theoretic way to do it, I think, is to find the center of mass of the surface, translate to the origin, then find the principal axes, transform again, and show the radii are all larger than 0. Sounds like a lot of work for a simple question, and I'm kind of stuck on the "center of mass" part, I don't feel the idea here is to integrate the formula over $\mathbb{R}^3$. I've also been thinking about calculating curvature or something with the Frenet-Serret formulas, but my geometry courses have been too long ago to know what I want here.

Is there a simple test to determine this? Something along the lines of "the eigenvalues of matrix foo are all positive" or "the value of bar is real" or I don't know what...

Thanks a lot for any insights you can give me!

1 Answers 1

6

There are two steps involved.

Step 1: translate to centre of mass

The goal of translation to the centre of mass is to remove the linear terms in your coefficients. You want to turn your expression into the form

$$ AX^2 + BY^2 + CZ^2 + D XY + E YZ + F XZ + G = 0$$

where $X = (x - x_0)$, $Y = (y - y_0)$ and $Z = (z-z_0)$, where $x_0, y_0, z_0$ are the coordinates for the centre of mass. To do so you can solve forward by completing the square, or solve backwards by plugging in the expressions for $X,Y,Z$ and matching coefficients.

Step 2: check for positivity of eigenvalues

Relative to the centre of mass, you can write your equation as

$$ Q(v) = -G $$

where $Q$ is the quadratic form

$$ Q(v) = AX^2 + BY^2 + \cdots + F XZ, \qquad v = (X,Y,Z) $$

Then recall that the level surfaces of a quadratic form define ellipsoids if and only if it is positive definite. So computing the eigenvalues (among the many ways to check that a quadratic form is positive definite) will get you the answer. (Of course, this will also require $G$ to be a negative number for there to be a solution.)

  • 0
    The matrix involved may be merely 3-by-3, but [Cholesky](http://math.stackexchange.com/questions/13282) remains a quicker way of verifying positive definiteness than explicitly computing eigenvalues. However, if you'll be doing more than merely verifying positive definiteness (e.g. you need the directions of the principal axes), then eigendecomposition is required.2010-12-13
  • 0
    Well, there are ways of verifying positive definiteness, but I'm almost sure the OP is supposed to be doing this by hand; Cholesky is surely an over-kill. He doesn't even need to compute the eigenvalues. All he need is the determinant and $AB - D^2$ both positive (since if $A,B,C$ not all positive you can trivially rule-out positive-definiteness...)2010-12-13
  • 0
    "The question is in homework where we have to fit this equation to a set of points, using QR-decomposition." - I supposed from this that this is a programming problem... unless he'd be doing least-squares via Gram-Schmidt, but I don't believe anybody's that cruel...2010-12-13
  • 0
    @J.M. Why not? I've given homework questions like that before. And the phrasing of the follow-up question suggests to me it is not programming. But as long as the OP does not specify, this debate is purely academic.2010-12-13
  • 0
    "purely academic" - That's a fair point, Willie. :) Anyway, I just gave a useful approach if somebody's going to be doing this by computer instead of doing it manually.2010-12-13
  • 0
    For a 3 by 3 matrix, Sylvester's criterion for checking for positive definiteness is pretty quick and painless.2010-12-13
  • 1
    Yes, this works out. Thanks a lot! I did the calculations in Matlab, but I still had to provide proof of my method in the report, obviously. The QR-decomposition was also with computer, I wouldn't want to do the least-squares fit of 150 points with Gram-Schmidt...2010-12-13
  • 0
    150 points? Oh... even I can't be that evil. :-) @Zaricuse: I knew there is a name for what I wrote about $AB-D^2$, thanks.2010-12-13