using fmincon in matlab

Question

Say i have a function f(X) which i want to minimize with constraints such that some other functions- A(X) = 0 and B(X) = 0 and 0 < C(X) < pi. There are many algorithms to do it, but to make my life easier, i want to use built in function fmincon() in matlab. So i read this documentation: http://www.mathworks.com/help/optim/ug/fmincon.html

But I don't understand how I should I pass the parameters to solve my problem in specific. How do I do it? Can I do it at all?

i can also put my C(X) function such that, C(X)>0 – Sahil Chaudhary Jun 06 '13 at 06:22 — Sahil Chaudhary, Jun 06 '13 at 06:22

score 1 · Accepted Answer · answered Jun 06 '13 at 06:33

1

Use the nonlcon parameter of fmincon (I'm assuming here your constraints are nonlinear?). Then A(X) and B(X) are fine but for C(X) it must be in the form c(X) < 0 so you'll need to break it into two constraints of that form.

I pulled this example of how to specify a function for nonlcon from elsewhere in the documentation:

function [c,ceq]=myNonlinearContraints(x)
%First deal with your nonlinear equalities
c(1) = A(X);
c(2) = B(X);
%Then your inequalities transformed to be in the form ceq < 0
ceq(1) = -C(X);
ceq(2) = C(X) - pi;

answered Jun 06 '13 at 06:33

Dan

45,079
17
88
157

quick question, when i pass the arguments, should i pass it like this: fmincon(...myNonlinearContraints)? – Sahil Chaudhary Jun 06 '13 at 06:49
You need to pass it as a function handle so you need to preface it with an `@` i.e. `fmincon(...,@myNonlinearContraints)`. Exactly as you pass your objective function. – Dan Jun 06 '13 at 06:58
got it.and I can pass only one of it? not more i suppose. Say I have defined myNonlinearConstraints() for B and C. How is my fmincon(...) going to look? thanks – Sahil Chaudhary Jun 06 '13 at 06:59
No put all the nonlinear contraints in the same function. The example function above encompases your A, B and C – Dan Jun 06 '13 at 07:03
@SahilChaudhary Sure, glad it helped :) btw if an answer does help you, you should consider giving it an upvote ;) – Dan Jun 06 '13 at 07:26

score 0 · Answer 2 · answered Jun 06 '13 at 06:35

See if the functions A,B are linear or nonlinear. That is maybe A(X) is simply a integral, then the interpretation should be that it is linear. It does actually make a difference if you supply a linear constraint as nonlinear.

If they are nonlinear, then create a

function [c,ceq] = nonlcon(X)

which gives out the equality constraint value (ceq) and inequality constraint value (c). Remember that inequality constraints are interpreted as

ineq(X) < 0

So you need to compute it that way.

C(X) seems to be a trigonometric function, so it will also be part of nonlcon function. This nonlcon, you will pass to fmincon as an argument. Nonlcon is called for a specific value of X and it returns the constraint value. Pass your lower and upper bounds if any and try the optimization for different initial points x0. For some problems, more than one solution can be found.

using fmincon in matlab

2 Answers2