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

xiaohsu2006

木蟲 (著名寫手)

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

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

小弟沒有高版本的，求高手幫忙，感激不盡！
反應(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ù)此樓

» 猜你喜歡

計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有135人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(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擬合問題求助已經(jīng)有7人回復(fù)
rmse，sse計(jì)算求助已經(jīng)有3人回復(fù)
如何用matlab對(duì)四維數(shù)據(jù)進(jìn)行非線性擬合？已經(jīng)有5人回復(fù)
通過實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù) 已經(jīng)有14人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒有遇到過這種情況？已經(jīng)有8人回復(fù)
用1stopt或者origin非線性擬合磁導(dǎo)率已經(jīng)有5人回復(fù)
matlab數(shù)據(jù)擬合問題求助已經(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í)，無法運(yùn)行sharedmodel ,求助高手已經(jīng)有7人回復(fù)

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

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

xiaohsu2006

木蟲 (著名寫手)

應(yīng)助: 23 (小學(xué)生)
金幣: 2520.4
散金: 207
紅花: 4
帖子: 1362
在線: 1137.1小時(shí)
蟲號(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è)代碼，最近用matlab擬合，但是初值取不好，代碼寫好了，發(fā)站內(nèi)短消息了，金幣隨后奉上，謝謝！

贊一下

回復(fù)此樓

8樓2014-05-16 14:14:51

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

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

xiaohsu2006

木蟲 (著名寫手)

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

沒人么？自己先頂一下！

回復(fù)此樓

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

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

木蟲 (著名寫手)

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

送紅花一朵

引用回帖:

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

贊一下

回復(fù)此樓

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

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

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

[考研] 319求調(diào)劑 +9	太容易1018 2026-04-01	9/450	2026-04-01 21:58 by 我的船我的海
[考研] 070300化學(xué)279求調(diào)劑 +15	哈哈哈^_^ 2026-03-31	17/850	2026-04-01 21:37 by 給你你注意休息
[考研] 332求調(diào)劑 +8	Lyy930824@ 2026-03-29	8/400	2026-04-01 18:40 by 千鈞澄玉宇啊
[考研] 一志愿北交大材料工程總分358 +7	cs0106 2026-04-01	8/400	2026-04-01 18:34 by 記事本2026
[考研] 324分 085600材料與化工 +17	呆鵝oor 2026-03-27	17/850	2026-04-01 17:51 by JYD2011
[考研] 材料專業(yè)調(diào)劑 +5	啦啦啦哭 2026-03-31	6/300	2026-04-01 16:48 by JourneyLucky
[碩博家園] 考研調(diào)劑 +5	駱駝男人 2026-04-01	5/250	2026-04-01 14:28 by syjjj0321
[考研] 324求調(diào)劑 +4	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 14:19 by ZXlzxl0425
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +15	Doleres 2026-03-31	16/800	2026-04-01 09:49 by lfj11
[考研] 349求調(diào)劑 +6	zwjjjjjj 2026-03-31	6/300	2026-04-01 09:16 by JourneyLucky
[考研] 初試301，代碼085701環(huán)境工程，本碩一致，四六級(jí)已過，有二區(qū)一作，共發(fā)表5篇論文 +3	axibli 2026-04-01	3/150	2026-04-01 08:43 by i_cooler
[考研] 322求調(diào)劑 +4	熹僖XX 2026-03-31	4/200	2026-04-01 08:21 by JourneyLucky
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +9	手機(jī)用戶 2026-03-26	9/450	2026-03-31 18:09 by 544594351
[考研] 學(xué)碩274求調(diào)劑 +17	Li李魚 2026-03-26	17/850	2026-03-31 15:19 by 客爾美德
[考研] 調(diào)劑求院校招收 +7	鶴鯨鴿 2026-03-28	7/350	2026-03-31 11:21 by oooqiao
[考研] 274求調(diào)劑 +6	xiao愛同學(xué) 2026-03-30	6/300	2026-03-31 10:04 by cal0306
[考研] 08工科求調(diào)劑286 +5	tgs_001 2026-03-28	5/250	2026-03-31 08:18 by 一只好果子?
[考研] 求調(diào)劑 +7	爭取九點(diǎn)睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點(diǎn)睡
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)