This paper suggests the use of wavelet multiresolution analysis (WMRA) and neural network for generation of artificial earthquake accelerograms from target spectrum. This procedure uses the learning capabilities of radial basis function (RBF) neural network to expand the knowledge of the inverse mapping from response spectrum to earthquake accelerogram. In the first step, WMRA is used to decompose earthquake accelerograms to several levels that each level covers a special range of frequencies, and then for every level a RBF neural network is trained to learn to relate the response spectrum to wavelet coefficients. Finally the generated accelerogram using inverse discrete wavelet transform is obtained. An example is presented to demonstrate the effectiveness of the method.