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
ì
ï
í
ï
î
k1 a1 x1
+
k1 b1 x2
=
k1 c1
k2 a2 x1
+
k2 b2 x2
=
k2 c2
k3 a3 x1
+
k3 b3 x2
=
k3 c3
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 .

kuva1.png
Figure 1: The triangle

We define
x = æ
ç
è
x1
x2
ö
÷
ø
,   c = æ
ç
ç
ç
è
k1 c1
k2 c2
k3 c3
ö
÷
÷
÷
ø
,   A = æ
ç
ç
ç
è
k1 a1
k1 b1
k2 a2
k2 b2
k3 a3
k3 b3
ö
÷
÷
÷
ø
The matrix equation Ax = C has in general no solution. The least squares solution x0 of the normal equation
AT A x = AT c
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?

kuva2.png

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 ?

kuva3.png
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.