Differential density matrix equations for a three-level atomic system driven by polychromatic fields are solved numerically. The polarization spectrum of the probe field is obtained. Ultra-narrow resonances in the polarization spectrum on the probe field frequency emerge at multiple harmonics and subharmonics.