ABSTRACT. This investigation is continuation of our studies of the geodetic (relativistic) rotation of the Solar system bodies (Eroshkin and Pashkevich, 2007) and (Eroshkin and Pashkevich, 2009). For each body (the Moon, the Sun, the major planets and Pluto) the files of the values of the components of the angular velocity of the geodetic rotation are constructed over the time span from AD1000 to AD3000 with one day spacing, by using DE422/LE422 ephemeris (Folkner, 2011), with respect to the proper coordinate systems of the bodies (Seidelmann et al., 2005). For the first time in the perturbing terms of the physical librations for the Moon and in Euler angles for other bodies of the Solar system the most essential terms of the geodetic rotation are found by means of the least squares method and spectral analysis methods.