how to use fmincon in a loop

3 次查看（过去 30 天）
Cong
Cong 2014-2-7
评论： Amit 2014-2-10
I have 5 unknown parameters but only 4 equations, but each parameter has an interval, so I prefer to use fmincon. For example, these are 4 equations:
2*x(3) + x(5) + 5*x(1) - 1=0;
x7*(4) + 3*x(3) + x(2) - 10=0;
x(3) + 2*x(5) + 5*x(1) - 3=0;
3*x(4) + 9*x(5) + x(2) - 10=0;
My idea is like this:
for x(1)=0.1:0.1:1
fmincon(@myfun, x0,[]...)
end
But I do not know how to let x(1) be changed in the loop.
Thank you for any suggestions in advance.
  4 个评论
显示 2更早的评论
Amit
Amit 2014-2-8
3*x(4) + 9*x(6) + x(2) - 10=0;
is this correct? x(6)? I thought you had only 5 parameters?
Cong
Cong 2014-2-8
Sorry, I just correct my mistake. You are right, only 5 parameters. This is just an example.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Amit
Amit 2014-2-8
A = [5 0 2 0 1;
0 1 3 7 0;
5 0 1 0 2;
0 1 0 3 9];
b = [-1;-10;-3;-10];
FF = @(x) norm(A*x+b);
xx = fmincon(FF,rand(5,1),[],[],[],[],zeros(5,1),ones(5,1))
  5 个评论
显示 3更早的评论
Cong
Cong 2014-2-10
Hi Amit, the equations that I posted are not my real problem which are more complex and non-linear. Here I just want to get an idea how to vary my x(5) while using 'fmincon'. I figure out that it is a problem of using anonymous function. Here is a simple example, but note that some contents are not provided.
for p=1:10
i=i+1;
x(i,:) = fmincon(@(X) myfun(X,p), x0, [], [], [], [], ...
[], [], [], options);
end
Thank you for the time.
Amit
Amit 2014-2-10
what is i there?
In a way, lets say you're trying to vary x(5) only, you can do it two way.
  1. set the upper and lower bounds for x(5) as value-eps and value + eps, this will contain the value of x(5) within the range.
  2. For you objective function, use only 4 input variable and pass the value of x(5) as a scalar as in fmincon(@(x) myfun(x,x5val). In this scenario, you'll have to write the function accordingly.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by