# A Numerical Test of Padé Approximation for Some Functions with Singularity

Authors: Hiroaki S. Yamada; Kensuke S. Ikeda
Journal: International Journal of Computational Mathematics
 (2014)

---

## Abstract

The aim of this study is to examine some numerical tests of Padé approximation for some typical functions with singularities such as simple pole, essential singularity, brunch cut, and natural boundary. As pointed out by Baker, it was shown that the simple pole and the essential singularity can be characterized by the poles of the Padé approximation. However, it was not fully clear how the Padé approximation works for the functions with the branch cut or the natural boundary. In the present paper, it is shown that the poles and zeros of the Padé approximated functions are alternately lined along the branch cut if the test function has branch cut, and poles are also distributed around the natural boundary for some lacunary power series and random power series which rigorously have a natural boundary on the unit circle. On the other hand, Froissart doublets due to numerical errors and/or external noise also appear around the unit circle in the Padé approximation. It is also shown that the residue calculus for the Padé approximated functions can be used to confirm the numerical accuracy of the Padé approximation and quasianalyticity of the random power series.


---

## Body

## 1. Introduction

Padé approximation was introduced in mathematics by Hermite and Padé and it has been used in physics for more than 40 years ago [1–3]. In particular, there have been several important examples in physics, to which the Padé approximation was applied, such as summation of the divergent Rayleigh-Schrodinger perturbation series in scattering theory [4], critical phenomena in statistical physics [5–9], denoising from noisy data of time-series [10–14], and detection of singularity of phase space trajectories of Hamiltonian dynamical systems [15–19].Mathematically, the Padé approximation can be used to estimate analyticity of functions. Indeed, the Padé approximation is usually superior to the truncated Taylor expansions when the original function contains any singularity. Let us consider a simple example. The functionf
(
z
)
=
(
1
+
2
z
)
/
(
1
+
z
) has brunch points at z
=
-
1 and z
=
-
1
/
2. The domain of convergence is |
z
|
<
1
/
2. Nevertheless, we can obtain an exact solution √
2
=
1.4142 for z
→
∞ when we apply [
3
∣
3
] diagonal Padé approximation to the function (See Section 3.3 for more details on this example) Although the mathematical validity of the Padé approximation has not been exactly proved yet, the Padé approximation is practically very useful to continue a singular function beyond the domain of convergence.Let us consider a critical phenomenon for Ising model as a simple example in the statistical physics [5, 7, 8]. We assume that at the critical point u
=
u
c the exact magnetic susceptibility χ has a singularity as
(1)
χ
~
u
-
u
c
-
γ
,
where u is a function of temperature and interactions and so on. In this case we sometimes use the logarithmic derivative of χ when we estimate the critical point u
=
u
c and the critical exponent γ as a pole-type singularity as
(2)
d
log
⁡
χ
d
u
=
χ
′
χ
=
γ
(
u
-
u
c
)
,
where χ
′ denotes the derivative with respect to the variable u. In the low temperature expansion for the magnetic susceptibility χ with some coefficients a
n, we assume that the approximated susceptibility χ
N and the logarithmic derivative can be obtained as follows:
(3)
χ
~
χ
N
=
∑
n
=
0
N
‍
a
n
u
n
,
χ
′
χ
~
∑
n
=
0
N
‍
b
n
u
n
.
Then we can estimate the critical point u
=
u
c and the critical exponent γ by the Padé approximation to the truncated expansion χ
N. Note that the singularity on |
u
|
=
u
c will be infinitely differentiable if the coefficients a
n fall off sufficiently rapidly.In general, analytic continuation over the singular point is possible along any other path in complex plane even if the function diverges at the singular point determining the radius of convergence, as seen in the above singularity of the Ising model. Therefore, the Padé approximation is useful to improve the convergence of the power series and approximate the exact solution.Furthermore, the Padé approximation has been used to investigate convergence of Fourier series [2, 3] and the breakdown of KAM curves in complex plane for Hamiltonian map systems, which is described as the analytic domains of Lindstedt series for standard map [15–19].In addition, we can see an interesting example concerning the Padé approximation in noisy data analysis [10–14]. The power series with finite random coefficient “almost always” has a natural boundary on the unit circle in the complex plane [20–24]. In a finite time-series, Froissart has shown that a natural boundary generated by the random time-series is approximated by doublets of poles and zeros (Froissart doublets) of the Padé approximated function, which are surrounding the vicinity of the unit circle. Taking advantage of the characteristic, the Padé approximation has been used in order to remove the noise and extract the true poles associated with damping modes from the observed noisy time-series.The main purpose of the present paper is to investigate whether the Padé approximation is numerically useful for detecting the singularity of some test functions. In particular, we examine the usefulness for the functions with a natural boundary such as a lacunary power series and a random power series [2, 13, 18].The organization of the paper is as follows. In Section2, we give a brief explanation of the Padé approximation and some important reminders in the numerical calculation. In Section 3, we present some numerical results of the Padé approximation for test functions with branch cut, essential singularity. We also try to apply the Padé approximation to an entire function. In Section 4, application of the Padé approximation to some lacunary series which is known to have a natural boundary is given. In Section 5, numerical results of the application to random power series with a natural boundary and some test functions with random noise are also shown. In Section 6, we discuss the residue calculus for the Padé approximated functions to confirm the numerical errors and quasianalyticity of the random power series. In the last section, we give summary and discussion.In AppendixA, the general result for Fibonacci generator used in Section 3.2 is given. In Appendix B, some theorems concerning zeros of polynomials are given. In Appendix C, some mathematical theorems for lacunary series, which are useful in reading the main text, are summarized. In Appendix D, exact Padé approximated functions to some lacunary power series with a natural boundary are given. Residue analysis for quasianalytic functions of Carleman class is given in Appendix E. Furthermore, some theorems concerning the random power series are given in Appendix F.


## 2. Padé Approximation

In this section, we give preliminary instructions to the Padé approximation and some important reminders on the numerical calculation.For a given functionf
(
z
), a truncated Taylor expansion f
[
N
]
(
z
) of order N about zero is given as
(4)
f
(
z
)
~
f
[
N
]
(
z
)
=
∑
n
=
0
N
‍
c
n
z
n
,
where c
n denotes the coefficients of the Taylor expansion. Padé approximation is more accurate approximation for f
(
z
) up to order O
(
z
N
) than the Taylor expansion. The Padé approximation is a rational function, namely, a ratio of two polynomials, which agrees with the highest possible order O
(
z
N
) with a truncated polynomial f
[
N
]
(
z
) as follows:
(5)
f
N
z
=
a
0
+
a
1
z
+
a
2
z
2
+
⋯
a
L
z
L
1
+
b
1
z
+
b
2
z
2
+
⋯
b
M
z
M
≡
P
L
(
z
)
Q
M
(
z
)
≡
f
[
L
∣
M
]
(
z
)
,
where P
L
(
z
) is a polynomial of degree less than or equal to L and Q
M
(
z
) is a polynomial of degree less than or equal to M. Note that b
0
=
1 (normalized) here. A unique approximation can be specified for all choices of M and L such that N
=
L
+
M when it exists. The coefficients {
a
n
}
,
{
b
n
} can be obtained from the condition that the first (
L
+
M
+
1
) terms vanish in the Taylor series in (4). The difference between the Padé approximation and the original function satisfies the following equation:
(6)
f
(
z
)
-
f
[
L
∣
M
]
(
z
)
=
O
(
z
L
+
M
+
1
)
.
In this paper, we sometimes use “[
L
∣
M
] Padé approximation” or “[
L
∣
M
] Padé approximated function” for f
[
L
∣
M
]
(
z
). Solving the problem in (6) is called a linear Teoplitz problem which is generally ill-conditioned. We used full LU decomposition for the Teoplitz matrix in the problem as well as iterative improvement in order to eliminate the ill-posed problem [25, 26]. In addition, hereafter, we use the diagonal Padé approximation, that is, L
=
M, of order M
≤
65 to estimate the singularity of the test functions because of the convergence and limitation due to the round-off errors and other sources of numerical errors. Here, the singularity of the function f
(
z
) is approximated by configuration of the poles and zeros of the [
M
∣
M
] order diagonal Padé approximated function. As mentioned in Introduction, in general, the Padé approximations are useful for representing unknown functions with possible poles. The application of the diagonal Padé approximation is insured for the functions with isolated singular points and rational-type functions. However, it is not fully clarified how the poles and zeros of the Padé approximated function describe brunch cut and natural boundary [27–29].Generally, the magnitude of the residues associated with the spurious poles is much smaller than that with the true poles, and they are close to machine precision. Very recently, Gonnet et al. suggested an efficient algorithm for the Padé approximations [30, 31]. The algorithm detects and eliminates the spurious pole-zero pairs caused by the rounding errors by means of singular value decomposition for the Teoplitz matrix.Before closing this section, we list up some important reminders when we numerically apply the Padé approximation to unknown functions as follows.(1)
More accurate calculation becomes possible by scaling the expansion variablez if there is a simple pole with large magnitude ρ(≫1). That is, we should change the order of radius of convergence into O
(
1
) by the scaling the expansion variable as z
→
z
/
ρ in order to keep the numerical accuracy. This procedure is effective when we apply the Padé approximation to exponentially decaying coefficients {
a
n
} with fluctuation.
(2)
Poles (i.e., roots ofP
M
(
z
)
=
0) and zeros (i.e., roots of Q
M
(
z
)
=
0) are sometimes cancelled (zero-pole ghost pairs). We can remove the effects of the ghost pairs and confirm the singularity of the functions by using the residue analysis of the Padé approximated function.
(3)
Poles and zeros of the Padé approximation to the truncated random power series accumulate around the unit circle as Froissart doublets. It is difficult to distinguish whether the poles of the Padé approximation originated from a natural boundary of the original function or the natural boundary generated by the numerical error and/or noise. Therefore, the numerical accuracy will be important to determine the coefficients of the Padé approximation.
(4)
In general, the denominators of the diagonal Padé approximated functions to the lacunary power series and the random power series become lacunary and random polynomials, respectively. Accordingly, the distribution of the poles and zeros of the approximated functions are similar to the distribution of the zeros corresponding to the original lacunary power polynomial and random power polynomial. In particular, it is well known that the zeros of the random polynomials uniformly distribute about the unit circle (see AppendixB).


## 3. Examples of Padé Approximation for Some Functions

In this section we investigate the configuration of the poles and zeros of the Padé approximation to some test functions with singularity.


### 3.1. Comparison of Padé Approximation with Taylor Expansion

First, we try the Padé approximation to the following test functionf
1
(
z
) with a brunch point at z
=
-
1:
(7)
f
1
(
z
)
=
log
⁡
(
1
+
z
)
z
.
The truncated Taylor expansion of the order N
=
4 around z
=
0 is
(8)
f
1
[
4
]
(
z
)
=
1
-
1
2
z
+
1
3
z
2
-
1
4
z
3
+
1
5
z
4
.
The [
2
∣
2
] Padé approximation is given as
(9)
f
1
[
2
∣
2
]
(
z
)
=
1
+
(
7
/
10
)
z
+
(
1
/
30
)
z
2
1
+
(
6
/
5
)
z
+
(
3
/
10
)
z
2
.Figure1 shows the approximated functions, the truncated Taylor series f
1
[
4
]
(
z
), and the original test function f
1
(
z
). The approximated function f
1
[
4
]
(
z
) by the truncated Taylor expansion converges only within |
z
|
<
1 and deviates from the exact function f
1
(
z
) for |
z
|
>
1. On the other hand, it follows that the Padé approximated function f
1
[
2
∣
2
]
(
z
) well approximates the original function f
1
(
z
) with very high precession even beyond the radius of convergence up to Re
z
=
x
~
10. Therefore, it is found that the divergent power series expansion (Taylor expansion) does still contain information about the original function outside the convergence radius, and rearranging the coefficients of the expansion into the Padé approximation recovers the information. As a result, the conversion from the Taylor form to the Padé form usually accelerates the convergence and often allows good accuracy even outside the radius of convergence of the power series.Figure 1
The Padé approximated functionf
1
[
2
∣
2
]
(
z
), the truncated Taylor series f
1
[
4
]
(
z
), and the original test function f
1
(
z
).


### 3.2. Padé Approximation for Fibonacci Generating Function

Let us consider the Fibonacci generating functionf
F
(
z
), where Fibonacci sequence {
F
n
} is encoded in the power series as the coefficients. The Fibonacci sequence {
F
n
} is given by the following recursion relation:
(10)
F
n
=
F
n
-
1
+
F
n
-
2
(
n
≥
2
)
,
where F
0
=
0 and F
1
=
1. Then the Fibonacci generating function f
F
(
z
) becomes
(11)
f
F
z
=
∑
n
=
0
∞
‍
F
n
z
n
=
z
1
-
z
-
z
2
=
1
5
1
1
-
ϕ
+
z
-
1
5
1
1
-
ϕ
-
z
,
where ϕ
+
≡
(
1
+
5
/
2
)
(
=
1.61803
…
) and ϕ
-
≡
(
1
-
5
/
2
)
(
=
-
0.61803
…
). The generating function has poles at z
=
ϕ
+ and z
=
ϕ
-. Generating functions by more general recursion relation is given in Appendix A.It should be noted that the diagonal Padé approximationf
F
[
N
/
2
∣
N
/
2
]
(
z
) for the truncated Fibonacci generation function f
F
[
N
]
(
z
)
=
∑
n
=
0
N
‍
F
n
z
n has the following form:
(12)
f
F
[
N
/
2
∣
N
/
2
]
(
z
)
=
z
1
-
z
-
z
2
+
O
(
z
N
+
1
)
,
for even number N
≥
2. This means that the diagonal Padé approximation can detect the exact poles of the generating function irrespective of the order.


### 3.3. Examples of Some Test Functions with Pole, Brunch Cut, and Essential Singularity

In this subsection, we use some test functions in applying of the Padé approximation:(13)
f
2
z
=
e
-
z
,
f
3
(
z
)
=
1
+
2
z
1
+
z
,
f
4
(
z
)
=
e
-
z
/
(
1
+
z
)
,
f
5
(
z
)
=
tan
z
4
.
Here, f
2
(
z
) has no singularity for |
z
|
<
∞, f
3
(
z
) has a brunch cut along a line on [
-
1
,
-
1
/
2
], and f
4
(
z
) has an essential singularity at z
=
-
1. f
5
(
z
) has eight poles at points on the unit circle z
=
exp
⁡
{
i
π
(
m
/
4
)
}
(
m
=
0,1
,
2
,
…
,
7
).First, let us apply Padé approximation tof
2
(
z
). In this case the explicit form of the Padé approximated function can be obtained in the following form [33]:
(14)
P
M
(
z
)
=
∑
k
=
0
M
‍
(
2
M
-
k
)
!
M
!
(
2
M
)
!
k
!
(
M
-
k
)
!
-
z
k
,
Q
M
(
z
)
=
∑
k
=
0
M
‍
(
2
M
-
k
)
!
M
!
(
2
M
)
!
k
!
(
M
-
k
)
!
z
k
.
Note that the coefficients of the numerator P
M
(
z
) have always alternative sign, and the zeros and poles of the Padé approximated function are symmetrical to the imaginary axis with each other because P
M
(
z
)
=
Q
M
(
-
z
). Figure 2(a) shows the numerical results in the complex z-plane. All poles are on the left half-plane Re
z
>
0, and all zeros are on the right half-plane Re
z
<
0. The poles and zeros of the Padé approximated functions for the regular function f
2
(
z
) go infinite and disappear as M
→
∞ because the function f
2
(
z
) is an entire function.Distributions of poles (●, ◯) and zeros (+, ×) of the [
M
∣
M
] diagonal Padé approximated functions for some test functions: (a) f
2
(
z
), (b) f
3
(
z
), and (c) f
4
(
z
).
(a)
(b)
(c)In Figure2(b), distribution of the zeros and poles of the Padé approximated function to f
3
(
z
) is shown. The poles and zeros make a line alternately between the two branch points of the f
3
(
z
); z
=
-
1 and z
=
-
1
/
2.Figure2(c) shows the distribution of the zeros and poles of the Padé approximation to f
4
(
z
). The Padé approximation clusters the poles and zeros at the singular point for f
4
(
z
). As M
→
∞ the poles and zeros approach the singular point z
=
-
1 that reflects the essential singularity of the original function f
4
(
z
).Figure3(a) shows the poles and zeros of the Padé approximated function to the test function f
5
(
z
). The poles and zeros are alternatively distributed in eight directions from the origin. It seems that the distance between the poles and/or zeros on the same line becomes small as they approach the location of the true poles. Some “spurious poles” appear around the unit circle with the increase of the order of the Padé approximation, as seen in Figure 3(b), which are irrelevant poles due to insufficient numerical accuracy. It is found that the numerical accuracy of the Padé approximation fails for the higher order of the Padé approximation. We discuss the spurious poles in Sections 4 and 5 again.Distributions of poles (◯) and zeros (×) of the [
M
∣
M
] diagonal Padé approximated functions for the test function f
5
(
z
). (a) The [
50
∣
50
] Padé approximation. (b) The [
75
∣
75
] Padé approximation. The unit circle is drawn to guide the eye.
(a)
(b)


## 3.1. Comparison of Padé Approximation with Taylor Expansion

First, we try the Padé approximation to the following test functionf
1
(
z
) with a brunch point at z
=
-
1:
(7)
f
1
(
z
)
=
log
⁡
(
1
+
z
)
z
.
The truncated Taylor expansion of the order N
=
4 around z
=
0 is
(8)
f
1
[
4
]
(
z
)
=
1
-
1
2
z
+
1
3
z
2
-
1
4
z
3
+
1
5
z
4
.
The [
2
∣
2
] Padé approximation is given as
(9)
f
1
[
2
∣
2
]
(
z
)
=
1
+
(
7
/
10
)
z
+
(
1
/
30
)
z
2
1
+
(
6
/
5
)
z
+
(
3
/
10
)
z
2
.Figure1 shows the approximated functions, the truncated Taylor series f
1
[
4
]
(
z
), and the original test function f
1
(
z
). The approximated function f
1
[
4
]
(
z
) by the truncated Taylor expansion converges only within |
z
|
<
1 and deviates from the exact function f
1
(
z
) for |
z
|
>
1. On the other hand, it follows that the Padé approximated function f
1
[
2
∣
2
]
(
z
) well approximates the original function f
1
(
z
) with very high precession even beyond the radius of convergence up to Re
z
=
x
~
10. Therefore, it is found that the divergent power series expansion (Taylor expansion) does still contain information about the original function outside the convergence radius, and rearranging the coefficients of the expansion into the Padé approximation recovers the information. As a result, the conversion from the Taylor form to the Padé form usually accelerates the convergence and often allows good accuracy even outside the radius of convergence of the power series.Figure 1
The Padé approximated functionf
1
[
2
∣
2
]
(
z
), the truncated Taylor series f
1
[
4
]
(
z
), and the original test function f
1
(
z
).


## 3.2. Padé Approximation for Fibonacci Generating Function

Let us consider the Fibonacci generating functionf
F
(
z
), where Fibonacci sequence {
F
n
} is encoded in the power series as the coefficients. The Fibonacci sequence {
F
n
} is given by the following recursion relation:
(10)
F
n
=
F
n
-
1
+
F
n
-
2
(
n
≥
2
)
,
where F
0
=
0 and F
1
=
1. Then the Fibonacci generating function f
F
(
z
) becomes
(11)
f
F
z
=
∑
n
=
0
∞
‍
F
n
z
n
=
z
1
-
z
-
z
2
=
1
5
1
1
-
ϕ
+
z
-
1
5
1
1
-
ϕ
-
z
,
where ϕ
+
≡
(
1
+
5
/
2
)
(
=
1.61803
…
) and ϕ
-
≡
(
1
-
5
/
2
)
(
=
-
0.61803
…
). The generating function has poles at z
=
ϕ
+ and z
=
ϕ
-. Generating functions by more general recursion relation is given in Appendix A.It should be noted that the diagonal Padé approximationf
F
[
N
/
2
∣
N
/
2
]
(
z
) for the truncated Fibonacci generation function f
F
[
N
]
(
z
)
=
∑
n
=
0
N
‍
F
n
z
n has the following form:
(12)
f
F
[
N
/
2
∣
N
/
2
]
(
z
)
=
z
1
-
z
-
z
2
+
O
(
z
N
+
1
)
,
for even number N
≥
2. This means that the diagonal Padé approximation can detect the exact poles of the generating function irrespective of the order.


## 3.3. Examples of Some Test Functions with Pole, Brunch Cut, and Essential Singularity

In this subsection, we use some test functions in applying of the Padé approximation:(13)
f
2
z
=
e
-
z
,
f
3
(
z
)
=
1
+
2
z
1
+
z
,
f
4
(
z
)
=
e
-
z
/
(
1
+
z
)
,
f
5
(
z
)
=
tan
z
4
.
Here, f
2
(
z
) has no singularity for |
z
|
<
∞, f
3
(
z
) has a brunch cut along a line on [
-
1
,
-
1
/
2
], and f
4
(
z
) has an essential singularity at z
=
-
1. f
5
(
z
) has eight poles at points on the unit circle z
=
exp
⁡
{
i
π
(
m
/
4
)
}
(
m
=
0,1
,
2
,
…
,
7
).First, let us apply Padé approximation tof
2
(
z
). In this case the explicit form of the Padé approximated function can be obtained in the following form [33]:
(14)
P
M
(
z
)
=
∑
k
=
0
M
‍
(
2
M
-
k
)
!
M
!
(
2
M
)
!
k
!
(
M
-
k
)
!
-
z
k
,
Q
M
(
z
)
=
∑
k
=
0
M
‍
(
2
M
-
k
)
!
M
!
(
2
M
)
!
k
!
(
M
-
k
)
!
z
k
.
Note that the coefficients of the numerator P
M
(
z
) have always alternative sign, and the zeros and poles of the Padé approximated function are symmetrical to the imaginary axis with each other because P
M
(
z
)
=
Q
M
(
-
z
). Figure 2(a) shows the numerical results in the complex z-plane. All poles are on the left half-plane Re
z
>
0, and all zeros are on the right half-plane Re
z
<
0. The poles and zeros of the Padé approximated functions for the regular function f
2
(
z
) go infinite and disappear as M
→
∞ because the function f
2
(
z
) is an entire function.Distributions of poles (●, ◯) and zeros (+, ×) of the [
M
∣
M
] diagonal Padé approximated functions for some test functions: (a) f
2
(
z
), (b) f
3
(
z
), and (c) f
4
(
z
).
(a)
(b)
(c)In Figure2(b), distribution of the zeros and poles of the Padé approximated function to f
3
(
z
) is shown. The poles and zeros make a line alternately between the two branch points of the f
3
(
z
); z
=
-
1 and z
=
-
1
/
2.Figure2(c) shows the distribution of the zeros and poles of the Padé approximation to f
4
(
z
). The Padé approximation clusters the poles and zeros at the singular point for f
4
(
z
). As M
→
∞ the poles and zeros approach the singular point z
=
-
1 that reflects the essential singularity of the original function f
4
(
z
).Figure3(a) shows the poles and zeros of the Padé approximated function to the test function f
5
(
z
). The poles and zeros are alternatively distributed in eight directions from the origin. It seems that the distance between the poles and/or zeros on the same line becomes small as they approach the location of the true poles. Some “spurious poles” appear around the unit circle with the increase of the order of the Padé approximation, as seen in Figure 3(b), which are irrelevant poles due to insufficient numerical accuracy. It is found that the numerical accuracy of the Padé approximation fails for the higher order of the Padé approximation. We discuss the spurious poles in Sections 4 and 5 again.Distributions of poles (◯) and zeros (×) of the [
M
∣
M
] diagonal Padé approximated functions for the test function f
5
(
z
). (a) The [
50
∣
50
] Padé approximation. (b) The [
75
∣
75
] Padé approximation. The unit circle is drawn to guide the eye.
(a)
(b)


## 4. Natural Boundary of Lacunary Power Series

In this section, we examine the applicability of the Padé approximation to investigate the analyticity of some well-known test functions with a natural boundary on the unit circle|
z
|
=
1. This will provide a preliminary information about what occurs in the Padé approximated functions for some test functions with a natural boundary. Indeed, we do know only a very few numerical examples which have a natural boundary and allow an exact diagonal Padé approximation.The following are famous lacunary power series with a natural boundary on the unit circle|
z
|
=
1, f
Jac
(
z
)
=
∑
n
=
0
∞
‍
z
2
n, f
Wie
(
z
)
=
∑
n
=
0
∞
‍
z
n
!, and f
Kro
(
z
)
=
∑
n
=
0
∞
‍
z
n
2, where the f
Jac
(
z
), f
Wie
(
z
), and f
Kro
(
z
) are called after Jacobi, Weierstrass, and Kronecker. Some theorems for the lacunary series with a natural boundary are given in Appendix C [20–22].Here, we use polar-formf
r
(
θ
) for the function f
(
z
) by changing the variable, that is, z
=
r
e
i
θ, in order to simply display the functions as
(15)
f
r
(
θ
)
=
f
(
z
=
r
e
i
θ
)
=
∑
n
=
0
∞
‍
c
n
r
e
i
θ
n
.
Then, note that the modulus r works as a convergence factor of the series because it well converges for r
<
1. Typically we take r
=
1 on the unit circle or r
=
0.98 inside the circle in the following numerical calculations.


### 4.1. Example  1: Jacobi Lacunary Series

We try to apply Padé approximation to the functionf
Jac
(
z
) with a natural boundary on the unit circle |
z
|
=
1. The Padé approximated function exactly has the following form:
(16)
f
Jac
2
N
z
~
f
Jac
2
N
-
1
∣
2
N
-
1
z
=
A
Jac
N
(
z
)
1
+
∑
k
=
0
N
-
2
z
2
k
-
z
2
N
-
1
,
where the explicit form of the numerator A
Jac
N
(
z
) is given in Appendix D. Accordingly, the poles of the [
2
N
-
1
∣
2
N
-
1
] Padé approximated function are given by roots of the polynomial
(17)
1
+
∑
k
=
0
N
-
2
‍
z
2
k
-
z
2
N
-
1
=
0
.
This is also just a lacunary polynomial. In Figure 4 the numerical result of the Padé approximation for f
Jac
(
z
) is shown. The poles and zeros are plotted for the [
64
∣
64
] Padé approximation in Figure 4(a). Inside the circle |
z
|
=
1 some cancellations of the ghost pairs appear. The poles and zeros accumulate around |
z
|
=
1 with the increase of order of the Padé approximation. In the case of the M
=
64, the poles accumulate around |
z
|
=
1 with making the zero-pole pairings. Figure 4(b) shows the Padé approximated functions in the polar-form with r
=
1. It well approximates the original function f
Jac
(
z
) when the order of the Padé approximation increases.(a) Distribution of poles (◯) and zeros (×) of the [
64
∣
64
] Padé approximation for the test function f
Jac
(
z
) with a natural boundary on |
z
|
=
1. The unit circle is drawn to guide the eye. (b) The Padé approximated functions f
Jac
[
32
∣
32
]
(
θ
), f
Jac
[
64
∣
64
]
(
θ
) and the exact function f
Jac
(
θ
) in the polar-form with r
=
1.0 (after [32]).
(a)
(b)It is also shown that the complex zeros of the polynomial (17) cluster near unit circle |
z
|
=
1 and distribute uniformly on the circle as F
N
→
∞ by Erdos-Turan-type theorem given in Appendices C and B [34–41].


### 4.2. Example  2: Fibonacci Lacunary Series

As a second example, we would like to apply Padé approximation to the following lacunary series:(18)
f
Fib
(
z
)
=
∑
n
=
0
∞
‍
z
F
n
,
where F
n is nth Fibonacci number. This function also has a natural boundary on |
z
|
=
1. The Padé approximated function exactly has the following form:
(19)
f
Fib
F
N
z
~
f
Fib
F
N
/
2
∣
F
N
/
2
z
=
A
Fib
F
N
(
z
)
1
+
z
F
N
-
4
-
z
F
N
-
2
.
The explicit form of the numerator A
Fib
F
N
(
z
) is given in Appendix D. The poles of the [
F
N
/
2
∣
F
N
/
2
] Padé approximated function are given by zeros of the lacunary polynomial
(20)
1
+
z
F
N
-
4
-
z
F
N
-
2
=
0
.In Figure5 the numerical result of the Padé approximation to f
Fib
(
z
) is shown. The poles and zeros are plotted for the [
55
∣
55
] Padé approximation in Figure 5(a). The poles and zeros accumulate around |
z
|
=
1 with the increase of the order of the Padé approximation. No pole appears inside the unit circle. The original function is also well approximated by the [
56
∣
56
] Padé approximation (see Figure 5(b)).(a) Distribution of poles (◯) and zeros (+) of the [
55
∣
55
] Padé approximated function for the test function f
Fib
(
z
) with a natural boundary on |
z
|
=
1. The unit circle is drawn to guide the eye. (b) The Padé approximated functions f
Fib
[
21
∣
21
]
(
θ
), f
Fib
[
34
∣
34
]
(
θ
), andf
Fib
[
55
∣
55
]
(
θ
) corresponding to N
=
F
9
=
55, N
=
F
10
=
89, and N
=
F
11
=
144, respectively, and the exact function f
Fib
(
θ
) in the polar form with r
=
1.0 (after [32]).
(a)
(b)


## 4.1. Example  1: Jacobi Lacunary Series

We try to apply Padé approximation to the functionf
Jac
(
z
) with a natural boundary on the unit circle |
z
|
=
1. The Padé approximated function exactly has the following form:
(16)
f
Jac
2
N
z
~
f
Jac
2
N
-
1
∣
2
N
-
1
z
=
A
Jac
N
(
z
)
1
+
∑
k
=
0
N
-
2
z
2
k
-
z
2
N
-
1
,
where the explicit form of the numerator A
Jac
N
(
z
) is given in Appendix D. Accordingly, the poles of the [
2
N
-
1
∣
2
N
-
1
] Padé approximated function are given by roots of the polynomial
(17)
1
+
∑
k
=
0
N
-
2
‍
z
2
k
-
z
2
N
-
1
=
0
.
This is also just a lacunary polynomial. In Figure 4 the numerical result of the Padé approximation for f
Jac
(
z
) is shown. The poles and zeros are plotted for the [
64
∣
64
] Padé approximation in Figure 4(a). Inside the circle |
z
|
=
1 some cancellations of the ghost pairs appear. The poles and zeros accumulate around |
z
|
=
1 with the increase of order of the Padé approximation. In the case of the M
=
64, the poles accumulate around |
z
|
=
1 with making the zero-pole pairings. Figure 4(b) shows the Padé approximated functions in the polar-form with r
=
1. It well approximates the original function f
Jac
(
z
) when the order of the Padé approximation increases.(a) Distribution of poles (◯) and zeros (×) of the [
64
∣
64
] Padé approximation for the test function f
Jac
(
z
) with a natural boundary on |
z
|
=
1. The unit circle is drawn to guide the eye. (b) The Padé approximated functions f
Jac
[
32
∣
32
]
(
θ
), f
Jac
[
64
∣
64
]
(
θ
) and the exact function f
Jac
(
θ
) in the polar-form with r
=
1.0 (after [32]).
(a)
(b)It is also shown that the complex zeros of the polynomial (17) cluster near unit circle |
z
|
=
1 and distribute uniformly on the circle as F
N
→
∞ by Erdos-Turan-type theorem given in Appendices C and B [34–41].


## 4.2. Example  2: Fibonacci Lacunary Series

As a second example, we would like to apply Padé approximation to the following lacunary series:(18)
f
Fib
(
z
)
=
∑
n
=
0
∞
‍
z
F
n
,
where F
n is nth Fibonacci number. This function also has a natural boundary on |
z
|
=
1. The Padé approximated function exactly has the following form:
(19)
f
Fib
F
N
z
~
f
Fib
F
N
/
2
∣
F
N
/
2
z
=
A
Fib
F
N
(
z
)
1
+
z
F
N
-
4
-
z
F
N
-
2
.
The explicit form of the numerator A
Fib
F
N
(
z
) is given in Appendix D. The poles of the [
F
N
/
2
∣
F
N
/
2
] Padé approximated function are given by zeros of the lacunary polynomial
(20)
1
+
z
F
N
-
4
-
z
F
N
-
2
=
0
.In Figure5 the numerical result of the Padé approximation to f
Fib
(
z
) is shown. The poles and zeros are plotted for the [
55
∣
55
] Padé approximation in Figure 5(a). The poles and zeros accumulate around |
z
|
=
1 with the increase of the order of the Padé approximation. No pole appears inside the unit circle. The original function is also well approximated by the [
56
∣
56
] Padé approximation (see Figure 5(b)).(a) Distribution of poles (◯) and zeros (+) of the [
55
∣
55
] Padé approximated function for the test function f
Fib
(
z
) with a natural boundary on |
z
|
=
1. The unit circle is drawn to guide the eye. (b) The Padé approximated functions f
Fib
[
21
∣
21
]
(
θ
), f
Fib
[
34
∣
34
]
(
θ
), andf
Fib
[
55
∣
55
]
(
θ
) corresponding to N
=
F
9
=
55, N
=
F
10
=
89, and N
=
F
11
=
144, respectively, and the exact function f
Fib
(
θ
) in the polar form with r
=
1.0 (after [32]).
(a)
(b)


## 5. Natural Boundary of Random Power Series and the Noise Effect on Padé Approximation

In this section, we apply Padé approximation to the random power series with a natural boundary with probability 1 and investigate how the approximation detect the singularity of the series. In addition, we examine the effect of noise on the coefficients of the power expansion for some test functions. Some related theorems for the natural boundary of the function generated by the random power series are given in AppendixF.


### 5.1. Random Power Series and Natural Boundary

Let us consider a random power series:(21)
f
noise
1
(
z
)
=
∑
n
=
0
∞
‍
ϵ
r
n
z
n
.
Here the coefficients, r
0
,
r
1
,
r
2
,
…, are i.i.d. random variables which take a value within r
n
∈
[
0,1
], and ϵ is the strength of the randomness. It is shown that, in general, the random power series has a natural boundary on the unit circle |
z
|
=
1 with probability one. Figure 6(a) shows distribution of poles and zeros of the [
50
∣
50
] Padé approximated function for f
noise
1
(
z
). Some pairs of poles and zeros are perfectly cancelled inside the circle |
z
|
=
1. On the other hand, almost all the poles and zeros of the Padé approximated function assemble around the circle |
z
|
=
1 and not cancelled. The pair of poles and zeros around the circle |
z
|
=
1 is called “Froissart doublets,” and it well corresponds to the natural boundary of f
noise
1
(
z
). The original function is also well approximated by the [
50
∣
50
] Padé approximation (see Figure 6(b)).(a) Distribution of poles (◯) and zeros (×) of the Padé approximated function f
noise
1
[
50
∣
50
]
(
z
) for a random power series f
noise
1
(
z
) with ϵ
=
1. The unit circle is drawn to guide the eye. (b) The Padé approximated function f
noise
1
[
50
∣
50
]
(
θ
) and the exact function f
Fib
(
θ
) in the polar form with r
=
1.0.
(a)
(b)Figure7 shows an example of the coefficients {
c
n
}
=
{
ϵ
r
n
} of the random power series and the coefficients {
a
n
} and {
b
n
} of the [
50
∣
50
] Padé approximated function. The fluctuation of the coefficient {
b
n
} that determines the poles of the Padé approximated function is smaller than that {
a
n
} of the numerator.(a) The coefficient{
c
n
}
=
{
ϵ
r
n
} of a truncated random power series f
noise
1
[
100
]
(
z
) with ϵ
=
0.1. (b) The coefficients {
a
n
} and {
b
n
} of the Padé approximated function f
noise
1
[
50
∣
50
]
(
z
) for f
noise
1
[
100
]
(
z
).
(a)
(b)Note that the truncated random series is a random polynomial. As for the random polynomial, it is well known that the distribution of the zeros converges on the unit circle when the order of the random polynomial increases (Erdos-Turan-type theorem) [34, 35, 40]. Accordingly, we can generally interpret that in the Padé approximated function to the random power series the distribution of poles and zeros also accumulates around the unit circle when the order of the Padé approximation increases. The dependence of the zeros of the random polynomial and the zeros and poles of the Padé approximation has been studied by Gilewicz and Kryakin [42] and Ding and Xiao [43].


### 5.2. Effect of Noise on a Function with a Simple Pole

In the following subsections, we investigate influences of noise on the Padé approximation for some constructed noisy test functions as follows;(22)
f
test
+
noise
2
(
z
)
=
f
test
(
z
)
+
f
noise
2
(
z
)
,
where f
test
(
z
)
=
∑
n
=
0
∞
‍
a
n
z
n and f
noise
2
(
z
)
=
∑
n
=
0
∞
‍
ϵ
n
z
n. The {
ϵ
n
} is i.i.d. random variables within [
-
ϵ
,
ϵ
], where ϵ is the noise strength. Essentially, f
noise
2
(
z
) is the same as the random power series f
noise
1
(
z
). First of all, in this subsection, we consider a truncated function with a simple pole. Note that if a
n
=
C (constant) and ϵ
=
0, that is, in noise-free case, f
pole
+
noise
(
z
)
=
C
∑
n
=
0
∞
‍
z
n
=
C
/
(
1
-
z
) with a simple pole at z
=
1. In [2] by Baker Jr., the noise effect is summarized as follows: the [
M
∣
M
] Padé approximation has an unstable zero at the distance of order ϵ
-
1 from the origin, and the other zeros make (
M
-
1
) Froissart doublets (zero-pole pairs).Next, we consider a function(23)
f
pole
2
+
noise
2
z
=
f
pole
2
z
+
f
noise
2
z
,
=
∑
n
=
0
∞
‍
(
1
2
n
+
ϵ
n
)
z
n
with the noise strength ϵ
<
1. Note that
(24)
f
pole
2
(
z
)
=
2
(
2
-
z
)
with a simple pole at z
=
2 to clearly show the shift of the poles of the approximated function due to the noisy series.Figure8 shows distribution of the poles and zeros of the [
10
∣
10
] Padé approximated functions. It clearly shows the pole shift by the noise effect. In the noise-free case (ϵ
=
0), a pole of the Padé approximation appears at z
=
2 and the other poles are cancelled with zeros (zero-pole ghost pairs). In a case when the relatively small noise (ϵ
=
0.01) is added, the poles and zeros move toward |
z
|
=
1 with making Froissart doublets, although a pole at z
=
2 is quite stable. It becomes impossible to detect the true pole at z
=
2 when the noise strength is relatively large (ϵ
=
0.1), not shown in Figure 8.Figure 8
Distribution of poles (◯, △) and zeros (×, +) of the [
10
∣
10
] Padé approximated function f
pole
2
+
noise
2
[
10
∣
10
]
(
z
) with a stable pole at z
=
2 for noise strength ϵ
=
0, ϵ
=
0.01. The unit circle is drawn to guide the eye.As a result, it is found that the locations of the ghost pairs are unstable for noise, and the residues for the poles are much smaller than one corresponding to the true pole. We can guess that the proximity of the nonmodal poles and zeros of the Padé approximated function can be understood in a sense that the poles due to the noise need zeros to cancel with each other asϵ
→
0.


### 5.3. Effect of Noise on a Function with a Branch Cut

We investigate the effect of the noise on functions with a branch cut. First, let us consider a function(25)
f
branch
1
(
z
)
=
3
+
z
1
+
z
with an algebraic branch points at z
=
-
1 and z
=
-
3 and with the branch cut in [
-
3
,
-
1
]. Distribution of the poles and zeros of the Padé approximated function f
branch
1
+
noise
2
[
10
∣
10
]
(
z
) is shown in Figure 9. In a case with relatively small noise (ϵ
=
0.01), some poles make a line on the branch cut, and some poles and zeros move toward the unit circle |
z
|
=
1. It is impossible to detect the branch cut when the noise strength is relatively large (ϵ
=
0.1).Figure 9
Distribution of poles (◯, △) and zeros (×, +) of the [
10
∣
10
] Padé approximated function f
branch
1
+
noise
2
[
10
∣
10
]
(
z
) with a brunch cut from z
=
-
∞ to z
=
0 for the noise strength ϵ
=
0, ϵ
=
0.01. The unit circle is drawn to guide the eye.Next, let us consider a function(26)
f
branch
2
(
z
)
=
log
⁡
(
6
5
-
z
)
with a logarithmic branch point at z
=
6
/
5 and with a brunch cut from z
=
6
/
5 to z
=
∞. The distribution of the poles and zeros of the Padé approximated function f
branch
2
+
noise
2
[
10
∣
10
]
(
z
) for the f
branch
2
(
z
) with the noisy perturbation is shown in Figure 10. Some poles and zeros are making a line alternatively on the branch cut in the noise-free case (ϵ
=
0). It assembles around the unit circle |
z
|
=
1 with making Froissart doublets when the noise with strength ϵ
=
0.01 is added.Figure 10
Distribution of poles (◯, △) and zeros (×, +) of the [
10
∣
10
] Padé approximated function f
branch
2
+
noise
2
[
10
∣
10
]
(
z
) with a brunch cut from z
=
6
/
5 to z
=
∞ for the noise strength ϵ
=
0, ϵ
=
0.01. The unit circle is drawn to guide the eye.


### 5.4. Effect of Noise on a Function with a Natural Boundary

Figure11 shows distribution of the poles and zeros of the [
50
∣
50
] Padé approximated function for
(27)
f
Jac
+
noise
(
z
)
=
f
Jac
(
z
)
+
f
noise
2
(
z
)
,
which has a natural boundary on |
z
|
=
1.Figure 11
Distribution of poles (●, ◯, △) and zeros (×
+, *) of the Padé approximated function f
Jac
+
noise
[
50
∣
50
]
(
z
) for the lacunary series f
Jac
+
noise
[
100
]
(
z
) with a natural boundary on |
z
|
=
1. The noise strengths are ϵ
=
0, ϵ
=
0.001, and ϵ
=
0.01, respectively. The unit circle is drawn to guide the eye.In the noise-free case, the pairs of poles and zeros of the Padé approximated function are perfectly cancelled inside the unit circle|
z
|
=
1. The other poles and zeros of the Padé approximated function assemble around the circle |
z
|
=
1 without cancellation. In the relatively small noise case (ϵ
=
0.01), the location of the poles is not significantly changed compared with the zeros shifted outside the unit circle due to the noise effect. And, again, the poles and zeros move toward |
z
|
=
1 with making Froissart doublets when the noise strength is relatively large (ϵ
=
0.1). It is closely related to a fact that fluctuation of the coefficients of the numerator of the Padé approximated function is much larger than those in the denominator, as seen in Padé approximation to the random power series in Figure 7. As a result, the singularity of the Padé approximated function for the function with a natural boundary is more sensitive to the noisy perturbation than that in the functions with the other type singularity such as simple poles and branch points.It is very difficult to effectively distinguish whether the poles of the Padé approximation originated from the natural boundary on|
z
|
=
1 of the original function f
Jac
(
z
) or from the other natural boundary on |
z
|
=
1 generated by noisy series f
noise
2
(
z
) or numerical errors. Actually the round-off error affects the distribution of the poles and zeros of the Padé approximated function. Accordingly, to determine the expansion coefficients c
n with adequate accuracy becomes very important in the numerical calculation. This is a drawback of the Padé approximation when we use it for functions with unknown singularities.


### 5.5. Numerical Accuracy and Spurious Poles

As we observed in the last subsection, the effect of rounding error and accuracy limit of computers work in the numerical results of the Padé approximation. As the result of accumulation of the round-off error, the “spurious poles” appear around the unit circle|
z
|
=
1 as the pole-zero pairs when the order of Padé approximation increases (we used a term “Froissart doublets” for the poles-zero pairs generated by random power series, conveniently, although we cannot numerically distinguish it from the spurious poles due to the round-off errors; in the next section, we will discuss the Froissart doublets again).However, we can roughly distinguish between true poles and the spurious poles by “residue analysis” of the Padé approximated function because the spurious poles-zero pairs are unstable for the change of the order. In this subsection, we try to investigate the residues of the Padé approximation for some test functions. Up to now, the residue analysis has been mainly used for performance comparison between the different algorithms of the Padé approximation of the same order [30, 31]. On the other hand, it seems that the study by using the information of the residue analysis is still rare in the Padé approximation [10, 12].Generally, the rational polynomials of the diagonal Padé approximation can be uniquely identified by the poles{
z
k
} and the corresponding residues A
k as follows:
(28)
Q
M
(
z
)
P
M
(
z
)
=
∑
k
M
‍
A
k
z
-
z
k
,
where the residues are given by
(29)
A
k
=
Q
M
(
z
k
)
∏
j
(
≠
k
)
M
(
z
k
-
z
j
)
.Here, we investigate the convergence property of the magnitude of residues|
A
k
| arranged in descending order.Figure12 shows the absolute value of the residues |
A
k
| of some Padé approximated functions for the test function f
5
(
z
), which are arranged in descending order (note that they are noise-free cases.) The distribution of the poles and zeros of the Padé approximated functions is given in Figure 3. In a case of M
=
50, the magnitude of the all residues |
A
k
| is larger than O
(
1
0
-
3
), which correspond to the relevant poles arranged radially in eight directions from the true poles. On the other hand, in a case of M
=
75, the spurious poles appear and distribute around the unit circle |
z
|
=
1 (see Figure 3(b)). It is found that the absolute values of the residues corresponding the spurious poles are several order of magnitude smaller than the relevant poles.Figure 12
Absolute values of the residues|
A
k
| of the [
50
∣
50
] and [
75
∣
75
] Padé approximated functions for the test function f
5
(
z
) without noise. The |
A
k
| are arranged in descending order.Distribution of poles and zeros of the Padé approximated functionf
branch
2
[
20
∣
20
]
(
z
) for the test function f
branch
2
(
z
) is shown in Figure 13. The stable poles and zeros are lined on [
6
/
5
,
∞
], and the spurious poles appear around |
z
|
=
1. The magnitude of the residues of the spurious poles is also enormously small compared with that of the stable poles remaining with the increase of the order of the Padé approximation.(a) Distribution of poles (◯) and zeros (×) of the [
20
∣
20
] Padé approximated function f
branch
2
[
20
∣
20
]
(
z
). The unit circle is drawn to guide the eye. (b) Absolute values of the residues |
A
k
| of the [
10
∣
10
] and [
20
∣
20
] Padé approximated functions for the test function f
branch
2
(
z
) without noise. The residues are arranged in descending order.
(a)
(b)Figure14 is also the result of the residues analysis for the Padé approximated function for the test function f
Jac
(
z
) with a natural boundary on the unit circle |
z
|
=
1. In the [
50
∣
50
] Padé approximated function, the magnitude of the residues |
A
k
| is shown in changing the noise strengths ϵ
=
0,0.01,0.1 corresponding to poles-zeros distribution in Figure 11.(a) Absolute values of the residues|
A
k
| of the [
50
∣
50
] Padé approximated functions to the noise added test function f
Jac
+
noise
(
z
) with the strength ϵ
=
0,0.01,0.1. (b) Absolute values of the residues |
A
k
| of the Padé approximated functions f
Jac
[
50
∣
50
]
(
z
), f
Jac
[
32
∣
32
]
(
z
), and f
Jac
[
64
∣
64
]
(
z
) for the test function f
Jac
(
z
) without the noise. The residues are arranged in descending order.
(a)
(b)In the small noise case (ϵ
=
0.01), the results of the residue analysis for f
Jac
+
noise
2
(
z
) is almost the same as the noise-free case (ϵ
=
0), and in the case with relatively strong noise (ϵ
=
0.1), the noise shifts the magnitude of the residues with larger value. In addition, the result of the residue analysis of the noise-free cases for some different orders of the Padé approximation is shown in Figure 14(b). We should have in mind that the order is important when we apply Padé approximation to the lacunary power series because we should not take the order of the approximation in the gap of the series.


## 5.1. Random Power Series and Natural Boundary

Let us consider a random power series:(21)
f
noise
1
(
z
)
=
∑
n
=
0
∞
‍
ϵ
r
n
z
n
.
Here the coefficients, r
0
,
r
1
,
r
2
,
…, are i.i.d. random variables which take a value within r
n
∈
[
0,1
], and ϵ is the strength of the randomness. It is shown that, in general, the random power series has a natural boundary on the unit circle |
z
|
=
1 with probability one. Figure 6(a) shows distribution of poles and zeros of the [
50
∣
50
] Padé approximated function for f
noise
1
(
z
). Some pairs of poles and zeros are perfectly cancelled inside the circle |
z
|
=
1. On the other hand, almost all the poles and zeros of the Padé approximated function assemble around the circle |
z
|
=
1 and not cancelled. The pair of poles and zeros around the circle |
z
|
=
1 is called “Froissart doublets,” and it well corresponds to the natural boundary of f
noise
1
(
z
). The original function is also well approximated by the [
50
∣
50
] Padé approximation (see Figure 6(b)).(a) Distribution of poles (◯) and zeros (×) of the Padé approximated function f
noise
1
[
50
∣
50
]
(
z
) for a random power series f
noise
1
(
z
) with ϵ
=
1. The unit circle is drawn to guide the eye. (b) The Padé approximated function f
noise
1
[
50
∣
50
]
(
θ
) and the exact function f
Fib
(
θ
) in the polar form with r
=
1.0.
(a)
(b)Figure7 shows an example of the coefficients {
c
n
}
=
{
ϵ
r
n
} of the random power series and the coefficients {
a
n
} and {
b
n
} of the [
50
∣
50
] Padé approximated function. The fluctuation of the coefficient {
b
n
} that determines the poles of the Padé approximated function is smaller than that {
a
n
} of the numerator.(a) The coefficient{
c
n
}
=
{
ϵ
r
n
} of a truncated random power series f
noise
1
[
100
]
(
z
) with ϵ
=
0.1. (b) The coefficients {
a
n
} and {
b
n
} of the Padé approximated function f
noise
1
[
50
∣
50
]
(
z
) for f
noise
1
[
100
]
(
z
).
(a)
(b)Note that the truncated random series is a random polynomial. As for the random polynomial, it is well known that the distribution of the zeros converges on the unit circle when the order of the random polynomial increases (Erdos-Turan-type theorem) [34, 35, 40]. Accordingly, we can generally interpret that in the Padé approximated function to the random power series the distribution of poles and zeros also accumulates around the unit circle when the order of the Padé approximation increases. The dependence of the zeros of the random polynomial and the zeros and poles of the Padé approximation has been studied by Gilewicz and Kryakin [42] and Ding and Xiao [43].


## 5.2. Effect of Noise on a Function with a Simple Pole

In the following subsections, we investigate influences of noise on the Padé approximation for some constructed noisy test functions as follows;(22)
f
test
+
noise
2
(
z
)
=
f
test
(
z
)
+
f
noise
2
(
z
)
,
where f
test
(
z
)
=
∑
n
=
0
∞
‍
a
n
z
n and f
noise
2
(
z
)
=
∑
n
=
0
∞
‍
ϵ
n
z
n. The {
ϵ
n
} is i.i.d. random variables within [
-
ϵ
,
ϵ
], where ϵ is the noise strength. Essentially, f
noise
2
(
z
) is the same as the random power series f
noise
1
(
z
). First of all, in this subsection, we consider a truncated function with a simple pole. Note that if a
n
=
C (constant) and ϵ
=
0, that is, in noise-free case, f
pole
+
noise
(
z
)
=
C
∑
n
=
0
∞
‍
z
n
=
C
/
(
1
-
z
) with a simple pole at z
=
1. In [2] by Baker Jr., the noise effect is summarized as follows: the [
M
∣
M
] Padé approximation has an unstable zero at the distance of order ϵ
-
1 from the origin, and the other zeros make (
M
-
1
) Froissart doublets (zero-pole pairs).Next, we consider a function(23)
f
pole
2
+
noise
2
z
=
f
pole
2
z
+
f
noise
2
z
,
=
∑
n
=
0
∞
‍
(
1
2
n
+
ϵ
n
)
z
n
with the noise strength ϵ
<
1. Note that
(24)
f
pole
2
(
z
)
=
2
(
2
-
z
)
with a simple pole at z
=
2 to clearly show the shift of the poles of the approximated function due to the noisy series.Figure8 shows distribution of the poles and zeros of the [
10
∣
10
] Padé approximated functions. It clearly shows the pole shift by the noise effect. In the noise-free case (ϵ
=
0), a pole of the Padé approximation appears at z
=
2 and the other poles are cancelled with zeros (zero-pole ghost pairs). In a case when the relatively small noise (ϵ
=
0.01) is added, the poles and zeros move toward |
z
|
=
1 with making Froissart doublets, although a pole at z
=
2 is quite stable. It becomes impossible to detect the true pole at z
=
2 when the noise strength is relatively large (ϵ
=
0.1), not shown in Figure 8.Figure 8
Distribution of poles (◯, △) and zeros (×, +) of the [
10
∣
10
] Padé approximated function f
pole
2
+
noise
2
[
10
∣
10
]
(
z
) with a stable pole at z
=
2 for noise strength ϵ
=
0, ϵ
=
0.01. The unit circle is drawn to guide the eye.As a result, it is found that the locations of the ghost pairs are unstable for noise, and the residues for the poles are much smaller than one corresponding to the true pole. We can guess that the proximity of the nonmodal poles and zeros of the Padé approximated function can be understood in a sense that the poles due to the noise need zeros to cancel with each other asϵ
→
0.


## 5.3. Effect of Noise on a Function with a Branch Cut

We investigate the effect of the noise on functions with a branch cut. First, let us consider a function(25)
f
branch
1
(
z
)
=
3
+
z
1
+
z
with an algebraic branch points at z
=
-
1 and z
=
-
3 and with the branch cut in [
-
3
,
-
1
]. Distribution of the poles and zeros of the Padé approximated function f
branch
1
+
noise
2
[
10
∣
10
]
(
z
) is shown in Figure 9. In a case with relatively small noise (ϵ
=
0.01), some poles make a line on the branch cut, and some poles and zeros move toward the unit circle |
z
|
=
1. It is impossible to detect the branch cut when the noise strength is relatively large (ϵ
=
0.1).Figure 9
Distribution of poles (◯, △) and zeros (×, +) of the [
10
∣
10
] Padé approximated function f
branch
1
+
noise
2
[
10
∣
10
]
(
z
) with a brunch cut from z
=
-
∞ to z
=
0 for the noise strength ϵ
=
0, ϵ
=
0.01. The unit circle is drawn to guide the eye.Next, let us consider a function(26)
f
branch
2
(
z
)
=
log
⁡
(
6
5
-
z
)
with a logarithmic branch point at z
=
6
/
5 and with a brunch cut from z
=
6
/
5 to z
=
∞. The distribution of the poles and zeros of the Padé approximated function f
branch
2
+
noise
2
[
10
∣
10
]
(
z
) for the f
branch
2
(
z
) with the noisy perturbation is shown in Figure 10. Some poles and zeros are making a line alternatively on the branch cut in the noise-free case (ϵ
=
0). It assembles around the unit circle |
z
|
=
1 with making Froissart doublets when the noise with strength ϵ
=
0.01 is added.Figure 10
Distribution of poles (◯, △) and zeros (×, +) of the [
10
∣
10
] Padé approximated function f
branch
2
+
noise
2
[
10
∣
10
]
(
z
) with a brunch cut from z
=
6
/
5 to z
=
∞ for the noise strength ϵ
=
0, ϵ
=
0.01. The unit circle is drawn to guide the eye.


## 5.4. Effect of Noise on a Function with a Natural Boundary

Figure11 shows distribution of the poles and zeros of the [
50
∣
50
] Padé approximated function for
(27)
f
Jac
+
noise
(
z
)
=
f
Jac
(
z
)
+
f
noise
2
(
z
)
,
which has a natural boundary on |
z
|
=
1.Figure 11
Distribution of poles (●, ◯, △) and zeros (×
+, *) of the Padé approximated function f
Jac
+
noise
[
50
∣
50
]
(
z
) for the lacunary series f
Jac
+
noise
[
100
]
(
z
) with a natural boundary on |
z
|
=
1. The noise strengths are ϵ
=
0, ϵ
=
0.001, and ϵ
=
0.01, respectively. The unit circle is drawn to guide the eye.In the noise-free case, the pairs of poles and zeros of the Padé approximated function are perfectly cancelled inside the unit circle|
z
|
=
1. The other poles and zeros of the Padé approximated function assemble around the circle |
z
|
=
1 without cancellation. In the relatively small noise case (ϵ
=
0.01), the location of the poles is not significantly changed compared with the zeros shifted outside the unit circle due to the noise effect. And, again, the poles and zeros move toward |
z
|
=
1 with making Froissart doublets when the noise strength is relatively large (ϵ
=
0.1). It is closely related to a fact that fluctuation of the coefficients of the numerator of the Padé approximated function is much larger than those in the denominator, as seen in Padé approximation to the random power series in Figure 7. As a result, the singularity of the Padé approximated function for the function with a natural boundary is more sensitive to the noisy perturbation than that in the functions with the other type singularity such as simple poles and branch points.It is very difficult to effectively distinguish whether the poles of the Padé approximation originated from the natural boundary on|
z
|
=
1 of the original function f
Jac
(
z
) or from the other natural boundary on |
z
|
=
1 generated by noisy series f
noise
2
(
z
) or numerical errors. Actually the round-off error affects the distribution of the poles and zeros of the Padé approximated function. Accordingly, to determine the expansion coefficients c
n with adequate accuracy becomes very important in the numerical calculation. This is a drawback of the Padé approximation when we use it for functions with unknown singularities.


## 5.5. Numerical Accuracy and Spurious Poles

As we observed in the last subsection, the effect of rounding error and accuracy limit of computers work in the numerical results of the Padé approximation. As the result of accumulation of the round-off error, the “spurious poles” appear around the unit circle|
z
|
=
1 as the pole-zero pairs when the order of Padé approximation increases (we used a term “Froissart doublets” for the poles-zero pairs generated by random power series, conveniently, although we cannot numerically distinguish it from the spurious poles due to the round-off errors; in the next section, we will discuss the Froissart doublets again).However, we can roughly distinguish between true poles and the spurious poles by “residue analysis” of the Padé approximated function because the spurious poles-zero pairs are unstable for the change of the order. In this subsection, we try to investigate the residues of the Padé approximation for some test functions. Up to now, the residue analysis has been mainly used for performance comparison between the different algorithms of the Padé approximation of the same order [30, 31]. On the other hand, it seems that the study by using the information of the residue analysis is still rare in the Padé approximation [10, 12].Generally, the rational polynomials of the diagonal Padé approximation can be uniquely identified by the poles{
z
k
} and the corresponding residues A
k as follows:
(28)
Q
M
(
z
)
P
M
(
z
)
=
∑
k
M
‍
A
k
z
-
z
k
,
where the residues are given by
(29)
A
k
=
Q
M
(
z
k
)
∏
j
(
≠
k
)
M
(
z
k
-
z
j
)
.Here, we investigate the convergence property of the magnitude of residues|
A
k
| arranged in descending order.Figure12 shows the absolute value of the residues |
A
k
| of some Padé approximated functions for the test function f
5
(
z
), which are arranged in descending order (note that they are noise-free cases.) The distribution of the poles and zeros of the Padé approximated functions is given in Figure 3. In a case of M
=
50, the magnitude of the all residues |
A
k
| is larger than O
(
1
0
-
3
), which correspond to the relevant poles arranged radially in eight directions from the true poles. On the other hand, in a case of M
=
75, the spurious poles appear and distribute around the unit circle |
z
|
=
1 (see Figure 3(b)). It is found that the absolute values of the residues corresponding the spurious poles are several order of magnitude smaller than the relevant poles.Figure 12
Absolute values of the residues|
A
k
| of the [
50
∣
50
] and [
75
∣
75
] Padé approximated functions for the test function f
5
(
z
) without noise. The |
A
k
| are arranged in descending order.Distribution of poles and zeros of the Padé approximated functionf
branch
2
[
20
∣
20
]
(
z
) for the test function f
branch
2
(
z
) is shown in Figure 13. The stable poles and zeros are lined on [
6
/
5
,
∞
], and the spurious poles appear around |
z
|
=
1. The magnitude of the residues of the spurious poles is also enormously small compared with that of the stable poles remaining with the increase of the order of the Padé approximation.(a) Distribution of poles (◯) and zeros (×) of the [
20
∣
20
] Padé approximated function f
branch
2
[
20
∣
20
]
(
z
). The unit circle is drawn to guide the eye. (b) Absolute values of the residues |
A
k
| of the [
10
∣
10
] and [
20
∣
20
] Padé approximated functions for the test function f
branch
2
(
z
) without noise. The residues are arranged in descending order.
(a)
(b)Figure14 is also the result of the residues analysis for the Padé approximated function for the test function f
Jac
(
z
) with a natural boundary on the unit circle |
z
|
=
1. In the [
50
∣
50
] Padé approximated function, the magnitude of the residues |
A
k
| is shown in changing the noise strengths ϵ
=
0,0.01,0.1 corresponding to poles-zeros distribution in Figure 11.(a) Absolute values of the residues|
A
k
| of the [
50
∣
50
] Padé approximated functions to the noise added test function f
Jac
+
noise
(
z
) with the strength ϵ
=
0,0.01,0.1. (b) Absolute values of the residues |
A
k
| of the Padé approximated functions f
Jac
[
50
∣
50
]
(
z
), f
Jac
[
32
∣
32
]
(
z
), and f
Jac
[
64
∣
64
]
(
z
) for the test function f
Jac
(
z
) without the noise. The residues are arranged in descending order.
(a)
(b)In the small noise case (ϵ
=
0.01), the results of the residue analysis for f
Jac
+
noise
2
(
z
) is almost the same as the noise-free case (ϵ
=
0), and in the case with relatively strong noise (ϵ
=
0.1), the noise shifts the magnitude of the residues with larger value. In addition, the result of the residue analysis of the noise-free cases for some different orders of the Padé approximation is shown in Figure 14(b). We should have in mind that the order is important when we apply Padé approximation to the lacunary power series because we should not take the order of the approximation in the gap of the series.


## 6. Froissart Doublets

The problem of constructing theZ-transform Z
(
z
) of a finite time-series is a standard problem in mathematics [10–14]. For example, it is shown that for a sum of oscillating damped signals, the Z-transform associated with the time-series can be characterized by a sum of the poles of the Padé approximated function. The position of each pole is simply linked to the damping factor and the frequency of each of the oscillators. Also, it is important to note that all these poles lie strictly outside the unit circle because it corresponds to the damping [10–13]. In addition, we will consider quasianalyticity property of the random power series by the residue analysis of the Padé approximation.


### 6.1. Noise Attractor

In signal processing, we can use the fact that the poles and zeros of the Padé approximated function to the noisy series distribute around the unit circle|
z
|
=
1 when we remove the noise from the observed data through Z-transform and/or Fourier transform of the data. Let a sequence {
s
0
,
s
1
,
…
,
s
n
,
…
} be a sample signal without noise. Then we define the Z-transform of the sequence as
(30)
Z
(
z
)
=
∑
n
=
0
N
‍
s
n
z
n
.
The function Z
(
z
) is analytic interior of |
z
|
<
1 if the number of signals N is finite [44]. Note that discrete Fourier transform is a special case of the Z-transform.Next, let us consider a signal sequence int
∈
[
0
,
T
] consisting of the superimposed damping oscillators as
(31)
s
k
=
∑
l
‍
A
l
e
i
ω
l
(
k
/
N
)
T
,
k
=
0,1
,
…
,
N
-
1
,
where A
l is the amplitude of the lth oscillator and ω
l
=
2
π
f
l
+
i
α
l. Here, f
l and α
l are the frequency and the damping factor of the lth oscillator. Then, the Z-transform is
(32)
Z
z
=
∑
n
=
0
∞
‍
s
n
z
n
=
∑
n
=
0
∞
‍
∑
l
‍
A
l
e
i
ω
l
(
n
/
N
)
T
z
n
=
∑
l
‍
A
l
1
-
z
z
l
,
where we take a limit n
→
∞ keeping T
/
N and z
l
≡
e
i
ω
l
(
T
/
N
). Accordingly, the singularity of Z
(
z
) appears as the poles at z
=
z
l
-
1
≡
e
-
i
ω
l
(
T
/
N
) outside the unit circle |
z
|
>
1, and the residue is Re
s
(
z
l
-
1
)
=
z
l
-
1
A
l.On the other hand, let us consider a noise-added sequence{
S
0
,
S
1
,
…
,
S
n
,
…
}. Then, the Froissart pointed out that there are two types of the poles: stable poles and unstable poles when we apply the diagonal Padé approximation to the unknown data set. In general, the Z-transform Z
(
z
)
=
∑
n
=
0
N
‍
S
n
z
n of the noisy sequence has a natural boundary on the unit circle |
z
|
=
1 with probability 1. In fact, the poles and zeros (Froissart doublets) of the Padé approximated function often distribute around the unit circle when the numerical error and/or noise are mixed into the Taylor series of the analytic functions, as seen in the last section. That is to say, we sometimes call the unit circle |
z
|
=
1
noise attractor in a sense that the poles and zeros are attracted to the circle as the Froissart doublets [45]. Accordingly, it is found that Padé approximated function for the function Z
(
z
) has stable poles associated with the damping modes and unstable spurious poles associated with the noisy fluctuation. After elimination of the spurious poles around the noise attractor from the noisy sequence we can reconstruct the noise-free sequence consisting of the stable poles located in the domain |
z
|
>
1. Another remarkable feature of the nonmodal poles is that the absolute values of the Cauchy residues associated with them are usually much smaller than those associated with true poles.


### 6.2. Random Power Series and Quasianalytic Function

Weierstrass defined the analytic function by direct analytic continuation of function. Then, apparently the analytic continuation is impossible beyond the natural boundary even if we can uniquely define the function and it is analytic outside the analytic domain. Borel and Gammel extended the narrow condition for the analyticity and gave a definition ofquasianalytic functions [46, 47]. Gammel conjectured the following for the random power series [10, 45].Gammel Conjecture (1973). The random power series belongs to the Borel class of quasianalytic functions as the following form:
(33)
f
Gammel
(
z
)
=
∑
k
=
0
∞
‍
B
k
1
-
w
k
z
,
where w
k
=
e
i
2
π
X
k and {
X
k
} are real numbers in the interval in X
k
∈
[
0,1
] and B
k decreases rapidly with k. Then, the natural boundary in the Weierstrass sense can be crossed.The function (33) is a simple example that poles are densely distributed on the unit circle. Then the convergence property of the sequence |
B
k
| is important for the analyticity of the function. Carleman proved that f
Gammel
(
z
) is quasianalytic if B
k satisfies the following condition:
(34)
|
B
k
|
<
C
e
-
k
1
+
c
,
c
>
0
.
This is Carleman class of quasianalytic functions. See Gammel’s paper [45] for the details. Moreover, Gammel and Nuttall proved that the quasianalytic functions can be exactly approximated by the Padé approximation [45].Gammel-Nuttall Theorem (1973). If B
k in (33) satisfies the condition (34) and |
ω
k
|
=
1, then the sequence of [
N
+
J
∣
N
] Padé approximation to the f
Gammel
(
z
) converges in measure to the function f
Gammel
(
z
) as N
→
∞ in any closed, bounded region of the complex plane, where J is a natural number that equals N or less.Is the Gammel Conjecture True? We try to examine the validity of the Gammel conjecture by applying residue analysis of the Padé approximated function to the random power series f
noise
2
(
z
). Figure 15 shows the absolute values of the residues |
A
k
| of the Padé approximated functions f
noise
[
45
∣
45
]
(
z
) for three different samples in descending order. |
A
k
| roughly exponentially decreases with respect to k as
(35)
A
k
~
exp
⁡
(
-
β
k
)
,
where β is the decay exponent. It shows exponential decay (or faster), and, on the surface, supports the Gammal conjecture.Figure 15
Absolute values of the residues|
A
k
| of the Padé approximated functions f
noise
2
[
45
∣
45
]
(
z
) for three samples of the truncated random power series of order N
=
90, and ϵ
=
1. The |
A
k
| are arranged in descending order.However, it is not nearly so simple. We should check the stability of the exponential-like decay of the magnitude of the residues by changing the order of the Padé approximation. Figure16 shows the result for the three different orders: M
=
15, M
=
45, and M
=
55. It expresses an indication that the decay exponent β does not converge to a positive certain value. It seems that the exponent behaves β
→
0 as a limit M
→
∞. On the other hand, if we directly apply the Padé approximation to the quasianalytic function f
Gammel
(
z
) with B
k
=
e
-
k, the exponent β is stable for changing the order of the Padé approximation (see Appendix E). These facts suggest that the random power series does not belong to Carleman class of quasianalytic functions although it has a natural boundary on the unit circle, and it has the form (33). As a result, we can say that no optimism is warranted on the Gammel conjecture.Figure 16
Absolute values of the residues|
A
k
| of the Padé approximated functions f
noise
2
[
45
∣
45
]
(
z
), f
noise
2
[
25
∣
25
]
(
z
), and f
noise
2
[
15
∣
15
]
(
z
) for a sample truncated random series, and ϵ
=
1. The |
A
k
| are arranged in descending order.How does the residue analyses of the Padé approximation for the analyticity and/or quasianalyticity of unknown function work? It is an interesting and future problem.


## 6.1. Noise Attractor

In signal processing, we can use the fact that the poles and zeros of the Padé approximated function to the noisy series distribute around the unit circle|
z
|
=
1 when we remove the noise from the observed data through Z-transform and/or Fourier transform of the data. Let a sequence {
s
0
,
s
1
,
…
,
s
n
,
…
} be a sample signal without noise. Then we define the Z-transform of the sequence as
(30)
Z
(
z
)
=
∑
n
=
0
N
‍
s
n
z
n
.
The function Z
(
z
) is analytic interior of |
z
|
<
1 if the number of signals N is finite [44]. Note that discrete Fourier transform is a special case of the Z-transform.Next, let us consider a signal sequence int
∈
[
0
,
T
] consisting of the superimposed damping oscillators as
(31)
s
k
=
∑
l
‍
A
l
e
i
ω
l
(
k
/
N
)
T
,
k
=
0,1
,
…
,
N
-
1
,
where A
l is the amplitude of the lth oscillator and ω
l
=
2
π
f
l
+
i
α
l. Here, f
l and α
l are the frequency and the damping factor of the lth oscillator. Then, the Z-transform is
(32)
Z
z
=
∑
n
=
0
∞
‍
s
n
z
n
=
∑
n
=
0
∞
‍
∑
l
‍
A
l
e
i
ω
l
(
n
/
N
)
T
z
n
=
∑
l
‍
A
l
1
-
z
z
l
,
where we take a limit n
→
∞ keeping T
/
N and z
l
≡
e
i
ω
l
(
T
/
N
). Accordingly, the singularity of Z
(
z
) appears as the poles at z
=
z
l
-
1
≡
e
-
i
ω
l
(
T
/
N
) outside the unit circle |
z
|
>
1, and the residue is Re
s
(
z
l
-
1
)
=
z
l
-
1
A
l.On the other hand, let us consider a noise-added sequence{
S
0
,
S
1
,
…
,
S
n
,
…
}. Then, the Froissart pointed out that there are two types of the poles: stable poles and unstable poles when we apply the diagonal Padé approximation to the unknown data set. In general, the Z-transform Z
(
z
)
=
∑
n
=
0
N
‍
S
n
z
n of the noisy sequence has a natural boundary on the unit circle |
z
|
=
1 with probability 1. In fact, the poles and zeros (Froissart doublets) of the Padé approximated function often distribute around the unit circle when the numerical error and/or noise are mixed into the Taylor series of the analytic functions, as seen in the last section. That is to say, we sometimes call the unit circle |
z
|
=
1
noise attractor in a sense that the poles and zeros are attracted to the circle as the Froissart doublets [45]. Accordingly, it is found that Padé approximated function for the function Z
(
z
) has stable poles associated with the damping modes and unstable spurious poles associated with the noisy fluctuation. After elimination of the spurious poles around the noise attractor from the noisy sequence we can reconstruct the noise-free sequence consisting of the stable poles located in the domain |
z
|
>
1. Another remarkable feature of the nonmodal poles is that the absolute values of the Cauchy residues associated with them are usually much smaller than those associated with true poles.


## 6.2. Random Power Series and Quasianalytic Function

Weierstrass defined the analytic function by direct analytic continuation of function. Then, apparently the analytic continuation is impossible beyond the natural boundary even if we can uniquely define the function and it is analytic outside the analytic domain. Borel and Gammel extended the narrow condition for the analyticity and gave a definition ofquasianalytic functions [46, 47]. Gammel conjectured the following for the random power series [10, 45].Gammel Conjecture (1973). The random power series belongs to the Borel class of quasianalytic functions as the following form:
(33)
f
Gammel
(
z
)
=
∑
k
=
0
∞
‍
B
k
1
-
w
k
z
,
where w
k
=
e
i
2
π
X
k and {
X
k
} are real numbers in the interval in X
k
∈
[
0,1
] and B
k decreases rapidly with k. Then, the natural boundary in the Weierstrass sense can be crossed.The function (33) is a simple example that poles are densely distributed on the unit circle. Then the convergence property of the sequence |
B
k
| is important for the analyticity of the function. Carleman proved that f
Gammel
(
z
) is quasianalytic if B
k satisfies the following condition:
(34)
|
B
k
|
<
C
e
-
k
1
+
c
,
c
>
0
.
This is Carleman class of quasianalytic functions. See Gammel’s paper [45] for the details. Moreover, Gammel and Nuttall proved that the quasianalytic functions can be exactly approximated by the Padé approximation [45].Gammel-Nuttall Theorem (1973). If B
k in (33) satisfies the condition (34) and |
ω
k
|
=
1, then the sequence of [
N
+
J
∣
N
] Padé approximation to the f
Gammel
(
z
) converges in measure to the function f
Gammel
(
z
) as N
→
∞ in any closed, bounded region of the complex plane, where J is a natural number that equals N or less.Is the Gammel Conjecture True? We try to examine the validity of the Gammel conjecture by applying residue analysis of the Padé approximated function to the random power series f
noise
2
(
z
). Figure 15 shows the absolute values of the residues |
A
k
| of the Padé approximated functions f
noise
[
45
∣
45
]
(
z
) for three different samples in descending order. |
A
k
| roughly exponentially decreases with respect to k as
(35)
A
k
~
exp
⁡
(
-
β
k
)
,
where β is the decay exponent. It shows exponential decay (or faster), and, on the surface, supports the Gammal conjecture.Figure 15
Absolute values of the residues|
A
k
| of the Padé approximated functions f
noise
2
[
45
∣
45
]
(
z
) for three samples of the truncated random power series of order N
=
90, and ϵ
=
1. The |
A
k
| are arranged in descending order.However, it is not nearly so simple. We should check the stability of the exponential-like decay of the magnitude of the residues by changing the order of the Padé approximation. Figure16 shows the result for the three different orders: M
=
15, M
=
45, and M
=
55. It expresses an indication that the decay exponent β does not converge to a positive certain value. It seems that the exponent behaves β
→
0 as a limit M
→
∞. On the other hand, if we directly apply the Padé approximation to the quasianalytic function f
Gammel
(
z
) with B
k
=
e
-
k, the exponent β is stable for changing the order of the Padé approximation (see Appendix E). These facts suggest that the random power series does not belong to Carleman class of quasianalytic functions although it has a natural boundary on the unit circle, and it has the form (33). As a result, we can say that no optimism is warranted on the Gammel conjecture.Figure 16
Absolute values of the residues|
A
k
| of the Padé approximated functions f
noise
2
[
45
∣
45
]
(
z
), f
noise
2
[
25
∣
25
]
(
z
), and f
noise
2
[
15
∣
15
]
(
z
) for a sample truncated random series, and ϵ
=
1. The |
A
k
| are arranged in descending order.How does the residue analyses of the Padé approximation for the analyticity and/or quasianalyticity of unknown function work? It is an interesting and future problem.


## 7. Summary and Discussion

In the present paper we numerically examined the effectiveness of the Padé approximation for some test functions with branch point, essential singularity, and natural boundary by watching the singularities of the Padé approximated functions. For the functions with a branch cut, the poles and zeros of the Padé approximated function are lined along the true branch cut. The poles and zeros are distributed around the true natural boundary if the original test function has a natural boundary. In addition, we gave the explicit Padé approximated functions for some lacunary power series which are useful to check the numerical result. It was shown that, in particular, the distribution of poles and zeros of the Padé approximated function for lacunary power series and the random power series accumulated around the unit circle when the order of the approximation increases.We often suffer from the difficulty to distinguish whether or not the poles of the Padé approximation are intrinsically originated from the natural boundary of the original power series because the numerical errors contained in the expansion coefficients also yield a false natural boundary. Therefore, the expansion coefficients with adequate numerical accuracy are necessary when we apply the Padé approximation to functions with unknown singularities.Furthermore, the residue calculus of the Padé approximated function is useful when we detect the singularity of the original power series from the asymptotic behavior of the truncated series. It is useful also for estimating the accuracy of the approximation. As a result, the residue calculus suggested that the random power series does not obey Gammel’s conjecture; that is, it does not belong to Borel class of the quasianalytic functions.We finally remark that the most serious problem to be improved is the numerical accuracy due to the limitation of the order in the Padé approximation when we use it for detecting unknown singularities of wave functions in quantum physics [32].



---

*Source: 587430-2014-11-20.xml*