Abstract. The aim of the given paper is the development of an approach for parametric identification of Hammerstein systems with piecewise linear nonlinearities, i.e., when the saturation-like function with unknown slopes is followed by a linear part with unknown parameters. It is shown here that by a simple input data rearrangement and by a following data partition the problem of identification of a nonlinear Hammerstein system could be reduced to the linear parametric estimation problem. Afterwards, estimates of the unknown parameters of linear regression models are calculated by processing respective particles of input-output data. A technique based on ordinary least squares is proposed here for the estimation of parameters of linear and nonlinear parts of the Hammerstein system, including the unknown threshold of the piecewise nonlinearity, too. The results of numerical simulation and identification obtained by processing observations of input-output signals of a discrete-time Hamm...