A numerical solution for the polarization of two-level atoms that interact with a polyharmonic field was obtained. An analytical solution for the particular case of a symmetrical position of carrier frequency relative to transition frequency is possible. The results showed that nonlinear effects occur in a polarization spectrum at the small amplitudes of comb components and small frequency distances between them. Consequently, it is necessary to take into account nonlinear effects in comb spectroscopy.