| 8 | 1/1 | 返回列表 |
| 查看: 3068 | 回復: 7 | |||
apolloking金蟲 (小有名氣)
|
[交流]
【求助】使用Matlab擬合反應(yīng)動力學方程問題 已有4人參與
|
|
求助大家一個問題 在固定床積分反應(yīng)器上測定催化劑催化性能 條件如下 給定不同的溫度 和空速 測定相應(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 也就是一個溫度 對于不同的空速有不同的轉(zhuǎn)化率 在試驗中消除了內(nèi)外擴散影響 反應(yīng)器按照等溫操作時理想置換積分反應(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) 文獻中說是用四階的龍格庫塔法沿床層積分上式 采用Nelder-Mead單純形法進行參數(shù)優(yōu)化 參數(shù)優(yōu)化目標函數(shù)如附圖所示 最終求得k0=7.32E9 E為64550J/mol 請問這個用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. |
金蟲 (正式寫手)
至尊木蟲 (著名寫手)
至尊木蟲 (著名寫手)
金蟲 (小有名氣)
至尊木蟲 (著名寫手)
木蟲 (正式寫手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 一志愿中國石油大學(華東) 本科齊魯工業(yè)大學 +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 278求調(diào)劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
[考研] 【考研調(diào)劑】化學專業(yè) 281分,一志愿四川大學,誠心求調(diào)劑 +8 | 吃吃吃才有意義 2026-03-19 | 8/400 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 316求調(diào)劑 +5 | 梁茜雯 2026-03-19 | 5/250 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 0703化學 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 277材料科學與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|