Motivated by image processing and numerical wavelet methods for Partial differential equations, we study the theoretical interactions between wavelets and the diffusion equations. Important properties of wavelets, such as the translation and scaling invariance, the p-vanishing-moment condition, and the atomic decomposition, are integrated into the diffusion process and lead to many interesting results.