| 查看: 1617 | 回复: 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 |
回复此楼» 猜你喜欢
只有每年这种时候来逛逛小木虫
已经有8人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
