The formulas for atomic displacements and Hamiltonian of a thin crystal film in phonon occupation number representation are obtained with the aid of Greens function theory. On the basis of these results, the formulas for thermal expansion coefficients of the thin crystal film are derived with the perturbation theory, and the numerical calculations are carried out. The results show that the thinner films have larger thermal expansion coefficients.