| 查看: 1821 | 回复: 4 | ||
qinzhen724新虫 (小有名气)
|
[求助]
Matlab中用fmincon求解最小值问题，调试了好多次，一直找不到可行解，求大神指导 已有2人参与
|
|
用fmincon求解最小值问题（1个目标函数、9个变量、6个线性等式约束），一直都是No feasible solution，不知道本身无解还是什么问题。 x0=[4000;4000;2000;500;100;10;10;10;10]; A=[]; b=[]; Aeq=[-88.840439250000020, 20.679950666666663, -2.158236236666667e+02, -3.601469919166667e+02, -36.694763666666690, 21.462819666666660, -13.261596999999995, 7.477251833333336, -1.023975406666667e+02; 0.012, 0, 0, 0.012, 0.012, 0, 0, 0, 0.012; 0, 0.002, 0.002, 0, 0.004, 0, 0.003, 0.002, 0; 0.016, 0, 0.016, 0.032, 0, 0, 0, 0, 0.016; 0, 0, 0, 0, 0, 0.028, 0.014, 0, 0; 0, 0, 0, 0, 0, 0, 0, 0.032, 0.032] beq=[4588048.12066667; 171.600000000000; 21.5666666666667; 241.933333333333; 2.10000000000000; 0.600000000000000] VLB=[0.001;0.001;0.001;0.001;0.001;0.001;0.001;0.001;0.001]; VUB=[6601.833333;4978.305556;4978.305556;5172.684028;2489.152778;34.625;69.25;8.65625;8.65625]; options = optimoptions(@fmincon,'Algorithm','interior-point','MaxIter',5000); [x,fval,exitflag,output]=fmincon(@Gibbsfun,x0,A,b,Aeq,beq,VLB,VUB,[],options) %%目标函数为myobjfun function G=myobjfun(x) deltaG1=-8.841665127258336e+04; deltaG2=-1.661952727160000e+05; deltaG3=5.504948170966668e+04; deltaG4=-1.666712325249982e+03; deltaG5=1.092642304030001e+05; deltaG6=-2.281491138470000e+05; deltaG7=1.169078584030000e+05; deltaG8=1.087887085851667e+05; deltaG9=4.941138012599998e+04; G=x(1)*(deltaG1+8314*log(x(1)/sum(x)*1.2))+x(2)*(deltaG2+8314*log(x(2)/sum(x)*1.2))... +x(3)*(deltaG3+8314*log(x(3)/sum(x)*1.2))+x(4)*(deltaG4+8314*log(x(4)/sum(x)*1.2))... +x(5)*(deltaG5+8314*log(x(5)/sum(x)*1.2))+x(6)*(deltaG6+8314*log(x(6)/sum(x)*1.2))... +x(7)*(deltaG7+8314*log(x(7)/sum(x)*1.2))+x(8)*(deltaG8+8314*log(x(8)/sum(x)*1.2))... +x(9)*(deltaG9+8314*log(x(9)/sum(x)*1.2)) 计算结果如下： No feasible solution found. fmincon stopped because the size of the current step is less than the default value of the step size tolerance but constraints are not satisfied to within the default value of the constraint tolerance. <stopping criteria details> x = 1.0e+03 * 4.666024308553983 3.857233354391819 2.387214185792911 0.000001000000005 0.000001000000001 0.015465337002847 0.000001000000003 0.008656249999525 0.006364985001793 fval = -1.006342210855910e+09 exitflag = -2 output = iterations: 70 funcCount: 710 constrviolation: 5.438285104570996e+06 stepsize: 9.313859771460607e-10 algorithm: 'interior-point' firstorderopt: 0.043644013275751 cgiterations: 76 message: 'No feasible solution found. |
回复此楼» 猜你喜欢
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有11人回复
-
闲聊
已经有3人回复
-
江苏省自然基金 什么时候出结果
已经有4人回复
-
如何查询2026年获批项目名称？
已经有3人回复
-
两块石头
已经有10人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
初秋的晨风
已经有6人回复
-
Ei源刊怎么投
已经有4人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2016-01-22 20:36:15
qinzhen724
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1526.9
- 帖子: 135
- 在线: 46.4小时
- 虫号: 2706419
- 注册: 2013-10-08
- 性别: GG
- 专业: 工程热物理相关交叉领域
3楼2016-01-27 09:10:11
wang_chuan
银虫 (正式写手)
- 应助: 64 (初中生)
- 金币: 1786
- 散金: 1900
- 红花: 9
- 帖子: 512
- 在线: 90.7小时
- 虫号: 953632
- 注册: 2010-02-07
- 性别: GG
- 专业: 人工智能与知识工程
4楼2016-01-27 22:30:02
qinzhen724
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1526.9
- 帖子: 135
- 在线: 46.4小时
- 虫号: 2706419
- 注册: 2013-10-08
- 性别: GG
- 专业: 工程热物理相关交叉领域
5楼2016-02-01 23:25:11