Subspace-based rational interpolation of analytic functions from phase data