| 1 | 1/1 | 返回列表 |
| 查看: 928 | 回復: 0 | |||
Justiceall新蟲 (初入文壇)
|
[求助]
求助 編出來的matlab規(guī)劃問題程序運行不了(線性約束方面的)
|
|
代碼如下: 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. >> 請各位大哥幫忙看看,我都改了一天了還是運行不出來啊。。。。。急死。。。就這些金幣了全給你們了。。。。。謝謝 |
找到一些相關的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277跪求調劑 +5 | 1915668 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 322求調劑 +3 | 宋明欣 2026-03-27 | 3/150 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 化學調劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 352分 化工與材料 +3 | 海納百川Ly 2026-03-27 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調劑 +7 | 邱gl 2026-03-27 | 7/350 |
|
|
[考研] 308求調劑 +6 | 墨墨漠 2026-03-27 | 6/300 |
|
|
[考研] 322求調劑 +4 | 我真的很想學習 2026-03-23 | 4/200 |
|
|
[考研] 一志愿鄭大085600,310分求調劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[碩博家園] 北京林業(yè)大學碩導招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[碩博家園] 招收生物學/細胞生物學調劑 +3 | IceGuo 2026-03-26 | 4/200 |
|
|
[考研] 調劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 材料考研求調劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
[考研] 271求調劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
[考研] 334分 一志愿武理 材料求調劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 263求調劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 293求調劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 308求調劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 求調劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 一志愿北京化工大學070300 學碩336求調劑 +5 | vv迷 2026-03-21 | 8/400 |
|