Padé-Approximation

Die Padé-Approximation bezeichnet in der Mathematik die beste Approximation einer Funktion durch rationale Funktionen.

Die Padé-Approximation ist benannt nach dem französischen Mathematiker Henri Padé, der sie 1892 bekannt machte, wobei allerdings der deutsche Mathematiker Georg Frobenius bereits 1881 seine diesbezüglichen Untersuchungen über die rationale Approximation von Potenzreihen veröffentlichte.

Die Padé-Approximation führt oft zu besseren Ergebnissen als die Approximation mittels Taylorreihen. Manchmal erhält man auch dann Approximationen, wenn die Taylorreihe nicht konvergiert. Daher wird sie häufig in Computerberechnungen verwendet. Auch im Gebiet der Diophantischen Approximation ist sie nützlich.

Definition

Sei $f$ eine Funktion und $m\geq 0$, $n\geq 1$ natürliche Zahlen, dann ist die Padé-Approximation der Ordnung $[m/n]$ die rationale Funktion

$R(x)={\frac {\sum _{j=0}^{m}a_{j}x^{j}}{1+\sum _{k=1}^{n}b_{k}x^{k}}}={\frac {a_{0}+a_{1}x+a_{2}x^{2}+\cdots +a_{m}x^{m}}{1+b_{1}x+b_{2}x^{2}+\cdots +b_{n}x^{n}}}$,

welche mit $f(x)$ in der höchstmöglichen Ordnung übereinstimmt, woraus folgt:

${\begin{array}{rcl}f(0)&=&R(0)\\f'(0)&=&R'(0)\\f''(0)&=&R''(0)\\&\vdots &\\f^{(m+n)}(0)&=&R^{(m+n)}(0).\end{array}}$

Eine äquivalente Definition lautet: Entwickelt man $R(x)$ in eine Maclaurinsche Reihe, d. h. in eine Taylorreihe um den Punkt 0, dann stimmen die ersten $m+n$ Terme von $R(x)$ und $f(x)$ überein. Daraus folgt für den Approximationsfehler

$f(x)-R(x)=c_{m+n+1}x^{m+n+1}+c_{m+n+2}x^{m+n+2}+\cdots$

Für jedes vorgegebene $m$ und $n$ ist die Padé-Approximation eindeutig, d. h. die Koeffizienten $a_{0},a_{1},\dots ,a_{m},b_{1},\dots ,b_{n}$ sind eindeutig aus den Koeffizienten $c_{0}\dots c_{m+n}$ der taylorartigen Reihe bestimmbar.

Im Nenner von $R(x)$ wurde der Anfangsterm ohne Beschränkung der Allgemeinheit $b_{0}=1$ gewählt. Andernfalls erhält man durch geeignetes Kürzen die genannte Form.

Die Padé-Approximation wird auch dargestellt als

$[m/n]_{f}(x).\,$

Berechnung

Zu einem gegebenen $x$ kann man die Padé-Approximation nach dem sogenannten „Epsilon-Verfahren“ des belgischen Mathematikers Peter Wynn, oder auch anderer Folgentransformationen berechnen. Dabei verwendet man die Teilsummen

$T_{N}(x)=c_{0}+c_{1}x+c_{2}x^{2}+\cdots +c_{N}x^{N}$

der Taylorreihe von $f$; die $c_{k}$ sind also gemäß

$c_{k}={\frac {f^{(k)}(0)}{k!}}$

durch $f$ bestimmt.

Bei der Funktion $f$ kann es sich auch um eine formale Potenzreihe handeln, so dass Padé-Approximationen auch auf die Summierung divergenter Reihen angewandt werden können.

Zur Berechnung der Padé-Approximation kann man den erweiterten euklidischen Algorithmus für den größten gemeinsamen Polynomteiler anwenden. Die Beziehung

$R(x)=P(x)/Q(x)=T_{m+n}(x){\text{ mod }}x^{m+n+1}$

ist äquivalent zur Existenz eines Faktors $K(x)$ derart, dass

$P(x)=Q(x)T_{m+n}(x)+K(x)x^{m+n+1}$.

Dies lässt sich als die Bézout-Gleichung eines Schrittes der Berechnung des größten gemeinsamen Polynomteilers interpretieren:

$T_{m+n}(x)$ und $x^{m+n+1}$.

Für die $[m/n]$-Approximation wendet man den erweiterten euklidischen Algorithmus an für

$r_{0}=x^{m+n+1},\;r_{1}=T_{m+n}(x)$

und stoppt, wenn $v_{k}$ vom Grade kleiner gleich $n$ ist. Dann stellt der Quotient der Polynome $P=r_{k},\;Q=v_{k}$ die $[m/n]$-Padé-Approximation $R=P/Q$ dar.

Riemann–Padé Zeta-Funktion

Zur Untersuchung von divergenten Reihen, etwa

$\sum _{z=1}^{\infty }f(z)$

kann es hilfreich sein, die Padé- oder rationale Zeta-Funktion einzuführen:

$\zeta _{R}(s)=\sum _{z=1}^{\infty }{\frac {R(z)}{z^{s}}}$,

wobei

$R(x)=[m/n]_{f}(x)\,$

die Padé-Approximation der Ordnung $[m/n]$ der Funktion $f$ ist. Der Wert für $s=0$ ist die Summe der divergenten Reihen. Die Funktionsgleichung für diese Zeta-Funktion lautet:

$\sum _{j=0}^{n}a_{j}\zeta _{R}(s-j)=\sum _{j=0}^{m}b_{j}\zeta _{0}(s-j),$

wobei $a_{j}$ und $b_{j}$ die Koeffizienten der Padé-Approximation sind. Der Index 0 steht für die Padé-Approximation der Ordnung [0/0] und ergibt so die Riemannsche ζ-Funktion.

DLog-Padé Methode

Mit Padé-Approximationen können kritische Punkte und Exponenten einer Funktion ermittelt werden. In der Thermodynamik heißt $x=r$ kritischer Punkt und $p$ der zugehörige kritische Exponent von $f$, wenn sich die Funktion $f$ in der Nähe eines Punktes $x=r$ wie $f(x)\sim \left|x-r\right|^{p}$ nicht-analytisch verhält. Sind hinreichend viele Terme der Reihenentwicklung von $f$ bekannt, dann ergeben sich näherungsweise die kritischen Punkte und die kritischen Exponenten aus den Polen und Residuen der Padé-Approximationen $\left[n/n+1\right]_{g}\left(x\right)$ mit $g={\frac {f'}{f}}$.

Verallgemeinerungen

Eine Padé-Approximation approximiert eine Funktion in einer Variablen. Eine Approximation in zwei Variablen heißt Chisholm-Approximation, in mehr als zwei Variablen Canterbury-Approximation (benannt nach Graves-Morris an der University of Kent).

Literatur

• G. A. Baker, Jr., P. Graves-Morris: Padé Approximants. Cambridge U.P., 1996, ISBN 0-521-45007-1.
• C. Brezinski, M. Redivo Zaglia: Extrapolation Methods. Theory and Practice. North-Holland, 1991, ISBN 0-444-88814-4.
• W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. Flannery: Numerical Recipes. The Art of Scientific Computing. 3. Auflage. Cambridge University Press, New York 2007, ISBN 978-0-521-88068-8, Section 5.12 Padé Approximants. (apps.nrbook.com)
• W. B. Gragg: The Pade Table and Its Relation to Certain Algorithms of Numerical Analysis. In: SIAM Review. Vol. 14, No. 1, 1972, S. 1–62.
• P. Wynn: Upon systems of recursions which obtain among the quotients of the Padé table. In: Numerische Mathematik. 8 (3), 1966, S. 264–269.

Weblinks

• Eric W Weisstein: Padé Approximant von MathWorld (abgerufen am 3. Juni 2014)
• Kapitel 3.5 (abgerufen am 4. Juni 2014)
• Beispiele von Padé-Approximationen (engl.) (abgerufen am 14. Februar 2016)
• Padé-Approximation für sin(x) (engl.) (Memento vom 1. März 2014 im Internet Archive) (abgerufen am 14. Februar 2016)