This paper focuses on the parameter identification of a class of nonlinear mechanical dynamic systems subject to feedback quantization errors.} A second-order integral method for the parameter identification of a nonlinear two-degrees-of-freedom robotic manipulator is studied. To obtain accurate parameter estimates, an optimal excitation trajectory subject to state constraints is proposed. The designated excitation trajectory is able to minimize estimation uncertainties, providing the controlled system trajectory reaches the desired optimal excitation trajectory. Moreover, for most of the nonlinear robotic systems, the output measurements are acquired from the encoders subject to apparent quantization error. Directly differentiating the position to estimate the velocity and acceleration trajectories is unrealistic and unsuitable due to the amplified quantization errors will cause considerable estimation errors. Instead, the finite Fourier series is introduced to perform optimal estimation of the position, velocity, and acceleration histories to capture the dominant frequencies of the system response. A weighted least-squares solution is further proposed to improve parameter identification precision. The comparative studies reveal that the integral method has superior parameter estimation performance than the direct difference model and the filtered regression model in the existing literature. Experimental simulations illustrate the effectiveness of the presented method.