Our work is devoted to the creation of a generalized dynamic model that describes libration processes in the motion of the natural satellite of the planet, which has been tested on the simulation of the moons rotational parameters(MRPs) of the Earth’s Moon and allows adapting description system of MRPs to other celestial objects, having a librational nature of rotation as well.In this case, it is assumed that an analytical approach was used to describe the rotational motion, which, on the one hand, is rather complicated for solving the equations of rotation and has a lower accuracy than the numerical one, but, on the other hand, the resulting analytical representation of MRPs provides a more convenient tool for analyzing the behavior of MRPs with changes in various parameters, determining the rotational dynamics of the satellite.It is with its help that the procedure of flexible computer simulation of the rotation process and the identification of those observational manifestations, which are determined primarily by the parameters of the figure of a celestial body - the laws of distribution of body mass, compression, non-sphericity.