  In this paper, we make a multi-epoch analysis of NGC 1052, one of the
prototypical LLAGN, using XMM-Newton, Suzaku and NuSTAR observations, taken
from 2001 to 2017. This is the first time that results from NuSTAR observations
are reported for NGC 1052. On the technical aspects, we found a
wavelength-dependent calibration issue between simultaneous XMM-Newton and
NuSTAR spectra. It is described by a change on the photon index of $\rm{
\Gamma_{NuSTAR}- \Gamma_{XMM-Newton}=0.17\pm0.04}$. We use ancillary Chandra
data to decontaminate the nuclear spectrum from circumnuclear contributors. We
find that two baseline models can fit the broad (0.5-50 keV) X-ray spectrum of
the source. One consists of a power-law like continuum which is absorbed by a
uniform absorber, reflection from neutral material, and a separate power-law
component in the soft band. The second model presents a clumpy absorber. The
reflection component is still present, but not the soft band power-law.
Instead, absorption by a warm absorber is necessary to fit the spectra. This is
the first time that a reflection component is established in this object,
thanks to high energy data from NuSTAR. This component is constant in flux and
shape, supporting the idea that is produced away from the central source
(probably the torus). We find flux, spectral slope and absorption variations on
timescales of months to years. We also find that a patchy-absober can explain
the behaviour of this source better as it is $\sim$ 200 times more likely than
the uniform absober while it yields to smaller intrinsic variations.
