Teorema de existencia y unicidad de ecuaciones diferenciales

Por Lizbeth Fernández Villegas

Introducción

En la entrada anterior trabajamos con la ecuación diferencial dy(x)dx=y(x) con condición inicial y(0)=1. Al identificar propiedades enunciadas en el teorema de punto fijo de Banach encontramos su solución. En esta ocasión repetiremos el proceso para demostrar que la solución a una ecuación diferencial general existe y es única.

Primeramente, veamos un concepto.

Definición. Función localmente Lipschitz continua en la segunda variable. Sea (a,b)R y sea ΩR tal que Ω es abierto. Si F:(a,b)×ΩR es una función que satisface que para cada x0(a,b) y y0Ω existen δ0>0 y c>0 tales que [x0δ0,x0+δ0](a,b),[y0δ0,y0+δ0]Ω y además que si x[x0δ0,x0+δ0] y si y1,y2[y0δ0,y0+δ0] entonces
|F(x,y1)F(x,y2)|c|y1y2|
diremos que F es localmente Lipschitz continua en la segunda variable.

Solución a la ecuación diferencial dy(x)dx=F(x,y(x))

Sea dy(x)dx=F(x,y(x)) una ecuación diferencial con condición inicial y(x0)=y0 donde:

  1. F es una función localmente Lipschitz continua en la segunda variable.
  2. y es una función, al menos derivable, de variable x que manda valores reales en valores reales.
  3. x0 es un punto donde la y buscada toma valor y0.

Plan para resolverla con el teorema de punto fijo de Banach: Propondremos un espacio métrico completo (X,d) de funciones entre las cuales deberá estar la y buscada y una contracción ϕ:XX cuyo punto fijo sea la solución de la ecuación diferencial.

Sean δ0>0 y c>0 para F localmente Lipschitz continua como en la definición. Se dejará como ejercicio al lector probar que F restringida en [x0δ0,x0+δ0]×[y0δ0,y0+δ0] es continua. Como este conjunto es compacto, se sigue que F está acotada en este conjunto. Por lo tanto existe M>0 tal que para toda (x,y)[x0δ0,x0+δ0]×[y0δ0,y0+δ0] se cumple
|f(x,y)|M.

Sea δ tal que 0<δ<min{1c,δ0M}

Considera X:={fC0([x0δ,x0+δ],R):d(f,y0)δM}
donde y0 representa, en este caso, a la función constante que arroja el valor y0. Nota que X es un espacio cerrado en el espacio métrico C0([x0δ,x0+δ],R) que recordemos, tiene la propiedad de ser completo. Por lo visto en la última proposición de la entrada Espacios métricos completos concluimos que X también es completo.

Propongamos la contracción ϕ deseada

Si fX satisface la ecuación diferencial entonces para todo x[x0δ,x0+δ] se sigue:

df(x)dx=F(x,f(x))x0xdf(t)dtdt=x0xF(t,f(t))dtf(x)f(x0)=x0xF(t,f(t))dtf(x)=f(x0)+x0xF(t,f(t))dtf(x)=y0+x0xF(t,f(t))dt

Como buscamos que esta solución sea punto fijo de una contracción ϕ en X entonces ϕ(f(x))=f(x). La última igualdad nos lleva a proponer:

ϕ(f(x)):=y0+x0xF(t,f(t))dt

Nota que ϕ(f(x)) pertenece a C0([x0δ,x0+δ],R). Probaremos que también pertenece a X. Si x[x0δ,x0+δ], tenemos dos casos:

Si x0x

|ϕ(f(x))y0|=|x0xF(t,f(t))dt|x0x|F(t,f(t))|dt(xx0)M=δM

Si x<x0

|ϕ(f(x))y0|=|x0xF(t,f(t))dt|=|xx0F(t,f(t))dt|=|xx0F(t,f(t))dt|xx0|F(t,f(t))|dt(x0x)M=δM

De ambos casos podemos concluir que d(f,y0)δM, por lo tanto ϕ(f)X.

ϕ es contracción en X

Sean f,gX. Considera I=[x0δ,x0+δ] entonces si xI, tenemos dos casos.

Si x0x.

|ϕ(f(x))ϕ(g(x))|=|y0+x0xF(t,f(t))dt(y0+x0xF(t,g(t))dt)|=|x0xF(t,f(t))dtF(t,g(t))dt|x0x|F(t,f(t))dtF(t,g(t))|dtx0xc|f(t)g(t)|dt(xx0)cd(f,g)(δc)d(f,g)

Si x<x0.

|ϕ(f(x))ϕ(g(x))|=|y0+x0xF(t,f(t))dt(y0+x0xF(t,g(t))dt)|=|x0xF(t,f(t))dtF(t,g(t))dt|=|xx0F(t,f(t))dtF(t,g(t))dt|=|xx0F(t,f(t))dtF(t,g(t))dt|xx0|F(t,f(t))dtF(t,g(t))|dtxx0c|f(t)g(t)|dt(x0x)cd(f,g)(δc)d(f,g)

Por lo tanto, la distancia entre ϕ(f) y ϕ(g) se puede estimar como

d(ϕ(f(x)),ϕ(g(x)))=SupxI{|ϕ(f(x))ϕ(g(x))|}SupxI{δcd(f,g)}=(δc)d(f,g)

Sea α:=δc entonces α<1, por lo tanto ϕ es contracción en X.

Lo que hemos visto en esta entrada demuestra el siguiente:

Teorema. Picard-Lindelöf. Sea F:(a,b)×ΩR una función continua y localmente Lipschitz continua en la segunda variable. Entonces, dados x0(a,b) y y0Ω existe δ>0 tal que la ecuación diferencial
dy(x)dx=F(x,y(x)),y(x0)=y0
tiene una única solución en el intervalo [x0δ,x0+δ].

Generalización en Rn

Si ΩRn y F tiene su contradominio en Rn entonces la definición y el teorema quedan como sigue:

Definición. Función localmente Lipschitz continua en la segunda variable. Sea (a,b)R y sea ΩRn tal que Ω es abierto. Si F:(a,b)×ΩRn es una función que satisface que para cada x0(a,b) y y0Ω existen δ0>0 y c>0 tales que [x0δ0,x0+δ0](a,b),B(y0,δ0)Ω y además que si |xx0|δ0 y si y1,y2B(y0,δ0) entonces
F(x,y1)F(x,y2)cy1y2
diremos que F es localmente Lipschitz continua en la segunda variable.

Teorema. Picard-Lindelöf. Sea F:(a,b)×ΩRn una función continua y localmente Lipschitz continua en la segunda variable. Entonces, dados x0(a,b) y y0Ω existe δ>0 tal que la ecuación diferencial
dy(x)dx=F(x,y(x)),y(x0)=y0
tiene una única solución en el intervalo [x0δ,x0+δ].

En este caso el espacio completo donde podemos encontrar la solución es

X:={fC0([x0δ,x0+δ],Rn):fy0δM}

Donde δ y M se definen de forma análoga a la demostración anterior.

Más adelante

Pasaremos a la siguiente sección de esta asignatura con temas de compacidad. Aunque ya se han usado algunos resultados para el caso del espacio métrico euclidiano, mostraremos cómo el concepto puede generalizarse en otros espacios a partir de la topología que la métrica induce en ellos.

Tarea moral

  1. Sean δ0>0 y c>0 para F localmente Lipschitz continua como en la definición. Prueba que F restringida en [x0δ0,x0+δ0]×[y0δ0,y0+δ0] es continua.
  2. Sea F:R×RR tal que F(x,y)=3y2/3.
    a) Prueba que F no es localmente Lipschitz continua en la segunda variable.
    b) Prueba que para cualesquiera α<0<β, la función
    fα,β(x)={(xα)3si x α,0si αxβ,(xβ)3si xβ.
    Es diferenciable en R y es solución de
    dy(x)dx=3y2/3,y(0)=0.
    Así, si F no es localmente Lipschitz continua en la segunda variable la ecuación puede tener una infinidad de soluciones.
  3. Sea F:R×RR tal que F(x,y)=y2.
    a) Prueba que F es localmente Lipschitz continua en la segunda variable.
    b) Para α0 considera la ecuación
    dy(x)dx=y2,y(0)=1α.
    Prueba que f(x)=1xα es su solución en algún intervalo que contiene a 0.
    c) ¿Cuál es el intervalo máximo para el que esta ecuación tiene solución?

Enlaces

Deja una respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Este sitio usa Akismet para reducir el spam. Aprende cómo se procesan los datos de tus comentarios.