This paper presents a simple procedure for constructing reduced-order models that match some Markov parameters and time moments of an original system as well as the energies of its impulse response and of other combinations of the system modes. In this way, the stability of the reduced model of a stable system is ensured.