Teoría de Sturm-Liouville

Teoría de Sturm-Liouville

En matemáticas, una ecuación de Sturm-Liouville, que toma su nombre de Jacques Charles François Sturm (1803-1855) y Joseph Liouville (1809-1882), es una ecuación diferencial lineal de segundo orden de la forma

(1)  -\frac{d}{dx}\left[p(x)\frac{dy}{ dx}\right]+q(x)y=\lambda w(x)y

donde las funciones \ p(x), q(x), w(x) están prescritas, y en el caso más simple son continuas en un intervalo finito cerrado [a,b]. El problema generalmente viene formulado con condiciones de frontera, es decir, valores específicos de y y/o \frac{dy}{dx} en los extremos a,b. La función w(x) es llamada función de densidad o función de peso.

El valor de λ no se especifica en la ecuación; el encontrar los valores λ donde exista una solución no trivial de la ecuación que satisfaga condiciones de frontera se denomina el problema de Sturm-Liouville (S-L).

Tales valores de λ son llamados valores propios o eigenvalores del problema de S-L que plantea (1) conjuntamente con las condiciones de frontera. Las soluciones correspondientes son las funciones propias o las eigenfunciones o los eigenvectores del problema. Bajo suposiciones normales en los coeficientes de las funciones p(x),q(x),w(x), éstas inducen operadores diferenciales hermíticos en algunas funciones definidas por las condiciones de frontera. La teoría resultante de la existencia y el comportamiento asintótico de los valores propios, la teoría cualitativa correspondiente de las funciones propias y sus funciones adecuadas completas se conoce como teoría de Sturm-Liouville. Esta teoría es importante en matemática aplicada, donde los problemas S-L ocurren muy comúnmente, particularmente al resolver ecuaciones diferenciales parciales con separación de variables.

Contenido

Teoría de Sturm-Liouville

Cuando las condiciones de frontera son regulares de la forma

(2)  y(a)\cos \alpha - p(a)y^{\prime}(a)\sin \alpha = 0

(3)  y(b)\cos \beta - p(b)y^{\prime}(b)\sin \beta = 0

donde p(x) es diferenciable, las funciones p(x),q(x),w(x) son continuas y las funciones \ p(x), w(x) son positivas sobre el intervalo [a,b], y los valores α,β están en el intervalo [0,π[ la teoría nos indica que

  • Los valores propios λn del problema de S-L, son valores reales y bien ordenados en el sentido de que \lambda_1 < \lambda_2 < \lambda_3 < \cdots < \lambda_n < \cdots \to \infty; .
  • A cada valor propio λn le corresponde una única función propia \ y_n(x) y \ y_n(x) tiene exactamente \ n-1 ceros en la frontera \ (a,b).
  • Las funciones propias son mutuamente ortogonales y satisfacen la relación de ortogonalidad

(4)  \int_{a}^{b}y_n(x)y_m(x)w(x)\,dx = 0 , m \ne n,

donde \ w(x) es la función de peso.

  • Un conjunto ortonormal puede ser formado si el conjunto de funciones propias satisface la relación de ortogonalidad

(5)  \int_{a}^{b}y_n(x)y_m(x)w(x)\,dx = \delta_{mn},

donde \ \delta_{mn} es la delta de Kronecker.

  • Los valores propios del problema de S-L pueden ser caracterizados por el cociente de Rayleigh

 \lambda_n = \frac{-p(x) y_{n}(x) y' _{n}(x)|_a^b + \int_a^b [p y'_{n}(x)^2 + q y_{n}(x)^2]\,dx}{\int_a^b y_{n}(x)^2 w(x)\, dx}

.

Forma de Sturm-Liouville

La ecuación diferencial

 - {d\over dx}\left[p(x){d\over dx}y(x)\right]+q(x)y(x)=\lambda w(x)y(x)

se dice que es de la forma de S-L o de la forma autoadjunta. Toda ecuación diferencial ordinaria lineal de segundo orden puede llevada a esta forma al multiplicarle por un factor integrante apropiado.

Ejemplos

\ x^2y''+xy'+(\lambda^2x^2-\nu^2)y=0

puede ser escrita en la forma de S-L así:

(xy')'+(\lambda^2 x-\nu^2/x)y=0.\,

  • La ecuación de Legendre:

(1-x^2)y''-2xy'+\nu(\nu+1)y=0\;\!

puede ser transformada fácilmente en una forma de Sturm-Lioville, ya que \scriptstyle (1 - x^2)' = -2x; así la ecuación de Legendre equivalente es:

[(1-x^2)y']'+\nu(\nu+1)y=0\;\!

  • Otro ejemplo simple es una ecuación diferencial de la forma:

\ x^3y''-xy'+2y=0

Si dividimos por \ x^3 tenemos:

y''-{x\over x^3}y'+{2\over x^3}y=0

Multiplicando por un factor integrante:

e^{\int -{x / x^3}\,dx}=e^{\int -{1 / x^2}\, dx}=e^{1 / x},

nos da

e^{1 / x}y''-{e^{1 / x} \over x^2} y'+ {2 e^{1 / x} \over x^3} y = 0

que puede ponerse fácilmente en la forma de S-L así:

D e^{1 / x} = -{e^{1 / x} \over x^2}

que es equivalente a decir:

(e^{1 / x}y')'+{2 e^{1 / x} \over x^3} y =0.

  • En general, dada una ecuación diferencial

P(x)y''+Q(x)y'+R(x)y=0\,

dividida para  \ P(x), multiplicada por un factor integrante tenemos la forma de S-L:

e^{\int {Q(x) / P(x)}\,dx},

Operadores diferenciales Sturm-Liouville

El operador lineal:

(6) L  u  =-{d\over dx}\left[p(x){du\over dx}\right]+q(x)u

puede ser vista como la transformación de una función \ u en otra función \ Lu. Se puede estudiar éste operador lineal en el contexto del análisis funcional. Si ponemos  \ w = 1 en la ecuación (1), podemos escribirla:

(7) L  u  = \lambda u \,.

Éste es precisamente un problema de valores propios; donde se trata de hallar valores propios λ y vectores propios \ u del operador \ L. Sin embargo, también se debe incluir las condiciones de frontera. Como ejemplo se dirá que vamos a evaluar el problema en el intervalo \ [0,1] y se pondra las condiciones de frontera \ u(0) = u(1) = 0.

La importancia de problemas de valores propios esta en el hecho que nos ayuda a resolver problemas asociados inhomogéneos:

L u  = f \,

en el intervalo \ [0,1]

u = 0 \,

en 0 y 1.

Aquí \ f es la función en el espacio \ L^2. Si una solución \ u existe y es única, se la puede escribir de la forma:

u = A f \,

porque la transformación de \ f a \ u debe ser lineal. Ahora se observa que el hallar los vectores propios y los valores propios de \ A es esencialmente lo mismo que hallar los vectores y valores propios de \ L. Efectivamente, si \ u es un vector propio de \ L con valores propios λ debe existir un \ u que también es vector propio de \ A con valores propios \frac{1}{u}.


Ejemplo

Se busca una función u(x) que resuelva el siguiente problema S-L:

(8)  L  u  = \frac{d^2u}{dx^2} = \lambda u

en donde las incógnitas son λ and u(x). Tomemos por ejemplo las condiciones de frontera

 u(0) = u(\pi) = 0 \,

y observemos que si k es cualquier entero, entonces la función

 u(x) = \sen kx \,

es una solución con eigenvalor λ = −k2. Sabemos que las soluciones de un problema S-L forman una base ortogonal, y de la teoría de las series de Fourier sabemos que este conjunto de funciones sinusoidales es una base ortoganal. Dado que las bases otorgonales por definición son máximas, concluimos que este problema S-L no tiene más eigenvectores.

Con base en esto resolvamos el problema inhomogéneo

  L u=x\,

con  x\in(0,\pi) y con las mismas condiciones de frontera. En este caso debemos poner f(x) = x en formas de serie de Fourier. El lector puede verificar, sea integrando ∫exp(ikx)x dx o consultando una tabla de transformadas de Fourier, que

L  u  =\sum_{k=1}^{\infty}-2\frac{(-1)^k}{k}\sin kx.

Esta serie de Fourier es problemática por sus malas propiedades de convergencia: no está claro a priori si converge puntualmente. Por el hecho (de la teoría de análisis de Fourier) que los coefientes son cuadrado-sumables, la serie converge en L2, lo cual basta para la presente discusión. Mencionamos para el lector interesado que podemos aplicar un resultado que dice que la serie de Fourier converge en cada punto de diferenciabilidad, y en los puntos de salto (por ejemplo (la función x, considerada como función periódica, tiene un salto en n) converge al promedio de los límites izquierdo y derecho.

Por lo tanto, por la fórmula (8) obtenemos que la solución es

u=\sum_{k=1}^{\infty}2\frac{(-1)^k}{k^3}\sin kx.

Podríamos haber encontrado la solución con antidiferenciación. Esta técnica da u=(x32x)/6, cuya serie de Fourier concuerda con la solución que encontramos. La técnica de antidiferención generalmente no es útil para las ecuaciones diferenciales de varias variables.

Aplicación a los modos de vibración normales

Ciertas ecuaciones en derivadas parciales pueden resolverse con la ayuda de la teoría de S-L. Supongamos que nos interesan los modos de vibración de una membrana delgada, sostenida en un marco rectangular, 0≤xL1, 0≤yL2. El desplazamiento vertical W(x, y, t) de la membrana es gobernada por la ecuación de onda:

\frac{\partial^2W}{\partial x^2}+\frac{\partial^2W}{\partial y^2} = \frac{1}{c^2}\frac{\partial^2W}{\partial t^2}.

El método de separación de variables sugiere buscar primero soluciones de la forma sencilla W = X(x) × Y(y) × T(t). Para una tal función W la ecuación en derivadas parciales puede escribirse como X"/X + Y"/Y = (1/c2)T"/T. Dado que los tres términos de esta ecuación son funciones de x,y,t por separado, deben ser constantes. Por ejemplo, el primer término da X"=λX donde λ es constante. Las condiciones de frontera ("sostenida en un marco rectangular") son W=0 when x=0, L1 ó y = 0, L2 lo cual define los problemas de S-L eigenvalue lo más sencillos posibles como en el ejemplo y dan las "soluciones de modos normales" para W con dependencia armónica del tiempo,

W_{mn}(x,y,t) = A_{mn}\sin\left(\frac{m\pi x}{L_1}\right)\sin\left(\frac{n\pi y}{L_2}\right)\cos\left(\omega_{mn}t\right)

donde m,n son enteros no nulos, Amn son constantes arbitrarias y

\omega^2_{mn} = c^2 \left(\frac{m^2\pi^2}{L_1^2}+\frac{n^2\pi^2}{L_2^2}\right).

Las funciones Wmn forman una base del espacio de Hilbert de soluciones (generalizadas) de la ecuación de onda; eso es, una solución arbitraria W puede descomponerse como suma de estos modos, que vibran a frecuencias individuales ωmn. Esta representación puede requerir una suma infinita convergente.

Operadores de Sturm-Liouville como operadores Hermíticos

Muchas de las propiedades de los operadores de S-L vienen del hecho que éstos son operadores hermíticos con respecto al producto interno:

 \langle u, v \rangle_w = \int_a^b w(x)u(x)v(x)dx.

Y así los valores propios de los operadores de S-L son reales y que las funciones propias corresponden a diferentes valores propios son ortogonales.

Representación de Soluciones y Cálculo Numérico

La ecuación de S-L (1) con condiciones de frontera puede resolverse en la práctica por una variedad de métodos numéricos. En casos difíciles puede ser necesario llevar a cabo los cálculos intermedios con varios cientos de dígitos decimales para lograr los eigenvalores correctamente a unas cuantas cifras.

1. Métodos de disparo.[1] [2] Estos métodos funcionan adivinando un valor de λ, resolviendo un problema en valores iniciales definido por las condiciones de frontera en un extremo, digamos a, del intervalo [a,b], comparando el valor que esta solución toma en el otro extremo b con la otra condicion de frontera deseada y finalmente aumentando o reduciendo λ segun sea necesario para corregir el valor original. Esta estrategia no es aplicable para encontrar eigenvalores complejos.

2. Método de Diferencias Finitas.

3. El método de Series de Potencias del Parámetro Espectral[3] (SPPS por sus siglas en inglés) aplica una generalización del seguiente hecho sobre las ecuaciones en derivadas ordinarias de segundo orden: si y es una solución que no se anula en ningún punto de [a,b], entonces la función

 y(x) \int_a^x \frac{dt}{p(t)y(t)^2}

es una solución de la misma ecuación y es linealmente independiente de y. Además, todas las soluciones son combinaciones lineales de estas dos soluciones. En el algoritmo SPPS, uno debe comenzar con un valor arbitrario λ0* (a menudo λ0*=0; no tiene que ser un eigenvalor) y con cualquier solución y0 de (1) con λ=λ0* la cual no se anule en [a,b]. (Discusión abajo sobre maneras de encontrar y0 y λ0* apropiados.) Dos sucesiones de funciones X(n)(t), X~(n)(t) en [a,b], que se llamarán integrales iteradas, se definen de manera recursiva como sigue. Primero, cuando n=0, se toman iguales a 1 idénticamente en [a,b]. Para obtener las siguientes funciones se multiplican alternativamente por 1/(py02) y wy02 y luego se integran, específicamente

(9)  X^{(n)}(t) = - \int_a^x  X^{(n-1)}(t) p(t)^{-1}  y_0(t)^{-2}\,dt para n impar,  \int_a^x  X^{(n-1)}(t)y_0(t)^{2} w(t) \,dt para n par,

(10)  \tilde X^{(n)}(t) = \int_a^x   \tilde X^{(n-1)}(t)y_0(t)^{2} w(t)\,dt para n impar,  - \int_a^x   \tilde X^{(n-1)}(t) p(t)^{-1}  y_0(t)^{-2} \,dt para n par,

cuando n>0. Las integrales iteradas así obtenidas se aplican ahora como coeficientes en las dos siguientes series de potencias en λ:

 u_0 = y_0 \sum_{k=0}^\infty (\lambda-\lambda_0^*)^k \tilde X^{(2k)} y  u_1 = y_0 \sum_{k=0}^\infty (\lambda-\lambda_0^*)^k X^{(2k+1)}  .

Entonces para cualquier λ (real o complejo), u0 y u1 son soluciones linealmente dependientes de la ecuación (1) correspondiente. (Las funciones p(x) y q(x) participan en esta construcción de forma implícita por su influencia en la elección de y0.)

Ahora se escogen coeficientes c0, c1 de manera que la combinación y=c0u0 + c1u1 satisfaga la primera condición de frontera (2). Esto es fácil de hacer porque X(n)(a)=0 y X~(n)(a)=0, para n>0. Los valores de X(n)(b) y X~(n)(b) dan los valores de u0(b) y u1(b) y de las derivadas u0'(b) y u1'(b), entonces la segunda condición de frontera (3) se convierte en una serie de potencias en λ. Para el trabajo numérico uno puede truncar esta serie a un número finito de términos y obtener un polinomio en λ calculable, cuyas raíces son aproximaciones de los eigenvalores que se buscan.

Cuando λ= λ0, esto se reduce a la construcción original descrita arriba para una solución linealmente independiente a una dada. Las representaciones (9),(10) también tienen aplicaciones teóricas en la teoría de S-L.[3]

Construcción de una solución que no se anula

El mismo método SPPS puede usarse para encontrar una solución inicial y0. Considérese la ecuación(p'y)'=μqy; i.e., q, w, y λ se reemplazan en (1) con 0, -q y μ respectivamente. Entonces la función constante 1 es una solución que no se anula, con respecto al eigenvalor μ0=0. Mientras no hay grantía que u0 ó u1 no se anule, la función compleja y0=u0+iu1 nunca se anulará porque dos soluciones linealmentes de una ecuación S-L regular no pueden anularse simultáneamente (como consecuencia de la unicidad de soluciones para problemas de valores iniciales). Este truco da una solución y0 de (1) para el valor λ0=0. En la práctica, si (1) tiene coeficientes reales, las soluciones basadas en y0 tendrán partes imaginarias muy pequeñas que se tendrán que descartar.

Véase también

Referencias

  1. J. D. Pryce, Numerical Solution of Sturm-Liouville Problems, Clarendon Press, Oxford, 1993.
  2. V. Ledoux, M. Van Daele, G. Vanden Berghe, "Efficient computation of high index Sturm-Liouville eigenvalues for problems in physics," Comput. Phys. Comm. 180, 2009, 532–554.
  3. a b V. V. Kravchenko, R. M. Porter, "Spectral parameter power series for Sturm-Liouville problems," Mathematical Methods in the Applied Sciences (MMAS) 33, 2010, 459-468

Wikimedia foundation. 2010.

Mira otros diccionarios:

  • Sturm — significa tormenta en idioma alemán. En dialecto suizo alemán también se puede emplear como equivalente de vértigo (mir ist sturm). También puede referirse a alguno de los siguientes conceptos: Contenido 1 Miscelánea 2 Apellido 2.1 Personajes …   Wikipedia Español

  • Joseph Liouville — Joseph Liouville. Nacimiento 24 de marzo de 1809 …   Wikipedia Español

  • Jacques Charles François Sturm — Nacimiento 29 de septiembre de 1803 …   Wikipedia Español

  • Método de separación de variables — El método de separación de variables se refiere a un procedimiento para encontrar una solución completa particular para ciertos problemas que involucran ecuaciones en derivadas parciales como serie cuyos términos son el producto de funciones que… …   Wikipedia Español

  • Polinomios ortogonales — Saltar a navegación, búsqueda Los polinomios ortogonales son conjuntos de polinomios que forman una base ortogonal son importantes porque aparecen en la teoría de ecuaciones diferenciales, muy especialmente en la teoría de Sturm Liouville, la… …   Wikipedia Español

  • Ecuación de onda — Saltar a navegación, búsqueda La ecuación de onda es una importante ecuación diferencial parcial lineal de segundo orden que describe la propagación de una variedad de ondas, como las ondas sonoras, las ondas de luz y las ondas en el agua. Es… …   Wikipedia Español

  • Modo normal — Varios modos normales de una red unidimensional. Un modo normal de un sistema oscilatorio es la frecuencia a la cual la estructura deformable oscilará al ser perturbada. Los modos normales son también llamados frecuencias naturales o frecuencias… …   Wikipedia Español

  • Armónicos esféricos — Saltar a navegación, búsqueda En matemáticas, los armónicos esféricos son funciones armónicas que representan la variación espacial de un conjunto ortogonal de soluciones de la ecuación de Laplace cuando la solución se expresa en coordenadas… …   Wikipedia Español

  • Muelle elástico — Saltar a navegación, búsqueda Para otros usos de este término, véase Muelle. «Resorte» redirige aquí. Para otras acepciones, véase Resorte (desambiguación). Muelles de tracción. Se conoce como …   Wikipedia Español

  • Resorte — Para otros usos de este término, véase Resorte (desambiguación). Muelles de tracción. Se conoce como resorte o muelle a un operador elástico capaz de almacenar energía y desprenderse de ella sin sufrir deformación permanente cuando cesan las… …   Wikipedia Español

Compartir el artículo y extractos

Link directo
Do a right-click on the link above
and select “Copy Link”