Advances in Numerical Analysis

Volume 2013 (2013), Article ID 263467, 8 pages

http://dx.doi.org/10.1155/2013/263467

Research Article

A New Extended Padé Approximation and Its Application

Z. Kalateh Bojdi,¹ S. Ahmadi-Asl,¹ and A. Aminataei²

¹Department of Mathematics, Birjand University, Birjand, Iran

²Faculty of Mathematics, K. N. Toosi University of Technology, P.O. Box 16315-1618, Tehran, Iran

Received 19 June 2013; Revised 16 September 2013; Accepted 7 October 2013

Academic Editor: Weimin Han

Copyright © 2013 Z. Kalateh Bojdi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We extend ordinary Padé approximation, which is based on a set of standard polynomials as , to a new extended Padé approximation (Müntz Padé approximation), based on the general basic function set    (the particular case of Müntz polynomials) using general Taylor series (based on fractional calculus) with error convergency. The importance of this extension is that the ordinary Padé approximation is a particular case of our extended Padé approximation. Also the parameterization ( is the corresponding parameter) of new extended Padé approximation is an important subject which, obtaining the optimal value of this parameter, can be a good question for a new research.

1. Introduction

Rational approximations of an arbitrary function are an important topic in numerical analysis due to their high applications in physical sciences, chemistry, engineering, and other applied sciences [1, 2]. The Padé approximation is a particular and classical type of rational fraction approximation. The idea of this approximation is to expand a function as a ratio of two power series and determining both the numerator and denominator coefficients using the coefficients of Taylor series expansion of a function [1]. The Padé approximation is the best approximation of a function by a rational function of a given order [1]. The technique was developed around 1890 by Henri Padé, but it goes back to George Freobenius who introduced the idea and investigated the features of rational approximations of power series. The Padé approximation is usually superior when functions contain poles, because the use of rational function allows them to be well represented [1]. The Padé approximation often gives better approximation of the function than truncating its Taylor series, and it may still work where the Taylor series does not converge [1]. For these reasons, Padé approximation is used extensively in computer calculations. The Padé approximation has also been used as an auxiliary function 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 it provides an approximation to the function throughout the whole complex plane, the study of Padé approximation is simultaneously a topic in mathematical approximation theory and analytic function theory. The generalized Padé approximation is given in [2]. For the connection of Padé approximation with continued fractions and orthogonal polynomials, see [2]. Also multivariate Padé approximation was done by [3]. Two versions of this approximation which concerned with interpolation problems are Hermite-Padé and Newton-Padé approximations and they are in fact rational interpolation problems [4–6]. In this paper, using the generalized Taylor series (based on fractional calculus), we extend ordinary Padé approximation to the general Müntz padé approximation. The paper is organized as follows. In Sections 1.1 to 1.3, we introduce an ordinary Padé approximation (and their properties), present the Müntz polynomials (rational Müntz functions and suitable theorems), and also generalize Taylor series (based on fractional calculus), respectively. Also Müntz Padé approximation is presented in Section 2. The uniqueness and convergence analysis results are given in Section 3. Some numerical tests are presented in Section 4. Two applications of this extended Padé approximation (Müntz Padé approximation) on the numerical solution of sound and vibration problem and microwave heating model in a slab using differential equations are given in Section 5. Finally we have monitored a brief conclusion in Section 6.

1.1. Ordinary Padé Approximation

Suppose that we are given a power series , represent a function , so that

This expansion is the fundamental starting point of any analysis using Padé approximation. A Padé approximation is a rational fraction as

which has a Maclaurin expansion which agrees with (1) as far as possible. Notice that, in (2), there are numerator coefficients and denominator coefficients. There is a more or less irrelevance common factor between them, and for definiteness we take . So there are independent numerator coefficients and independent denominator coefficients, making unknown coefficients in all. This number suggests that normally the out to fit the power series (1) through the orders . For the notation of formal power series, we have

and returning to (3) and cross-multiplying, we find that

Equating the coefficients of from (4), we find

If , we define for consistency. Since , (5) becomes a set of linear equations for the unknown denominator coefficients and we obtain

from which the may be found. The numerator coefficients follow immediately from (4) by equating the coefficients of as

Thus (4) and (7) normally determine the Padé numerator and denominator which are called Padé equations and we have constructed an Padé approximation which agrees with through order . If the given power series converge to the same function for with , then a sequence of Padé approximation may converge for , where is a domain longer than [1].

Remark 1. In general, Padé approximation does not exist (for an arbitrary function with respect to particular and ) [1].

Theorem 2. If the Padé approximation exists, then it will be unique [2].

1.2. Müntz Polynomials

An interesting generalization of Weierstrass theorem (1885) (density of polynomials set in , equipped with a supremum norm) goes back to the German mathematician Herman Müntz. He performed this generalization by a suitable tool which is named Müntz polynomials. For presenting his theorem, first we introduce a definition and a remark.

Definition 3. Suppose ; then, the set is named as a set of Müntz or Müntz polynomials [6].

Remark 4. Let be any sequences of distinct real numbers and ; then,

is dense in [6].

Theorem 5 (generalized Müntz theorem). If , then the Müntz polynomials with respect to these sequences are dense in , if and only if [6].

Remark 6. In this paper, by considering , which is a particular case of Müntz polynomials, we introduce a new and extended Müntz Padé approximation. Also, it is clear that by considering , we can obtain the classical Weierstrass theorem [6].

1.3. Generalized Taylor Series

In calculus, Taylor’s theorem gives an approximation of a -times differentiable function around a given point by a th order Taylor polynomial as

with the error

and for analytical polynomials, we have

Taylor’s theorem also generalizes to multivariate and vector valued functions. But the generalized version of this theorem from fractional calculus view or the theory of derivatives of arbitrary order was done by [7], because of its high applications in solving fractional ordinary differential equations integral ordinary differential equations and obtaining an operational matrix, and so forth. For this purpose, we need some definitions in fractional calculus.

Definition 7. The Caputo fractional derivative of of order with is defined as

for , , .

From Caputo fractional derivative, we have

where and . However, the Caputo fractional derivative is a linear operator which means

For more literature review of fractional calculus, see [8].

Theorem 8. Suppose that , for , where ; then, one has

where , for all , and , (-times) [7].

If we consider , then we obtain a fractional Maclaurin series and in a similar manner as the arbitrary function that has an infinite Caputo differentiable at point named as analytical functions on in a fractional sense if

2. Müntz Padé Approximations

Now by using the general Taylor series (based on fractional calculus), we generalize the classical and ordinary Padé approximations to a new and extended Müntz Padé approximation. For this goal, we present the Müntz Padé approximation as a ratio of two Müntz polynomials constructed from the coefficients of generalized Taylor series expansion of a function. Also in a similar approach from ordinary Padé approximation, we prove the uniqueness of Müntz Padé approximation. Now in a similar manner from ordinary Padé approximation, suppose that is an analytical function (from fractional calculus) in the neighborhood of , and then we can write

and also Müntz Padé approximation is defined as a ratio of two Müntz polynomials as

and then we suppose that (Müntz Padé approximation) has a Maclaurin series (17) or

By cross multiplying, we find that

Equating the coefficients , from (20) we obtain a system of equations similar to (5) and we can obtain coefficients. Also immediately by equating the coefficients of , we can obtain coefficients in a recursion formula as (7). Thus by assumption of existence of Müntz Padé approximation, for establishing the uniqueness of this approximation, we must prove that the matrix of system of (5) is nonsingular, but we perform the proof of the uniqueness of Müntz Padé approximation by a different approach in the next section.

Remark 9. Let ; then, we can obtain the ordinary Padé approximation.

3. Uniqueness and Convergence Analysis of Müntz Padé Approximation

In this section, we present the uniqueness and convergence analysis of Müntz Padé approximation in some theorems.

Theorem 10 (uniqueness). When Müntz Padé approximation exists, then it is unique for any formal power series.

Proof. Assume that there are two such Müntz Padé approximations and , where the degrees of and are less than or equal to and that of and are less than or equal to . Then by (19), we must have

since both approximate the same series. If we multiply (21) by , we obtain

but the left-hand side of (22) is a polynomial of degree at most and thus is identically zero. Since neither nor is identically zero, we conclude

Since, by definition, both and and and are relatively prime and ; we have shown that the two supposedly different Müntz Padé approximants are the same.

Theorem 11 (convergency). Let be analytic in (in fractional sense); then, an infinite subsequence of Müntz Padé approximants converges to uniformly in .

Proof. By hypothesis, is analytic in (in fractional sense) and consequently within a large interval, , with . Let

The Müntz Padé approximation is given by

If a subsequence of coefficients are zero, then are truncated fractional Maclaurin expansions which converge to uniformly in . So we assume that no infinite subsequence of vanishes and consider which is well defined for all sufficiently large , because

the sequence of Müntz Padé approximation given by (25) converges uniformly in unless there exists a sequence of values of for which within . Thus either a subsequence of the second row converges uniformly or else for some and all , . In the latter case

contradicting with (24), so the proof is completed.

4. The Test Experiments

Now in this section, in the two subsections, we show the applications of the Müntz Padé approximation in the functional approximation (see Section 4.1) and fractional calculus fields. The advantage of using the Müntz Padé approximation is shown for numerical approximation of fractional differential equations in Section 4.2.

4.1. Müntz Padé Approximation and Functional Approximation

In this section, we obtain the Müntz Padé approximations of , for different values of . Figure 1 shows the Müntz Padé approximations of and .

Figure 1: The Müntz Padé approximation of and exact function of .

Also the Müntz Padé approximation of , and is shown in Figure 2.

Figure 2: Comparison between the Müntz Padé approximation of and exact function of .

4.2. Müntz Padé Approximation and Fractional Calculus

This section is devoted to presentation of some numerical simulations obtained by applying the collocation method and based on a new extended Padé approximation (Müntz Padé approximation). The algorithm for numerical approximation of solutions to the initial value problems for the fractional differential equations was implemented by MATLAB. In the case of nonlinear equations, the MATLAB function was used for solving the nonlinear system. In the case that the exact solution to a problem is known, the dependence of approximation errors on the discretization parameter is estimated in 2-norm as

where is the approximated solution corresponding to the discretization parameter .

Experiment 1. We start with a simple nonlinear problem [8]

where we have a nonlinear and nonsmooth right-hand side. The solution has a smooth derivative of order [8]. The analytical solution subject to the initial condition is given by

Now we approximate the exact solution of (29) or , with the new and extended Padé approximation (Müntz Padé approximation) as

and substituting it into (29), we obtain

where

For the collocation points, we use the first roots of the Jacobi polynomial [9, 10], and then after enforcing the initial condition , we obtain a system of nonlinear algebraic equations and we use the MATLAB function for solving the nonlinear system. Thus, substituting the collocation points into (32) yields

Now from (34) and its initial condition, we have algebraic equations of unknown coefficients. Thus for obtaining the unknown coefficients, we must eliminate one arbitrary equation from these equations. But because of the necessity of holding the boundary conditions, we eliminate the last equation from (34). Finally, replacing the last equation of (34) by the equation of initial condition, we obtain a system of equations of unknowns . By implementing the method as presented, for and also for different parameters of and , we obtain the approximate solutions. The ordinary and Müntz Padé coefficients approximation of this problem are shown in Figures 3 and 4, respectively. We have observed that this method (the new extended Padé approximation (Müntz Padé approximation)) is very efficient for numerical approximation of the fractional ordinary differential equations. Also, closer look at the results of the Müntz Padé approximation scheme reveals that in this method of solution the coefficients decrease faster than the classical case. This is the advantage on the application of the Müntz Padé approximation to the fractional ordinary differential equations.

Figure 3: ordinary Padé approximation coefficients of , , and of Experiment 1.

Figure 4: new Padé approximation coefficients of , , , and of Experiment 1.

Experiment 2. Consider the nonlinear fractional integro-differential equation [11]

where

The exact solution of (35) is . We implement the collocation method based on the new Padé approximation of , , and . The obtained results of our method are presented in Table 1.

Table 1: The comparison between the collocation method based on Müntz Padé approximation and the exact solution of , , , of Experiment 2.

Exact
Approximate
Error

0.0
0.0000000000
0.0000000000
0.000000000

0.1
0.0953101798
0.0953103008
0.000000121

0.2
0.18232155679
0.1823218407
0.000000284

0.3
0.26236426446
0.2623649824
0.000000718

0.4
0.33647223662
0.3364726886
0.000000452

0.5
0.40546510810
0.40546528510
0.000000177

0.6
0.47000362924
0.4700043522
0.000000723

0.7
0.53062825106
0.5306290280
0.000000777

0.8
0.58778666490
0.58778755290
0.000000888

0.9
0.64185388617
0.64185410017
0.000000214

1.0
0.69314718055
0.69314789955
0.000000719

5. Application of Müntz Padé Approximation in Vibration and Electromagnetic Radiation Problems

In this section, we present two applications of Müntz Padé approximation for numerical approximation of some applicable ordinary differential equations.

Experiment 3. Consider the following model problem [12–15]:

where [12]. This model has high application in the theory of sound and vibration [12]. The exact solution to this initial value problem has the form

Now for , , and and using and Müntz Padé approximations, we obtain the following results shown in Table 2.

Also the obtained results for different kinds of Müntz Padé approximations are shown in Figure 5.

Table 2: Comparison between the and Müntz Padé approximations of , , and of Experiment 3.

Error of
Error of

0.0
0.0012
0.0009

0.1
0.0010
−0.0004

0.2
0.0013
0.0006

0.3
0.0022
−0.0002

0.4
0.0014
0.0007

0.5
0.0012
0.0003

Figure 5: A: the obtained results of Müntz Padé approximation of experiment 3. B: the obtained results of Müntz Padé approximation of Experiment 3.

Experiment 4. Consider the following model problem [16]:

with

This equation has high application in the theory of electromagnetic radiation and describes the steady state reaction-diffusion equations with source term that arise in modeling microwave heating in an infinite slab with isothermal walls [16]. Also and represent the thermal absorptivity and electric field decay rate parameters, respectively. The exact solution to this boundary value problem has the form

where , are Bessel functions of the first kind and , are Bessel functions of the second kind. Now for , and and using and Müntz Padé approximations, we obtain the following results shown in Table 3.

Table 3: Comparison between the and Müntz Padé approximations of and of Experiment 4.

Error of
Error of

0.0
0.0002
0.0001

0.1
−0.0033
−0.0011

0.2
0.0014
0.0009

0.3
0.0003
0.0001

0.4
0.0026
−0.0012

0.5
0.0023
0.0021

6. Conclusion

In this paper, using the general Taylor series (based on fractional calculus), we extend the ordinary Padé approximation to the general Müntz Padé approximation. The importance of this extension is that the ordinary Padé approximation is a particular case of our Müntz Padé approximation (). We have applied the method in the application of functional approximation, fractional exponent, and vibration and electromagnetic radiation model problems and have obtained the results with a good order of accuracy. Also the uniqueness results and error analysis have been presented completely. In addition, the test experiments have been presented for showing the applicability and validity of the new Müntz Padé approximation.

Acknowledgments

The authors are grateful to the reviewers and the editor Professor Dr. Weimin Han for their helpful comments and suggestions which indeed improved the quality of this paper.

References

G. A. Baker and P. Graves-Morris, Padé Approximants, Addison-Wesley, 1981.

G. A. Baker, Essentials of Padé Approximants, Academic Press, New York, NY, USA, 1975.

A. Cuyt, “Multivariate Padé approximants revisited,” BIT, vol. 26, no. 1, pp. 71–79, 1986.

L. Wuytack, “On the osculatory rational interpolation problem,” Mathematics and Computers in Simulation, vol. 29, pp. 837–843, 1975.

A. A. M. Cuyt and B. M. Verdonk, “General order Newton-Padé approximants for multivariate functions,” Numerische Mathematik, vol. 43, no. 2, pp. 293–307, 1984.

P. Borwein and T. Erdélyi, Polynomials and Polynomials Inequalities, vol. 161 of Graduate Texts in Mathematics, Springer, Berlin, Germany, 1995.

Z. M. Odibat and N. T. Shawagfeh, “Generalized Taylor's formula,” Applied Mathematics and Computation, vol. 186, no. 1, pp. 286–293, 2007.

K. Diethelm, The Analysis of Fractional Differential Equations, Springer, Berlin, Germany, 2010.

D. Funaro, Polynomial Approximations of Differential Equations, Springer, 1992.

J. S. Hesthaven, S. Gottlieb, and D. Gottlieb, Spectral Methods for Time-Dependent Problems, Cambridge University, 2009.

P. Mokhtary and F. Ghoreishi, “The -convergence of the Legendre spectral tau matrix formulation for nonlinear fractional integro differential equations,” Numerical Algorithms, vol. 58, no. 4, pp. 475–496, 2011.

I. V. Andrianov and J. Awrejcewicz, “Analysis of jump phenomena using Padé approximations,” Journal of Sound and Vibration, vol. 260, no. 3, pp. 577–588, 2003.

I. Andrianov and J. Awrejcewicz, “Solutions in the Fourier series form, Gibbs phenomena and Padé approximants,” Journal of Sound and Vibration, vol. 245, no. 4, pp. 753–756, 2001.

I. V. Andrianov and J. Awrejcewicz, “Iterative determination of homoclinic orbit parameters and Padé approximants,” Journal of Sound and Vibration, vol. 240, no. 2, pp. 394–397, 2001.

G. Kudra and J. Awrejcewicz, “Tangents hyperbolic approximations of the spatial model of friction coupled with rolling resistance,” International Journal of Bifurcation and Chaos, vol. 21, pp. 2905–2917, 2011.

O. D. Makinde, “Solving microwave heating model in a slab using Hermite-Padé approximation technique,” Applied Thermal Engineering, vol. 27, pp. 599–603, 2007.