24小时热门版块排行榜    

查看: 787  |  回复: 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的回帖

sdkjdx_lixue

新虫 (初入文坛)

【答案】应助回帖

怎么没有人回帖啊,求求 各位大神帮忙啊!
2楼2014-05-31 10:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

中国地质

新虫 (初入文坛)

【答案】应助回帖

我也不知道问题出在哪了
3楼2014-05-31 10:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 changpingxu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 14:15 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+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 南海小哥
信息提示
请填处理意见