24小时热门版块排行榜    

查看: 785  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

changpingxu

木虫 (小有名气)

[求助] Matlab非线性规划编程求助

以下代码有错,但不知错在哪,望高手指出。
函数表达式请见附件。

A=[];B=[];
W0=rand(23,1);
LB=zeros(23,1);UB=ones(23,1);%下限都为0上限为1
Aeq=ones(1,23);Beq=[1];%线性等式矩阵A为系数B为右端常数
[w,fval]=fmincon('FUN1',W0,A,B,Aeq,Beq,LB,UB)

Ci=[ 0.4404 0 0.1180 0 0 1.0000 1.0000 0 1.0000 0.5476 1.0000 1.0000 0 1.0000 1.0000 1.0000 1.0000 0.7048 1.0000 1.0000 1.0000 1.0000 0;0 0.0062 0 0.4370    0.9310 0.1845 0 0.8640 0.4438 0.5476 1.0000 1.0000 1.0000 0 0 0.5476 0.5476 0.7048 0.5476 0 0.5476 0 0.8407;0.8450 0.7967 0.7998 0.2257 1.0000 0 0.6592    1.0000 0.4438 0 0.2945 0 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000;1.0000 1.0000 1.0000 1.0000 0.6049 0 0.6592 0.8640 0 1.0000 0 0.5476 0 0 0 0.5476 0.5476 0    0.5476 0 0.5476 0 0.4603];
function f=FUN1(w);
f=sqrt((Ci(1,1)*w1-w1)^2+(Ci(1,2)*w2-w2)^2+(Ci(1,3)*w3-w3)^2+(Ci(1,4)*w4)^2+(Ci(1,5)*w5)^2+(Ci(1,6)*w6)^2+(Ci(1,7)*w7-w7)^2+(Ci(1,8)*w8)^2+(Ci(1,9)*w9)^2+(Ci(1,10)*w10)^2+(Ci(1,11)*w11)^2+(Ci(1,12)*w12)^2+(Ci(1,13)*w13-w13)^2+(Ci(1,14)*w14-w14)^2+(Ci(1,15)*w15-w15)^2+(Ci(1,16)*w16-w16)^2+(Ci(1,17)*w17-w17)^2+(Ci(1,18)*w18-w18)^2+(Ci(1,19)*w19-w19)^2+(Ci(1,20)*w20-w20)^2+(Ci(1,21)*w21-w21)^2+(Ci(1,22)*w22-w22)^2+(Ci(1,23)*w23-w23)^2)-sqrt((Ci(1,1)*w1)^2+(Ci(1,2)*w2)^2+(Ci(1,3)*w3)^2+(Ci(1,4)*w4-w4)^2+(Ci(1,5)*w5-w5)^2+(Ci(1,6)*w6-w6)^2+(Ci(1,7)*w7)^2+(Ci(1,8)*w8-w8)^2+(Ci(1,9)*w9-w9)^2+(Ci(1,10)*w10-w10)^2+(Ci(1,11)*w11-w11)^2+(Ci(1,12)*w12-w12)^2+(Ci(1,13)*w13)^2+(Ci(1,14)*w14)^2+(Ci(1,15)*w15)^2+(Ci(1,16)*w16)^2+(Ci(1,17)*w17)^2+(Ci(1,18)*w18)^2+(Ci(1,19)*w19)^2+(Ci(1,20)*w20)^2+(Ci(1,21)*w21)^2+(Ci(1,22)*w22)^2+(Ci(1,23)*w23)^2)+sqrt((Ci(2,1)*w1-w1)^2+(Ci(2,2)*w2-w2)^2+(Ci(2,3)*w3-w3)^2+(Ci(2,4)*w4)^2+(Ci(2,5)*w5)^2+(Ci(2,6)*w6)^2+(Ci(2,7)*w7-w7)^2+(Ci(2,8)*w8)^2+(Ci(2,9)*w9)^2+(Ci(2,10)*w10)^2+(Ci(2,11)*w11)^2+(Ci(2,12)*w12)^2+(Ci(2,13)*w13-w13)^2+(Ci(2,14)*w14-w14)^2+(Ci(2,15)*w15-w15)^2+(Ci(2,16)*w16-w16)^2+(Ci(2,17)*w17-w17)^2+(Ci(2,18)*w18-w18)^2+(Ci(2,19)*w19-w19)^2+(Ci(2,20)*w20-w20)^2+(Ci(2,21)*w21-w21)^2+(Ci(2,22)*w22-w22)^2+(Ci(2,23)*w23-w23)^2)-sqrt((Ci(2,1)*w1)^2+(Ci(1,2)*w2)^2+(Ci(2,3)*w3)^2+(Ci(2,4)*w4-w4)^2+(Ci(2,5)*w5-w5)^2+(Ci(2,6)*w6-w6)^2+(Ci(2,7)*w7)^2+(Ci(2,8)*w8-w8)^2+(Ci(2,9)*w9-w9)^2+(Ci(2,10)*w10-w10)^2+(Ci(1,11)*w11-w11)^2+(Ci(2,12)*w12-w12)^2+(Ci(2,13)*w13)^2+(Ci(2,14)*w14)^2+(Ci(2,15)*w15)^2+(Ci(2,16)*w16)^2+(Ci(2,17)*w17)^2+(Ci(2,18)*w18)^2+(Ci(2,19)*w19)^2+(Ci(2,20)*w20)^2+(Ci(2,21)*w21)^2+(Ci(2,22)*w22)^2+(Ci(2,23)*w23)^2)+sqrt((Ci(3,1)*w1-w1)^2+(Ci(3,2)*w2-w2)^2+(Ci(3,3)*w3-w3)^2+(Ci(3,4)*w4)^2+(Ci(3,5)*w5)^2+(Ci(3,6)*w6)^2+(Ci(3,7)*w7-w7)^2+(Ci(3,8)*w8)^2+(Ci(3,9)*w9)^2+(Ci(3,10)*w10)^2+(Ci(3,11)*w11)^2+(Ci(3,12)*w12)^2+(Ci(3,13)*w13-w13)^2+(Ci(3,14)*w14-w14)^2+(Ci(3,15)*w15-w15)^2+(Ci(3,16)*w16-w16)^2+(Ci(3,17)*w17-w17)^2+(Ci(3,18)*w18-w18)^2+(Ci(3,19)*w19-w19)^2+(Ci(3,20)*w20-w20)^2+(Ci(3,21)*w21-w21)^2+(Ci(3,22)*w22-w22)^2+(Ci(3,23)*w23-w23)^2)-sqrt((Ci(3,1)*w1)^2+(Ci(1,2)*w2)^2+(Ci(3,3)*w3)^2+(Ci(3,4)*w4-w4)^2+(Ci(3,5)*w5-w5)^2+(Ci(3,6)*w6-w6)^2+(Ci(3,7)*w7)^2+(Ci(3,8)*w8-w8)^2+(Ci(3,9)*w9-w9)^2+(Ci(3,10)*w10-w10)^2+(Ci(1,11)*w11-w11)^2+(Ci(3,12)*w12-w12)^2+(Ci(3,13)*w13)^2+(Ci(3,14)*w14)^2+(Ci(3,15)*w15)^2+(Ci(3,16)*w16)^2+(Ci(3,17)*w17)^2+(Ci(3,18)*w18)^2+(Ci(3,19)*w19)^2+(Ci(3,20)*w20)^2+(Ci(3,21)*w21)^2+(Ci(3,22)*w22)^2+(Ci(3,23)*w23)^2)+sqrt((Ci(4,1)*w1-w1)^2+(Ci(4,2)*w2-w2)^2+(Ci(4,3)*w3-w3)^2+(Ci(4,4)*w4)^2+(Ci(4,5)*w5)^2+(Ci(4,6)*w6)^2+(Ci(4,7)*w7-w7)^2+(Ci(4,8)*w8)^2+(Ci(4,9)*w9)^2+(Ci(4,10)*w10)^2+(Ci(4,11)*w11)^2+(Ci(4,12)*w12)^2+(Ci(4,13)*w13-w13)^2+(Ci(4,14)*w14-w14)^2+(Ci(4,15)*w15-w15)^2+(Ci(4,16)*w16-w16)^2+(Ci(4,17)*w17-w17)^2+(Ci(4,18)*w18-w18)^2+(Ci(4,19)*w19-w19)^2+(Ci(4,20)*w20-w20)^2+(Ci(4,21)*w21-w21)^2+(Ci(4,22)*w22-w22)^2+(Ci(4,23)*w23-w23)^2)-sqrt((Ci(4,1)*w1)^2+(Ci(1,2)*w2)^2+(Ci(4,3)*w3)^2+(Ci(4,4)*w4-w4)^2+(Ci(4,5)*w5-w5)^2+(Ci(4,6)*w6-w6)^2+(Ci(4,7)*w7)^2+(Ci(4,8)*w8-w8)^2+(Ci(4,9)*w9-w9)^2+(Ci(4,10)*w10-w10)^2+(Ci(1,11)*w11-w11)^2+(Ci(4,12)*w12-w12)^2+(Ci(4,13)*w13)^2+(Ci(4,14)*w14)^2+(Ci(4,15)*w15)^2+(Ci(4,16)*w16)^2+(Ci(4,17)*w17)^2+(Ci(4,18)*w18)^2+(Ci(4,19)*w19)^2+(Ci(4,20)*w20)^2+(Ci(4,21)*w21)^2+(Ci(4,22)*w22)^2+(Ci(4,23)*w23)^2);%定义函数

Matlab非线性规划编程求助
目标函数表达式.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

中国地质

新虫 (初入文坛)

【答案】应助回帖

我也不知道问题出在哪了
3楼2014-05-31 10:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

sdkjdx_lixue

新虫 (初入文坛)

【答案】应助回帖

怎么没有人回帖啊,求求 各位大神帮忙啊!
2楼2014-05-31 10:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:39 by G6APbkg8SA6w
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:38 by G6APbkg8SA6w
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 09:17 by G6APbkg8SA6w
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见