24小时热门版块排行榜    

查看: 3809  |  回复: 1

laoyeye2

木虫 (小有名气)

[求助] matlab中用ga 多目标优化求解最小值

求助:
%目标函数
function f = my_first_multi2(x)
f=28.44-4.02583*x(1)-1.69425*x(2)+1.61083*x(3)+3.73125*x(4)+0.040375*x(1)*x(2)+0.055625*x(1)*x(3)-0.10937*x(1)*x(4)-0.008375*x(2)*x(3)-0.050375*x(2)*x(4)-0.10062*x(3)*x(4)+0.21031*x(1)^2+0.026363*x(2)^2 +0.0028125*x(3)^2+0.026563*x(4)^2;

%约束函数
function [c,ceq]=ga_con1(x)
c=[-2.11833+0.955422*x(1)+0.105*x(2)-0.81042*x(3)+0.3075*x(4)+0.080625*x(1)*x(2)-0.51688*x(1)*x(3)+0.28563*x(1)*x(4)-0.15338*x(2)*x(3)+0.10562*x(2)*x(4)-0.58938*x(3)*x(4)-0.014896*x(1)^2-0.000254*x(2)^2 +0.67885*x(3)^2+0.1476*x(4)^2-5];
ceq=[];

求解代码:
FitnessFunction=@ my_first_multi2;
lb=[4 15 4 5];
ub=[8 35 8 9];
[x,fval] = ga(FitnessFunction,4,[],[],[],[],lb,ub,@ga_con1)

最后出现的问题:Optimization terminated: no feasible point found.

x =

     4    15     4     5


大神帮我看看,问题出在哪里了,万分感谢。
回复此楼

» 猜你喜欢

世上无难事。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laoyeye2

木虫 (小有名气)

首先非常感谢小木虫里面的“文刀星 ”,他的回复很好的解决了我的问题。以下为回复:
遇到这种问题的时候,第一,看看自己的约束函数有没有写错,如果是约束错了,当然就找不到满足约束的点啦。你随便把一个x=[4 15 4 5];代入你的约束函数计算,会发现,本来你的约束只有一个,但是返回的却是两个约束值!!!当然就是约束有问题了。问题的产生在于你写约束函数的格式或者说习惯很不好造成的。c=[...]的写法很不好,只要里面符号之间多了一个空格,那就多了一个约束,而且不是你想要的约束,你可以查查看,在0.000254*x(2)^2 +0.67885*x(3)^2的+号前多了一个空格,使得你的约束由一变二。最好这样写:function [c,ceq]=ga_con1(x)c=zeros(1,1);%多个约束的时候c=zeros(N,1),然后c(1)=...;c(2)=...;等等c(1)=-2.11833+0.955422*x(1)+0.105*x(2)-... 0.81042*x(3)+0.3075*x(4)+0.080625*x(1)*x(2)-0.51688*x(1)*x(3)+... 0.28563*x(1)*x(4)-0.15338*x(2)*x(3)+0.10562*x(2)*x(4)-... 0.58938*x(3)*x(4)-0.014896*x(1)^2-0.000254*x(2)^2+... 0.67885*x(3)^2+0.1476*x(4)^2-5;ceq=[];式子太长的时候,可以在符号的后面加三个点,然后回车在下面写,这样好看些,便于检查。记得是在符号的后面加,在前面加的话,某些情况也会出现由一生二的情况!!
世上无难事。
2楼2015-12-02 13:36:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 laoyeye2 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 0703化学调剂,求各位老师收留 +9 秋有木北 2026-03-14 9/450 2026-03-18 19:49 by macy2011
[教师之家] 焦虑 +8 水冰月月野兔 2026-03-13 12/600 2026-03-18 15:27 by 咪呜喵呜
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +7 rare12345 2026-03-18 7/350 2026-03-18 14:31 by laoshidan
[考研] 085601材料工程专硕求调剂 +6 慕寒mio 2026-03-16 6/300 2026-03-18 14:26 by 007_lilei
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
信息提示
请填处理意见