 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2254 | 回復(fù): 6 | ||
Brono博博鐵蟲 (初入文壇)
|
[求助]
動(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é)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有177人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
2樓2020-01-29 14:48:56
lvdq44
木蟲 (著名寫手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 16186.2
- 紅花: 3
- 帖子: 1156
- 在線: 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
獨(dú)孤神宇: 金幣+5, 鼓勵(lì)交流 2020-02-01 15:43:48
Brono博博: 金幣+7, ★★★很有幫助, 謝謝您！ 2020-02-02 21:28:17
|
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ù)方程組；另一種是解非線性微分方程組。看你用的是哪一種實(shí)驗(yàn)反應(yīng)器。4.給待求參數(shù)的初值時(shí)，要避免直接給“零”值。需要時(shí)，給個(gè)非常小的正值，例如1.0e－6。上面這些，與你用何種編程語言無關(guān)。 發(fā)自小木蟲Android客戶端 |
» 本帖已獲得的紅花（最新10朵）
3樓2020-02-01 13:49:21
Brono博博
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 114.5
- 帖子: 26
- 在線: 8.6小時(shí)
- 蟲號(hào): 7939131
- 注冊(cè): 2018-02-06
- 性別: GG
- 專業(yè): 化學(xué)反應(yīng)工程
送紅花一朵 |
謝謝您的指導(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ù)，感謝！！ |
4樓2020-02-02 22:00:26
Brono博博
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 114.5
- 帖子: 26
- 在線: 8.6小時(shí)
- 蟲號(hào): 7939131
- 注冊(cè): 2018-02-06
- 性別: GG
- 專業(yè): 化學(xué)反應(yīng)工程
5樓2020-02-02 22:02:24
獨(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ì)算模擬
6樓2020-02-02 22:06:39
Brono博博
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 114.5
- 帖子: 26
- 在線: 8.6小時(shí)
- 蟲號(hào): 7939131
- 注冊(cè): 2018-02-06
- 性別: GG
- 專業(yè): 化學(xué)反應(yīng)工程
7樓2020-02-05 21:29:55
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Brono博博 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料調(diào)劑 總分330 +3 | 池池丶 2026-03-03 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
|
[考研] 304分材料專碩求調(diào)劑 +8 | qiuzhigril 2026-03-03 | 10/500 |
|
|
|
[考研] 0856材料與化工，270求調(diào)劑 +14 | YXCT 2026-03-01 | 16/800 |
|
|
|
[考研] 347求調(diào)劑 +5 | 啊歐歐歐 2026-03-03 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿西工大材料學(xué)，英一數(shù)二總分321分，求調(diào)劑。 +3 | zz05zz 2026-03-02 | 5/250 |
|
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研]
材料工程專碩283求調(diào)劑
 5+8
|
，！？ 2026-03-02 | 10/500 |
|
|
|
[考研] 302材料工程求調(diào)劑 +5 | Doleres 2026-03-01 | 6/300 |
|
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。，。，。，。i 2026-03-02 | 3/150 |
|
|
|
[考研] 281求調(diào)劑 +5 | 2026計(jì)算機(jī)_誠(chéng)心 2026-03-01 | 8/400 |
|
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|