To develop a mathematical model of rotational motion of an artificial satellite about its center of mass under the action of various forces (magnetic, Lorentz, etc.) caused by the geomagnetic field, it is necessary to know the induction of the Earth's magnetic field (EMF) as a function of the radius vector of a given point in the near-Earth space. Because the EMF possesses a complex structure and the above-mentioned functional dependence is unavailable in explicit analytic form, a set of approximate models of the EMF should be used. The simplest such model-aright dipole (aligned with the axis of rotation)-does not enable one to reveal in detail the influence of diurnal EMF rotation on the rotational motion of a satellite. The next EMF approximation-an inclined magnetic dipole-does not suffer from the above-mentioned drawback. However, it is shown that not all corrections to the magnetic induction of the EMF of the same order of magnitude are taken into account in the course of transformation from the model of aligned dipole to the model of inclined dipole. So, to develop the EMF model accurately accounting for the absence of axial symmetry of the EMF with respect to the axis of diurnal rotation of the Earth, in general, the effect of the quadrupole component of the geomagnetic potential on the EMF induction (and, probably, even the components of higher orders) must be taken into consideration. By using the International Geomagnetic Reference Field IGRF-2000, the multipole models of the EMF, corresponding to quadrupole, octupole, and higher-order approximations, were constructed and studied in this work. The EMF potential is expressed in terms of its multipole tensors. As a result, projections of the induction and induction gradient of EMF in the center of mass of the satellite onto the axes of the orbital coordinate system can be written in convenient and concise form. The expressions for the first four multipole tensors through the known geomagnetic constants are found. A method for estimating the reliability of these models is put forward, and the regions of applicability of the quadrupole and octupole models are drawn on the plane of orbital parameters.