a_tree

應(yīng)助: 0 (幼兒園)
金幣: 213.5
帖子: 32
在線: 16.4小時(shí)
蟲(chóng)號(hào): 1097439
注冊(cè): 2010-09-13
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

[求助] 求1stOpt高手解決 “函數(shù)表達(dá)式有誤” 已有1人參與

我想求解一個(gè)大型非線性方程組，第一次接觸1stOpt，所以求解了一個(gè)2元方程組練手，但是一直出現(xiàn)錯(cuò)誤，請(qǐng)大家?guī)蛶臀�！只�?個(gè)幣了。代碼如下：
Constant n=2;
Constant M(1:n,1:n)=[1,2,
                  3,4 ];
Constant K(1:n,1:n)=[2,4,
                  0,3];
Constant NGB(1:n)=[0.2,0.5];
parameter x(1:n);
Function For (i = 1 : n)  ((Sum(j=1:n)(M[i,j]*x[j]*x[j]+K[i,j]*x[j] ))-NGB=0);

1@C9)E46{@V6N3[3Z`JXX$C.png

回復(fù)此樓

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有90人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

1stopt程序，求修改已經(jīng)有4人回復(fù)
用含有定積分的隱函數(shù)擬合數(shù)據(jù)，matlab或者origin 已經(jīng)有18人回復(fù)
親們，求助一下，這個(gè)非線性方程組怎么解？已經(jīng)有11人回復(fù)
通過(guò)實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù) 已經(jīng)有14人回復(fù)

1樓 2015-05-13 17:06:35

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

回帖支持 ( 顯示支持度最高的前 50 名 )

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

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

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

下面寫(xiě)法，6.9里沒(méi)問(wèn)題，其它版本就不知道了：

CODE:

Constant n=2;

Constant M(1:n,1:n)=[1,2,

                      3,4 ];

Constant K(1:n,1:n)=[2,4,

                    0,3];

Constant NGB(1:n)=[0.2,0.5];

ConstStr For(i=1:n)(f[i]=(Sum(j=1:n)(M[i,j]*x[j]*x[j]+K[i,j]*x[j] ))-NGB[i]);

parameter x(1:n);

Function For(i=1:n)(f[i]);

x1: -0.167448164499229
x2: 0.119566309158215

贊一下(1人)

回復(fù)此樓

2樓2015-05-13 17:46:12

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

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

★ ★ ★ ★ ★ ★
a_tree: 金幣+6, ★★★很有幫助 2015-05-13 20:05:52

抱歉，筆誤，應(yīng)該是6.0里而非6.9

贊一下(1人)

回復(fù)此樓

3樓2015-05-13 17:47:05

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

普通回帖

a_tree

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

應(yīng)助: 0 (幼兒園)
金幣: 213.5
帖子: 32
在線: 16.4小時(shí)
蟲(chóng)號(hào): 1097439
注冊(cè): 2010-09-13
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

3樓: Originally posted by dingd at 2015-05-13 17:47:05
抱歉，筆誤，應(yīng)該是6.0里而非6.9

非常感謝！您的代碼我在1.5免安裝版上試了一下，提示錯(cuò)誤。我把M,K,NGB的數(shù)據(jù)貼上來(lái)，您能幫我跑一下嗎？

贊一下

回復(fù)此樓

4樓2015-05-13 18:03:43

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

a_tree

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

應(yīng)助: 0 (幼兒園)
金幣: 213.5
帖子: 32
在線: 16.4小時(shí)
蟲(chóng)號(hào): 1097439
注冊(cè): 2010-09-13
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

3樓: Originally posted by dingd at 2015-05-13 17:47:05
抱歉，筆誤，應(yīng)該是6.0里而非6.9

你好 M K NGB 在附件里。n=81。想求解這個(gè)非線性方程組很久了，一直很難得到全局收斂的解。這對(duì)我非常非常重要。非常期待您的幫助！

贊一下

回復(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 : K.txt

2015-05-13 20:07:34, 23.9 K

附件 2 : M.txt

2015-05-13 20:07:35, 23.3 K

附件 3 : NGB.txt

2015-05-13 20:07:35, 522 bytes

5樓2015-05-13 20:09:24

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

a_tree

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

應(yīng)助: 0 (幼兒園)
金幣: 213.5
帖子: 32
在線: 16.4小時(shí)
蟲(chóng)號(hào): 1097439
注冊(cè): 2010-09-13
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

經(jīng)過(guò)反復(fù)檢查和修改，終于找到為什么會(huì)出錯(cuò)了，把For(i=1:n)與((Sum(j=1:n)(M[i,j]*x[j]*x[j]+K[i,j]*x[j] ))-NGB=0);之間的空格去掉就OK了，希望能幫到和我遇到同樣問(wèn)題的人。
Constant n=2;
Constant M(1:n,1:n)=[1,2,
3,4 ];
Constant K(1:n,1:n)=[2,4,
0,3];
Constant NGB(1:n)=[0.2,0.5];
parameter x(1:n);
Function For (i = 1 : n)((Sum(j=1:n)(M[i,j]*x[j]*x[j]+K[i,j]*x[j] ))-NGB=0);

贊一下

回復(fù)此樓

6樓2015-05-13 21:28:49

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

云過(guò)半山

金蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1142.5
散金: 732
紅花: 2
帖子: 102
在線: 141.2小時(shí)
蟲(chóng)號(hào): 1990166
注冊(cè): 2012-09-10
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by dingd at 2015-05-13 17:47:05
抱歉，筆誤，應(yīng)該是6.0里而非6.9

請(qǐng)問(wèn) 在網(wǎng)上只能下到1stopts1.5的版本這個(gè)6版本怎么弄到？購(gòu)買么？

贊一下

回復(fù)此樓

也無(wú)風(fēng)雨也無(wú)晴

7樓2015-05-15 11:15:40

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

引用回帖:

7樓: Originally posted by 云過(guò)半山 at 2015-05-15 11:15:40
請(qǐng)問(wèn) 在網(wǎng)上只能下到1stopts1.5的版本這個(gè)6版本怎么弄到？購(gòu)買么？...

是的。

回復(fù)此樓

8樓2015-05-15 13:21:52

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

a_tree

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

應(yīng)助: 0 (幼兒園)
金幣: 213.5
帖子: 32
在線: 16.4小時(shí)
蟲(chóng)號(hào): 1097439
注冊(cè): 2010-09-13
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

8樓: Originally posted by dingd at 2015-05-15 13:21:52
是的。...

你好，請(qǐng)問(wèn)在1stOpt里怎么設(shè)置最大迭代數(shù)？每次計(jì)算的“計(jì)算終止原因”都是“達(dá)到最大迭代數(shù)”，迭代次數(shù)都顯示為5001.

贊一下

回復(fù)此樓

9樓2015-05-16 11:36:27

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

a_tree

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

應(yīng)助: 0 (幼兒園)
金幣: 213.5
帖子: 32
在線: 16.4小時(shí)
蟲(chóng)號(hào): 1097439
注冊(cè): 2010-09-13
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

8樓: Originally posted by dingd at 2015-05-15 13:21:52
是的。...

6.0版多少元可以買到？

贊一下

回復(fù)此樓

10樓2015-05-16 11:38:04

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

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

[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +4	披星河 2026-03-30	4/200	2026-03-31 04:15 by fmesaito
[考博] 材料專業(yè)申博 +3	杜雨婷dyt 2026-03-29	3/150	2026-03-31 02:04 by nxgogo
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 286分調(diào)劑 +10	Faune 2026-03-30	11/550	2026-03-30 20:35 by 啊李999
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25	稻妻小編 2026-03-26	25/1250	2026-03-30 20:11 by 滴滴上岸呀
[考研] 311求調(diào)劑 +7	藍(lán)月亮亮 2026-03-30	7/350	2026-03-30 19:24 by 無(wú)際的草原
[考研] 289求調(diào)劑 +16	新時(shí)代材料 2026-03-27	16/800	2026-03-30 19:04 by Wang200018
[考研] 297求調(diào)劑 +17	田洪有 2026-03-26	18/900	2026-03-30 18:32 by nothing投稿中
[考研] 生物技術(shù)與工程 +7	1294608413 2026-03-25	8/400	2026-03-30 11:36 by 唐沐兒
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無(wú)際的草原
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 086000調(diào)劑 +3	7901117076 2026-03-26	3/150	2026-03-27 21:34 by Jianing_Mi
[考研] 一志愿西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 321求調(diào)劑 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 085602化學(xué)工程求調(diào)劑。 +4	平樂(lè)樂(lè)樂(lè) 2026-03-26	4/200	2026-03-26 17:57 by fmesaito
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚(yú)的貓
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky