We show that, in spite of rather common opinion, correlations of observables on subsystems of a composite quantum system can be represented as correlations of classical Gaussian variables. We restrict our model to the finite dimensional case (which is important, e.g., in quantum information theory). Here quantum correlations are represented by integrals with Gaussian densities which can be directly calculated. In particular, the EPRBell correlations for polarizations of entangled photons can be represented as correlations of Gaussian random variables. The tricky point of our construction is the necessity to introduce the Gaussian noise ("vacuum fluctuations") to obtain positively defined covariance matrices for "prequantum random variables." Quantum mechanics can be considered as a method to proceed by subtracting the influence of this Gaussian background noise.