There exist many methods for facial recognition, such as eigenface, templates, artificial neural networks, etc., based on the given facial sample data (patterns). When an input facial image (target) involve simple geometrical transformations and illumination, the performance of these methods are not very satisfactory. In this thesis, following Li et al., we propose a new face recognition system, which can eliminate translation, rotation, scaling, and prospective transformations of facial images automatically, and can also eliminate illumination. According to facial features, we use this method to find the best transformation and the closet illumination, and then to eliminate them for identification by the best matching between a target and the patterns. Finally, we use the least squares method to recognize the target. This method is validated by numerical examples.