In this paper, a computer-aided method for reducing the order of a digital control system is developed. This is done by minimizing a weighted mean square error between the frequency responses of the reduced and the original systems. The method is demonstrated by a numerical example and the results are compared with those of the continued fractions method developed by Shin and Wu and the optimization method of model reduction developed by Luus.