The least square problem of the convolution result and real seismic data can be considered as the solution of a huge rarefactional matrix equation, which can be solved by singular value decomposition.