26/01/2020, 22:10:34 pm *
Bienvenido(a), Visitante. Por favor, ingresa o regístrate.
¿Perdiste tu email de activación?

Ingresar con nombre de usuario, contraseña y duración de la sesión
Noticias: Homenaje a aladan
 
 
Páginas: [1]   Ir Abajo
  Imprimir  
Autor Tema: Un sistema de ecuaciones diferenciales que parece fácil... pero no lo es  (Leído 1088 veces)
0 Usuarios y 1 Visitante están viendo este tema.
manooooh
Pleno*
*****

Karma: +1/-0
Desconectado Desconectado

Sexo: Masculino
Argentina Argentina

Mensajes: 2.547


Ver Perfil
« : 18/11/2018, 05:32:25 am »

Hola!

La cuestión es esta.

Estaba realizando ejercicios de sistemas de ecuaciones diferenciales ordinarias de primer orden (¡un tipo de ecuaciones extremadamente reducido!) cuando topé con este. Observemos el ejercicio:

Hallar la familia de líneas de campo de \[\vec f(x,y,z)=(y,z,x).\] El enunciado es relativamente claro: se pide resolver el siguiente sistema: \[\begin{cases}x'=y\\y'=z\\z'=x.\end{cases}\]
La verdad es que uno piensa que la solución de este sistema se puede expresar como funciones elementales; y esto es cierto. Lo que no es cierto es que este sistema tenga una solución sencilla y corta de términos; tiene tantos que uno podría perderse entre las variables y funciones.

A mi me resulta asombroso. La dificultad no sobrepasa las cuentas que hay que hacer, pero me parece algo tan impresionante que de tres ecuaciones tenga una solución rara. ¿A ustedes también?

¿Existen otros casos parecidos, como por ejemplo la integral "inocente" [texx]\int{\sqrt{\tan x}\,\mathrm dx}[/texx] o el UTF pero menos conocidos, como este sistema?



Tengo entendido que esto sucede porque las componentes no están "alineadas"; si el campo hubiese sido [texx]\vec f(x,y,z)=(x,y,z)[/texx] el problema hubiese sido mucho más fácil (¡mucho!). ¿Tiene alguna otra explicación? ¿Siempre lo que está "entrelazado" es lo que cuesta más "desenredar"?



Moraleja. No confíes en las primeras impresiones...

:risa:

Saludos
En línea
Masacroso
Pleno*
*****

Karma: +2/-0
Desconectado Desconectado

España España

Mensajes: 1.735


Ver Perfil
« Respuesta #1 : 18/11/2018, 13:41:37 pm »

Si llamas [texx]v:=(x,y,z)[/texx] entonces el sistema se puede escribir como [texx]v'=Av[/texx] donde

[texx]\displaystyle A=\begin{bmatrix}0&1&0\\0&0&1\\1&0&0\end{bmatrix}[/texx]

Entonces las soluciones vienen dadas por [texx]u=e^{tA}v(0)[/texx]. Como [texx]A^3=I[/texx] entonces se puede saber la forma de los coeficientes de la matriz [texx]e^{tA}[/texx]: [texx][e^{tA}]_{i,j}=\sum_{k=0}^\infty\frac{t^{3k+c}}{(3k+c)!}[/texx] para algún [texx]c\in\{0,1,2\}[/texx], que depende de [texx](i,j)[/texx], es decir, de la posición del coeficiente en la matriz.
En línea
manooooh
Pleno*
*****

Karma: +1/-0
Desconectado Desconectado

Sexo: Masculino
Argentina Argentina

Mensajes: 2.547


Ver Perfil
« Respuesta #2 : 18/11/2018, 14:10:49 pm »

Hola Masacroso

Interesante forma de verlo.

¿Por qué las soluciones tienen un [texx]v(0)[/texx]?

Si uno escribe [texx]e^{tA}v(0)[/texx] a mi parecer no hay una sumatoria que de lugar al resto de términos de la solución. ¿Esa expresión tiene incorporada la idea de sumar varios términos?

Saludos y gracias
En línea
Masacroso
Pleno*
*****

Karma: +2/-0
Desconectado Desconectado

España España

Mensajes: 1.735


Ver Perfil
« Respuesta #3 : 18/11/2018, 14:52:24 pm »

Hola Masacroso

Interesante forma de verlo.

¿Por qué las soluciones tienen un [texx]v(0)[/texx]?

Si uno escribe [texx]e^{tA}v(0)[/texx] a mi parecer no hay una sumatoria que de lugar al resto de términos de la solución. ¿Esa expresión tiene incorporada la idea de sumar varios términos?

Saludos y gracias

Bueno, olvida el [texx]v(0)[/texx], representa en verdad el vector [texx](x(0),y(0),z(0))[/texx], que es el valor que se suele dar para resolver problemas de valor inicial. La solución general es de la forma [texx]u(t)=e^{tA}w[/texx] para algún vector [texx]w\in\Bbb R^3[/texx].

La forma de [texx]e^{tA}[/texx] (cortesía de Wolfram Mathematica) es

[texx] e^{t A}=\left(
\begin{array}{ccc}
 \frac{1}{3} \left(2 e^{-t/2}
   \cos \left(\frac{\sqrt{3}
   t}{2}\right)+e^t\right) &
   \frac{1}{3} e^{-t/2}
   \left(-\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^{3
   t/2}+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{e^t}{3}-\frac{1}{3}
   e^{-t/2} \left(\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) \\
 \frac{e^t}{3}-\frac{1}{3}
   e^{-t/2} \left(\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{1}{3} \left(2
   e^{-t/2} \cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^t\right) &
   \frac{1}{3} e^{-t/2}
   \left(-\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^{3
   t/2}+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) \\
 \frac{1}{3} e^{-t/2}
   \left(-\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^{3
   t/2}+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{e^t}{3}-\frac{1}{3}
   e^{-t/2} \left(\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{1}{3} \left(2
   e^{-t/2} \cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^t\right) \\
\end{array}
\right)[/texx]

Los coeficientes matriciales también pueden expresarse con las series que he dejado antes. Es decir

[texx]\displaystyle {\frac{1}{3} \left(e^t+2 e^{-t/2}\cos\left(\frac{\sqrt{3}t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k}}{(3k)!}\\
\frac{1}{3} e^{-t/2} \left(e^{3 t/2}+\sqrt{3} \sin\left(\frac{\sqrt{3} t}{2}\right)-\cos\left(\frac{\sqrt{3} t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k+1}}{(3k+1)!}\\
\frac{e^t}{3}-\frac{1}{3} e^{-t/2} \left(\sqrt{3} \sin\left(\frac{\sqrt{3} t}{2}\right)+\cos\left(\frac{\sqrt{3} t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k+2}}{(3k+2)!}
}[/texx]

quedando

[texx]\displaystyle e^{tA}=\left(
\begin{array}{ccc}
 \sum _{k=0}^{\infty } \frac{t^{3 k}}{(3 k)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+1}}{(3 k+1)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+2}}{(3 k+2)!} \\
 \sum _{k=0}^{\infty } \frac{t^{3 k+2}}{(3 k+2)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k}}{(3 k)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+1}}{(3 k+1)!} \\
 \sum _{k=0}^{\infty } \frac{t^{3 k+1}}{(3 k+1)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+2}}{(3 k+2)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k}}{(3 k)!} \\
\end{array}
\right)[/texx]
En línea
manooooh
Pleno*
*****

Karma: +1/-0
Desconectado Desconectado

Sexo: Masculino
Argentina Argentina

Mensajes: 2.547


Ver Perfil
« Respuesta #4 : 18/11/2018, 15:29:24 pm »

Hola

Bueno, olvida el [texx]v(0)[/texx], representa en verdad el vector [texx](x(0),y(0),z(0))[/texx], que es el valor que se suele dar para resolver problemas de valor inicial. La solución general es de la forma [texx]u(t)=e^{tA}w[/texx] para algún vector [texx]w\in\Bbb R^3[/texx].

La forma de [texx]e^{tA}[/texx] (cortesía de Wolfram Mathematica) es

[texx]\displaystyle e^{t A}=\left(
\begin{array}{ccc}
 \frac{1}{3} \left(2 e^{-t/2}
   \cos \left(\frac{\sqrt{3}
   t}{2}\right)+e^t\right) &
   \frac{1}{3} e^{-t/2}
   \left(-\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^{3
   t/2}+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{e^t}{3}-\frac{1}{3}
   e^{-t/2} \left(\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) \\
 \frac{e^t}{3}-\frac{1}{3}
   e^{-t/2} \left(\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{1}{3} \left(2
   e^{-t/2} \cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^t\right) &
   \frac{1}{3} e^{-t/2}
   \left(-\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^{3
   t/2}+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) \\
 \frac{1}{3} e^{-t/2}
   \left(-\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^{3
   t/2}+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{e^t}{3}-\frac{1}{3}
   e^{-t/2} \left(\cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+\sqrt{3} \sin
   \left(\frac{\sqrt{3}
   t}{2}\right)\right) &
   \frac{1}{3} \left(2
   e^{-t/2} \cos
   \left(\frac{\sqrt{3}
   t}{2}\right)+e^t\right) \\
\end{array}
\right)[/texx]

Los coeficientes matriciales también pueden expresarse con las series que he dejado antes. Es decir

[texx]\displaystyle {\frac{1}{3} \left(e^t+2 e^{-t/2}\cos\left(\frac{\sqrt{3}t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k}}{(3k)!}\\
\frac{1}{3} e^{-t/2} \left(e^{3 t/2}+\sqrt{3} \sin\left(\frac{\sqrt{3} t}{2}\right)-\cos\left(\frac{\sqrt{3} t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k+1}}{(3k+1)!}\\
\frac{e^t}{3}-\frac{1}{3} e^{-t/2} \left(\sqrt{3} \sin\left(\frac{\sqrt{3} t}{2}\right)+\cos\left(\frac{\sqrt{3} t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k+2}}{(3k+2)!}
}[/texx]

Vale, creo que voy entendiendo.

¿Cómo se dedujo que, por ejemplo, [texx]\frac{1}{3} \left(e^t+2 e^{-t/2}\cos\left(\frac{\sqrt{3}t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k}}{(3k)!}[/texx]?

Saludos
En línea
Masacroso
Pleno*
*****

Karma: +2/-0
Desconectado Desconectado

España España

Mensajes: 1.735


Ver Perfil
« Respuesta #5 : 18/11/2018, 15:43:57 pm »

Vale, creo que voy entendiendo.

¿Cómo se dedujo que, por ejemplo, [texx]\frac{1}{3} \left(e^t+2 e^{-t/2}\cos\left(\frac{\sqrt{3}t}{2}\right)\right)=\sum_{k=0}^\infty\frac{t^{3k}}{(3k)!}[/texx]?

Saludos

Lo dedujo Wolfram Mathematica. Sabía que esa expresión era equivalente a una de las series antes descrita pero no sabía cual así que usé el programa para hacer una serie de Maclaurin truncada para saber cuál de las tres era.
En línea
manooooh
Pleno*
*****

Karma: +1/-0
Desconectado Desconectado

Sexo: Masculino
Argentina Argentina

Mensajes: 2.547


Ver Perfil
« Respuesta #6 : 18/11/2018, 16:39:03 pm »

Hola

Lo dedujo Wolfram Mathematica. Sabía que esa expresión era equivalente a una de las series antes descrita pero no sabía cual así que usé el programa para hacer una serie de Maclaurin truncada para saber cuál de las tres era.

[texx]\displaystyle e^{tA}=\left(
\begin{array}{ccc}
 \sum _{k=0}^{\infty } \frac{t^{3 k}}{(3 k)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+1}}{(3 k+1)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+2}}{(3 k+2)!} \\
 \sum _{k=0}^{\infty } \frac{t^{3 k+2}}{(3 k+2)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k}}{(3 k)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+1}}{(3 k+1)!} \\
 \sum _{k=0}^{\infty } \frac{t^{3 k+1}}{(3 k+1)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k+2}}{(3 k+2)!} & \sum
   _{k=0}^{\infty } \frac{t^{3 k}}{(3 k)!} \\
\end{array}
\right)[/texx]

Ahh, a esto me refería. ¡Gracias!



¿Conocen otros ejemplos que a primera vista parezcan simples de resolver pero luego no?

Saludos
En línea
geómetracat
Moderador Global
Pleno*
*

Karma: +0/-0
Desconectado Desconectado

Sexo: Masculino
España España

Mensajes: 892



Ver Perfil
« Respuesta #7 : 19/11/2018, 07:05:31 am »

Me parece que os estáis complicando la vida.
La observación clave es que ese sistema es equivalente a la ecuación diferencial homogénea [texx]x'''=x[/texx]. Una vez uno se da cuenta de ésto el problema se vuelve muy fácil, pues la ecuación tiene polinomio característico [texx]\lambda^3 - 1 =0[/texx], por tanto sus raíces son las raíces cúbicas de la unidad y ya puedes escribir la solución para [texx]x[/texx] como combinación lineal de exponenciales. Para sacar [texx]y[/texx] y [texx]z[/texx] derivas y listo.
En línea

La ecuación más bonita de las matemáticas: [texx]d^2=0[/texx]
Masacroso
Pleno*
*****

Karma: +2/-0
Desconectado Desconectado

España España

Mensajes: 1.735


Ver Perfil
« Respuesta #8 : 19/11/2018, 11:55:30 am »

Me parece que os estáis complicando la vida.
La observación clave es que ese sistema es equivalente a la ecuación diferencial homogénea [texx]x'''=x[/texx]. Una vez uno se da cuenta de ésto el problema se vuelve muy fácil, pues la ecuación tiene polinomio característico [texx]\lambda^3 - 1 =0[/texx], por tanto sus raíces son las raíces cúbicas de la unidad y ya puedes escribir la solución para [texx]x[/texx] como combinación lineal de exponenciales. Para sacar [texx]y[/texx] y [texx]z[/texx] derivas y listo.

Ajá... Lo cierto es que tengo un conocimiento muy limitado de ecuaciones diferenciales, fundamentalmente teórico además, así que tiendo a escribir esas cosas que son al final poco prácticas y bastante aparatosas.

En algún momento dedicaré mi tiempo a ampliar mi conocimiento en esa área... aunque realmente las ecuaciones diferenciales me parecen bastante aburridas y tediosas pero lo cierto es que aparecen en todas partes así que hay que conocer mejor el tema.
En línea
geómetracat
Moderador Global
Pleno*
*

Karma: +0/-0
Desconectado Desconectado

Sexo: Masculino
España España

Mensajes: 892



Ver Perfil
« Respuesta #9 : 19/11/2018, 14:43:26 pm »

Ajá... Lo cierto es que tengo un conocimiento muy limitado de ecuaciones diferenciales, fundamentalmente teórico además, así que tiendo a escribir esas cosas que son al final poco prácticas y bastante aparatosas.

En algún momento dedicaré mi tiempo a ampliar mi conocimiento en esa área... aunque realmente las ecuaciones diferenciales me parecen bastante aburridas y tediosas pero lo cierto es que aparecen en todas partes así que hay que conocer mejor el tema.

Te entiendo, a mí resolver ecuaciones diferenciales siempre me pareció aburridote. Pasa como con las integrales: solamente hay métodos generales de resolución para unos tipos muy concretos de EDOs, y si la que tienes no se ajusta a esos tipos tienes que buscarte la vida como sea (lo más probable es que ni siquera tengan solución expresable en términos de funciones elementales).
De todas formas, lo que es la teoría que hay detrás es muy bonita, y vale la pena estudiarla (problemas de existencia y unicidad, teoría cualitativa, sistemas dinámicos, etc).
En línea

La ecuación más bonita de las matemáticas: [texx]d^2=0[/texx]
manooooh
Pleno*
*****

Karma: +1/-0
Desconectado Desconectado

Sexo: Masculino
Argentina Argentina

Mensajes: 2.547


Ver Perfil
« Respuesta #10 : 19/11/2018, 19:23:15 pm »

Hola

Me gustó el enfoque de Masacroso por la novedad que para mí representa(ba) :guiño:.

geómetracat, es claro que todo se reduce a resolver [texx]x'''-x'=0[/texx] e ir "subiendo" en cada ecuación inconclusa que dejemos, pero lo sorprendente (al menos para mí) es la complejidad de las cuentas; del resto no tengamos ninguna duda de que es mecánico.

Te entiendo, a mí resolver ecuaciones diferenciales siempre me pareció aburridote. Pasa como con las integrales: solamente hay métodos generales de resolución para unos tipos muy concretos de EDOs, y si la que tienes no se ajusta a esos tipos tienes que buscarte la vida como sea (lo más probable es que ni siquera tengan solución expresable en términos de funciones elementales).

Yo considero que resolver tan pocas EDOs o integrales no es aburrido sino frustrante, por la baja calidad de conocimientos que los estudiantes de grado aprendemos. Me encantaría tomar un curso dedicado a EDOs o integrales, pero lamentablemente aquí en Argentina o no hay, o dejaron de existir hace años (tampoco ir a competir internacionalmente, porque allí se valora la rapidez y no la capacidad de la persona en resolver un problema).

Gracias por sus aportes y opiniones :sonrisa:.

Saludos
En línea
Abdulai
Moderador Global
Pleno*
*

Karma: +0/-0
Conectado Conectado

Sexo: Masculino
Argentina Argentina

Mensajes: 2.280


Ver Perfil
« Respuesta #11 : 20/11/2018, 09:07:27 am »

....
¿Existen otros casos parecidos, como por ejemplo la integral "inocente" [texx]\int{\sqrt{\tan x}\,\mathrm dx}[/texx] o el UTF pero menos conocidos, como este sistema?
...

Algo que en su momento me dejó ojiplático fue la región de convergencia de la sucesión  [texx]z_{n+1}=z_n^2+c[/texx]  en el plano complejo.

Lo que menos imaginaba es que algo tan sencillo pudiera ser asi.



En línea
Páginas: [1]   Ir Arriba
  Imprimir  
 
Ir a:  

Impulsado por MySQL Impulsado por PHP Powered by SMF 1.1.4 | SMF © 2006, Simple Machines LLC XHTML 1.0 válido! CSS válido!