The effect of preservation of the residual polarization of backscattered light in the case of multiply scattered disordered media illumination by a linearly polarized plane wave is examined using the path-integral approach and Monte Carlo simulation. Disordered ensembles of non-interacting dielectric particles are considered as the model of scattering media. The influence of the anisotropy parameter of the scattering system on the degree of residual polarization is analysed. Experimental results obtained for various scattering systems at different wavelengths of illuminating light are in satisfactory agreement with the results of theoretical analysis and Monte Carlo simulation. The dependence of statistical properties of the polarization states of backscattered field partial components, such as probability distributions of ellipticity, on the anisotropy parameter is studied.