We have reproduced some sunspot wave signatures obtained from spectropolarimetric observations through 3D MHD numericalsimulations. The results of the simulations arecompared with the oscillations observed simultaneously at different heights from the SiI lambda10827Å line, HeI lambda10830Å line, the CaII H core and the FeI blends at the wings of the CaII H line. The simulations show a remarkable agreement with the observations, and we have used them to quantify the energy contribution of the magneto-acoustic waves to the chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is 5-10 times lower than the amount needed to balance the chromospheric radiative losses.