The localization of the hybrid modes of phonons and photons in polar matter is investigated in the presence of random scatterers theoretically. We employ the self-consistent generalized Born-Huang approach to derive effective equations describing the phonon-polariton fields. Based on these equations, the density of states and various localization properties are exploited in two-dimensional systems both analytically and numerically within the framework of the Anderson model with a non-Hermitian effective Hamiltonian. Consequently, it is shown that the disorder effect brings some intriguing features which include the appearance of the localized states in the polariton bottleneck in the energy spectrum and the collapse of the energy gap. In addition, an analysis is given of the polariton level-spacing distribution.