24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

199107011523

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 234.8
帖子: 93
在线: 31.9小时
虫号: 3441617
注册: 2014-09-25
专业: 工程热物理相关交叉领域

[求助] 优化分析已有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

回复此楼

» 猜你喜欢

药学硕士，第一、第二作者已发表6 篇 SCI，药理方向及相关方向2026年/2027年博士申请已经有4人回复
26年博士申请自荐-电催化已经有3人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有6人回复
收到国自然专家邀请后几年才会有本子送过来评已经有4人回复
考博已经有5人回复
26年申博自荐-计算机视觉已经有4人回复
药化及相关博士的申请已经有3人回复
一篇MDPI论文改变了学习工作和生活已经有4人回复

1楼 2016-06-21 18:15:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

199107011523

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 234.8
帖子: 93
在线: 31.9小时
虫号: 3441617
注册: 2014-09-25
专业: 工程热物理相关交叉领域

引用回帖:

5楼: Originally posted by fungarwai at 2016-06-22 18:45:12
我看到的是这样

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 info ...

恩恩，就是这样子的，可是，t的结果为什么和约束条件这么大相径庭呢