Based on the theory of surface lattice dynamics, the surface phonon spectrums along three symmetrical directions of \bar \varGamma \bar L, \bar L\bar M and \bar \varGamma \bar M are simulated for the W(100) surface by using the modified analytic embedded atom method. The polarization vectors at different symmetrical points are also calculated. According to the criterion and marking method of surface mode, the surface modes along different symmetrical directions are drawn, the distribution range and mode coupling of surface modes are discussed as well. The vibration frequencies of surface modes calculated by us have been compared to available experimental datum and some theoretical values correspondingly. The results display that the present results are general agreement with the referenced experimental or theoretical results. Based on the calculated polarization vector, the surface vibration states are constructed for the atomic layers in the neighboring surface. And the polarization and local features of the surface modes along different symmetrical directions are analyzed. The results show that there are some coupling phenomena between surface mode dispersion, such as avoid crossing and independence crossing. The avoid crossing is found between the surface-mode branch S₁ and the surface-mode branch S₂ near \bar \zeta _y = 0.32 along \bar L\bar M direction. In the region, going from \bar L to \bar M, S₁ changes from y polarization to z polarization, and S₂ changes from z polarization to y polarization. The independence crossings exist between surface-mode branch S₁ and surface-mode branch S₂ at \bar \zeta _x = 0.5 along \bar \varGamma \bar L direction, and surface-mode branch S₂ and surface-mode branch S₃ at \bar \zeta _x = 0.5 along \bar L\bar M direction, respectively. Before and after the crossings, the polarization and local features of the surface modes have not changed. Inspection of the polarization vectors, the coupling phenomena are iconically demonstrated.