The thermal problem of functionally graded materials (FGM) under linear heat source is studied by a hybrid numerical method. The accuracy of the analytical method and the efficiency of the finite element method are taken into account. The volume fraction of FGM in the thickness direction can be changed by changing the gradient parameters. Based on the weighted residual method, the heat conduction equation under the third boundary condition is established. The temperature distribution of FGM under the action of linear heat source is obtained by Fourier transform. The results show that the closer to the heat source it is, the greater the influence of the heat source is and the influence of the heat source is local. The temperature change trend of the observation points is consistent with the heat source, showing a linear change. The results also show that the higher the value of gradient parameter is, the higher the temperature of location point is. The temperature distribution of observation points is positively correlated with gradient parameter. When the gradient parameter value exceeds a certain value, it has a little effect on the temperature change in the model and the heat conduction in the model tends to be pure metal heat conduction, the optimal gradient parameters combined the thermal insulation property of ceramics and the high strength toughness of metals are obtained.