This paper is aimed at combining wavelet multiresolution analysis and nonstationary Kanai-Tajimi model for the simulation of earthquake accelerograms. The proposed approach decomposes earthquake accelerograms using wavelet multiresolution analysis for the simulation of earthquake accelerograms. This study is on the basis of some Iranian earthquake records, namely Naghan 1977, Tabas 1978, Manjil 1990 and Bam 2003. The obtained results indicate that the simulated records preserve the significant properties of the actual accelerograms. In order to investigate the efficiency of the model, the spectral response curves obtained from the simulated accelerograms have been compared with those from the actual records. The results revealed that there is a good agreement between the response spectra of simulated and actual records.