This paper deals with the problem of model reduction by moment matching for linear differential inclusions. The problem is formally formulated and the notions of momentset, perturbed moment trajectory, approximate reduced order model and robust reduced order model are introduced. Two sets of results are presented. The first part of the paper deals with robustness of the reduced order models with respect to input perturbations. Exploiting this result an enhanced model reduction scheme for linear differential equations is presented. In the second part of the paper we focus on the problem of model reduction by moment matching for time-varying systems driven by time-varying signal generators. Finally, these two sets of results are used to solve the problem of model reduction for linear differential inclusions. The results are illustrated by means of numerical examples.