Abstract: A recursive matrix method is adopted for the computation of dyadic Green＇ s functions in horizontally stratified media with arbitrary number of layers. Three matrix equation groups for computing the coefficients of the Sommerfeld integrals are obtained according to the continuity condition of electric and magnetic fields across the interface between different layers, which are in correspondence with the TM wave produced by a vertical unit electric dipole and the TE or TM wave produced by a horizontal unit electric dipole, respectively. All the equation groups can be quickly solved via the recursive method. The dyadic Green＇s functions with source point and field point be- ing in any layer can be conveniently obtained by merely changing the position of the elements within the source term of the equation groups. The expression of the dyadic Green＇s functions provided by this paper is terse in form and does not have exponentially increasing term, so it does not overflow during computation.

Key words: dyadic Green＇s functions, stratified media, recursive matrix methed, Sommerfeld integrals