# An analog of the Sinc approximation for periodic functions

In this paper, we propose an interpolation formula for periodic functions. This formula can be regarded as an analog of the Sinc approximation, which is an interpolation formula for functions defined on the entire infinite interval. Theoretical error analysis and numerical examples show that the proposed formula converges exponentially for analytic periodic functions.

04/06/2020

### On an optimal interpolation formula in K_2(P_2) space

The paper is devoted to the construction of an optimal interpolation for...
04/17/2018

### On indefinite sums weighted by periodic sequences

For any integer q≥ 2 we provide a formula to express indefinite sums of ...
05/10/2020

### New conformal map for the trapezoidal formula for infinite integrals of unilateral rapidly decreasing functions

While the trapezoidal formula can attain exponential convergence when ap...
08/28/2019

### A Method for Dimensionally Adaptive Sparse Trigonometric Interpolation of Periodic Functions

We present a method for dimensionally adaptive sparse trigonometric inte...
12/30/2018

### New conformal map for the Sinc approximation for exponentially decaying functions over the semi-infinite interval

The Sinc approximation has shown its high efficiency for numerical metho...
04/19/2021

### A linear barycentric rational interpolant on starlike domains

When an approximant is accurate on the interval, it is only natural to t...
02/01/2016

### Holonomic Tools for Basic Hypergeometric Functions

With the exception of q-hypergeometric summation, the use of computer al...

## 1 Introduction

The trapezoidal formula with equal mesh is an elementary numerical integration formula, and it works very well for the integral of an analytic function over the entire infinite interval and the integral of an analytic periodic function over an interval of its period length, that is, the formula

 ∫∞−∞f(x)dx≃h∞∑k=−∞f(kh) (1)

is efficient if is an analytic function on the entire infinite interval , and

 ∫a0f(x)dx≃h{12f(0)+N−1∑k=1f(kh)+12f(a)}(h=aN) (2)

is efficient if is an analytic periodic function of period . Especially, the formula (1) over the entire infinite interval is widely used in science and engineering combined with the technique of a variable transform [9].

By the way, there exists an interpolation formula with equal mesh over the entire infinite interval. It is the Sinc approximation [7]

 f(x)≃Lh[f](x)=∞∑k=−∞f(kh)sin[(π/h)(x−kh)](π/h)(x−kh). (3)

It is an interpolation formula with sampling points , that is,

 Lh[f](kh)=f(kh)(k∈Z).

it works very well especially if is an analytic function over the entire infinite interval , and we obtain the trapezoidal formula over the entire infinite interval by integrating (3) formally term by term, that is,

 ∫∞−∞f(x)dx≃ ∫∞−∞Lh[f](x)dx = ∞∑k=−∞f(kh)∫∞−∞sin[(π/h)(x−kh)](π/h)(x−kh)dx = h∞∑k=−∞f(kh).

As the trapezoidal formula for numerical integration over the entire infinite interval , the Sinc approximation is applied to various numerical computations combined with the technique of variable transforms such as indefinite integrals [3]

, initial value problems of ordinary differential equations

[5], boundary value problems of ordinary differential equations [8] and integral equations [4, 6].

Taking into account that

1. the trapezoidal formula works very well for analytic functions over the entire infinite interval and analytic periodic functions, and

2. the Sinc approximation is an interpolation formula over the entire infinite interval which generates the trapezoidal formula over the entire infinite interval,

it is a natural question whether there exists an interpolation formula for periodic functions with equal mesh which generates the trapezoidal formula for periodic functions, that is,

 L[f](kh)=f(kh)(k=0,1,2,…,N−1;h=aN), (4) where a is the period of f(x), and ∫a0L[f](x)dx=h{12f(0)+N−1∑k=1f(kh)+12f(a)}. (5)

This paper gives the answer “yes” to this question, that is, there exists an interpolation formula for periodic functions which satisfies the properties (4) and (5).

The remainder of this paper are structured as follows. In Section 2, we propose an interpolation formula for periodic functions which we desire. In Section 3, we show a theorem on the exponential convergence of the proposed method for periodic analytic functions. In Section 4, we show the effectiveness of the proposed method by some numerical examples. In Section 5, we show concluding remarks.

## 2 Sinc approximation for periodic functions

Let be a periodic function of period . We propose the following interpolation formula for .

 f(x)≃ LN[f](x) (6) = ⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩12N0+12N0∑k=0f(ka2N0+1)sin[(2N0+1)πa(x−ka2N0+1)]sin[πa(x−ka2N0+1)]if N=2N0+1 (odd),12N02N0−1∑k=0f(ka2N0)cot[πa(x−ka2N0)]sin[2N0πa(x−ka2N0)]if N=2N0 (even).

We can easily check that

1. is an interpolation of with sampling points , that is,

 LN[f](kaN)=f(kaN)(k=0,1,2,…,N−1).
2. is a periodic function of period , that is,

 LN[f](x+a)=LN[f](x).
3. The trapezoidal formula for over with an equal mesh is obtained by integrating , that is,

 aN{12f(0)+N−1∑k=1f(kaN)+12f(a)}=∫a0LN[f](x)dx.
4. The limit as and with fixed gives the Sinc approximation

 LN[f](x)→Lh[f](x)=∞∑k=−∞f(kh)sin[(π/h)(x−kh)](π/h)(x−kh).

The first and second properties are satisfied also by the Sinc approximation. We call the approximation (6) the Sinc approximation for periodic functions.

### Derivation

We show here how to find the Sinc approximation for periodic functions (6).

We begin with the Fourier transform of the Sinc approximation over the entire infinite interval

 F[Lh[f]](ξ)=∫∞−∞Lh[f](x)e−2πiξxdx=h∞∑k=−∞f(kh)e−2πikhχ(−π/h,π/h)(ξ), (7) where χ(−π/h,π/h)(ξ)is the characteristic function of the interval (−π/h,π/h) χ(−π/h,π/h)(ξ)={1ξ∈(−π/h,π/h)0ξ∉[−π/h,π/h].

We compare (7) with the Fourier transform of

 F[f](ξ)=∫∞−∞f(x)e−2πiξxdx. (8)

We find that is non-zero only in the interval and, in (7), the integral on the right-hand side of (8) is replaced with its approximation by the trapezoidal formula of equal mesh .

Let be a periodic function of period and expand it to the Fourier series

 f(x)=∑n∈Zˆfnexp(i2nπax), (9) where ˆfn=1a∫∞−∞f(x)exp(−i2nπax)dx (10)

From the above observation of the Sinc approximation over the entire infinite interval, we expect that we can obtain the Sinc approximation of the periodic function by truncating the infinite sum on the right-hand side of (8) and replace the integral on the right-hand side of (10) by its approximation by the trapezoidal formula.

First, we derive the Sinc approximation of the case that the number of sampling points is odd. We truncate the right-hand side of (9) as

 f(x)≃N0∑n=−N0ˆfnexp(i2nπax). (11)

Then, we approximate the integral on the right-hand side of (10) by its approximation by the -point trapezoidal formula

 ˆfn≃ ˆf(2N0+1)n (12) = 1aa2N0+12N0∑k=0f(ka2N0+1)exp(−i2nπaka2N0+1) = 12N0+12N0∑k=0f(ka2N0+1)exp(−i2nkπ2N0+1).

Substituting (12) into (11), we have

 f(x)≃ N0∑n=−N0ˆf(2N0+1)nexp(i2nπax) = 12N0+12N0∑k=0f(ka2N0+1)N0∑n=−N0exp[i2nπa(x−ka2N0+1)].

Using the formula

 N0∑n=−N0eiκn=sin[(N0+1/2)κ]sin(κ/2),

we obtain the Sinc approximation (6) for .

Second, we derive the Sinc approximation for the case that the number of sampling points is even. We truncate the Fourier series (9) as

 f(x)≃ N0−1∑n=−N0−1ˆfnexp(i2nπax)+12ˆf2N0exp(i2(2N0)πax) (13) +12ˆf−2N0exp(i2⋅(−2N0)πax).

Then, we approximate the integral on the right-hand side of (10) by its approximation by the -point trapezoidal formula

 ˆfn≃ ˆf(2N0)n (14) = 1aa2N02N0−1∑k=0f(ka2N0)exp(−i2nπaka2N0) = 12N02N0−1∑k=0f(ka2N0)exp(−i2nkπ2N0).

Substituting (14) into (13), we have

 f(x)≃ N0−1∑n=−N0+1ˆf(2N0)nexp(i2nπax)+12ˆf(2N0)N0exp(i2N0πax) (15) +12ˆf(2N0)−N0exp(−i2N0πax) = 12N02N0−1∑k=0f(ka2N0)⎧⎨⎩N0−1∑n=−N0+1exp[i2nπa(x−ka2N0)] +12exp[i2N0πa(x−ka2N0)]+12exp[−i2N0πa(x−ka2N0)]}.

Using the formula

 N0−1∑n=−N0+1eiκn+12eiκN0+12e−iκN0=cot(κ2)sin(N0κ),

we obtain the Sinc approximation (6) for .

## 3 Convergence theorem

The following theorem shows that the Sinc approximation for periodic functions converges exponentially as the number of sampling points increases if the function is an analytic periodic function.

###### Theorem 1

We suppose that is a periodic function of period and analytic in a complex domain containing the strip . Then we have the inequality

 |f(x)−LN[f](x)|≦Cd,a∥f∥dexp(−πdaN)(∀x∈R), (16)

where

 ∥f∥d=max|Imz|=d|f(z)| (17)

and is a positive constant depending and only.

We compare this theorem with the following theorem for the Sinc approximation over the entire infinite interval (Theorem 3.1.3 in [7]).

###### Theorem 2

Let be the following complex domain defined for by

 Dd(ϵ)={z∈C||Rez|<1/ϵ,|Imz|

We suppose that a function is analytic in the strip domain

 Dd={z∈C||Imz|

and satisfies

 N(f,Dd)≡limϵ→0∫∂Dd(ϵ)|f(z)||dz|<∞.

Then, we have

 |f(x)−Lh[f](x)|≦1πdN(f,Dd)exp(−πd/h)1−exp(−2πd/h)(∀x∈R).

This theorem shows that the error of the Sinc approximation is of order as the mesh goes to zero. On the other hand, the error of the approximation (6) for periodic analytic functions is of order where is the mesh size . Therefore, the error of the Sinc approximation for periodic functions is of the same order as that of the Sinc approximation over the entire infinite interval.

### Proof of Theorem 1

We prove Theorem 1 only for . We can prove the theorem for similarly. Since

 f(x)−L2N0+1[f](x)=N0∑n=−N0(ˆfn−ˆf(2N0+1)n)exp(i2nπax)−∑|n|>N0ˆfnexp(i2nπax),

we have

 |f(x)−L2N0+1[f](x)|≦N0∑n=−N0|ˆfn−ˆf(2N0+1)n|+∑|n|>N0|ˆfn|. (18)

Each term of the first term on the right-hand side of (18) is the error of the approximation of the integral

by the -point trapezoidal formula of the equal mesh , and, from the theorem in §4.6.5 of [1], we have

 |ˆfn−ˆf(2N0+1)n|≦2∥f∥dexp(2πda|n|)exp(−(2πd/a)(2N0+1))1−exp(−2πd/a).

Then, we have

 |ˆfn−ˆf(2N0+1)n|≦ 2∥f∥dexp(−(2πd/a)(2N0+1))1−exp(−2πd/a)N0∑n=−N0exp(2πda|n|) ≦ 2∥f∥dexp(−(2πd/a)N0)[1−exp(−2πd/a)]2.

As for the second term on the right-hand side of (18), since

 |fn|≦∥f∥dexp(−2|n|πda)(n≠0),

we have

 ∑|n|>N0|ˆfn|≦ ∥f∥d∑|n|>N0exp(−2|n|πda) ≦ 2∥f∥dexp(−(2πd/a)(N0+1))1−exp(−2πd/a).

Therefore, we have (16) for .

## 4 Numerical examples

We evaluated the errors of the -point Sinc approximation and the -point Chebyshev interpolation

 ϵN=max0≦x≦a|f(x)−LN[f](x)|, (19)

where is the Sinc approximation or the Chebyshev interpolation, of the periodic functions with period

 (1)exp(sinx),(2)12+cosx,(3)15+cosx. (20)

The maximum on the right-hand side of (19) was actually evaluated on equidistant points in the interval . All the computations were performed using programs coded in C++ and multi-precision arithmetics with decimal digit precision by the software exflib [2]. Figure 1 shows the errors of the two methods applied to the functions (20). From these figures, the Sinc approximation is superior to the Chebyshev interpolation and converges exponentially as the number of sampling points increases.

We evaluated the decay rate of the error of the Sinc approximation using the fit command of the software gnuplot. Table 1

shows the results. We can also estimate the decay rate of the error

by using Theorem 1 for the functions (2) and (3) in (20). Since these functions have poles respectively at

 (2) z=−ilog(−2±√3)=∓1.3170…i+(2n+1)π(n∈Z), (3) z=−ilog(−5±√212)=∓1.5668…i+(2n+1)π(n∈Z),

the value in Theorem 1 is estimated as for the function (2) and for the function (3). Using them, we can estimate theoretically the decay rates of the error of the Sinc approximation which are shown in Table 1. From this table, the theoretical error estimate make a good agreement with the numerical result especially for the function (2).

## 5 Summary

In this paper, we proposed an approximation of periodic functions, which can be regarded as an analog of the Sinc approximation for functions over the entire infinite interval in the sense that it has equidistant sampling points and it gives the trapezoidal formula with equal mesh by integrating it. From theoretical error estimate and numerical examples, it converges exponentially as the number of sampling points increases if the approximated function is analytic.

As an application of the proposed approximation, it is expected that it can be applied to the boundary element method for boundary value problems of partial differential equations.

## References

• [1] P. J. Davis and P. Rabinowitz (1984) Methods of numerical integration, second ed.. Academic Press, San Diego. Cited by: §3.
• [2] H. Fujiwara Exflib information. Cited by: §4.
• [3] M. Muhammad and M. Mori (2003) Double exponential formula for numerical indefinite integration. J. Comput. Appl. Math. 161, pp. 431–448. Cited by: §1.
• [4] M. Muhammad, A. Nurmuhammad, M. Mori, and M. Sugihara (2005) Numerical solution of integral equations by means of the sinc collocation method based on the double exponential transformation. J. Comput. Appl. Math. 177, pp. 269–286. Cited by: §1.
• [5] A. Nurmuhammad, M. Muhammad, and M. Mori (2005) Numerical solution of initial value problems based on the double exponential transformation. Publ. RIMS, Kyoto Univ. 41, pp. 937–948. Cited by: §1.
• [6] T. Okayama, T. Matsuo, and M. Sugihara (2011) Improvement of a sinc-collocation method for fredholm integral equations of the second kind. BIT Numer. Math. 51, pp. 339–366. Cited by: §1.
• [7] F. Stenger (1993) Numerical methods based on sinc and analytic functions. Springer-Verlag, New York. Cited by: §1, §3.
• [8] M. Sugihara (2002) Double exponential transformation in the sinc-collocation method for two-point boundary value problems. J. Comput. Appl. Maht. 149, pp. 239–250. Cited by: §1.
• [9] H. Takahasi and M. Mori (1978) Double exponential formulas for numerical integration. Publ. Res. Inst. Math. Sci., Kyoto Univ. 339, pp. 721–741. Cited by: §1.