A two-stage subspace method is proposed to identify MDOF nonlinear system. The method gives more accurate estimations and is more robust. The method gives a more suitable estimation of the order of system. Numerical and experimental examples are given to verify relevant analysis.