In this paper, we show that, for the discrete multivariable systems, the solution of the Lyapunov equation in the controllable canonical form transforms a companion matrix to another companion matrix. We also show that the first (last)

columns of the inverse of this solution matrix, where the block size is

, give the left (right) matrix fraction decomposition of a strictly positive real function.