Firstly, based on the meshing principle and the transformation relationship between generating coordinate system of NC grinding, the mathematical model of tooth surface of face gears can be derived, followed by a division of the tooth surface grids of face gear in 5×9, and a calculation of the theoretical coordinate values of nodes. Then three coordinate measuring machines (CMM) are used to measure the coordinate values of the grinding tooth surface joints, and tooth surface deviations of grid nodes can be obtained by analyzing and decomposing the tooth surface error. With the identification equation of the tooth surface error having been established, the sequential quadratic programming (SQP) method has been proposed to optimize the machining parameters of the machine tool, thus completing the correction of the tooth surface error of the opposite gears. The result of tooth surface error correction shows that the error of tooth surface in the working area tends to be stable, and the corrected tooth surface error is obviously improved, but the transition area near the tooth surface needs further improvement.