In this article we shall use two special classes of reproducing kernel Hilbert spaces (which originate in the work of de Branges [dB) and de Branges-Rovnyak [dBRl), respectively) to solve matrix versions of a number of classical interpolation problems. Enroute we shall reinterpret de Branges' characterization of the first of these spaces, when it is finite dimensional, in terms of matrix equations of the Liapunov and Stein type and shall subsequently...