 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 計算模擬 » 程序代碼 » 動力學參數(shù)擬合問題求助
| 5 | 1/1 | 返回列表 |
| 查看: 2267 | 回復: 6 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
Brono博博鐵蟲 (初入文壇)
|
[求助]
動力學參數(shù)擬合問題求助 已有3人參與
|
|
|
目前在做反應動力學擬合，可是實驗室無人會，只能自己摸索，目前知道了動力學方程和一些數(shù)據(jù)，想擬合出參數(shù)，求各位大佬幫忙matlab或1stop擬合一下，有搞動力學的同學可以一起討論一下，真是太難了，麻煩各位了�。。∠旅媸俏页醪綄懙�1stopt程序，結果跑出來的數(shù)據(jù)根本不能用，想請教下大家擬合度低究竟是模型參數(shù)太多的原因還是算法有問題，請大神幫幫忙��！ 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 |
回復此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
幫我的英語口語老師找學生
已經(jīng)有4人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有201人回復
-
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復
-
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復
-
哈爾濱理工大學物理系招收物理學考研調劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調劑
已經(jīng)有2人回復
Brono博博
鐵蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 114.5
- 帖子: 26
- 在線: 8.6小時
- 蟲號: 7939131
- 注冊: 2018-02-06
- 性別: GG
- 專業(yè): 化學反應工程
5樓2020-02-02 22:02:24
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
2樓2020-01-29 14:48:56
lvdq44
木蟲 (著名寫手)
- 應助: 9 (幼兒園)
- 金幣: 16211.2
- 紅花: 3
- 帖子: 1161
- 在線: 878.5小時
- 蟲號: 587512
- 注冊: 2008-08-22
- 性別: GG
- 專業(yè): 化學反應工程
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
獨孤神宇: 金幣+5, 鼓勵交流 2020-02-01 15:43:48
Brono博博: 金幣+7, ★★★很有幫助, 謝謝您！ 2020-02-02 21:28:17
獨孤神宇: 金幣+5, 鼓勵交流 2020-02-01 15:43:48
Brono博博: 金幣+7, ★★★很有幫助, 謝謝您！ 2020-02-02 21:28:17
|
1.首先要確定你所有的微分方程是不是都寫對了？比如括號成對等這問題。2.你的微分方程組等號的左邊項應該是反應速率。而這個反應速率在求解動力學參數(shù)時，應該是一個實驗測到的已知數(shù)或者是一個可以解得的數(shù)值。用無梯度反應器可以直接測到反應速率。用管式等溫反應器則需要對微分方程進行積分。這種情況下，積分反應器的進出口條件都是實驗測到的，是已知的。而參數(shù)估計要做的就是:找到一組參數(shù)，它能使微分方程從已知的反應器進口條件積分到出口處的結果，與實測的出口組分濃度之間的誤差為極小。3.所以，上述情況，一種是解非線性代數(shù)方程組；另一種是解非線性微分方程組�？茨阌玫氖悄囊环N實驗反應器。4.給待求參數(shù)的初值時，要避免直接給“零”值。需要時，給個非常小的正值，例如1.0e－6。上面這些，與你用何種編程語言無關。 發(fā)自小木蟲Android客戶端 |
» 本帖已獲得的紅花（最新10朵）
3樓2020-02-01 13:49:21
Brono博博
鐵蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 114.5
- 帖子: 26
- 在線: 8.6小時
- 蟲號: 7939131
- 注冊: 2018-02-06
- 性別: GG
- 專業(yè): 化學反應工程
送紅花一朵 |
謝謝您的指導！還是有個問題想請教您，我的反應動力學模型是基于羥醛縮合加氫過程集成一步的反應模型，即整個反應體系包含的是兩個獨立的反應步驟，即先進行羥醛縮合其次對縮合產(chǎn)物進行加氫，以往的研究基本上是分開段分別進行動力學的研究，這樣研究的結果較好，根據(jù)經(jīng)驗羥醛縮合步為速率控制步驟，那么我這種集成研究是否會對整個體系的動力學研究有影響呢？ 對您提出問題的回復：我是在Parr高壓在線取樣反應釜上進行的動力學研究，反應速率是對微分方程積分進而求得的，所以求解的是非線性微分方程組，但結果求出來的擬合線基本都是直線，這與實測各物質濃度變化規(guī)律根本不符，所以我想除了模型不合適外還可能是擬合計算方法有著較大的紕漏（因為數(shù)據(jù)擬合這方面從來沒接觸過）。最后還是謝謝您的回復，感謝�。� |
4樓2020-02-02 22:00:26
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調劑 +9 | 王曉陽- 2026-03-09 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 復試調劑 +5 | 呼呼？~+123456 2026-03-08 | 6/300 |
|
|
|
[考研] 0817一志愿蘇州大學280復試調劑 +4 | kk扛 2026-03-07 | 4/200 |
|
|
|
[考研] 086000生物與醫(yī)藥319分求調劑 +4 | Tolkien 2026-03-07 | 4/200 |
|
|
|
[考研] 一志愿東北大學080500，英一數(shù)二343分求調劑 +7 | Seele_v5 2026-03-04 | 7/350 |
|
|
|
[考研] 085701環(huán)境工程專碩求調劑 +5 | 幼稚園大班打手 2026-03-08 | 5/250 |
|
|
|
[考研] 290 材料與化工求調劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
|
[考研] 一志愿武理314求調劑 +4 | （￣～￣;） 2026-03-08 | 5/250 |
|
|
|
[考研] 346分材料求調劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
|
[考研] 求調劑，不管什么專業(yè)，我是可塑造的人才一枚，希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|
|
[考研] 一志愿中國石油大學（華東） 本科齊魯工業(yè)大學 求調劑 +3 | snw石 2026-03-07 | 3/150 |
|
|
|
[考研] 一志愿南京大學資源與環(huán)境專業(yè)283分求調劑 +8 | 追云?? 2026-03-04 | 8/400 |
|
|
|
[考研] 求調劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
|
[考研] 求調劑 +4 | 呼呼？~+123456 2026-03-06 | 4/200 |
|
|
|
[考研] 268求調劑 +4 | 劉合華 2026-03-05 | 4/200 |
|
|
|
[考研] 328求調劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
|
[考研] 復試調劑 +5 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 334求調劑 +3 | Trying] 2026-03-03 | 4/200 |
|
|
|
[考研] 320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗 +5 | 錘某人 2026-03-04 | 8/400 |
|
|
|
[考研] 一志愿武漢理工大學-085602-總分296分-求調劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|