Based on the characteristic matrix of a lamellar structure, a matrix method about thermal wave transmission and reflection coefficients is presented, and a calculating formula for photothermal radiometry signal of the multilayer media is given. The results of two three-ply samples calculated by this method are in good agreement with experimental data.