Abstract: In this paper, a fast recursive algorithm with two-grade QR decomposition for multivarible transfer function matrix identification is presented. By using the algorithm, repeatedly processing of the measured data of subsystems identification algorithms can be avoided, and without expanding of parameter vector dimension, a parameter estimation that minimizes the overall cost function can be obtained by using this algorithm. Combined with the fast HOUSEHOLDER transformation, this algorithm has better numerical stability and can enhance the estimation accuracy.