Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Distance Between Lines in R^3 (fwd)
Replies: 15   Last Post: Sep 13, 2013 1:25 PM

 Messages: [ Previous | Next ]
 Thomas Nordhaus Posts: 433 Registered: 12/13/04
Distance between lines in R^n (Was: Distance Between Lines in R^3
(fwd))

Posted: Jul 27, 2013 6:40 AM

Am 26.07.2013 16:26, schrieb Bart Goddard:
> William Elliot <marsh@panix.com> wrote in
> news:Pine.NEB.4.64.1307252033010.24715@panix3.panix.com:
>

>> How do we find the shortest distance between two lines L, L' in R^3
>
> Each line has a direction vector. The cross product n of the two
> direction vectors is perpendicular to both lines. Pick a point
> on each line, P and Q. The projection of the vector from P to Q
> onto the vector n is the distance.

That's nice and simple. However the situation in R^3 seems to be special:

Let L: s -> v0 + s*v, L': t -> w0 + t*w. Suppose v and w aren't parallel
then there are exactly two unit vectors n1 and n2 that are perpendicular
to v and w. Furthermore n1 = -n2. You can compute them easily via the
vector (cross) product in R3: n1,n2 = +- (vxw)/norm(vxw).
Then the distance D(L,L') between the lines L and L' is given by

(A): D(L,L') = |<vxw,v0-w0>/(norm(vxw))|, <,> the scalar product.

Again, what makes this special is that in R^3 you only have one line
through the origin which is perpendicular to the plane spanned by v and
w. In R^n this is a hyperplane of dimension n-2.

On the other hand, using variational principles it is easy to compute
D(L,L') for all n>=2 just using the euclidean structure, i.e. scalar
products. I got the following general result (using Matlab-Notation for
matrices and vectors):

Let L: s -> v0 + s*v, L': t -> w0 + t*w be non-parallel lines.
Let M=[<w,w>,<v,w> ; <v,w>,<v,v>]/(<v,v><w,w> - <v,w>^2),
[S;T] = M*[<v,w0-v0>;<w,v0-w0>].

Then D(L,L') = norm(v0+S*v-(w0+T*w)).

(I derived this in the accompanying thread). This even works for n=2,
giving D(L,L')=0.

************
Question: Is there a generalization of the vector-product solution (A)
to higher dimensional euclidean space R^n, n>3 using some sort of
generalized vector-product?
************

--
Thomas Nordhaus

Date Subject Author
7/25/13 William Elliot
7/26/13 quasi
7/26/13 Thomas Nordhaus
7/26/13 Thomas Nordhaus
7/28/13 Ken.Pledger@vuw.ac.nz
7/28/13 Thomas Nordhaus
7/28/13 Virgil
7/28/13 Ken.Pledger@vuw.ac.nz
7/28/13 fom
7/29/13 Thomas Nordhaus
7/29/13 Peter Percival
7/26/13 Bart Goddard
7/27/13 Thomas Nordhaus
7/29/13 Tucsondrew@me.com
9/11/13 Brian Q. Hutchings
9/13/13 Brian Q. Hutchings