24小时热门版块排行榜    

查看: 1616  |  回复: 3

laoyeye2

木虫 (小有名气)

[求助] matlab 多目标优化求解最小值 已有1人参与

matlab 多目标优化求解最小值
1.jpg


matlab 多目标优化求解最小值-1
2.jpg

以下为求解过程。。。。。。。。。。。。。。。。。。。。。。。。。
目标函数:
function f = my_first_multi2(x)
f(1)=121.06621+5.47522*x(1)-28.46988*x(2)-13.5757*x(3)-0.54759*x(1)*x(2)-0.38818*x(1)*x(3)+1.32485*x(2)*x(3) ...
    +0.14034*x(1)^2+1.87034*x(2)^2+0.52653*x(3)^2;
约束函数:
function [c,ceq]=ga_con1(x)
c=[161.23586-10.31686*x(1)-24.81371*x(2)-9.95421*x(3)+0.91842*x(1)*x(2)-0.098083*x(1)*x(3) ...
    +0.377*x(2)*x(3)+0.24078*x(1)^2+1.18925*x(2)^2+0.81175*x(3)^2-12;
    512.14149-54.38052*x(1)+59.73721*x(2)+70.37891*x(3)+5.09733*x(1)*x(2)-1.81717*x(1)*x(3)-6.75775*x(2)*x(3) ...
    1.50459*x(1)^2-7.99343*x(2)^2-1.05542*x(3)^2-530];
ceq=[];
求解:
FitnessFunction=@ my_first_multi2;
lb=[8 4 6];
ub=[14 6 8];
[x,fval] = ga(FitnessFunction,3,[],[],[],[],lb,ub,@ga_con1)



最后出现的问题:Optimization terminated: no feasible point found.
x =
  8     4     6
fval =
  23.0374
回复此楼
世上无难事。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

文刀星

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
laoyeye2: 金币+20, ★★★★★最佳答案, 谢谢。 2015-09-21 16:42:41
内容已删除
2楼2015-09-21 09:30:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

文刀星

木虫 (小有名气)

不知道为什么,按理说你写的约束方程中,得出的c应该是个列向量,可是实际运行的时候却是一个2*2的矩阵;矩阵中的元素分别是...前后的式子计算出来的,所以约束函数中的约束就已经不是你原来所想要的约束了。你可以试试看,我很好奇,为什么会发生这样的情况。
比如
x=[8,4,6];
c=[161.23586-10.31686*x(1)-24.81371*x(2)-9.95421*x(3)+0.91842*x(1)*x(2)-0.098083*x(1)*x(3) ...
    +0.377*x(2)*x(3)+0.24078*x(1)^2+1.18925*x(2)^2+0.81175*x(3)^2-12;
    512.14149-54.38052*x(1)+59.73721*x(2)+70.37891*x(3)+5.09733*x(1)*x(2)-1.81717*x(1)*x(3)-6.75775*x(2)*x(3) ...
    +1.50459*x(1)^2-7.99343*x(2)^2-1.05542*x(3)^2-530];
运行后得到的结果是
c =

  -55.5977   60.7089
  652.0240 -599.5962
但是我们想要的结果应该是
c =

  5.1113
  52.4278
这是什么情况?
3楼2015-09-21 09:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

文刀星

木虫 (小有名气)

啊!!!为什么把+号放到...之前就好了呢?
x=[8,4,6];
c=[161.23586-10.31686*x(1)-24.81371*x(2)-9.95421*x(3)+0.91842*x(1)*x(2)-0.098083*x(1)*x(3) +...
    0.377*x(2)*x(3)+0.24078*x(1)^2+1.18925*x(2)^2+0.81175*x(3)^2-12;
    512.14149-54.38052*x(1)+59.73721*x(2)+70.37891*x(3)+5.09733*x(1)*x(2)-1.81717*x(1)*x(3)-6.75775*x(2)*x(3)+ ...
    1.50459*x(1)^2-7.99343*x(2)^2-1.05542*x(3)^2-530];

c =

  5.1113
  52.4278
4楼2015-09-21 09:58:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 laoyeye2 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见