24小时热门版块排行榜    

查看: 969  |  回复: 0

Justiceall

新虫 (初入文坛)

[求助] 求助 编出来的matlab规划问题程序运行不了(线性约束方面的)

代码如下:
function profit
[x,fval]=ga(@fitness,5,[],[],[],[],0,1,@nonl)
function y=fitness(x)
y=(-(7.42-(1.56*x(1)+258.98*x(2)+1.26*x(3)+3.32*x(4)+5.00*x(5))));
function [c,ce,A,B,X,C,D,E,F,O,S,G,H,I,J,L,e,r,m]=nonl(x)
A=[18,0,0,0,0;0,0,0,0,0;0,0,8,0,0;0,0,0,21,0;0,0,0,0,9]; %s=R-M的对角阵
B=[8,0,8,21,9];
s=B';
X=[x(1),x(2),x(3),x(4),x(5)];
C=[118,0,95,111,78];
r=C';
D=[100,0,87,90,69];
m=D';
E=[0,0,0,0,0];  %组分油烯烃含量
O=E';
F=[0,0,0,0,0];  %组分油烯烃含量平方
S=F';
G=[0,0,0,0,0];  %组分油芳烃含量
H=G';
I=[0,0,0,0,0];  %组分油芳烃含量平方
J=I';
L=[1,1,1,1,1];
e=L';
R=r*X+0.03224*(r*A*X-((r*X)*(s*X))/(e*X))+0.00101*((S*X)-((O*X).^2)/(e*X));
M=m*X+0.0445*(m*A*X-((m*X)*(s*X))/(e*X))+0.00081*((S*X)-((O*X).^2)/(e*X))+((-0.00645)/10000*(e*X))/(((J*X)-((H*X).^2)/(e*X)).^2);
c(1)=-(r*X+0.03224*(r*A*X-((r*X)*(s*X))/(e*X))+0.00101*((S*X)-((O*X).^2)/(e*X)))+90;
c(2)=(-(R+M)/2)+85;
c(3)=((x(1)*(50.5).^1.25)+(x(2)*(0.00667).^1.25)+(x(3)*(39).^1.25)+(x(4)*(5.87).^1.25)+(x(5)*(75).^1.25)).^0.8-72;
ce=[];
matlab检测出的错误是:
>> profit
Error using  *
Inner matrix dimensions must agree.

Error in profit>nonl (line 24)
R=r*X+0.03224*(r*A*X-((r*X)*(s*X))/(e*X))+0.00101*((S*X)-((O*X).^2)/(e*X));

Error in validate>@(x)nonlcon(x,NonconFcnArgs{:}) (line 142)
    nonlcon = @(x) nonlcon(x,NonconFcnArgs{:});

Error in constrValidate (line 20)
        [cineq,ceq] = nonlcon(Iterate.x');

Error in gacommon (line 132)
[LinearConstr, Iterate,nineqcstr,neqcstr,ncstr] = constrValidate(NonconFcn, ...

Error in ga (line 322)
[x,fval,exitFlag,output,population,scores,FitnessFcn,nvars,Aineq,bineq,Aeq,beq,lb,ub, ...

Error in profit (line 2)
[x,fval]=ga(@fitness,5,[],[],[],[],[],[],@nonl)

Caused by:
    Failure in initial user-supplied nonlinear constraint function evaluation.

>>
请各位大哥帮忙看看,我都改了一天了还是运行不出来啊。。。。。急死。。。就这些金币了全给你们了。。。。。谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 Justiceall 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 10/500 2026-11-02 06:16 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 05:41 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:37 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-02 05:34 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 03:13 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 24/1200 2026-11-02 02:39 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
信息提示
请填处理意见