Periodic trapezoidal Fourier aliasing (source code)

= Periodic trapezoidal Fourier aliasing

For a two-periodic Fourier series sampled at $2N$ equally spaced points, discrete averaging retains exactly the modes divisible by $2N$:
$$
I_N(h)-I(h)=\sum_{j\ne0}\widehat h_{2Nj}.
$$
If $|\widehat h_n|\leq Mc^{|n|}$ with $0<c<1$, then
$$
|I_N-I|\leq\frac{2Mc^{2N}}{1-c^{2N}},
$$
which is exponentially small.