Solution of Wiener-Hopf and Fredholm integral equations by fast Hilbert and Fourier transforms
Abstract
We present numerical methods based on the fast Fourier transform (FFT) to solve convolution integral equations on a semi-infinite interval (Wiener-Hopf equation) or on a finite interval (Fredholm equation). We extend and improve a FFT-based method for the Wiener-Hopf equation due to Henery, expressing it in terms of the Hilbert transform, and computing the latter in a more sophisticated way with sinc functions. We then generalise our method to the Fredholm equation reformulating it as two coupled Wiener-Hopf equations and solving them iteratively. We provide numerical tests and open-source code.
- Publication:
-
arXiv e-prints
- Pub Date:
- June 2021
- DOI:
- 10.48550/arXiv.2106.05326
- arXiv:
- arXiv:2106.05326
- Bibcode:
- 2021arXiv210605326G
- Keywords:
-
- Mathematics - Numerical Analysis;
- 45B05;
- 65R20;
- 65T99 (Primary) 45E10;
- 44A15 (Secondary);
- G.1.9
- E-Print:
- 30 pages, 19 figures