troyzz

木蟲 (初入文壇)

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

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

在網(wǎng)上看了一遍關(guān)于配料比優(yōu)化的論文，完全按照論文內(nèi)容試著編了一下，可是無法運(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)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 程序語(yǔ)言.docx

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

» 猜你喜歡

309求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)

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

給小老板認(rèn)識(shí)的老外寫郵件，請(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)問題，煩請(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)，怎么寫到txt里？已經(jīng)有14人回復(fù)
有會(huì)看fortran程序的，能幫忙看看這個(gè)程序嗎？已經(jīng)有10人回復(fù)
內(nèi)循環(huán)問題，求大蝦幫忙已經(jīng)有3人回復(fù)
JCA審稿意見返回，編輯要求很嚴(yán)厲，語(yǔ)言的問題，想找一個(gè)公司幫忙改，大家推薦下已經(jīng)有12人回復(fù)
求幫忙，關(guān)于含有bessel函數(shù)的方程的編程求解問題已經(jīng)有11人回復(fù)
寫了一個(gè)fortran90的小程序，編譯通不過，請(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ù)
【求助】幫忙改寫展寬程序已經(jīng)有11人回復(fù)
【求助】幫忙看一下，問題程序如何修改已經(jīng)有3人回復(fù)

Fighting~~

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

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

ausir

木蟲 (正式寫手)

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

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

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

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

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

贊一下

回復(fù)此樓

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

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

ausir

木蟲 (正式寫手)

應(yīng)助: 18 (小學(xué)生)
金幣: 2009.2
紅花: 4
帖子: 500
在線: 178.9小時(shí)
蟲號(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è)問題：

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 矩陣中的寫法，1000 寫作 1 000， ...

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

贊一下

回復(fù)此樓

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

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

troyzz

木蟲 (初入文壇)

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

引用回帖:

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

原來是復(fù)制過程中出現(xiàn)很多錯(cuò)誤，我又仔細(xì)的看了一下那篇論文，并參考了您的意見，自己演算了一遍，附上修改后的語(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

木蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1733.3
帖子: 42
在線: 64.2小時(shí)
蟲號(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ā)表

[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8	柳溪邊 2026-03-26	8/400	2026-03-29 20:42 by 唐沐兒
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運(yùn)的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 0856求調(diào)劑 +7	楒桉 2026-03-28	7/350	2026-03-29 08:28 by fmesaito
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5	niko- 2026-03-27	5/250	2026-03-28 14:25 by 唐沐兒
[考研] 311求調(diào)劑 +9	lin0039 2026-03-26	9/450	2026-03-28 13:05 by 唐沐兒
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語(yǔ)已過四六級(jí)，六級(jí)521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 322求調(diào)劑 +4	我真的很想學(xué)習(xí) 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 333求調(diào)劑 +6	wfh030413@ 2026-03-23	6/300	2026-03-26 22:45 by 學(xué)員8dgXkO
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106