An improved estimate of PSWF approximation and approximation by Mathieu functions

In this paper, an error estimate of spectral approximations by prolate spheroidal wave functions (PSWFs) with explicit dependence on the bandwidth parameter and optimal order of convergence is derived, which improves the existing result in [Chen et al., Spectral methods based on prolate spheroidal wave functions for hyperbolic PDEs, SIAM J. Numer. Anal. 43 (5) (2005) 1912-1933]. The underlying argument is applied to analyze spectral approximations of periodic functions by Mathieu functions, which leads to new estimates featured with explicit dependence on the intrinsic parameter.     

Power series expansions for spheroidal wave functions with small arguments

Power series expansions for the angular spheroidal wave functions of the first kind S_mn(c,η), with small arguments c, are derived for general integer values of m and n. The various evaluated expansion coefficients can also be used in the calculation of the corresponding angular functions of the second kind, as well as for the radial functions of any kind. Only the prolate functions are considered explicitly, but corresponding formulas for the oblate ones are obtained immediately.     

Minimax approximation by rational fractions of the inverse polynomial type

Rational fractions of the formR(x)/(c ₁ +c ₂ x +c ₃ x ² + ...) ^r are useful for approximating decay type functions over infinite and semi-infinite domains. A procedure is given which produces the optimal coefficients with no more effort than for linear approximations. No initial guess is needed for the values of the coefficients nor for the maximum error of approximation.     

The Eigenvalues of the Angular Spheroidal Wave Equation

Approximate relations are obtained between the eigenvalues λ and the ellipticity parameter c² of the angular spheroidal wave equation. Although based on WKBJ methods and the assumption that λ is large, the relations are useful throughout the complex c²-plane. They are exact at c² = 0, and reproduce the standard asymptotic formulas for λ when c² is large. At intermediate values of c², they provide approximations for the square-root branch points of the multivalued function λ(c²) in the complex c²-plane at which adjacent eigenvalues of the same class become equal in pairs. These branch points lie on an infinite sequence of distorted circular rings. Their exact locations have been computed for the first four rings for angular wavenumbers m = 0,…,4.     

A weak Galerkin mixed finite element method for the Helmholtz equation with large wave numbers

In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H¹ and L² norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.     

Nonstationary queues with Interrrupted Poisson arrivals and unreliable/repairable servers

A queueing model having a nonstationary Interrupted Poisson arrival process (IPP(t)),s time-dependent exponential unreliable/repairable servers and finite capacityc is introduced, and an approximation method for analysis of it is developed and tested. Approximations are developed for the time-dependent queue length moments and the system viewpoint waiting time distributions and moments. The approximation involves state-space partitioning and numerically integrating partial-moment differential equations (PMDEs). Surrogate distribution approximations (SDA's) are used to close the system of PMDEs. The approximations allow for analysis using only (s + 1)(s + 6) differential equations for the queue length moments rather than the 2(c + 1)(s +1) equations required by the classic method of numerically integrating the full set of Kolmogorov-forward equations. Effectively hours of cpu time are reduced to minutes for even modest capacity systems. Approximations for waiting time distributions and moments are developed.This research was partially funded by National Science Foundation grant ECS-8404409.     

Formulation of Quantum Scattering Theory in Terms of Proper Differentials (Stationary Wave Packets)

Constructing the basic operators of scattering theory on and off the mass shell in terms of spatially bounded stationary wave packets or proper differentials is described. For this, we use a technique based on a certain scheme for discretizing the continuum. Finite-dimensional approximations for the Green's functions and T-matrix, which are first found here, are immediately constructed for any energy using a single simple diagonalization of the Hamiltonian matrix in an L ₂-type complete basis. We show that the developed approach leads to a convenient finite-dimensional representation of the scattering operators in the basis of the wave functions of a harmonic oscillator. The method allows an immediate extension to the problem of three and more bodies.     

Best simultaneous L1 approximations

Three possible definitions are proposed for best simultaneous L₁ approximation to n continuous real-valued functions, and the relation between best simultaneous approximations and best L₁ approximations to the arithmetic mean of the n functions is discussed.     

Near-best L1 approximations on circular and elliptical contours

Polynomial approximations are obtained to analytic functions on circular and elliptical contours by forming partial sums of order n of their expansions in Taylor series and Chebyshev series of the second kind, respectively. It is proved that the resulting approximations converge in the L₁ norm as n → ∞, and that they are near-best L₁ approximations within relative distances of the order of log n. Practical implications of the results are discussed, and they are shown to provide a theoretical basis for polynomial approximation methods for the evaluation of indefinite integrals on contours.     

On the Prolate Spheroidal Wave Functions and Hardy’s Uncertainty Principle

We prove a weak version of Hardy’s uncertainty principle using properties of the prolate spheroidal wave functions. We describe the eigenvalues of the sum of a time limiting operator and a band limiting operator acting on \(L^2(\mathbb {R})\) . A weak version of Hardy’s uncertainty principle follows from the asymptotic behavior of the largest eigenvalue as the time limit and the band limit approach infinity. An asymptotic formula for this eigenvalue is obtained from its well-known counterpart for the prolate integral operator.     

Wavelets Based on Prolate Spheroidal Wave Functions

The article is concerned with a particular multiresolution analysis (MRA) composed of Paley–Wiener spaces. Their usual wavelet basis consisting of sinc functions is replaced by one based on prolate spheroidal wave functions (PSWFs) which have much better time localization than the sinc function. The new wavelets preserve the high energy concentration in both the time and frequency domain inherited from PSWFs. Since the size of the energy concentration interval of PSWFs is one of the most important parameters in some applications, we modify the wavelets at different scales to retain a constant energy concentration interval. This requires a slight modification of the dilation relations, but leads to locally positive kernels. Convergence and other related properties, such as Gibbs phenomenon, of the associated approximations are discussed. A computationally friendly sampling technique is exploited to calculate the expansion coefficients. Several numerical examples are provided to illustrate the theory.     

On the evaluation of prolate spheroidal wave functions and associated quadrature rules

As demonstrated by Slepian et al. in a sequence of classical papers (see Slepian (1983) [33], Slepian and Pollak (1961) [34], Landau and Pollak (1961) [18], Slepian and Pollak (1964) [35], Slepian (1965) [36]), prolate spheroidal wave functions (PSWFs) provide a natural and efficient tool for computing with bandlimited functions defined on an interval. Recently, PSWFs have been becoming increasingly popular in various areas in which such functions occur – this includes physics (e.g. wave phenomena, fluid dynamics), engineering (signal processing, filter design), etc.To use PSWFs as a computational tool, one needs fast and accurate numerical algorithms for the evaluation of PSWFs and related quantities, as well as for the construction of corresponding quadrature rules, interpolation formulas, etc. During the last 15 years, substantial progress has been made in the design of such algorithms – see, for example, Xiao et al. (2001) [40] (see also Bowkamp (1947) [6], Slepian and Pollak (1961) [34], Landau and Pollak (1961) [18], Slepian and Pollak (1964) [35] for some classical results).The complexity of many of the existing algorithms, however, is at least quadratic in the band limit c. For example, the evaluation of the nth eigenvalue of the prolate integral operator requires $O(c^2+n^2)$ operations (see e.g. Xiao et al. (2001) [40]); the construction of accurate quadrature rules for the integration (and associated interpolation) of bandlimited functions with band limit c requires $O(c^3)$ operations (see e.g. Cheng et al. (1999) [8]). Therefore, while the existing algorithms are satisfactory for moderate values of c (e.g. $c⩽{10}^3$), they tend to be relatively slow when c is large (e.g. $c⩾{10}^4$).In this paper, we describe several numerical algorithms for the evaluation of PSWFs and related quantities, and design a class of PSWF-based quadratures for the integration of bandlimited functions. While the analysis is somewhat involved and will be published separately (currently, it can be found in Osipov and Rokhlin (2012) [27]), the resulting numerical algorithms are quite simple and efficient in practice. For example, the evaluation of the nth eigenvalue of the prolate integral operator requires $O(n+c\cdot \log c)$ operations; the construction of accurate quadrature rules for the integration (and associated interpolation) of bandlimited functions with band limit c requires $O(c)$ operations. All algorithms described in this paper produce results essentially to machine precision. Our results are illustrated via several numerical experiments.     

Existence, uniqueness and asymptotic stability of time periodic traveling waves for a periodic Lotka–Volterra competition system with diffusion

We study the existence, uniqueness, and asymptotic stability of time periodic traveling wave solutions to a periodic diffusive Lotka–Volterra competition system. Under certain conditions, we prove that there exists a maximal wave speed c^? such that for each wave speed c?c^?, there is a time periodic traveling wave connecting two semi-trivial periodic solutions of the corresponding kinetic system. It is shown that such a traveling wave is unique modulo translation and is monotone with respect to its co-moving frame coordinate. We also show that the traveling wave solutions with wave speed c<c^? are asymptotically stable in certain sense. In addition, we establish the nonexistence of time periodic traveling waves for nonzero speed c>c^?.     

THE ASYMPTOTIC THEORY OF INITIAL VALUEPROBLEMS FOR SEMILINEAR WAVE EQUATIONS IN THREE SPACE DIMENSIONS

This paper deals with the asymptotic theory of initial value problems for semilinear waveequations in three space dimensions. The well-posedness and validity of formal approximations ona long time scale of order |ε|^-1 are discussed in the classical sense of C^2 This result describes aceu-ratively the approximations of solutions. At the end of this paper an application of the asymptotictheory is given to analyze a special model for a perturbed wave equation,     

Traveling waves for a four-compartment lattice epidemic system with exposed class and standard incidence

This paper is considering the problem of traveling wave solutions (TWS) for a susceptible-exposed-infectious-recovered (SEIR) epidemic model with discrete diffusion. The threshold condition for the existence and nonexistence of TWS is obtained. More specifically, such kind of solutions are governed by the threshold number ?₀. We can find a critical wave speed c^? if ?₀ > 1, by employing the Schauder's fixed point theorem, limiting argument and two-sided Laplace transform, we confirm that there exists TWS for c > c^?, while there exists no TWS for c < c^?. We also obtain the nonexistence of TWS for ?₀ ≤ 1. At last, we give some biological explanations from the epidemiological perspective.     

Generalized Prolate Spheroidal Wave Functions: Spectral Analysis and Approximation of Almost Band-Limited Functions

In this work, we first give various explicit and local estimates of the eigenfunctions of a perturbed Jacobi differential operator. These eigenfunctions generalize the famous classical prolate spheroidal wave functions (PSWFs), founded in 1960s by Slepian and his co-authors and corresponding to the case \(\alpha =\beta =0.\) They also generalize the new PSWFs introduced and studied recently in Wang and Zhang (Appl Comput Harmon Anal 29:303–329, 2010), denoted by GPSWFs and corresponding to the case \(\alpha =\beta .\) The main content of this work is devoted to the previous interesting special case \(\alpha =\beta >- 1.\) In particular, we give further computational improvements, as well as some useful explicit and local estimates of the GPSWFs. More importantly, by using the concept of a restricted Paley–Wiener space, we relate the GPSWFs to the solutions of a generalized energy maximisation problem. As a consequence, many desirable spectral properties of the self-adjoint compact integral operator associated with the GPSWFs are deduced from the rich literature of the PSWFs. In particular, we show that the GPSWFs are well adapted for the spectral approximation of the classical c-band-limited as well as almost c-band-limited functions. Finally, we provide the reader with some numerical examples that illustrate the different results of this work.     

Spectral analysis of the finite Hankel transform and circular prolate spheroidal wave functions

In this paper, we develop two practical methods for the computation of the eigenvalues as well as the eigenfunctions of the finite Hankel transform operator. These different eigenfunctions are called circular prolate spheroidal wave functions (CPSWFs). This work is motivated by the potential applications of the CPSWFs as well as the development of practical methods for computing their values. Also, in this work, we should prove that the CPSWFs form an orthonormal basis of the space of Hankel band-limited functions, an orthogonal basis of L²([0,1]) and an orthonormal system of L²([0,+∞[). Our computation of the CPSWFs and their associated eigenvalues is done by the use of two different methods. The first method is based on a suitable matrix representation of the finite Hankel transform operator. The second method is based on the use of an efficient quadrature method based on a special family of orthogonal polynomials. Also, we give two Maple programs that implement the previous two methods. Finally, we present some numerical results that illustrate the results of this work.     

Constrained Lp-Approximation by Generalized n-Convex Functions Induced by ECT-Systems

The problem of finding a best L_p-approximation (1 ≤ p < ∞) to a function in L_p from a special subcone of generalized n-convex functions induced by an ECT-system is considered. Tchebycheff splines with a countably infinite number of knots are introduced and best approximations are characterized in terms of local best approximations by these splines. Various properties of best approximations and their uniqueness in L₁ are investigated. Some special results for generalized monotone and convex cases are obtained.     

Approximations of pseudo-Boolean functions; applications to game theory

This paper studies the approximation of pseudo-Boolean functions by linear functions and more generally by functions of (at most) a specified degree. Here a pseudo-Boolean function means a real valued function defined on {0,1} ⁿ , and its degree is that of the unique multilinear polynomial that expresses it; linear functions are those of degree at most one. The approximation consists in choosing among all linear functions the one which is closest to a given function, where distance is measured by the Euclidean metric onR ²ⁿ . A characterization of the best linear approximation is obtained in terms of the average value of the function and its first derivatives. This leads to an explicit formula for computing the approximation from the polynomial expression of the given function. These results are later generalized to handle approximations of higher degrees, and further results are obtained regarding the interaction of approximations of different degrees. For the linear case, a certain constrained version of the approximation problem is also studied. Special attention is given to some important properties of pseudo-Boolean functions and the extent to which they are preserved in the approximation. A separate section points out the relevance of linear approximations to game theory and shows that the well known Banzhaf power index and Shapley value are obtained as best linear approximations of the game (each in a suitably defined sense).Supported by the Air Force Office of Scientific Research (under grant number AFOSR 89-0512 and AFOSR 90-0008 to Rutgers University), as well as the National Science Foundation (under grant number DMS 89-06870).     

Some solutions (linear in the spatial variables) and generalized Chapman–Enskog functions basic to the linearized Boltzmann equations for a binary mixture of rigid spheres

A Legendre expansion of the (matrix) scattering kernel relevant to the (vector- valued) linearized Boltzmann equation for a binary mixture of rigid spheres is used to define twelve solutions that are linear in the spatial variables {x, y, z}. The twelve (asymptotic) solutions are expressed in terms of three vector-valued functions A ⁽¹⁾(c), A⁽²⁾(c), and B(c). These functions are generalizations of the Chapman–Enskog functions used to define asymptotic solutions and viscosity and heat conduction coefficients for the case of a single-species gas. To provide evidence that the three Chapman–Enskog vectors exist as solutions of the defining linear integral equations, numerical results developed in terms of expansions based on Hermite cubic splines and a collocation scheme are reported for two binary mixtures (Ne-Ar and He-Xe) with various molar concentrations.