A reduced–order model for a multi–input single–output or a single–input multi–output continuous system is proposed. The model is especially suitable for the case where accuracy in the high-frequency region is important. It also satisfies basic conditions for reduced-order models, such as uniqueness, stability, and a simple calculation requirement. Chained aggregation with the minimization of truncation errors introduces an internal system representation characterized by a system matrix of Hessenberg form and a diagonalized covariance matrix.