Abstarct: Ground response evaluation of is one of the most common and important issues in seismic geotechnics. Ground response analysis is used to predict ground surface movements and editing the response spectra to determine soil dynamic stresses and strains to assess geotechnical hazards (different types of ground instabilities) and calculate earthquake induced forces that can cause ground and retaining structures instability. In the classical solution of the wave equation, which uses the results of ground response analysis, the constant shear wave velocity is considered. With the help of auxiliary methods such as linear equivalence method it is tried to model the nonlinear essence of shear wave velocity, which this method is also unable to accurately model nonlinear soil behavior. Wavelet Analysis is relatively one of the new and exciting achievements of pure mathematics, which is baxsed on decades of research in harmonic analysis, and today has found important applications in many disciplines of science and engineering and new possibilities have been provided for understanding its mathematical aspects as well as enhancing its applications. In this method, the wave equation is solved with the help of the Haar wavelet by considering the shear wave velocity as a nonlinear function of depth. For this aim, two soil profiles of Gilroy and Treasure Islan regions with different input earthquakes have been studied. In order to investigate the effect of nonlinear soil behavior, the analysis is performed in two stages. The first step is to consider multilxayer soil profiles and the second one to consider single laxyer profiles. Finally the results of the two steps are compared with those of the Haar wavelet analysis and this comparison indicates that this method is able to solve nonlinear wave equation with any arbitrary function for shear wave velocity.