The metal-metal contact nonlinearity can degrade the performance of communication system, which has become a hot issue in the electrical contact field. The key problem is to establish the relationship between the parameters of rough surface and the real contact area under the preload. The surface roughness parameters are dependent on the metal material and processing technology etc. The real contact area is generally calculated by the finite element method (FEM). However, there exist some difficulties for FEM such as exact model of rough surface and tedious calculation in the engineering. Aimed at this problem, we propose a new method to predict the real contact area of fractal rough surface, and establish the least-square support vector machines (SVM) regression model and optimize the hyper-parameters by the combination of the coupled simulated annealing and grid searching algorithms. The model can predict the real contact area of metal-metal contact for different working conditions, and address the problem to establish the relationship between the parameters of rough surfaces and the real contact area under the preload. The work lays good foundation for analysis of metal-metal contact nonlinearity.