CYYT

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4098.2
散金: 10
紅花: 3
帖子: 161
在線: 195.7小時(shí)
蟲號(hào): 2244949
注冊(cè): 2013-01-16
專業(yè): 能源化工

[求助] 求高版本 1stopt跑一下已有1人參與

Title "fun";
Variable f1,f2,f3,f4,T,k,r1;
Function r1=ln(k1)*exp(E/-8.314/T)*(f1*f3*f3-f4/k)/[1+f1*exp(a1-b1*(1/T-1/508.9))+f2*exp(a2-b2*(1/T-1/508.9))+f3*exp(a3-b3*(1/T-1/508.9))]^3;
data;
0.055401453 0.06568836 0.310940986 0.03388052 505.85 0.375502591 3.637487987
0.062438972 0.073957325 0.350203427 0.03983632 504.45 0.400287073 3.875227293
0.06951785 0.082275676 0.389613788 0.046072104 504.85 0.393030011 3.838005558
0.075342296 0.089100948 0.422011831 0.050927616 504.35 0.402123899 3.969231143
0.062428462 0.074049173 0.349970255 0.0432256 516.85 0.229853036 3.827809524
0.062421901 0.074093139 0.349849554 0.04253792 523.55 0.172092587 3.774803171
0.06116489 0.072336052 0.343314957 0.027342157 490.85 0.75820361 2.327431006
0.061168402 0.072258127 0.343479618 0.021084269 482.65 1.132812726 1.789155527
0.028574399 0.034715226 0.350085136 0.022047021 515.55 0.243329481 1.770600166
0.027986893 0.033932423 0.343158611 0.022257254 493.65 0.662970349 1.680596038
0.067324173 0.03064251 0.285081748 0.03654528 528.45 0.139877688 3.992327182
0.067340041 0.03060705 0.285279281 0.042628301 514.05 0.259951812 4.380735361

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有294人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

1樓 2015-11-01 12:50:07

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
CYYT: 金幣+10, ★★★很有幫助, 謝謝你了。。 2015-11-02 18:51:59

把中括號(hào)“[]”改成小括號(hào)“()”低版本的也能運(yùn)行吧。
公式仔細(xì)檢查下，似乎有過擬合現(xiàn)象，參數(shù)值不唯一：
1：
均方差(RMSE):0.21698029314672
殘差平方和(SSE):0.564965371368438
相關(guān)系數(shù)(R): 0.975158853698942
相關(guān)系數(shù)之平方(R^2): 0.950934789947434
確定系數(shù)(DC): 0.950779667454345
卡方系數(shù)(Chi-Square): 0.130254031694643
F統(tǒng)計(jì)(F-Statistic): 11.6096040370822

參數(shù) 最佳估算
---------- -------------
k1 0.999999988673665
e -95862.8129665484
a1 -867833.168533956
b1 43.9156147666947
a2 -4982.47558760389
b2 0
a3 -0.00833194402515029
b3 1656.35069912403

2：
均方差(RMSE):0.216980293146719
殘差平方和(SSE):0.564965371368435
相關(guān)系數(shù)(R): 0.975158853730298
相關(guān)系數(shù)之平方(R^2): 0.950934790008589
確定系數(shù)(DC): 0.950779667454345
卡方系數(shù)(Chi-Square): 0.130254031979546
F統(tǒng)計(jì)(F-Statistic): 11.6096040370823

參數(shù) 最佳估算
---------- -------------
k1 0.999999988673666
e -95862.8132743037
a1 -17440154.537724
b1 4235845.00626166
a2 -1395.58639416365
b2 -1613072.39980911
a3 -0.00833195730676999
b3 1656.35070431295

贊一下

回復(fù)此樓

2樓2015-11-01 19:01:37

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時(shí)
蟲號(hào): 3929201
注冊(cè): 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計(jì)算模擬

是的，公式里只能有小括號(hào)

贊一下

回復(fù)此樓

一團(tuán)和氣

3樓2015-11-02 15:13:29

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

CYYT

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4098.2
散金: 10
紅花: 3
帖子: 161
在線: 195.7小時(shí)
蟲號(hào): 2244949
注冊(cè): 2013-01-16
專業(yè): 能源化工

引用回帖:

3樓: Originally posted by lipenggg at 2015-11-02 15:13:29
是的，公式里只能有小括號(hào)

恩恩，謝謝提醒。。

回復(fù)此樓

4樓2015-11-02 18:53:01

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

CYYT

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4098.2
散金: 10
紅花: 3
帖子: 161
在線: 195.7小時(shí)
蟲號(hào): 2244949
注冊(cè): 2013-01-16
專業(yè): 能源化工

引用回帖:

3樓: Originally posted by lipenggg at 2015-11-02 15:13:29
是的，公式里只能有小括號(hào)

不知道為什么，我擬合的結(jié)果里出現(xiàn)了“nan”，求解答啊。。

捕獲.JPG

贊一下

回復(fù)此樓

5樓2015-11-02 18:55:37

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

CYYT

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4098.2
散金: 10
紅花: 3
帖子: 161
在線: 195.7小時(shí)
蟲號(hào): 2244949
注冊(cè): 2013-01-16
專業(yè): 能源化工

引用回帖:

2樓: Originally posted by dingd at 2015-11-01 19:01:37
把中括號(hào)“[]”改成小括號(hào)“()”低版本的也能運(yùn)行吧。
公式仔細(xì)檢查下，似乎有過擬合現(xiàn)象，參數(shù)值不唯一：
1：
均方差(RMSE):0.21698029314672
殘差平方和(SSE):0.564965371368438
相關(guān)系數(shù)(R): 0.97515885369 ...

請(qǐng)問，過擬合是什么意思啊，我是小白，不太清楚啊。。

贊一下

回復(fù)此樓

6樓2015-11-02 19:01:15

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

dingd

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

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

引用回帖:

6樓: Originally posted by CYYT at 2015-11-02 19:01:15
請(qǐng)問，過擬合是什么意思啊，我是小白，不太清楚啊。。...

簡(jiǎn)單舉例，如果擬合方程是： y=a+b*c*x
b和c會(huì)有無窮多組組合，這就是過擬合，公式簡(jiǎn)化成： y=a+b*x參數(shù)就穩(wěn)定唯一了。你的公式或許存在類似問題而導(dǎo)致目標(biāo)函數(shù)值一樣但參數(shù)組值不同，自己檢查。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

CYYT

7樓2015-11-02 20:06:39

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時(shí)
蟲號(hào): 3929201
注冊(cè): 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計(jì)算模擬

出現(xiàn)nan，公式仔細(xì)檢查一下。估計(jì)是公式有問題，軟件不識(shí)別。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

CYYT

一團(tuán)和氣

8樓2015-11-03 10:05:27

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時(shí)
蟲號(hào): 3929201
注冊(cè): 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計(jì)算模擬

r1=ln(k1)*exp(E/-8.314/T)*(f1*f3*f3-f4/k)/(1+f1*exp(a1-b1*(1/T-1/508.9))+f2*exp(a2-b2*(1/T-1/508.9))+f3*exp(a3-b3*(1/T-1/508.9)))^3;

贊一下

回復(fù)此樓

一團(tuán)和氣

9樓2015-11-03 10:06:01

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

CYYT

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4098.2
散金: 10
紅花: 3
帖子: 161
在線: 195.7小時(shí)
蟲號(hào): 2244949
注冊(cè): 2013-01-16
專業(yè): 能源化工

送紅花一朵

引用回帖:

7樓: Originally posted by dingd at 2015-11-02 20:06:39
簡(jiǎn)單舉例，如果擬合方程是： y=a+b*c*x
b和c會(huì)有無窮多組組合，這就是過擬合，公式簡(jiǎn)化成： y=a+b*x參數(shù)就穩(wěn)定唯一了。你的公式或許存在類似問題而導(dǎo)致目標(biāo)函數(shù)值一樣但參數(shù)組值不同，自己檢查。...

哦哦，明白了，謝謝你的指點(diǎn)，送個(gè)小花！

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

10樓2015-11-03 16:52:24

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

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

[考研] 化學(xué)調(diào)劑 +3	yzysaa 2026-03-21	3/150	2026-03-21 16:40 by macy2011
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 求調(diào)劑 +4	十三加油 2026-03-21	4/200	2026-03-21 16:20 by barlinike
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	11/550	2026-03-21 15:43 by ColorlessPI
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤(rùn)采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽(yáng)
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽(yáng)w
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu