The identification of Hammerstein models relates to a computerized method for identifying Hammerstein models in which the linear dynamic part is modeled by a space-state model and the static nonlinear part is modeled using a radial basis function neural network (RBFNN), and in which a particle swarm optimization (PSO) algorithm is used to estimate the neural network parameters and a numerical algorithm for subspace state-space system identification (N4SID) is used for estimation of parameters of the linear part.