08.08.2013 Views

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 ...

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!