  The Dirac-{\delta} distribution may be realized through sequences of
convlutions, the latter being also regarded as approximation to the identity.
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 in signal and
image analysis, as well as in system identification.
