domingo, 26 de abril de 2015

ecuaciones diferenciales



Valores propios m�ltiples: ejemplo

Ejemplo. Considere el sistema
displaymath87

1.
Encuentre la soluci�n general.
2.
Encuentre la soluci�n que satisface la condici�n inicialdisplaymath89
3.
Dibuje algunas soluciones en el fase-plano incluyendo la soluci�n encontrada en 2.
Respuesta. El coeficiente de la matriz del sistema es
displaymath91
Para encontrar los valores propios para considerar el polinomio caracter�stico
displaymath93
Puesto que  tex2html_wrap_inline95  , tenemos un valor propio repetido igual a 3. Encontremos el vector propio asociado  tex2html_wrap_inline97  . Fije
displaymath99
Entonces debemos tener  tex2html_wrap_inline101  a que traduce
displaymath103
Esto reduce a . Por lo tanto podemos tomar
displaymath107
Buscamos despu�s el segundo vector  tex2html_wrap_inline109  . La ecuaci�n que da este vector es  tex2html_wrap_inline111  cu�l traduce al sistema algebraico
displaymath113
donde
displaymath115
Las dos ecuaciones reducen claramente a la ecuaci�n = 1 o = 1 + , donde se puede elegir para ser cualquier n�mero. Tan si tomamos = 0 por ejemplo, conseguimos
displaymath125
Por lo tanto las dos soluciones independientes est�n
displaymath127
La soluci�n general entonces estar�
displaymath129

Para encontrar la soluci�n que satisface la condici�n inicial
displaymath89
debemos tener
displaymath133
Esto implica  tex2html_wrap_inline135  y  tex2html_wrap_inline137  . Por lo tanto la soluci�n est�
displaymath139

El plano de la fase con algunas soluciones se da en el cuadro abajo:











Valores propios Complejos


Considere el sistema homog�neo linear
displaymath221
El polinomio caracter�stico es
displaymath223
En esta secci�n, consideramos el caso cuando la ecuaci�n cuadr�tica antedicha tiene ra�ces complejas (que es si  tex2html_wrap_inline225  ). Las ra�ces (valores propios) son
displaymath227
donde
displaymath229
En este caso, la dificultad miente con la definici�n de
displaymath231
Para conseguir alrededor de esta dificultad utilizamos el f�rmula de Euler
displaymath233
Por lo tanto, tenemos
displaymath235
En este caso, el vector propio asociado  tex2html_wrap_inline237  a la voluntad tiene componentes complejos.

Ejemplo. Encuentre los valores propios y los vectores propios de la matriz
displaymath239
Respuesta. El polinomio caracter�stico es
displaymath241
Sus ra�ces son
displaymath243
Sistema  tex2html_wrap_inline245  . El vector propio asociado es dado por la ecuaci�n  tex2html_wrap_inline249  . Sistema
displaymath251
La ecuaci�n  tex2html_wrap_inline249  traduce a
displaymath255
Puesto que  tex2html_wrap_inline257  , entonces las dos ecuaciones son iguales (que deben haber esperado, usted ven porqu�?). Por lo tanto tenemos  tex2html_wrap_inline259  que implica que es un vector propio
displaymath261
Lo dejamos al lector para demostrar que para el valor propio  tex2html_wrap_inline263  , el vector propio est�
displaymath265

Vayamos de nuevo al sistema
displaymath221
con valores propios complejos  tex2html_wrap_inline269  . Observe eso si , donde
displaymath273
es un vector propio asociado a  tex2html_wrap_inline275  , entonces el vector
displaymath277
(donde  tex2html_wrap_inline279  est� la conjugaci�n de v) es un vector propio asociado a  tex2html_wrap_inline283  . Por otra parte, hemos visto eso
displaymath285
son las soluciones. Observe que estas soluciones son funciones complejas. Para encontrar soluciones verdaderas, utilizamos las observaciones antedichas. Sistema
displaymath287
entonces tenemos
displaymath289
cu�l da
displaymath291
Tenemos semejantemente
displaymath293
Poniendo todo juntos conseguimos
displaymath295
Esto implica claramente  tex2html_wrap_inline297  donde
displaymath299
Es f�cil ver que tenemos
displaymath301
Puesto que la suma y la diferencia de soluciones conducen a otra soluci�n, entonces ambas  tex2html_wrap_inline303  y  tex2html_wrap_inline305  son soluciones del sistema. �stas son soluciones verdaderas. Es muy f�cil llegar hecho de que son linear independientes. Resumamos la t�cnica antedicha.

Resumen (del caso complejo). Considere el sistema
displaymath221

 tex2html_wrap_inline309
Anote el polinomio caracter�sticodisplaymath223
y encuentre sus ra�ces
displaymath313
estamos asumiendo eso  tex2html_wrap_inline225  . Observe que en este paso, usted necesita saber  tex2html_wrap_inline317  y  tex2html_wrap_inline319  . El error com�n es olvidarse de dividirse por 2.
 tex2html_wrap_inline309
Encuentre un vector propio asociado al valor propio  tex2html_wrap_inline275  . Anote el vector propio comodisplaymath327

 tex2html_wrap_inline309
Dos soluciones independientes son dadas linear por los f�rmulasdisplaymath299

 tex2html_wrap_inline309
La soluci�n general esdisplaymath335
donde  tex2html_wrap_inline337  y  tex2html_wrap_inline339  est�n los n�meros arbitrarios. Observe que en este caso, tenemos
displaymath341
Ejemplo. Considere el oscilador arm�nico
displaymath343
Encuentre la soluci�n general usando la t�cnica del sistema.

Respuesta. Primero reescribimos la segunda ecuaci�n de la orden en el sistema
displaymath345
El coeficiente de la matriz de este sistema es
displaymath239
Hemos encontrado ya los valores propios y los vectores propios de esta matriz. Los valores propios est�n de hecho
displaymath349
Por lo tanto tenemos
displaymath351
El vector propio asociado a  tex2html_wrap_inline237  es
displaymath355
Anotamos despu�s las dos linear soluciones independientes
displaymath357
y
displaymath359
La soluci�n general del sistema equivalente es
displaymath335
o
displaymath423

Debajo de nosotros dibujamos algunas soluciones. El aviso c�mo las soluciones tuercen en espiral y ti�en en el origen (v�ase la discusi�n abajo)


Puesto que estamos buscando la soluci�n general de la ecuaci�n diferencial, consideramos solamente el primer componente. Por lo tanto tenemos
displaymath363
Usted puede desear comprobar que el segundo componente sea justo el derivado de y .
Debajo de nosotros dibujamos algunas soluciones para la ecuaci�n diferencial


An�lisis cualitativo de sistemas con valores propios complejos. 

Recuerde que en este caso, la soluci�n general est� dada cerca
displaymath341

El comportamiento de las soluciones en el plano de la fase depende de la parte real  tex2html_wrap_inline317  . De hecho, tenemos tres casos:


 tex2html_wrap_inline309
el caso:  tex2html_wrap_inline14  . Las soluciones tienden al origen (cuando  tex2html_wrap_inline437  ) mientras que tuercen en espiral. En este caso, el punto del equilibrio se llama un fregadero espiral .

 tex2html_wrap_inline309
El caso:  tex2html_wrap_inline441  Las soluciones estallan o consiguen lejos del origen (cuando  tex2html_wrap_inline437  ) mientras que tuercen en espiral. En este caso, el punto del equilibrio se llama una fuente espiral .

 tex2html_wrap_inline309
El caso:  tex2html_wrap_inline16  Las soluciones son peri�dicas. Esto significa que la trayectoria es curvas o ciclos cerrados. En este caso, el punto del equilibrio se llama un centro .
 

No hay comentarios:

Publicar un comentario