A new method for obtaining generalized q-Markov Cover models for discrete-time SISO systems is proposed and is based on the inverse solution of the Lyapunov equation. The reduced order models not only match the Markov parameters and high-frequency power moments like the conventional q-Markov Cover techniques but also match time moments and low-frequency power moments of the original system.