A method for performing quantum state tomography for quadrupole nuclei is presented in this paper. First, it is shown that upon appropriate phase cycling, the nuclear-magnetic-resonance (NMR) intensities of quadrupole nuclei depend only on diagonal elements of the density-matrix. Thus, a method for obtaining the density-matrix elements, which consists of dragging off-diagonal elements into the main diagonal using fine phase-controlled selective radio-frequency pulses, was derived. The use of the method is exemplified through ${}^{23}\mathrm{Na}$ NMR (nuclear spin $I=3∕2$) in a lyotropic liquid crystal at room temperature, in three applications: (a) the tomography of pseudopure states, (b) the tomography of the quadrupole free evolution of the density matrix, and (c) the unitary state evolution of each qubit in the system over the Bloch sphere upon the application of a Hadamard gate. Further applications in the context of pure NMR and in the context of quantum information processing, as well as generalizations for higher spins, are discussed.