Abstract: This talk introduces that the matrix fraction description (MFD), given by a pair fA(z) B(z)g of matrix polynomials of z, for a linear time-invariant system may not be unique even if A(z) is monic, A(z) and B(z) have no common left factor, and the composite matrix consisting of the highest-order terms of A(z) and B(z) are full row rank. The orders of all possible MFDs of a given system are completely characterized. Testing criteria for determining whether a matrix pair is an MFD of the system are derived, which involve rank tests of certain Toeplitz matrices derived from either the impulse response or output correlation functions of the system. A decision procedure is devised that generates sequentially all MFDs for a given system. Identification algorithms are introduced that estimate all MFDs of a given system from its input output data or output data only. The results are then extended to cover ARMAX systems.