The invention provides an artificial wave response spectrum fitting calculation method based on multi-phase spectrum iteration, and aims at solving the problems existing in artificial seismic waves ofa traditional frequency domain method fitting envelope design response spectrum, when Fourier amplitude spectrum iteration is carried out, R * 10 groups of random phase spectrums are generated and R* 10 artificial waves are generated in the R-th iteration, and the R * 10 groups of random phase spectrums are calculated. Iteration is carried out by using an error between the average value of the reaction spectrums corresponding to the artificial waves and the designed reaction spectrum, so that the condition that the final fitting result of the random phase spectrum cannot be converged is avoided. In the process of designing the reaction spectrum through frequency domain method fitting, a multi-phase spectrum matrix is mainly used, a single and invariant phase spectrum is not used any more, a complex Fourier phase spectrum is used for each iteration, and the influence of the random phase spectrum of the artificial wave on the convergence precision of a final fitting result is also reduced.