It is possible to obtain either exact realizations or useful approximations of linear systems or matrix-vector products arising in many different applications, by synthesizing them in the form of repeated or multi-channel filtering operations in fractional Fourier domains, resulting in much more efficient implementations with acceptable decreases in accuracy. By varying the number and configuration of filter blocks, which may take the form of arbitrary flow graphs, it is possible to trade off between accuracy and efficiency in the desired manner. The proposed scheme constitutes a systematic way of exploiting the information inherent in the regularity or structure of a given linear system or matrix, even when that structure is not readily apparent.