The dynamics of spin 3/2 systems is analyzed using the density matrix theory of relaxation. By using the superoperator formalism, an algebraic formulation of the density matrix's evolution is obtained, in which the contributions from free relaxation and RF application are easily factored out. As an intermediate step, an exact form for the propagator of the density matrix for a spin 3/2 system, in the presence of static quadrupolar coupling, inhomogeneous static magnetic field, and relaxation is demonstrated. Using this algebraic formulation, exact expressions for the behavior of the density matrix in the classical one-, two-, and three-pulse experiments are derived. These theoretical formulas are then used to illustrate the bias introduced on the measured relaxation parameters by the presence of large spatial variations in the B0 and B1 fields. The theoretical predictions are easily evaluated through simple matrix algebra and the results agree very well with the experimental observations. This approach could prove useful for the characterization of the spatial variations of the signal intensity in multiple quantum-filtered sodium MRI experiments.