I'm working on a model predictive control problem in which there are constraints both on the input of the system and on the state variables. Since solvers like quadprog or mpcActiveSetSolve deal with the problem J = x'Hx + F'x, allowing to specify constraints only for the x variable (A*x <= b or the equality ones), which is in my case the input of the system for all the steps of the horizon (u bar as in the Cost_Function_LA.png), I don't know how to specify to the solver the constraints I put on the variables of the state.