A generalized harmonic balance method for forced non-linear oscillations: the subharmonic cases