Computing dominant poles of power system multivariable transfer functions

This paper describes a numerical linear algebraic algorithm to compute the dominant poles of multi-input-multi-output (MIMO) high-order transfer functions. The results presented are related to the study of electromechanical oscillations in large electrical power systems, but the algorithm is completely general. The computed dominant poles may then be used to build modal equivalents for MIMO transfer functions of large linear systems, among other applications.