Ecuaciones Integrales - Daniela Araya
Ecuaciones Integrales - Daniela Araya
Ecuaciones Integrales - Daniela Araya
UNA INTRODUCCION
Trabajo de Graduacion presentado a la Facultad de Ciencias, en cumplimiento parcial de los requisitos exigidos para optar al grado de Licenciado
en Educacion Matematicas y Computacion.
Profesor Informante
Profesor Gua
Profesor Informante
Director
Agradecimientos
Al nalizar esta linda etapa de mi vida, miro hacia atras y veo a todas
las personas que me acompa~naron en este lindo y difcil sendero, personas
que por diferentes motivos no estan a mi lado y personas que siguen a mi
lado, que han tenido paciencia y el amor de ense~narme y soportar mi difcil
caracter, mi familia, mi novio Carlos y mi amigo Sebastian.
Agradezco tambien a todos los profesores que me hicieron clases, ya que cada uno de ellos me ense~no diferentes cosas, en especial al profesor Carlos
Lizama que con su paciencia y consejos nos ayudo a realizar nuestro trabajo.
Tambien agradezco a la familia de Sebastian por los ricos almuerzos y onces
que compartimos juntos. Por ultimo a Dios que me ha dado todo para poder
estar en esta etapa.
Indice general
1. Ecuaciones de Fredholm
1.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . .
48
5
2.2.1. Descomposicion de Nucleos . . . . . . . . . . . . . . . 54
2.3. Funciones Propias . . . . . . . . . . . . . . . . . . . . . . . . . 58
2.4. Teora en Lp . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
Introducci
on
La teora de Ecuaciones Integrales fue ampliamente desarrollada por el
matematico sueco Erik Ivar Fredholm quien desde temprana edad tuvo acceso a una de las mejores escuelas de Estocolmo, donde obtuvo con honores
su Bachillerato el 16 de Mayo de 1885. Despues de un a~no de estudiar en el
Royal Technological Institute ingresa a la Universidad de Uppsala donde obtiene el grado de Magister en Ciencias el 28 de Mayo de 1888. Para obtener el
grado de Doctor en Ciencias trabaja bajo la tutela del profesor Mittag-Leer
de la Universidad de Estocolmo, esto lo logra inscribiendo el Doctorado en
Uppsala, pero realizando sus estudios en Estocolmo. Recibe el grado de Doctor el 31 de Mayo de 1898. Su tesis doctoral abordo el estudio de Ecuaciones
Diferenciales Parciales que fue motivado por un problema de equilibrio en
elasticidad dentro del area de la Fsica. Esto dio origen a una de las investigaciones mas trascendentales para la vida de Fredholm y para la matematica
de la primera cuarta parte del siglo 20, a saber, Ecuaciones Integrales y
la Teora Espectral. Esta
fue utilizada por matematicos de renombre como
Schwarz, Neumann y Poincare. Cabe destacar que Hilbert complemento el
trabajo de Fredholm incluyendo la teora de valores propios para la Ecuacion
Integral de Fredholm.
El presente trabajo de titulacion comprende dos captulos, el primero
se titula Ecuaciones de Fredholm el cual presenta las ecuaciones generales
de Fredholm a resolver y analiza en profundidad la existencia de la(s) solucion(es) segun el tipo de nucleo que la compongan,estos son: nucleos separa-
7
Introduccion; en esta se presentan las Ecuaciones Integrales de primer y
segundo tipo, ademas se dene el Operador Integral involucrado en las ecuaciones anteriores, as como tambien se da un peque~no ejemplo particular que
ilustra las funciones que pertenecen al rango de este operador.
Caso Nucleos Separables; en esta seccion se estudia los nucleos separables
y de estos se desprende un caso particular: nucleos resolventes. Ademas, se
analizan las condiciones necesarias para la existencia de la(s) solucion(es) de
la ecuacion matricial analoga de la Ecuacion Integral de segundo tipo, esta
herramienta fundamental se encuentra en el Teorema 1.2.9. La seccion naliza con el estudio del Teorema Ecuacion de Fredholm y su contraparte para
la Ecuacion Integral de segundo tipo.
Caso Nucleos no Separables; en esta seccion se dene la funcion norma que
es una herramienta necesaria para estudiar un nuevo metodo para resolver la
Ecuacion Integral de segundo tipo, cuando esta se compone de un nucleo no
separable. Este metodo recibe el nombre de metodo de Schmidt. Se muestran
tres formas de lograr la descomposicion que propone. El metodo de Schmidt
requiere de la serie de Neumann para encontrar el inverso de un operador
(I
operador cK0 para que dicha serie converja. Asimismo, se estudian las condiciones para la existencia y unicidad de las soluciones solo cuando los nucleos
son peque~nos en el sentido de la norma que se considere.
El captulo dos, Teora de Fredholm en espacios de Lebesgue, esta dirigido
al estudio de las condiciones que debe cumplir el nucleo -el cual pertenece
a estos tipos de espacios- que compone la ecuacion Integral de Fredholm de
segundo tipo. Las secciones que componen este captulo son:
8
Funciones absolutamente Integrables; se dene el espacio de Lebesgue L1 y
las funciones que lo comprenden, ademas buscamos una cota adecuada para
Funciones de Cuadrado Integrable; se realiza un estudio analogo a la seccion anterior para el espacio de Lebesgue L2 . Por otra parte se analizan
caractersticas especiales de este espacio que permiten descomponer nucleos,
as como lo propone el metodo de Schmidt. Se presenta el Teorema de Fischer-
Cap
tulo 1
Ecuaciones de Fredholm
1.1.
Introducci
on
Denimos dos ecuaciones integrales que llamaremos ecuaciones de Fredholm de primer y segundo tipo respectivamente.
Z b
a
u(x) = f (x) +
Z b
a
x 2 [a; b]
(1.1)
x 2 [a; b]:
(1.2)
Z b
a
K0 u = f y u = f + K0 u:
9
(1.3)
10
Ejemplo 1.1.1.
Considerar el nucleo k (x; y ) = ex+y y las ecuaciones
Z
ex+y u(y ) dy = x
Z
u(x) = x +
(1.4)
(1.5)
Observamos que
Z
K0 u(x) =
Z
donde C =
x +y
u(y ) dy = e
ey u(y ) dy = Cex ;
u(x) = x +
entonces
u0 (x) = 1 + ex
Z
0
ey u(y ) dy = 1 + (u(x)
x)
u(x) = x + Cex ;
donde C es una constante arbitraria.
11
1.2.
N
ucleos Separables
1.2.1. N
ucleos Separables
Denicion 1.2.1. Decimos que un nucleo k(x; y) es separable si este se
puede escribir como
k (x; y ) =
n
X
i=1
yn
= xn
y
=
n
X
i=1
fi (x)gi (y ):
xn
x
yn
es separable. En efecto,
y
+ xn 2 y + xn 3 y 2 + xn 4 y 3 + ::: + xy n
+ yn
xn i y i 1 :
Supondremos que las funciones fi y gi en la denicion 1.2.1 son linealmente independientes. As el rango del operador K0 es n-dimensional y consiste en la combinacion lineal de funciones fi . Esta armacion se obtiene
como sigue:
(K0 u)(x) =
=
=
Z b
a
Z b
a
n
X
i=1
k (x; y )u(y ) dy
n
X
i=1
fi (x)
fi (x)gi (y ) u(y ) dy
Z b
a
gi (y )u(y ) dy:
12
Deniendo los coecientes ui como los productos internos
Z b
ui := hgi ; ui =
gi (y )u(y ) dy;
n
X
i=1
fi (x)ui :
(1.6)
(1.7)
existe una restriccion para la existencia de su solucion, esto es, h(x) debe
ser combinacion lineal de las funciones fi (x); ademas si esto ocurre, existen
innitas soluciones que satisfacen la ecuacion.
En el caso de la ecuacion de segundo tipo
u(x) = h(x) +
Z b
a
(1.8)
u(x) = h(x) +
j =1
fj (x)uj ;
(1.9)
donde solo resta encontrar los coecientes uj . Para esto denimos
hi := hgi ; hi =
Z b
a
gi (x)h(x) dx;
(1.10)
Kij := hgi ; fj i =
Z b
a
gi (x)u(x) dx =
Z b
a
gi (x)h(x) dx +
n Z b
X
j =1
gi (x)fj (x) dx uj
(1.11)
13
o usando la notacion (1.10) y (1.11):
ui = h i +
n
X
j =1
Kij uj :
(1.12)
Kij uj = hi
(1.13)
u(x) = 1 +
(1 + x + y + xy )
1=2
u(y ) dy:
k (x; y ) = (1 + x)
1=2
1=2
(1 + y )
y denimos
f1 (x) := (1 + x)
1=2
y g1 (y ) := (1 + y )
1=2
Luego
Z
h1 =
p1
1+x
dx = 2( 2
1) y K11 =
1
dx = ln 2:
1+x
14
As obtenemos
2( 2 1)
:
u1 = h1 + K11 u1 =
(1 2 ln 2)
As, por (1.9) una solucion es
2( 2
u(x) = 1 + p
x + 1(1
1)
:
2 ln 2)
u(x) = h(x) +
En este caso la integral
1=2
(xy )
u(y ) dy:
1
dx;
0 x
diverge y no podemos formar u(x) utilizando (1.12). Ahora bien, esto no
K11 =
u(y ) = y
h(x) = x
para
1=2
Z
0
py dy:
1.2.2. N
ucleo Resolvente
En esta seccion, reemplazaremos el operador K0 por operadores cK0 ,
donde c
cK )u = h:
Mu = h:
(1.14)
15
Denotaremos por D(c) al determinante de M .En lo que sigue supondremos
que D(c) es distinto de cero, as M sera invertible.
Sea C la matriz de los cofactores de M , ver [2], entonces tenemos que
C t M = MC t = ID(c):
Dado que existe M 1 , si denotamos M 1 := R(c) obtenemos que
R(c)
= C t =D(c);
uj =
n
X
i=1
Rji hi :
u(x) = h(x) +
= h(x) +
Ya que por (1.10) hi =
Z b
a
n
X
cfj (x)uj
j =1
"
n
X
j =1
cfj (x)
n
X
i=1
Rji hi :
gi (y )h(y ) dy , se obtiene
u(x) = h(x) +
= h(x) +
= h(x) +
n
X
"
cfj (x)
n
X
j =1
i=1
Z b
n
n
XX
a
j =1 i=1
"
Z b
n X
n
X
a
j =1 i=1
Rji
Z b
a
gi (y )h(y ) dy
16
n X
n
X
j =1 i=1
u(x) = h(x) +
Z b
a
(1.15)
u = f + cR(c):
(1.16)
u = (I + cR0 )h;
as el operador I + cR0 coincide justamente con el operador inverso de I cK0
ya que u = (I
cK0 ) 1 h.
Ejemplo 1.2.8. Para ilustrar como obtener un nucleo resolvente utilizaremos la siguiente ecuacion
Z
u(x) = h(x) + c
con nucleo separable k (x; y ) = e
x2
x2 y 2
K11 =
(1.17)
x2 e y 2 .
y a g1 (y ) = e
de la matriz K es
Z
u(y ) dy;
2x2
dx =
y2 .
1p p
2 ;
2
17
y as K =
1p p
2 . Luego
2
M = 1
2
=
y entonces
M
para todo c distinto de
1p p
2
2 p
p
c 2
;
2
p2 p ;
c 2
(1.17) es
pp
c 2
R(x; y ; c) =
x2
y2
:
2
Ahora bien, si h(x) = x entonces la solucion para la ecuacion (1.17) es de
acuerdo a (1.15) lo siguiente
Z
u(x) = x +
= x+c
= x+c
= x:
Pero >que ocurre si c =
2
2
pp
c 2
pp
c 2
pp
c 2
x2
x2
x2
e
Z
y2
y2
y dy
y dy
0
18
Derivando (1.17) con respecto a x, y suponiendo que h es una funcion derivable, se tiene que
u0 (x) = h0 (x)
= h0 (x)
x2
2xce
y2
1
h(x)):
2x(u(x)
u(y ) dy
h(x) obtenemos
g 0 (x) = 2xg (x);
g ( x) = e
= e
2t dt+k
x 2 +k
= Ae
x2
; (A = ek ):
u(x) = h(x) + Ae
x2
u(x) = x + Ae
x2
Ku = f;
19
sin saber si esta tiene alguna solucion o no, especulamos la posibilidad de
que pueda tener dos soluciones distintas u1 y u2 , entonces
Ku1 = f
Ku2 = f:
Restando las ecuaciones se obtiene:
K (u1
Denimos ' := u1
u2 ) = 0:
u = f + Ku;
obteniendo as la ecuacion K' = ', lo que equivale a (I K )' = 0. Deniendo A := I
K se obtiene
A' = 0;
(1.18)
20
ecuacion u = f + Ku:
K )u = f
u = (I
K ) 1 f:
X
i
fi (x)fi (y ):
Entonces la ecuacion
K0 ' = '
es igual a:
Z b
a
esto es:
Z b
(1.19)
21
o
XZ b
a
que es equivalente a
X
i
fi (x)
Z b
a
fi (y )'(y ) dy = '(x)
de donde obtenemos
X
i
fi (x)
Z b
a
fi (y )'(y ) dy
'(x) = 0
f i ( x)
Z b
a
fi (y )'(y ) dy + '(x) = 0:
esto es,
Z bX
a
o
X
i
que es equivalente a
fi (x)hfi ; 'i' dx +
Z b
a
'2 dx = 0
X
i
hfi; 'i
+ h'; 'i = 0:
Como cada termino del lado izquierdo es positivo, se tiene que estos terminos
deben ser iguales a 0, en particular h'; 'i = 0. Si ' es continua, entonces
22
Observacion 1.2.11. El producto interno entre vectores u y v 2 Rn sera denotado por el smbolo (u; v ); esto es
(u; v ) =
n
X
i=1
ui vi :
+ 2 xn
+ xn ) + xn (xn
+ 2x2n
+ 2x2n + 2xn 1 xn ) +
+ 2 x2 x 3 + 2 x4 x5 + 2 x6 x7 + : : : + 2 xn 2 xn 1 :
23
Veremos ahora que (x; Mx) = 0 solo cuando x = 0. En efecto, si (x; Mx) = 0
entonces
(x; Mx) = (2x21 + 2x22 + 2x1 x2 ) + (2x23 + 2x24 + 2x3 x4 ) + : : :
+ (2x2n
+ 2 x6 x7 + : : : + 2 xn 2 x n
=0
+ 2 xn 2 xn
+ x2n 1 ) + x2n = 0
+ xn )2
+ xn 1 )2 + x21 + x2n = 0:
M = 0;
es decir,
( ; Mi ) = 0 para cada i = 1; : : : ; n:
24
Como Mi son las las de la matriz M t , lo anterior se puede reescribir como
M t = 0:
En resumen se tiene:
(1.20)
de la ecuacion M t
u = u0 +
k
X
I =1
25
2 Ker(M t):
; fi =
Z b
(K0 u)(x) =
Z b
a
Z b
a
u(x)(K0 v )(x) dx =
Z b
a
(1.21)
Ahora bien;
Z b
a
u(x)(K0 v )(x) dx =
Z b
a
u(x)
Z b
a
k (x; y )v (y ) dy dx =
Z bZ b
a
u(x)(K0 v )(x) dx =
Z bZ b
a
26
As, de (1.21) tenemos
Z bZ b
a
Z b
a
k (y; x)u(y ) dy
(K0 u)(x) =
Z b
u(x) = f (x) + 2
(1.22)
'(x) = 2
(1.23)
= cx1=2
Z
donde c = 2
y 1=2 '(y ) dy: Esto muestra que existe una solucion distinta de
cero de la ecuacion (I
(x) =
(y )k (x; y ) dy:
27
En efecto:
Z
( x) =
k (x; y ) (y ) dy
k t (x; y ) (y ) dy
(yx)1=2 (y ) dy
k (x; y ) (y ) dy:
f (x) (x) dx = 0 o
f (x)x1=2 dx = 0:
Observacion 1.2.19. Si resolvemos directamente la ecuacion (1.22), podemos ver como se origina la condicion de Fredholm. En efecto, la ecuacion
(1.22) es equivalente a
u(x) f (x)
p
=
2 x
Ahora, derivando
1 1
u0 (x) = f 0 (x) + 2 p
2 x
esto es:
o
1
u0 (x) = f 0 (x) + p
x
u0 (x) = f 0 (x) +
y 1=2 u(y ) dy
u(x) f (x)
p
;
2 x
1
(u(x)
2x
f (x));
28
lo que es equivalente a:
u0 (x)
Sea g (x) = u(x)
1
(u(x)
2x
f 0 (x) =
g 0 ( x) =
de donde
px:
g ( x) =
Luego u(x) =
px + f (x) = f (x) + 2 Z
esto es
(xy )1=2 ( y + f (y )) dy
o
f (x)):
y dy +
lo que es igual a
o
Z
p
pyf (y) dy:
0=2 x
1
Z
0
pyf (y) dy = 0
o equivalentemente
Z
pyf (y) dy = Z
(y )f (y ) dy = ( ; f ) = 0;
29
1.3.
N
ucleos No Separables
1.3.1. Normas
Denicion 1.3.1. Si en un espacio vectorial E se puede denir una funcion
kvk > 0 si v 6= 0
2.
kcvk = jcjkvk
3.
ku + vk kuk + kvk,
para todo u; v
funcion
2E
yc
2 K, entonces el espacio E
k k la llamamos norma.
se dice normado y la
kKvk B (K )kvk
30
hecho, se puede demostrar que lo es cuando se encuentra una representacion
de este numero) y la denotamos por kK k. Entonces:
kKvk kK kkvk;
con la igualdad para algun v 6= 0.
kK k = miax
X
j
jKij j:
j(Kv)ij =
X
Kij vj
j
X
j
jKij jjvj j:
Como kv k es mas grande que cada uno de los valores jvj j, obtenemos
Ri (K ) =
X
j
jKij j:
sgn(x) :=
8
>
>
>
<
>
>
>
:
1;
x > 0;
0;
x = 0;
1; x < 0.
31
Tomando vj = sgn(Kmj ) se tiene kv k = 1 y
(Kv )m =
X
j
jKmj j
kK k = miax
Observacion 1.3.3. Para cada n
X
j
jKij j:
2 N, kK nk kK kn.
En efecto: Para
kKK k kK kkK k:
De lo anterior obtenemos, por induccion
kK nk kK kn:
Observar que en general kK n k es estrictamente mas peque~na que kK kn
como veremos en el ejemplo (1.3.14).
32
Demostracion. Es claro que kuk > 0 y que kcv k = jcjkv k. Ademas como
sup ju(x) + v (x)j
kK k = sup
0
Z b
jk(x; y)j dy
(1.24)
u = h + S0 u + E0 u
esto es,
(I
E0 )u = h + S0 u:
E0 ) 1 .
Entonces se tiene:
u = h + S0 u
(1.25)
33
donde h = (I
E0 ) 1 h, S0 = (I
E0 ) 1 S0 .
s(x; y ) =
Luego
n
X
i=1
fi (x)gi (y ):
X
s (x; y ) =
fi (x)gi (y );
donde f = (I
i
i=1
E ) fi .
1
Por lo tanto, para resolver la ecuacion integral (1.2) inicial, solo tenemos
que resolver la ecuacion separable (1.25).
E0 ) puede
ser calculada formalmente por medio de series. Estas series son las llamadas
series de Neumann :
(I
E0 )
1 m
;
n
n
con m = 1; : : : ; n:
34
Denimos
1
2
km (x) = k x; m
Consideremos la funcion
um (y ) =
8
<
:
1; y 2
0; y 2
=
m 1
n ;
m 1
n ;
1
:
n
m ;
n
m .
n
y , por
k (x; y ) =
n
X
m=1
km (x)um (y ) + e(x; y )
u(x) = h(x) +
= h(x) +
n
X
m=1
n
X
m=1
km (x)um (y )u(y ) dy
km (x)
Z (m=n)
(m 1)=n
u(y ) dy:
1
km (x) =
2
Z
k (x; y )e
imy
dy:
k (x; y ) =
n
X
m= n
km (x)um (y ) + e(x; y )
! k(x; y)
35
Notemos que la aproximacion obtenida es diferenciable en y y se tiene que
la ecuacion de segundo tipo se puede escribir como
Z
u(x) = h(x) +
= h(x) +
n
X
kk (x)eimy u(y ) dy
0
n
Z
n
X
1
2
k (x; s)e
ims
Z
ds
p(x; y ; ") =
X
n
kn (x; ")y n ;
u(x) = h(x) +
= h(x) +
X
n
p(x; y ; ")u(y ) dy
Z
kn (x; ")
y n u(y ) dy:
1.3.4. Iteraci
on por medio de Series de Neumann
Sea c
36
a traves de la formula
un+1 = f + cKun ;
(1.26)
un = f + cKun
= f + cK (f + cKun 2 ) = f + cKf + c2 k 2 un
un = f + cKf + c2 K 2 f + c3 K 3 f + ::: + cn 1 K n 1 f + cn K n u0 :
Lo anterior sugiere que:
u=f+
1
X
n=1
cn K n f:
(1.27)
1.
Demostracion. Sea SN =
kf k < 1 se tiene
kSM
donde kf k
SN k =
M
X
n=N
1
X
n=0
N
X1
n=0
M
X
n n
c K f
n =N
M
X
n=N
jcj kK kkf k
n
jcjnkK kn ! 0, cuando N
M
X
n=N
jcjnkK knkf k;
n=0
37
cK ) 1 .
1
X
n=0
cn K n coincide con
cK )
1
X
n=0
cn K n =
Demostracion. En efecto
1
X
(I cK )
cn K n =
n=0
1
X
1
X
n=0
cn K n (I
cn K n (I
n=0
c K
n=0
= I+
Analogamente
1
X
1
X
n=1
1
X
n=0
cn K n
cK ) = I:
cn+1 K n+1
1
X
n=0
cn+1 K n+1 = I:
cK ) = I:
u = f+
1
X
n=1
cn K n f = I +
1
X
n=1
cn K n f
Por otra parte, de (1.16) se tiene que u = (I + cR(c))f , luego cR(c) coincide
1
X
exactamente con la serie
cn K n .
n=1
38
jcj kK k < 1:
2
1
X
n=1
jcj kK k kf k
n
1
X
n=1
Proposicion 1.3.13. Para que la serie de Neumann sea convergente es suciente que
K=@
0
2
1
3
1
A;
K2 = @
6 11
1
A
jcjkK k =
2 1 2
< 1:
39
1 1
K=@
0 1
1
A
1 n
Kn = @
0 1
1
A
(I
cK )
B
B
B
@
1+
1
X
n=1
1
X
nc
n=1
1
X
cn
1+
1
C
C
C;
A
n=1
(I
cK )
B
=@ 1
c (1
c
1
c)2
C
A:
jcj < 1=
40
kK k
1 = max
kk
X
i
jKij j;
es norma.
Ahora comparemos las cotas encontradas a traves de suponer que jcjkK k1 < 1
y jcj2 kK 2 k1 < 1. Dada la matriz
K=@
0
2
1
3
1
A;
suciente que jcjkK k < 1,(ver Proposicion 1.3.8) en tal caso no se ocuparon
las propiedades particulares de la norma asociada a una matriz.
41
k um
n
X
j =0
cj K j f satisface
un k < ";
para cada m > n: Ya que el espacio vectorial de las funciones con la norma denida es completo (ver[4]), entonces la serie de Cauchy es una serie
convergente, es decir, existe una funcion u tal que
kun
uk ! 0; cuando n ! 1
n
X
2 N, entonces la suma
k=1
Z
0
zamos dos casos en que la norma del operador asociado a K0 es menor que
1 y, as la serie de Neumann converge uniformemente a la solucion u.
1. k (x; y ) = sin(xy ): Como 0 < x < 1 entonces xy < y de donde sin(xy ) <
sin(y ). De esto se tiene que:
kK k = sup
0
Z
0
j sin(xy)j dy
Z
0
sin y dy = 1
cos(1) < 1:
42
2. k (x; y ) = 2jx
kK k
= sup
2jx
Z
sup 2
x
y j dy
Z
jxj dy +
jyj dy
= sup 2x + 1 3:
x
kK k
2
0
= sup
x
= sup 4
x
0
1
jx
4jx
z jjz
zj
Z z
0
y j dz dy
jz
y j dy +
Z
z
jz
y j dy dz;
donde
Z z
0
jz yj dy +
jz yj dy =
Z z
0
z y dy +
1
z + y dy = z + + z 2 ;
2
por lo tanto:
kK k
2
0
=
=
=
=
jx zj z + 21 + z2
sup 4
x
0
Z z
Z x
1
1
2
2
z + + z dz
(x z )
z + + z dz + (z x)
2
2
x
0
2
3
4
2
3
4
3x 2x + x
1 2x + 3 x 2x + x
sup 4
+
2
6
x
2
3
4
1 2x + 3 x 2x + x
2
sup 4
= :
6
3
x
1
u(x) = f (x) + c
ex y u(y ) dy
43
encontraremos los valores de c para que la serie de Neumann converja. Para
esto resolvemos la ecuacion directamente.
Derivando con respecto a x obtenemos:
Z
0
f (x):
1 = e f (x)
1
X
n=1
cn ;
por tanto para que esta serie converja se debe tener que jcj < 1.
Ahora bien aplicando la denicion de norma sobre el operador integral K0
se tiene que:
kK k = sup
0
je j dy = sup e [
x y
e ] = sup e
1
0
1
e
1
0; 3679.
e
Al comparar las dos formas que hemos utilizado para encontrar el valor de c,
As kK0 k < e. Suponiendo jcjkK0 k < 1 se tiene que jcj <
observamos que no hay una diferencia signicativa para este valor, es decir,
se conserva la condicion de que jcj < 1, para que la ecuacion inicial posea
solucion.
Observacion 1.3.21. .
1. Si f es una funcion continua y K0 u es continua para cualquier u entonces la solucion u es continua.
44
2. Si K0 es un operador continuo y f es un funcion discontinua entonces
que f .
K0 )
u = (I
K0 ) 1 f
kuk k(I
K0 )
kkf k < 1:
45
' 6= 0:
K' = ';
kK k < 1 obtenemos:
0
donde se tiene una contradiccion que nace del supuesto u1 6= u2 ; as se tiene
que ' = 0, por tanto u1 = u2 . As, la ecuacion u = f + Ku tiene una unica
solucion.
Veamos el siguiente ejemplo que ilustra la importancia de que la norma
de ' sea nita.
Dado el nucleo
k (x; y ) :=
8
<
:
yx y ; 0 y x 1
0x<y1
0;
Calculamos la norma de K0 :
kK k
0
= sup
x
= sup
x
Dado que x
por:
Z x
0
Z x
yx
jk(x; y)j dy +
Z
x
jk(x; y)j dy
dy:
kK k sup
0
Z x
0
y dy =
1
< 1:
2
46
As, kK0 k < 1. Sea '(x) = xx
entonces su norma es
Anteriormente vimos que si k (x; y ) = s(x; y )+ e(x; y ) donde s(x; y ) es separable y e(x; y ) es peque~no, entonces la ecuacion u = f + K0 u puede ser reducida
a un sistema de ecuaciones lineales algebraicas. De esto se desprende que el
espacio de las funciones con norma nita kuk <
C0 (u)(x) =
Z b
a
Z b
a
kC k = sup
0
Z b
a
Z b
a
47
Luego:
kK k = kC
0
kE k.
1
Observacion 1.3.25. Existe la posibilidad de que e1 posea innitas discontinuidades en su dominio, pero aun as kE1 k < 1.
k (x; y ) = jx
c(x; y ) = (jx
yj
1
2
y j + ")
1
2
Cap
tulo 2
Teor
a de Fredholm en espacios
de Lebesgue
48
49
2.1.
kuk
Z b
a
ju(x)j dx < 1:
(2.1)
Observacion 2.1.2. Se verica que kuk1 = 0 s y solo s u(x) = 0 para cada
kK k sup
0
Z b
a
50
Demostracion. En efecto,
kK uk
0
Z b
j(Ku)(x)jdx =
a
Z bZ b
a a
Z b
a
= sup
y
jk(x; y)jju(y)j dy dx =
ju(y)j sup
Z b
a
ess sup
Z b
a
ju(y)j
jk(x; y)j dx dy
jk(x; y)j dx dy
Observaci
on 2.1.5. Si y
Z
b
Z b Z b
dx
k
(
x;
y
)
u
(
y
)
dy
a
a
Z b
Z b
jk(x; y)j dx, donde ess sup denota supremo esencial (ver [1]).
A n de probar que esta solucion pertenece al espacio L1 , formamos la sucesion un de la siguiente forma:
un =
n
X
i=0
ci K0i f:
Ejemplo 2.1.6. Para ilustrar como esta nueva norma ampla los tipos de
funciones f y k que podemos ocupar, analicemos la siguiente ecuacion
u(x) = x
1=2
(x
1)
1=3
u(y ) dy:
(2.2)
51
En primer lugar observamos que
kf k =
sup jx
x2[0;1]
1=2
j=
1
sup p ;
x2[0;1] x
kK k =
0
sup j(1
x2[0;1]
x)
1=3
j=
1
sup p
;
3
x
x2[0;1] 1
kf k
kK k
0
jx
1=2
sup
y2[0;1]
j dx = 2;
j(1
x)
1=3
j dx = 32 :
Ya que ambas normas son nitas concluimos que la ecuacion (2.2) tiene
solucion y esta pertenece a L1 (0; 1).
2.2.
kuk
Z b
a
u (x) dx < 1:
2
(2.3)
kf k =
Z b
a
jf (x)jdx
1=2
52
Proposicion 2.2.4.
kK k
Z b Z b
jk(x; y)j
1=2
dydx
kK uk
0
2
2
Z b
a
Z b
a
(K0 u) (x) dx =
2
ahora si B (K0 ) =
2
2
2
Z b
a
Z b Z b
a
dx = kuk
2
2
kk(x; )k
2
2
2
dx =
k (x; y )u(y ) dy
2
2
Z b
a
kk(x; )k
2
2
Z b
a
hk(x; ); ui dx
dx;
dx entonces
kK uk B (K )kuk :
0
Observacion 2.2.5. Notemos que B (K0 ) < 1 implica kK0 k < 1, as la
ecuacion u = f + K0 u tiene solucion. Esta solucion se encuentra en L2 [a; b],
puesto que este espacio es completo, ver [4].
Los siguientes teoremas complementan la demostracion del Teorema FischerRiesz que es una herramienta util para construir funciones de cuadrado integrable a partir de vectores de dimension innita, el cual nos permite llevar
la teora de Fredholm a este tipo de espacios vectoriales.
53
Teorema 2.2.6 (Beppo Levi). Sea (fn (x)) una sucesion de funciones, tal
que
f (x) es integrable, y
Z
f (x) dx =
XZ
fn (x) dx:
fn (x)
L2 tal que
1
X
n=1
g=
1
X
n=1
fn un ;
1
X
i=0
fi ui (x) y L
kSm
S n k L1 =
Z X
m
f
u
(
x
)
i i
i=n+1
1=2
= L
= L1=2
donde
n+1
X
i=1
fi2
kSm
m
X
dx
S n k L2 = L
Z
1=2
m
X
i=n+1 j =n+1
12 dx
m
X
i=n+1
1=2
2
Z X
m
@
fi ui (x)
i=n+1
m
X
fi ui (x);
fi fj hui ; uj i = L1=2
n+1
X
i=1
11=2
j =n+1
dxA
fj uj (x)i
fi2 ;
54
convergente con lmite g en L1 . De esto se tiene que Sn es una suma parcial
acotada, con cota k . Ahora, si hacemos
Z
Sn (x) dx =
2
y como Sn2
n
X
i=1
!2
fi ui
k 2 dx < 1
dx <
2.2.1. Descomposicion de N
ucleos
En el metodo de Schmidt se supone que un nucleo k (x; y ) se puede descomponer como
x; y 2 [a; b];
k (x; y ) =
jnj<N
kn (x)einy +
jnjN
kn (x)einy ;
1
donde kn (x) =
k (x; y )e iny dy son los coecientes de Fourier de la
2
funcion y ! k (x; y ) con x jo. (ver [1]).
El uso de los coecientes de Fourier tiene una motivacion especial. Si aproximamos una funcion f , primero por medio de una combinacion lineal arbitraria y luego utilizando coecientes de Fourier, obtenemos respectivamente:
f ( x) =
N
X
n=1
cn un (x) + E (x);
f (x) =
N
X
n=1
fn un (x) + e(x);
55
donde un es cualquier sistema ortonormal completo. Si hacemos la diferencia
entre los terminos anteriores, encontramos
E (x) =
N
X
n=1
( fn
kE k
2
2
=k
N
X
n=1
2
2
=
=
=
N
X
N
X
n=1
n=1
fn un ya que
hf
fm um ; (fn cn )un i = hf
fm um ;
fn un
c n un i
X
X
X
X
X
X
hf; fnuni hf; cnuni h fmum; fnuni + h fmum; cnuni
XX
X
X
XX
fm cn hum ; un i
fm fn hum ; un i +
fn hf; un i
cn hf; un i
X
fn hf; un i
m n
X
fn2
n
cn hf; un i
fn cn = 0:
kE k
2
2
= kek +
2
2
N
X
n=1
(fn
cn )2 :
As, concluimos que el error se minimiza cuando fn = cn . Por lo tanto el utilizar coecientes de Fourier para aproximar una funcion en L2 [a; b] optimiza
el tama~no del error.
X
n2Z
hf; uniun:
56
El smbolo =2 representa una igualdad con la norma del espacio L2 [a; b],
L
es decir,
n
X
n
f; ui ui
! 0;
n ! 1:
(2.4)
kf k
2 L [a; b] se tiene
2
X
n2Z
jfij ;
(2.5)
n
X
n
2
f; ui ui
= hf
n
X
n
X
hf; uiiui; f
hf; uiiuii
que es equivalente a
hf; f i
o
2hf; f
kf k
n
X
n
X
hf; uiiuii + hf
n
X
n
fi +
2
n
X
n
hf; uiiui; f
fi = k f k
2
n
X
n
hf; gi =
X
n2Z
fi gi ;
kf + g k
= kf k2 + 2hf; g i + kg k2 :
n
X
n
fi2 :
hf; uj iuj i
57
Por otro lado ocupando la identidad (2.5) se tiene que
kf + gk
n2Z
Consideremos y
(fi + gi )2 =
hf; gi =
fi2 + 2
n2Z
X
n2Z
X
n2Z
fi gi +
X
n2Z
gi2
fi gi :
k (x; y ) =2
L
X
n2Z
kn (x)un (y );
donde kn (x) son los coecientes de Fourier de k (x; y ). Estos son:
Z b
a
k (x; y ) =
Deniendo sn (x; y ) =
mos
jnj<N
jnj<N
kn (x)un (y ) +
jnj>N
kn (x)un (y ):
kn (x)un (y ) y en (x; y ) =
jnj>N
(2.6)
kn (x)un (y ) obtene-
Z b
a
Armacion: en ! 0 cuando n ! 1:
En efecto, como k (x; y ) =
X
n2Z
kn (x)un (y ), entonces
en (x; y ) = k (x; y )
jnjN
kn (x)un (y )
58
de donde que
n!1
sn k = 0;
kEnk
Z bZ b
a
en (x; y )2 dx dy;
2.3.
Funciones Propias
Denicion 2.3.1. Diremos que una sucesion (an )n2N en un espacio E con
producto interno <; > converge debilmente a un lmite a
cribiremos an ! a, si
w
2 E , lo cual es-
cuando n ! 1:
E
59
Teorema 2.3.4. Sea K0 el operador integral de la ecuacion (1.2) y supongamos que K0 es simetrico, entonces K0 posee una funcion propia.
Por denicion de supremo se tiene que existe una sucesion (un ) A tal que
1
n
jkj
jQ(un)j;
jkj
1
n
jQ(un)j jkj:
sucesion, (unk ), de la sucesion (un ) tal que unk ! ' para algun ' 2 A.
w
jQ(unk )j !w jQ(')j; n ! 1:
As, por unicidad de lmite se tiene que:
jQ(')j = jkj:
60
Por tanto, jQ(u)j jQ(')j para cada u 2 A.
Dado y
Q
pero
Q
y
ky k2
y
kyk2
y
y
kyk2 ; K0 kyk2
y
2
ky k2
. Entonces u
2 A y por la parte
jQ(')j
1
=
h
y; K0 y i = jQ(y )j
2
kyk2
De esto se tiene que jQ(y )j jQ(')j para todo y 2 L2 [a; b], y as Q posee un
maximo en ', lo cual prueba la armacion.
q (t) = Q(' + ty ) =
=
2
2
2
2
2
2
q (t) =
2
2
2
2
61
Derivando la funcion q obtenemos
q 0 (t) =
(2hK0 '; y i + 2thy; K0 y i)k' + ty k22 h' + ty; K0 (' + ty )i(2h'; y i + 2tky k22 )
(k'k22 + 2th'; y i + t2 ky k22 )2
2h'; K0 'ih'; y i
= 2hK0 '; y i
2h'h'; K0 'i; y i;
concluimos que
hK '
0
para todo y
'h'; K0 'i; y i = 0;
K0 ' = 'h'; K0 'i donde ' es la funcion propia que buscabamos y h'; K0 'i
es su valor propio.
2.4.
Teor
a en
Lp
kukp =
Z b
a
ju(x)j
dx
1=p
< 1;
p 1:
62
Se denota por Lp [a; b] al espacio vectorial de todas las funciones que
cumplan con la Denicion 2.4.1, donde u = v s y solo s u(x) = v (x) para
cada x que no este en un conjunto de medida nula.
ju(x)v(x)j dx
Z b
a
ju(x)j
dx
1 1
+ = 1 con p; q
p q
1, entonces
1=p Z b
1=q
jv(x)j
dx
u(x)
1 dx
=
Z b
n=m Z b
u(x) dx
a
" Z
b
a
u(x)m dx
a
1=m #n
<
(b
(m
n)=m
dx
a)m=(m
m=(m n)
n)
ju(x) + v(x)j
dx
1=p
Z b
a
ju(x)j
1=p Z b
dx
jv(x)j
dx
1=p
63
Las demostraciones de las Proposiciones 2.4.2 y 2.4.5 se encuentran en
[6].
kk
por
k kp
Proposicion 2.4.7.
kK ukp
Z b
kk(x; )k
p
p dx
1=p
Demostracion. En efecto,
kK ukp
=
Z b
j(K u)(x)j
dx
1=p
a
p
Z b Z b
k (x; y )u(y ) dy
a
a
Z b Z b
a
Z b
a
jk(x; y)u(y)j dy
hk(x; ); ui
dx
1=p
!1=p
dx
p
!1=p
dx
kK ukp
0
Z b
a
kukq
Z b
a
p
q dx
kk(x; )k
1=p
p
p dx
1=p
k kp es una
64
Demostracion. En efecto,
Z b
a
(u(x)v (x)) dx
1=p
Z b
u(x) v (x) dx
1=p
1=p
u(x) v (x) dx
kupkpkvpkq :
(u(x)v (x)) dx
1=p
kukppkvkpq < 1:
Observacion 2.4.9. .
1. De la proposicion anterior se obtiene que el operador K0 u pertenece a
1=p
p
; observamos que para
k (x; ) p dx
a
con f (x) en Lp [a; b] y k (x; y ) en Lq [a; b] solo
Z b
k
p!1
p!1
Bibliograf
a
65