domingo, 26 de abril de 2015

ecuaciones diferenciales



Ecuaciones de segundo orden y sistemas

Hemos visto antes c�mo solucionar un cierto tipo de ecuaciones de segundo orden. De hecho, hemos visto solamente c�mo solucionar ecuaciones lineales. �Qu� ocurre con las ecuaciones no lineales? Desafortunadamente muchos de problemas de la vida real son modelados por ecuaciones no lineales. Aqu� demostraremos c�mo se puede pasar una ecuaci�n de segundo orden a sistema. La t�cnica desarrollada para el sistema se puede entonces utilizar para estudiar las ecuaciones de segundo orden
Sea cualquier ecuaci�n diferencial de orden dos (en la forma expl�cita)
displaymath73
Introduzcamos la funci�n
displaymath75
Entonces tenemos
displaymath77
Poniendo todo juntos conseguimos
displaymath79
Es muy f�cil ver que (t) es soluci�n a la segunda ecuaci�n de la orden si y solamente si ( y , v) es soluci�n al sistema. El plano de fase del sistema no es pues importante como si nos ocup�ramos solamente del sistema asociado. �sta es otra raz�n por la que debemos tambi�n prestar la atenci�n al gr�fico de y contra .

Oscilador Arm�nico 
En f�sica, a una masa con un muelle se llama oscilador arm�nico y es uno de los modelos m�s importantes de la ciencia. Este modelo tambi�n aparece en la teor�a del circuito (circuitos de RLC) y en la f�sica de part�culas.
Una descripci�n general del sistema es
 tex2html_wrap_inline91
la masa del objeto unido al resorte;
 tex2html_wrap_inline91
la constante del resorte  tex2html_wrap_inline97  (que es un resultado directo de la ley de Hooke);
 tex2html_wrap_inline91
el coeficiente de "damping"  tex2html_wrap_inline101  que puede ser proporcional al vector de la velocidad o tener una forma muy complicada (no lineal).
De llas leyes de la mec�nica newtoniana, podemos deducir la ecuaci�n diferencial que modela el oscilador arm�nico
displaymath103
en ausencia de una fuerza externa que act�a en el objeto. Recuerde que (t) denota la posici�n del objeto en el tiempo . Estamos asumiendo claramente que el movimiento es lineal (a lo largo de una l�nea). Antes de que anotemos el sistema asociado, reescribimos la ecuaci�n en la forma expl�cita
displaymath109
El sistema asociado es
displaymath111

Osciladores Arm�nicos sin "damping" 
�stos son los osciladores arm�nicos para los cuales  tex2html_wrap_inline113  . En este caso, la ecuaci�n diferencial se reduce a
displaymath115
o equivalentemente
displaymath117
o
displaymath109
El sistema asociado es
displaymath121

Ejemplo. Considere los osciladores arm�nicos
1.
= 1,  tex2html_wrap_inline125  y  tex2html_wrap_inline127  ;
2.
= 1,  tex2html_wrap_inline131  y  tex2html_wrap_inline133  ;
3.
= 1,  tex2html_wrap_inline113  y  tex2html_wrap_inline139  .
Los planos de la fase con algunas soluciones son
1.
2.
3.










Introducci�n y motivaci�n

Las ecuaciones diferenciales son muy provechosas en muchas �reas de la ciencia. Pero la mayor�a de problemas verdaderos interesantes de la vida implican m�s de una funci�n desconocida. Por lo tanto, el uso del sistema de ecuaciones diferenciales es muy �til. Sin p�rdida de generalidad, nos concentraremos en los sistemas de dos ecuaciones diferenciales

displaymath80

Como motivaci�n vamos a considerar una isla con dos especies: Conejos y zorros. Uno desempe�a el papel del depredador mientras que el otro el papel de presa. Si estamos interesados modelar los crecimientos de las poblaciones de ambas especies, tenemos que tener presente que si, por ejemplo, la poblaci�n del zorro aumenta, entonces la poblaci�n del conejo se ver� afectada. El �ndice del cambio de la poblaci�n de un tipo depender� tan de la poblaci�n real del otro tipo. Por ejemplo, en ausencia de la poblaci�n del conejo, la poblaci�n del zorro disminuir�  para hacer frente a cierta extinci�n. Algo que la mayor�a de nosotros quisieran evitar. Un modelo para este problema de la Despredador-Presa fue desarrollado por Lotka (en 1925) y Volterra (en 1926) y se conoce como el sistema de Lotka-Volterra
displaymath82
donde (t) mide la poblaci�n de conejos, (t) mide la poblaci�n de zorros, y todas la constantes implicadas  tex2html_wrap_inline88  son n�meros positivos. Observe que son el �ndice de crecimiento de la presa, y el �ndice de mortalidad del depredador.  tex2html_wrap_inline94  y  tex2html_wrap_inline96  son las medidas entre efecto de la interacci�n de los conejos y los zorros.

Observe que en el sistema de Lotka-Volterra, falta la variable t. Esta clase de sistema se llama Sistemas aut�nomos y se escriben
displaymath100

Plano De la Fase 

Vayamos de nuevo al caso general
displaymath80
Una soluci�n a este sistema es el par de las funciones ((t), (t)) que satisfacen ambas ecuaciones diferenciales del sistema. Cuando variamos t, conseguimos un sistema de puntos en el plano XY que, en la f�sica, generalmente llamamos una trayectoria. El objeto m�vil tiene los coordenadas (x(t), (t)) en el tiempo . La velocidad a la trayectoria en el tiempo ser�
displaymath114
Observe que no necesitamos saber la soluci�n ((t), (t)) para determinar el vector de la velocidad en el tiempo . De hecho, tenemos
displaymath120
mientras sabemos , y . En detalle, podemos dibujar todos los vectores de la velocidad por todas partes en el plano para los sistemas aut�nomos. Esto se conoce como el campo de vectores.
El campo de vectores se debe entender como el an�logo del campo de la direcci�n para las ecuaciones diferenciales.
Ejemplo. Dibuje el campo de vectores para el problema presa-despredador
displaymath128
Por ejemplo, para el punto ( x , y) = (2.1), tenemos
displaymath132
Si ponemos m�s vectores de la velocidad, conseguimos
Est� claro que la longitud de los vectores de la velocidad puede afectar nuestra comprensi�n de las soluciones. Si somos interesados solamente por la direcci�n del movimiento no su velocidad, despu�s es natural fijar una longitud para estos vectores
Compare los dos cuadros.

Observaci�n. Observe que para el problema presa-depredador, podemos estar interesados descubrir si una de las dos especies se extingue. En decir podemos estar interesados estudiar las funciones (t) y y(t) por separado. Es decir dos m�s gr�ficos se asocian naturalmente a un sistema. Por el ejemplo antedicho, el gr�fico de la soluci�n que satisface la condici�n
displaymath138
se dibuja abajo
y los gr�ficos de x (t) y de y (t) son
y

Soluciones de Equilibrio. 
Una soluci�n de equilibrio (o soluci�n cr�tica ) de un Sistema Aut�nomo es la trayectoria de un objeto m�vil que no se mueve. O sea el vector de  velocidad es b�sicamente igual a 0. Considere el Sistema aut�nomo
displaymath100
Las soluciones del equilibrio son las soluciones algebraicas del sistema
displaymath146
Las soluciones del equilibrio tambi�n se llaman los puntos del equilibrio puesto que la trayectoria asociada es exactamente puntos en el plano de fase.
Ejemplo. Encuentre los puntos del equilibrio del sistema de Duffin
displaymath148
Respuesta. Consideramos el sistema
displaymath150
Tapando = 0 en la segunda ecuaci�n, conseguimos
displaymath154
lo cu�l da = 0, o  tex2html_wrap_inline158  . Por lo tanto los puntos del equilibrio son
displaymath160

No hay comentarios:

Publicar un comentario