The Dirac-δ distribution may be realized through sequences of convolutions, the latter being also regarded as approximation to the identity (Stein-Weiss). The present study proposes the so called pre-orthogonal adaptive Fourier decomposition (POAFD) method to realize fast approximation to the identity. The type of sparse representation method has potential applications to signal and image analysis, as well as to system identification.