Elucidating the evolution law of the elastic properties of the matrix phase is of great significance for the control of steel properties and quality during continuous casting and subsequent heat treatment. In this paper, thermal expansion experiments and ab initio calculations are used to study the elastic properties of the interstitial free (IF) steel matrix phase in different magnetic states and crystal structures. The results show that the bulk modulus B and the tetragonal shear elastic constant C' for the entire temperature range decrease with increasing temperature, but C-44 is the opposite. While from paramagnetic (PM) to ferromagnetic (FM) state, C'(C-44) have changed similar to 188% (similar to 27%), B increases by similar to 55% during the crystal structure change (fcc -> bcc). With the FM to PM state, the Zener anisotropy parameter increases sharply, and Young's modulus decreases significantly in the [001] direction; the maximum difference is similar to 76 GPa. The evolution rate of average Young's modulus in single bcc-phase FM (fcc-phase PM) range reaches similar to 5.5(similar to 5.6) x 10(-2) GPa K-1. The research provides an effective method for ab initio calculation of the elastic properties of interstitial free and ultra-low carbon steels at high temperature, also furnishing a basis for the application of ab initio calculations to the high temperature performance of steel materials.