We consider a class of discrete differential operators acting on multidimensional Haar wavelet basis with the aim of finding wavelet approximate solutions of partial differential problems. Although these operators depend on the interpolating method used for the Haar wavelets regularization and the scale dimension space, they can be easily used to define the space of approximate wavelet solutions.