An algorithm for the computation of polyhedral invariant sets for closed-loop linear MPC systems