| 查看: 722 | 回复: 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 |
回复此楼» 猜你喜欢
求个博导看看
已经有8人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有7人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
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
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