Fyskode Learning

Métodos topológicos para EDO no lineales · Posgrado inicial · 10 horas

Banach, contracción, derivada de Fréchet y función implícita

Puntos fijos con unicidad, prolongación de soluciones y ramas locales mediante diferenciación en espacios normados.

El teorema del punto fijo de Banach produce más información que un resultado puramente topológico: entrega existencia, unicidad, convergencia de la iteración y una estimación explícita del error. Su precio es una hipótesis métrica fuerte, la contracción. En ecuaciones diferenciales esa hipótesis suele valer sólo en un intervalo corto, una bola pequeña o después de elegir una norma ponderada.

La derivada de Fréchet permite reconocer contracciones locales y estudiar cómo cambia una solución cuando varían los datos. Si la linealización respecto de la incógnita es invertible, el teorema de la función implícita construye una rama local única. Cuando la linealización pierde invertibilidad, el método deja de aplicar y aparecen las preguntas de resonancia y bifurcación de la unidad siguiente.

Contracciones en espacios completos

Sea (X,d)(X,d) completo y T:XXT:X\to X tal que existe 0<q<10<q<1 con

d(Tx,Ty)qd(x,y)d(Tx,Ty)\le q\,d(x,y)

para todos x,yx,y. Si xn+1=T(xn)x_{n+1}=T(x_n), entonces

d(xn+1,xn)qnd(x1,x0).d(x_{n+1},x_n)\le q^n d(x_1,x_0).

La desigualdad triangular da, para m>nm>n,

d(xm,xn)k=nm1qkd(x1,x0)qn1qd(x1,x0).d(x_m,x_n) \le\sum_{k=n}^{m-1}q^k d(x_1,x_0) \le\frac{q^n}{1-q}d(x_1,x_0).

La sucesión es de Cauchy y converge a cierto xx^*. Por continuidad, Tx=xTx^*=x^*. Si yy^* fuera otro punto fijo,

d(x,y)=d(Tx,Ty)qd(x,y),d(x^*,y^*)=d(Tx^*,Ty^*)\le qd(x^*,y^*),

lo cual obliga a x=yx^*=y^*. Además se obtienen los controles

d(xn,x)qn1qd(x1,x0),d(x_n,x^*)\le\frac{q^n}{1-q}d(x_1,x_0), d(xn,x)q1qd(xn,xn1).d(x_n,x^*)\le\frac{q}{1-q}d(x_n,x_{n-1}).

El primero es una estimación a priori y el segundo una estimación a posteriori. Declarar convergencia sin cuantificar qq omite el dato que vuelve útil al teorema.

Iteración gráfica de un mapa hacia un punto fijo
La telaraña converge si el mapa contrae en la región pertinente. La imagen ilustra la iteración, mientras que la prueba requiere una cota uniforme de la constante Lipschitz menor que uno.

Atlas reproducible de sistemas dinámicos Funciones: plot_return_floquet, animate_return_cobweb

Descargar .py

Picard, elección del intervalo y prolongación

Para

x=f(t,x),x(t0)=x0,x'=f(t,x),\qquad x(t_0)=x_0,

definimos en C([t0h,t0+h],Rm)C([t_0-h,t_0+h],\mathbb R^m)

(Tx)(t)=x0+t0tf(s,x(s))ds.(Tx)(t)=x_0+\int_{t_0}^{t}f(s,x(s))\,ds.

Supongamos que f(t,x)M|f(t,x)|\le M y que ff es Lipschitz en xx con constante LL sobre el rectángulo

tt0a,xx0b.|t-t_0|\le a,\qquad |x-x_0|\le b.

Si hah\le a, MhbMh\le b y Lh<1Lh<1, entonces TT conserva la bola funcional xx0b\|x-x_0\|_\infty\le b y

TxTyLhxy.\|Tx-Ty\|_\infty\le Lh\|x-y\|_\infty.

Banach da una solución única. La restricción Lh<1Lh<1 no significa que la solución deje de existir después de hh. Se reinicia el argumento desde el extremo y se prolonga mientras la gráfica permanezca en un compacto donde ff sea regular. El teorema de prolongación se resume así: una solución maximal definida en [t0,b)[t_0,b) con b<b<\infty debe abandonar todo compacto del dominio cuando tbt\uparrow b. Si permanece acotada y lejos de la frontera del dominio, puede extenderse, contradiciendo maximalidad.

Una norma de Bielecki,

xα=maxteαtt0x(t),\|x\|_\alpha=\max_t e^{-\alpha|t-t_0|}|x(t)|,

puede reducir la constante del operador integral. Sin embargo, cambiar de norma no crea cotas globales sobre una no linealidad que crece sin control; la invariancia del conjunto sigue necesitando una estimación.

El tratamiento de Amster presenta Banach como origen simultáneo de Picard y de varios resultados locales posteriores, con atención a la elección del espacio y la bola. Amster, 2021 Pablo Amster (2021) Métodos topológicos en el estudio de las ecuaciones diferenciales no lineales 1.ª ed., Papirhos Textos 12 · Instituto de Matemáticas, UNAM Ubicación consultada: cap. 4, §§4.1–4.5: contracción, diferenciación de Fréchet y función implícita Abrir fuente

Diferenciación de Fréchet

Sean X,YX,Y espacios normados y F:UXYF:U\subset X\to Y. Decimos que FF es diferenciable en sentido de Fréchet en xx si existe un operador lineal continuo DF(x):XYDF(x):X\to Y tal que

limhX0F(x+h)F(x)DF(x)hYhX=0.\lim_{\|h\|_X\to0} \frac{\|F(x+h)-F(x)-DF(x)h\|_Y}{\|h\|_X}=0.

La aproximación debe ser uniforme respecto de la dirección h/hh/\|h\|. Poseer derivadas direccionales en todas las direcciones no basta. Esta uniformidad permite la regla de la cadena y el teorema de la función implícita.

Consideremos el operador de Nemytskii

N(u)(t)=g(t,u(t))N(u)(t)=g(t,u(t))

de C([0,T])C([0,T]) en sí mismo. Si gg posee derivada continua gxg_x sobre una banda compacta, entonces

DN(u)v(t)=gx(t,u(t))v(t).DN(u)v(t)=g_x(t,u(t))v(t).

Para justificarlo, escribimos mediante el teorema del valor medio integral

g(t,u+v)g(t,u)gx(t,u)v=01[gx(t,u+sv)gx(t,u)]vds.g(t,u+v)-g(t,u)-g_x(t,u)v =\int_0^1[g_x(t,u+sv)-g_x(t,u)]v\,ds.

La continuidad uniforme de gxg_x en la banda hace que el cociente por v\|v\|_\infty tienda a cero uniformemente en tt.

Si KK es lineal continuo, la derivada de KNK\circ N es KDN(u)K\circ DN(u). Para

F(u)=uKN(u),F(u)=u-KN(u),

resulta

DF(u)v=vK[gx(,u())v].DF(u)v=v-K[g_x(\cdot,u(\cdot))v].

Esta fórmula transforma la invertibilidad abstracta de DF(u)DF(u) en un problema lineal de contorno.

Teorema de la función implícita

Sea F:UX×RpYF:U\subset X\times\mathbb R^p\to Y de clase C1C^1, y supongamos

F(u0,λ0)=0.F(u_0,\lambda_0)=0.

Si la derivada parcial

DuF(u0,λ0):XYD_uF(u_0,\lambda_0):X\to Y

es un isomorfismo lineal continuo, existen vecindades de u0u_0 y λ0\lambda_0 y una función única u=u(λ)u=u(\lambda) de clase C1C^1 tal que

F(u(λ),λ)=0.F(u(\lambda),\lambda)=0.

Al derivar la identidad se obtiene

DuF(u(λ),λ)u(λ)+DλF(u(λ),λ)=0,D_uF(u(\lambda),\lambda)u'(\lambda) +D_\lambda F(u(\lambda),\lambda)=0,

y en el punto base

u(λ0)=[DuF(u0,λ0)]1DλF(u0,λ0).u'(\lambda_0) =-[D_uF(u_0,\lambda_0)]^{-1} D_\lambda F(u_0,\lambda_0).

La prueba reduce el problema a una contracción dependiente del parámetro. Sea A=DuF(u0,λ0)A=D_uF(u_0,\lambda_0). Resolver F(u,λ)=0F(u,\lambda)=0 equivale a buscar un punto fijo de

Tλ(u)=uA1F(u,λ).T_\lambda(u)=u-A^{-1}F(u,\lambda).

En una bola pequeña, DuTλ=IA1DuFD_uT_\lambda=I-A^{-1}D_uF tiene norma menor que uno porque se anula en (u0,λ0)(u_0,\lambda_0) y depende continuamente de (u,λ)(u,\lambda). También se verifica que TλT_\lambda conserva la bola. Banach produce el punto fijo único y su dependencia regular. Las clases abiertas de la UBA desarrollan esta construcción para hacer explícita la relación entre invertibilidad y unicidad local. Amster, 2020 Pablo Amster (2020) Análisis no lineal: métodos topológicos Curso de 24 clases con notas y videos · Departamento de Matemática, Universidad de Buenos Aires Ubicación consultada: clases sobre derivada de Fréchet y teorema de la función implícita Abrir fuente

Ejemplo trabajado: una rama de soluciones de Dirichlet

Estudiemos

u=λsinu,u(0)=u(1)=0.-u''=\lambda\sin u, \qquad u(0)=u(1)=0.

Sea

X=C2([0,1])C0([0,1]),Y=C([0,1]),X=C^2([0,1])\cap C_0([0,1]), \qquad Y=C([0,1]),

con la norma usual de C2C^2 en XX. Definimos

F(u,λ)=uλsinu.F(u,\lambda)=-u''-\lambda\sin u.

Se cumple F(0,λ)=0F(0,\lambda)=0 para todo λ\lambda: existe la rama trivial. La derivada respecto de uu es

DuF(0,λ)v=vλv.D_uF(0,\lambda)v=-v''-\lambda v.

Para λ0<π2\lambda_0<\pi^2, este operador es invertible. En efecto, si

vλ0v=0,v(0)=v(1)=0,-v''-\lambda_0v=0,\qquad v(0)=v(1)=0,

multiplicar por vv e integrar da

01v2dt=λ001v2dt.\int_0^1|v'|^2dt =\lambda_0\int_0^1|v|^2dt.

La desigualdad de Poincaré,

01v2dtπ201v2dt,\int_0^1|v'|^2dt\ge\pi^2\int_0^1|v|^2dt,

obliga a v=0v=0. La teoría lineal de Dirichlet da sobreyectividad. Por función implícita, cerca de (0,λ0)(0,\lambda_0) la única solución es la trivial.

En λ=π2\lambda=\pi^2, el núcleo contiene sin(πt)\sin(\pi t) y la hipótesis falla. No se concluye que exista una bifurcación sólo por esa falla; se concluye que la unicidad local ya no está garantizada por este teorema. La reducción y las condiciones adicionales necesarias pertenecen al análisis de bifurcación.

Para ver una rama no constante, modifiquemos el problema a

u=λsinu+h(t),-u''=\lambda\sin u+h(t),

y supongamos conocida una solución u0u_0 en λ0\lambda_0 para la cual vλ0cos(u0)v-v''-\lambda_0\cos(u_0)v es invertible. Entonces

u(λ0)=[d2dt2λ0cosu0]1sinu0.u'(\lambda_0) =\big[-\tfrac{d^2}{dt^2}-\lambda_0\cos u_0\big]^{-1}\sin u_0.

La expresión determina la tangente de la rama y su cálculo requiere resolver un problema lineal, no volver a resolver el problema no lineal completo.

Procedimiento funcional

  1. Elige un espacio completo cuya norma controle las cantidades presentes en la ecuación.
  2. Construye un conjunto cerrado preservado por el operador.
  3. Calcula una constante de contracción explícita y verifica que sea menor que uno.
  4. Registra una cota del error de iteración; no uses sólo la diferencia visual entre curvas.
  5. Para prolongar, demuestra que la solución permanece en compactos del dominio.
  6. Si hay parámetros, formula F(u,λ)=0F(u,\lambda)=0 y calcula la derivada de Fréchet con su resto uniforme.
  7. Traduce la invertibilidad de DuFD_uF a un problema lineal y comprueba núcleo y sobreyectividad.
  8. Distingue la rama local garantizada de cualquier afirmación global o bifurcación posterior.

Ejercicios

  1. Demuestra las dos estimaciones de error del teorema de Banach y compara cuál puede calcularse durante una iteración.
  2. Aplica Picard a x=t+x2x'=t+x^2, x(0)=0x(0)=0, y determina valores concretos de hh, bb, MM y LL que hagan invariante la bola.
  3. Prueba que una solución maximal de x=f(t,x)x'=f(t,x) puede prolongarse si su gráfica permanece en un compacto contenido en el dominio de ff.
  4. Demuestra la diferenciabilidad de Fréchet del operador N(u)=euN(u)=e^u sobre una bola de C[0,1]C[0,1] y estima el resto cuadrático.
  5. Calcula la derivada del operador integral F(u,λ)=uλK(u3)hF(u,\lambda)=u-\lambda K(u^3)-h y formula su problema lineal asociado.
  6. Verifica mediante series de Fourier que vλv-v''-\lambda v con Dirichlet es invertible exactamente cuando λn2π2\lambda\ne n^2\pi^2.
  7. Implementa la iteración de Picard y una corrección de Newton para un problema escalar; compara residuos y explica qué hipótesis locales justifican cada convergencia.

Banach aporta control métrico y la derivada de Fréchet convierte la ecuación no lineal en una familia lineal bien definida. El teorema implícito combina ambos: la inversa de la linealización crea una contracción local. Su alcance es preciso —rama única en una vecindad— y su fallo identifica, sin resolver todavía, los parámetros donde pueden surgir resonancia o bifurcación.

Fuentes consultadas

Obras citadas en el desarrollo; los localizadores indican los capítulos o secciones consultados.
  1. Pablo Amster (2021). Métodos topológicos en el estudio de las ecuaciones diferenciales no lineales. 1.ª ed., Papirhos Textos 12 Instituto de Matemáticas, UNAM. ISBN 978-607-30-5203-0.

    cap. 4, §§4.1–4.5: Banach, derivada de Fréchet y función implícita · Contracciones, prolongación, diferenciación de Fréchet y función implícita.
  2. Pablo Amster (2020). Análisis no lineal: métodos topológicos. Curso de 24 clases con notas y videos Departamento de Matemática, Universidad de Buenos Aires.

    curso abierto, clases sobre Banach y teorema de la función implícita · Construcción funcional y dependencia respecto de parámetros.
  3. Gerald Teschl (2012). Ordinary Differential Equations and Dynamical Systems. Graduate Studies in Mathematics, vol. 140; edición en línea autorizada American Mathematical Society. ISBN 978-0-8218-8328-0.

    caps. 2 y 5, existencia local, dependencia de datos y contorno · Apoyo para Picard, variación de constantes y ecuaciones lineales.
  4. Bjorn Poonen y Haynes Miller (2009). Honors Differential Equations, 18.034. Notas, problemas y soluciones; CC BY-NC-SA 4.0 MIT OpenCourseWare.

    notas abiertas sobre existencia, unicidad y problemas no lineales · Ejemplos complementarios de iteración y estimaciones.