DFT calculation of the generalized and Drazin inverse of a polynomial matrix

A new algorithm is presented for the determination of the (a) generalized inverse and (b) Drazin inverse, of a polynomial matrix. The proposed algorithms are based on the discrete Fourier transform and thus are computationally fast in contrast to other known algorithms. The above algorithms are implemented in the Mathematica programming language and illustrated via examples.