xiao__guan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 60
帖子: 15
在線: 13.8小時
蟲號: 2086322
注冊: 2012-10-25
專業(yè): 有色金屬冶金

[求助] 1stopt 計算問題

命令如下：

Title "Type your title here";
constant k=1863;
constant m1=14.86;
constant m2=66.80;
constant m3=1.46;
constant m4=17.66;
Parameters x(1:20)[0,1];
Function
0=exp((13836+25.44*k)/8.314/k)*x1^3*x2-x5;
0=exp((35288+204.59*k)/8.314/k)*x1^12*x2^7-x6;
0=exp((18946+25.44*k)/8.314/k)*x1*x2-x7;
0=exp((18566+23.68*k)/8.314/k)*x1*x2^2-x8;
0=exp((-4752.67-75.14*k)/8.314/k)*x1*x2^6-x9;
0=exp((107800+13.44*k)/8.314/k)*x1^3*x3^2-x10;
0=exp((122731+8.85*k)/8.314/k)*x1^2*x3-x11;
0=exp((302587+66.24*k)/8.314/k)*x1^3*x3^2-x12;
0=exp((86231+0.995*k)/8.314/k)*x1*x3-x13;
0=exp((-19512+28.19*k)/8.314/k)*x2^3*x3^2-x14;
0=exp((354932-417.63*k)/8.314/k)*x1^2*x2*x3-x15;
0=exp((1252543-277.06*k)/8.314/k)*x1*x2*x3^3-x16;
0=exp((203172+21.38*k)/8.314/k)*x1^3*x4^2-x17;
0=exp((75158+21.38*k)/8.314/k)*x1*x4-x18;
0=exp((313417-282.66*k)/8.314/k)*x1*x2*x3-x19;
0=exp((2279-13.54*k)/8.314/k)*x2*x4-x20;
0=(0.5*x1+3*x5+12*x6+x7+x8+x9+3*x10+2*x11+3*x12+x13+2*x15+x16+3*x17+x18+x19)*(m1/40+m2/102+m3/64+m4/84)-m1/40;
0=(x2+x5+7*x6+x7+2*x8+6*x9+3*x14+x15+x16+x20)*(m1/40+m2/102+m3/64+m4/84)-m2/102;
0=(x3+x10+x11+2*x12+x13+2*x14+x15+2*x16+x19)*(m1/40+m2/102+m3/64+m4/84)-m3/64;
0=(x4+2*x17+x18+x19+x20)*(m1/40+m2/102+m3/64+m4/84)-m4/84;
0=(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20)-1;

這個程序計算出的結(jié)果：
迭代數(shù): 10740
計算用時(時:分:秒:毫秒): 00:00:22:238
計算中止原因: 達(dá)到收斂判定標(biāo)準(zhǔn)
優(yōu)化算法: 標(biāo)準(zhǔn)簡面體爬山法 + 通用全局優(yōu)化法
目標(biāo)函數(shù)值: 0.337300657640182
x1: 0.00280473183833806
x2: 0.389556948473223
x3: 2.47955133843895E-9
x4: 0.028798627469645
x5: 4.47794033080389E-7
x6: 7.19960278061123E-12
x7: 0.079174458763718
x8: 0.0243537011949764
x9: 7.21467656804983E-12
x10: 0.0181130086301488
x11: 2.49572142641952E-10
x12: 3.93442470085121E-11
x13: 2.20909868564788E-9
x14: 2.29566803469628E-11
x15: 4.50226794231978E-11
x16: 7.36016488772663E-9
x17: 0.000119112732787977
x18: 0.13534118437126
x19: 7.8257263248991E-12
x20: 0.00255012081464932

請問其中的目標(biāo)函數(shù)值是什么意思？
另外為什么計算出的數(shù)值相加不符合第21個方程等于1 的條件？
能否用其他方法進(jìn)一步計算其數(shù)值？
能否用其他方法估算計算的初值使其可以應(yīng)用于matlab或其他需要賦初值的軟件中？

回復(fù)此樓

» 猜你喜歡

材料求調(diào)劑一志愿哈工大總分298分，前三科223分已經(jīng)有5人回復(fù)
085600 286分材料求調(diào)劑已經(jīng)有6人回復(fù)
086502化學(xué)工程342求調(diào)劑已經(jīng)有4人回復(fù)
材料求調(diào)劑一志愿哈工大324 已經(jīng)有7人回復(fù)
328求調(diào)劑已經(jīng)有7人回復(fù)
340求調(diào)劑已經(jīng)有5人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑推薦材料 304 已經(jīng)有15人回復(fù)
291求調(diào)劑已經(jīng)有19人回復(fù)
331環(huán)境科學(xué)與工程求調(diào)劑已經(jīng)有3人回復(fù)

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

1stopt 擬合計算已經(jīng)有6人回復(fù)
畢不了業(yè)了，發(fā)帖求助Origin8.0多元非線性擬合的問題已經(jīng)有9人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計）已經(jīng)有11人回復(fù)
1stOPT1.5版每次仿真結(jié)果不一樣，求助！急！已經(jīng)有3人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
動力學(xué)方程參數(shù)估計方法已經(jīng)有14人回復(fù)
1stopt運行程序，總是出現(xiàn)字符函數(shù)錯誤的提示，請高手幫助已經(jīng)有5人回復(fù)
大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時候有沒有遇到過這種情況？已經(jīng)有8人回復(fù)
matlab數(shù)據(jù)擬合問題求助已經(jīng)有7人回復(fù)
origin或1stopt計算積分面積已經(jīng)有4人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
Gulp 擬合勢函數(shù) 初值設(shè)置的問題已經(jīng)有12人回復(fù)
origin進(jìn)行非線性擬合如何選擇函數(shù)或者自定義合適的函數(shù) 已經(jīng)有19人回復(fù)
1stOpt求解常微分方程邊值問題已經(jīng)有5人回復(fù)
求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
為什么同樣的數(shù)據(jù)，origin 和excel做成的擬合公式不同？已經(jīng)有13人回復(fù)
請教下已知曲線以及曲線上的各個散點如何求出這條曲線的方程表達(dá)式已經(jīng)有10人回復(fù)
用1stopt進(jìn)行數(shù)據(jù)擬合時，無法運行sharedmodel ,求助高手已經(jīng)有7人回復(fù)

1樓 2013-09-03 18:12:20

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
xiao__guan: 金幣+2, ★有幫助 2013-09-03 20:57:12

把解代回原方程的殘差和，表示解的精度。
21個方程20個未知數(shù)，很難滿足所有方程，解只能是近似的。

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2013-09-03 18:41:55

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

xiao__guan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 60
帖子: 15
在線: 13.8小時
蟲號: 2086322
注冊: 2012-10-25
專業(yè): 有色金屬冶金

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2013-09-03 18:41:55
把解代回原方程的殘差和，表示解的精度。
21個方程20個未知數(shù)，很難滿足所有方程，解只能是近似的。

那函數(shù)值是越小越好嗎？
算了很多次，目標(biāo)函數(shù)值都在0.34 左右，是否表示該方法已經(jīng)無法繼續(xù)優(yōu)化？？

贊一下

回復(fù)此樓

3樓2013-09-03 20:56:59

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

★ ★ ★
xiao__guan: 金幣+3, ★★★很有幫助 2013-09-16 16:36:32

引用回帖:

3樓: Originally posted by xiao__guan at 2013-09-03 20:56:59
那函數(shù)值是越小越好嗎？
算了很多次，目標(biāo)函數(shù)值都在0.34 左右，是否表示該方法已經(jīng)無法繼續(xù)優(yōu)化？？...

理想狀態(tài)下為0。21個方程20個未知數(shù)，很少有完美解的。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

4樓2013-09-03 20:59:38

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

xiao__guan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 60
帖子: 15
在線: 13.8小時
蟲號: 2086322
注冊: 2012-10-25
專業(yè): 有色金屬冶金

引用回帖:

4樓: Originally posted by 月只藍(lán) at 2013-09-03 20:59:38
理想狀態(tài)下為0。21個方程20個未知數(shù)，很少有完美解的。...

即使將最后一個方程刪去，解依然變化不大啊

贊一下

回復(fù)此樓

5樓2013-09-04 08:26:23

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

dingd

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

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+2, 辛苦了 2013-09-04 13:06:47

去掉最后一個方程，5.0算的結(jié)果：

目標(biāo)函數(shù)值(最小): 1.0648202469962E-21
x1: 0.00437883851908341
x2: 0.34215847834201
x3: 2.38789778950921E-7
x4: 0.0197863611761941
x5: 1.49674494864709E-6
x6: 1.98891513976249E-11
x7: 0.108569800331441
x8: 0.0293322423655484
x9: 6.15537849483763E-10
x10: 7.05738130338043E-12
x11: 3.66667811347899E-8
x12: 4.21253375613354E-9
x13: 3.08443362662997E-7
x14: 1.14635795502336E-11
x15: 2.10731677895587E-11
x16: 0.00905621804209686
x17: 0.000213969107469129
x18: 0.145174955071669
x19: 6.36062068449683E-12
x20: 0.00153890350917598

贊一下

回復(fù)此樓

6樓2013-09-04 08:47:26

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

★
jjdg: 金幣+1, 感謝建議 2013-09-04 13:07:10

引用回帖:

5樓: Originally posted by xiao__guan at 2013-09-04 08:26:23
即使將最后一個方程刪去，解依然變化不大啊...

參考專家顧問 dingd 的結(jié)果，他的結(jié)果應(yīng)該是最好的。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

7樓2013-09-04 10:46:16

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

xiao__guan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 60
帖子: 15
在線: 13.8小時
蟲號: 2086322
注冊: 2012-10-25
專業(yè): 有色金屬冶金

引用回帖:

7樓: Originally posted by 月只藍(lán) at 2013-09-04 10:46:16
參考專家顧問 dingd 的結(jié)果，他的結(jié)果應(yīng)該是最好的。...

那請問如何限制matlab求方程時解的范圍？類似于Parameters x(1:20)[0,1]這個命令？

贊一下

回復(fù)此樓

8樓2013-09-04 12:24:49

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

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

[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +6	辰沐5211314 2026-03-26	6/300	2026-03-28 03:31 by fmesaito
[考研] 265求調(diào)劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 277跪求調(diào)劑 +5	1915668 2026-03-27	8/400	2026-03-27 19:53 by WYUMater
[考研] 286求調(diào)劑 +8	PolarBear11 2026-03-26	8/400	2026-03-27 18:05 by yu221
[考研] 化學(xué)調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考研] 322求調(diào)劑 +4	舊吢 2026-03-24	4/200	2026-03-27 15:38 by 不吃魚的貓
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 299求調(diào)劑 +6	嗯嗯嗯嗯2 2026-03-27	6/300	2026-03-27 12:43 by 果果媽咪
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 材料調(diào)劑 5+4	想要一壺桃花水 2026-03-25	10/500	2026-03-26 19:56 by 不吃魚的貓
[考研] 085600 材料與化工 329分求調(diào)劑 +9	Mr. Z 2026-03-25	9/450	2026-03-26 10:36 by baoball
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪�。浚�
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte