OBJETIVO
Llevar a cabo una anlisis de sistemas de ecuaciones diferenciales ordinarias
planteadas para la resolucin de un problema de cristalizacin de nitrato de
potasio dentro de un laboratorio de fisicoqumica, mediante la aproximacin
de datos usando el mtodo de Runge-Kutta de cuarto orden para sistemas de
ecuaciones, se comparan los resultados obtenidos con el modelo matemtico
trabajado con el algoritmo 5.7 de la pgina del doctor Burden.
MARCO TEORICO
Teorema de Picard-Lindelf
El teorema de Picard-Lindelf es un resultado matemtico de gran importancia dentro del
estudio de las ecuaciones diferenciales ordinarias (EDOs). Establece bajo qu condiciones
puede asegurarse la existencia y unicidad de solucin de una EDO dado un problema de
Cauchy (problema de valor inicial).
Teorema
El teorema debe su nombre al matemtico francs Charles mile Picard y al toplogo fins
Ernst Leonard Lindelf, ste ltimo enunci la teora de Picard tras su muerte.
Enunciado general
Sea
f ( t , x ) : R x R n Rn
x
Lipschitz respecto de
donde
(interprtese
es abierto, una funcin continua y localmente
f (t , x )
como la forma estndar de una EDO n-
dimensional de primer orden). Entonces, dado
I =[ t 0 , t 0+ ] R , R
intervalo cerrado
( t 0 , x 0 )
, podemos encontrar un
donde existe una nica solucin del
problema de Cauchy:
Que cumple que los pares
( t , x (t ) ) , t I
De hecho, ste puede ser encontrado de manera explcita, en la demostracin se dan detalles de ello.
Un enunciado ms restrictivo
El resultado anterior exige los requisitos mnimos que debe cumplir una funcin si queremos aplicar el
teorema. Aadiendo ms condiciones al enunciado original, podemos dar este otro ms sencillo: "Sea
n
f ( t , x ) : R x R R
nica solucin x (t)
Definida
una funcin Lipschitz. Entonces, dados
( t 0 , x 0 ) [ a , b ] x Rn
" existe una
del problema de valor inicial.
t [ a , b ] .
Observacin
Es importante observar que el teorema de Picard slo nos garantiza la existencia y
unicidad local de la solucin de una EDO. Es decir, ms all del intervalo proporcionado
por el teorema (dado que su demostracin es constructiva) no podemos decir nada, en
principio, del comportamiento de la solucin del problema de valor inicial. Es posible
complementar el teorema sealando que existe un intervalo abierto, que llamaremos
intervalo maximal en el cual puede garantizarse que la solucin existe y es nica; fuera de
este intervalo, el teorema de Picard no puede aplicarse.
Demostracin
Sea
) x B (x
)
C a ,b= I a (t
0
b
0
t I a ( t 0 )=[ t 0 ,t 0 + ]
el cilindro compacto donde
)= x b , x + b
B b (x
[ 0
] . Sea
0
0
mxima pendiente en mdulo. Y finalmente sea
est definida, esto es
M =f , s decir, el valor de
la constante de Lipschtitz de
respecto la segunda variable.
Definimos el siguiente operador entre funciones continuas, el operador de Picard, como
sigue:
definido como:
Vamos a imponer que est bien definido, es decir, que su imagen sea una funcin que
tome valores en
B b ( x 0) , es decir, que la norma de
sea menor que
El
imposicin, por lo que deber ser que
b .
ltimo
paso
es
Veamos ahora que el operador de Picard es contractivo bajo ciertas hiptesis sobre que
ms adelante podrn ser omitidas.
Dadas dos funciones
Pero como
queremos:
es Lipschitz respecto la segunda variable tenemos que:
Esto es contractivo si
o equivalentemente para tener igualdad si
Por lo tanto como el operador de Picard es un operador entre espacios de Banach (en
particular espacios mtricos inducidos por la norma) y contractivo, por el teorema del
punto fijo de Banach, existe una nica funcin
tal que
es decir,
solucin del problema de valor inicial definida en
donde debe satisfacer las
condiciones dadas, es decir,
Mtodos de Runge Kutta de Cuarto Orden
Un procedimiento de Runge-Kutta de Cuarto orden cosiste en determinar parmetros de
modo que la frmula
(1)
Donde
(2)
Concuerda con un polinomio de Taylor de grado cuatro. Esto da como resultado un sistema
de 11 ecuaciones con 13 incgnitas. El conjunto de valores usado con ms frecuencia para
los parmetros produce el siguiente resultado:
(3)
Mientras que las otras frmulas de cuarto orden se deducen con facilidad, el agorado
resumido es muy usado y reconocido como una invaluable herramienta de clculo, se
denomina el mtodo de Runge-Kutta de cuarto orden o mtodo clsico de Runge-Kutta. Se
aconseja tener cuidado con las formulas ya que k 2 depende de k1, k3 depende de k2 y k4
depende de k3.
Mtodos de Runge Kutta de Cuarto Orden para Sistemas
de Orden Superior
Para un sistema de la forma
(4)
Se parece a
(5)
Donde:
(6)
COMENTARIO
CONCLUSIONES
BIBLIOGRAFIA
Universidad de Salamanca. [En lnea] [Citado el: 4 de Octubre de 2015.]
http://campus.usal.es/~mpg/Personales/PersonalMAGL/Docencia/TeoriaTema1MM.pdf
.
Burden, Richard L. 2011. Analisis Numerico. [ed.] Michelle Julet. Novena Edicion.
s.l. : Cengage Learning, 2011. pgs. 259-355.
R., Molero. Caminos UPM. [En lnea] [Citado el: 4 de Octubre de 2015.]
http://www2.caminos.upm.es/Departamentos/matematicas/Fdistancia/PIE/Analisis
%20matematico/Temas/C08_Existencia_Soluciones.pdf.
UNAM. Universidad Autonoma De Mexico. [En lnea] [Citado el: 4 de Octubre de
2015.]
https://www.uam.es/personal_pdi/ciencias/fquiros/Numerico2_03_04/capitulo2.pdf.
2014. Wikipedia. Teorema de Picard-Lindelf. [En lnea] 29 de Diciembre de 2014.
[Citado el: 04 de Octubre de 2015.] https://es.wikipedia.org/wiki/Teorema_de_PicardLindel%C3%B6f.
Zill, Dennis G. 2009. Ecuaciones Diferenciales con problemas de valores en la
frontera. [trad.] Anal Garcia. Septima Edicion. s.l. : Cengage Learning, 2009. pgs.
339-362.
