24小时热门版块排行榜    

查看: 946  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 函评 +3 布布和一二 2026-06-13 3/150 2026-06-13 23:35 by 芋圆香飘飘
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] T^T +5 q1zg6x4lf3 2026-06-12 5/250 2026-06-13 18:53 by lkmi35jdzt
[硕博家园] O(∩_∩)O +5 jk31rb6vyz 2026-06-12 7/350 2026-06-13 18:45 by lkmi35jdzt
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[考博] T^T +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 14:45 by 80j16ujcu8
[考研] ╭(╯ε╰)╮ +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 14:22 by 80j16ujcu8
[考博] T^T +4 jk31rb6vyz 2026-06-12 4/200 2026-06-13 12:09 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[考博] →_→ +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 08:49 by t4hz8376wy
[找工作] *囍 +3 m584r2z4j3 2026-06-12 3/150 2026-06-13 08:29 by t4hz8376wy
[论文投稿] →_→ +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 08:13 by t4hz8376wy
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[教师之家] 重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶 +6 yexuqing 2026-06-07 7/350 2026-06-12 16:46 by wwwkkk83
[基金申请] 前几天时间戳更新了 +8 vito刘 2026-06-11 13/650 2026-06-12 14:17 by BlakeReary
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
请填处理意见