In the present thesis we study the unitary dynamics and the thermalization properties of free-fermion-like Hamiltonians after a sudden quantum quench in presence of disorder. With analytical and numerical arguments, we show that the existence of a stationary state and its description with a generalized Gibbs ensemble (GGE) depend crucially on the observable considered (local versus extensive, one-body versus many-body) and on the localization properties of the final Hamiltonian. We then show an extension of the Wang-Landau algorithm which allows the computation of weighted distributions associated to quantum quenches, like the diagonal and the GGE ensemble expectation-value distributions. We present results on three one-dimensional models, the Anderson model, a disordered one-dimensional fermionic chain with long-range hopping, and the disordered Ising/XY spin chain.