In der Analysis wird die Sattelpunktsnäherung verwendet, um Integrale der Form
![{\displaystyle I=\lim _{N\to \infty }\int \limits _{-\infty }^{\infty }e^{-Nf(x)}\;\mathrm {d} x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1594531afd0fbc5288280e71fbaa3af2ce569c13)
näherungsweise zu berechnen. Die Methode stammt von Pierre Simon de Laplace (1774) und wird manchmal nach ihm benannt. Sie ist Teil der asymptotischen Analyse.
Falls die Funktion
analytisch ist und ein globales Minimum bei
besitzt, so erhält man:
![{\displaystyle I=\lim _{N\to \infty }e^{-Nf(x_{0})}{\sqrt {\frac {2\pi }{Nf''(x_{0})}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6fddc35a677aef115549e2737c2e24668a4069f9)
mit
.
Die zweite Ableitung ist positiv, da hier ein Minimum vorliegt. Das Ergebnis gilt asymptotisch, das heißt für
gegen Unendlich.
Dabei können auch endliche Integrationsgrenzen
vorliegen.
Die Verallgemeinerung der Sattelpunktnäherung in die komplexe Zahlenebene wird auch Sattelpunktmethode genannt. Aus ihr erklärt sich die Benennung nach einem Sattelpunkt.
Es kann auch ein anderes Vorzeichen im Exponenten betrachtet werden:
Mit anderem Vorzeichen gilt für
![{\displaystyle I=\lim _{N\to \infty }\int \limits _{-\infty }^{\infty }e^{Nf(x)}\;\mathrm {d} x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95b6e84381db45e24453be0440c38aaad7e5f943)
falls bei
ein globales Maximum vorliegt asymptotisch:
mit
.
Da hier ein Maximum vorliegt, ist die zweite Ableitung negativ.
Begründung
Betrachtet wird der erste Fall (Minimum bei
), die Argumentation im zweiten Fall ist analog.
Für große
wird die Exponentialfunktion außerhalb der Umgebung von
beliebig klein. Deshalb wird
um
in eine Taylorreihe entwickelt:
.
(Wegen des globalen Minimums bei
ist
)
Einsetzen ins Integral liefert
.
Die Größe
ist also der Grenzwert
des Produkts aus
und dem nichtelementaren Integral. Letzteres ist eng mit dem gaußschen Fehlerintegral
bzw. der Gauß-Verteilung verwandt. Das Integral
ist von der Form
, wobei
gilt. Es ist insbesondere
, da das Minimum bei
eine positive zweite Ableitung bedingt, was im Folgenden wichtig sein wird:
Für alle
mit
lässt sich folgende Relation (bspw. über Substitution
) zeigen:
Weiterhin beeinflusst eine Verschiebung von
um die Konstante
nicht den Wert von
, da sich das Integral durch die lineare Substitution
mit
leicht in das obige Integral
überführen lässt, dessen Wert bereits bekannt ist.
Man erhält also mit
:
Somit folgt für
(asymptotisch):
Anwendungen
Die Sattelpunktsnäherung und Sattelpunktmethoden findet verschiedene Anwendungen in der theoretischen Physik, unter anderem in der statistischen Physik im Grenzfall großer Systeme, in der Quantenfeldtheorie bei der Auswertung von Pfadintegralen oder in der Optik.
Eine Anwendung ist die Stirlingformel
![{\displaystyle N!\approx {\sqrt {2\pi N}}N^{N}e^{-N}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1412e81d4e14368c7ca9cf935b1ead543c4796f4)
für große
.
Aus der Definition der Gammafunktion folgt
![{\displaystyle N!=\Gamma (N+1)=\int _{0}^{\infty }e^{-x}x^{N}\,\mathrm {d} x.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ad9ccc8fb8473aec405aa4f0b3851d04f8e6d9c)
Mit der Variablentransformation
(so dass ::
) erhält man:
![{\displaystyle {\begin{aligned}N!&=\int _{0}^{\infty }e^{-Nz}\left(Nz\right)^{N}N\,\mathrm {d} z\\&=N^{N+1}\int _{0}^{\infty }e^{-Nz}z^{N}\,\mathrm {d} z\\&=N^{N+1}\int _{0}^{\infty }e^{-Nz}e^{N\ln z}\,\mathrm {d} z\\&=N^{N+1}\int _{0}^{\infty }e^{N(\ln z-z)}\,\mathrm {d} z.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf18dfccf08504f584cc99810683321cbcdd9aa4)
Nun kann man die Sattelpunktnäherung in der zweiten Form (für Maxima) anwenden mit
![{\displaystyle f\left(z\right)=\ln {z}-z}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe9182e53d3d3a89aea0275106e38634189b7e7d)
mit den Ableitungen
![{\displaystyle f'(z)={\frac {1}{z}}-1,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/00112267c3ad76ff7d06aff310a11eec60876ab2)
![{\displaystyle f''(z)=-{\frac {1}{z^{2}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2c5391be065933b63dead23dca21fab1074180af)
Das Maximum von f liegt bei
mit dem Wert der zweiten Ableitung −1. Man erhält mit der Sattelpunktnäherung:
![{\displaystyle N!\approx N^{N+1}{\sqrt {\frac {2\pi }{N}}}e^{-N}={\sqrt {2\pi N}}N^{N}e^{-N}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2ca8e7def238dc434a58a2c191b3d151c34b08bb)
Verallgemeinerung
Die Sattelpunktnäherung wird bei Betrachtung im Komplexen in der Methode des steilsten Abstiegs (englisch: Method of steepest descent) bzw. der Methode der stationären Phase (englisch: Method of stationary phase) verallgemeinert (allgemein Sattelpunktmethode). Ziel ist die asymptotische Auswertung von geschlossenen Wegintegralen in der komplexen Zahlenebene (
)
![{\displaystyle \displaystyle \int _{C}f(z)e^{\lambda g(z)}\mathrm {d} z}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c6e5d6f7fca50bcffbe8f5518ffd77fe7b11069)
für große reelle
. Dabei deformiert man im Komplexen den Integrationsweg so, dass ein stationärer Punkt (Nullstelle der ersten Ableitung von g)
von
auf dem Integrationsweg liegt und geht dann ähnlich wie oben vor (unter zusätzlicher Anwendung des Cauchyschen Integralsatzes). In der Version der Methode des steilsten Abstiegs legt man den Integrationsweg bei
so, dass der Realteil u von g dort ein Maximum hat. Da der Realteil u von g eine harmonische Funktion ist, können
und
nicht dasselbe Vorzeichen haben: es liegt ein Sattelpunkt vor und man legt den Integrationsweg längs des Wegs des „steilsten Abstiegs“. Daher der Name der Methode.
Bei der Methode der stationären Phase werden speziell Integrale betrachtet, bei denen der Exponent der Exponentialfunktion längs des Weges imaginär ist:
![{\displaystyle \displaystyle \int _{C}f(z)e^{i\lambda u(z)}\mathrm {d} z}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a8b591e36a57df9c3cc6a390c7bdc5d515d2c246)
Mit einer reellen Funktion u und großem
.
Die Methode wurde zuerst von Peter Debye 1909 zur Abschätzung von Besselfunktionen veröffentlicht, aber auch schon von Bernhard Riemann benutzt.[1]
Literatur
- Jon Matthews, R. L. Walker: Mathematical Methods of Physics, Addison-Wesley, Kapitel 3.6 (Saddle Point Methods)
- Philip M. Morse, Herman Feshbach: Methods of Theoretical Physics, Band 1, McGraw Hill, 1953, Kapitel 4.6
- Arthur Erdélyi: Asymptotic Expansions, Dover 1956
- Nicolaas Govert de Bruijn: Asymptotic Methods in Analysis, Wiley 1961
Einzelnachweise
- ↑ Sattelpunktmethode. In: Guido Walz (Hrsg.): Lexikon der Mathematik. 1. Auflage. Spektrum Akademischer Verlag, Mannheim/Heidelberg 2000, ISBN 3-8274-0439-8.