Maybe convolution on a discrete basis can be represented as a matrix multiplication with a circulant matrix?