| 5 | 1/1 | 返回列表 |
| 查看: 556 | 回復: 4 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
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 ] |

木蟲 (正式寫手)
木蟲 (正式寫手)
|
這里面有幾個問題: 1. 把 myfun 這個用戶自定義函數(shù)單獨存在一個.m文件中。后面在 fgoalattain中調(diào)用; 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)自動停止運算,沒有結(jié)果。 5. 最后的關(guān)鍵,第12行的函數(shù)中,根本就沒有提到你自己定義的函數(shù) myfun,而是使用了沒有任何定義的 xy。根本不會有任何運算。 6. 關(guān)于fgoalattain函數(shù)的調(diào)用語法問題,在myfun 前面要加一個 @符號或使用 'myfun' 。運算結(jié)果一樣。都無法取得附件中中所要求的結(jié)果。 以上看法供參考。 |
木蟲 (初入文壇)
|
原來是復制過程中出現(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ā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工 求調(diào)劑 +15 | enenenhui 2026-03-13 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 3/150 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5 | Doleres 2026-03-19 | 5/250 |
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 化學工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 考研化學學碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 0703 物理化學調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 080500,材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|