The hyperfine structure of energy levels of muonic molecules tdµ, tpµ and dpµ is calculated on the basis of stochastic variational method. The basis wave functions are taken in the Gaussian form. The matrix elements of the Hamiltonian are calculated analytically. Vacuum polarization, relativistic and nuclear structure corrections are taken into account to increase the accuracy. For numerical calculation, a computer code is written in the MATLAB system. Numerical values of energy levels of hyperfine structure in muonic molecules tdµ, tpµ and dpµ are obtained.