版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

Brono博博

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 114.5
帖子: 26
在線: 8.6小時(shí)
蟲號(hào): 7939131
注冊(cè): 2018-02-06
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

[求助] 動(dòng)力學(xué)參數(shù)擬合問題求助已有3人參與

目前在做反應(yīng)動(dòng)力學(xué)擬合，可是實(shí)驗(yàn)室無人會(huì)，只能自己摸索，目前知道了動(dòng)力學(xué)方程和一些數(shù)據(jù)，想擬合出參數(shù)，求各位大佬幫忙matlab或1stop擬合一下，有搞動(dòng)力學(xué)的同學(xué)可以一起討論一下，真是太難了，麻煩各位了�。。∠旅媸俏页醪綄懙�1stopt程序，結(jié)果跑出來的數(shù)據(jù)根本不能用，想請(qǐng)教下大家擬合度低究竟是模型參數(shù)太多的原因還是算法有問題，請(qǐng)大神幫幫忙！！
Parameters k(1:10),k1=[0,],k2=[0,],k3=[0,],k4=[0,],k5=[0,],k6=[0,],k7=[0,],k8=[0,],k9=[0,],k10=[0,];
variables t,c(1:5);
ODEFunction
c1'=k1*c2*(c2+c3+c4)/(1+k2*(c2*(c2+c3+c4))^0.5+k3*c2+k4*c3+k5*c4+k6*c5)^2+k7*(c2*(c2+c3+c3))^0.5/(1+k2*(c2*(c2+c3+c3))^0.5+k3*c2+k4*c3+k5*c4+k6*c5)/(1+k8)^2;
c2'=k1*c2*(c2+c3+c4)/(1+k2*(c2*(c2+c3+c3))^0.5+k3*c2+k4*c3+k5*c4+k6*c5)^2-k9*c2/(1+k2*(c2*(c2+c3+c3))^0.5+k3*c2+k4*c3+k5*c4+k6*c5)/(1+k8)^2;
c3'=(k9*c2-k10*c3)/(1+k2*(c2*(c2+c3+c3))^0.5+k3*c2+k4*c3+k5*c4+k6*c5)/(1+k8)^2;
c4'=k10*c3/(1+k2*(c2*(c2+c3+c3))^0.5+k3*c2+k4*c3+k5*c4+k6*c5)/(1+k8)^2;
c5'=k7*(c2*(c2+c3+c3))^0.5/(1+k2*(c2*(c2+c3+c3))^0.5+k3*c2+k4*c3+k5*c4+k6*c5)/(1+k8)^2;
data;//
//t c1 c2 c3 c4 c5
10 5.0249 0.0000 0.9794 0.0000 0.0000
20 4.8894 0.0000 1.3043 0.0000 0.0000
30 4.7004 0.0000 1.5203 0.0000 0.0000
45 4.4588 0.0000 1.7503 0.0000 0.0000
60 4.3103 0.0000 1.7336 0.0000 0.0000
75 4.2412 0.0000 1.7511 0.0000 0.0000
90 4.0658 0.0000 1.7148 0.0000 0.0000
105 3.9587 0.0000 1.7531 0.0000 0.0000
125 3.7984 0.1646 1.7507 0.1039 0.0000
140 3.5823 0.2002 1.6003 0.1406 0.0000
160 3.2773 0.2529 1.4905 0.2211 0.0000
175 2.9786 0.3404 1.4087 0.3433 0.0000
190 2.6945 0.4994 1.3894 0.6631 0.0000
210 2.3301 0.7709 0.9335 0.8718 0.0000
230 2.0823 0.8962 0.7938 1.1992 0.0000
250 1.3216 1.0824 0.3107 1.3042 0.0448
280 0.5662 1.7420 0.1108 1.4905 0.1773
300 0.2618 1.8755 0.0000 1.4112 0.3343
320 0.0000 2.1430 0.0000 0.7001 0.6788
340 0.0000 2.1758 0.0000 0.6303 1.3191
360 0.0000 2.1846 0.0000 0.4146 1.7942

回復(fù)此樓

» 猜你喜歡

幫我的英語口語老師找學(xué)生已經(jīng)有4人回復(fù)
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有51人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有1人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有2人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)

1樓 2020-01-28 20:52:46

已閱回復(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
月只藍(lán): 金幣+1, 鼓勵(lì)交流 2020-01-29 17:57:01
獨(dú)孤神宇: 金幣+1, 鼓勵(lì)交流 2020-02-01 15:43:58
Brono博博: 金幣+3, ★有幫助 2020-02-02 22:00:04

自己先檢查下模型公式是否正確無誤？數(shù)據(jù)是否也無問題？

贊一下(1人)

回復(fù)此樓

2樓2020-01-29 14:48:56

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

lvdq44

木蟲 (著名寫手)

應(yīng)助: 9 (幼兒園)
金幣: 16216.2
紅花: 3
帖子: 1162
在線: 878.5小時(shí)
蟲號(hào): 587512
注冊(cè): 2008-08-22
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
獨(dú)孤神宇: 金幣+5, 鼓勵(lì)交流 2020-02-01 15:43:48
Brono博博: 金幣+7, ★★★很有幫助, 謝謝您！ 2020-02-02 21:28:17

引用回帖:

2樓: Originally posted by dingd at 2020-01-29 14:48:56
自己先檢查下模型公式是否正確無誤？數(shù)據(jù)是否也無問題？

1.首先要確定你所有的微分方程是不是都寫對(duì)了？比如括號(hào)成對(duì)等這問題。2.你的微分方程組等號(hào)的左邊項(xiàng)應(yīng)該是反應(yīng)速率。而這個(gè)反應(yīng)速率在求解動(dòng)力學(xué)參數(shù)時(shí)，應(yīng)該是一個(gè)實(shí)驗(yàn)測(cè)到的已知數(shù)或者是一個(gè)可以解得的數(shù)值。用無梯度反應(yīng)器可以直接測(cè)到反應(yīng)速率。用管式等溫反應(yīng)器則需要對(duì)微分方程進(jìn)行積分。這種情況下，積分反應(yīng)器的進(jìn)出口條件都是實(shí)驗(yàn)測(cè)到的，是已知的。而參數(shù)估計(jì)要做的就是:找到一組參數(shù)，它能使微分方程從已知的反應(yīng)器進(jìn)口條件積分到出口處的結(jié)果，與實(shí)測(cè)的出口組分濃度之間的誤差為極小。3.所以，上述情況，一種是解非線性代數(shù)方程組；另一種是解非線性微分方程組�？茨阌玫氖悄囊环N實(shí)驗(yàn)反應(yīng)器。4.給待求參數(shù)的初值時(shí)，要避免直接給“零”值。需要時(shí)，給個(gè)非常小的正值，例如1.0e－6。上面這些，與你用何種編程語言無關(guān)。

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

贊一下(3人)

回復(fù)此樓

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

Brono博博

3樓2020-02-01 13:49:21

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

Brono博博

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 114.5
帖子: 26
在線: 8.6小時(shí)
蟲號(hào): 7939131
注冊(cè): 2018-02-06
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

送紅花一朵

引用回帖:

3樓: Originally posted by lvdq44 at 2020-02-01 13:49:21
1.首先要確定你所有的微分方程是不是都寫對(duì)了？比如括號(hào)成對(duì)等這問題。2.你的微分方程組等號(hào)的左邊項(xiàng)應(yīng)該是反應(yīng)速率。而這個(gè)反應(yīng)速率在求解動(dòng)力學(xué)參數(shù)時(shí)，應(yīng)該是一個(gè)實(shí)驗(yàn)測(cè)到的已知數(shù)或者是一個(gè)可以解得的數(shù)值。用 ...

謝謝您的指導(dǎo)！還是有個(gè)問題想請(qǐng)教您，我的反應(yīng)動(dòng)力學(xué)模型是基于羥醛縮合加氫過程集成一步的反應(yīng)模型，即整個(gè)反應(yīng)體系包含的是兩個(gè)獨(dú)立的反應(yīng)步驟，即先進(jìn)行羥醛縮合其次對(duì)縮合產(chǎn)物進(jìn)行加氫，以往的研究基本上是分開段分別進(jìn)行動(dòng)力學(xué)的研究，這樣研究的結(jié)果較好，根據(jù)經(jīng)驗(yàn)羥醛縮合步為速率控制步驟，那么我這種集成研究是否會(huì)對(duì)整個(gè)體系的動(dòng)力學(xué)研究有影響呢？
對(duì)您提出問題的回復(fù)：我是在Parr高壓在線取樣反應(yīng)釜上進(jìn)行的動(dòng)力學(xué)研究，反應(yīng)速率是對(duì)微分方程積分進(jìn)而求得的，所以求解的是非線性微分方程組，但結(jié)果求出來的擬合線基本都是直線，這與實(shí)測(cè)各物質(zhì)濃度變化規(guī)律根本不符，所以我想除了模型不合適外還可能是擬合計(jì)算方法有著較大的紕漏（因?yàn)閿?shù)據(jù)擬合這方面從來沒接觸過）。最后還是謝謝您的回復(fù)，感謝！！

贊一下

回復(fù)此樓

4樓2020-02-02 22:00:26

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

Brono博博

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 114.5
帖子: 26
在線: 8.6小時(shí)
蟲號(hào): 7939131
注冊(cè): 2018-02-06
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

2樓: Originally posted by dingd at 2020-01-29 14:48:56
自己先檢查下模型公式是否正確無誤？數(shù)據(jù)是否也無問題？

模型公式和數(shù)據(jù)應(yīng)該是沒問題的

贊一下

回復(fù)此樓

5樓2020-02-02 22:02:24

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★ ★ ★
Brono博博: 金幣+7, ★★★很有幫助 2020-02-05 21:29:18

Parameters k(1:10)；
ParameterDomain = [0,];

參數(shù)比較多，微分方程擬合計(jì)算比較慢，一般迭代幾次就可以大致預(yù)判結(jié)果了。多算幾次試試。。

贊一下(1人)

回復(fù)此樓

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

Brono博博

數(shù)值計(jì)算

6樓2020-02-02 22:06:39

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

Brono博博

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 114.5
帖子: 26
在線: 8.6小時(shí)
蟲號(hào): 7939131
注冊(cè): 2018-02-06
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

送紅花一朵

引用回帖:

6樓: Originally posted by 獨(dú)孤神宇 at 2020-02-02 22:06:39
Parameters k(1:10)；
ParameterDomain = ;

參數(shù)比較多，微分方程擬合計(jì)算比較慢，一般迭代幾次就可以大致預(yù)判結(jié)果了。多算幾次試試。。

好的謝謝您！

回復(fù)此樓

7樓2020-02-05 21:29:55

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 Brono博博的主題更新

返回列表

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

[考研] 材料與化工085600調(diào)劑求老師收留 +4	jiaanl 2026-03-11	4/200	2026-03-11 19:38 by peike
[考研] 材料工程專碩，一志愿中國(guó)礦業(yè)大學(xué)，總分314，求調(diào)劑 +4	無懈可擊的巨人 2026-03-10	4/200	2026-03-10 22:45 by 星空星月
[考研] 化工0817調(diào)劑 +8	燦若星晨 2026-03-10	8/400	2026-03-10 22:44 by 星空星月
[考研] 332材料求調(diào)劑 +6	zjy101327 2026-03-05	8/400	2026-03-10 17:55 by zjy101327
[考研] 311求調(diào)劑 +3	牛乳糖的卡卡 2026-03-10	3/150	2026-03-10 16:19 by 球場(chǎng)大飛機(jī)
[考研] 0703化學(xué)調(diào)劑 +3	三dd. 2026-03-10	3/150	2026-03-10 15:45 by peike
[考研] 0860求調(diào)劑（272分） +3	lllllcsjsj 2026-03-05	4/200	2026-03-10 15:29 by circleffyy
[考研] 08工科 +5	li李樂成 2026-03-06	5/250	2026-03-10 14:51 by 在風(fēng)落中
[考研] 材料專碩調(diào)劑 +7	慕辰123 2026-03-05	10/500	2026-03-10 14:22 by peike
[考研] 278求調(diào)劑 +8	Gale1314 2026-03-06	8/400	2026-03-10 14:18 by Linda Hu
[考研] 320求調(diào)劑 +4	魏zy 2026-03-08	4/200	2026-03-09 16:14 by ruiyingmiao
[考研] 296求調(diào)劑 +4	Xinyu Wu311 2026-03-09	4/200	2026-03-09 11:51 by 勇敢太監(jiān)王公公
[考研] 求0856調(diào)劑 +5	小力氣珂珂 2026-03-08	5/250	2026-03-09 06:42 by houyaoxu
[考研] 0817找調(diào)劑 +6	kk扛 2026-03-08	6/300	2026-03-09 06:38 by houyaoxu
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級(jí)有競(jìng)賽求調(diào)劑 +5	galaxary 2026-03-07	7/350	2026-03-08 04:37 by wang_dand
[考研] 301求調(diào)劑 +5	一二LV 2026-03-07	5/250	2026-03-07 22:20 by 18137688336
[考研] 276求調(diào)劑本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4	路lyh123 2026-03-06	6/300	2026-03-07 11:37 by shdliugang
[考研] 334求調(diào)劑 +6	Trying] 2026-03-05	8/400	2026-03-06 00:16 by wutongshun
[考研] 求材料調(diào)劑 +4	berdmond 2026-03-05	4/200	2026-03-05 19:45 by 黑衣饅頭人
[考研] 304求調(diào)劑 +4	曼殊2266 2026-03-05	4/200	2026-03-05 17:10 by zhukairuo

24小時(shí)熱門版塊排行榜

[求助] 動(dòng)力學(xué)參數(shù)擬合問題求助 已有3人參與

» 猜你喜歡

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

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

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

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

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

[求助] 動(dòng)力學(xué)參數(shù)擬合問題求助已有3人參與