System approximation by matching the impulse response energies