Abstract: Aiming at the growth of ablative Rayleigh-Taylor instability with two perturbations, the evolutions of the amplitudes of high-order harmonics excited by two-mode coupling under different preheating conditions are studied by using a high-precision numerical simulation method. When the fundamental modes are a long-wavelength and a short-wavelength mode, the long-wavelength modes of the excited harmonics are dominant, while the development of short-wavelength modes are obviously suppressed; when the fundamental modes are two short-wavelength modes, many fast-growing and long-wavelength modes are excited, and the growth of short-wavelength modes are in the form of small oscillation. By comparing the two different two-mode coupling cases, it is found that the long-wavelength structures are dominant in the weakly nonlinear stage. Especially, in the two short-wavelength modes coupling case, the bubbles and spikes show long-wavelength structures which are different from the two fundamental modes. By further comparing the three preheating ablative effects, it is found that the higher the preheat degree is, the more the coupled harmonics growth will be weakened. It is of great significance to control the development of ablative Rayleigh-Taylor instability in inertial confinement fusion engineering.