The full vibrational spectra especially those high-lying vibrational energies in the dissociation region of four specific electronic states 1(3)Δ(g), 33Σ(+)(g), 1(3)Σ-(g) and b(3)Π(u) have been obtained by using the improved variational algebraic method (VAM). The analytical potential energy functions (APEFs) of these electronic states are also determined with corresponding adjustable parameter λ by using the variational algebraic energy consistent method (VAECM) based on the VAM vibrational spectra. The full vibrational energies, vibrational spectroscopic constants, force constants f(n), and expansion coefficients a(n) of the VAECM potential are also tabulated for each electronic state in this study. The results show that the VAECM analytical potentials are superior to some other widely used analytical ones, and do not have the unphysical tiny barriers existing in the precious AECM potentials.