| 5 | 1/1 | 返回列表 |
| 查看: 557 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
troyzz木蟲 (初入文壇)
|
[求助]
Matlab二次規(guī)劃語言,高分求大神幫忙修改
|
||
|
在網(wǎng)上看了一遍關(guān)于配料比優(yōu)化的論文,完全按照論文內(nèi)容試著編了一下,可是無法運行,求大神幫忙改一下。論文就不附上了,有涉嫌侵權(quán)嫌疑,把從里面摘抄的一部分附上。 1 function f=myfun(X) 2 f=0.35*X(1)+0.002*X(2)+0.02*X(3)+0.025*X(4) +0.035*X(5)+6*X(6); 3 goal=170; 4 weight=170; 5 X0=[400;170;100;550;1200;6]; 6 b=[0;0;0;0;0;0;0;-500;600]; 7 Aeq=[0.3226 1 0.4808 0.3774 0.3704 0]; 8 beq=990; 9 lb=[300;130;50;500;1000;315]; 10 ub=[500;210;200;1 000;1400;12]; 11 A=[0.2 -1 0.2 0 0 0; -0.4 1 -0.4 0 0 0; 0 0 0 -0.65 0.35 0; 0 0 0 0.6 01.4 0; -14.6 52.71 -14.6 0 0 0; 0.01 0 0 0 0 -1; -0.02 0 0 0 0 1; -1 0 -1 0 0 0; 1 0 1 0 0 0]; 12 [X,fval,attainfactor,exitflag]=fgoalattain(xy,X0,goal,weight,A,b,Aeq,beq,lb,ub)[ Last edited by troyzz on 2013-5-18 at 12:07 ] |

木蟲 (初入文壇)
|
原來是復(fù)制過程中出現(xiàn)很多錯誤,我又仔細(xì)的看了一下那篇論文,并參考了您的意見,自己演算了一遍,附上修改后的語言還是不行,其實我對fgoalattain不熟悉。 function f=myfun(X) f=0.35*X(1)+0.002*X(2)+0.02*X(3)+0.025*X(4) +0.035*X(5)+6*X(6); goal=170; weight=170; X0=[400;170;100;550;1200;6]; b=[0;0;0;0;0;0;0;-500;600]; Aeq=[0.3226 1 0.4808 0.3774 0.3704 0]; beq=990; lb=[300;130;50;500;1000;3.5]; ub=[500;210;200;1000;1400;12]; A=[0.2 -1 0.2 0 0 0; -0.4 1 -0.4 0 0 0; 0 0 0 -0.65 0.35 0; 0 0 0 0.6 -0.4 0; -14.6 52.708 -14.6 0 0 0; 0.01 0 0 0 0 -1; -0.02 0 0 0 0 1; -1 0 -1 0 0 0; 1 0 1 0 0 0]; [X,fval,attainfactor,exitflag]=fgoalattain(@myfun,X0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon) |

木蟲 (正式寫手)
|
這里面有幾個問題: 1. 把 myfun 這個用戶自定義函數(shù)單獨存在一個.m文件中。后面在 fgoalattain中調(diào)用; 2. ub 的最后一項比 lb 對應(yīng)的數(shù)值小,ub中是12, lb中是315; 3. ub 矩陣中的寫法,1000 寫作 1 000,這個導(dǎo)致識別錯誤; 4. 矩陣A的第4行第5列的數(shù)值到底是多少? 01.4? 可以按照1.4來計算,但是計算迭代到700次時依然無法收斂,系統(tǒng)自動停止運算,沒有結(jié)果。 5. 最后的關(guān)鍵,第12行的函數(shù)中,根本就沒有提到你自己定義的函數(shù) myfun,而是使用了沒有任何定義的 xy。根本不會有任何運算。 6. 關(guān)于fgoalattain函數(shù)的調(diào)用語法問題,在myfun 前面要加一個 @符號或使用 'myfun' 。運算結(jié)果一樣。都無法取得附件中中所要求的結(jié)果。 以上看法供參考。 |
木蟲 (正式寫手)
木蟲 (初入文壇)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +8 | S.H1 2026-03-18 | 8/400 |
|
|
[考研] 290求調(diào)劑 +4 | ^O^乜 2026-03-19 | 4/200 |
|
|
[考研] 08工學(xué)調(diào)劑 +4 | 用戶573181 2026-03-20 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|