| 查看: 1607 | 回复: 0 | ||
chenxh97银虫 (著名写手)
|
[求助]
matlab最优化求解fmincon寻求帮助
|
|
运行 [x,fval,exitflag]=fmincon(@opt_fun,x0,A,[100 90 80],Aeq,[20 20],lb,ub,@opt_con_EEVEN)结果有错误，结果应该是：x=[20 0 0 20 0 0]，但是运行的结果与初始值相同x0=[0 0 20 0 0 20]. 运行结果与相关值如下所示： >> [x,fval,exitflag]=fmincon(@opt_fun_EEVEN,x0,A,[100 90 80],Aeq,[20 20],lb,ub,@opt_con_EEVEN) Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 5.741181e-17. > In /usr/local/MATLAB/R2014a/toolbox/optim/optim/private/backsolveSys.p>backsolveSys at 17 In /usr/local/MATLAB/R2014a/toolbox/optim/optim/private/solveKKTsystem.p>solveKKTsystem at 15 In /usr/local/MATLAB/R2014a/toolbox/optim/optim/private/computeTrialStep.p>computeTrialStep at 66 In /usr/local/MATLAB/R2014a/toolbox/optim/optim/barrier.p>barrier at 357 In fmincon at 818 Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 8.881786e-18. > In /usr/local/MATLAB/R2014a/toolbox/optim/optim/private/backsolveSys.p>backsolveSys at 17 In /usr/local/MATLAB/R2014a/toolbox/optim/optim/private/solveKKTsystem.p>solveKKTsystem at 15 In /usr/local/MATLAB/R2014a/toolbox/optim/optim/private/computeTrialStep.p>computeTrialStep at 66 In /usr/local/MATLAB/R2014a/toolbox/optim/optim/barrier.p>barrier at 357 In fmincon at 818 Local minimum found that satisfies the constraints. Optimization completed because the objective function is non-decreasing in feasible directions, to within the default value of the function tolerance, and constraints are satisfied to within the default value of the constraint tolerance. <stopping criteria details> x = 0.0000 0.0000 20.0000 0.0000 0.0000 20.0000 fval = 6.4000e+03 exitflag = 1 其中 function y=opt_fun(x) global cc; y=sum(cc.*x); function [c,ceq]=opt_con_EEVEN(x) c=[]; ceq(1)=x(1)*(x(1)-20); ceq(2)=x(2)*(x(2)-20); ceq(3)=x(3)*(x(3)-20); ceq(4)=x(4)*(x(4)-20); ceq(5)=x(5)*(x(5)-20); ceq(6)=x(6)*(x(6)-20); A = 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 Aeq = 1 0 1 0 1 0 0 1 0 1 0 1 lb = 0 0 0 0 0 0 ub = 100 100 100 100 100 100 x0=0 0 20 0 0 20 cc = 1 1 20 20 300 300 |
回复此楼» 猜你喜欢
生命口会评
已经有5人回复
-
关于如何从代码看上不上会
已经有22人回复
-
大龄残疾硕士的一点执念
已经有25人回复
-
今年E04面上
已经有22人回复
-
chemdraw
已经有6人回复
-
祈祷青基必中
已经有10人回复
-
b口会评
已经有6人回复
-
信息学部会评时间
已经有3人回复
-
在职考研
已经有7人回复
-
MSER送审了还被拒稿
已经有5人回复