looking for a solution for an optimization problem

2020 2 5
4 个回答
回答已采纳
更新时间：2020 2 9
3 次查看（30 天）

0 个投票

hi, I am looking for a solution for an optimization problem,
I have a set of non-linear equations and I am looking for positive results, what is the method to follow? an example if possible.
thanks a lot

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

 采纳的回答

Alan Weiss
Alan Weiss 2020-2-5

0 个投票

You might find this example relevant. If you have a recent MATLAB version, then look at this example.
Alan Weiss
MATLAB mathematical toolbox documentation

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Zinedine Gueddal
Zinedine Gueddal 2020-2-5
Thanks, excuse me because I’m a beginner, for exemple: If i look for solving this 3 equations: (X(1)^2)-1=0; (x(2)^2)-4=0; (x(3)^2)-9=0; and i have contrainte: X(1:3)>=0; I maked program:confuneq.m: Function [c, ceq] =confuneq(x) c=-x(1:3); ceq=(X(1)^2)-1=0 (x(2)^2)-4=0 (x(3)^2)-9=0; And exemple.m: X0=[0,0,0]; options=optimoptions('lsqnonlin',FiniteDifferenceType',central'); X=lsqnonlin(@confuneq, x0, [], [], options); The x result [0 0 0]! Spite not satisfate [1 2 3] result why?
Alan Weiss
Alan Weiss 2020-2-5
The problem is your initial point has gradient 0. Take almost any nonzero initial point.
Alan Weiss
MATLAB mathematical toolbox documentation
Zinedine Gueddal
Zinedine Gueddal 2020-2-5
Thanks a lot, is there an other program using fmincon for solving this?
Zinedine Gueddal
Zinedine Gueddal 2020-2-5
I changed the initial value for x0 =[0. 1,0.1,0.1] but the result x =[9.8900e-16,0.1,0.1] not [1 2 3], can you help me please

请先登录，再进行评论。

更多回答（3 个）

Alan Weiss
Alan Weiss 2020-2-6

0 个投票

I think that you need to be more careful when trying to follow the example. EITHER use fmincon OR use lsqnonlin, not both. The fmincon workflow:
fun = @(x)0;
The nonlinear constraint function is nlcon at the end of this example.
Set lower bounds of zero.
lb = zeros(3,1);
Run fmincon starting from [0.1;0.1;0.1];
x0 = 0.1*ones(3,1);
[x,fval,eflag] = fmincon(fun,x0,[],[],[],[],lb,[],@nlcon)
This finds the solution you want.
Alan Weiss
MATLAB mathematical toolbox documentation
function [c,ceq] = nlcon(x)
c = [];
ceq = [x(1)^2 - 1;
x(2)^2 - 4;
x(3)^2 - 9];
end
Zinedine Gueddal
Zinedine Gueddal 2020-2-8

0 个投票

Thank you very much, it is for the first time I see this excellent idea (flag instead of exitflag)!
Zinedine Gueddal
Zinedine Gueddal 2020-2-8

0 个投票

I tried this program for fmincon on its market and I found positive results, but I tried another program but I found results at x near zero, how can I find x and not negative ?

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Alan Weiss
Alan Weiss 2020-2-9
I'm sorry, but I do not understand what you are asking. Are you saying that different input data yields different results, and some of those results are not satisfactory? Or are you saying that you tried running another program and it gave better results? Or are you asking something else?
Alan Weiss
MATLAB mathematical toolbox documentation

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Systems of Nonlinear Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by