The algorithm of dispersion function of Rayleigh waves´ in porous layered half-space

The dispersion functions of Rayleigh waves are theoretically deduced in layered half-space, which consist of elastic and porous layers. Different with the Rayleigh waves of elastic half-space, the dispersion function changes when the porous layer exists, and relates with the position of porous layers. The dispersion functions are given in cases of the porous layers locating at different positions. Furthermore, by the resolution of the original transfer matrix, an optimized square matrix of 15 orders is deduced in order to eliminate the cubic terms of e, this greatly improves the computable frequency range of Rayleigh waves´ dispersion.