The onset of matrix steady-state cracking stands as a pivotal mechanical characteristic in fiber reinforced ceramic matrix composites (FRCMCs), garnering substantial attention and investigations. Experimentally, it has been demonstrated that increasing interphase layer thickness may cause non-monotonic changes in matrix cracking stress. However, the existing models can hardly elucidate this phenomenon thoroughly due to the neglect of interphase thickness. This paper presents a comprehensive analytical model for the matrix cracking incorporating interphase, the Poisson effect, Coulomb friction, fiber asperities, residual thermal stress (RTS), and their coupling effects, along with a modified criterion for interfacial debonding that accounts for the presence of the axial RTS. Based on the proposed model, three distinct cracking domains, i.e., perfectly bonded, debonding with and without interfacial separation, have been identified with the critical conditions deduced analytically. Thereby the mechanism of the non-monotonic influence of interphase thickness is thoroughly revealed as the transition of cracking modes. Meanwhile, the role of interphase on the matrix cracking is systematically studied, and the results indicate that interphase has a notable effect through relieving axial RTS, adjusting interfacial friction, altering interfacial shear modulus, and influencing debonding toughness. The outcomes of this study offer valuable guidance for the interphase design of FRCMCs.