troyzz

應(yīng)助: 0 (幼兒園)
金幣: 1733.3
帖子: 42
在線: 64.2小時(shí)
蟲(chóng)號(hào): 1525201
注冊(cè): 2011-12-05
性別: GG
專業(yè): 玻璃材料

[求助] Matlab二次規(guī)劃語(yǔ)言，高分求大神幫忙修改

在網(wǎng)上看了一遍關(guān)于配料比優(yōu)化的論文，完全按照論文內(nèi)容試著編了一下，可是無(wú)法運(yùn)行，求大神幫忙改一下。論文就不附上了，有涉嫌侵權(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ù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 程序語(yǔ)言.docx

2013-05-17 16:43:10, 31.2 K

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

給小老板認(rèn)識(shí)的老外寫(xiě)郵件，請(qǐng)求幫忙修改投稿語(yǔ)言，求個(gè)模板，多謝多謝已經(jīng)有7人回復(fù)
幫忙用Mathematica求解一個(gè)微分方程組已經(jīng)有6人回復(fù)
幫忙看一下這條語(yǔ)句是什么意思已經(jīng)有7人回復(fù)
linux終端運(yùn)行程序出錯(cuò)，求大神幫忙已經(jīng)有6人回復(fù)
求高手幫忙修改一段c++程序已經(jīng)有4人回復(fù)
初學(xué)matlab 求解方程大家?guī)蛶兔Π�。�！?/a> 已經(jīng)有11人回復(fù)

Perl 生物信息學(xué)編程求助，運(yùn)行出現(xiàn)問(wèn)題，煩請(qǐng)幫忙查找一下錯(cuò)誤，謝謝！已經(jīng)有3人回復(fù)
求高人幫忙看一下我這個(gè)Fortran程序的錯(cuò)誤出在哪里？已經(jīng)有12人回復(fù)
請(qǐng)幫忙用matlab求矩陣特征值已經(jīng)有8人回復(fù)
幫忙解釋一個(gè)read語(yǔ)句，謝謝已經(jīng)有4人回復(fù)
求幫忙啊~ matlab一多項(xiàng)式結(jié)果太長(zhǎng)，怎么寫(xiě)到txt里？已經(jīng)有14人回復(fù)
有會(huì)看fortran程序的，能幫忙看看這個(gè)程序嗎？已經(jīng)有10人回復(fù)
內(nèi)循環(huán)問(wèn)題，求大蝦幫忙已經(jīng)有3人回復(fù)
JCA審稿意見(jiàn)返回，編輯要求很嚴(yán)厲，語(yǔ)言的問(wèn)題，想找一個(gè)公司幫忙改，大家推薦下已經(jīng)有12人回復(fù)
求幫忙，關(guān)于含有bessel函數(shù)的方程的編程求解問(wèn)題已經(jīng)有11人回復(fù)
寫(xiě)了一個(gè)fortran90的小程序，編譯通不過(guò)，請(qǐng)大俠幫忙已經(jīng)有59人回復(fù)
【求助】50金幣求大家?guī)兔Ω囊欢纬绦?/a> 已經(jīng)有13人回復(fù)

【求助】幫忙看看哪錯(cuò)了error: invalid lvalue in assignment 已經(jīng)有6人回復(fù)
【求助】f90子程序怎么編譯錯(cuò)誤？請(qǐng)大家?guī)兔�，謝謝已經(jīng)有3人回復(fù)
【求助】用fortran怎么描述下面這樣的情況，我是新手，請(qǐng)大家?guī)兔?/a> 已經(jīng)有8人回復(fù)

【求助】幫忙看一個(gè)提示，子程序還是undefined reference to 已經(jīng)有4人回復(fù)
【求助】幫忙改寫(xiě)展寬程序已經(jīng)有11人回復(fù)
【求助】幫忙看一下，問(wèn)題程序如何修改已經(jīng)有3人回復(fù)

Fighting~~

1樓 2013-05-17 16:43:29

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ausir

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 18 (小學(xué)生)
金幣: 2009.2
紅花: 4
帖子: 500
在線: 178.9小時(shí)
蟲(chóng)號(hào): 1939287
注冊(cè): 2012-08-15
性別: GG
專業(yè): 集成電路制造與封裝

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

這里面有幾個(gè)問(wèn)題：

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

贊一下

回復(fù)此樓

2樓2013-05-18 13:35:58

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ausir

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 18 (小學(xué)生)
金幣: 2009.2
紅花: 4
帖子: 500
在線: 178.9小時(shí)
蟲(chóng)號(hào): 1939287
注冊(cè): 2012-08-15
性別: GG
專業(yè): 集成電路制造與封裝

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★
troyzz: 金幣+8, ★有幫助, 謝謝你的解答 2013-05-19 11:53:01
xzhdty: 謝謝參與 2013-05-19 12:00:29

引用回帖:

2樓: Originally posted by ausir at 2013-05-18 13:35:58
這里面有幾個(gè)問(wèn)題：

1. 把 myfun 這個(gè)用戶自定義函數(shù)單獨(dú)存在一個(gè).m文件中。后面在 fgoalattain中調(diào)用；
2. ub 的最后一項(xiàng)比 lb 對(duì)應(yīng)的數(shù)值小，ub中是12， lb中是315；
3. ub 矩陣中的寫(xiě)法，1000 寫(xiě)作 1 000， ...

補(bǔ)充說(shuō)明一下：我在計(jì)算時(shí)，把Ub的最后一項(xiàng)設(shè)置為1200，這樣就大于lb中的315，能夠進(jìn)行運(yùn)算。但是無(wú)法保證是否是真實(shí)的需要，請(qǐng)樓主斟酌。

贊一下

回復(fù)此樓

3樓2013-05-18 13:38:06

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

troyzz

木蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1733.3
帖子: 42
在線: 64.2小時(shí)
蟲(chóng)號(hào): 1525201
注冊(cè): 2011-12-05
性別: GG
專業(yè): 玻璃材料

引用回帖:

3樓: Originally posted by ausir at 2013-05-18 13:38:06
補(bǔ)充說(shuō)明一下：我在計(jì)算時(shí)，把Ub的最后一項(xiàng)設(shè)置為1200，這樣就大于lb中的315，能夠進(jìn)行運(yùn)算。但是無(wú)法保證是否是真實(shí)的需要，請(qǐng)樓主斟酌。...

原來(lái)是復(fù)制過(guò)程中出現(xiàn)很多錯(cuò)誤，我又仔細(xì)的看了一下那篇論文，并參考了您的意見(jiàn)，自己演算了一遍，附上修改后的語(yǔ)言還是不行，其實(shí)我對(duì)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ù)此樓

Fighting~~

4樓2013-05-18 15:01:59

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

troyzz

木蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1733.3
帖子: 42
在線: 64.2小時(shí)
蟲(chóng)號(hào): 1525201
注冊(cè): 2011-12-05
性別: GG
專業(yè): 玻璃材料

引用回帖:

我已經(jīng)解決了，把程序分成了兩個(gè)，一、二行保存為myfun.m，組成嵌入式，謝謝你的解答。

贊一下

回復(fù)此樓

Fighting~~

5樓2013-05-19 11:52:02

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3	一瓶苯甲酸 2026-03-14	3/150	2026-03-19 15:17 by 盡舜堯1
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +11	yangfz 2026-03-17	11/550	2026-03-19 15:06 by houyaoxu
[考研] 化學(xué)求調(diào)劑 +3	臨澤境llllll 2026-03-17	4/200	2026-03-19 13:59 by houyaoxu
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲(chóng)200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3	Ncdx123456 2026-03-19	3/150	2026-03-19 13:18 by houyaoxu
[考研] 346求調(diào)劑[0856] +3	WayneLim327 2026-03-16	6/300	2026-03-19 11:21 by WayneLim327
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無(wú)際的草原
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無(wú)際的草原
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu