In this paper, we study Lagrange-Sylvester interpolation of rational matrix functions which are analytic at infinity and propose a new interpolation algorithm based on the recent subspace-based identification methods. As an application, we consider the problem of system identification with interpolation constraints.