On the calculation of highly oscillatory integrals with an algebraic singularity

In this paper we explore the evaluation of highly oscillatory integrals containing an algebraic singularity based on the change of variable t = x^−r, where r is a positive real number, and the analytical continuation. Robust numerical results demonstrate the accuracy and effectiveness of the proposed approach.     

Uniform approximations to Cauchy principal value integrals of oscillatory functions

This paper presents an interpolatory type integration rule for the numerical evaluation of Cauchy principal value integrals of oscillatory integrands , where -1<τ<1, for a given smooth function f(x). The proposed method is constructed by interpolating f(x) at practical Chebyshev points and subtracting out the singularity. A numerically stable procedure is obtained and the corresponding algorithm can be implemented by fast Fourier transform. The validity of the method has been demonstrated by several numerical experiments.     

一种非规则积分区域的二重高振荡函数积分方法

在传统L ev in方法与新F ilon型方法的基础上,本文提出了一种求解非规则区域下的二重高振荡函数数值积分方法,通过利用L ev in匹配法将二重积分化为一重积分,并避免了对复杂的m om en ts的求解,能提高计算的效率,且有很高的求积精度.     

Note on the homotopy perturbation method for multivariate vector-value oscillatory integrals

This paper considers a homotopy perturbation method for approximating multivariate vector-value highly oscillatory integrals. The asymptotic formulae of the integrals and the asymptotic order of the asymptotic method are presented. Numerical examples show the efficiency of the approximation method.     

Numerical analysis of a fast integration method for highly oscillatory functions

The integration of systems containing Bessel functions is a central point in many practical problems in physics, chemistry and engineering. This paper presents a new numerical analysis for the collocation method presented by Levin for and gives more accurate error analysis about the integration of systems containing Bessel functions. The effectiveness and accuracy of the quadrature is tested for Bessel functions with large arguments. AMS subject classification (2000)  65D32, 65D30     

Decay estimates for oscillatory integrals with polynomial phase in terms of p and p or in terms of p and p

We obtain estimates for certain oscillatory integrals with polynomial (degree n) phase, p(t). These estimates are stated in terms of differences between the roots, real or complex, of p⁽ⁿ⁻³⁾(t) and p⁽ⁿ⁻²⁾(t) or between p⁽ⁿ⁻²⁾(t) and p⁽ⁿ⁻¹⁾(t). The sharpness of these results is also explored. This result is a partial generalization of the results found in [J. Math. Anal. Appl. 280 (2003) 424].     

Asymptotic approximation of a highly oscillatory integral with application to the canonical catastrophe integrals

We consider the highly oscillatory integral $F(w):={\int }_{-\infty }^{\infty }{e}^{iw(t^{K+2}+e^{i\theta }{t}^p)}g(t)dt$ for large positive values of w, , K and p positive integers with $1\le p\le K$, and $g(t)$ an entire function. The standard saddle point method is complicated and we use here a simplified version of this method introduced by López et al. We derive an asymptotic approximation of this integral when $w\to +\infty$ for general values of K and p in terms of elementary functions, and determine the Stokes lines. For $p\ne 1$, the asymptotic behavior of this integral may be classified in four different regions according to the even/odd character of the couple of parameters K and p; the special case $p=1$ requires a separate analysis. As an important application, we consider the family of canonical catastrophe integrals ${\mathrm{\Psi }}_K(x_1,x_2,\text{…},x_K)$ for large values of one of its variables, say $x_p$, and bounded values of the remaining ones. This family of integrals may be written in the form $F(w)$ for appropriate values of the parameters w, θ and the function $g(t)$. Then, we derive an asymptotic approximation of the family of canonical catastrophe integrals for large $|x_p|$. The approximations are accompanied by several numerical experiments. The asymptotic formulas presented here fill up a gap in the NIST Handbook of Mathematical Functions by Olver et al.     

Numerical approximation of vector-valued highly oscillatory integrals

We present a method for the efficient approximation of integrals with highly oscillatory vector-valued kernels, such as integrals involving Airy functions or Bessel functions. We construct a vector-valued version of the asymptotic expansion, which allows us to determine the asymptotic order of a Levin-type method. Levin-type methods are constructed using collocation, and choosing a basis based on the asymptotic expansion results in an approximation with significantly higher asymptotic order. AMS subject classification (2000)  65D30     

estimates for a class of oscillatory integrals

We obtain a simpler proof of Theorem 3.1 of The complete mapping properties of some oscillatory integrals in several dimensions, by G. Sampson and P. Szeptycki (Canad. Math. J. 53 (5) (2001), 1031-1056).

     

A quadrature method for a class of strongly oscillatory infinite integrals

A constructive process is presented which leads to a class of quadrature formulas with any preassigned compound precision for the numerical integration of rapidly oscillating functions on (0, ) of the forme ^–x F(x, x), where is a large parameter andF(·,y) is periodic of period unity iny.1980 Mathematics Subject Classification, Primary 41A55, Scondary 41A60     

Moment-free numerical integration of highly oscillatory functions

^** Email: s.olver{at}damtp.cam.ac.uk The aim of this paper is to derive new methods for numericallyapproximating the integral of a highly oscillatory function.We begin with a review of the asymptotic and Filon-type methodsdeveloped by Iserles and Nørsett. Using a method developedby Levin as a point of departure, we construct a new methodthat utilizes the same information as a Filon-type method, andobtains the same asymptotic order, while not requiring the computationof moments. We also show that a special case of this methodhas the property that the asymptotic order increases with theaddition of sample points within the interval of integration,unlike all the preceding methods whose orders depend only onthe endpoints.     

一些高振荡积分、高振荡积分方程的高性能计算

电磁、声波散射问题的研究涉及一类数学物理问题, 此类问题具有深刻的理论价值和重要的应用背景, 亟待解决. 高振荡微分、积分方程是刻画这些问题的重要的数学模型, 其数值计算存在许多挑战性研究课题. 本文从积分方程解法角度出发, 综述了求解这类高振荡问题的一些最新进展, 特别是针对广义Fourier 变换、Bessel 变换的高效算法、高振荡核Volterra 积分方程的数值解法作了详细介绍. 这些数值方法共有特点是振荡频率越高算法精度愈高, 且可望为电磁计算的研究提供一些新的高效算法.     

Efficient quadrature for highly oscillatory integrals involving critical points

This paper based on the Levin collocation method and Levin-type method together with composite two-point Gauss–Legendre quadrature presents efficient quadrature for integral transformations of highly oscillatory functions with critical points. The effectiveness and accuracy of the quadrature are tested.     

Efficient evaluation of highly oscillatory acoustic scattering surface integrals

We consider the approximation of some highly oscillatory weakly singular surface integrals, arising from boundary integral methods with smooth global basis functions for solving problems of high frequency acoustic scattering by three-dimensional convex obstacles, described globally in spherical coordinates. As the frequency of the incident wave increases, the performance of standard quadrature schemes deteriorates. Naive application of asymptotic schemes also fails due to the weak singularity. We propose here a new scheme based on a combination of an asymptotic approach and exact treatment of singularities in an appropriate coordinate system. For the case of a spherical scatterer we demonstrate via error analysis and numerical results that, provided the observation point is sufficiently far from the shadow boundary, a high level of accuracy can be achieved with a minimal computational cost.     

An automatic integration procedure for infinite range integrals involving oscillatory kernels

Let the real functionsK(x) andL(x) be such thatM(x)=K(x)+iL(x)=e^ix g(x), whereg(x) is infinitely differentiable for all largex and is non-oscillatory at infinity. We develop an efficient automatic quadrature procedure for numerically computing the integrals _a K(t)f(t) and _a L(t)f(t)dt, where the functionf(t) is smooth and nonoscillatory at infinity. One such example for which we also provide numerical results is that for whichK(x)=J (x) andL(x)=Y (x), whereJ (x) andY (x) are the Bessel functions of order . The procedure involves the use of an automatic scheme for Fourier integrals and the modified W-transformation which is used for computing oscillatory infinite integrals.     

Efficient numerical methods for Cauchy principal value integrals with highly oscillatory integrands

Numerical Algorithms - This paper focus on the numerical evaluation of the Cauchy principal value integrals with oscillatory integrands where α, β &gt; ??1,??1...