A new method for the approximation of bilinear systems is proposed. The reduction scheme applies to both stable and unstable bilinear systems. The technique uses generalized input normal representations to retain the dominant part of the original system. The algorithm is evaluated on a synchronous induction generator and is shown to lead to acceptable reduced approximations of the original system. A frequency weighting is also introduced in the reduction scheme to further improve the approximation.