EN
Results of Monte Carlo simulations of a 2D system of charged donors are presented. They enable to study the effects related to a spatial correlation of donor charges located on a random donor matrix. A qualitative difference between DX^{+} and DX¯ models is observed. In the first case, strong temperature dependence of the correlations and a "freezing-like" behaviour is found. The origin of the freezing is traced to the random distribution of donor sites by comparing the system with a liquid-like model where the charges may assume arbitrary positions within a plane. In the second case only nearest-neighbour correlations are observed. The simulations have direct application in analysis of the behaviour of the DX centers in the GaAs planarly doped with Si.