| 1 | 1/1 | 返回列表 |
| 查看: 927 | 回復(fù): 0 | ||
Justiceall新蟲 (初入文壇)
|
[求助]
求助 編出來(lái)的matlab規(guī)劃問(wèn)題程序運(yùn)行不了(線性約束方面的)
|
|
代碼如下: 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的對(duì)角陣 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檢測(cè)出的錯(cuò)誤是: >> 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. >> 請(qǐng)各位大哥幫忙看看,我都改了一天了還是運(yùn)行不出來(lái)啊。。。。。急死。。。就這些金幣了全給你們了。。。。。謝謝 |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +3 | 辰沐5211314 2026-03-26 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-27 | 3/150 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑!。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|