In this paper, the reduction method uses the concepts of stability-equation and important poles to find the denominator of the reduced model. Then the numerator of the reduced model is found by complex curve fitting. This method tends to simultaneously guarantee a stable reduced model from a stable system and obtain a satisfactory result, since it considers the distribution of important poles. Examples are presented to illustrate this advantage.