Optimization of two-point Padé-type approximants

In this paper, we first give characterization theorems for the best two-point Padé-type approximants (2PTAs) in the uniform norm. Secondly, we consider sequences of 2PTAs in a domain of the complex plane from the viewpoint of the asymptotic degree of convergence, and we also give conditions for geometric convergence.     

A note on the determinant formulas computation of generalized inverse matrix Padé approximation

In this paper, the computation of two special determinants which appear in the construction of a generalized inverse matrix Padé approximation of type [n/2k] (described in [Linear Algebra Appl. 322 (2001) 141]) for a given power series is investigated. Here a common computational approach of determinant can not be used. The main tool to be used to do the two special determinants is the well-known Schur complement theorem.     

Padé iteration method for regularization

In this study we present iterative regularization methods using rational approximations, in particular, Padé approximants, which work well for ill-posed problems. We prove that the (k, j)-Padé method is a convergent and order optimal iterative regularization method in using the discrepancy principle of Morozov. Furthermore, we present a hybrid Padé method, compare it with other well-known methods and found that it is faster than the Landweber method. It is worth mentioning that this study is a completion of the paper [A. Kirsche, C. Böckmann, Rational approximations for ill-conditioned equation systems, Appl. Math. Comput. 171 (2005) 385–397] where this method was treated to solve ill-conditioned equation systems.     

Padé approximation for a multivariate Markov transform

Methods of Padé approximation are used to analyse a multivariate Markov transform which has been recently introduced by the authors. The first main result is a characterization of the rationality of the Markov transform via Hankel determinants. The second main result is a cubature formula for a special class of measures.     

Hermite–Padé approximation and simultaneous quadrature formulas

We study Hermite–Padé approximation of the so-called Nikishin systems of functions. In particular, the set of multi-indices for which normality is known to take place is considerably enlarged as well as the sequences of multi-indices for which convergence of the corresponding simultaneous rational approximants takes place. These results are applied to the study of the convergence properties of simultaneous quadrature rules of a given function with respect to different weights.     

Discrete non-local boundary conditions for split-step Padé approximations of the one-way Helmholtz equation

This paper deals with the efficient numerical solution of the two-dimensional one-way Helmholtz equation posed on an unbounded domain. In this case, one has to introduce artificial boundary conditions to confine the computational domain. The main topic of this work is the construction of the so-called discrete transparent boundary conditions for state-of-the-art parabolic equation methods, namely a split-step discretization of the high-order parabolic approximation and the split-step Padé algorithm of Collins. Finally, several numerical examples arising in optics and underwater acoustics illustrate the efficiency and accuracy of our approach.     

Theoretical and numerical considerations about Padé approximants for the matrix logarithm

We show that for a vast class of matrix Lie groups, which includes the orthogonal and the symplectic, diagonal Padé approximants of log((1+x)/(1−x)) are structure preserving. The conditioning of these approximants is analyzed. We also present a new algorithm for the Briggs–Padé method, based on a strategy for reducing the number of square roots in the inverse scaling and squaring procedure.     

Multipoint Padé approximants to complex Cauchy transforms with polar singularities

We study diagonal multipoint Padé approximants to functions of the form

where R is a rational function and λ is a complex measure with compact regular support included in , whose argument has bounded variation on the support. Assuming that interpolation sets are such that their normalized counting measures converge sufficiently fast in the weak-star sense to some conjugate-symmetric distribution σ, we show that the counting measures of poles of the approximants converge to , the balayage of σ onto the support of λ, in the weak^* sense, that the approximants themselves converge in capacity to F outside the support of λ, and that the poles of R attract at least as many poles of the approximants as their multiplicity and not much more.     

Convergence of Multipoint Padé-type Approximants

Let μ be a finite positive Borel measure whose support is a compact subset K of the real line and let I be the convex hull of K. Let r denote a rational function with real coefficients whose poles lie in C\I and r(∞)=0. We consider multipoint rational interpolants of the function where some poles are fixed and others are left free. We show that if the interpolation points and the fixed poles are chosen conveniently then the sequence of multipoint rational approximants converges geometrically to f in the chordal metric on compact subsets of C\I.     

Padé approximation of Stieltjes series

Using Nuttall's compact formula for the [n, n − 1] Pad'e approximant, the authors show that there is a natural connection between the Padé approximants of a series of Stieltjes and orthogonal polynomials. In particular, we obtain the precise error formulas. The [n, n − 1] Padé approximant in this case is just a Gaussian quadrature of the Stieltjes integral. Hence, analysis of the error is now possible and under very mild conditions it is shown that the [n, n + j], j − 1, Padé approximants converge to the Stieltjes integral.     

Moment methods in two point Padé approximation

In the separable Hilbert space (H, ·, ·) the following “operator moment problem” is solved: given a complex sequence (c_k)_{k ε Z} generated by a meromorphic function f, find T ε B(H) and u₀ ε H such that T^ku₀, u₀ = c_k (k ε Z). If the sequence (c_k)_{k ε Z} is “normal,” an adapted form of Vorobyev's method of moments yields a sequence of two point Padé approximants to f. A sufficient condition for convergence of this sequence of approximants is given.     

B-splines and Hermite–Padé approximants to the exponential function

This paper is the continuation of a work initiated in [P. Sablonnière, An algorithm for the computation of Hermite–Padé approximations to the exponential function: divided differences and Hermite–Padé forms. Numer. Algorithms 33 (2003) 443–452] about the computation of Hermite–Padé forms (HPF) and associated Hermite–Padé approximants (HPA) to the exponential function. We present an alternative algorithm for their computation, based on the representation of HPF in terms of integral remainders with B-splines as Peano kernels. Using the good properties of discrete B-splines, this algorithm gives rise to a great variety of representations of HPF of higher orders in terms of HPF of lower orders, and in particular of classical Padé forms. We give some examples illustrating this algorithm, in particular, another way of constructing quadratic HPF already described by different authors. Finally, we briefly study a family of cubic HPF.     

On the convergence of bivariate two-point Padé-type approximants

Some choices of denominators are given which ensure the geometrical convergence of certain convergence of bivariate two-point Padé-type approximants to functions being holomorphic on certain domains     

A numerical solution of the elasticity problem of settled of the wronkler ground with variable coefficients

In this paper, we have given numerical solution of the elasticity problem of settled on the wronkler ground with variable coefficient. The approximation solution of boundary value problem which is pertinent to this has been converted to integral equations, and then by using the successive approximation methods, has been reached. In addition to this, the approximation solution of the problem was put into Padé series form. We applied these methods to an example which is the elasticity problem of unit length homogeny beam, which is a special form of boundary value problem. First we calculate the successive approximation of the given boundary value problem then transform it into Padé series form, which give an arbitrary order for solving differential equation numerically.     

Rational approximation with sign-sensitive weight

For metrics of sign-sensitive weight, a generalization of the Jackson theorem and an analog of the Dolzhenko theorem on the estimate of the variation of a function in terms of its least rational deviations are obtained. In the same metric, two-sided estimates for rational deviations of the function signx are given.Translated fromMatematickeskie Zametki, Vol. 60, No. 5, pp. 715–725, November, 1996.     

Analytic approximation of the blow‐up time for nonlinear differential equations by the ADM–Padé technique

We present a new approach to calculate analytic approximations of blow‐up solutions and their critical blow‐up times. Our approach applies the Adomian decomposition–Padé method to quickly and easily compute the critical blow‐up times, which comprises the Adomian decomposition method combined with the Padé approximants technique. We validate our new approach with a variety of numerical examples, including nonlinear ODEs, systems of nonlinear ODEs, and nonlinear PDEs. Furthermore, our new method is shown to be more convenient than prior art that relies on compound discretized algorithms. Copyright © 2013 John Wiley & Sons, Ltd.     

Shear stress analysis in axially symmetric flow using the Padé approximation technique

This paper evaluates the shear stress on the surface of a circular cylinder in axially symmetric flow within the boundary layer. An extended solution is obtained by established techniques, e.g. asymptotic series or Pohlhausen method. A Padé approximation technique is applied to the extended solution obtained. The region of validity of the solution is extended and the results are then compared and examined against known approximate Pohlhausen and series solutions. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Quadratic Hermite–Padé polynomials associated with the exponential function

The asymptotic behavior of quadratic Hermite–Padé polynomials associated with the exponential function is studied for n→∞. These polynomials are defined by the relation

(*)

p_n(z)+q_n(z)e^z+r_n(z)e^2z=O(z³ⁿ⁺²) as z→0,

where O(·) denotes Landau's symbol. In the investigation analytic expressions are proved for the asymptotics of the polynomials, for the asymptotics of the remainder term in (*), and also for the arcs on which the zeros of the polynomials and of the remainder term cluster if the independent variable z is rescaled in an appropriate way. The asymptotic expressions are defined with the help of an algebraic function of third degree and its associated Riemann surface. Among other possible applications, the results form the basis for the investigation of the convergence of quadratic Hermite–Padé approximants, which will be done in a follow-up paper.     

A Padé family of iterations for the matrix sign function and related problems

In this paper we consider the Pad'e family of iterations for computing the matrix sign function and the Padé family of iterations for computing the matrix p‐sector function. We prove that all the iterations of the Padé family for the matrix sign function have a common convergence region. It completes a similar result of Kenney and Laub for half of the Padé family. We show that the iterations of the Padé family for the matrix p‐sector function are well defined in an analogous common region, depending on p. For this purpose we proved that the Padé approximants to the function (1?z)^?σ, 0<σ<1, are a quotient of hypergeometric functions whose poles we have localized. Furthermore we proved that the coefficients of the power expansion of a certain analytic function form a positive sequence and in a special case this sequence has the log‐concavity property. Copyright © 2011 John Wiley & Sons, Ltd.     

Copositive approximation by rational functions with prescribed numerator degree

The paper proves that, if f（x） ∈ L^p[-1,1],1≤p〈∞ ,changes sign I times in （-1, 1）,then there exists a real rational function r（x） ∈ Rn^（2μ-1）l which is eopositive with f（x）, such that the following Jackson type estimate ｜｜f-r｜｜p≤Cδl^2μωφ（f,1/n）p holds, where μ is a natural number ≥3/2＋1/p, and Cδ is a positive constant depending only on δ.