Sturm-Liouville-Problem

Ein klassisches Sturm-Liouville-Problem (nach Joseph Liouville und Charles-François Sturm) ist folgendes Eigenwertproblem aus der Analysis: Finde alle komplexen Zahlen λ, für die die Differentialgleichung

$-\left( p \cdot \psi' \right)' + q \cdot \psi = \lambda \cdot w \cdot \psi$

auf dem Intervall (a,b) eine Lösung besitzt, die den Randbedingungen

$\cos(\alpha) \psi(a) + \sin(\alpha) p(a)\psi'(a) = 0,\quad \cos(\beta) \psi(b) + \sin(\beta) p(b)\psi'(b) = 0$

genügt.

Führt man den linearen Operator der Form

$L = \frac{1}{w} \left( -\frac{d}{dx} \, p\, \frac{d}{dx} +q \right)$

ein, den Sturm-Liouville-Operator, so kann die Eigenwertgleichung Lψ = λψ mithilfe von Methoden aus der Funktionalanalysis (Spektraltheorie) im Hilbertraum der bezüglich der Gewichtsfunktion w quadratintegrierbaren Funktionen behandelt werden.

Ist das Intervall kompakt und sind die Koeffizienten w,p ^{− 1},q integrierbar, so spricht man von einem regulären Sturm-Liouville-Problem. Ist das Intervall unbeschränkt oder sind die Koeffizienten nur lokal integrierbar, so spricht man von einem singulären Sturm-Liouville-Problem.

Reguläre Sturm-Liouville-Probleme

Die Eigenwertgleichung

$-(p \cdot \psi')' + q \cdot \psi = \lambda \cdot w \cdot \psi$

mit integrierbaren reellen Funktionen w(x) > 0,p(x) ^{− 1} > 0,q(x), zusammen mit Randbedingungen der Form

$\cos(\alpha) \psi(a) + \sin(\alpha) p(a)\psi'(a) = 0,\quad \cos(\beta) \psi(b) + \sin(\beta) p(b)\psi'(b) = 0, \qquad \alpha,\beta\in [0,\pi),$

nennt man ein reguläres Sturm-Liouville-Problem über dem Intervall [a,b], wenn dieses Intervall endlich ist.

Im Fall ψ(a) = ψ(b) = 0 spricht man von Dirichlet-Randbedingungen und im Fall ψ'(a) = ψ'(b) = 0 von Neumann-Randbedingungen.

Für das reguläre Sturm-Liouville-Problem gilt, dass es eine abzählbare Folge von reellen Eigenwerten gibt, die gegen $+\infty$ divergiert:

$\lambda_1 < \lambda_2 < \lambda_3 < \cdots < \lambda_n < \cdots \to \infty.$

Die Eigenwerte verhalten sich asymptotisch wie (Weyl Asymptotik):

$\lambda_n= \pi^2\left(\int_a^b \sqrt{\frac{w(x)}{p(x)}} dx\right)^{-2} n^2 + O(n).$

Die zugehörigen Eigenfunktionen ψ_n bilden eine Orthonormalbasis im Hilbertraum L²([a,b],w(x)dx) der bezüglich der Gewichtsfunktion w quadratintegrierbaren Funktionen.

Der Sturm'sche Oszillationssatz besagt, dass im Fall von Dirichlet-Randbedingungen die n-te Eigenfunktion ψ_n genau n − 1 Nullstellen im Intervall (a,b) besitzt. Insbesondere kann die tiefste Eigenfunktion positiv gewählt werden.

Beispiel

Ein einfaches Beispiel ist die Differentialgleichung

− ψ'' = λψ

auf dem Intervall [0,π], zusammen mit den Dirichlet-Randbedingungen

ψ(0) = ψ(π) = 0.

Die Eigenwerte sind

λ_n = n²,

und die normierten Eigenfunktionen sind

$\psi_n(x) = \sqrt{\frac{2}{\pi}} \sin(n\, x).$

Die zugehörige Eigenfunktionsentwicklung ist die Fourierreihe.

Mathematische Theorie

Der geeignete mathematische Rahmen ist der Hilbertraum L²([a,b];w(x)dx) mit dem Skalarprodukt

$\langle f, g\rangle := \int_{a}^{b} \overline{f(x)} g(x)w(x)\,{\rm d}x$.

In diesem Raum ist L ein selbstadjungierter Operator, wenn er auf der Menge der (im Sinne der schwachen Ableitung) differenzierbaren Funktionen, die die Randbedingungen erfüllen, definiert wird:

$\mathfrak{D}(L) = \{ f \in L^2([a,b]; w(x){\rm d}x) : f,pf' \in AC[a,b], \, L f \in L^2([a,b]; w(x){\rm d}x), \, \cos(\alpha) f(a) + \sin(\alpha) p(a) f'(a) = \cos(\beta) f(b) + \sin(\beta) p(b) f'(b) = 0\}.$

Hierbei bezeichnet AC[a,b] die Menge der auf [a,b] absolut stetigen Funktionen. Da L ein unbeschränkter Operator ist, betrachtet man die Resolvente

$(L - z)^{-1}, \qquad z \in\mathbb{C}$,

wobei z kein Eigenwert sein darf. Es stellt sich heraus, dass die Resolvente ein Integraloperator mit stetigem Kern (die Green'sche Funktion des Randwertproblems) ist. Somit ist die Resolvente ein kompakter Operator, und die Existenz einer abzählbaren Folge von Eigenfunktionen folgt aus dem Spektralsatz für kompakte Operatoren.

Der Zusammenhang zwischen den Eigenwerten von L und der Resolvente folgt, da (L − z) ^{− 1}ψ = αψ äquivalent ist zu Lψ = λψ mit λ = (z + α ^{− 1}).

Singuläre Sturm-Liouville-Probleme

Sind obige Bedingungen nicht erfüllt, so spricht man von einem singulären Sturm-Liouville-Problem. Das Spektrum besteht dann im Allgemeinen nicht mehr nur aus Eigenwerten und besitzt auch einen kontinuierlichen Anteil. Es gibt weiterhin verallgemeinerte Eigenfunktionen, und die zugehörige Eigenfunktionsentwicklung ist eine Integraltransformation (vergleiche Fouriertransformation anstelle von Fourierreihe).

Wechseln p oder w das Vorzeichen, so spricht man von einem indefiniten Sturm-Liouville-Problem

Anwendung

• Der Fall p(x) = w(x) = 1 entspricht der eindimensionalen zeitunabhängigen Schrödingergleichung.
• Der Separationsansatz zur Lösung partieller Differentialgleichungen führt auf Sturm-Liouville-Probleme.

Weblinks

Kai Gehrs: Sturm-Liouville-Probleme. (PDF; 152 kB)

Literatur

• Gerald Teschl: Ordinary Differential Equations and Dynamical Systems. American Mathematical Society, Providence 2011 (freie Onlineversion).
• Wolfgang Walter: Gewöhnliche Differentialgleichungen. 7. Auflage. Springer, Berlin 2000, ISBN 3-540-67642-2.
• Joachim Weidmann: Lineare Operatoren in Hilberträumen / Teil 2. Anwendungen. Teubner, Stuttgart; Leipzig; Wiesbaden 2003, ISBN 3-519-02237-0.