La PizarraMatemática.

Grado en Matemáticas

Problemas de contorno y valores propios

Dos datos en un mismo instante fijan una solución regular de segundo orden. Si los ponemos en extremos distintos, pueden ser compatibles, incompatibles o dejar libertad. Los valores propios señalan precisamente los parámetros que permiten una solución homogénea no nula. Los calcularemos y aprenderemos por qué sus modos son ortogonales y cómo responden a un forzamiento.

1. Dos condiciones en extremos distintos pueden dejar un sistema singular

Una ecuación lineal regular de segundo orden tiene dos constantes libres. y(a),y′(a) las fijan de forma única por el teorema inicial. En cambio, y(a),y(b) pueden producir dos ecuaciones dependientes o incompatibles para esas constantes.

Se comprueba el sistema resultante: contar condiciones no basta para afirmar existencia o unicidad. Si ambas condiciones se reducen a la misma igualdad, queda libertad; si exigen valores contradictorios de la misma combinación, no habrá solución.

El mismo intervalo admite infinitas o ninguna

y″+4y=0, y(0)=0.

  1. Familia

    y=Acos(2x)+Bsen(2x); el primer dato fija A=0.
  2. Segundo extremo cero

    Pedir y(π/2)=0 no restringe B, porque senπ=0.
  3. Segundo extremo uno

    Pedir y(π/2)=1 exigiría 0=1.
Resultado

Dos condiciones pueden dejar infinitas soluciones o ninguna.

Lo comprobamos

El intervalo y la ecuación permanecieron iguales; cambió el valor que debía cumplirse en el segundo extremo.

Un problema de contorno que sí fija ambas constantes

y″=0, y(0)=−1,y(2)=3.

  1. Integrar

    y=Ax+B.
  2. Primer dato

    B=−1.
  3. Segundo

    2A−1=3, así A=2.
Resultado

y=2x−1 es la única solución.

Lo comprobamos

Las dos ecuaciones para las constantes fueron independientes y compatibles.

2. Dirichlet, Neumann y Robin cambian el problema y sus valores propios

Dirichlet fija el valor de y; Neumann fija su derivada normal; Robin combina ambos. En [0,L], la normal exterior izquierda apunta hacia menos x y la derecha hacia más x: las derivadas normales son −y′(0) y y′(L). Condiciones homogéneas tienen lado derecho cero.

La derivada normal no tiene el mismo signo en ambos extremos
Las normales apuntan hacia fuera02∂ₙy=−y′(0)∂ₙy=y′(2)El signo depende del extremo.No cambiamos el eje al leer la normal.
Intervalo de cero a dos con normal exterior hacia la izquierda en cero y hacia la derecha en dos. Por eso la derivada normal izquierda es menos la derivada usual en cero, y la derecha es la derivada usual en dos. La variable espacial sigue creciendo hacia la derecha.

El operador no queda especificado solo por la expresión diferencial: también forman parte sus condiciones. Cambiarlas puede añadir un modo constante o incluso valores propios negativos. Los signos de Robin deben escribirse con la convención normal elegida.

Un modo constante cambia de admisible al cambiar el borde

−y″=λy en [0,2].

  1. Neumann homogénea

    y′(0)=y′(2)=0 admite y=1 con λ=0.
  2. Dirichlet homogénea

    y(0)=y(2)=0 excluye cualquier constante no nula.
  3. Conclusión

    La expresión −y″ es la misma, pero no el espectro.

Resultado

Las condiciones son parte del problema de valores propios.

Lo comprobamos

La derivada normal izquierda tiene signo contrario, aunque para dato cero ambas escrituras de Neumann coincidan.

Robin permite ver por qué importan los signos

y(0)−y′(0)=0 y −y(2)+y′(2)=0.

  1. Lectura normal

    La izquierda es y+∂ₙy=0 y la derecha −y+∂ₙy=0.
  2. Probar

    y=eˣ cumple las dos, porque y′=y.
  3. Ecuación

    −y″=−y, así λ=−1 es un valor propio con estas condiciones.

Resultado

No toda elección de borde obliga a valores propios positivos.

Lo comprobamos

Es un problema separado de Robin con coeficientes reales; no contradice la positividad que veremos para Dirichlet con energía no negativa.

3. Un valor propio requiere una solución distinta de cero

En L[y]=λwy con condiciones homogéneas, y=0 siempre sirve. λ es valor propio únicamente si también existe alguna solución no nula. Esa función se llama autofunción y describe un modo permitido por la ecuación y el borde simultáneamente.

Multiplicar una autofunción por una constante no nula no cambia el valor propio. Por eso buscamos direcciones de funciones y después elegimos una normalización. Multiplicar por cero produce la solución cero, que no es una autofunción.

Un parámetro que no es propio para Dirichlet

−y″=0, y(0)=y(2)=0.

  1. General

    y=Ax+B.
  2. Extremos

    B=0 y 2A=0.
  3. Resultado

    Solo queda y=0.
Resultado

λ=0 no es valor propio de este problema.

Lo comprobamos

La existencia de la solución cero no satisface la condición de modo no trivial.

Una familia de múltiplos de un mismo modo

−y″=(π²/4)y en [0,2], extremos cero.

  1. Candidata

    y=sen(πx/2).
  2. Verificar

    Su segunda derivada es −(π²/4)y y se anula en cero y dos.

  3. Múltiplos

    A y también sirve para cada A no nulo.

Resultado

π²/4 sí es valor propio y el modo puede escalarse.

Lo comprobamos

No aparecen valores propios diferentes por cambiar solo la amplitud de la misma autofunción.

4. Para el modelo Dirichlet hay que estudiar los tres signos de λ

En −y″=λy sobre (0,L), con extremos cero, si λ=k²>0 usamos A cos(kx)+B sen(kx). El primer extremo da A=0; una función no nula exige sen(kL)=0.

Así k=nπ/L, n≥1, y λₙ=(nπ/L)².
Los extremos fijos solo permiten algunas frecuencias
−y″=λy ; y(0)=y(2)=0n=1 ; λ=π²/402n=2 ; λ=π²02n=3 ; λ=9π²/402Más n: más oscilaciones y mayor λ.
Tres autofunciones seno de n pi x medios en el intervalo de longitud dos, con n uno dos tres. Los valores propios son pi cuadrado cuartos, pi cuadrado y nueve pi cuadrado cuartos. Todas se anulan en ambos extremos y el número de ceros interiores crece con n.

Si λ=0, las funciones afines quedan anuladas por los extremos. Si λ=−k²<0, la base es cosh(kx),senh(kx); el primer extremo elimina cosh y el segundo elimina senh porque kL>0. Con esto no se omite ningún signo. Cambiar la segunda condición cambia la ecuación que deben satisfacer las frecuencias.

Los tres primeros modos del dibujo

L=2 y Dirichlet en ambos extremos.

  1. Frecuencias

    k=nπ/2, n=1,2,3,….
  2. Valores

    π²/4, π², 9π²/4,….

  3. Funciones

    sen(nπx/2); el modo n tiene n−1 ceros interiores.

Resultado

El borde permite solo frecuencias discretas.

Lo comprobamos

Todas las curvas del dibujo se anulan en cero y dos, pero oscilan con ritmos distintos.

Dirichlet a un lado y Neumann al otro

y(0)=0,y′(2)=0.

  1. Caso positivo

    y=Bsen(kx); y′(2)=Bkcos(2k)=0.
  2. Frecuencias

    k=(2n+1)π/4, n=0,1,2,….
  3. Otros signos

    Una afín con y(0)=0 y y′(2)=0 es nula; una senh tiene derivada kcosh(2k) no nula salvo amplitud cero.

Resultado

Los valores son ((2n+1)π/4)².

Lo comprobamos

El primer valor es π²/16, diferente del π²/4 del borde Dirichlet doble.

5. Sturm–Liouville incorpora coeficientes, peso y condiciones autoadjuntas

La forma regular es L[y]=−(py′)′+qy=λwy en [a,b].
Suponemos p C¹ positiva y q,w reales continuas, con w>0.

Añadimos condiciones homogéneas separadas reales, no triviales, que anulen p(u v′−u′v) al evaluar los bordes.

El producto ponderado es ⟨u,v⟩w=∫w u conjugado(v).

La positividad de w garantiza norma positiva para una función no nula. En ese espacio, el operador asociado es T[y]=L[y]/w con el dominio de condiciones elegido. La integración por partes mostrará su simetría; estas hipótesis sustentan el teorema espectral regular.

El peso cambia la escala de los valores propios

−2y″=3λy en [0,2], extremos cero.

  1. Datos

    p=2,q=0,w=3, todos regulares con p,w positivos.
  2. Normalizar ecuación

    −y″=(3λ/2)y.
  3. Frecuencias

    3λ/2=(nπ/2)², así λₙ=n²π²/6.
Resultado

Las autofunciones son los mismos senos, pero los valores propios cambiaron.

Lo comprobamos

Confundir el peso con uno habría dado el espectro equivocado.

Un peso que no cumple la hipótesis de norma positiva

En [0,2], intenta usar w=−1.

  1. Forma

    La expresión ∫w|u|² sería negativa para toda u no nula.

  2. Prueba simple

    Con u=1, da −2.
  3. Consecuencia

    No es el producto interior positivo del marco regular que acabamos de definir.

Resultado

No se aplica ese teorema con sus hipótesis sin cambiar la formulación.

Lo comprobamos

Que un símbolo aparezca multiplicando λ no basta para que sea un peso positivo admisible.

6. La identidad de Green produce realidad y ortogonalidad

Si L[u]=λwu y L[v]=μwv, multiplicamos la primera por v, la segunda por u, restamos e integramos.

Integrar por partes da (λ−μ)∫wuv=[p(u v′−u′v)] entre a y b.

Las condiciones autoadjuntas anulan el borde.

Para λ,μ reales distintos, el producto ponderado es cero. Para una autofunción compleja usamos su conjugada en la identidad: (λ−conjugado(λ))∫w|u|²=0. Como la norma es positiva, λ debe ser real. Ni la realidad ni la ortogonalidad se deducen solo de que un ejemplo use senos.

Comprobar dos modos ortogonales en el modelo sencillo

u=sen(πx/2),v=sen(πx), intervalo [0,2].

  1. Producto

    uv=[cos(πx/2)−cos(3πx/2)]/2.
  2. Integrar

    Las primitivas son múltiplos de los senos de esos ángulos.

  3. Extremos

    En cero y dos, los senos valen cero: ∫uv=0.

Resultado

Los valores propios π²/4 y π² son distintos y los modos son ortogonales.

Lo comprobamos

La identidad trigonométrica comprueba el resultado general con una cuenta concreta.

Un peso no constante participa realmente en la cancelación

p=w=eˣ,q=0, con extremos cero en [0,2].

  1. Cambio

    Escribe y=e^(−x/2)v; −y″−y′=λy se vuelve −v″+(1/4)v=λv.
  2. Modos

    yₙ=e^(−x/2)sen(nπx/2), λₙ=(nπ/2)²+1/4.
  3. Producto ponderado

    eˣyₘyₙ=sen(mπx/2)sen(nπx/2), cuya integral es cero para m≠n.

Resultado

La ortogonalidad correcta usa el peso eˣ.

Lo comprobamos

Las dos exponenciales de las autofunciones se cancelan con el peso; no se elimina este sin una razón.

7. Normalizar modos permite expandir funciones con una noción precisa de convergencia

El teorema espectral regular, bajo las hipótesis anteriores y bordes separados autoadjuntos, da una sucesión real creciente de valores propios simples que tiende a +∞. Sus autofunciones pueden normalizarse y forman una base completa en L² con peso w.

Normalizamos dividiendo por √(∫w|yₙ|²). Los coeficientes de f son sus productos interiores con los modos normalizados. La convergencia general es en norma L²w. Convergencia puntual, uniformidad y derivar la serie requieren hipótesis adicionales sobre datos y regularidad.

Calcular una normalización, no solo nombrarla

Dirichlet en [0,3], yₙ=sen(nπx/3), w=1.

  1. Cuadrado

    sen²θ=(1−cos(2θ))/2.
  2. Norma

    Integrar da ∫₀³yₙ²=3/2, porque el término seno de borde se anula.

  3. Dividir

    ϕₙ=√(2/3)sen(nπx/3).
Resultado

Cada ϕₙ tiene norma uno.

Lo comprobamos

La amplitud original de un modo era libre; la normalización fija una escala útil para coeficientes.

El peso cambia también los coeficientes de expansión

En [0,2], w=3 y f=2sen(πx/2)−3sen(πx).

  1. Modos unitarios

    ϕₙ=sen(nπx/2)/√3, ya que la integral sin peso de cada seno al cuadrado vale uno.

  2. Coeficientes

    c₁=2√3, c₂=−3√3 y los demás cero por ortogonalidad.
  3. Reconstruir

    c₁ϕ₁+c₂ϕ₂ devuelve exactamente f.

Resultado

Los coeficientes dependen de la normalización y del producto interior.

Lo comprobamos

No se confunden los números de una combinación en modos sin normalizar con los de una base ortonormal.

8. La energía comprueba el signo y aproxima el primer valor propio

Con bordes Dirichlet, multiplicar L[y]=λwy por y e integrar da λ=∫[p(y′)²+qy²]/∫wy². El borde se anula y el denominador es positivo. Si q≥0, el numerador de una autofunción no nula es positivo: si fuese cero, y′ sería cero y los extremos obligarían y=0.

La caracterización variacional del primer valor propio es el mínimo de este cociente de Rayleigh entre funciones de prueba admisibles no nulas con extremos cero y energía finita. Una prueba concreta ofrece una cota superior. Con otros bordes hay que conservar sus términos, y con q negativo no se deduce positividad.

Una aproximación de energía cercana al primer valor

−y″ en [0,2], extremos cero; prueba v=x(2−x).

  1. Numerador

    ∫₀²(v′)²=∫₀²(2−2x)²=8/3.
  2. Denominador

    ∫₀²v²=16/15.
  3. Cociente

    R[v]=(8/3)/(16/15)=5/2, así λ₁≤2,5.
Resultado

El valor exacto π²/4, aproximadamente 2,467, queda por debajo de esa cota.

Lo comprobamos

La función de prueba cumple los dos extremos; no cualquier función produce una cota del mismo problema.

La positividad tiene condiciones que no se omiten

Compara dos cambios del problema.

  1. Potencial negativo

    En [0,π], Dirichlet y L=−y″−2y, el modo sen x tiene λ=−1.
  2. Otro borde

    Con q=0 y los Robin de la segunda sección, eˣ también tiene λ=−1.

  3. Razón

    En el primero q no era no negativo; en el segundo no desaparece el término de energía que exige Dirichlet.

Resultado

No hay una regla de positividad independiente del operador y sus condiciones.

Lo comprobamos

Ambos ejemplos respetan su ecuación, pero no las hipótesis de la conclusión positiva anterior.

9. En un valor propio, el forzamiento debe evitar el modo resonante

Para (T−λ)y=g, con T=L/w y λ propio, la identidad autoadjunta exige ∫w g ϕ=0 para cada modo ϕ del núcleo.

Si escribimos en cambio L[y]−λwy=f, entonces g=f/w y la condición es ∫fϕ=0.

Se especifica la formulación antes de colocar el peso.

En este marco regular autoadjunto, esa compatibilidad es también suficiente para el forzamiento apropiado. La solución no será única: se puede añadir cualquier modo homogéneo del núcleo. Imponer que y sea ortogonal al núcleo selecciona una representante. Fuera del espectro, el problema homogéneo solo tiene cero y el problema forzado tiene solución única.

Un forzamiento incompatible

−y″−(π²/4)y=sen(πx/2), en (0,2), extremos cero.

  1. Modo resonante

    ϕ=sen(πx/2) está en el núcleo.
  2. Multiplicar e integrar

    La autoadjunción deja cero del lado izquierdo.

  3. Derecha

    ∫₀²ϕ²=1, distinto de cero.
Resultado

No existe solución con esos bordes.

Lo comprobamos

No es un fallo de nuestro método de cálculo: la identidad obliga a una contradicción.

Un forzamiento compatible deja libertad homogénea

Mismo operador, lado derecho sen(πx).

  1. Ortogonalidad

    El segundo seno es ortogonal al primero, así cumple la compatibilidad.

  2. Particular

    yₚ=4sen(πx)/(3π²), porque π²−π²/4=3π²/4.
  3. Todas

    y=yₚ+Csen(πx/2).
    Si pedimos ∫yϕ=0, entonces C=0.
Resultado

Hay soluciones, pero hace falta una normalización adicional para unicidad.

Lo comprobamos

Sustituir en el operador y comprobar los dos extremos verifica cada miembro de la familia.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Espectro con extremos cero

Resuelve −y″=λy, y(0)=y(π)=0.

  1. Objetivo

    La solución cero existe para todo λ; buscamos soluciones no nulas.

  2. Signos no positivos

    λ=0 da una afín anulada por los extremos; λ=−k² da cosh y senh cuyos coeficientes también se anulan.

  3. Caso positivo

    λ=k²,k>0: y=Acos(kx)+Bsen(kx).
  4. Extremo cero

    A=0 y se requiere B no nulo.
  5. Extremo π

    Bsen(kπ)=0 fuerza k=n entero positivo, así λₙ=n².
  6. Modos

    yₙ=sen(nx), n≥1, y sus múltiplos no nulos.

    La segunda derivada y los dos extremos verifican el resultado; fuera de esos parámetros solo queda cero.

Resultado

λₙ=n²; yₙ=sin(nx).

Lo comprobamos

−(sen nx)″=n²sen nx y ambos extremos son cero para n≥1. Los casos no positivos se han examinado por separado.

EJEMPLO RESUELTO 2

Un problema de contorno compatible

Resuelve y″=0, y(0)=1,y(1)=3.

  1. Integrar

    y″=0 da y′=A.
  2. Segunda integración

    y=Ax+B.
  3. Dato en cero

    B=1.
  4. Dato en uno

    A+1=3, por tanto A=2.
  5. Compatibilidad

    Las dos condiciones determinan constantes de forma independiente.

  6. Verificar

    y=2x+1 tiene segunda derivada cero, valor uno en cero y tres en uno. Este problema es compatible y tiene solución única.

Resultado

y=2x+1.

Lo comprobamos

La recta une exactamente los dos datos. El rótulo de la fuente que lo llamaba incompatible se corrige porque el cálculo demuestra una solución única.

Los símbolos que hemos usado

[a,b]
Intervalo espacial finito con a<b donde se plantea el problema regular.
Dirichlet
Condición que prescribe el valor de y en un extremo.
Neumann y Robin
Condiciones sobre la derivada normal o una combinación lineal de valor y derivada.
λ
Parámetro espectral; es valor propio si hay una solución no nula de la ecuación y condiciones homogéneas.
y_n
Autofunción asociada al valor propio λ_n, definida salvo multiplicación por una constante no nula.
L[y]=−(py′)′+qy
Operador diferencial de Sturm–Liouville con p>0 y coeficientes regulares.
⟨u,v⟩_w
Producto escalar ponderado ∫_a^b w u conjugado(v), con peso w>0.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar