Orthogonale Polynome
Da die Laguerre-Polynome für
und/oder
divergent sind, bilden sie keinen Prähilbertraum und keinen Hilbertraum. Deshalb wird eine Gewichtsfunktion eingeführt, welche die Lösung der Differentialgleichung ungeändert lässt und welche dafür sorgt, dass die Laguerre-Polynome quadratintegrierbar werden. Unter diesen Voraussetzungen bilden die Eigenfunktionen
eine Orthonormalbasis im Hilbertraum
der quadratintegrierbaren Funktionen mit der Gewichtsfunktion
. Demzufolge gilt
![{\displaystyle \langle L_{n},L_{m}\rangle =\int _{0}^{\infty }\mathrm {e} ^{-x}L_{n}(x)L_{m}(x)\mathrm {d} x=\delta _{nm}.}](//wikimedia.org/api/rest_v1/media/math/render/svg/3fa64d401d6bbddf38b8fc411849d8ce271fc864)
Hierbei bedeutet
das Kronecker-Delta.
- Beweis
Teil 1: Zunächst wird gezeigt, dass die Laguerre-Polynome mit dem Gewicht
orthogonal sind, für
gilt demnach ![{\displaystyle \langle L_{n},L_{m}\rangle =\int _{0}^{\infty }\mathrm {e} ^{-x}L_{n}(x)L_{m}(x)\mathrm {d} x=0.}](//wikimedia.org/api/rest_v1/media/math/render/svg/63c31855263fe02d76b3eb1c8d3d9e3d51a6dea6)
Mit dem Sturm-Liouville-Operator
ergeben sich für die Laguerre-Polynome
folgende Ausgangsgleichungen:
- (1)
![{\displaystyle \quad {\mathcal {L}}L_{n}=-\mathrm {e} ^{x}{\frac {\mathrm {d} }{\mathrm {d} x}}\left(x{\mathrm {e} }^{-x}{\frac {\mathrm {d} L_{n}}{\mathrm {d} x}}\right)=nL_{n}}](//wikimedia.org/api/rest_v1/media/math/render/svg/509989ba5368f584571e0fec4c1b00562a2a7022)
und
- (2)
.
Wird Gleichung (1) von links mit
multipliziert und von Gleichung (2), welche ebenfalls von links mit
multipliziert wird, subtrahiert, so ergeben sich die beiden Gleichungen:
- (3)
![{\displaystyle \quad L_{n}{\mathcal {L}}L_{m}-L_{m}{\mathcal {L}}L_{n}=-L_{n}\mathrm {e} ^{x}{\frac {\mathrm {d} }{\mathrm {d} x}}\left(x{\mathrm {e} }^{-x}{\frac {\mathrm {d} L_{m}}{\mathrm {d} x}}\right)+L_{m}\mathrm {e} ^{x}{\frac {\mathrm {d} }{\mathrm {d} x}}\left(x{\mathrm {e} }^{-x}{\frac {\mathrm {d} L_{n}}{\mathrm {d} x}}\right)}](//wikimedia.org/api/rest_v1/media/math/render/svg/aa80b586eece46edc9e55ddae19359581acf39b4)
und
- (4)
.
Zunächst wird Gleichung (3) zusammengefasst. Mit der Produktregel für Ableitungen, der Term
bleibt hierbei unberücksichtigt, ergeben sich folgende Darstellungen
![{\displaystyle \textstyle L_{n}{\frac {\mathrm {d} }{\mathrm {d} x}}\left(x{\mathrm {e} }^{-x}{\frac {\mathrm {d} L_{m}}{\mathrm {d} x}}\right)={\frac {\mathrm {d} }{\mathrm {d} x}}\left(x{\mathrm {e} }^{-x}L_{n}{\frac {\mathrm {d} L_{m}}{\mathrm {d} x}}\right)-\left(x{\mathrm {e} }^{-x}{\frac {\mathrm {d} L_{m}}{\mathrm {d} x}}\right){\frac {\mathrm {d} L_{n}}{\mathrm {d} x}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/f496422df663404c100bc9f360e5faaa5a867b68)
und
.
Auf diese Weise wird erkennbar, dass der zweite Term in beiden Ableitungen gleich ist und bei der Differenzenbildung verschwindet, also:
- (5)
![{\displaystyle {\begin{aligned}\quad L_{n}{\mathcal {L}}L_{m}-L_{m}{\mathcal {L}}L_{n}&=-\mathrm {e} ^{x}{\frac {\mathrm {d} }{\mathrm {d} x}}\left(x{\mathrm {e} }^{-x}L_{n}{\frac {\mathrm {d} L_{m}}{\mathrm {d} x}}\right)+\mathrm {e} ^{x}{\frac {\mathrm {d} }{\mathrm {d} x}}\left(x{\mathrm {e} }^{-x}L_{m}{\frac {\mathrm {d} L_{n}}{\mathrm {d} x}}\right)\\\\&=-\mathrm {e} ^{x}{\frac {\mathrm {d} }{\mathrm {d} x}}{\bigg (}x{\mathrm {e} }^{-x}\left(L_{n}{\frac {\mathrm {d} L_{m}}{\mathrm {d} x}}-L_{m}{\frac {\mathrm {d} L_{n}}{\mathrm {d} x}}\right){\bigg )}\\\\&=-\mathrm {e} ^{x}{\frac {\mathrm {d} }{\mathrm {d} x}}{\bigg (}x{\mathrm {e} }^{-x}W(L_{n},L_{m}){\bigg )},\\\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/309a1dc014c6adebf117d63147d03415fa59f9f2)
wobei
die Wronski-Determinante der Funktionen
bedeutet.
Zur Berechnung der Wronski-Determinante mittels der Abelschen Identität wird die Differentialgleichung
oder
betrachtet, so dass eine hebbare Singularität bei
entsteht. Die Koeffizientenmatrix des Fundamentalsystems lautet dann
und deren Spur ist
. Somit lautet die Abelsche Identität:
.
Da
und
linear unabhängig sind, ist
– bei genauer Betrachtung ist
– und es ergibt sich folgendes Resultat:
![{\displaystyle {\begin{aligned}W(L_{n},L_{m})(x)&=W(L_{n},L_{m})(0)\exp \left(\int _{0}^{x}{\bigg (}1-{\frac {1}{\xi }}{\bigg )}\mathrm {d} \xi \right)=W(L_{n},L_{m})(0)\exp {\Bigg (}{\bigg [}\xi -\ln \xi {\bigg ]}_{0}^{x}{\Bigg )}\\&=\lim _{\xi \to x}W(L_{n},L_{m})(0)\exp {\Big (}\xi -\ln \xi {\Big )}-\lim _{\xi \to 0}W(L_{n},L_{m})(0)\exp {\Big (}\xi -\ln \xi {\Big )}\\&=\lim _{\xi \to x}W(L_{n},L_{m})(0){\frac {\exp(\xi )}{\exp(\ln \xi )}}-\lim _{\xi \to 0}W(L_{n},L_{m})(0){\frac {\exp(\xi )}{\exp(\ln \xi )}}\\&=\lim _{\xi \to x}W(L_{n},L_{m})(0){\frac {\mathrm {e} ^{\xi }}{\xi }}+\lim _{\xi \to 0}W(L_{n},L_{m})(0){\frac {\mathrm {e} ^{\xi }}{\xi }}\\&=W(L_{n},L_{m})(0){\frac {\mathrm {e} ^{x}}{x}}+\lim _{\xi \to 0}W(L_{n},L_{m})(0){\frac {\mathrm {e} ^{\xi }}{\xi }}+C.\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/564ae29e9266412c3a404c26fcb0e43f8ed32baf)
Die Integrationskonstante wird
gewählt und Gleichung (5) wird mit
multipliziert, so dass folgt:
![{\displaystyle {\begin{aligned}\mathrm {e} ^{-x}{\big (}L_{n}{\mathcal {L}}L_{m}-L_{m}{\mathcal {L}}L_{n}{\big )}&=-{\frac {\mathrm {d} }{\mathrm {d} x}}{\bigg (}x{\mathrm {e} }^{-x}W(L_{n},L_{m})(0){\frac {\mathrm {e} ^{x}}{x}}{\bigg )}\\&=-{\frac {\mathrm {d} }{\mathrm {d} x}}{\bigg (}W(L_{n},L_{m})(0){\bigg )}\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/5a78c3c794389848af9fed5cecfbcbdf2955343a)
Nach Umformen und Trennung der Variablen lautet die Gleichung nun:
![{\displaystyle -\mathrm {e} ^{-x}{\big (}L_{n}{\mathcal {L}}L_{m}-L_{m}{\mathcal {L}}L_{n}{\big )}\mathrm {d} x=\mathrm {d} {\bigg (}W(L_{n},L_{m})(0){\bigg )}}](//wikimedia.org/api/rest_v1/media/math/render/svg/c9070f062761797fbcd83fb28064a12e210c6623)
Auf beiden Seiten der Gleichung stehen nun eindimensionale Pfaffsche Formen und da
eine konstante Funktion ist, gilt
. Für die Berechnung der verbleibenden Pfaffschen Form ist eine geeignete Parametrisierung
zu wählen. Das Integral lautet nun:
.[1]
Demnach verschwindet das Integral längs dem Intervall
, so dass unter Verwendung von Gleichung (4) gilt:
![{\displaystyle 0=(m-n)\int _{0}^{\infty }\mathrm {e} ^{-x}L_{m}L_{n}\mathrm {d} t}](//wikimedia.org/api/rest_v1/media/math/render/svg/e60fe505902dec9afc10978de6b5b55864d001f2)
Diese Bedingung kann nur erfüllt werden, wenn:
.
Teil 2: Im Folgenden wird gezeigt, dass die Laguerre-Polynome mit dem Gewicht
beschränkt sind,[2] für
gilt demnach
, oder abkürzend
.
Für den Beweis wird einerseits die Reihendarstellung
und anderseits die Rodrigues-Formel
benutzt. Es gilt:
.
Für
mit
ergibt sich:
.
Wird nun für
das Laguerre-Polynom zerlegt, so folgt:
![{\displaystyle \langle L_{n},L_{n}\rangle =\sum _{k=0}^{n-1}{\frac {(-1)^{k}}{k!}}\int _{0}^{\infty }x^{k}{\frac {1}{n!}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}\mathrm {d} x+{\frac {(-1)^{n}}{n!}}\int _{0}^{\infty }x^{n}{\frac {1}{n!}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}\mathrm {d} x.}](//wikimedia.org/api/rest_v1/media/math/render/svg/8fff008ddd656b19907138b8c5d8396d07a74ade)
Durch diese Zerlegung wird der Grad des Polynoms in der Summe um 1 reduziert und in der Folge gilt
, wie in Teil 1 gezeigt. Es verbleibt somit lediglich der zweite Term, der mit partieller Integration berechnet wird, also:
![{\displaystyle {\begin{aligned}\langle L_{n},L_{n}\rangle &={\frac {(-1)^{n}}{n!}}\int _{0}^{\infty }x^{n}{\frac {1}{n!}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}\mathrm {d} x\\&={\frac {(-1)^{n}}{n!}}{\bigg [}x^{n}{\frac {1}{n!}}{\frac {\mathrm {d} ^{(n-1)}}{\mathrm {d} x^{(n-1)}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}{\bigg ]}_{0}^{\infty }-n{\frac {(-1)^{n}}{n!}}\int _{0}^{\infty }x^{(n-1)}{\frac {1}{n!}}{\frac {\mathrm {d} ^{(n-1)}}{\mathrm {d} x^{(n-1)}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}\mathrm {d} x\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/8905a83663b8b0b8f61eced14ee8409c30ebb852)
Die Stammfunktion wurde mithilfe der Produktregel berechnet und es ergibt sich im Grenzwert
. Dasselbe Resultat wird im Grenzwert
erhalten. Da dieses Ergebnis für alle
partiellen Integrationen gilt, folgt:
![{\displaystyle {\begin{aligned}\langle L_{n},L_{n}\rangle &=(-1)^{1}n{\frac {(-1)^{n}}{n!}}\int _{0}^{\infty }x^{(n-1)}{\frac {1}{n!}}{\frac {\mathrm {d} ^{(n-1)}}{\mathrm {d} x^{(n-1)}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}\mathrm {d} x\\&=(-1)^{2}n(n-1){\frac {(-1)^{n}}{n!}}\int _{0}^{\infty }x^{(n-2)}{\frac {1}{n!}}{\frac {\mathrm {d} ^{(n-2)}}{\mathrm {d} x^{(n-2)}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}\mathrm {d} x\\&\;\;\vdots \\&=(-1)^{n}n!{\frac {(-1)^{n}}{n!}}\int _{0}^{\infty }x^{(n-n)}{\frac {1}{n!}}{\frac {\mathrm {d} ^{(n-n)}}{\mathrm {d} x^{(n-n)}}}{\bigg (}\mathrm {e} ^{-x}x^{n}{\bigg )}\mathrm {d} x\\&={\frac {(-1)^{2n}}{n!}}\int _{0}^{\infty }\mathrm {e} ^{-x}x^{n}\mathrm {d} x\\&={\frac {1}{n!}}\int _{0}^{\infty }\mathrm {e} ^{-x}x^{n}\mathrm {d} x\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/6a32796d20050ae17a1f5b2caac140c163394578)
Mittels weiterer
-facher partieller Integration oder Integrationstabelle folgt
und somit:
.
Aus Teil 1 und Teil 2 ergibt sich:
![{\displaystyle \langle L_{n},L_{m}\rangle =\int _{0}^{\infty }\mathrm {e} ^{-x}L_{n}(x)L_{m}(x)\mathrm {d} x=\delta _{nm}.}](//wikimedia.org/api/rest_v1/media/math/render/svg/3fa64d401d6bbddf38b8fc411849d8ce271fc864)
![{\displaystyle \Box }](//wikimedia.org/api/rest_v1/media/math/render/svg/029b77f09ebeaf7528fc831fe57848be51f2240b)