ABSTRACTA new technique for frequency limited model order reduction of discrete time second-order systems is presented. Discrete time frequency limited Gramians (DFLGs) and corresponding discrete algebraic Lyapunov equations are developed. An efficient technique for the computation of DFLGs and their Cholesky factors is presented. Computed DFLGs are partitioned to obtain position and velocity Gramians. These Gramians are balanced with different combinations to obtain various balanced transformations that yield Hankel singular values (HSVs) for order reduction. Frequency limited discrete time balanced truncation framework is proposed and truncation based on magnitudes of HSVs is applied to obtain the reduced order model. Moreover, stability conditions for reduced order models are stated. Results of the proposed technique are compared with infinite Gramians balancing scheme in order to certify the usefulness of the presented technique for frequency limited applications.