MateWiki - Contribuciones del usuario [es]

Archivo:0-2000 e).jpg

2014-03-04T22:39:05Z

Noemí Ortiz:

Archivo:0-500 e).jpg

2014-03-04T22:38:31Z

Noemí Ortiz:

Archivo:0-100 e).jpg

2014-03-04T22:37:56Z

Noemí Ortiz:

Archivo:0-2000 d).jpg

2014-03-04T22:37:30Z

Noemí Ortiz:

Archivo:0-500 d).jpg

2014-03-04T22:37:01Z

Noemí Ortiz:

Archivo:0-100 d).jpg

2014-03-04T22:36:25Z

Noemí Ortiz:

Archivo:0-2000 c).jpg

2014-03-04T22:35:58Z

Noemí Ortiz:

Archivo:0-500 c).jpg

2014-03-04T22:35:30Z

Noemí Ortiz:

Archivo:0-100 c).jpg

2014-03-04T22:34:54Z

Noemí Ortiz:

Archivo:0-2000 b).jpg

2014-03-04T22:34:26Z

Noemí Ortiz:

Archivo:0-500 b).jpg

2014-03-04T22:33:55Z

Noemí Ortiz:

Archivo:0-100 b).jpg

2014-03-04T22:33:27Z

Noemí Ortiz:

Archivo:0-2000 a).jpg

2014-03-04T22:32:47Z

Noemí Ortiz:

Archivo:0-500 a).jpg

2014-03-04T22:32:10Z

Noemí Ortiz:

Orbita plana de una luna entorno a un planeta

2014-03-04T22:31:38Z

Noemí Ortiz:

{{ TrabajoED | Orbita plana de una luna entorno a un planeta. Grupo 13-B | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED13/14|Curso 2013-14]] | Noemí Ortiz, Alicia Chacón, Miguel Sánchez, Mónica Gómez, Cristina Jiménez}}

El trabajo realizado consiste en el estudio de la órbita lunar alrededor de un planeta.
Para ello disponemos de dos ecuaciones diferenciales de segundo grado, que describen dicha trayectoria, y que deberemos convertir en un sistema de cuatro ecuaciones diferenciales no lineales de primer orden, para poder aplicar los métodos de cálculo numérico correspondientes.

\[\left\{\begin{matrix}\ x''=-G \frac{mx}{(x^{2}+y^{2})^{3/2}}\ , & t\in [0,T]\\ y''=-G \frac{my}{(x^{2}+y^{2})^{3/2}}\ , & \end{matrix}\right.\]

==Reducción al sistema==
Al disponer de una parametrización de dos variables, x(t) e y(t), para transformarlo en el sistema anteriormente citado, hemos tomado dichas variables y sus derivadas, asignando a cada una, una nueva variable.
:<math> x(t)=y_1(t) </math>
:<math> x'(t)=y_2(t) </math>
:<math> y(t)=y_3(t) </math>
:<math> y'(t)=y_4(t) </math>
De manera que nuestras ecuaciones iniciales quedarían definidas de la siguiente forma:
:<math> y'_2= x''= -G \frac{mx}{(x^{2}+y^{2})^{3/2}}</math>
:<math> y'_4= y''= -G \frac{my}{(x^{2}+y^{2})^{3/2}}</math>
Y el sistema final con el que trabajaríamos sería el siguiente:
\[\left\{\begin{matrix}\ y'_1(t)=y_2(t) \\ y'_2(t)= -G \frac{my_1}{(y_1^{2}+y_2^{2})^{3/2}}\\ y'_3(t)= y_4(t)\\ y'_4(t)= -G \frac{my_3}{(y_1^{2}+y_3^{2})^{3/2}} & \end{matrix}\right.\]

A continuación explicamos los dos métodos de resolución empleados para este sistema.

== Método de Euler modificado ==
Se llama método de Euler al método numérico consistente en ir incrementando paso a paso la variable independiente y hallando su imagen con la derivada correspondiente.
Por otra parte, el método de Euler modificado se basa en la misma idea pero hace un refinamiento en la aproximación, tomando un promedio entre ciertas pendientes.
Este método es equivalente a un Runge-Kutta de orden 2.

\[y_{n+1}= y_{n}+h\cdot \frac{[(f(x_{n},y_{n}+f(x_{n+1},y_{n+1})]}{2}\]
Donde
\[y_{n+1}= y_{n}+h\cdot f(x_{n},y_{n})\]
Dicho método en código Matlab es el siguiente:
{{matlab|codigo=
%Médoto Euler modificado
clear all
t0=0;
tf=100;
%longitud de paso
h=0.1;
N=(tf-t0)/h;
t=t0:h:tf;
% definimos a b c d
a=1; b=0; c=0; d=1;
y0=[a ;b;c;d];
y(:,1)=y0;
yy=y0;
Gm=1;
%Bucle

for n=1:N
r=(yy(1)^2+yy(3)^2)^(3/2);
K1=[yy(2); (-Gm*yy(1)/r);yy(4); (-Gm*yy(3)/r)];
r1=[((yy(1)+h*K1(1))^2+(yy(3)+h*K1(3))^2)^(3/2)];
K2=[yy(2)+h*K1(2); -(Gm*(yy(1)+h*K1(1)))/r1; yy(4)+h*K1(4);-(Gm*(yy(3)+h*K1(3)))/r1];
yy=yy+(h/2)*(K1+K2);
y(:,n+1)=yy;
end

plot(y(1,:),y(3,:))}}

Además hemos ido variando la longitud de paso para demostrar que la precisión del método es menor cuanto mayor es la longitud de paso.
A continuación se muestran las gráficas correspondientes a dichas variaciones.

[[Archivo:Euler.jpeg|thumb|left|400px|Gráfica Euler para h=0.1]]
[[Archivo:Euler10.jpg|thumb|left|400px|Gráfica Euler para h=1]]

== Método Runge-Kutta ==
Los métodos de Runge-Kutta son un conjunto de métodos iterativos (implícitos y explícitos) para la aproximación de soluciones de ecuaciones diferenciales ordinarias.
En nuestro caso utilizaremos el Runge-Kutta de orden 4 dado por la siguiente ecuación:
:<math>y_{i+1} = y_i + {1 \over 6}h\left ( k_1 + 2k_2 + 2k_3 + k_4 \right ) </math>
Donde

:<math>
\begin{cases}
k_1 & = f \left( x_i, y_i \right) \\
k_2 & = f \left( x_i + {1 \over 2}h, y_i + {1 \over 2}k_1 h \right) \\
k_3 & = f \left( x_i + {1 \over 2}h, y_i + {1 \over 2}k_2 h \right) \\
k_4 & = f \left( x_i + h, y_i + k_3h \right) \\
\end{cases}
</math>

A continuación mostramos el código Matlab utilizado para la descripción de la trayectoria utilizando este método:

{{matlab|codigo=
%metodo Runge-Kutta de cuarto orden

clear all
t0=0;
tf=2000;
%longitud de paso
h=0.1;
N=(tf-t0)/h;
t=t0:h:tf;
%Definimos a, b, c y d
a=1; b=0; c=0; d=1;
yo=[a;b;c;d]
y(:,1)=yo;
yy=yo;
Gm=1.1;

%Empezamos el bucle
for n=1:N
r=(yy(1)^2+yy(3)^2)^(3/2);
K1=[yy(2); (-Gm*yy(1)/r);yy(4); (-Gm*yy(3)/r)];
r1=[((yy(1)+h*K1(1))^2+(yy(3)+h*K1(3))^2)^(3/2)];
K2=[yy(2)+(h/2)*K1(2); -(Gm*(yy(1)+(h/2)*K1(1)))/r1; yy(4)+(h/2)*K1(4);-(Gm*(yy(3)+(h/2)*K1(3)))/r1];
K3=[yy(2)+(h/2)*K2(2); -(Gm*(yy(1)+(h/2)*K2(1)))/r1; yy(4)+(h/2)*K2(4);-(Gm*(yy(3)+(h/2)*K2(3)))/r1];
K4=[yy(2)+(h)*K3(2); -(Gm*(yy(1)+(h)*K3(1)))/r1; yy(4)+(h)*K3(4);-(Gm*(yy(3)+(h)*K3(3)))/r1];
yy=yy+(h/6)*(K1+2*(K2)+2*K3+K4);
y(:,n+1)=yy;
end

plot(y(1,:),y(3,:))}}

Y al igual que en el método anterior, hemos obtenido diferentes gráficas en función de los diferentes valores que hemos asignado a la longitud de paso y a los extremos del intervalo:

[[Archivo:0-100_a).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,100] para a=1, b=c=0, d=1, Gm=1]]
[[Archivo:0-500_a).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,500] para a=1, b=c=0, d=1, Gm=1]]
[[Archivo:0-2000_a).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,2000] para a=1, b=c=0, d=1, Gm=1]]
[[Archivo:0-100_b).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,100] para a=1.3, b=c=0, d=1, Gm=1]]
[[Archivo:0-500_b).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,500] para a=1.3, b=c=0, d=1, Gm=1]]
[[Archivo:0-2000_b).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,2000] para a=1.3, b=c=0, d=1, Gm=1]]
[[Archivo:0-100_c).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,100] para a=1.3, b=c=0, d=1.1, Gm=1]]
[[Archivo:0-500_c).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,500] para a=1.3, b=c=0, d=1.1, Gm=1]]
[[Archivo:0-2000_c).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,2000] para a=1.3, b=c=0, d=1.1, Gm=1]]
[[Archivo:0-100_d).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,100] para a=1, b=c=0, d=1.4, Gm=1]]
[[Archivo:0-500_d).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,500] para a=1, b=c=0, d=1.4, Gm=1]]
[[Archivo:0-2000_d).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,2000] para a=1, b=c=0, d=1.4, Gm=1]]
[[Archivo:0-100_e).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,100] para a=1, b=c=0, d=1, Gm=1.1]]
[[Archivo:0-500_e).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,500] para a=1, b=c=0, d=1, Gm=1.1]]
[[Archivo:0-2000_e).jpg|thumb|left|400px|Gráfica Runge-kutta para h=1 en el intervalo [0,2000] para a=1, b=c=0, d=1, Gm=1.1]]

Archivo:0-100 a).jpg

2014-03-04T22:23:53Z

Noemí Ortiz:

Orbita plana de una luna entorno a un planeta

2014-03-04T22:23:04Z

Noemí Ortiz: