In order to understand the relation between the harmonic spectra efficiency of H+2 and its isotope molecule with laser wavelength, the relation between the harmonic intensity of H+2 and D+2 with the wavelength in the range of 600nm~1600nm was theoretically studied by solving 2-D time-dependent Schrdinger equation. It is shown that the intensity of harmonic spectrum is decreased as the wavelength increases. In shorter wavelength region, the decrease rate of harmonic intensity of H+2 is greater than that of D+2. In longer wavelength region, the decrease rate of harmonic intensity of H+2 is smaller than that of D+2. Furthermore, driven by lower laser intensity, the harmonic yield of H+2 is always higher than that of D+2. Driven by stronger laser intensity, the harmonic yield of H+2 is lower than that of D+2 in shorter wavelength region; while, it is higher than that of D+2 in longer wavelength region. Theoretical analyses show that the extension of nuclear distance and charge resonance enhanced ionization play the important role in the change of harmonic yield of H+2 and D+2. The results are helpful for molecular harmonic control.