troyzz

應(yīng)助: 0 (幼兒園)
金幣: 1733.3
帖子: 42
在線(xiàn): 64.2小時(shí)
蟲(chóng)號(hào): 1525201
注冊(cè): 2011-12-05
性別: GG
專(zhuān)業(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)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(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分人工智能專(zhuān)業(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)返回，編輯要求很?chē)?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
在線(xiàn): 178.9小時(shí)
蟲(chóng)號(hào): 1939287
注冊(cè): 2012-08-15
性別: GG
專(zhuān)業(yè): 集成電路制造與封裝

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

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

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

1. 把 myfun 這個(gè)用戶(hù)自定義函數(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
在線(xiàn): 178.9小時(shí)
蟲(chóng)號(hào): 1939287
注冊(cè): 2012-08-15
性別: GG
專(zhuān)業(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è)用戶(hù)自定義函數(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
在線(xiàn): 64.2小時(shí)
蟲(chóng)號(hào): 1525201
注冊(cè): 2011-12-05
性別: GG
專(zhuān)業(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
在線(xiàn): 64.2小時(shí)
蟲(chóng)號(hào): 1525201
注冊(cè): 2011-12-05
性別: GG
專(zhuān)業(yè): 玻璃材料

引用回帖:

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

贊一下

回復(fù)此樓

Fighting~~

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

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

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

[考研] 生物學(xué)調(diào)劑招人！��！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋?xiě)?/a>
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	6/300	2026-03-19 14:16 by 申子申申
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線(xiàn)與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 085601專(zhuān)碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu