Este trabajo contrasta métodos numéricos con métodos algebraicos aplicados ambos a la resolución de un sistema de ecuaciones diferenciales lineales 3-dimensionales con matrices antisimétricas definidas en un cuerpo diferencial no constante. Al mismo sistema se aplican métodos y algorítmos propios de la Teoría de Galois Diferencial, lo que permite resolverlo algebraicamente y métodos numéricos, en particular métodos de la familia de Runge - Kutta. Por último, se calculan los errores absolutos y relativos entre las soluciones Liouvillianas, obtenidas mediante la resolución algebraica y las soluciones obtenidas aplicando métodos numéricos.
Palabras claves: Teoría de Galois diferencial, métodos de la familia de Runge - Kutta, soluciones Liouvillianas, sistemas de ecuaciones diferenciales, matrices antisimétricas, cuerpo diferencial no constante.
This work contrasts numerical methods with algebraic methods. These methods are applied to solve a three dimensional linear differential system with skew symmetric matrices defined in a non- constant differential field. Algorithms and methods of Differential Galois Theory, are used to provide an algebraic solution, while numerical methods, in particular, methods from Runge - Kutta family, are applied to the same system. Finally, the absolute and relative errors between Liouvillians solution are calculated comparing the solutions obtained by means of algebraic methods and by means of numerical methods.
Este trabajo es una versión mejorada de [12], la cual corresponde a la tesis de maestría del segundo autor, dirigida por el primer y el tercer autor en la Universidad del Norte.
El objetivo de este trabajo es contrastar mé todos numéricos con métodos algebraicos en la resolución de sistemas de ecuaciones diferenciales lineales 3-dimensionales con matrices antisimétricas definidas en un cuerpo diferencial no constante.
Para lograr este objetivo se hizo un estudio de las formas de resolución de sistemas de ecuaciones diferenciales lineales cuando la matriz de coeficientes es constante y se muestra que en el caso que los coeficientes del sistema de ecuaciones diferenciales pertenezca a un cuerpo diferencial no constante, el método de solución útil cuando los coeficientes son constantes, sólo es aplicable a un tipo muy restringido de matrices con entradas en un cuerpo diferencial no constante.
Ante esta realidad se hace uso de la Teoría de Galois Diferencial, [1,15]. El sistema original se somete a una serie de transformaciones a trav és de las cuales el sistema se convierte en una ecuación diferencial lineal reducida a la cual puede aplicarse el algoritmo implementado por Jerry Kovacic, [2]. Para el desarrollo de estas transformaciones se aplicaron conceptos expuestos por Emili Picard, Ernest Vessiot y Jean Gaston Darboux, en trabajos relativos a estos temas. Es un hecho conocido que só lo una clase bastante pequeña de ecuaciones diferenciales o de sistemas de las mismas tiene solución mediante métodos analíticos, ante esta circunstancia los métodos numéricos proveen una herramienta poderosa y adecuada que permite una aproximación a la solución bajo ciertas condiciones que debe cumplir el sistema. El estudio del sistema objetivo permitió comprobar que es posible resolverlo mediante la aplicación de métodos numéricos. Al sistema de ecuaciones se le aplicó métodos de la familia de Runge-Kutta, en particular mé todos de cuarto orden que en Matlab [14] se especifican como ODE45, [5,6,7,8,11,13,16]. Recientemente se han aplicado técnicas algebraicas y numéricas usando la algebrización Hamiltoniana para resolver un sistema de ecuaciones no-lineales de primer orden, ver [4].
Un sólido rígido es un conjunto de puntos del espacio que conservan la distancia entre ellos, bajo la acción de fuerzas aplicadas. Cualquier movimiento de un sólido rígido, es un movimiento rígido. Los movimientos rígidos pueden ser traslaciones o rotaciones. En una traslación, dos puntos cualesquiera del conjunto de puntos se mueven siguiendo una trayectoria paralela. Una rotación se realiza alrededor de un eje que bien puede pasar por puntos del sólido o por puntos que no forman parte del mismo, dos puntos cualesquiera que se encuentren sobre el mismo vector, cuyo origen es un punto del eje de rotación se mueven con la misma velocidad angular, en general, el movimiento de un sólido rígido es una combinación de rotación y traslación. Todo movimiento rígido es una isometría y por tanto puede ser descrito mediante transformaciones algebraicas, si V y W son espacios vectoriales, cada uno con un producto interno y dV es la métrica asociada al producto interno en V y dW es la métrica asociada al producto interno en W, una isometría es una función de V en W tal que
|
Un movimiento de un sólido rígido en el espacio tridimensional se puede describir mediante un grupo de isometrías uniparamétricas, eligiendo como parámetro el tiempo. Sea una función analítica, I un intervalo en R, la longitud de arco entre los puntos a y b de I se define como
|
El vector tangente unitario de la curva en R3 cuando esta ha sido parametrizada por su longitud de arco es,
|
La función de curvatura de es
|
El vector normal principal de es
|
y el vector binormal de es
|
Además, existe una función escalar llamada la torsión de , tal que . Si e1; e2; e3 son vectores unitarios ortogonales dos a dos, y p R3, un sistema ortonormal de coordenadas con origen en p se simboliza como: (p; e1; e2; e3) y es un sistema de referencia en p.
Los vectores , , , son vectores unitarios, ortogonales entre si, y junto con forman un sistema de referencia para la curva conocido como sistema de referencia de Frenet para la curva. El siguiente teorema enuncia que es posible especificar una curva si se conoce su función de curvatura y su torsión:
Teorema Fundamental de Curvas: Dadas dos funciones diferenciables existe una curva regular tal que es la longitud de arco, es la curvatura y y es la torsión de . Además, si otra curva satisface las mismas condiciones, esta difiere de por un movimiento rígido, esto es, existe una transformacion lineal y ortogonal con determinante positivo y un vector tal que , ver [10].
El sistema diferencial lineal objeto de nuestro interés es:
|
con matriz de coeficientes antisimétrica definida en un cuerpo diferencial no constante, este sistema modela el movimiento de un sólido rígido en el espacio vectorial o en un espacio afín. Estamos interesados en mostrar dos enfoques, uno numérico, otro analitico, en la resolución del sistema en cuestión.
La teoría de Galois establece que es posible encontrar las raices de un polinomio , definido sobre un cuerpo , realizando operaciones aritméticas y radicaciones con sus coeficientes, siempre y cuando el grupo de Galois del polinomio sea resoluble, esto es, exista una cadena finita de subgrupos normales, . El grupo de Galois de un polinomio está formado por los homomorfismos que al actuar sobre las raices del polinomio producen también raices.
La teoría de Galois diferencial establece proposiciones análogas para las ecuaciones diferenciales: permite decidir cuándo una ecuación diferencial se puede resolver a partir de sus coeficientes encontrando el correspondiente grupo resoluble.
Para determinar este grupo, es necesario desarrollar un proceso de algebrización de la ecuación diferencial, el cual consiste en expresar las funciones coeficientes de la ecuación diferencial como funciones racionales, el proceso de algebrización es posible mediante el procedimiento llamado cambio hamiltoniano de variable, [1]. Un cambio de variable es un cambio Hamiltoniano de variable, si y sólo si es una solución del sistema autónomo clásico Hamiltoniano con un grado de libertad, , para algún , elemento de un cuerpo diferencial con coeficientes en Si un cambio de variable es Hamiltoniano, existe tal que y por tanto, define una derivación para la nueva variable. Mediante el cambio hamiltoniano de variable las funciones transcendentes del sistema diferencial objetivo se transforman en funciones racionales, lo que permite la algebrización del proceso de solución del sistema diferencial, en efecto:
|
|
además,
|
entonces, el sistema se transforma en:
|
Mediante el uso de las coordenadas simétricas de Darboux, [10], para la esfera compleja : , el sistema algebrizado se transforma en una ecuación lineal diferencial reducida.
|
|
Si se deriva a con respecto a , se tiene que
|
(1) |
Teniendo en cuenta que:
|
La expresión (1) se transforma en una Riccati:
|
(2) |
Ahora, mediante el cambio , donde es una función de , por determinar, se tranformará la ecuación (2) en una ecuación diferencial de segundo orden
|
(3) |
donde . Para que (3) sea una ecuación diferencial lineal de segundo orden es necesario que , por tanto lo que finalmente conduce a:
|
(4) |
Ahora, la ecuación (4) puede transformarse en una EDLR [1]
|
El proceso de algebrización del sistema diferencial lineal ha permitido expresar el coeficiente de , , como una función racional, en efecto
|
.
Al transformar el sistema de ecuaciones diferenciales lineales con coeficientes no constantes en una EDL o en una EDLR con funciones racionales en , como coeficientes, es posible aplicar el Algoritmo de Kovacic para obtener una base de soluciones [3]. El Algoritmo de Kovacic clasifica las EDLR' s en cuatro tipos, éstos están relacionados con los ceros de , los de y el orden de en infinito, definido como: , además se define el conjunto como , donde . Para la ecuación que nos ocupa, y tienen 12 y 14 ceros, respectivamente, dos de los ceros de son de multiplicidad algebraica 2 el resto son de orden o multiplicidad algebraíca 1, todos los ceros de son de orden par, y es 2, por tanto, si existe una solución para la EDLR debe ser de la forma donde y se construyen con los pasos del algoritmo, ver [2]. La ecuación:
|
cae entonces en el caso 1 del Algoritmo de Kovacic y de existir la solución, tendrá la forma y será una solución liouvilliana [9].
Solucionar a mano, mediante el Algoritmo de Kovacic, EDLR es, en general, muy laborioso, normalmente la EDLR se utiliza para conocer el tipo de ecuación y la estructura de su correspondiente grupo de Galois Diferencial. En la práctica, para solucionar la ecuación diferencial, se aplica el programa"kovacicsols", implementado en Maple, a la EDL Las soluciones que se obtienen al aplicar "kovacicsols" a la ecuación son:
|
al algoritmo proveer dos soluciones se tiene una base de soluciones, y cualquier solución particular será de la forma:
|
La base de soluciones de la EDLR será:
|
La base de soluciones de la Riccati asociada a la EDL, será:
|
Aplicando el hecho que y que
|
es posible reconstruir la curva.
En las últimas décadas, Matlab [14] se ha convertido en un lenguage cada vez más importante para la computación científica, por ello en esta sección presentamos resultados numéricos, resueltos por medio de un código computacional generado en el lenguaje de Matlab.
Para que sea posible solucionar mediante métodos numéricos un sistema diferencial, es necesario que la función vectorial sea continua, acotada, y satisfaga una condición de Lipschitz con constante de Lipschitz en una región de . A continuación se prueba que la función vectorial,
|
cumple estas condiciones, por tanto, si es un punto cualquiera de , entonces, el problema de valores iniciales
|
tiene solución única sobre el intervalo y la gráfica de esta solución se encuentra en .
La función es de Lipschitz en
La función , definida mediante,
|
es de Lipschitz en su segunda variable, , en efecto, sea
|
|
Dado que
|
Si
|
, alcanza su único máximo relativo en y , se concluye que Sea entonces,
|
La norma infinito de la matriz está acotada por luego
|
la función satisface una condición de Lipschitz en la variable con constante de Lipschitz en una región
|
y cada una de las funciones son funciones acotadas en
El punto es un punto de , entonces el problema de valores iniciales
|
tiene solución única en
Para encontrar la solución numérica del sistema diferencial
|
y de la forma algebrizada del mismo, se usará el método de Runge-Kutta explícito de Dormand y Prince, DOPRI 5(4).
En la figura 1 se presenta la gráfica de la solución del sistema exponencial.![]() |
Figura 1: Solución del sistema exponencial. |
![]() |
Figura 2: Sol numérica versus Sol analítica S. Exp. |
![]() |
Figura 3: Error absoluto entre y . |
![]() |
Figura 4: Error absoluto entre y . |
![]() |
Figura 5: Error absoluto entre y . |
![]() |
Figura 6: Error relativo entre y . |
![]() |
Figura 7: Error relativo entre y . |
![]() |
Figura 8: Error relativo entre y . |
![]() |
Figura 9: Solución sistema algebrizado. |
![]() |
Figura 10: Sol. numérica versus Sol. analítica. S. Algebrizado. |
![]() |
Figura 11: Error absoluto entre y . |
![]() |
Figura 12: Error absoluto entre y . |
![]() |
Figura 13: Error absoluto entre y . |
![]() |
Figura 14: Error relativo entre y . |
![]() |
Figura 15: Error relativo entre y . |
![]() |
Figura 16: Error relativo entre y . |
Las gráficas que se obtienen al correr el programa que se diseñó, ponen de manifiesto una diferencia notable entre las soluciones numérica y analítica del sistema diferencial cuando el mismo está en forma exponencial.
Se observa que la diferencia existente entre las soluciones numérica y analítica, cuando se trabaja con la forma algebrizada del sistema diferencial, es imperceptible, gráficamente no hay distinción clara, el orden del error entre las dos soluciones es pequeño.
Al menos para este sistema, es evidente la necesidad de algebrizar previamente el sistema para que la solución numérica obtenida sea la adecuada.
El proceso que permite la algebrización del sistema diferencial hace uso de procesos proyectivos y de las coordenadas simétricas de Darboux, mediante estos procesos se obtuvieron las transformaciones expuestas. Este trabajo, al igual que [4], pone en evidencia que la algebrización Hamiltoniana es una herramienta útil, junto con los métodos numéricos, puesto que permiten minimizar los errores.
Los autores agradecen a la oficina de la DIDI de la Universidad del Norte, por el apoyo brindado para la publicación de este artículo. Los autores agradecen también a la Universidad Simón Bolívar. Finalmente los autores agradecen a los árbitros anónimos por su ayuda para mejorar este documento.
[16] T. Sauer, Análisis numérico. Pearson, Segunda edición, 2013.
Published on 03/01/18
Accepted on 07/05/17
Submitted on 27/03/17
Volume 34, Issue 1, 2018
DOI: 10.23967/j.rimni.2017.7.008
Licence: CC BY-NC-SA license
Are you one of the authors of this document?