SummaryIntensive research in the field of mathematical modeling of pneumatic servo drives has shown that their mathematical models are nonlinear in which many important details cannot be included in the model. Owing the influence of the combination of heat coefficient, unknown discharge coefficient, and change of temperature, it was supposed that parameters of the pneumatic cylinder are random (stochastic parameters). On the other side, it has been well known that the nonlinear model can be approximated by a linear model with time‐varying parameters. Due to the aforementioned reasons, it can be assumed that the pneumatic cylinder model is a linear stochastic model with variable parameters. In practical conditions, in measurements, there are rare, inconsistent observations with the largest part of population of observations (outliers). Therefore, synthesis of robust algorithms is of primary interest. In this paper, the robust recursive algorithm for output error models with time‐varying parameters is proposed. The convergence property of the proposed robust algorithm is analyzed using the methodology of an associated ordinary differential equation system. Because ad hoc selection of model orders leads to overparameterization or parsimony problem, the robust Akaike's criterion is proposed to overcome these problems. By determining the least favorable probability density for a given class of probability distribution represents a base for design of the robust version of Akaike's criterion. The behavior of the proposed robust identification algorithm is considered through intensive simulations that demonstrate the superiority of the robust algorithm in relation to the linear algorithms (derived under an assumption that the stochastic disturbance has a Gaussian distribution). The good practical values of the proposed robust algorithm to identification of the pneumatic cylinder are illustrated by experimental results. Copyright © 2016 John Wiley & Sons, Ltd.