Structure-preserving model order reduction for K-power bilinear systems via Laguerre functions

This paper presents a series of structure-preserving model order reduction algorithms for K-power bilinear systems via Laguerre functions. The method first aims to rewritten the K-power bilinear system as a general bilinear system and calculate the approximate low-rank factors of the cross Gramian of the bilinear system by combining the idea of Laguerre functions expansion of the matrix exponential function. After that, the approximate balanced system of the K-power bilinear system is constructed by the corresponding projection transformation of each subsystem. In order to achieve the purpose of model order reduction, the states with smaller singular values are then truncated, so as to further obtain the reduced order model. For this approach, there is a disadvantage that unstable systems may be generated although the original one is stable. To alleviate the inadequacies of this approach, we have improved the model reduction procedure, which is based upon the dominant subspace projection method. In addition, we also carried out a correlation analysis on the stability of improved algorithms. Finally, numerical experiments are employed to substantiate the effectiveness of the presented algorithms.