Thread Subject: nonlinear constraint help

Hello to all,

I have two maybe very dumb questions but I cannot get my head around these: I have a problem subject to 18 linear constraints and 1 nonlinear constraint. In this last one appear in part variables present also in the linear ones and other variables.
 Questions:
1. from the nonlinear constraint to recall the variables present in the linear constraints, I'm calling them with x(position), given the example here http://www.mathworks.com/access/helpdesk/help/toolbox/optim/index.html?/access/helpdesk/help/toolbox/optim/ug/brhkghv-7.html&http://www.mathworks.com/matlabcentral/newsreader/view_thread/256555 is this correct? the name of the solution in my function where I call fmincon is x.

2. the second question is if I can add more variables to be solved for in the nonlinear constraint, but these would be only in this constraint (they are multiplied by the other variables present also in the linear constraints so they should be identified, isn't it?). I have seen the documentation on how to pass parameters in the nonlinear constraints but not variables, is this possible?

Thanks a lot for your help in advance,

Best Regards,

Simone

I am sorry, I do not understand what you are asking. Perhaps you can
include an example of the type of parameters you would like to include,
and an example of the syntax you use to call your optimization function.

Alan Weiss
MATLAB mathematical toolbox documentation

Simone wrote:
> Hello to all,
>
> I have two maybe very dumb questions but I cannot get my head around these: I have a problem subject to 18 linear constraints and 1 nonlinear constraint. In this last one appear in part variables present also in the linear ones and other variables.
> Questions:
> 1. from the nonlinear constraint to recall the variables present in the linear constraints, I'm calling them with x(position), given the example here http://www.mathworks.com/access/helpdesk/help/toolbox/optim/index.html?/access/helpdesk/help/toolbox/optim/ug/brhkghv-7.html&http://www.mathworks.com/matlabcentral/newsreader/view_thread/256555 is this correct? the name of the solution in my function where I call fmincon is x.
>
> 2. the second question is if I can add more variables to be solved for in the nonlinear constraint, but these would be only in this constraint (they are multiplied by the other variables present also in the linear constraints so they should be identified, isn't it?). I have seen the documentation on how to pass parameters in the nonlinear constraints but not variables, is this possible?
>
> Thanks a lot for your help in advance,
>
> Best Regards,
>
> Simone

"Simone " <simonematool@gmail.com> wrote in message <h8ufpv$2iq$1@fred.mathworks.com>...
> Hello to all,
>
> I have two maybe very dumb questions but I cannot get my head around these: I have a problem subject to 18 linear constraints and 1 nonlinear constraint. In this last one appear in part variables present also in the linear ones and other variables.
> Questions:
> 1. from the nonlinear constraint to recall the variables present in the linear constraints, I'm calling them with x(position), given the example here http://www.mathworks.com/access/helpdesk/help/toolbox/optim/index.html?/access/helpdesk/help/toolbox/optim/ug/brhkghv-7.html&http://www.mathworks.com/matlabcentral/newsreader/view_thread/256555 is this correct?
====================

Yes.


>
> 2. the second question is if I can add more variables to be solved for in the nonlinear constraint, but these would be only in this constraint
==============

Yes.

 (they are multiplied by the other variables present also in the linear constraints so they should be identified, isn't it?). I have seen the documentation on how to pass parameters in the nonlinear constraints but not variables, is this possible?
================

If your variables decompose as [x,y] where x and y are vectors, then neither your objective function f(x,y) nor your constraints Ci(x,y) have to have an explicit dependence on both x and y, even if they both take x and y as input arguments.


However, if f(x,y) is independent of y, i.e. can theoretically be written f(x), that gives the algorithm more freedom in finding a y that satisfies Ci(x,y)>=0. If it's important to distinguish between more than one possible solution for y, you must of course add more constraints....

Dear Alan and all,

thanks a lot for the answer and sorry for my little clarity, I am gonna try again.
So the idea is that I have a series of linear constraints called, as usual, with the A= matrix, b=vector, Aeq=matrix and beq= vector. I never specified in my use of MATLAB (which until recently has been mostly with linear programming) the name of the variables that satisfy these constraints, then I have a nonlinear constraint in which are both the variables present in the linear constraints and other objects (that I am interested in as a solution) in something like this:
function [c, ceq] = myconstr(x,l0,l)
c = [l(1,1)./l0.*x(2).*x(2*n+2)+...+l(1,n)./l0.*x(n+1).*x(3*n+1) -x(2*n+2).*x(n+2)- x(3*n+1).*x(2*n+1) ]
ceq=[]
l0 and l are data given. I assumed that the variables of the system in the linear constraints are all in the vector x as given in the example from the Mathworks website (and I'm calling also my solution x as per below command). In the linear part I have 2*n+1 variables in the constraints so I started calling these new variables I would like to pass onto the system with x(2*n+2) until x(3*n+1), is it possible to do this?
The following is the way I call the optimization function:
[x, fval, exitflag]=fmincon(@myobj,xs,A,b,Aeq,beq,lb,ub,@myconstr);

Hoping this is clearer and can help you to help me,

Hoping you have a great weekend,

Best Regards,

Simone




Alan Weiss <aweiss@mathworks.com> wrote in message <h8vuq3$5ng$2@fred.mathworks.com>...
> I am sorry, I do not understand what you are asking. Perhaps you can
> include an example of the type of parameters you would like to include,
> and an example of the syntax you use to call your optimization function.
>
> Alan Weiss
> MATLAB mathematical toolbox documentation
>
> Simone wrote:
> > Hello to all,
> >
> > I have two maybe very dumb questions but I cannot get my head around these: I have a problem subject to 18 linear constraints and 1 nonlinear constraint. In this last one appear in part variables present also in the linear ones and other variables.
> > Questions:
> > 1. from the nonlinear constraint to recall the variables present in the linear constraints, I'm calling them with x(position), given the example here http://www.mathworks.com/access/helpdesk/help/toolbox/optim/index.html?/access/helpdesk/help/toolbox/optim/ug/brhkghv-7.html&http://www.mathworks.com/matlabcentral/newsreader/view_thread/256555 is this correct? the name of the solution in my function where I call fmincon is x.
> >
> > 2. the second question is if I can add more variables to be solved for in the nonlinear constraint, but these would be only in this constraint (they are multiplied by the other variables present also in the linear constraints so they should be identified, isn't it?). I have seen the documentation on how to pass parameters in the nonlinear constraints but not variables, is this possible?
> >
> > Thanks a lot for your help in advance,
> >
> > Best Regards,
> >
> > Simone

Thanks a lot Matt for your answer, as you probably saw I was asnswering Alan already so only when I had already posted my reply I saw yours, sorry for the double post but I thank you for your answer, but then, if adding more variables can be done, do you think can it be that MATLAB gives me an error saying that the index exceeds matrix dimension in the nonlinear constraint because I add more variables in the nonlinear constraint or am I doing something else wrong??? Like the number of columns of A should be equal to the total number of variables, does it matter???

Thanks a lot again,

Best Regards,

Simone




"Matt " <xys@whatever.com> wrote in message <h90802$61u$1@fred.mathworks.com>...
> "Simone " <simonematool@gmail.com> wrote in message <h8ufpv$2iq$1@fred.mathworks.com>...
> > Hello to all,
> >
> > I have two maybe very dumb questions but I cannot get my head around these: I have a problem subject to 18 linear constraints and 1 nonlinear constraint. In this last one appear in part variables present also in the linear ones and other variables.
> > Questions:
> > 1. from the nonlinear constraint to recall the variables present in the linear constraints, I'm calling them with x(position), given the example here http://www.mathworks.com/access/helpdesk/help/toolbox/optim/index.html?/access/helpdesk/help/toolbox/optim/ug/brhkghv-7.html&http://www.mathworks.com/matlabcentral/newsreader/view_thread/256555 is this correct?
> ====================
>
> Yes.
>
>
> >
> > 2. the second question is if I can add more variables to be solved for in the nonlinear constraint, but these would be only in this constraint
> ==============
>
> Yes.
>
> (they are multiplied by the other variables present also in the linear constraints so they should be identified, isn't it?). I have seen the documentation on how to pass parameters in the nonlinear constraints but not variables, is this possible?
> ================
>
> If your variables decompose as [x,y] where x and y are vectors, then neither your objective function f(x,y) nor your constraints Ci(x,y) have to have an explicit dependence on both x and y, even if they both take x and y as input arguments.
>
>
> However, if f(x,y) is independent of y, i.e. can theoretically be written f(x), that gives the algorithm more freedom in finding a y that satisfies Ci(x,y)>=0. If it's important to distinguish between more than one possible solution for y, you must of course add more constraints....

Simone, I am still not 100% sure I understand you, but I will make a
guess at what you are asking.

You might be asking how to include the l0 and l matrices in your
constraint function. You should include them exactly as described in the
documentation:

function [c, ceq] = myconstr(x,l0,l)
c = [l(1,1)./l0.*x(2).*x(2*n+2)+...+l(1,n)./l0.*x(n+1).*x(3*n+1)
-x(2*n+2).*x(n+2)- x(3*n+1).*x(2*n+1) ]
ceq=[]

Load l and l0 into your MATLAB workspace, and define
constr = @(x)myconstr(x,l0,l)

Use constr as follows:

[x, fval, exitflag]=fmincon(@myobj,xs,A,b,Aeq,beq,lb,ub,constr);

If this is not what you were asking, please ask again, using more
examples of exactly what you are trying to do and your syntax.

Alan Weiss
MATLAB mathematical toolbox documentation

Simone wrote:
> Dear Alan and all,
>
> thanks a lot for the answer and sorry for my little clarity, I am gonna try again.
> So the idea is that I have a series of linear constraints called, as usual, with the A= matrix, b=vector, Aeq=matrix and beq= vector. I never specified in my use of MATLAB (which until recently has been mostly with linear programming) the name of the variables that satisfy these constraints, then I have a nonlinear constraint in which are both the variables present in the linear constraints and other objects (that I am interested in as a solution) in something like this:
> function [c, ceq] = myconstr(x,l0,l)
> c = [l(1,1)./l0.*x(2).*x(2*n+2)+...+l(1,n)./l0.*x(n+1).*x(3*n+1) -x(2*n+2).*x(n+2)- x(3*n+1).*x(2*n+1) ]
> ceq=[]
> l0 and l are data given. I assumed that the variables of the system in the linear constraints are all in the vector x as given in the example from the Mathworks website (and I'm calling also my solution x as per below command). In the linear part I have 2*n+1 variables in the constraints so I started calling these new variables I would like to pass onto the system with x(2*n+2) until x(3*n+1), is it possible to do this?
> The following is the way I call the optimization function:
> [x, fval, exitflag]=fmincon(@myobj,xs,A,b,Aeq,beq,lb,ub,@myconstr);
>
> Hoping this is clearer and can help you to help me,
>
> Hoping you have a great weekend,
>
> Best Regards,
>
> Simone
>
>
>
>
> Alan Weiss <aweiss@mathworks.com> wrote in message <h8vuq3$5ng$2@fred.mathworks.com>...
>> I am sorry, I do not understand what you are asking. Perhaps you can
>> include an example of the type of parameters you would like to include,
>> and an example of the syntax you use to call your optimization function.
>>
>> Alan Weiss
>> MATLAB mathematical toolbox documentation
>>
>> Simone wrote:
>>> Hello to all,
>>>
>>> I have two maybe very dumb questions but I cannot get my head around these: I have a problem subject to 18 linear constraints and 1 nonlinear constraint. In this last one appear in part variables present also in the linear ones and other variables.
>>> Questions:
>>> 1. from the nonlinear constraint to recall the variables present in the linear constraints, I'm calling them with x(position), given the example here http://www.mathworks.com/access/helpdesk/help/toolbox/optim/index.html?/access/helpdesk/help/toolbox/optim/ug/brhkghv-7.html&http://www.mathworks.com/matlabcentral/newsreader/view_thread/256555 is this correct? the name of the solution in my function where I call fmincon is x.
>>>
>>> 2. the second question is if I can add more variables to be solved for in the nonlinear constraint, but these would be only in this constraint (they are multiplied by the other variables present also in the linear constraints so they should be identified, isn't it?). I have seen the documentation on how to pass parameters in the nonlinear constraints but not variables, is this possible?
>>>
>>> Thanks a lot for your help in advance,
>>>
>>> Best Regards,
>>>
>>> Simone

"Simone " <simonematool@gmail.com> wrote in message <h908vv$b5f$1@fred.mathworks.com>...
> Dear Alan and all,
>
> thanks a lot for the answer and sorry for my little clarity, I am gonna try again.
> So the idea is that I have a series of linear constraints called, as usual, with the A= matrix, b=vector, Aeq=matrix and beq= vector. I never specified in my use of MATLAB (which until recently has been mostly with linear programming) the name of the variables that satisfy these constraints, then I have a nonlinear constraint in which are both the variables present in the linear constraints and other objects (that I am interested in as a solution) in something like this:
> function [c, ceq] = myconstr(x,l0,l)
> c = [l(1,1)./l0.*x(2).*x(2*n+2)+...+l(1,n)./l0.*x(n+1).*x(3*n+1) -x(2*n+2).*x(n+2)- x(3*n+1).*x(2*n+1) ]
> ceq=[]
> l0 and l are data given. I assumed that the variables of the system in the linear constraints are all in the vector x as given in the example from the Mathworks website (and I'm calling also my solution x as per below command). In the linear part I have 2*n+1 variables in the constraints so I started calling these new variables I would like to pass onto the system with x(2*n+2) until x(3*n+1), is it possible to do this?
=========================

No.

You must define both myobj(X) and myconstr(X) so that the input vector X is of length 3*n+1 and contains all variables that ever appear in the problem. It is not necessary that myobj() and myconstr actually use all 3*n+1 components of X. It is okay if they just use the first 2*n+1 of them.

You must also modify A and Aeq so that they have 3*n+1 columns. The columns 2*n+2 through 3*n+1 should be zero.

Similarly, you must extend ub and lb to be of length 3*n+1...

Dear Matt,
thanks a lot again, I will try!
Best regards,
Simone


"Matt " <xys@whatever.com> wrote in message <h90h6h$59s$1@fred.mathworks.com>...
> "Simone " <simonematool@gmail.com> wrote in message <h908vv$b5f$1@fred.mathworks.com>...
> > Dear Alan and all,
> >
> > thanks a lot for the answer and sorry for my little clarity, I am gonna try again.
> > So the idea is that I have a series of linear constraints called, as usual, with the A= matrix, b=vector, Aeq=matrix and beq= vector. I never specified in my use of MATLAB (which until recently has been mostly with linear programming) the name of the variables that satisfy these constraints, then I have a nonlinear constraint in which are both the variables present in the linear constraints and other objects (that I am interested in as a solution) in something like this:
> > function [c, ceq] = myconstr(x,l0,l)
> > c = [l(1,1)./l0.*x(2).*x(2*n+2)+...+l(1,n)./l0.*x(n+1).*x(3*n+1) -x(2*n+2).*x(n+2)- x(3*n+1).*x(2*n+1) ]
> > ceq=[]
> > l0 and l are data given. I assumed that the variables of the system in the linear constraints are all in the vector x as given in the example from the Mathworks website (and I'm calling also my solution x as per below command). In the linear part I have 2*n+1 variables in the constraints so I started calling these new variables I would like to pass onto the system with x(2*n+2) until x(3*n+1), is it possible to do this?
> =========================
>
> No.
>
> You must define both myobj(X) and myconstr(X) so that the input vector X is of length 3*n+1 and contains all variables that ever appear in the problem. It is not necessary that myobj() and myconstr actually use all 3*n+1 components of X. It is okay if they just use the first 2*n+1 of them.
>
> You must also modify A and Aeq so that they have 3*n+1 columns. The columns 2*n+2 through 3*n+1 should be zero.
>
> Similarly, you must extend ub and lb to be of length 3*n+1...