Optimisation des Q-matrices pour la résolution de contraintes linéaires