You can implement the constraints on psi as a linear inequality constraint. If you check the documentation for your solver (e.g., doc fmincon), you should see the definitions of the A, Aeq, B, Beq matrices, which you can use to implement the constraint. For example, if your constraint is that
psi < b
You can convert this constraint to an inequality:
psi - b < 0
implement this as:
A = [1 -1]; B = 0;
%assuming your variables are in the order [psi; b].
Similarly, for psi > b:
A = [1 -1]; B = 0;
Note that you will need to pad out A and B with zeros if you have additional parameters to optimize.
