| 查看: 822 | 回复: 7 | ||
199107011523新虫 (小有名气)
|
[求助]
优化分析 已有1人参与
|
|
已知t1=67;t2=110;t3=82;t4=73;t5=33;t6=0.5;t7=0.54;69≤t4≤75;31≤t5≤35;3%≤t7-t6≤5%； y1=2t5-5t5^2-3/t5^2+(a-1)t5-at3 y2=2t5-5t5^2-3/t5^2+3t4 y3=3t4+5t5^2-3/t5^3+(a-1)t-at1 z1=y1/(t1+3*t2); z2=y2/(3*t3+2*t2); z3=y3/(t1+t5); a=t7/(t7-t6); 10≤a≤50; Z=z1+z2+z3; F=(y1+y2)/y3; 求min f=Z/F maxF minZ 以min f=Z/F为例，编写m.文件如下，运行有警告提示，希望得到高手指点： function f=fun567(t) a=t(7)/(t(7)-t(6)); y1=2*t(5)-5*t(5)^2-3/t(5)^2+(a-1)*t(2)-a*t(3); y2=2*t(5)-5*t(5)^2-3/t(5)^2+3*t(4); y3=3*t(4)+5*t(5)^2-3/t(5)^3+(a-1)*t(2)-a*t(1); z1=y1/(t(1)+3*t(2)); z2=y2/(3*t(3)+2*t(2)); z3=y3/(t(1)+t(5)); Z=z1+z2+z3; F=(y1+y2)/y3; f=Z/F function [c,ceq]=mycon3_1(t) c(1)=t(4)-8; c(2)=5-t(4); c(3)=t(5)-7; c(4)=3-t(5); c(5)=t(7)/(t(7)-t(6))-50; c(6)=10-t(6)/(t(7)-t(6)); c(7)=t(7)-t(6)-0.06; c(8)=0.03-t(7)+t(6); ceq=[]; 调用函数如下： t=[67;110;82;73;33;0.5;0.54]; [t,fval]=fmincon(@(t) fun567(t),t,[],[],Aeq,beq,[],[],@mycon3_1)@laosam280 |
回复此楼» 猜你喜欢
如何实现卤原子转化
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
考博自荐
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
-
申请博士
已经有3人回复
199107011523
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 234.8
- 帖子: 93
- 在线: 31.9小时
- 虫号: 3441617
- 注册: 2014-09-25
- 专业: 工程热物理相关交叉领域
2楼2016-06-22 10:18:08
3楼2016-06-22 17:57:54
199107011523
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 234.8
- 帖子: 93
- 在线: 31.9小时
- 虫号: 3441617
- 注册: 2014-09-25
- 专业: 工程热物理相关交叉领域
4楼2016-06-22 18:32:19
|
我看到的是这样 Warning: The default trust-region-reflective algorithm does not solve problems with the constraints you have specified. FMINCON will use the active-set algorithm instead. For information on applicable algorithms, see Choosing the Algorithm in the documentation. > In fmincon at 500 In test at 4 Warning: Your current settings will run a different algorithm (interior-point) in a future release. > In fmincon at 505 In test at 4 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> Active inequalities (to within options.TolCon = 1e-06): lower upper ineqlin ineqnonlin 2 5 t = 1.0e+08 * -0.8936 -0.2746 1.2088 0.0000 0.0000 0.0000 0.0000 fval = -3.4171 你的c里面那些值怎么取得跟题目不一样 c(1)=t(4)-75; c(2)=69-t(4); c(3)=t(5)-35; c(4)=31-t(5); c(7)=t(7)-t(6)-0.05; |
5楼2016-06-22 18:45:12
199107011523
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 234.8
- 帖子: 93
- 在线: 31.9小时
- 虫号: 3441617
- 注册: 2014-09-25
- 专业: 工程热物理相关交叉领域
6楼2016-06-22 19:40:41
|
这些值改过来了吗？ c(1)=t(4)-75; c(2)=69-t(4); c(3)=t(5)-35; c(4)=31-t(5); c(7)=t(7)-t(6)-0.05; Warning: The default trust-region-reflective algorithm does not solve problems with the constraints you have specified. FMINCON will use the active-set algorithm instead. For information on applicable algorithms, see Choosing the Algorithm in the documentation. > In fmincon at 500 In test at 4 Warning: Your current settings will run a different algorithm (interior-point) in a future release. > In fmincon at 505 In test at 4 Solver stopped prematurely. fmincon stopped because it exceeded the function evaluation limit, options.MaxFunEvals = 700 (the default value). t = -32.1395 183.1122 70.3997 74.3837 32.1395 2.1947 2.2426 fval = -2.1166e+10 |
7楼2016-06-22 20:18:00
199107011523
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 234.8
- 帖子: 93
- 在线: 31.9小时
- 虫号: 3441617
- 注册: 2014-09-25
- 专业: 工程热物理相关交叉领域
|
是的…就是这个过程中t1,t2,t6,t7的变化太大…我原来设定约束条件Aeq=[1,0,0,0,0,0,0;0,1,0,0,0,0,0];beq=[67,110]做过分析，发现这个过程中，t6,t7的影响占据主要因素，变化和初值差距太大了…和预期想的结果不一样的 发自小木虫Android客户端 |
8楼2016-06-22 21:01:23
10