In this paper, we propose the discrete model reduction method of bounded real transfer functions. From the discrete bounded real lemma, we obtain the two riccati equations and define the disrete bounded real balancing using solutions of these two riccati equations. And we get the reduced order discrete model from the GSPA of full order model. Especially, when free parameter of GSPA is .+-.1, we show that the reduced order discrete model retains minimality, stability, and bounded real and BR-balancing properties. And we derive the .inf.-norm error bound between full order model and reduced order model. Finally to illustrate the validity of proposed method, we give an example.