Model reduction by matching Markov parameters, time moments, and impulse-response energies