Two frequency-domain reduction (approximation) methods are introduced, for the purpose of applying them to high-degree transfer functions of single-input single-output linear time-invariant systems and obtain corresponding adequate reduced order models. These two reduction methods (i.e. Pade and Routh) are coupled with fitting of time moments and Markov parameters matching techniques to better evaluate the adequacy of the obtained reduced order model. The mentioned methods have been applied successfully to the 6th-order transfer function of a practical hydro power system.