A matrix formalism is introduced here for the analytical calculation of the thermo-elastic response of one-dimensional optically inhomogeneous layered structures exposed to the periodic laser heating. The multilayer structure consists of thermally and mechanically homogeneous solid layers which can be interleaved by liquid or gaseous interlayers. The method enables the modelling of photoacoustic and photothermal measurements in a variety of engineering applications.