In this paper, the time-resolved degenerate four-wave mixing in two-level system is treated rigorously with the density-matrix method. The phenomena related to the saturation effect are discussed. The physical origin of these phenomena is analysed.