Based on the Weibull distribution function, a rock meso-heterogeneity model was introduced to numerically simulate the Brazilian splitting test with discontinuous deformation analysis (DDA) method. The heterogeneities of the elastic modulus, the Poisson's ratio and the strength parameters were considered in this model. The simulation results of the heterogeneous and the homogeneous Brazilian disks show that the heterogeneity model can derive the macroscopic non-linear mechanical response of rock with meso-scale linear constitutive relationship and leads to lower rock macroscopic strength than the homogeneeus modle. And after comparing the simulation result of considering the heterogeneities of the elastic modulus and the Poisson's ratio with that of considering the heterogeneities of the elastic modulus, the Poisson's ratio and the strength parameters, we found the macroscopic strength of the disk with taking the heterogeneity of the strength parameters into consideration was lower and the cracks generated were more scattered. This study is beneficial for the understanding of the rock heterogeneous and non-linear characteristics and their numerical analysis.