Archivo de la etiqueta: Teorema de Picard

Ecuaciones Diferenciales I: Teorema de Existencia y Unicidad – Ecuación Integral, Funciones Lipschitzianas y Lema de Gronwall

Por Omar González Franco

Estudié matemáticas, la locura de la razón.
– Benjamin Moser

Introducción

A lo largo de esta primera unidad hemos estudiado una variedad de ecuaciones diferenciales ordinarias de primer orden y hemos desarrollado distintas técnicas para resolver cada tipo de ecuación. Vimos que una sola ecuación puede tener infinitas soluciones y sólo cuando le imponemos una condición inicial es como podremos obtener una solución particular de esa ecuación diferencial. Ahora bien, si la solución existe, entonces debe ser única pero, ¿es siempre cierto esto?.

Ya presentamos el teorema de existencia y unicidad para ecuaciones diferenciales de primer orden y el teorema de existencia y unicidad para el caso de ecuaciones diferenciales lineales de primer orden, nuestro objetivo ahora es tener un teorema de existencia y unicidad general que pueda aplicarse a cualquier ecuación diferencial ordinaria de primer orden.

Este teorema, conocido como teorema de existencia y unicidad de Picard – Lindelöf contiene las hipótesis suficientes para garantizar que si existe una solución a un problema de valor inicial (PVI), entonces dicha solución es única.

Cabe mencionar que es posible enunciar un teorema de existencia y unicidad de tipo global y uno de tipo local. En el caso de tipo global el intervalo de existencia de la solución se conoce a priori, mientras que en uno de tipo local se asegura que existe un intervalo, en un principio desconocido, donde el PVI tiene solución única. En este curso demostraremos el resultado de tipo global y veremos que el de tipo local es consecuencia del global, además de que puedes encontrar la demostración al teorema de tipo local en la sección de videos.

Demostrar el teorema de existencia y unicidad de Picard – Lindelöf no es tarea fácil, primero será necesario desarrollar una teoría preliminar en la que estableceremos algunos conceptos nuevos y, así mismo, haremos un breve repaso sobre conceptos que conocemos y que nos serán de utilidad para demostrar dicho teorema. Esta teoría preliminar la desarrollaremos a lo largo de esta y la siguiente entrada para finalmente demostrar el teorema en la última entrada de esta primera unidad.

Comenzaremos enunciando el teorema de existencia y unicidad de Picard – Lindelöf para tenerlo presente, a pesar de que quizá algunas cosas no queden claras, el objetivo de esta teoría preliminar será comprender lo que nos quiere decir este teorema, además de brindarnos las herramientas necesarias para demostrarlo.

Bien, ¡comencemos!.

Teorema de Existencia y Unicidad de Picard-Lindelöf

El teorema global de existencia y unicidad para ecuaciones diferenciales ordinarias de primer orden es el siguiente.

Una observación importante es que el punto (x0,y0) puede estar en la frontera de la banda vertical U=[a,b]×R, es decir, puede ser de la forma (a,y0) o (b,y0).

Podemos notar que en el enunciado se hace mención de términos que aún no conocemos, como lo son función lipschitziana e Iterantes de Picard, así que necesitamos definirlos.

Este teorema corresponde al resultado global en el que el intervalo es una banda vertical U=[a,b]×R, en el caso local se considera una región limitada definida como

R={(x,y)R2:|xx1|a,|yy1|b,a,bR}

y la solución esta definida en el intervalo δ=[x0h,x0+h] para cierta hR. Una vez demostrado el resultado global retomaremos el caso local.

En esta teoría preliminar veremos que el PVI (1) puede ser equivalente a resolver una ecuación integral, estudiaremos las funciones lipschitzianas de una y dos variables, demostraremos algunas proposiciones al respecto, demostraremos el lema de Gronwall, repasaremos algunos conceptos importantes sobre sucesiones, series y convergencia, definiremos las iteraciones de Picard y veremos algunos ejemplos. Una vez desarrollada esta teoría pasaremos a la demostración del teorema de existencia y unicidad.

Para comenzar, veamos que el PVI (1) se puede escribir de forma equivalente como una ecuación integral cuando la función f es continua.

Ecuación integral equivalente a un PVI

Un PVI como (1) se puede escribir de forma equivalente como una ecuación integral en el caso en el que la función f sea continua. Evidentemente este no es un curso ecuaciones integrales, pero para entender esta equivalencia definiremos lo que es una ecuación integral.

Teniendo en cuenta esta definición demostremos nuestro primer teorema de esta teoría preliminar el cual refleja el hecho de que un PVI como (1) es equivalente a resolver una ecuación integral.

Demostración:

Supongamos que y:δR, con gráfica contenida en U, es solución del PVI, entonces cumple que

dydx=f(x,y);y(x0)=y0

Como y es solución de la ecuación diferencial en el intervalo δ, entonces debe ser continua en el mismo intervalo, así tenemos que f y y son continuas y por tanto dydx y la función

g:δR,tg(t)=f(t,y(t))

también son continuas, de esta manera podemos integrar la ecuación dydx=f(x,y) para cualquier xδ.

x0xdydx(t)dt=x0xf(t,y(t))dt

Aplicando el teorema fundamental del cálculo (regla de Barrow) en el lado izquierdo de la ecuación, tenemos

y(x)y(x0)=x0xf(t,y(t))dty(x)=y(x0)+x0xf(t,y(t))dty(x)=y0+x0xf(t,y(t))dt

obteniendo así que y(x) verifica la ecuación integral (3).

Ahora supongamos que y(x) es una función continua en δ y que satisface la ecuación integral

y(x)=y0+x0xf(t,y(t))dt

Derivemos esta expresión.

dydx=ddx(y0+x0xf(t,y(t))dt)=dy0dx+ddx(x0xf(t,y(t))dt)=0+f(x,y)=f(x,y)

Donde se ha aplicado el teorema fundamental del cálculo. Con este resultado vemos que se ha recuperado la ecuación diferencial dydx=f(x,y), mostrando así que y(x) es solución a la ecuación diferencial y además

y(x0)=y0+x0x0f(t,y(t))dt=y0+0=y0

es decir, se satisface la condición inicial y(x0)=y0, de esta manera queda demostrado que y(x) es solución del PVI.

◻

Este resultado es muy útil en muchos resultados sobre ecuaciones diferenciales y nos será de utilidad para motivar, más adelante, la introducción a las llamadas iterantes de Picard.

Continuando con nuestra teoría preliminar, un concepto sumamente importante que estudiaremos a continuación es el de funciones lipschitzianas.

Funciones Lipschitzianas

Como estamos trabajando con la ecuación diferencial

dydx=f(x,y)

la función f es una función de dos variables, así que nos interesa estudiar las funciones lipschitzianas de dos variables, sin embargo es probable que este sea un concepto nuevo y para que sea más intuitivo entenderlo presentaremos la definición de función lipschitziana para el caso de una función de una variable y realizaremos algunos ejemplos sencillos para posteriormente definir la función lipschitziana en el caso de dos variables.

Con esta definición observamos que si x1x2 el cociente

f(x1)f(x2)x1x2

corresponde a la pendiente de la recta secante a la gráfica de f que pasa por los puntos (x1,f(x1)) y (x2,f(x2)), de esta forma la condición de Lipschitz indica que todas estas pendientes están acotadas, es decir, existe una constante L>0, tal que

|f(x1)f(x2)x1x2|L

para cada x1,x2I, con x1x2.

Recta secante que une a los puntos (x1,f(x1)) y (x2,f(x2)).

No entraremos es muchos detalles para el caso de una función de una variable, pero cabe mencionar que cualquier función lipschitziana es uniformemente continua, ya que dado ε>0 basta tomar δ=εL y la condición de Lipschitz (4) para que se verifique que

|x1x2|<δ|f(x1)f(x2)|<ε

Como ejemplo mostremos que toda recta es una función lipschitziana.

Ejemplo: Mostrar que la función

f(x)=mx+b

es una función lipschitziana, con L=|m|.

Solución: Queremos probar que se cumple (4). Vemos que

|f(x1)f(x2)|=|mx1+b(mx2+b)|=|mx1+bmx2b|=|mx1mx2|=|m||x1x2|=L|x1x2|

En donde consideramos que L=|m|. En este caso se da la igualdad

|f(x1)f(x2)|=L|x1x2|

probando así que la función f(x)=mx+b es una función Lipschitziana.

◻

Hay funciones uniformemente continuas que no son lipschitzianas, un ejemplo puede ser la función f:[0,1]R definida como f(x)=x, esta función es uniformemente continua pero no lipschitziana. Mostremos este hecho.

Ejemplo: Mostrar que la función f:[0,1]R, definida como f(x)=x no es lipschitziana.

Solución: Vamos a suponer que la función f(x)=x es lipschitziana y lleguemos a una contradicción. Si f(x)=x fuera lipschitziana debería satisfacer que

|f(x1)f(x2)|L|x1x2|

x1,x2[0,1] y para alguna L0. Vemos que

|f(x)f(0)|=|x0|L|x0|

es decir, x[0,1] (x es positiva),

xLx

Si x(0,1] (x0), entonces

xxL1xL

Este último resultado nos dice que la función 1x esta acotada por L para x(0,1], sin embargo si tomamos el límite x0 por la derecha obtenemos

limx0+1x=!

Hemos llegado a una contradicción y todo ocurrió de suponer que la función f(x)=x era lipschitziana. Por lo tanto, a pesar de ser uniformemente continua, f(x)=x no es lipschitziana.

◻

Un resultado más que no demostraremos es el siguiente teorema.

Hay funciones lipschitzianas que no son derivables, por ejemplo la función f:RR definida por f(x)=|x|.

Podemos decir, entonces, que la condición de Lipschitz es una condición intermedia entre continuidad uniforme y la existencia de derivada acotada.

Con esto en mente, ahora definamos lo que es una función lipschitziana para el caso en el que la función f es de dos variables. Para este caso, la condición de Lipschitz sólo afectará a una de las variables, concretamente a la segunda, importante considerar este hecho.

La relación (5) es lo que se pide que se cumpla en la tercer hipótesis del teorema de Picard – Lindelöf.

Enunciemos dos proposiciones importantes con respecto a las funciones lipschitzianas de dos variables que nos serán de utilidad a la hora de demostrar el teorema de Picard – Lindelöf.

Demostración: Sea f(x,y) una función lipschitziana respecto de la variable y y supongamos que existe su derivada parcial con respecto a dicha variable fy. Por definición, para (x,y)U se tiene que

(7)fy(x,y)limh0f(x,y+h)f(x,y)h

Dado un δ>0 y para h suficientemente pequeño |h|<δ, el punto (x,y+h) pertenece a U. Sea L una constante de Lipschitz de f respecto de y en U. De acuerdo a la definición de la condición de Lipschitz se verifica que

(8)|f(x,y+h)f(x,y)|L|y+hy|=L|h|

Usando (7) y (8) tenemos lo siguiente.

|fy(x,y)|=|limh0f(x,y+h)f(x,y)h|=limh0|f(x,y+h)f(x,y)h|limh0L|h||h|=L

Esto es,

|fy(x,y)|L

lo que significa que fy esta acotada en U por la constante de Lipschitz L.

◻

Ahora revisemos el resultado recíproco de la proposición anterior en donde es necesario que U sea un conjunto convexo.

Demostración: Para demostrar esta proposición haremos uso del teorema del valor medio para funciones de una variable, de aquí la necesidad de que U sea convexo.

Por hipótesis, fy esta acotada en U, sea L>0, tal que

(9)|fy(x,y)|L

para cada (x,y)U, y sean (x,y1),(x,y2)U con y1<y2. Como U es convexo tenemos garantizado que para cada y tal que y1<y<y2 el punto (x,y) pertenece a U, pues dicho punto pertenece al segmento que une los puntos (x,y1) y (x,y2), con estos resultados la función

gx:[y1,y2]R,gx(y)=f(x,y)

está bien definida y es derivable

gx(y)=fy(x,y)

para cada y[y1,y2]. Por el teorema del valor medio, existe y tal que y1<y<y2 y tal que

gx(y1)gx(y2)=gx(y)(y1y2)

es decir,

f(x,y1)f(x,y2)=fy(x,y)(y1y2)

Esta igualdad también la podemos escribir como

(10)|f(x,y1)f(x,y2)|=|fy(x,y)||y1y2|

Por la desigualdad (9), tenemos

(11)|fy(x,y)||y1y2|L|y1y2|

De los resultados (10) y (11) concluimos que

|f(x,y1)f(x,y2)|L|y1y2|

lo que prueba que f es una función lipschitziana con respecto de la segunda variable.

◻

Esta proposición es bastante útil, pues basta verificar que la derivada fy de f=f(x,y) esta acotada en un conjunto convexo U para concluir que f es una función lipschitziana respecto de la segunda variable. Realicemos un ejemplo.

Ejemplo: Sea U=[1,1]×R. Mostrar que la función f:UR definida como

f(x,y)=|x|sin2(y)

es una función lipschitziana respecto de la segunda variable.

Solución: Es claro que el conjunto U es convexo y que existe la derivada de f con respecto a y dada por

fy=2|x|sin(y)cos(y)

Como

|sin(y)cos(y)|1

yR y |x|<1,x[1,1], notamos que

2|x||sin(y)cos(y)|2

Esto es,

|fy|2

esto muestra que la derivada de f esta acotada, por la proposición anterior concluimos que la función f es lipschitziana y podemos tomar como constante de Lipchitz el valor L=2.

◻

En este ejemplo vimos que el valor L=2 es una cota de |fy|, sin embargo cualquier número mayor a 2 cumple también la desigualdad y por tanto también puede ser una constante de Lipschitz en U. En general, una buena constante de Lipschitz puede ser

(12)L=sup(x,y)U|fy(x,y)|

De ambas proposiciones podemos realizar la siguiente caracterización de Lipschitz, bastante útil en la práctica.

En este corolario unimos los resultados de las dos proposiciones anteriores.

Con esto concluimos el estudio de las funciones lipschitzianas, es importante tener presente este último corolario ya que será de suma relevancia en la demostración del teorema de Picard.

Para concluir con esta entrada presentaremos una herramienta más que nos será de mucha utilidad a la hora de demostrar el teorema de Picard – Lindelöf, en particular nos ayudará a probar la unicidad de la solución al PVI (1). Revisemos el Lema de Gronwall.

Lema de Gronwall

Este resultado fue desarrollado por Thomas Hakon Grönwall en 1919.

Demostración: Definamos la función

(15)g(x)=x0xf(t)dt

Notemos que

g(x0)=0ydgdx=f(x)

En términos de g(x) y dgdx la desigualdad (13) se puede escribir de la siguiente forma.

0dgdxα+βg(x)

de donde,

(16)dgdxβg(x)α

Multipliquemos ambos lados de la desigualdad por eβ(xx0).

eβ(xx0)(dgdxβg(x))eβ(xx0)α(17)eβ(xx0)dgdxβeβ(xx0)g(x)αeβ(xx0)

Identificamos que el lado izquierdo de la última desigualdad corresponde a la derivada del producto de las funciones eβ(xx0) y g(x), en efecto

ddx(g(x)eβ(xx0))=dgdxeβ(xx0)+g(x)(βeβ(xx0))=eβ(xx0)dgdxβeβ(xx0)g(x)

Sustituimos en la desigualdad (17).

(18)ddx(g(x)eβ(xx0))αeβ(xx0)

Integremos de x0 a x.

x0xddt(g(t)eβ(tx0))dtαx0xeβ(tx0)dtg(x)eβ(xx0)g(x0)eβ(x0x0)α[1β(eβ(xx0)eβ(x0x0))]

pero

g(x0)=x0x0f(t)dt=0yeβ(x0x0)=1

Así,

(19)g(x)eβ(xx0)αβ(eβ(xx0)1)

Multipliquemos ambos lados de la desigualdad por eβ(xx0).

g(x)αβeβ(xx0)(eβ(xx0)1)=αβ(1eβ(xx0))=αβ(eβ(xx0)1)

es decir,

(20)g(x)αβ(eβ(xx0)1)

De la desigualdad original (13) sabemos que

0f(x)α+βx0xf(t)dt0f(x)α+βg(x)

de donde,

(21)f(x)αβg(x)

De los resultados (20) y (21), tenemos

f(x)αβg(x)αβ(eβ(xx0)1)

lo que nos interesa es la desigualdad

f(x)αβαβ(eβ(xx0)1)

haciendo un poco de álgebra obtenemos lo siguiente.

f(x)αβαβ(eβ(xx0)1)f(x)αβαβ(eβ(xx0)1)f(x)α+βαβ(eβ(xx0)1)f(x)α+α(eβ(xx0)1)f(x)α+αeβ(xx0)αf(x)αeβ(xx0)

Por lo tanto,

f(x)αeβ(xx0)

Con esto queda demostrado que si se cumple la desigualdad (13), entonces f(x)αeβ(xx0), xI.

◻

Usando el lema de Gronwall podemos demostrar el siguiente corolario de manera inmediata.

Demostración: Debido a que se cumplen todas las hipótesis del lema de Gronwall sabemos que xI

0f(x)αeβ(xx0)

Pero si α=0, entonces

0f(x)0

de donde se deduce que f(x)=0, xI.

◻

Con esto concluimos la primer entrada sobre la teoría preliminar que necesitamos conocer para poder demostrar el teorema de existencia y unicidad de Picard – Lindelöf.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  1. Probar que la función f:RR, f(x)=c es una función lipschitziana
  1. Probar que la función f:RR, f(x)=|x| es lipschitziana, con L=1
  1. Probar que la función f:RR, f(x)=x2 no es una función lipschitziana.
    Hint: Suponer que lo es, es decir |f(x2)f(x1)|L|x2x1| y considerar la definición de derivada limx2x1|f(x2)f(x1)||x2x1|=|f(x1)| para llegar a una contradicción.

En los siguientes ejercicios se puede usar la definición de función lipschitziana respecto de la segunda variable o las proposiciones vistas.

  1. Probar que la función f:UR con U={(x,y):0x1,yR} definida como f(x,y)=ycos(x) es una función lipschitziana respecto de la segunda variable, con L=1.
  1. Probar que la función f:UR con U={(x,y):1x2,yR} definida como f(x,y)=2xy+exsin(x) es una función lipschitziana respecto de la segunda variable, con L=2.

Más adelante…

En esta entrada conocimos el teorema de existencia y unicidad de Picard – Lindelöf para ecuaciones diferenciales ordinarias de primer orden. Vimos que el PVI (1) es equivalente a resolver la ecuación integral (3), definimos a las funciones lipschitzianas de dos variables, demostramos algunos resultados al respecto y concluimos con la demostración del lema de Gronwall. Todos estos resultados los aplicaremos más adelante en la demostración del teorema de Picard – Lindelöf.

En la siguiente entrada continuaremos desarrollando esta teoría preliminar. Definiremos el concepto de aproximaciones sucesivas, mejor conocidas como iterantes de Picard, haremos un breve repaso sobre convergencia de series y sucesiones de funciones, presentaremos el resultado local del teorema de existencia y unicidad y resolveremos un ejercicio al respecto.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»