El teorema de Lambert se deriva de la ecuación de Kepler . Cuando un cuerpo se mueve de un punto P1 a otro P2 de su trayectoria, emplea un tiempo t 2 -t 1 . Este intervalo de tiempo depende de dos parámetros:
la suma r 1 +r 2 de las longitudes de los radio vectores que unen el el centro de fuerzas situado en el foco con los puntos P1 y P2
la longitud de la cuerda C, o la distancia entre los puntos P1 y P2
Trayectorias elípticas
Diferencia de tiempos t 2 -t 1
Los instantes t 1 y t 2 se calculan a partir de la ecuación de Kepler para trayectorias elípticas
t
1
=
P
2 π
(
E
1
− ε sin
E
1
)
t
2
=
P
2 π
(
E
2
− ε sin
E
2
)
t
2
−
t
1
=
P
2 π
(
E
2
−
E
1
− ε (
sin
E
2
− sin
E
1
)
)
Donde E 1 es el ángulo que forma el vector posición OP1 con el eje X y E 2 es el ángulo que forma el vector OP2 con el eje X. ε <1, es la excentricidad de la elipse.
Utilizamos las relaciones trigonométrics
{
sin ( A + B ) = sin A cos B + sin B cos A
sin ( A − B ) = sin A cos B − sin B cos A
sin ( A + B ) − sin ( A − B ) = 2 sin B cos A
Llamando u=A+B y v=A-B
A =
u + v
2
, B =
u − v
2
sin u − sin v = 2 cos (
u + v
2
) sin (
u − v
2
)
El resultado es
t
2
−
t
1
=
P
2 π
{
E
2
−
E
1
− 2 ε cos (
E
2
+
E
1
2
) sin (
E
2
−
E
1
2
)
}
Distancia C entre los dos posiciones P1 y P2
Utilizamos la ecuación de una elipse en términos del ángulo E
r = a ( 1 − ε cos E )
Calculamos las coordenadas del punto P, x=r cosθ , y=r sinθ en términos del ángulo E
De la relación entre θ y E
tan (
θ
2
) =
1 + ε
1 − ε
tan (
E
2
)
Hemos obtenido
cos θ =
cos E − ε
1 − ε cos E
, sin θ =
1 −
cos
2
θ
=
1 −
ε
2
sin E
1 − ε cos E
Las coordenadas x e y del punto P son
{
x = a (
cos E − ε
)
y = a
1 −
ε
2
sin E
La longitud de la cuerda C es
C
2
=
(
x
2
−
x
1
)
2
+
(
y
2
−
y
1
)
2
C
2
=
a
2
(
cos
E
2
− cos
E
1
)
2
+
a
2
(
1 −
ε
2
)
(
sin
E
2
− sin
E
1
)
2
C
2
= 2
a
2
− 2
a
2
cos (
E
2
−
E
1
) −
a
2
ε
2
(
sin
E
2
− sin
E
1
)
2
C
2
= 4
a
2
sin
2
(
E
2
−
E
1
2
) − 4
a
2
ε
2
cos
2
(
E
2
+
E
1
2
)
sin
2
(
E
2
−
E
1
2
)
C
2
= 4
a
2
sin
2
(
E
2
−
E
1
2
) {
1 −
ε
2
cos
2
(
E
2
+
E
1
2
)
}
Suma de las longitudes de los radio vectores, r 1 +r 2
r
1
+
r
2
= a ( 1 − ε cos
E
1
) + a ( 1 − ε cos
E
2
) = 2 a (
1 −
1
2
ε (
cos
E
1
+ cos
E
2
)
)
Utilizamos las relaciones trigonométricas
{
cos ( A + B ) = cos A cos B − sin A sin B
cos ( A − B ) = cos A cos B + sin A sin B
cos ( A + B ) + cos ( A − B ) = 2 cos A cos B
Llamando u=A+B y v=A-B
A =
u + v
2
, B =
u − v
2
cos u + cos v = 2 cos (
u + v
2
) cos (
u − v
2
)
El resultado es
r
1
+
r
2
= 2 a {
1 − ε cos (
E
2
+
E
1
2
) cos (
E
2
−
E
1
2
)
}
Relaciones
Los resultados parciales obtenidos son
{
t
2
−
t
1
=
P
2 π
{
E
2
−
E
1
− 2 ε cos (
E
2
+
E
1
2
) sin (
E
2
−
E
1
2
)
}
C
2
= 4
a
2
{
1 −
ε
2
cos
2
(
E
2
+
E
1
2
)
}
sin
2
(
E
2
−
E
1
2
)
r
1
+
r
2
= 2 a {
1 − ε cos (
E
2
+
E
1
2
) cos (
E
2
−
E
1
2
)
}
Para simplificar las expresiones, llamamos
x = ε cos (
E
2
+
E
1
2
) , y =
E
2
−
E
1
2
{
t
2
−
t
1
=
P
2 π
(
2 y − 2 x sin y
)
C = 2 a
1 −
x
2
sin y
r
1
+
r
2
= 2 a (
1 − x cos y
)
Sumamos y restamos las dos últimas ecuaciones
{
r
1
+
r
2
+ C = 2 a (
1 − x cos y +
1 −
x
2
sin y
)
r
1
+
r
2
− C = 2 a (
1 − x cos y −
1 −
x
2
sin y
)
Hacemos un nuevo cambio de variable
x = cos (
α + β
2
) , y =
α − β
2
{
t
2
−
t
1
=
P
2 π
(
α − β − 2 cos (
α + δ
2
) sin (
α − β
2
)
)
r
1
+
r
2
+ C = 2 a (
1 − cos (
α + β
2
) cos (
α − β
2
) + sin (
α + β
2
) sin (
α − β
2
)
)
r
1
+
r
2
− C = 2 a (
1 − cos (
α + β
2
) cos (
α − β
2
) − sin (
α + β
2
) sin (
α − β
2
)
)
Utilizando la relaciones trigonométricas
{
cos ( A + B ) = cos A cos B − sin A sin B
cos ( A − B ) = cos A cos B + sin A sin B
cos ( A − B ) − cos ( A + B ) = 2 sin A sin B
Llamando u=A+B y v=A-B
A =
u + v
2
, B =
u − v
2
cos v − cos u = 2 sin (
u + v
2
) sin (
u − v
2
)
El resultado es
{
r
1
+
r
2
+ C = 2 a (
1 −
1
2
(
cos β + cos α
) +
1
2
(
cos β − cos α
)
)
r
1
+
r
2
− C = 2 a (
1 −
1
2
(
cos β + cos α
) −
1
2
(
cos β − cos α
)
)
{
r
1
+
r
2
+ C = 2 a (
1 − cos α
) = 4 a
sin
2
α
2
r
1
+
r
2
− C = 2 a (
1 − cos β
) = 4 a
sin
2
β
2
Se ha empleado la relación
cos ( 2 θ ) =
cos
2
θ −
sin
2
θ
Conocidos la longitud de los radio vectores, r 1 y r 2 y la longitud de la cuerda C , los ángulos α y β valen
{
sin
α
2
=
1
2
r
1
+
r
2
+ C
a
sin
β
2
=
1
2
r
1
+
r
2
− C
a
Para calcular el intervalo de tiempo en función de los ángulos α y β utilizamos la relación trigonométrica, sinu -sinv . El resultado final es
t
2
−
t
1
=
P
2 π
(
α − β − (
sin α − sin β
)
)
El intervalo de tiempo t 2 -t 1 depende de los ángulos α y β , es decir, de r 1 +r 2 y la longitud C
Ejemplo
La posición P1 es θ 1 =0
La posición P2 es θ 2 =π/2
La ecuación de la elipse es
r =
d
1 + ε cos θ
=
a (
1 −
ε
2
)
1 + ε cos θ
La longitud de los radio vectores es
{
θ
1
= 0,
r
1
= a (
1 − ε
)
θ
2
=
π
2
,
r
2
= a (
1 −
ε
2
)
La distancia C entre los dos puntos P1 y P2 es
C =
r
1
2
+
r
2
2
= a (
1 − ε
)
1 +
(
1 + ε
)
2
= a (
1 − ε
)
2 + 2 ε +
ε
2
De otro modo. La relación entre los ángulos θ y E es
tan (
E
2
) =
1 − ε
1 + ε
tan (
θ
2
)
{
θ
1
= 0,
E
1
= 0
θ
2
=
π
2
, tan (
E
2
2
) =
1 − ε
1 + ε
La longitud de la cuerda C vale
x = ε cos (
E
2
2
) , y =
E
2
2
C = 2 a
1 −
x
2
sin y
Teniendo en cuenta las relaciones trigonométricas
cos θ =
1
1 +
tan
2
θ
, sin θ =
tan θ
1 +
tan
2
θ
El resultado es
x = ε cos (
E
2
2
) = ε
1 + ε
2
sin (
E
2
2
) =
1 − ε
2
C = 2 a
1 −
ε
2
1 + ε
2
1 − ε
2
= a
(
2 −
ε
2
−
ε
3
) (
1 − ε
)
= a (
1 − ε
)
2 + 2 ε +
ε
2
Conocidos r 2 +r 1 y C , los ángulos α y β son
{
C = a (
1 − ε
)
2 + 2 ε +
ε
2
r
1
+
r
2
= a (
1 − ε
) (
2 + ε
)
{
sin
α
2
=
1
2
1 − ε
2 + ε +
2 + 2 ε +
ε
2
sin
β
2
=
1
2
1 − ε
2 + ε −
2 + 2 ε +
ε
2
El intervalo de tiempo es
t
2
−
t
1
=
P
2 π
(
α − β − (
sin α − sin β
)
)
Calculamos cosα y cosβ , empleando las relaciones trigonométricas
cos α =
cos
2
(
α
2
) −
sin
2
(
α
2
) = 1 − 2
sin
2
(
α
2
)
cos α = 1 −
1
2
(
1 − ε
) (
2 + ε +
2 + 2 ε +
ε
2
)
cos β = 1 −
1
2
(
1 − ε
) (
2 + ε −
2 + 2 ε +
ε
2
)
El producto cosα ·cosβ es
cos α cos β =
− 1 + ε +
ε
2
+
ε
3
2
Calculamos sinα y sinβ , empleando las relaciones trigonométricas
sin α = 2 sin (
α
2
) cos (
α
2
) = 2 sin (
α
2
)
1 −
sin
2
(
α
2
)
sin α =
1 − ε
2 + ε +
2 + 2 ε +
ε
2
1 −
1
4
(
1 − ε
) (
2 + ε +
2 + 2 ε +
ε
2
)
=
1 − ε
2 + ε +
2 + 2 ε +
ε
2
−
1
4
(
1 − ε
)
(
2 + ε +
2 + 2 ε +
ε
2
)
2
sin α =
1 −
ε
2
2
ε
2
+ ε + 1 + ε
2 + 2 ε +
ε
2
sin β =
1 −
ε
2
2
ε
2
+ ε + 1 − ε
2 + 2 ε +
ε
2
El producto sinα ·sinβ es
sin α sin β =
1 −
ε
2
2
(
ε
2
+ ε + 1
)
2
−
ε
2
(
2 + 2 ε +
ε
2
)
=
1 −
ε
2
2
(
1 + ε
) =
1 + ε −
ε
2
−
ε
3
2
El coseno de la diferencia de dos ángulos es
cos (
α − β
) = cos α cos β + sin α sin β
cos (
α − β
) = ε
α − β = arccos (
ε
)
Elevamos al cuadrado la diferencia, sinα -sinβ
(
sin α − sin β
)
2
=
sin
2
α +
sin
2
β − 2 sin α sin β =
(
1 −
ε
2
2
) (
ε
2
+ ε + 1 + ε
2 + 2 ε +
ε
2
) + (
1 −
ε
2
2
) (
ε
2
+ ε + 1 − ε
2 + 2 ε +
ε
2
) − 2
1 + ε −
ε
2
−
ε
3
2
=
(
1 −
ε
2
) (
ε
2
+ ε + 1
) − (
1 + ε −
ε
2
−
ε
3
) =
ε
2
(
1 −
ε
2
)
sin α − sin β = ε
1 −
ε
2
El resultado final es
t
2
−
t
1
=
P
2 π
(
arccos ε − ε
1 −
ε
2
)
Trayectorias parabólicas
La trayectoria parabólica es una trayectoria elíptica cuyo semieje mayor a tiende a infinito. Entonces los ángulos α y β son pequeños y se pueden aproximar
{
α
2
≈
1
2
r
1
+
r
2
+ C
a
β
2
≈
1
2
r
1
+
r
2
− C
a
Desarrollamos en serie la función sin(x ) tomando los dos primeros términos
>> syms x;
>> taylor(sin(x))
ans =x^5/120 - x^3/6 + x
sin x ≈ x −
x
3
6
La diferencia de tiempos es aproximadamente
t
2
−
t
1
=
P
2 π
(
α − β − (
sin α − sin β
)
) ≈
P
2 π
(
α − β − (
α −
α
3
6
− β +
β
3
6
)
)
t
2
−
t
1
≈
P
2 π
(
α
3
6
−
β
3
6
)
El periodo o el tiempo que tarda el cuerpo en describir una trayectoria elíptica de semieje a es
P
2
= 4
π
2
a
3
G M
El intervalo de tiempos es independiente del semieje a de la elipse, cuando a es grande
t
2
−
t
1
≈
a
3 / 2
6
G M
{
(
r
1
+
r
2
+ C
a
)
3 / 2
−
(
r
1
+
r
2
− C
a
)
3 / 2
}
t
2
−
t
1
≈
1
6
G M
{
(
r
1
+
r
2
+ C
)
3 / 2
−
(
r
1
+
r
2
− C
)
3 / 2
}
Trayectorias hiperbólicas
Diferencia de tiempos t 2 -t 1
Los instantes t 1 y t 2 se calculan a partir de la ecuación de Kepler para trayectorias hiperbólicas
t
1
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
(
ε sinh
F
1
−
F
1
)
t
2
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
(
ε sinh
F
2
−
F
2
)
t
2
−
t
1
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
(
ε (
sinh
F
2
− sinh
F
1
) − (
F
2
−
F
1
)
)
Donde F 1 es el ángulo que forma el vector posición OP1 con el eje X y F 2 es el ángulo que forma el vector OP2 con el eje X. ε >1, es la excentricidad de la hipérbola.
Utilizamos las relaciones trigonométrics
{
sinh ( A + B ) = sinh A cosh B + sinh B cosh A
sinh ( A − B ) = sinh A cosh B − sinh B cosh A
sinh ( A + B ) − sinh ( A − B ) = 2 sinh B cosh A
Llamando u=A+B y v=A-B
A =
u + v
2
, B =
u − v
2
sinh u − sinh v = 2 cosh (
u + v
2
) sinh (
u − v
2
)
El resultado es
t
2
−
t
1
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
{
2 ε cosh (
F
2
+
F
1
2
) sinh (
F
2
−
F
1
2
) − (
F
2
−
F
1
)
}
Distancia C entre los dos posiciones P1 y P2
Utilizamos la ecuación de una hipérbola en términos del ángulo E
r = a ( ε cosh F − 1 )
Calculamos las coordenadas del punto P, x=r cosθ , y=r sinθ en términos del ángulo F
De la relación entre θ y F
tan (
θ
2
) =
ε + 1
ε − 1
tan (
F
2
)
Hemos obtenido
cos θ =
ε − cosh F
ε cosh F − 1
, sin θ =
1 −
cos
2
θ
=
1 −
ε
2
sinh F
ε cosh F − 1
Las coordenadas x e y del punto P son
{
x = a (
ε − cosh F
)
y = a
1 −
ε
2
sinh F
Hemos empleado la relación
cosh
2
x −
sinh
2
x = 1
La longitud de la cuerda C es
C
2
=
(
x
2
−
x
1
)
2
+
(
y
2
−
y
1
)
2
C
2
=
a
2
(
cosh
F
2
− cosh
F
1
)
2
+
a
2
(
ε
2
− 1
)
(
sinh
F
2
− sinh
F
1
)
2
C
2
= 4
a
2
sinh
2
(
F
2
+
F
1
2
)
sinh
2
(
F
2
−
F
1
2
) + 4
a
2
(
ε
2
− 1
)
cosh
2
(
F
2
+
F
1
2
)
sinh
2
(
F
2
−
F
1
2
)
C
2
= 4
a
2
sinh
2
(
F
2
+
F
1
2
) {
sinh
2
(
F
2
−
F
1
2
) + (
ε
2
− 1
)
cosh
2
(
F
2
+
F
1
2
)
}
C
2
= 4
a
2
sinh
2
(
F
2
+
F
1
2
) {
cosh
2
(
F
2
+
F
1
2
) − 1
}
Hemos empleado la relación
cosh u − cosh v = 2 sinh (
u + v
2
) sinh (
u − v
2
)
cuya demostración es similar a sinhu -sinhv
Suma de las longitudes de los radio vectores, r 1 +r 2
r
1
+
r
2
= a ( ε cosh
F
1
− 1 ) + a ( ε cosh
F
2
− 1 ) = 2 a (
1
2
ε (
cosh
F
1
+ cosh
F
2
) − 1
)
Utilizando la relación
cosh u + cosh v = 2 cosh (
u + v
2
) cosh (
u − v
2
)
El resultado es
r
1
+
r
2
= 2 a {
ε cosh (
F
2
+
F
1
2
) cosh (
F
2
−
F
1
2
) − 1
}
Relaciones
Los resultados parciales obtenidos son
{
t
2
−
t
1
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
{
2 ε cosh (
F
2
+
F
1
2
) sinh (
F
2
−
F
1
2
) − (
F
2
−
F
1
)
}
C
2
= 4
a
2
sinh
2
(
F
2
+
F
1
2
) {
cosh
2
(
F
2
+
F
1
2
) − 1
}
r
1
+
r
2
= 2 a {
ε cosh (
F
2
+
F
1
2
) cosh (
F
2
−
F
1
2
) − 1
}
Para simplificar las expresiones, llamamos
x = ε cosh (
F
2
+
F
1
2
) , y =
F
2
−
F
1
2
{
t
2
−
t
1
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
(
2 x sinh y − 2 y
)
C = 2 a
x
2
− 1
sinh y
r
1
+
r
2
= 2 a (
x cosh y − 1
)
Sumamos y restamos las dos últimas ecuaciones
{
r
1
+
r
2
+ C = 2 a (
x cosh y − 1 +
x
2
− 1
sinh y
)
r
1
+
r
2
+ C = 2 a (
x cosh y − 1 −
x
2
− 1
sinh y
)
Hacemos un nuevo cambio de variable
x = cosh (
α + β
2
) , y =
α − β
2
{
t
2
−
t
1
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
(
2 cosh (
α + δ
2
) sinh (
α − β
2
) − (
α − β
)
)
r
1
+
r
2
+ C = 2 a (
cosh (
α + β
2
) cosh (
α − β
2
) + sinh (
α + β
2
) sinh (
α − β
2
) − 1
)
r
1
+
r
2
− C = 2 a (
cosh (
α + β
2
) cosh (
α − β
2
) − sinh (
α + β
2
) sinh (
α − β
2
) − 1
)
que se simplifican, teniendo en cuenta las relaciones coshu +coshv y coshu -coshv
{
r
1
+
r
2
+ C = 2 a (
1
2
(
cosh α + cosh β
) +
1
2
(
cosh α − cosh β
) − 1
)
r
1
+
r
2
− C = 2 a (
1
2
(
cosh α + cosh β
) −
1
2
(
cosh α − cosh β
) − 1
)
El resultado es
{
r
1
+
r
2
+ C = 2 a (
cosh α − 1
) = 4 a
sinh
2
(
α
2
)
r
1
+
r
2
− C = 2 a (
cosh β − 1
) = 4 a
sinh
2
(
β
2
)
Se ha utilizado la relación
cosh ( 2 x ) =
cosh
2
x +
sinh
2
x
Conocidos la longitud de los radio vectores, r 1 y r 2 y la longitud de la cuerda C , los ángulos α y β valen
{
sinh
α
2
=
1
2
r
1
+
r
2
+ C
a
sinh
β
2
=
1
2
r
1
+
r
2
− C
a
Para calcular el intervalo de tiempo en función de los ángulos α y β utilizamos la relación sinhu -sinhv . El resultado final es
t
2
−
t
1
=
L
3
(
ε
2
− 1
)
3 / 2
G
2
M
2
m
3
(
sinh α − sinh β − (
α − β
)
)
El intervalo de tiempo t 2 -t 1 depende de los ángulos α y β , es decir, de r 1 +r 2 y la longitud C
Referencias
Victor G. Szebehely, Hans Mark. Adventures in Celestial Mechanics . Wiley-VCH (2004), pp. 115-123