Un problema resonante para un sistema - Universidad de Buenos ...
Un problema resonante para un sistema - Universidad de Buenos ...
Un problema resonante para un sistema - Universidad de Buenos ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>Un</strong> <strong>problema</strong> <strong>resonante</strong><br />
<strong>para</strong> <strong>un</strong> <strong>sistema</strong><br />
con condiciones periódicas<br />
Pablo Amster - Pablo De Nápoli<br />
p<strong>de</strong>napo@dm.uba.ar<br />
<strong>Un</strong>iversidad <strong>de</strong> <strong>Buenos</strong> Aires y CONICET
Introducción<br />
Estudiamos el <strong>sistema</strong> <strong>de</strong> ecuaciones no lineales <strong>de</strong><br />
seg<strong>un</strong>do or<strong>de</strong>n <strong>para</strong> <strong>un</strong>a f<strong>un</strong>ción vectorial<br />
u : [0, 2π] → R N que satisface<br />
u ′′ + m 2 u + g(u) = p(t) 0 < t < 2π (1)<br />
con condiciones <strong>de</strong> frontera periódicas:<br />
u(0) = u(2π), u ′ (0) = u ′ (2π). (2)<br />
Aquí m = 0 es <strong>un</strong> entero, p ∈ L 2 (0, 2π), y la no<br />
linealidad g es continua y acotada. Entonces, (1-2) es<br />
<strong>un</strong> <strong>problema</strong> <strong>resonante</strong>, ya que el núcleo <strong>de</strong>l operador<br />
Lmu := u ′′ + m 2 u actuando en el espacio <strong>de</strong><br />
f<strong>un</strong>ciones 2π-periodicas es no-trivial.
Esta situación se conoce en la literatura como<br />
resonancia en <strong>un</strong> a<strong>un</strong>tovalor <strong>de</strong> or<strong>de</strong>n superior: <strong>de</strong><br />
hecho, si consi<strong>de</strong>ramos el <strong>problema</strong> <strong>de</strong> autovalores<br />
−u ′′ = λu con condiciones períodicas, los autovalores<br />
son λm = m 2 ∈ N0.<br />
Recor<strong>de</strong>mos que en el caso m = 0 <strong>para</strong> <strong>un</strong>a ecuación<br />
escalar, existe <strong>un</strong>a solución si se asumen las<br />
condiciones <strong>de</strong> Lan<strong>de</strong>sman-Lazer conditions, que<br />
fueron introducidas en el trabajo:<br />
E. Lan<strong>de</strong>sman and A. Lazer: Nonlinear perturbations<br />
of linear elliptic bo<strong>un</strong>dary value problems at<br />
resonance, J. Math. Mech. 19 (1970), 609-623.<br />
<strong>para</strong> <strong>un</strong>a ecuación elítpica <strong>resonante</strong> escalar <strong>de</strong><br />
seg<strong>un</strong>do or<strong>de</strong>n, bajo condiciones <strong>de</strong> Dirichlet.
<strong>Un</strong> survey sobre condiciones <strong>de</strong> Lan<strong>de</strong>sman-Lazer:<br />
J. Mawhin: Lan<strong>de</strong>sman-Lazer conditions for<br />
bo<strong>un</strong>dary value problems: A nonlinear version of<br />
resonance. Bol. <strong>de</strong> la Sociedad Española <strong>de</strong><br />
Mat.Aplicada 16 (2000), 45-65.<br />
Grosso modo, las condiciones anteriores dicen que si<br />
p (el promedio <strong>de</strong> p) se encuentra entre los límites en<br />
±∞ <strong>de</strong> la nolinealidad g, entonces el <strong>problema</strong><br />
admite por lo menos <strong>un</strong>a solución. Esta condición<br />
pue<strong>de</strong> consi<strong>de</strong>rarse como <strong>un</strong>a condición <strong>de</strong> grado<br />
sobre la esfera S0 = {−1, 1}, en el siguiente sentido:<br />
<strong>para</strong> v = ±1, <strong>de</strong>finimos g±1 = g(±∞), entonces la<br />
f<strong>un</strong>ción θ : S0 → S0 dada por θ(v) = gv−p<br />
|gv−p| cambia <strong>de</strong><br />
signo, y en consecuencia, tiene grado no nulo.
De este modo, el siguiente resultado, que es <strong>un</strong>a<br />
adaptación <strong>de</strong> <strong>un</strong> teorema <strong>de</strong> Nirenberg <strong>para</strong> <strong>sistema</strong>s<br />
elípticos, que aparece en el artículo:<br />
L. Nirenberg: Generalized <strong>de</strong>gree and nonlinear<br />
problems, Contributions to nonlinear f<strong>un</strong>ctional<br />
analysis, Ed. E. H. Zarantonello, Aca<strong>de</strong>mic Press<br />
New York (1971), 1-9.<br />
pue<strong>de</strong> ser consi<strong>de</strong>rado como <strong>un</strong>a extensión natural <strong>de</strong>l<br />
teorema <strong>de</strong> Lan<strong>de</strong>sman y Lazer:
Teorema (Niremberg)<br />
Supongamos que existen los límites radiales<br />
gv := límr→+∞ g(rv) exist <strong>un</strong>iformly respecto a<br />
v ∈ S N−1 , la esfera <strong>un</strong>idad <strong>de</strong> R N . Entonces (1-2) con<br />
m = 0 posee al menos <strong>un</strong>a solución T -periodica, si se<br />
satisfacen las siguientes condiciones:<br />
• gv = p := 1<br />
T<br />
T<br />
0 p(t)dt <strong>para</strong> todo v ∈ SN−1 .<br />
• El grado <strong>de</strong> la aplicación θ : S N−1 → S N−1 dada<br />
por<br />
θ(v) = gv − p<br />
|gv − p|<br />
es diferente <strong>de</strong> 0.
Observamos que el promedio p pue<strong>de</strong> ser visto como<br />
la proyección ortogonal <strong>de</strong>l término forzante p sobre<br />
el núcleo <strong>de</strong>l operador lineal L0, que es el conj<strong>un</strong>to <strong>de</strong><br />
las f<strong>un</strong>ciones constantes, que se i<strong>de</strong>ntifica<br />
naturalmente con R N .
En contraste con el caso anterior, si m = 0 el núcleo<br />
es 2N-dimensional, más precisamente:<br />
Ker(Lm) = cos(mt)α + sin(mt)β : (α,β) ∈ R 2N .<br />
Entonces, <strong>un</strong>o pue<strong>de</strong> esperar que <strong>un</strong>a condición <strong>de</strong>l<br />
tipo <strong>de</strong> Lan<strong>de</strong>sman-Lazer correspondiente a este caso,<br />
<strong>de</strong>be po<strong>de</strong>r expresarse en términos <strong>de</strong> la proyección<br />
ortogonal <strong>de</strong> p sobre Vm o, equivalentemente, en<br />
términos <strong>de</strong> los m-ésimos coeficientes <strong>de</strong> Fourier <strong>de</strong><br />
p. Para N = 1, esto fue <strong>de</strong>mostrado en el artículo:<br />
A. Lazer, D. Leach, Bo<strong>un</strong><strong>de</strong>d perturbations of forced<br />
harmonic oscillators at resonance, Ann. Mat. Pura<br />
Appl. 82 (1969), pp. 49-68.
Teorema (Lazer-Leach)<br />
Sea N = 1 y supongamos que g ∈ C(R) posee<br />
límites en at infinito. A<strong>de</strong>más, sean αp y βp los<br />
coeficientes <strong>de</strong> Fourier m-ésimos <strong>de</strong> p. Entonces, si<br />
<br />
α2 p + β2 p < 2<br />
|g(+∞) − g(−∞)| , (3)<br />
π<br />
el <strong>problema</strong> (1-2) admite al menos <strong>un</strong>a solución<br />
2π-periodica.<br />
En esta com<strong>un</strong>icación, presentaremos <strong>un</strong>a<br />
generalización <strong>de</strong>l teorema <strong>de</strong> Lazer-Leach <strong>para</strong><br />
N > 1. Es conveniente observar que alg<strong>un</strong>as<br />
dificulta<strong>de</strong>s son <strong>de</strong> esperarse cuando <strong>un</strong>o intenta<br />
exten<strong>de</strong>r este resultado a <strong>un</strong> <strong>sistema</strong> <strong>de</strong> ecuaciones.
Por ejemplo, <strong>para</strong> el caso escalar, no es necesario<br />
suponer que los límites radiales en g(±∞) existan;<br />
Sin embargo, el argumento no se generaliza a <strong>un</strong><br />
sitema. Cuando m = 0, <strong>un</strong> interesante ejemplo fue<br />
dado en el artículo:<br />
R. Ortega and L. Sánchez: Periodic solutions of<br />
forced oscillators with several <strong>de</strong>grees of freedom,<br />
Bull. London Math. Soc. 34 (2002), 308-318.<br />
mostrando que la existencia <strong>de</strong> los límites radiales <strong>de</strong><br />
g es en algún sentido necesaria. Más presisamente, los<br />
autores presenta <strong>un</strong> <strong>sistema</strong> <strong>para</strong> el cual no existen<br />
soluciones periódicas, a<strong>un</strong>que las siguientes<br />
condiciones se cumplen <strong>para</strong> algún R > 0:
• g(u) = p <strong>para</strong> |u| ≥ R.<br />
• El grado <strong>de</strong> la aplicación θR : S N−1 → S N−1<br />
dada por<br />
es diferente <strong>de</strong> 0.<br />
θR(v) =<br />
g(Rv) − p<br />
|g(Rv) − p|<br />
A pesar <strong>de</strong> este ejemplo, mostraremos que la hipótesis<br />
<strong>de</strong> la existencia <strong>de</strong> los límites radiales pue<strong>de</strong><br />
reemplazarse por <strong>un</strong>a condición más débil (ver la<br />
condición (C1) más a<strong>de</strong>lante).
Aplicando métodos <strong>de</strong> grado topológico, obtendremos<br />
soluciones <strong>de</strong> (1-2) bajo condiciones apropiadas <strong>de</strong><br />
tipo Lazer-Leach type. En particular, si la nolinealidad<br />
g tiene límites radiales <strong>un</strong>iformes en infinito, estas<br />
condiciones involucran los m-ésimos coeficientes <strong>de</strong><br />
Fourier <strong>de</strong> <strong>un</strong>a extensión a<strong>de</strong>cuada <strong>de</strong> g la esfera en<br />
infinito. Sin embargo, a diferencia <strong>de</strong>l resultado <strong>de</strong><br />
Nirenberg, nuestra condición (C1) no asume que<br />
todos los límites radiales existan: asumiremos en<br />
cambio la existencia <strong>de</strong> límites superiores, y sólo en<br />
ciertas direcciones específicas.
Primero necesitamos introducir alg<strong>un</strong>as<br />
notaciones.Sea:<br />
H := u ∈ H 1 ([0, 2π], R N ) : u(0) = u(2π) <br />
equipado con la norma usual u := uH1, y sea<br />
D := u ∈ H ∩ H 2 ([0, 2π], R N ) : u ′ (0) = u ′ (2π) .<br />
El operador Lm : D → L 2 ([0, 2π], R N ) se <strong>de</strong>fine<br />
como antes, y su núcleo Vm se pue<strong>de</strong> <strong>de</strong>scribir como:<br />
don<strong>de</strong><br />
Vm := Ker(Lm) = {uw : w = (α,β) ∈ R 2N },
Los coeficientes <strong>de</strong> Fourier m-ésimos <strong>de</strong><br />
u ∈ L 1 ([0, 2π], R N ) serán <strong>de</strong>notados respectivamente<br />
por αu and βu, i.e.<br />
αu = 1<br />
π<br />
βu = 1<br />
π<br />
2π<br />
0<br />
2π<br />
0<br />
cos(mt)u(t) dt,<br />
sin(mt)u(t) dt.<br />
A<strong>de</strong>más si w(u) = (αu,βu), entonces la proyección<br />
ortogonal P <strong>de</strong>l espacio L 2 ([0, 2π], R N ) sobre Vm se<br />
pue<strong>de</strong> <strong>de</strong>finir como Pu = J(w(u)) = u w(u). Es<br />
inmediato que el rango <strong>de</strong> Lm es el complemento<br />
ortogonal <strong>de</strong> Vm con respecto al producto escalar <strong>de</strong>
Entonces, po<strong>de</strong>mos <strong>de</strong>finir <strong>un</strong> inverso a la <strong>de</strong>recha<br />
(compacto) K : R(L) → H <strong>de</strong>l operador Lm, dado<br />
por Kϕ = u, don<strong>de</strong> u ∈ D es la única solución <strong>de</strong>l<br />
<strong>problema</strong>:<br />
u ′′ + m 2 u = ϕ<br />
Pu = 0.<br />
A<strong>de</strong>más, tenemos la siguiente estimación estándar:<br />
Lema 1 Existe <strong>un</strong>a constante c tal que:<br />
<strong>para</strong> todo u ∈ D.<br />
u − PuH 2 ≤ cLm(u)L2
La condición (C1)<br />
En lo sucesivo, vamos a asumir que se verifica la<br />
siguiente condición:<br />
(C1) Existe <strong>un</strong> cubrimiento por abiertos {Uj}j=1...,K<br />
<strong>de</strong> la esfera <strong>un</strong>idad S 2N−1 ⊂ R 2N , y vectores<br />
wj = (α j ,β j ) ∈ S 2N−1 tales que <strong>para</strong> cada w ∈ Uj el<br />
límite<br />
gw,j(t) := lím sup〈g(suw(t)),uwj<br />
s→+∞<br />
(t)〉<br />
es semi-continuou superiormente en w <strong>para</strong> casi todo<br />
t, don<strong>de</strong> 〈·, ·〉 <strong>de</strong>nota el producto interno usual en R N .
Teorema principal<br />
Supongamos que se verifica la condición (C1) y que:<br />
• (C2) Para cada w ∈ S2N−1 existe j ∈ {1,...,K}<br />
tal que<br />
1<br />
π<br />
2π<br />
0<br />
g w,j(t) dt < 〈αp,α j 〉 + 〈βp,β j 〉.<br />
• (C3) Para cada R ≫ 0 el grado <strong>de</strong> Brouwer<br />
<strong>de</strong>gB(G,BR(0), 0) es diferente <strong>de</strong> cero, don<strong>de</strong><br />
G : R 2N → R 2N es la aplicación <strong>de</strong>finida por:<br />
G(w) = w(p) − J −1 P(g ◦ uw).<br />
Entonces el <strong>problema</strong> (1-2) admite al menos <strong>un</strong>a<br />
solución.
Se sigue <strong>de</strong> las <strong>de</strong>finiciones anteriores que G pue<strong>de</strong><br />
expresarse en términos <strong>de</strong> los coeficientes <strong>de</strong> Fourier<br />
m-ésimos <strong>de</strong> la f<strong>un</strong>ción ϕ(t) = g(uw(t)), <strong>de</strong> hecho:<br />
G(w) = (αp − αg◦uw ,βp − βg◦uw )<br />
El teorema anterior tiene <strong>un</strong>a consecuencia inmediata<br />
si asumimos que los líimites radiales<br />
gv = líms→+∞ g(sv) existen <strong>un</strong>iformemente <strong>para</strong><br />
v ∈ S N−1 . De hecho, en este caso, po<strong>de</strong>mos <strong>de</strong>finir<br />
<strong>para</strong> cada t ∈ [0, 2π] y cada w ∈ S 2N−1 el limite<br />
gw(t) := lím<br />
s→+∞ g(suw(t)). (4)
A<strong>un</strong>que uw(t) podría eventualmente ser 0 <strong>para</strong> <strong>un</strong><br />
conj<strong>un</strong>to finito <strong>de</strong> valores <strong>de</strong> t, este “conj<strong>un</strong>to<br />
singular” <strong>de</strong> valores <strong>de</strong> t no juega ningún rol cuando<br />
usamos el teoremas <strong>de</strong> convergencia <strong>de</strong> Lebesgue <strong>para</strong><br />
la integral. Por otro lado, si uw(t) = 0, entonces gw(t)<br />
es continua como f<strong>un</strong>ction <strong>de</strong> w; thus, entonces es<br />
facil ver que las conditiones (C1) y (C2) <strong>de</strong>l teorema<br />
anterior son satifechas.
En este marco, el teorema anterior pue<strong>de</strong> ser<br />
formulado como el resultado <strong>de</strong> Nirenberg’s, en<br />
términos <strong>de</strong> <strong>un</strong>a condición sobre la extensión <strong>de</strong> g a la<br />
esfera en infinito o, siendo más precisos, en términos<br />
<strong>de</strong> los m-ésimos coeficientes <strong>de</strong> Fourier <strong>de</strong> dicha<br />
extensión. De hecho, <strong>para</strong> w ∈ S 2N−1<br />
lím<br />
s→+∞ G(sw) = (αp − αgw ,βp − βgw ) = 0,<br />
Entonces, la aplicación θ : S 2N−1 → S 2N−1 dada por<br />
θ(w) = (αp − αgw ,βp − βgw )<br />
|(αp − αgw ,βp − βgw )|<br />
está bien <strong>de</strong>finida. A partir <strong>de</strong> las propieda<strong>de</strong>s <strong>de</strong>l
Corolario:<br />
Supongamos que los límites radiales gv existen<br />
<strong>un</strong>iformemente <strong>para</strong> v ∈ S N−1 , y <strong>para</strong> cada<br />
w ∈ S 2N−1 <strong>de</strong>finamos la f<strong>un</strong>ción gw(t) por (4).<br />
A<strong>de</strong>más, supongamos que:<br />
1. (αgw ,βgw ) = (αp,βp) <strong>para</strong> todo w ∈ S 2N−1 .<br />
2. <strong>de</strong>g(θ) = 0.<br />
Entonces el <strong>problema</strong> (1-2) admite al menos <strong>un</strong>a<br />
solución.
En el caso particular N = 1, si w = (α,β) ∈ S 1<br />
tenemos que uw(t) = cos(mt − ω), don<strong>de</strong> α = cos(ω)<br />
y β = sin(ω). Se sigue que:<br />
don<strong>de</strong> I + ω<br />
I− ω<br />
g(suw(t)) →<br />
g(+∞) if t ∈ I + ω<br />
g(−∞) if t ∈ I − ω ,<br />
= {t ∈ [0, 2π] : cos(mt − ω) > 0} and<br />
= {t ∈ [0, 2π] : cos(mt − ω) < 0}. En<br />
consecuencia,<br />
gw(t) = g(+∞)χ I + ω (t) + g(−∞)χ I − ω (t),<br />
salvo <strong>para</strong> finitos valores <strong>de</strong> t. Se sigue que,
I ± ω<br />
I ± ω<br />
Por lo tanto,<br />
cos(mt) dt = ±2 cos(ω) = ±2α,<br />
sin(mt) dt = ±2 sin(ω) = ±2β.<br />
lím<br />
s→+∞ G(sw) = (αp,βp) − 2<br />
(g(+∞) −g(−∞))(α,β)<br />
π<br />
a partir <strong>de</strong> lo que se recupera el resultado original<br />
Lazer y Leach.
El corollario anteriror <strong>para</strong> N = 1 permite <strong>un</strong>a<br />
interpretación natural <strong>de</strong>l teorema <strong>de</strong> Lazer-Leach en<br />
términos <strong>de</strong> <strong>un</strong>a integral compleja. De lo anterior<br />
resulta que el grado <strong>de</strong> la aplicación θ : S 1 → S 1 dada<br />
por líms→+∞ G(sw)<br />
|G(sw)|<br />
es equivalente al ínidice en<br />
z0 = αp + iβp <strong>de</strong> la curva γ : [0, 2π] → C <strong>de</strong>finida por<br />
γ(t) = 2<br />
π (g(+∞) − g(−∞))eit .<br />
De (3), resulta que |z0| < |γ(t)|, y que<br />
I(γ,z0) = 1<br />
2πi<br />
2π<br />
0<br />
γ ′ (t)<br />
γ(t) − z0<br />
dt = ±1.
I<strong>de</strong>a <strong>de</strong> la prueba:<br />
Para λ ∈ [0, 1], <strong>de</strong>finimos el operador <strong>de</strong> the Fredholm<br />
Fλ : H → H dado por Fλu = u − Tλu, don<strong>de</strong> el<br />
operador compacto Tλ se <strong>de</strong>fine por:<br />
Tλu = Pu+P(p−g(u))+λK p−g(u)−P(p−g(u)) .<br />
Para λ > 0, u es <strong>un</strong> cero <strong>de</strong> Fλ si y sólo si u ∈ D, y<br />
Lm(u) = λ(p − g(u)).<br />
En consecuencia, buscamos <strong>un</strong> cero <strong>de</strong> F1. Se sigue<br />
<strong>de</strong> las hipótesis que, <strong>para</strong> algún R ≫ 0 a<strong>de</strong>cuado, la<br />
f<strong>un</strong>ción Fλ no se anula en ∂Ω, con Ω = BR(0). Por lo<br />
tanto, el grado <strong>de</strong> Leray-Schau<strong>de</strong>r <strong>de</strong> Fλ en 0 está<br />
<strong>de</strong>finido y, <strong>de</strong>gLS(F1, Ω, 0) = <strong>de</strong>gLS(F0, Ω, 0).<br />
A<strong>de</strong>más, utilizando la <strong>de</strong>finición <strong>de</strong>l grado y el hecho<br />
<strong>de</strong> que
concluimos que:<br />
F0u = u − P(u + p − g(u)),<br />
<strong>de</strong>gLS(F0, Ω, 0) = <strong>de</strong>gB(F0|Vm , Ω ∩ Vm, 0).<br />
Pero, <strong>para</strong> u ∈ Vm es claro que F0u = −P(p − g(u)).<br />
Escribiendo u = uw con w ∈ R 2N <strong>de</strong>ducimos que<br />
J −1 F0J(w) = − (αp,βp) − J −1 P(g(uw)) = −G(w)<br />
Salvo <strong>un</strong> factor (−1) N , el grado <strong>de</strong> F1 en 0 sobre Ω<br />
pue<strong>de</strong> i<strong>de</strong>ntificarse con el grado <strong>de</strong> Brouwer <strong>de</strong> G en 0<br />
<strong>para</strong> <strong>un</strong>a bola <strong>de</strong> R 2N <strong>de</strong> radio gran<strong>de</strong>. Luego, el<br />
resultado se sigue <strong>de</strong> la hipótesis (C2).
<strong>Un</strong> <strong>sistema</strong> débilmente acoplado<br />
Como <strong>un</strong>a aplicación <strong>de</strong>l teorema principal,<br />
consi<strong>de</strong>remos el <strong>sistema</strong><br />
u ′′<br />
i + m 2 ui + ˜gi(ui) + hi(u) = pi(t) i = 1,...,N,<br />
don<strong>de</strong> ˜gi tiene limites en infinito, y hi(u) → 0<br />
<strong>un</strong>iformemente cuando |ui| → +∞. Conviene<br />
observar que, en este caso, los límites radiales<br />
g = ˜g + h no necesariamente existen <strong>para</strong> aquellos<br />
v ∈ S N−1 tales que vi = 0 <strong>para</strong> algún i, pues hi(sv)<br />
no necesariamente converge cuando s → +∞. Sin<br />
embargo, la condición (C1) se verifica, así como la<br />
condición (C1) en el teorema principal.
Luego <strong>de</strong> alg<strong>un</strong>os cálculos, se sigue que si<br />
<br />
(αp) 2 i + (βp) 2 i<br />
< 2<br />
π |˜gi(+∞) − ˜gi(−∞)|<br />
<strong>para</strong> i = 1,...,N, entonces G es homotopicamente<br />
equivalente a la aplicación<br />
(αp,βp) − 2<br />
π<br />
N<br />
[˜gi(+∞) − ˜gi(−∞)](αiei + βieN+i)<br />
i=1<br />
en ∂BR <strong>para</strong> algún R ≫ 0. y utilizando las<br />
propieda<strong>de</strong>s <strong>de</strong>l grado vemos que:<br />
<strong>de</strong>gB(G,BR, 0) = <strong>de</strong>gB(T,BR, 0) = ±1.