Fourier expansions employing polyharmonic–Neumann eigenfunctions have demonstrated improved convergence over those using the classical trigonometric system, due to the rapid decay of their Fourier coefficients. Building on this insight, we investigate interpolations on a finite interval that are exact for polyharmonic–Neumann eigenfunctions and exhibit similar benefits. Furthermore, we enhance the convergence of these interpolations by incorporating the concept of quasi-periodicity, wherein the basis functions are periodic over a slightly extended interval. We demonstrate that those interpolations achieve significantly better convergence rates away from the endpoints of the approximation interval and offer increased accuracy over the entire interval. We establish these properties for a specific case of polyharmonic–Neumann eigenfunctions known as the modified Fourier system. For other basis functions, we provide supporting evidence through numerical experiments. While the latter methods display superior convergence rates, we demonstrate that interpolations using the modified Fourier basis offer distinct advantages. Firstly, they permit explicit representations via the inverses of certain Vandermonde matrices, whereas other interpolation methods require approximate computations of the eigenvalues and eigenfunctions involved. Secondly, these matrix inverses can be efficiently computed for numerical applications. Thirdly, the introduction of quasi-periodicity improves the convergence rates, making them comparable to those of other interpolation techniques.