This paper extends the polynomial function to double logarithmic function, constructing a class of multi parameters iterative function, and uses this function to calculate the fractal interpolated surface for given interpolation points, and establishes the iterative function mathematical model of multi parameters fractal interpolation. In order to verify the effectiveness and reliability of this proposed model algorithm, this paper uses MATLAB numerical simulation method to calculate, and programs the Newton iterative function of multi parameters fractal interpolation surface, finally gets calculation nephogram of multi parameters fractal interpolation curved surface through calculating. Finally, using iterative method reduces the surface grid size, increasing the smoothness of the surface, so the surface is closer to the actual surface. It provides a new computer method research of fractal interpolation function.