24小时热门版块排行榜

返回列表

changpingxu

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 3692.5
帖子: 85
在线: 119.8小时
虫号: 2523733
注册: 2013-06-27
性别: GG
专业: 工程管理

[求助] 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);%定义函数

目标函数表达式.jpg

回复此楼