當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

xiaohsu2006

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 23 (小學(xué)生)
金幣: 2520.4
散金: 207
紅花: 4
帖子: 1362
在線: 1137.1小時(shí)
蟲(chóng)號(hào): 1012350
注冊(cè): 2010-05-06
性別: GG
專業(yè): 多相流熱物理學(xué)

[求助] 求1stopt擬合動(dòng)力學(xué)參數(shù)

小弟沒(méi)有高版本的，求高手幫忙，感激不盡！
反應(yīng)動(dòng)力學(xué)：
A+B - 2.5C+3.5D+0.5F+G
A+D - C+F
3C+D - E+A

r1=k1*CB
r2=k2*CD
r3=k3*CD

dCAdt=-r1-r2+r3
dCBdt=-r1
dCCdt=2.5*r1+r2-3*r3
dCDdt=3.5*r1-r2-r3
dCEdt=r3
dCFdt=0.5r1+r2
dCG=r1

初始條件 CA=4.0799, CB=0.0707，CC=CD=CE=CF=CG=0
1
t/s CC CD CE CF
30 2.35E-04 3.36E-03 1.15E-04 4.53E-03
60 1.78E-03 8.62E-03 9.35E-04 6.27E-03
120 2.92E-03 1.20E-02 2.27E-03 6.74E-03
180 4.43E-03 1.54E-02 2.93E-03 6.99E-03
2

t/s CC CD CE CF
30 7.33E-04 5.87E-03 4.50E-04 5.01E-03
60 3.07E-03 1.36E-02 2.67E-03 6.25E-03
120 5.07E-03 1.66E-02 4.42E-03 6.62E-03
180 5.86E-03 1.71E-02 5.30E-03 7.29E-03
3

t/s CC CD CE CF
30 1.67E-03 9.37E-03 1.08E-03 4.77E-03
60 5.55E-03 1.68E-02 4.33E-03 5.91E-03
120 8.29E-03 1.77E-02 6.02E-03 7.50E-03
180 1.00E-02 1.68E-02 6.28E-03 8.12E-03
分別擬合出三個(gè)溫度下的反應(yīng)速率常數(shù)k1,k2,k3.

回復(fù)此樓

» 猜你喜歡

中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有197人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(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)有2人回復(fù)

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

拜求！matlab擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有14人回復(fù)
最小二乘解算擬合方程已經(jīng)有12人回復(fù)
求助用matlab擬合動(dòng)力學(xué)方程已經(jīng)有9人回復(fù)
matlab 擬合方法已經(jīng)有4人回復(fù)
matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù) 已經(jīng)有8人回復(fù)
1stopt擬合問(wèn)題求助已經(jīng)有7人回復(fù)
rmse，sse計(jì)算求助已經(jīng)有3人回復(fù)
如何用matlab對(duì)四維數(shù)據(jù)進(jìn)行非線性擬合？已經(jīng)有5人回復(fù)
通過(guò)實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù) 已經(jīng)有14人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒(méi)有遇到過(guò)這種情況？已經(jīng)有8人回復(fù)
用1stopt或者origin非線性擬合磁導(dǎo)率已經(jīng)有5人回復(fù)
matlab數(shù)據(jù)擬合問(wèn)題求助已經(jīng)有7人回復(fù)
matlab求解動(dòng)力學(xué) 方程中的參數(shù) 已經(jīng)有13人回復(fù)
動(dòng)力學(xué)參數(shù)擬合已經(jīng)有26人回復(fù)
matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

MATLAB擬合工具如何引入誤差棒？已經(jīng)有11人回復(fù)
學(xué)1stopt要多久？已經(jīng)有12人回復(fù)
求可以做曲面擬合的軟件已經(jīng)有10人回復(fù)
【求助】擬合動(dòng)力學(xué)方程求助已經(jīng)有13人回復(fù)
用1stopt進(jìn)行數(shù)據(jù)擬合時(shí)，無(wú)法運(yùn)行sharedmodel ,求助高手已經(jīng)有7人回復(fù)

1樓 2014-03-21 16:53:16

已閱回復(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

引用回帖:

6樓: Originally posted by xiaohsu2006 at 2014-04-12 19:02:41
你好，能不能幫我再跑個(gè)代碼，寫(xiě)好了！麻煩你了，先謝過(guò)！
ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD;
Variable t,CC,CD,CE,CF;
InitialODEValue t=0, CA=4079.9, CB=70.7,CC=0, CD=0, CE=0, CF=0, CG=0;
ODEFunct ...

均方差(RMSE):2.06284664311107
殘差平方和(SSE):68.0853803679137
相關(guān)系數(shù)(R): 0.858287409589619
相關(guān)系數(shù)之平方(R^2): 0.736657277460058
確定系數(shù)(DC): -2.62363520042484
F統(tǒng)計(jì)(F-Statistic): 0.302815847593867

參數(shù) 最佳估算
-------------------- -------------
k1 0.000571141146818135
k2 0.00237051259796401
k3 0.00354131349016907

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) cc 計(jì)算 cc 目標(biāo) cd 計(jì)算 cd 目標(biāo) ce 計(jì)算 ce 目標(biāo) cf 計(jì)算 cf
1 30 0.235008523 2.51033673349595 3.363239986 3.85109149173353 0.115121201 0.211249848393764 4.531327956 0.74194373907573
2 60 1.780010877 4.10425865484917 8.616751392 7.01089022534024 0.934615204 0.793820390023803 6.268066726 1.72224301149238
3 120 2.923654401 5.14738702699264 12.0012815 11.6919857722601 2.266319353 2.81501817652468 6.737997984 4.22595972144955
4 180 4.433302192 4.10534301010109 15.40323797 14.746973106494 2.929517334 5.64773674450559 6.994715676 7.23413221207872

贊一下

回復(fù)此樓

7樓2014-04-12 19:52:02

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

查看全部 8 個(gè)回答

xiaohsu2006

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 23 (小學(xué)生)
金幣: 2520.4
散金: 207
紅花: 4
帖子: 1362
在線: 1137.1小時(shí)
蟲(chóng)號(hào): 1012350
注冊(cè): 2010-05-06
性別: GG
專業(yè): 多相流熱物理學(xué)

沒(méi)人么？自己先頂一下！

回復(fù)此樓

2樓2014-03-24 10:23:41

已閱回復(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)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
xiaohsu2006: 金幣+90, ★★★★★最佳答案 2014-03-24 11:33:21

代碼及第一組數(shù)據(jù)結(jié)果：

CODE:

ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD;

Variable t,CC,CD,CE,CF;

InitialODEValue t=0, CA=4.0799, CB=0.0707,CC=0, CD=0, CE=0, CF=0, CG=0;

ODEFunction CA'=-r1-r2+r3;

        CB'=-r1;

        CC'=2.5*r1+r2-3*r3;

        CD'=3.5*r1-r2-r3;

        CE'=r3;

        CF'=0.5*r1+r2;

        CG'=r1;

Data;

//t/s        CC        CD        CE        CF

30        2.35E-04        3.36E-03        1.15E-04        4.53E-03

60        1.78E-03        8.62E-03        9.35E-04        6.27E-03

120        2.92E-03        1.20E-02        2.27E-03        6.74E-03

180        4.43E-03        1.54E-02        2.93E-03        6.99E-03

均方差(RMSE):0.00206339294816877
殘差平方和(SSE):6.81214473368416E-5
相關(guān)系數(shù)(R): 0.857496435844223
相關(guān)系數(shù)之平方(R^2): 0.735300137485547
確定系數(shù)(DC): -2.62619250526438
F統(tǒng)計(jì)(F-Statistic): 0.302254888086507

參數(shù) 最佳估算
-------------------- -------------
k1 0.000571090248522173
k2 0.00237024445720707
k3 0.00354220354969495

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) cc 計(jì)算 cc 目標(biāo) cd 計(jì)算 cd 目標(biāo) ce 計(jì)算 ce 目標(biāo) cf 計(jì)算 cf
1 30 0.000235 0.00250994279834437 0.00336 0.00385071637240616 0.000115 0.000211282929925092 0.00453 0.000741861291302684
2 60 0.00178 0.0041032640107171 0.00862 0.00701015349688254 0.000935 0.000793940277216935 0.00627 0.00172202530354545
3 120 0.00292 0.00514480375607497 0.012 0.0116906247177385 0.00227 0.00281539943089081 0.00674 0.00422532664962184
4 180 0.00443 0.00410080160481611 0.0154 0.0147450856032192 0.00293 0.00564846729296284 0.00699 0.00723295135707574

贊一下

回復(fù)此樓

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

xiaohsu2006

3樓2014-03-24 10:57:31

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

xiaohsu2006

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 23 (小學(xué)生)
金幣: 2520.4
散金: 207
紅花: 4
帖子: 1362
在線: 1137.1小時(shí)
蟲(chóng)號(hào): 1012350
注冊(cè): 2010-05-06
性別: GG
專業(yè): 多相流熱物理學(xué)

送紅花一朵

引用回帖:

3樓: Originally posted by dingd at 2014-03-24 10:57:31
代碼及第一組數(shù)據(jù)結(jié)果：

ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD;
Variable t,CC,CD,CE,CF;
InitialODEValue t=0, CA=4.0799, CB=0.0707,CC=0, CD=0, CE=0, CF=0, CG=0;
ODEFunction CA'=-r1-r2+r3;
C ...

非常感謝您！金幣給你了，能不能幫忙再跑下后面兩個(gè)，我再追加金幣！謝謝了！

贊一下

回復(fù)此樓

4樓2014-03-24 11:35:56

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

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

[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 化學(xué)308分調(diào)劑 +7	你好明天你好 2026-03-30	7/350	2026-03-30 11:02 by 1172367218
[考研] 085600材料與化工調(diào)劑 +6	kikiki7 2026-03-30	6/300	2026-03-30 10:57 by 星空星月
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +6	@taotao 2026-03-26	7/350	2026-03-30 10:43 by 我是小康
[考研] 281求調(diào)劑 +4	亞克西good 2026-03-26	6/300	2026-03-29 21:06 by 啊李999
[考研] 332求調(diào)劑 +4	@MZB382400 2026-03-28	4/200	2026-03-28 21:02 by 唐沐兒
[考研] 求調(diào)劑 +6	蘆lty 2026-03-25	7/350	2026-03-28 13:13 by 唐沐兒
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +5	dongfang59 2026-03-27	5/250	2026-03-28 04:53 by wxiongid
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito
[考研] 求調(diào)劑推薦材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 0856調(diào)劑 +5	求求讓我有書(shū)讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 343求調(diào)劑 +4	贈(zèng)我一本書(shū) 2026-03-23	4/200	2026-03-27 00:40 by wxiongid
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚(yú)的貓
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243