Hi!
I'm trying to solve a problem and am currently stuck, so I came here in the hope that someone could help me.
I have a 10x3 matrix (call it "A") containing position data for points in space and I need to find out if these lie on a circle using QR-factorization (and, later on, the least square method). It's a matlab excercise, so the calculations themselves aren't a concern, but rather the theory behind them.
I've been staring at the problem for a long time now, but I can't quite figure it out.
I realize that the points must lie on a plane, i.e. a subspace R2 to the space R3 in which the data points are scattered.
The equation of a circle is (x - a)^2 + (y - b)^2 = r. Or, if the circle is centered around origo, x^2 + y^2 = r (but the given data points don't appear to be).
I know the principle of QR-factorization, and have managed to QR-factorize A. When QR-factoring A into QR, Q forms an orthonormal basis for col(A), which I think is the clue to determining if the points form a circle. But that's about as far as I get.
I can't quite see the full connection to QR-factorization and I'm not sure how to write the equation of the circle so that the least square approximation is applicable.
Can anyone please give me some help?


LinkBack URL
About LinkBacks