In this paper,we study the qualitative behaviour of satellite systems using bifurcation diagrams, Poincaré section, Lyapunov exponents, dissipation, equilibrium points, Kaplan–Yorke dimension etc. Bifurcation diagrams with respect to the known parameters of satellite systems are analysed. Poincaré sections with different sowing axes of the satellite are drawn. Eigenvalues of Jacobian matrices for the satellite system at different equilibrium points are calculated to justify the unstable regions. Lyapunov exponents are estimated. From these studies, chaosin satellite system has been established. Solution of equations of motion of the satellite system are drawn in the form of three-dimensional, two-dimensional and time series phase portraits. Phase portraits and time series display the chaotic nature of the considered system.