The phonon dispersion unfolding method is useful for obtaining hidden Bloch symmetries and comparing theoretical results with experiment spectrums (e.g. inelastic neutron scattering, inelastic X-ray scattering, Raman ). In this paper, we propose a method to unfold phonon dispersions. The main advantage of this method is the ability to handle systems with heavy breaking of spatial translational symmetry. Its validity is tested by pure diamond, diamond with Si substitution, and diamond with C vacancies.Since the discovery of effective energy bands in semiconductors [1,2] and alloys [3], unfolding methods for treating the electronic band structures of the weakly periodic solid systems have been rapidly developed [4][5][6][7][8][9][10][11][12][13][14]. Moreover, it has also stimulated the study of phonon dispersion unfolding method [14][15][16][17]. Like electronic energy bands, the phonon dispersions are also affected by translational symmetry (TS) breaking. Most usually, one has to artificially use a supercell to calculate the phonon dispersions of a TS-broken system. However, the theoretically calculated phonon dispersions of a supercell is much denser than that of the primitive cell. As a result, it can not be directly compared with the phonon dispersions detected in inelastic neutron (or X-ray) scattering experiments. Thus, as a further step, one needs to unfold the phonon dispersions into the Brillouin zone (BZ) of the primitive cell.