An efficient algorithm for calculation of the Luenberger canonical form

A new algorithm is presented to obtain the Luenberger canonical form for multivariable systems. A distinct feature of the method is that the canonical form is obtained directly and, if necessary, the similarity transformation can be computed. There is a substantial reduction in the amount of computation compared to Luenberger´s method. The reduced computations along with Gaussian techniques lend greater inherent accuracy and the ability to refine the solution with additional computations. An example is presented to illustrate the technique.