We study subhalo populations surrounding massive dark matter haloes by using three AMIGA simulations each having different mass and spatial resolutions. Our analysis shows that the slope of the subhalo mass function has a value 0.9, which agrees with previous studies. The dependence of mass functions on redshift is the same for subhaloes and main haloes. In all simulations, combined subhalo masses are about 0.1-0.2 of main halo masses and this mass fraction increases slightly with redshift and the mass of the main halo. The distribution of mass fractions for subhaloes is close to Gaussian at z = 0 and differs slightly at earlier epochs. Spatial distribution of subhaloes as measured in units of virial radius Rvir of the main halo does not depend on redshift and follows r1/3 rule. Spatial distribution of all haloes surrounding main haloes continues up to 3 times Rvir with equal slope but lower amplitude. Beyond 16 times Rvir, the average distribution of haloes becomes uniform.