The Embedded Atom Method is applied to study the properties of surface phonons of transition metal Cu. We derive the formula of surface lattice dynamical matrix from the se-miempirical functions. As an example, we calculate the surface phonon dispersions on Cu (100) surface; the calculated Raileigh phonon dispersion curve on Cu(100) along the direction ΓX is compared with that measured by electon-energy-loss-spectroscopy, and satifactory agreement is obtained.