Order reduction of bilinear MIMO dynamical systems using new block Krylov subspaces

In this paper we study numerical methods for the model-order reduction of large-scale bilinear multi-input multi-output systems. A new projection method is proposed. The projection subspace is the union of some new block Krylov subspaces. We show that the reduced-order bilinear system constructed by the new method can match a desired number of moments of multivariable transfer functions corresponding to the kernels of Volterra series representation of the original system. Some numerical examples are presented to illustrate the effectiveness of the proposed method.