The main goal of this work is to realize a PEMFC model that can be used efficiently for the global modelling of the fuel cell system. The modelling method proposed in the paper is an approach from an, empirical point of view that allows a PEMFC model of "black-box" class to be developed. Moving least squares (MLS) have therefore been employed to approximate the cell voltage characteristics V, using an experimental dataset measured in determinate conditions. The MLS approach appears to present a good balance of response surface accuracy, smoothness, robustness, and ease of use. This kind of numerical model offers good perspectives for the systems identification, the simulation of the systems, the design and the optimization of process control, etc. The results prove that the method is suitable for predicting and describing the fuel cell behaviour in all the points of the approximation domain. The proposed model can be included in a numerical application to optimize the operation of an existing fuel cell system. (C) 2007 Elsevier B.V. All rights reserved.