Padé approximant explained

In mathematics, a Padé approximant is the "best" approximation of a function near a specific point by a rational function of given order. Under this technique, the approximant's power series agrees with the power series of the function it is approximating. The technique was developed around 1890 by Henri Padé, but goes back to Georg Frobenius, who introduced the idea and investigated the features of rational approximations of power series.

The Padé approximant often gives better approximation of the function than truncating its Taylor series, and it may still work where the Taylor series does not converge. For these reasons Padé approximants are used extensively in computer calculations. They have also been used as auxiliary functions in Diophantine approximation and transcendental number theory, though for sharp results ad hoc methods—in some sense inspired by the Padé theory—typically replace them. Since a Padé approximant is a rational function, an artificial singular point may occur as an approximation, but this can be avoided by Borel–Padé analysis.

The reason the Padé approximant tends to be a better approximation than a truncating Taylor series is clear from the viewpoint of the multi-point summation method. Since there are many cases in which the asymptotic expansion at infinity becomes 0 or a constant, it can be interpreted as the "incomplete two-point Padé approximation", in which the ordinary Padé approximation improves the method truncating a Taylor series.

Definition

Given a function and two integers and, the Padé approximant of order is the rational function

R(x) = \frac = \frac,which agrees with to the highest possible order, which amounts to\begin f(0) &= R(0), \\ f'(0) &= R'(0), \\ f(0) &= R(0), \\ &\mathrel \\ f^(0) &= R^(0).\end

Equivalently, if

R(x)

is expanded in a Maclaurin series (Taylor series at 0), its first

m+n

terms would equal the first

m+n

terms of

f(x)

, and thusf(x) - R(x) = c_ x^ + c_ x^ + \cdots

When it exists, the Padé approximant is unique as a formal power series for the given m and n.[1]

The Padé approximant defined above is also denoted as[m/n]_f(x).

Computation

For given, Padé approximants can be computed by Wynn's epsilon algorithm[2] and also other sequence transformations[3] from the partial sumsT_N(x) = c_0 + c_1 x + c_2 x^2 + \cdots + c_N x^N of the Taylor series of, i.e., we havec_k = \frac. can also be a formal power series, and, hence, Padé approximants can also be applied to the summation of divergent series.

One way to compute a Padé approximant is via the extended Euclidean algorithm for the polynomial greatest common divisor.[4] The relationR(x) = P(x)/Q(x) = T_(x) \bmod x^is equivalent to the existence of some factor

K(x)

such thatP(x) = Q(x)T_(x) + K(x)x^,which can be interpreted as the Bézout identity of one step in the computation of the extended greatest common divisor of the polynomials

Tm+n(x)

and

xm+n+1

.

Recall that, to compute the greatest common divisor of two polynomials p and q, one computes via long division the remainder sequencer_0=p,\;r_1=q,\quad r_= q_k r_k +r_, with

\degrk+1<\degrk

, until

rk+1=0

. For the Bézout identities of the extended greatest common divisor one computes simultaneously the two polynomial sequencesu_0=1,\;v_0=0,\quad u_1=0,\;v_1=1,\quad u_=u_-q_k u_k,\;v_=v_-q_kv_kto obtain in each step the Bézout identityr_k(x) = u_k(x) p(x) + v_k(x) q(x).

For the approximant, one thus carries out the extended Euclidean algorithm forr_0=x^,\;r_1=T_(x)and stops it at the last instant that

vk

has degree or smaller.

Then the polynomials

P=rk,Q=vk

give the Padé approximant. If one were to compute all steps of the extended greatest common divisor computation, one would obtain an anti-diagonal of the Padé table.

Riemann–Padé zeta function

To study the resummation of a divergent series, say \sum_^\infty f(z), it can be useful to introduce the Padé or simply rational zeta function as \zeta_R(s) = \sum_^\infty \frac, where R(x) = [m/n]_f(x)is the Padé approximation of order of the function . The zeta regularization value at is taken to be the sum of the divergent series.

The functional equation for this Padé zeta function is \sum_^n a_j \zeta_R(s-j) = \sum_^m b_j \zeta_0(s-j), where and are the coefficients in the Padé approximation. The subscript '0' means that the Padé is of order [0/0] and hence, we have the Riemann zeta function.

DLog Padé method

Padé approximants can be used to extract critical points and exponents of functions.[5] [6] In thermodynamics, if a function behaves in a non-analytic way near a point like

f(x)\sim|x-r|p

, one calls a critical point and the associated critical exponent of . If sufficient terms of the series expansion of are known, one can approximately extract the critical points and the critical exponents from respectively the poles and residues of the Padé approximants

[n/n+1]g(x)

, where

g=f'/f

.

Generalizations

A Padé approximant approximates a function in one variable. An approximant in two variables is called a Chisholm approximant (after J. S. R. Chisholm),[7] in multiple variables a Canterbury approximant (after Graves-Morris at the University of Kent).[8]

Two-points Padé approximant

The conventional Padé approximation is determined to reproduce the Maclaurin expansion up to a given order. Therefore, the approximation at the value apart from the expansion point may be poor. This is avoided by the 2-point Padé approximation, which is a type of multipoint summation method.[9] At

x=0

, consider a case that a function

f(x)

which is expressed by asymptotic behavior

f0(x)

:f \sim f_0(x) + o\big(f_0(x)\big), \quad x \to 0, and at

x\toinfty

, additional asymptotic behavior

finfty(x)

:f(x) \sim f_\infty(x) + o\big(f_\infty(x)\big), \quad x \to \infty.

By selecting the major behavior of

f0(x),finfty(x)

, approximate functions

F(x)

such that simultaneously reproduce asymptotic behavior by developing the Padé approximation can be found in various cases. As a result, at the point

x\toinfty

, where the accuracy of the approximation may be the worst in the ordinary Padé approximation, good accuracy of the 2-point Padé approximant is guaranteed. Therefore, the 2-point Padé approximant can be a method that gives a good approximation globally for

x=0\siminfty

.

In cases where

f0(x),finfty(x)

are expressed by polynomials or series of negative powers, exponential function, logarithmic function or

xlnx

, we can apply 2-point Padé approximant to

f(x)

. There is a method of using this to give an approximate solution of a differential equation with high accuracy. Also, for the nontrivial zeros of the Riemann zeta function, the first nontrivial zero can be estimated with some accuracy from the asymptotic behavior on the real axis.

Multi-point Padé approximant

A further extension of the 2-point Padé approximant is the multi-point Padé approximant. This method treats singularity points

x=xj(j=1,2,3,...,N)

of a function

f(x)

which is to be approximated. Consider the cases when singularities of a function are expressed with index

nj

byf(x) \sim \frac, \quad x \to x_j.

Besides the 2-point Padé approximant, which includes information at

x=0,x\toinfty

, this method approximates to reduce the property of diverging at

x\simxj

. As a result, since the information of the peculiarity of the function is captured, the approximation of a function

f(x)

can be performed with higher accuracy.

Examples

: \sin(x) \approx \frac

Notes and References

  1. Web site: Padé Approximant . Wolfram MathWorld.
  2. Theorem 1 in On the Convergence and Stability of the Epsilon Algorithm. Peter . Wynn . SIAM Journal on Numerical Analysis . 3 . Mar 1966 . 91–122. 2949688. 1 . 1966SJNA....3...91W . 10.1137/0703007.
  3. Extrapolation algorithms and Padé approximations . C. . Brezenski . Applied Numerical Mathematics . 20 . 1996 . 299–318 . 10.1016/0168-9274(95)00110-7 . 3 . 10.1.1.20.9528.
  4. Book: Polynomial and Matrix computations - Volume 1. Fundamental Algorithms. Dario . Bini. Victor . Pan . Birkhäuser. Progress in Theoretical Computer Science . 1994. 978-0-8176-3786-6 . Problem 5.2b and Algorithm 5.2 (p. 46) .
  5. Adler . Joan . Series expansions . Computers in Physics . 1994 . 8 . 3 . 287 . 10.1063/1.168493 . 1994ComPh...8..287A . free .
  6. Baker . G. A. Jr. . Padé approximant . Scholarpedia . 2012 . 7 . 6 . 9756 . 10.4249/scholarpedia.9756 . 2012SchpJ...7.9756B . free .
  7. Chisholm. J. S. R.. 1973. Rational approximants defined from double power series. Mathematics of Computation. 27. 124. 841–848. 10.1090/S0025-5718-1973-0382928-6. 0025-5718. free.
  8. Graves-Morris. P.R. . Roberts. D.E. . 1975 . Calculation of Canterbury approximants . Computer Physics Communications . 10. 4 . 234–244 . 10.1016/0010-4655(75)90068-5 . 1975CoPhC..10..234G.
  9. Book: Ueoka, Yoshiki . Introduction to multipoints summation method Modern applied mathematics that connects here and the infinite beyond: From Taylor expansion to application of differential equations.