Variational and numerical study of soliton switching in two and three core couplers with higher-order perturbations