On bi-orthogonal systems of trigonometric functions and quadrature formulas for periodic integrands

In this paper, quadrature formulas with an arbitrary number of nodes and exactly integrating trigonometric polynomials up to degree as high as possible are constructed in order to approximate 2π-periodic weighted integrals. For this purpose, certain bi-orthogonal systems of trigonometric functions are introduced and their most relevant properties studied. Some illustrative numerical examples are also given. The paper completes the results previously given by Szeg? in Magy Tud Akad Mat Kut Intez Közl 8:255–273, 1963 and by some of the authors in Annales Mathematicae et Informaticae 32:5–44, 2005.     

Approximating exponential and logarithmic functions using polynomial interpolation

This article takes a closer look at the problem of approximating the exponential and logarithmic functions using polynomials. Either as an alternative to or a precursor to Taylor polynomial approximations at the precalculus level, interpolating polynomials are considered. A measure of error is given and the behaviour of the error function is analysed. The results of interpolating polynomials are compared with those of Taylor polynomials.     

The Bernoulli Spline and Approximation by Trigonometric Blending Polynomials

Using an extension of the notion of Bernoulli spline to a multivariate setting, a Jackson-Favard estimate is derived for approximation of 2π-periodic test functions by trigonometric blending polynomials. Techniques involved in the proof include properties of periodic distributions and of their Fourier transforms. The usual Jackson-Favard estimates from the literature can be derived from our result if limits of test functions are considered.     

Operator adapted spectral element methods I: harmonic and generalized harmonic polynomials

Summary. The paper presents results on the approximation of functions which solve an elliptic differential equation by operator adapted systems of functions. Compared with standard polynomials, these operator adapted systems have superior local approximation properties. First, the case of Laplace's equation and harmonic polynomials as operator adapted functions is analyzed and rates of convergence in a Sobolev space setting are given for the approximation with harmonic polynomials. Special attention is paid to the approximation of singular functions that arise typically in corners. These results for harmonic polynomials are extended to general elliptic equations with analytic coefficients by means of the theory of Bergman and Vekua; the approximation results for Laplace's equation hold true verbatim, if harmonic polynomials are replaced with generalized harmonic polynomials. The Partition of Unity Method is used in a numerical example to construct an operator adapted spectral method for Laplace's equation that is based on approximating with harmonic polynomials locally. Received May 26, 1997 / Revised version received September 21, 1998 / Published online September 7, 1999     

Fourier and barycentric formulae for equidistant Hermite trigonometric interpolation

We consider the Hermite trigonometric interpolation problem of order 1 for equidistant nodes, i.e., the problem of finding a trigonometric polynomial t that interpolates the values of a function and of its derivative at equidistant points. We give a formula for the Fourier coefficients of t in terms of those of the two classical trigonometric polynomials interpolating the values and those of the derivative separately. This formula yields the coefficients with a single FFT. It also gives an aliasing formula for the error in the coefficients which, on its turn, yields error bounds and convergence results for differentiable as well as analytic functions. We then consider the Lagrangian formula and eliminate the unstable factor by switching to the barycentric formula. We also give simplified formulae for even and odd functions, as well as consequent formulae for Hermite interpolation between Chebyshev points.     

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets is presented for solving 2-D potential problems defined inside or outside of a circular boundary in this paper. In this approach, an equivalent variational form of the corresponding boundary integral equation for the potential problem is used; the trigonometric wavelets are employed as trial and test functions of the variational formulation. The analytical formulae of the matrix entries indicate that most of the matrix entries are naturally zero without any truncation technique and the system matrix is a block diagonal matrix. Each block consists of four circular submatrices. Hence the memory spaces and computational complexity of the system matrix are linear scale. This approach could be easily coupled into domain decomposition method based on variational formulation. Finally, the error estimates of the approximation solutions are given and some test examples are presented.     

A multiscale Galerkin method for the hypersingular integral equation reduced by the harmonic equation

The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.     

Statistical Approximation by Double Poisson–Cauchy Singular Integral Operators

In this paper we show that it is possible to approximate a continuous and 2π-periodic function on the disk centered at origin with radius π by means of double Poisson–Cauchy singular integral operators which do not need to be positive in general. Our results cover not only the classical approximation but also the statistical approximation process.     

Approximation of Periodic Functions of High Smoothness by Interpolation Trigonometric Polynomials in the Metric of L1

We establish an asymptotically exact estimate for the error of approximation of ²-periodic functions of high smoothness by interpolation trigonometric polynomials in the metric of L ₁.     

ϕ-Strong approximation of functions by trigonometric polynomials

Mathematical Notes - The rate of ?-strong approximation of periodic functions by trigonometric polynomials constructed on the basis of interpolating polynomials with equidistant nodes is...     

Trigonometric wavelet method for some elliptic boundary value problems

In this paper, we apply trigonometric wavelet method to investigate the numerical solution of elliptic boundary value problem in an exterior unit disk. The simple computation formulae of the entries in the stiffness matrix are obtained. It shows that we only need to compute 2(j²+1) elements of one ²j+2ײj+2 stiffness matrix. Moreover, the error estimates of the approximation solutions are given and some test examples are presented.     

Approximate solutions and error estimates for a Stokes boundary value problem

The aim of this paper is the numerical treatment of a boundary value problem for the system of Stokes’ equations. For this we extend the method of approximate approximations to boundary value problems. This method was introduced by Maz’ya (DFG-Kolloquium des DFG-Forschungsschwerpunktes Randelementmethoden, 1991) and has been used until now for the approximation of smooth functions defined on the whole space and for the approximation of volume potentials. In the present paper we develop an approximation procedure for the solution of the interior Dirichlet problem for the system of Stokes’ equations in two dimensions. The procedure is based on potential theoretical considerations in connection with a boundary integral equations method and consists of three approximation steps as follows. In a first step the unknown source density in the potential representation of the solution is replaced by approximate approximations. In a second step the decay behavior of the generating functions is used to gain a suitable approximation for the potential kernel, and in a third step Nyström’s method leads to a linear algebraic system for the approximate source density. For every step a convergence analysis is established and corresponding error estimates are given.     

A SPECIAL MODULUS OF CONTINUITY AND THE K-FUNCTIONAL

We consider the questions connected with the approximation of a real continuous 1-periodic functions and give a new proof of the equivalence of the special Boman-Shapiro modulus of continuity with Peetre's K-functional. We also prove Jackson's inequality for the approximation by trigonometric polynomials.     

An approximation for the nonlinear filtering problem,with error bound †

For the standard continuous-time nonlinear filtering problem an approximation approach is derived. The approximate filter is given by the solution to an appropriate discrete-time approximating filtering problem that can be explicitly solved by a finite-dimensional procedure. Furthermore an explicit upper bound for the approximation error is derived. The approximating problem is obtained by first approximating the signal and then using measure transformation to express the original observation process in terms of the approximating signal     

Defect correction from a Galerkin viewpoint

Summary We consider the numerical solution of systems of nonlinear two point boundary value problems by Galerkin's method. An initial solution is computed with piecewise linear approximating functions and this is then improved by using higher-order piecewise polynomials to compute defect corrections. This technique, including numerical integration, is justified by typical Galerkin arguments and properties of piecewise polynomials rather than the traditional asymptotic error expansions of finite difference methods.     

GUARANTEED BOUNDS FOR THE SOLUTION OF THE WAVE EQUATION

Abstract

Numerical solution of the wave equation in the form of close lower and upper bounds provides a secure a posteriori error estimate that can be used for efficient accuracy control. The method considered in this paper uses some monotone properties of the differential operator in the wave equation to construct bounds for the solution in the form of trigonometric polynomials of x. Aspects of the numerical implementation, the accuracy of the computed bounds and some numerical examples are discussed.     

Best approximations of classes of functions defined by means of the modulus of continuity

This paper is devoted to an exact solution of problems of best approximation in the uniform and integral metrics of classes of periodic functions representable as a convolution of a kernel not increasing the oscillation with functions having a given convex upwards majorant of the modulus of continuity. The approximating sets are taken to be the trigonometric polynomials in the case of the uniform and integral metrics, and convolutions of the kernel defining the class with polynomial splines in the case of the integral metric.Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 44, No. 5, pp. 579–589, May, 1992.     

Sharp inequalities for trigonometric polynomials with respect to integral functionals

The problem on sharp inequalities for linear operators on the set of trigonometric polynomials with respect to integral functionals \(\int_0^{2\pi } {\phi \left( {\left| {f\left( x \right)} \right|} \right)dx}\) is discussed. A solution of the problem on trigonometric polynomials with given leading harmonic of least deviation from zero with respect to such functionals over the set of all functions φ defined, nonnegative, and nondecreasing on the semiaxis [0,+∞) is given.     

Interpolation Methods for the Evaluation of 2π-Periodic Finite Baire Measure

We discuss the definition and effectiveness of a Padé-type approximation to 2π-periodic finite Baire measures on [-π,π]. In the first two sections we recall the definitions and basic properties of the Padré-type approximants to harmonic functions in the unit disk and to L ^p -functions on the unit circle. Section 3 deals with the extension of these definitions and properties to a finite 2π-periodic Baire measure. Finally, section 4 is devoted to a study of the convergence of a sequence of such approximants, in the weak star topology of measures.     

Numerical integration in Galerkin meshless methods, applied to elliptic Neumann problem with non-constant coefficients

In this paper, we explore the effect of numerical integration on the Galerkin meshless method used to approximate the solution of an elliptic partial differential equation with non-constant coefficients with Neumann boundary conditions. We considered Galerkin meshless methods with shape functions that reproduce polynomials of degree k?≥?1. We have obtained an estimate for the energy norm of the error in the approximate solution under the presence of numerical integration. This result has been established under the assumption that the numerical integration rule satisfies a certain discrete Green’s formula, which is not problem dependent, i.e., does not depend on the non-constant coefficients of the problem. We have also derived numerical integration rules satisfying the discrete Green’s formula.