This paper reports the eigenstructure of a set of first-order hyperbolic partial differential equations for modeling waves in solids with a trigonal 32 symmetry. The governing equations include the equation of motion and partial differentiation of the elastic constitutive relation with respect to time. The result is a set of nine, first-order, fully coupled, hyperbolic partial differential equations with velocity and stress components as the unknowns. Shown in the vector form, the model equations have three coefficient matrices. The wave physics are fully described by the eigenvalues and eigenvectors of these matrices; i.e., the nontrivial eigenvalues are the wave speeds, and a part of the corresponding left eigenvectors represents wave polarization. For a wave moving in a certain direction, three wave speeds can be identified by calculating the eigenvalues of the coefficient matrix in a rotated coordinate system. In this process, without using the plane-wave solution, we recover the Christoffel matrix and thus validate the formulation. To demonstrate this approach, two- and three-dimensional slowness profiles of quartz are calculated. Wave polarization vectors for wave propagation in several compression directions as well as noncompression directions are discussed.