A multi-stage least squares method for identifying Hammerstein model nonlinear systems