Complexes of localized states are numerically analyzed in two dynamical systems: directly driven nonlinear Schro¨dinger equation (NLS) and double sine-Gordon equation (2SG). Both systems have a wide range of physical applications. Our numerical approach is based on the numerical continuation with respect to the control parameters of the quiescent (stationary) solutions and stability and bifurcation analysis of the linearized eigenvalue problem. Multisoliton complexes of the NLS equation are studied in the undamped and the weak damping regimes. We show that in the weak damping case the directly driven NLS equation holds stable and unstable multi-soliton complexes. The results are confirmed by means of direct numerical simulations of the time-dependent NLS equation. Properties of the multi-fluxon solutions of 2SG equation are studied depending on the parameter of the second harmonic. We show that the second harmonic changes properties and increases the complexity of coexisting static fluxons of 2SG equation. Results are discussed within the frame of the long Josephson junction model.