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

Sumamos y restamos las dos últimas ecuaciones

{ r 1 + r 2 +C=2a( 1xcosy+ 1 x 2 siny ) r 1 + r 2 C=2a( 1xcosy 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( 1cos( α+β 2 )cos( αβ 2 )+sin( α+β 2 )sin( αβ 2 ) ) r 1 + r 2 C=2a( 1cos( α+β 2 )cos( αβ 2 )sin( α+β 2 )sin( αβ 2 ) )

Utilizando la relaciones trigonométricas

{ cos(A+B)=cosAcosBsinAsinB cos(AB)=cosAcosB+sinAsinB cos(AB)cos(A+B)=2sinAsinB

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

A= u+v 2 ,B= uv 2 cosvcosu=2sin( u+v 2 )sin( uv 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( 1cosα )=4a sin 2 α 2 r 1 + r 2 C=2a( 1cosβ )=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 )=12 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

sinxx 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 ( 2xsinhy2y ) C=2a x 2 1 sinhy r 1 + r 2 =2a( xcoshy1 )

Sumamos y restamos las dos últimas ecuaciones

{ r 1 + r 2 +C=2a( xcoshy1+ x 2 1 sinhy ) r 1 + r 2 +C=2a( xcoshy1 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