MATLAB Answers

0

Updating bounds in fmincon

Asked by Michelle on 13 Aug 2013
Hello,
I am trying to update bounds on a parameter I am optimizing. I have two parameters, b and psi. psi's bounds are a function of b so they need to be updated. How do I keep changing the bounds while also optimizing for both b and psi?
This is the inequality I need to define:
π-2tan^(-1)(2r/b)<Ψ<3π/2-2tan^(-1)(2r/b)

  0 Comments

Sign in to comment.

Tags

1 Answer

Answer by Matt Kindig on 13 Aug 2013
Edited by Matt Kindig on 13 Aug 2013

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.

  1 Comment

Matt Kindig on 13 Aug 2013
Now that you've edited the question, I see that your constraints are nonlinear. Thus, instead of using the linear inequality method I've described below, you will need to implement this using a custom nonlinear constraint function, with two inequality expressions. See the documentation for fmincon for details.

Sign in to comment.