“…are the transfer functions of the original and reduced-order models, i = √ −1, and || · || 2 denotes the spectral matrix norm, see [43]. Moreover, the balanced truncation method preserves stability in (39). It should, however, be noticed that the preservation of other properties, in particular, passivity cannot, in general, be guaranteed for this method.…”