For a given graph G of order n with m edges, and a real symmetric matrix associated to the graph, M(G)∈Rn×n, the interlacing graph reduction problem is to find a graph Gr of order r<n such that the eigenvalues of M(Gr) interlace the eigenvalues of M(G). Graph contractions over partitions of the vertices are widely used as a combinatorial graph reduction tool. In this study, we derive a graph reduction interlacing theorem based on subspace mappings and the minmax theory. We then define a class of edge-matching graph contractions and show how two types of edge-matching contractions provide Laplacian and normalized Laplacian interlacing. An O(mn) algorithm is provided for finding a normalized Laplacian interlacing contraction and an O(n2+nm) algorithm is provided for finding a Laplacian interlacing contraction.