This paper proposes a frequency transfer function synthesis for approximating a high-order model with resonance to a low-order model in the frequency domain. The presented model approximation method is based on minimizing the error function weighted by the numerator polynomial of approximated models, which is used of the RLS(Recursive Least Square) technique to estimate the coefficient vector of approximated models. The proposed method provides better fitting in a low frequency and peak resonance. And an example is given to illustrate feasibilities of the suggested schemes.