The study of acoustic properties within particle systems is of great significance for inverting system structures, designing soundproofing, or sound-guiding materials. To investigate the impact of shape-induced anisotropy on sound velocity, this study employs the discrete element method to simulate the propagation of sound waves in systems of ellipsoids with various shapes and calculates the corresponding sound velocities. The results indicate that the shape of the particles has a significant impact on sound velocity. In the ellipsoidal system defined by + + = 1, 1/alpha indicates that the major axis of the ellipsoid is alpha times the length of the minor axis. The sound velocity varies with alpha: it decreases as alpha increases within the range alpha is an element of [1/7.5, 1/4.5], then increases with alpha in the range alpha is an element of [1/4.5, 1/1.45], and finally decreases again as alpha increases in the range alpha is an element of [1/1.45, 1]. The sound velocity is the lowest for spherical particle systems, and it is highest when alpha = 1/7.5. The differences between the maximum and minimum velocities of compression waves and shear waves are 16 % and 24 %, respectively, indicating that shape-induced anisotropy has a significant and noticeable effect on sound velocity. Additionally, elastic moduli distributions across various shaped systems are computed and analyzed. The fourth-order anisotropic expansion of elastic moduli effectively elucidates the correlation between sound