| 5 | 1/1 | 返回列表 |
| 查看: 555 | 回復: 4 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
troyzz木蟲 (初入文壇)
|
[求助]
Matlab二次規(guī)劃語言,高分求大神幫忙修改
|
|
|
在網(wǎng)上看了一遍關于配料比優(yōu)化的論文,完全按照論文內容試著編了一下,可是無法運行,求大神幫忙改一下。論文就不附上了,有涉嫌侵權嫌疑,把從里面摘抄的一部分附上。 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 ] |

木蟲 (初入文壇)

木蟲 (正式寫手)
|
這里面有幾個問題: 1. 把 myfun 這個用戶自定義函數(shù)單獨存在一個.m文件中。后面在 fgoalattain中調用; 2. ub 的最后一項比 lb 對應的數(shù)值小,ub中是12, lb中是315; 3. ub 矩陣中的寫法,1000 寫作 1 000,這個導致識別錯誤; 4. 矩陣A的第4行第5列的數(shù)值到底是多少? 01.4? 可以按照1.4來計算,但是計算迭代到700次時依然無法收斂,系統(tǒng)自動停止運算,沒有結果。 5. 最后的關鍵,第12行的函數(shù)中,根本就沒有提到你自己定義的函數(shù) myfun,而是使用了沒有任何定義的 xy。根本不會有任何運算。 6. 關于fgoalattain函數(shù)的調用語法問題,在myfun 前面要加一個 @符號或使用 'myfun' 。運算結果一樣。都無法取得附件中中所要求的結果。 以上看法供參考。 |
木蟲 (正式寫手)
木蟲 (初入文壇)
|
原來是復制過程中出現(xiàn)很多錯誤,我又仔細的看了一下那篇論文,并參考了您的意見,自己演算了一遍,附上修改后的語言還是不行,其實我對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) |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279分求調劑 一志愿211 +3 | chaojifeixia 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工 求調劑 +14 | enenenhui 2026-03-13 | 15/750 |
|
|
[考研] 288求調劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 招收調劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|
[考研] 求調劑,一志愿:南京航空航天大學大學 ,080500材料科學與工程學碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 0703化學調劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 0703化學 305求調劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 295求調劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 298-一志愿中國農業(yè)大學-求調劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
[考研] 收復試調劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 0703化學求調劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[基金申請] 被我言中:新模板不強調格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 材料專碩326求調劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 085601求調劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 0703 物理化學調劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 070305求調劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 080500,材料學碩302分求調劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 297求調劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|