Abstract:A multi-color study of the variability of the magnetic cataclysmic variable BY Cam is presented. The observations were obtained at the Korean 1.8 m and Ukrainian 2.6 m, 1.2 m and 38 cm telescopes in [2003][2004][2005] 56 observational runs cover 189 hours. The variations of the mean brightness in different colors are correlated with a slope R/ V = 1 29(4), where the number in brackets denotes the error estimates in the last digits. For individual runs, this slope is much smaller ranging from 0.98(3) to 1.24(3), with a mean value of 1.11(1). Near the maximum, the slope becomes smaller for some nights, indicating more "blue" spectral energy distribution, whereas the night-to-night variability has an "infrared" character. For the simultaneous UBVRI photometry, the slopes increase with wavelength from U/ R = 0 23(1) to I/ R = 1 18(1). Such wavelength dependence is the opposite of that observed in non-magnetic cataclysmic variables, in agreement with the model of cyclotron emission. The principal component analysis shows two components of variablitity with different spectral energy distributions (with a third at the limit of detection), which possibly correspond to different regions of emission. The highest peak in the scalegram analysis corresponds to the 200 min spin variability, its quarter and to the 30 min and 8 min QPOs. The amplitudes of these components are dependent on wavelength and luminosity state. The light curves were fitted by a statistically optimal trigonometrical polynomial (up to 4th order) to take into account a 4-hump structure. The dependences of these parameters on the phase of the beat period and on mean brightness are discussed.