Brain functions are related to neuronal networks of different sizes and distribution, and neuronal networks of different sizes oscillate at different frequencies. Thus the synchronization of neuronal networks is often reflected by cross-frequency interaction. The description of this cross-frequency interaction is therefore a crucial issue in understanding the modulation mechanisms between neuronal populations. A number of different kinds of interaction between frequencies have been reported. In this paper, we develop a general harmonic wavelet transform based bicoherence using a phase randomization method. This allows us to measure the comodulation of oscillations between different frequency bands in neuronal populations. The performance of the method is evaluated by a simulation study. The results show that the improved wavelet bicoherence method can detect a reliable phase coupling value, and also identify zero bicoherence for waves that are not phase-coupled. Spurious bicoherences can be effectively eliminated through the phase randomization method. Finally, this method is applied to electrocorticogram data recorded from rats during transitions between slow-wave sleep, rapid-eye movement sleep and waking. The phase coupling in rapid-eye movement sleep is statistically lower than that during slow-wave sleep, and slightly less than those in the wakeful state. The degree of phase coupling in rapid-eye movement sleep after slow-wave sleep is greater than in rapid-eye movement sleep prior to waking. This method could be applied to investigate the cross-frequency interactions in other physiological signals.