In this paper, we propose a method to model nonlinear systems using polynomial nonlinear state space equations. Obtaining good initial estimates is a major problem in nonlinear modelling. It is solved here by identifying first the best linear approximation of the system under test. The proposed identification procedure is successfully applied to measurements of two physical systems. (C) 2010 Elsevier Ltd. All rights reserved.