An experimental study has been performed on vortex pairing under fundamental and subharmonic forcing with controlled initial phase differences through hot-wire measurements and a multi-smoke wire flow visualization. For the range of St$_{D}$ < 0.6, vortex pairing was controlled by means of fundamental and subharmonic forcing with varying initial phase differences. Much larger mixing rate was achieved by two-frequency forcing with a proper phase difference than one frequency forcing. As St$_{D}$ decreased, vortex pairing was limited to a narrow region of the initial phase difference between two disturbances and higher amplitudes of the fundamental and its subharmonic at the nozzle exit were required for more stable pairing. As the amplitude of the subharmonic at the nozzle exit increased for fixed St$_{D}$ and fundamental amplitude, the distribution of the subharmonic mode against the variation of the initial phase difference changed from a sine function form into a cusp-like form. Thus, vortex pairing can be controlled more precisely for the former case. For St$_{D}$ > 0.6, non-pairing advection of vortices due to the improper phase difference was sometimes observed in several fundamental forcing amplitudes when only the fundamental was applied. However, when its subharmonic was added, vortex pairing readily occurred. As the initial amplitude of this subharmonic increased, the position of vortex pairing moved upstream. This was thought to be due to the fact that the variation of the initial phase difference between the fundamental and its subharmonic has less effects on vortex pairing in the region of fundamental-only vortex pairing.pairing.