Archivo de la categoría: Sin clasificar

Teorema de la Función Inversa

Por Angélica Amellali Mercado Aguilar

Teorema de la Función Inversa (sistema fi:RnR)

Teorema 1. Sea URn un abierto y sean
f1:URfn:UR
con derivadas parciales continuas. Considerar las ecuaciones
f1(x1,x2,,xn)=y1f2(x1,x2,,xn)=y2fn(x1,x2,,xn)=yn Tratamos de resolver las n-ecuaciones para x1,x2,xncomo funciones de y1,y2,yn.
La condición de existencia para la solución en una vecindad del punto x0 es que el determinante de la matriz Df(x0) y f=(f1,f2,fn) sean distintos de cero.

La condición de existencia para la solución en una vecindad del punto x0 es que el determinante de la matriz Df(x0) y f=(fi,f2,fn) sean distintos de cero. Explicitamente:

[(f1,f2,,fn)(x1,x2,,xn)|x=x0=J(f)(x0)=|f1x1(x0)f1x1(x0)fnx1(x0)fnxn(x0)|0]

entonces el sistema anterior se puede resolver de manera ‘unica como x=g(y) para x cerca de x0 y y cerca de y0 ◻

Nota. La cuestión de existencia se responde por medio del teorema general de la función implícita aplicado a las funciones yifi(x1,x2,,xn) con las incognitas x1,x2,,xn.

Ejemplo. El problema de factorizar un polinomio xn+an1xn1++a0 en factores lineales es, en cierto sentido un problema de función inversa. Los coeficientes ai son funciones conocidas de las n raices rj. ¿Se podran expresar las raices como funciones de los coeficientes en alguna región?. Con n=3 , aplicar el teorema de la función inversa a este problema y enunciar la conclusión acerca de la posibilidad de hacer lo planteado.

Solución. Para el caso n=3 tenemos que podemos factorizar el polinomio de la siguiente forma
x3+a2x2+a1x+a0=(xr1)(xr2)(xr3)
desarrolando el lado derecho tenemos que
(xr1)(xr2)(xr3)=x3r3x2r2x2+r2r3xr1x2+xr1r3+r1r2xr1r2r3
que se puede escribir
x3+x2(r3r2r1)+x(r2r3+r1r3+r1r2)r1r2r3
igualando las expresiones
x3+a2x2+a1x+a0=x3+x2(r3r2r1)+x(r2r3+r1r3+r1r2)r1r2r3
por lo tanto igualando coeficientes
a0=r1r2r2a1=r2r3+r1r3+r1r2a2=r1r2r3

Al sistema anterior le aplicamos el teorema de la función implicita para comprobar si las raices se pueden expresar en términos de los coeficientes, para ello calculamos el determinante de jacobiano del sistema que en este caso es

J=(a0r1a0r2a0r3a1r1a1r2a2r3a2r1a2r2a3r3)=(r2r3r1r3r1r2r3+r2r3+r1r2+r1111)

de esta manera el determinante del jacobiano es
det(r1r3r1r3r1r2r3+r2r3+r1r2+r1111)=|r1r3r1r3r1r2r3+r2r3+r1r2+r1111|=
(r2r3)×|r3+r1r2+r111|(r1r3)×|r3+r2r2+r111|+(r1r2)×|r3+r2r3+r131|= (r2r3)×(r2r3)+(r1r3)×(r1r3)(r1r2)×(r1r2)=r2r3r2+r2r3r3+r1r3r1r1r3r3r1r2r1+r1r2r2

que se puede escribir
=r3r1r1r3r1r3r3r2r1+r3r2r3r2r1r1+r2r1r3+r2r2r1r2r2r3
=(r3r1r3r2r2r1+r2r2)r1(r3r1r3r2r2r1+r2r2)r3=(r3r1r3r2r2r1+r2r2)(r1r3)
=((r3r2)r1(r3r2)r2)(r1r3)=(r3r2)(r1r2)(r1r3)

Este último término no es cero si el polinomio tiene raices distintas. Así el teorema de la función inversa muestra que las raices se pueden hallar como funciones de los coeficientes en alguna vecindad de cualquier punto en el que las raices sean distintas. Esto es, si las rices r1, r2, r3 de x3+a2x2+a1x+a0 son todas diferentes, entonces hay vecindades V de (r1, r2, r3) y W de (a0, a1, a2) tales que las raices en V son funciones de los coeficientes en W.

Funciones de RnRm

Definición 1. Una función f de Rn en Rm denotada f:RnRm, es una relación que asigna a cada vector del espacio Rn un único vector del espacio Rm\Si f es una función de Rn en Rm, entonces f se expresa
f=(f1,f2,,fm) en donde fk  k=1,,m es la k-ésima función componente y fk:RnR k=1,,m

Definición 2. Si ARn, la imagen bajo la función f de Rn en Rm se denota f(A), y se define
f(A)=f(x)Rm | xA

Definición 3. El dominio de una función f de Rn en Rm es la intersección de los dominios de las funciones componentes fk es decir
Domf=k=1mDomfk=Domf1Domf2Domf3Domfm

Ejemplo. Encontrar el dominio y la imagen de la recta y=3x para la función f:R2R2 dada por f(x,y)=(4x+2y5,2x+y5)
\item[Solución] En este caso f1=(4x+2y5)  Domf1=R2
f2=(2x+y5)  Domf2=R2
por lo tanto
Domf=Domf1Domf=R2R2=R2
Para la imagen de la recta y=3x procedemos de la siguiente manera
f(x,y)=(4x+2y5,2x+y5)=(x,y)  f(x,3x)=(4x+2(3x)5,2x+(3x)5)=(x,y)  
x=4x+2(3x)5yy=2x+(3x)5  x=2xyy=xy=x2
por lo tanto la imagen de la recta y=3x sera:
f(3x)={(x,y)R2 | y=x2}

Definición 4. Sean f:ARnRm y DRm. Definimos la imagen inversa de D bajo f, que denotamos f1(D), como el conjunto dado por:
f1(D)={x^A | f(x^)D}

Definición 5. Sean f:ARnRm y DRm, BA. Definimos la imagen directa de B bajo f, que denotamos f(B), como el conjunto dado por:
f(B)={f(x^)Rm | x^B}

Proposición 1. Sean f:ARnRm, Aα, B,CA y Dα, D, ERm, con αI I un conjunto de indices. Pruebe que:
1.DE  f1(D)f1(E)2.f1(αIDα)=αIf1(Dα) 3.f1(αIDα)=αIf1(Dα)4.f1(Dc)=(f1(D))c5.BC  f(B)f(C)6.f(αIAα)=αIf(Aα)7.f(αIAα)αIf(Aα)8.f(A)f(B)f(AB)9.Bf1(f(B))10.f(f1)(D)D

Demostración. (1).DE  f1(D)f1(E)
xf1(D)  f(x)D
 DE f(x)E   xf1(E) (2).f1(αIDα)=αIf1(Dα)
xf1(αIDα)  f(x)αIDα
  f(x)Dαip.a. iI   xf1(Dαi)p.a. iI
  x(αIf1(Dα)
(3).f1(αIDα)=αIf1(Dα)
xf1(αIDα)  f(x)αIDα
  f(x)Dαi iI   xf1(Dαi) iI
  x(αIf1(Dα)
(4).f1(Dc)=(f1(D))c
xf1(Dc)  f(x)Dc
  f(x)D
  xf1(D)
  x(f1(D))c
(5).BC  f(B)f(C)
f(x)f(B)  xB
 BC xC   f(x)f(C) (6).f(αIAα)=αIf(Aα)
f(x)f(αIAα)  xαIAα
  xAαp.a α I   f(x) f(Aα) p.a α I   f(x)αIf(Aα)

(7).f(αIAα)αIf(Aα) f(x)f(αIAα)  xαIAα   xAα α I
  f(x) f(Aα)  α I
  f(x)αIf(Aα)
(8).f(A)f(B)f(AB)
f(A)f(B)=f(A)(f(B))c
  f(x)f(A)f(B)
  f(x)f(A)(f(B))c
  f(x)f(A)yf(x)(f(B))c
  xAyf(x)f(B)
 xAyxB  xAyxBc
  xABc
  xAB
  f(x)f(AB)
(9).Bf1(f(B))xB  f(x)f(B)  xf1(f(B))$


(10).f(f1(D))D
f(x)f(f1(D))  xf1(D)
  f(x)D ◻

Diferenciación

Por Angélica Amellali Mercado Aguilar

Diferenciación de funciones f:RnRm

Definición. Considere la función f:ARnRm definida en un conjunto abierto A de Rn y sea x0A. Se dice que esta función es diferenciable si

f(x0+h)=f(x0)+f(x0)h+r(h)
cumple
limh0r(h)|h|=0^

Ejemplo. Compruebe que la función f:R2R3 definida por
f(x,y)=(exy,x2+y,2x3y2) es diferenciable en (1,3)
Solución En este caso
limh0r(h)|h|=
lim(h1,h2)(0,0)f(1+h1,3+h2)f(1,3)((3e3,e3)(h1h2),(2,1)(h1h2),(54,12)(h1h2))|(h1,h2)|
=lim(h1,h2)(0,0)(e(1+h1)(3+h2),(1+h1)2+(3+h2),2(1+h1)3(3+h2)2)(e3,4,18)|(h1,h2)|
((3e3,e3)(h1h2),(2,1)(h1h2),(54,12)(h1h2))|(h1,h2)|
=(lim(h1,h2)(0,0)e(1+h1)(3+h2)e33e3h1e3h2|(h1,h2)|,lim(h1,h2)(0,0)(1+h1)2+(3+h2)42h1h2|(h1,h2)|,
lim(h1,h2)(0,0)2(1+h1)3(3+h2)21854h112h2|(h1,h2)|)
=(0,0,0)
por lo que la función es diferenciable.

En el ejemplo anterior se tiene que
lim(h1,h2)(0,0)f(1+h1,3+h2)f(1,3)((3e3,e3)(h1h2),(2,1)(h1h2),(54,12)(h1h2))|(h1,h2)|
se puede expresar
lim(h1,h2)(0,0)f(1+h1,3+h2)f(1,3)(3e3e3215412)[h1h2]h12+h22
lo que nos lleva a la siguiente definición.

Definición. A la matriz de m×n se le llama Matriz Jacobiana de la función f en x0 y se le denota Jf(x0).


Definición. Sea f:RnRm definida en el abierto Ω de Rn y x0Ω. Se dice que f es diferenciable en x0Ω si y solo si existe una matriz T de m×n tal que limh0f(x0+h)f(x0)T(x0)h|h|=0donde T(x0) es la matriz jacobiana denotada por Jf(x0) ó Df(x0). En notación matricial:
T(x0)h=(f1(x0)x1f1(x0)xnfm(x0)x1fm(x0)xn)(h1 hn)

En términos ϵδ se tiene que si 0<|h|<δ entonces |f(x0+h)f(x0)T(x0)h||h|<ϵ

Teorema 1. Supónga que f:ΩRnRm es diferenciable en x0Rn. Entonces la matriz T es única

Demostración. Supongamos que existen T1 y T2 que cumplen limh0f(x0+h)f(x0)T1(x0)h|h|=0ylimh0f(x0+h)f(x0)T2(x0)h|h|=0

limh0f(x0+h)f(x0)T1(x0)h|h|f(x0+h)f(x0)T2(x0)h|h|=0
limh0T2(x0)hT1(x0)h|h|=0
Sea x un vector unitario en la dirección del vector h y hacemos h=tx con tR

0=limh0T2(x0)hT1(x0)h|h|=limt0T2(x0)txT1(x0)tx|tx|=limt0tT2(x0)xT1(x0)x|t||x|
0=T2(x0)xT1(x0)xT2(x0)x=T1(x0)xT2(x0)=T1(x0)T2=T1 ◻

Operadores: Divergencia, Rotacional y Laplaciano

Considere la función f:R3R3 dada por
f(x,y,z)=(f1(x,y,z),f2(x,y,z),f3(x,y,z))
cuya matriz jacobiana es
(f1xf1yf1zf2xf2yf2zf3xf3yf3z )
Con los elementos de esta matriz se forman importantes combinaciones que son la divergencia y el rotacional, conocidos también como invariantes de primer orden de esta matriz. La razón por la que se llaman invariantes es porque el valor de dichas combinaciones no se altera al efectuar un cambio de coordenadas.

Definición. Dada la matriz Jacobiana (f1xf1yf1zf2xf2yf2zf3xf3yf3z )
Se define la divergencia de f como
div (f)=f1x+f2y+f3z
Se puede ver que es igual a la suma de los elementos de la diagonal principal, es decir, que constituye la traza de la matriz jacobiana de f. La defnición de divergencia puede darse también mediante el operador (nabla)
f=(x,y,z)(f1,f2,f3)=f1x+f2y+f3z

Dada la matriz Jacobiana
(f1xf1yf1z f2xf2yf2z f3xf3yf3z )
se define el rotacional como
rot (f)=(f3yf2z,f3xf1z,f2xf1y)
Y podemos ver que las componentes del rotacional están definidas por las diferencias de los elementos situados simétricamente con respecto a la diagonal principal de la matriz Jacobiana.
La defnición de rotacional puede darse también mediante el operador (nabla)
×f=|i^j^k^xyzf1f2f3|=(f3yf2z,f3xf1z,f2xf1y)

Definición. Sea f:RnR definida en el abierto Ω de Rn tal que f es de clase C2 en Ω. La expresión
2f=2fx2+2fy2
es llamada Laplaciano de f. La ecuación
2fx2+2fy2=0
es llamada la ecuación de Laplace. Las funciones f de clase C2 que cumplen la ecuación de Laplace se llaman funciones Armónicas.

Ejercicio. Sean f,g:AR3R3 dos funciones diferenciables en una región AR3 y φ:R3R
Pruebe que

(a) div (f+g)=div f+div g

(b) div (φf)=φ div (f)+grad (φ)f
(c) rot (f+g)=rot (f)+rot g
(d) Sean ϕ,ψ:AR2R dos funciones ψ=ψ(x,y), ϕ=ϕ(x,y) tales que ϕx=ψy,  ϕy=ψx
Demuestre que ϕ, ψ son armónicas.

Solución. Para el inciso (a) se tiene
div (f+g)=(f+g)
=(x,y,z)(f+g)
=x(f+g)+y(f+g)+z(f+g)
=xf+xg+yf+yg+zf+zg
=xf+yf+zf+xg+yg+zg
=f+g

Para el inciso (a) se tiene
div (f+g)=(f+g)
=(x,y,z)(f+g)
=x(f+g)+y(f+g)+z(f+g)
=xf+xg+yf+yg+zf+zg
=xf+yf+zf+xg+yg+zg
=f+g

Para el inciso (c) se tiene
×(f+g)=(x,y,z)×(f1+g1,f2+g2,f3+g3)
=((f3+g3)y(f2+g2)z,(f3+g3)x(f1+g1)z,(f2+g2)x(f1+g1)y)
=[(f3yf2z)+(g3yg2z),(f3xf1z)+(g3xg1z),(f2xf1y)+(g2xg1y)]
=×f+×g

Para el inciso (d) se tiene
ϕx=ψy,    ϕy=ψx
por lo que
2ϕx2=x(ψy)
=2ψxy
por otro lado
2ϕy2=y(ψx)
=2ψyx
por lo tanto
2ϕx2+2ϕy2=2ψxy2ψyx=0
Analogamente se tiene
2ψx2=x(ϕy)
=2ϕxy
por otro lado
2ψy2=y(ϕx)
=2ϕyx
por lo tanto
2ψx2+2ψy2=2ϕxy+2ϕyx=0
como las funciones ψ, ϕ satisfacen la ecuación de Laplace entonces ambas funciones son armónicas.

Convergencia e integración

Por Lizbeth Fernández Villegas

Introducción

Así como ya hicimos comparaciones de continuidad o diferenciabilidad del límite de una sucesión de funciones a partir de sus términos, en esta ocasión lo haremos con funciones integrables.

Partimos de una sucesión de funciones donde para cada nN,fn:[a,b]R,a,bR. Supón además que (fn)nN converge puntualmente a una función f en [a,b].

Si cada una de las funciones fn son integrables, ¿será f también integrable?

¿La sucesión de integrales converge? ¿Su límite coincide con la integral del límite? Veamos el siguiente:

Ejemplo.

Considera el conjunto Q[0,1]. Como es numerable, podemos identificarlo como Q[0,1]={xn:nN}. Para cada nN definimos X{xn} como la función característica dada por:

X{xn}={1si x=xn0si xxn

Función X{xn}.

Ahora, para cada nN definimos fn=X{x1}++X{xn}.

Función fn.

Entonces la función fn es integrable en [a,b] y la sucesión (fn)nN converge puntualmente a la función:

XQ[0,1]={1si xQ[0,1]0si xQ[0,1]

Pero XQ[0,1] no es integrable en [0,1]. Por lo tanto la convergencia puntual podría no bastar para que el límite sea integrable. ¿Y si la convergencia es uniforme?

Proposición: Sea (fn)nN una sucesión de funciones integrables en [a,b] que converge uniformemente a una función f en [a,b]. Entonces f es integrable y
abf=limnabfn.

Demostración:
Para cada nN sea εn:=supaxb|fn(x)f(x)|.
Entonces fnεnffn+εn, de modo que las integrales superior e inferior de f satisfacen:
ab(fnεn)dxfdxfdxab(fn+εn)dx

Entonces 0fdxfdxab(fn+εn)dxab(fnεn)dx=2εn[ba].
Dado que εn0 porque (fn)nNf de manera uniforme, se sigue que f=f. Por lo tanto f es integrable.

Podemos ver también que
|abfdxabfndx|εn[ba]0 lo que demuestra que
abf=limnabfn.

Es importante mencionar que la convergencia uniforme no es una condición necesaria para que se de esta igualdad. Veamos el siguiente:

Ejemplo.

Para cada nN sea fn(x)=xn con x[0,1]. En la entrada Convergencia uniforme y continuidad mostramos que la sucesión (xn)nN converge puntualmente a la función:

f(x)={0 si 0x<11 si x=1

Si calculamos las integrales tenemos que para cada nN:

01xndx=1n+10=01f(x)dx.

Por lo tanto
01f(x)dx=limn01fn(x)dx.

Las condiciones de este ejemplo pueden generalizarse. Antes conozcamos algunas definiciones:

Definición. Sucesión uniformemente acotada: Sea (fn)nN una sucesión de funciones con fn:ARR,nN. Diremos que es uniformemente acotada en A si existe M>0R tal que |fn(x)|M para cualquier xA y cualquier nN.

Definición. Sucesión acotadamente convergente: Una sucesión de funciones (fn)nN con fn:ARR,nN es acotadamente convergente en A si converge puntualmente y es uniformemente acotada en A.

Proposición: Sea (fn)nN una sucesión acotadamente convergente en [a,b] donde cada función es integrable en [a,b], y que la función límite f es integrable en [a,b]. Supongamos también que existe una partición P de [a,b], a saber, P={x0,x1,,xm}, tal que la sucesión (fn)nN es uniformemente convergente hacia f en cada subintervalo [c,d][a,b] que no contenga ninguno de los puntos xkP. Entonces:

limnabfn(x)dx=abf(x)dx.

Demostración:
Dado que f es acotada y (fn)nN es uniformemente acotada, existe M>0 tal que |f(x)|M para cada x[a,b] y para cualquier nN.
Sea ε>0 tal que 2ε<P,
sea h=ε2m, donde m es el múmero de subintervalos de P,
considera una nueva partición P de [a,b] dada por:
P={x0,x0+h,x1h,x1+h,,xm1h,xm1+h,xmh,xm}

Nota que la función |ffn| es integrable en [a,b] y es acotada por 2M. Consideremos la integral de esta función en cada uno de los intervalos de la nueva partición P.

Por un lado, consideremos la suma de las integrales de |ffn| tomadas sobre los intervalos que sí tienen algún punto de P, es decir los intervalos
[x0,x0+h],[x1h,x1+h],,[xm1h,xm1+h],[xmh,xm].

La suma está dada por:

x0x0+h|ffn|(x)dx+x1hx1+h|ffn|(x)dx++xm1hxm1+h|ffn|(x)dx+xmhxm|ffn|(x)dx2M(x0+hx0)+2M(x1+h(x1h))++2M(xm1+h(xm1h))+2M(xm(xmh))=2Mh+2M(2h)++2M(2h)+2Mh=2M(2hm)=2Mε

El subconjunto restante de [a,b] lo llamaremos S. Está formado por un número finito de intervalos cerrados en los que (fn)nN converge uniformemente hacia f (pues no tiene ningún punto de P). Por consiguiente, existe NN tal que para cada xS, si nN se cumple que
|f(x)fn(x)|<ε

De modo que la suma de las integrales de |ffn| sobre los intervalos de S es a lo sumo ε(ba), luego para cada nN:

ab|f(x)fn(x)|dx(2M+ba)ε0.

Esto demuestra que abfn(x)dxabf(x)dx cuando n.

En la última sección de Análisis Matemático I hablaremos de la integral de Riemann-Stieltjes, que es un concepto que generaliza la integral de Riemann. La proposición vista aquí se puede expresar como sigue:

Proposición. Sucesión de funciones Riemann-Stieltjes: Sea α monótona en [a,b]. Supón que para cada nN,fnR(α) en [a,b]. Si (fn)nN converge uniformemente a f en [a,b] entonces fR(α) en [a,b] y:

abfdα=limnabfndα.

Más adelante…

Hablaremos de series de funciones y del límite de ellas. Así conoceremos el concepto de convergencia uniforme pero ahora en sumas infinitas.

Tarea moral

  1. Sea fn como en el primer ejemplo. Prueba que en efecto la sucesión (fn)nN no converge uniformemente a la función:
    XQ[0,1]={1si xQ[0,1]0si xQ[0,1]
  2. Sea (fn)nN una sucesión de funciones acotadas con fn:ARR,nN, tal que converge uniformemente a una función f:AR. Demuestra que (fn)nN es uniformemente acotada en A.
  3. Regresa luego de ver la integral de Riemann-Stieljes y demuestra la última proposición de esta sección.

Enlaces:

Contracciones

Por Lizbeth Fernández Villegas

Introducción

Cuando los puntos de un espacio métrico son enviados al mismo espacio a través de una función, conviene saber si habrá algún punto que se envíe a sí mismo, es decir, que se conserve fijo. Las próximas entradas nos mostrarán cuándo esa situación ocurre y resultados interesantes derivados de ello. Comencemos con la primera:

Definición. Función contracción: Sea (X,d) un espacio métrico y ϕ:XX una función. Diremos que ϕ es una contracción si existe α(0,1) tal que para cualesquiera x,yX se cumple que:
d(ϕ(x),ϕ(y))αd(x,y)

Podemos pensar entonces, que una función contracción, justamente hace que los puntos sean más cercanos entre sí de lo que eran originalmente.

Representación de una función contracción.

Nota que una contracción es también una función Lipschitz continua con constante de Lipschitz c<1. Este concepto se vio en la entrada Más conceptos de continuidad. Demos paso a otra:

Definición. Punto fijo: Sea X un espacio métrico y xX. Decimos que x es punto fijo de la función ϕ:XX si ϕ(x)=x.

Representación de un punto fijo.

Para ejemplificar estas ideas, veamos dos funciones que son contracciones y cómo existe un punto fijo en los casos a mencionar:

Ejemplos. f(x)=x2, con α=12.

Considera f:RR tal que f(x)=x2 en el espacio euclidiano. Sean x,yR. Sucede que:

d(f(x),f(y))=|f(x)f(y)|=|x2y2|=|12(xy)|=12|xy|=12d(x,y)

Por lo tanto d(f(x),f(y))12d(x,y) lo que demuestra que f es una contracción con α=12.

La siguiente imagen representa la diferencia de las distancias antes y después de aplicar la función en dos puntos x y y. Basta con observar las proyecciones de la gráfica de f en los ejes coordenados.

f(x)=x2

Ahora busquemos un punto fijo:

f(x)=xx2=xx=2x0=2xx0=x

Es decir, 0 es el único punto fijo de f.

A continuación, vamos a construir una sucesión de la siguiente manera:

Toma cualquier x0X

x1:=f(x0)=x02

x2:=f(x1)=x022=x022

x3:=f(x2)=x0222=x023
.
.
.
xk:=f(xk1)=x02k

Entonces la sucesión se define como (xn)nN donde xn=x02n.
Nota que tiende a 0 en R.

En las siguientes gráficas podemos observar el comportamiento de la sucesión:

Sea x0X. Mostramos la gráfica de la función f(x)=x2 y la función identidad I(x)=x.
Señalamos los términos x0 y x1=f(x0)=x02 y la distancia entre f(x0) y f(x1) vistos como proyecciones de las gráficas de los puntos sobre los ejes del plano cartesiano:

Términos x0 y x1

Si continuamos, generamos el punto x2=f(x1). Gráficamente también es visible que las distancias entre dos puntos disminuyen en el eje vertical al continuar con las iteraciones.


Términos x0,x1 y x2

Podemos observar que los puntos convergen a 0 que recordemos, es también el punto fijo de f.

Veamos otro caso:

Ejemplo. f(x)=x2+6, con α=12.

Considera f:RR tal que f(x)=x2+6 en el espacio euclidiano. Sean x,yR. Sucede que:

|f(x)f(y)|=|x2+6(y2+6)|=|12(xy)|=12|xy|

De modo que d(f(x),f(y))12d(x,y) lo cual prueba que f es una contracción con α=12.

Busquemos puntos fijos:

f(x)=xx2+6=x6=x212=x

Entonces 12 es el único punto fijo de f.

El siguiente gráfico nos confirma estos resultados para la sucesión generada a partir de un punto x0X donde para cada nN,xn=f(xn1).

Queda como ejercicio al lector demostrar que (xn)nN12 en R.

Esto da pie para enunciar el:

Teorema de punto fijo de Banach. Sea (X,d) un espacio métrico completo y sea ϕ:XX una contracción, entonces:

  1. Para cada x0X la sucesión (ϕn(x0))nN es de Cauchy y, en consecuencia (ϕn(x0))nN converge a un punto xX. ϕn representa la composición ϕϕnveces
  2. El punto x descrito es punto fijo de ϕ.
  3. El punto fijo es único.
  4. Podemos estimar la distancia de ϕn(x0) a x usando la desigualdad:
    d(ϕn(x0),x)αn1αd(x0,ϕ(x0)).

Por lo pronto demostremos que si una contracción tiene un punto fijo entonces este es único.

Sean x,yX tales que ϕ(x)=x y ϕ(y)=y. Como ϕ es una contracción se tiene que:

d(x,y)=d(ϕ(x),ϕ(y))αd(x,y)

Como α<1 se sigue que:
αd(x,y)d(x,y)
Por lo tanto d(x,y)=d(x,y), y en consecuencia x=y.

Más adelante…

Continuaremos con la demostración del teorema de punto fijo de Banach. En la siguiente entrada comprobaremos que la sucesión (ϕn(x0))nN es de Cauchy.

Tarea moral

  1. Sea f:RR tal que f(x)=x2+6 en el espacio euclidiano. Sea x0R, prueba que la sucesión (fn(x0))nN converge a 12.
  2. Sea f:[a,b][a,b],a,bR una función continua. Demuestra que tiene al menos un punto fijo.
  3. Da un ejemplo de una función continua f:[a,b][a,b],a,bR con una infinidad de puntos fijos.
  4. Prueba que si f:RR y para cada xR,|f(x)|M<1 entonces f es una contracción.
  5. Da un ejemplo de un espacio métrico completo y una función ϕ:XX que satisface que para todo xyX,d(ϕ(x),ϕ(y))<d(x,y) pero que no tenga ningún punto fijo.

Enlaces

Versión cuatro del Teorema de la Función Implícita

Por Angélica Amellali Mercado Aguilar

Ejemplo. Se da el nivel cero de una función diferenciable F:R4R y un punto P perteneciente a este nivel. Diga en cada caso si en los alrededores del punto p es posible ver la gráfica de F como la gráfica de una función diferenciable del tipo

a) u=u(x,y,z) b) z=z(x,y,u)c) y=y(x,u,z) d) x=x(y,z,u) para x2+y2+z2+uu=4 en p=(1,1,1,1)

Solución. En este caso para todos los incisos podemos definir f(x,y,z,u)=x2+y2+z2+uu4=0 y para el inciso a, se tiene
Fu=2u|(1,1,1,1)=20 por lo tanto es posible ver a la gráfica de F como una función diferenciable del tipo u=u(x,y,z) y sus derivadas parciales seran:

ux(1,1,1,1)=Fx|(1,1,1,1)Fu|(1,1,1,1)=2x2u=1
uy(1,1,1,1)=Fy|(1,1,1,1)Fu|(1,1,1,1)=2y2u=1
uz(1,1,1,1)=Fz|(1,1,1,1)Fu|(1,1,1,1)=2z2u=1

para el inciso b, se tiene

Fz=2z|(1,1,1,1)=20 por lo tanto es posible ver a la gráfica de F como una función diferenciable del tipo z=z(x,y,u) y sus derivadas parciales seran: zx(1,1,1,1)=Fx|(1,1,1,1)Fz|(1,1,1,1)=2x2z=1 zy(1,1,1,1)=Fy|(1,1,1,1)Fz|(1,1,1,1)=2y2z=1 zz(1,1,1,1)=Fz|(1,1,1,1)Fz|(1,1,1,1)=2u2z=1

para el inciso c, se tiene

Fy=2y|(1,1,1,1)=20 por lo tanto es posible ver a la gráfica de F como una función diferenciable del tipo y=y(x,z,u) y sus derivadas parciales seran: yx(1,1,1,1)=Fx|(1,1,1,1)Fy|(1,1,1,1)=2x2y=1 yz(1,1,1,1)=Fz|(1,1,1,1)Fy|(1,1,1,1)=2z2y=1 yu(1,1,1,1)=Fu|(1,1,1,1)Fy|(1,1,1,1)=2u2y=1 para el inciso d, se tiene Fx=2x|(1,1,1,1)=20 por lo tanto es posible ver a la gráfica de F como una función diferenciable del tipo x=x(y,z,u) y sus derivadas parciales seran:
xy(1,1,1,1)=Fy|(1,1,1,1)Fx|(1,1,1,1)=2x2y=1
xz(1,1,1,1)=Fz|(1,1,1,1)Fx|(1,1,1,1)=2z2y=1
xu(1,1,1,1)=Fu|(1,1,1,1)Fx|(1,1,1,1)=2u2y=1

Teorema de la Función Implicita (version (4))

Consideremos ahora el sistema
au+bvk1x=0 cu+dvk2y=0
con a,b,c,d,k1,k2 constantes. Nos preguntamos cuando podemos resolver el sistema para u y v en términos de x y y.
Si escribimos el sistema como
au+bv=k1x cu+dv=k2y
y sabemos que este sistema tiene solución si det|abcd|0 en tal caso escribimos
u=1det|abcd|(k1dxk2by),     v=1det|abcd|(k2ayk1cx).
Esta solución no cambiaria si consideramos
au+bv=f1(x,y) cu+dy=f2(x,y)

donde f1 y f2 son funciones dadas de x y y. La posibilidad de despejar las variables u y v en términos de x y y recae sobre los coeficientes de estas variables en las ecuaciones dadas.

Ahora si consideramos ecuaciones no lineales en u y v escribimos el sistema como
g1(u,v)=f1(x,y) g2(u,v)=f2(x,y)

nos preguntamos cuando del sistema podemos despejar a uy v en términos de x y y. Mas generalmente, consideramos el problema siguiente, dadas las funciones F y G de las variables u,v,x,y nos preguntamos cuando de las expresiones

F(x,y,u,v)=0 G(x,y,u,v)=0

podemos despejar a u y v en términos de x y y en caso de
ser posible diremos que las funciones u=φ1(x,y) y
v=φ2(x,y) son funciones implícitas dadas. Se espera que
n funciones u=φ1(x,y) y
v=φ2(x,y) en
F(x,y,φ1(x,y),φ2(x,y) G(x,y,φ1(x,y),φ2(x,y)
con (x,y) en alguna vecindad V. Suponiendo que existen φ1 y φ2 veamos sus derivadas

Fxxx+Fyyx+Fuux+Fvvx=0Fuux+Fvvx=Fx

Gxxx+Gyyx+Guux+Gvvx=0    Guux+Gvvx=Gx

Lo anterior se puede ver como un sistema de 2 ecuaciones con 2 incógnitas ux y vx. Aquí se ve que para que el sistema tenga solución

det|FuFvFuGv|0 en (P) (el det Jacobiano) y según la regla de Cramer

ux=det|FxFvGxGv|det|FuFvFuGv|,    vx=det|FuFxGuGx|det|FuFvFuGv| (con los dos det Jacobianos).

Análogamente si derivamos con respecto a y obtenemos
Fuuy+Fvvy=Fy ,Guuy+Gvvy=Gy
de donde

uy=det|FyFvGyGv|det|FuFvFuGv|,   vy=det|FuFyGuGy|det|FuFvFuGv| (con los dos det Jacobianos).

Al determinante det|FuFvGuGv| lo llamamos Jacobiano y lo denotamos por (F,G)(u,v).

Teorema de la Función Implícita (Versión 4)

Teorema 1. Considere las funciones z1=F(x,y,u,v) y z2=G(x,y,u,v). Sea P=(x,y,u,v)R4 un punto tal que F(P)=G(P)=0.
Suponga que en una bola BR4 de centro P las funciones F y G tienen (sus cuatro) derivadas parciales continuas. Si el Jacobiano (F,G)(u,v)(P)0 entonces las expresiones
F(x,y,u,v)=0 y G(x,y,u,v)=0 definen funciones (implícitas) u=φ1(x,y) y v=φ2(x,y) definidas en una vecindad v de (x,y) las cuales tienen derivadas parciales continuas en v que se pueden calcular como se menciona arriba.

Demostración. Dado que det|FuFvFuGv|0 entonces Fu(p), Fv(p), Gu(p), Gv(p) no son cero al mismo tiempo, podemos suponer sin pérdida de generalidad que Gv(p)0. Entonces la función z1=G(x,y,u,v) satisface las hipótesis del T.F.I y en una bola abierta con centro p, v se puede escribir como v=ψ(x,y,u). Hacemos ahora H(x,y,u)=F(x,y,u,ψ(x,y,u)) y tenemos que Hu=Fxxu+Fyyu+Fuuu+Fvψu=Fu+Fvψu por otro lado
ψu=GuGv por lo tanto Hu=Fu+Fvψu=Fu+Fv(GuGv)=FuGvFvGuGv0por lo tanto para H(x,y,u)=0 tenemos que existe una función u=φ1(x,y) y por lo tanto v=ψ(x,y,u)=ψ(x,y,φ1(x,y,u))=φ2(x,y) y por tanto u,v se pueden expresar en términos de x,y en una vecindad de p ◻

Ejemplo. Analizar la solubilidad del sistema
eu+ev=x+ye ueu+vev=xye

Solución. En este caso definimos
F(x,y,u,v)=eu+evxye=0 G(x,y,u,v)=ueu+vevxye=0
por lo que el sistema tendrá solución si det|FuFvFuGv|0
En este caso
det|FuFvFuGv|=det|euevueu+eeuvev+ev|=eu(vev+ev)ev(ueu+eu)=veu+vuev+u0
por lo tanto u y v se pueden ver en términos de x,y se pueden calcular sus parciales en u=0, v=1, x=1, y=1 que es este caso dan

ux=det|1yeevvev+ev|veu+vuev+u=(vev+ev)+evyeveu+vuev+u|(1,1,1,1)=2ee2e=2e vx=det|euueu+eu1ye|veu+vuev+u=yeue+ueu+euveu+vuev+u|(1,1,1,1)=e1e=1e1 uy=det|exeevvev+ev|veu+vuev+u=e(vev+ev)+evxeveu+vuev+u|(1,1,1,1)=e2+e2e2e=e vy=det|euueu+euexe|veu+vuev+u=euxe+e(ueu+eu)veu+vuev+u|(1,1,1,1)=eee=0