Воздействие полихроматического поля на двухуровневую среду рассмотрено в асимптотическом пределе малых амплитуд. Найденные поправки третьего порядка позволяют описывать нелинейный отклик на внешнее поле и межмодовое взаимодействие. Полученные спектры поляризации сравнивались с результатами численного решения уравнения матрицы плотности в приближении вращающейся волны. Результаты исследований могут быть использованы в нелинейной комб-спектроскопии.