To investigate the polarization characteristics of snowflake-like ice crystals under terahertz waves (1.5 THz, 2.524 THz, 3.437 THz), the Discrete Dipole Approximation (DDA) method is employed to calculate their linear polarization degree, circular polarization degree, and the fraction of circularly polarized light intensity. These results are compared with those of hexagonal plate-like ice crystals of the same equivalent radius. The findings indicate that the polarization characteristics of snowflake-like ice crystals are contingent upon the particle′s equivalent radius, shape, and the frequency of the electromagnetic wave. For particles of identical size, the polarization characteristics of snowflake-like ice crystal particles are significantly more affected at a frequency of 3.437 THz. Within the same frequency band, the larger the size of the snowflake-like ice crystal particles, the more pronounced the variation in polarization characteristics with respect to the scattering angle, with 90°emerging as a pivotal boundary for the scattering angle and the alteration in polarization characteristics. When hexagonal plate-shaped ice crystals are used as substitutes for snowflake-shaped ones, the degree of circular polarization would be significantly underestimated at 1.5 THz, whereas the degree of linear polarization would be overestimated at 3.437 THz. These research outcomes are of significant importance for gaining a deeper understanding of the mechanisms underlying the interaction between terahertz waves and ice crystal particles.