| 8 | 1/1 | 返回列表 |
| 查看: 3069 | 回復(fù): 7 | |||
apolloking金蟲(chóng) (小有名氣)
|
[交流]
【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問(wèn)題 已有4人參與
|
|
求助大家一個(gè)問(wèn)題 在固定床積分反應(yīng)器上測(cè)定催化劑催化性能 條件如下 給定不同的溫度 和空速 測(cè)定相應(yīng)的轉(zhuǎn)化率 如 溫度 空速 轉(zhuǎn)化率 350 5000 0.224067 375 5000 0.308978 400 5000 0.418433 425 5000 0.583276 450 5000 0.712133 475 5000 0.834063 500 5000 0.927528 525 5000 0.971279 550 5000 0.989578 350 10000 0.167333 375 10000 0.26823 400 10000 0.359565 425 10000 0.492762 450 10000 0.605102 475 10000 0.706325 500 10000 0.870149 525 10000 0.936043 550 10000 0.972126 350 15000 0.161258 375 15000 0.229717 400 15000 0.305041 425 15000 0.445914 450 15000 0.562024 475 15000 0.702027 500 15000 0.847299 525 15000 0.901723 550 15000 0.956323 350 20000 0.103751 375 20000 0.143094 400 20000 0.227758 425 20000 0.297666 450 20000 0.456384 475 20000 0.592814 500 20000 0.743624 525 20000 0.825191 550 20000 0.897261 350 25000 0.089297 375 25000 0.149173 400 25000 0.213542 425 25000 0.304452 450 25000 0.433972 475 25000 0.563372 500 25000 0.703383 525 25000 0.801953 550 25000 0.870735 也就是一個(gè)溫度 對(duì)于不同的空速有不同的轉(zhuǎn)化率 在試驗(yàn)中消除了內(nèi)外擴(kuò)散影響 反應(yīng)器按照等溫操作時(shí)理想置換積分反應(yīng)器處理 ra=k0*exp(-E/R/T)*CA0/FA0*(1-x) ra=dx/dW 兩式聯(lián)立 dx/dW=k0*exp(-E/R/T)*CA0/FA0*(1-x) 文獻(xiàn)中說(shuō)是用四階的龍格庫(kù)塔法沿床層積分上式 采用Nelder-Mead單純形法進(jìn)行參數(shù)優(yōu)化 參數(shù)優(yōu)化目標(biāo)函數(shù)如附圖所示 最終求得k0=7.32E9 E為64550J/mol 請(qǐng)問(wèn)這個(gè)用matlab應(yīng)該怎么去做啊? 謝謝 [ Last edited by apolloking on 2009-10-27 at 12:36 ] |
|
There are multiple things involved: 1. Numerical Integration of dX/dW. Check 'help ode15s' 2. Optimization. Simplex method is coded in MATLAB function 'fminsearch'. however, I would suggest you try 'nlinfit' in your case. 3. Pose your problem properly. I would say this is the hardest step and it might take a while. Anyway, 1 and 2 are good starting points in my opinion. |
金蟲(chóng) (正式寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
金蟲(chóng) (小有名氣)
至尊木蟲(chóng) (著名寫(xiě)手)
|
如果用無(wú)梯度反應(yīng)器測(cè)定動(dòng)力學(xué),做粒度試驗(yàn),粒度減小而收率不變則消除了內(nèi)擴(kuò)散;做轉(zhuǎn)速試驗(yàn),轉(zhuǎn)速增加而收率不變則消除了外擴(kuò)散影響。在上述的轉(zhuǎn)速和粒度下完成測(cè)試就可以了。 替樓主回復(fù)了。 [ Last edited by zhdgzhdg on 2010-9-4 at 07:13 ] |
木蟲(chóng) (正式寫(xiě)手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +3 | 順順順mr 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +19 | 不知道叫什么! 2026-03-15 | 21/1050 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考研] 266求調(diào)劑 +5 | 陽(yáng)陽(yáng)哇塞 2026-03-14 | 10/500 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|