One can obtain either exact realizations or useful approximations of linear systems or matrix-vector products that arise in many different applications by implementing them in the form of multistage or multichannel fractional Fourier-domain filters, resulting in space-bandwidth-efficient systems with acceptable decreases in accuracy. Varying the number and the configuration of filters enables one to trade off between accuracy and efficiency in a flexible manner. The proposed scheme constitutes a systematic way of exploiting the regularity or structure of a given linear system or matrix, even when that structure is not readily apparent.