1 NORMALIZING THE LEAST-SQUARES SOLUTION OF LINES
1.1 Description of the problem
Let
a1, b1, c1, k1,
a2, b2, c2, k2,
a3, b3, c3, k3,
be real numbers such that
a12 + b12 = a22 + b22 = a32 + b32 = 1. |
|
We call d1, d2 and d3 the straight lines of equations
Note that in the Java Sketches the scaling numbers k1, k2 and k3 are denoted by a, b and c.
The lines form a triangle of vertices S1 = d2Çd3, S2 = d3Çd1 and S3 = d1Çd2,
see Figure .
We define
x = |
æ ç
è
|
|
ö ÷
ø
|
, c = |
æ ç ç
ç è
|
|
ö ÷ ÷
÷ ø
|
, A = |
æ ç ç
ç è
|
|
ö ÷ ÷
÷ ø
|
|
|
The matrix equation Ax = C has in general no solution. The least squares solution
x0
of the normal equation
is such that
F(x) = F(x1,x2) = ||A x - c||2 = (xT AT - cT)(Ax - c) |
|
is minimum for x = x0.
We denote by d(x,L) the distance from the point x to line L.
If one computes F(x) one gets:
F(x) = k12 (d(x,d1))2+k22 (d(x,d2))2+k32 (d(x,d3))2 |
|
I Case k1 = k2 = k3
The function to be minimized is
F(x) = (d(x,d1))2+(d(x,d2))2+(d(x,d3))2 |
|
Let us call l1 = S2 S3, l2 = S3 S1 and l3 = S1 S2,
the lengths of the sides of the triangle defined by the lines d1, d2 and d3.
The point x0 such that F(x) is minimum for x = x0
is the barycenter of S1, S2 and S3 with masses l12, l22 and l32:
x0 = |
l12
l12+l22+l32
|
S1+ |
l22
l12+l22+l32
|
S2+ |
l32
l12+l22+l32
|
S3 |
|
Problems
1.
If the triangle defined by the lines d1, d2 and d3 is isosceles (let us say S1 S2 = S1 S3),
on what line do you expect to find x0?
2.
If the triangle is equilateral where is x0?
3.
The distances of x0 to the sides of the triangle are proportional to the lengths of the
corresponding sides to a certain power r:
|
d(x0,d1)
l1r
|
= |
d(x0,d2)
l2r
|
= |
d(x0,d3)
l3r
|
|
|
Guess what is the value of r?
Note: If r = -1, the point x0 would be the center of gravity;
if r = 0, the point x0 would be the center of the inner circle.
4.
If the triangle is rectangular, e.g. d2 ^d3, on what line do you find x0
(see Figure ?
|
|
Figure 2: Rectangular triangle
|
|
H K2
S1 S3
|
= |
H K3
S2 S1
|
since the triangles are similar |
|
II Case k1, k2 k3 any real numbers
For this general case we already made experiments!
File translated from
TEX
by
TTH,
version 3.01.
On 16 Nov 2003, 15:48.