Abstract
The problem of computing oscillatory integrals with general oscillators is considered. We employ a Filon-type method, where the interpolation basis functions are chosen in such a way that the moments are in terms of elementary functions and the oscillator only. This allows us to evaluate the moments rapidly and easily without needing to engage hypergeometric functions. The proposed basis functions form a Chebyshev set for any oscillator function even if it has some stationary points in the integration interval. This property enables us to employ the Filon-type method without needing any information about the stationary points if any. Interpolation by the proposed basis functions at the Fekete points (which are known as nearly optimal interpolation points), when combined with the idea of splines, leads to a reliable convergent method for computing the oscillatory integrals. Our numerical experiments show that the proposed method is more efficient than the earlier ones with the same advantages.
Similar content being viewed by others
References
Bos, L., De Marchi, S., Sommariva, A., Vianello, M.: Computing multivariate Fekete and Leja points by numerical linear algebra. SIAM J. Numer. Anal. 48(5), 1984–1999 (2010)
Bos, L., Levenberg, N.: On the calculation of approximate Fekete points: the univariate case. Electron. Trans. Numer. Anal. 30, 377–397 (2008)
Deaño, A., Huybrechs, D.: Complex Gaussian quadrature of oscillatory integrals. Numer. Math. 112, 197–219 (2009)
Domínguez, V., Graham, I., Kim, T.: Filon-Clenshaw-Curtis rules for highly oscillatory integrals with algebraic singularities and stationary points. SIAM J. Numer. Anal. 51(3), 1542–1566 (2013)
Domínguez, V., Graham, I., Smyshlyaev, V.: Stability and error estimates for Filon-Clenshaw-Curtis rules for highly oscillatory integrals. IMA J. Numer. Anal. 31, 1253–1280 (2011)
Fejér, L.: Bestimmung derjenigen Abszissen eines Intervalles, für welche die Quadratsumme der Grundfunktionen der Lagrangeschen Interpolation im Intervalle [−1, 1] ein möglichst kleines Maximum besitzt. Ann. Scuola Norm. Sup. Pisa Sci. Fis. Mt. Ser. II 1, 263–276 (1932)
Huybrechs, D., Olver, S.: Superinterpolation in highly oscillatory quadrature. Found. Comput. Math. 12, 203–228 (2012)
Huybrechs, D., Vandewalle, S.: On the evaluation of highly oscillatory integrals by analytic continuation. SIAM J. Numer. Anal. 44, 1026–1048 (2006)
Iserles, A., Nørsett, S.: Efficient quadrature of highly oscillatory integrals using derivatives. Proc. R. Soc. Lond. A 461, 1383–1399 (2005)
Majidian, H., Babolian, E., Malekifard, E.: A hybrid algorithm for approximating Fekete points of 1d interpolation problems with arbitrary basis functions. Submitted
Olver, S.: Moment-free numerical approximation of highly oscillatory integrals with stationary points. Eur. J. Appl. Math. 18, 435–447 (2007)
Phillips, G.: Estimate of the maximum error in best polynomial approximations. Comput. J. 11(1), 110–111 (1968)
Shampine, L.: Integrating oscillatory functions in Matlab, II. Electron. Trans. Numer. Anal. 39, 403–413 (2012)
Sommariva, A., Vianello, M.: Computing approximate Fekete points by QR factorizations of vandermonde matrices. Comput. Math. Appl. 57(8), 1324–1336 (2009)
Sommariva, A., Vianello, M.: Approximate Fekete points for weighted polynomial interpolation. Electron. Trans. Numer. Anal. 37, 1–22 (2010)
Xiang, S.: Efficient Filon-type methods for \({\int }^{b}_{a}f(x)e^{\mathrm {i}\omega g(x)}\,\mathrm {d} x\). Numer. Math. 105, 633–658 (2007)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Majidian, H. Automatic computing of oscillatory integrals. Numer Algor 77, 867–884 (2018). https://doi.org/10.1007/s11075-017-0343-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-017-0343-8