Abstract:The study of trigonometric series has started at the beginning of the nineteenth century. Joseph Fourier made the important observation that almost every function of a closed interval can be decomposed into the sum of sine and cosine functions. This technique to develop a function into a trigonometric series was published for the first time in 1822 by Joseph Fourier. The resulting series is nowadays called Fourier series. Since Fourier’s time, many different approaches to understand the concept of Fourier series have been discovered, each of which emphasizes different aspects of the topic. Some of the more powerful and elegant approaches are based on mathematical ideas and tools that were not available at the time Fourier completed his original work. Although the original motivation was to solve the heat equation for a metal plate, it later became obvious that the same technique could be applied to a wide variety of mathematical and physical problems and has many applications in electrical engineering, vibration analysis, acoustics, optics, signal treatment, image processing, etc. Despite the importance of Fourier series, the method used until now to compute them via computer algebra systems (CAS) is essentially based on the same principle as in Fourier’s time, i.e. by the evaluation of certain integrals. Unfortunately this technique is not completely successful for many functions. Although numeric values of the Fourier coefficients might be available, symbolic values are often not accessible. Modern CAS like Maple or Mathematica can compute such integrals in many cases for a given n∈Z\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${n \in \mathbb{Z}}$$\end{document}. However if one is interested in the Fourier coefficients for all n∈Z\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${n \in \mathbb{Z}}$$\end{document}, then n is considered as a given symbolic variable and such integrals can be computed only in few cases. In this paper we introduce an algorithmic approach to compute those Fourier coefficients, involving differential equations of a particular form, and recurrence equations. This approach extrapolates the computation of the Fourier series for functions for which the computation of Fourier coefficients via the definition is out of reach for current CAS.A holonomic recurrence equation for an, i.e. a recurrence equation which is linear, homogeneous and has polynomial coefficients, can be written in operator notation as L(an) = 0. The operator L can be interpreted as a non-commutative polynomial via the commutator rule Nn − nN = N, N denoting the shift operator Nan = an+1. In the last section we show how our algorithm can be used to factorize such recurrence operators in certain cases.

Numerical Methods for Fast Nonlinear Fourier Transformation, Part I: Exponential Runge-Kutta and Linear Multistep Methods

Fast Numerical Nonlinear Fourier Transforms

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

Rapid variable-step computation of dynamic convolutions and Volterra-type integro-differential equations: RK45 Fehlberg, RK4

Algorithmic Approach for Formal Fourier Series

Nonlinear Stability and Convergence of ERKN Integrators for Solving Nonlinear Multi-Frequency Highly Oscillatory Second-Order ODEs with Applications to Semi-Linear Wave Equations

Symmetric integrators with improved uniform error bounds and long-time conservations for the nonlinear Klein-Gordon equation in the nonrelativistic limit regime

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

Exponentially accurate spectral and spectral element methods for fractional ODEs

Fast and Efficient Numerical Methods for an Extended Black-Scholes Model

Fast Numerical Approximation of Linear, Second-Order Hyperbolic Problems Using Model Order Reduction and the Laplace Transform

A Fast Discrete Spectral Method for Stochastic Partial Differential Equations

High-order exponential integrators for the Riesz space-fractional telegraph equation

Accelerating the convergence of Newton's method for nonlinear elliptic PDEs using Fourier neural operators

Fast Fourier Solvers for the Tensor Product High-Order FEM for a Poisson Type Equation

Simplified explicit exponential Runge-Kutta methods without order reduction

Low-synchronization Arnoldi Methods for the Matrix Exponential with Application to Exponential Integrators

Numerical Conservation of Energy, Momentum and Actions for Extended RKN Methods when Applied to Nonlinear Wave Equations Via Spatial Spectral Semi-Discretizations

Computation of highly oscillatory integrals using a Fourier extension approximation

Fast IMEX Time Integration of Nonlinear Stiff Fractional Differential Equations