The inverse thermoelastic problem of identification of the variable properties of a functionally graded rectangle is studied. Unsteady vibrations are excited by applying mechanical and thermal loads to the upper side of the rectangle. To solve the direct problem in Laplace transforms, the method of separation of variables and the shooting method for harmonics are used. Transformants are inverted by expanding the origin in terms of shifted Legendre polynomials. The method proposed for solving the direct problem is verified by comparison with a finite element solution. The influence of the laws of change of variable characteristics on the boundary physical fields is analyzed. The displacement components give additional information on the mechanical loading, and the temperature measured on the upper side of the rectangle over a certain time interval – on the thermal loading. Assuming that the additional information admits expansion in Fourier series, the two-dimensional inverse problem is reduced to one-dimensional problems for various harmonics. The solution of the obtained nonlinear inverse problems is carried out on the basis of an iterative process, at each stage of which, in order to find corrections for thermomechanical characteristics, systems of Fredholm integral equations of the 1st kind are solved. The possibility of simultaneous reconstruction of several characteristics is investigated. The results of computational experiments on the phased reconstruction of thermomechanical characteristics are presented. The influence of the thermomechanical coupling parameter on the results of the thermal stress coefficient reconstruction was clarified.