An efficient algorithm for partial fraction expansion of the linear matrix pencil inverse

Abstract A new algorithm for computations of matrix partial fractions representing the inverse of linear matrix pencil is based on an appropriate expression in matrix form of the Pascal triangle. It concerns singular and nonsingular systems and starts with the inverse of regular matrix linear pencil M(s) = sA 0 - A where only A 0 is singular, or both A 0 and A are singular. Nonsingular systems are considered as a particular case of singular systems. The presented algorithm of the matrix partial fraction expansion is suitable to determine the matrix transfer function, and is computer oriented because all manipulations can be performed on matrices with constant entries only.