Numerical integrating of highly oscillating functions: effective stable algorithms in case of linear phase

04/08/2021
by   Leonid A. Sevastianov, et al.
0

A practical and simple stable method for calculating Fourier integrals is proposed, effective both at low and at high frequencies. An approach based on the fruitful idea of Levin, to use of the collocation method to approximate the slowly oscillating part of the antiderivative of the desired integral, allows reducing the calculation of the integral of a highly oscillating function (with a linear phase) to solving a system of linear algebraic equations with a three-diagonal triangular or five-diagonal band Hermitian matrix. The choice of Gauss-Lobatto grid nodes as collocation points makes it possible to use the properties of discrete "orthogonality" of Chebyshev differentiation matrices in physical and spectral spaces. This is realized by increasing the efficiency of the numerical algorithm for solving the problem. The system pre-conditioning procedure leads to significantly less cumbersome and more economical calculation formulas. To avoid possible numerical instability of the algorithm, we proceed to the solution of a normal system of linear algebraic equations.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
11/26/2017

Solution of a System of Linear Equations in an Integral Ring

A modified Gauss's algorithm for solving a system of linear equations in...
research
09/20/2023

Preconditioning for time-harmonic Maxwell's equations using the Laguerre transform

A method of numerically solving the Maxwell equations is considered for ...
research
04/01/2021

Mixed precision recursive block diagonalization for bivariate functions of matrices

Various numerical linear algebra problems can be formulated as evaluatin...
research
07/14/2022

On the numerical solution of Volterra integral equations on equispaced nodes

In the present paper, a Nystrom-type method for second kind Volterra int...
research
01/13/2021

A three-dimensional Laguerre one-way wave equation solver

A finite difference algorithm based on the integral Laguerre transform i...
research
11/01/2020

Highly accurate decoupled doubling algorithm for large-scale M-matrix algebraic Riccati equations

We consider the numerical solution of large-scale M-matrix algebraic Ric...
research
02/25/2023

Automated tuning for the parameters of linear solvers

Robust iterative methods for solving systems of linear algebraic equatio...

Please sign up or login with your details

Forgot password? Click here to reset