1) Augment your parameter vector X with a.
Xaug = [X; a];
2) Define a-b(x,a) < 0 as a nonlinear constraint.
3) Define the bounds of a.
lba = 0;
uba = inf;
4) Solve the problem using fmincon (or another optimization routine).
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!