El teorema de Lambert

Trayectorias elípticas

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 =2a{ 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π ( 2y−2xsiny ) C=2a 1− x 2 siny r 1 + r 2 =2a( 1−xcosy )

Sumamos y restamos las dos últimas ecuaciones

{ r 1 + r 2 +C=2a( 1−xcosy+ 1− x 2 siny ) r 1 + r 2 −C=2a( 1−xcosy− 1− x 2 siny )

Hacemos un nuevo cambio de variable

x=cos( α+β 2 ), y= α−β 2 { t 2 − t 1 = P 2π ( α−β−2cos( α+δ 2 )sin( α−β 2 ) ) r 1 + r 2 +C=2a( 1−cos( α+β 2 )cos( α−β 2 )+sin( α+β 2 )sin( α−β 2 ) ) r 1 + r 2 −C=2a( 1−cos( α+β 2 )cos( α−β 2 )−sin( α+β 2 )sin( α−β 2 ) )

Utilizando la relaciones trigonométricas

{ cos(A+B)=cosAcosB−sinAsinB cos(A−B)=cosAcosB+sinAsinB cos(A−B)−cos(A+B)=2sinAsinB

Llamando u=A+B y v=A-B

A= u+v 2 , B= u−v 2 cosv−cosu=2sin( u+v 2 )sin( u−v 2 )

El resultado es

{ r 1 + r 2 +C=2a( 1− 1 2 ( cosβ+cosα )+ 1 2 ( cosβ−cosα ) ) r 1 + r 2 −C=2a( 1− 1 2 ( cosβ+cosα )− 1 2 ( cosβ−cosα ) ) { r 1 + r 2 +C=2a( 1−cosα )=4a sin 2 α 2 r 1 + r 2 −C=2a( 1−cosβ )=4a sin 2 β 2

Se ha empleado la relación

cos(2θ)= cos 2 θ− sin 2 θ

Conocidos la longitud de los radio vectores, r1 y r2 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 t2-t1 depende de los ángulos α y β, es decir, de r1+r2 y la longitud C

Ejemplo

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=2a 1− x 2 siny

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=2a 1− ε 2 1+ε 2 1−ε 2 =a ( 2− ε 2 − ε 3 )( 1−ε ) =a( 1−ε ) 2+2ε+ ε 2

Conocidos r2+r1 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α=2sin( α 2 )cos( α 2 )=2sin( α 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 β−2sinα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

sinx≈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 GM

El intervalo de tiempos es independiente del semieje a de la elipse, cuando a es grande

t 2 − t 1 ≈ a 3/2 6 GM { ( r 1 + r 2 +C a ) 3/2 − ( r 1 + r 2 −C a ) 3/2 } t 2 − t 1 ≈ 1 6 GM { ( r 1 + r 2 +C ) 3/2 − ( r 1 + r 2 −C ) 3/2 }

Trayectorias hiperbólicas

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 =2a{ ε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 ( 2xsinhy−2y ) C=2a x 2 −1 sinhy r 1 + r 2 =2a( xcoshy−1 )

Sumamos y restamos las dos últimas ecuaciones

{ r 1 + r 2 +C=2a( xcoshy−1+ x 2 −1 sinhy ) r 1 + r 2 +C=2a( xcoshy−1− x 2 −1 sinhy )

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 ( 2cosh( α+δ 2 )sinh( α−β 2 )−( α−β ) ) r 1 + r 2 +C=2a( cosh( α+β 2 )cosh( α−β 2 )+sinh( α+β 2 )sinh( α−β 2 )−1 ) r 1 + r 2 −C=2a( 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=2a( 1 2 ( coshα+coshβ )+ 1 2 ( coshα−coshβ )−1 ) r 1 + r 2 −C=2a( 1 2 ( coshα+coshβ )− 1 2 ( coshα−coshβ )−1 )

El resultado es

{ r 1 + r 2 +C=2a( coshα−1 )=4a sinh 2 ( α 2 ) r 1 + r 2 −C=2a( coshβ−1 )=4a sinh 2 ( β 2 )

Se ha utilizado la relación

cosh(2x)= cosh 2 x+ sinh 2 x

Conocidos la longitud de los radio vectores, r1 y r2 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 t2-t1 depende de los ángulos α y β, es decir, de r1+r2 y la longitud C

Referencias

Victor G. Szebehely, Hans Mark. Adventures in Celestial Mechanics. Wiley-VCH (2004), pp. 115-123