In this paper, a new model reduction method is proposed to obtain a reduced order model in the frequency domain. The method is developed based on the second-order plus dead time modeling technique. The initial value of the reduced model parameters can be obtained using this method coinciding four point(0, -$pi$/2, -$pi$, -3$pi$/2) on the Nyquist curve. The optimal parameters of the reduced model is obtained through calculation procedure with three steps. It is shown that Nyquist curves and unit step responses of the reduced models of numerical examples closely agree with those of original models.