 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4120 | 回復: 26 | ||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
[交流]
動力學參數(shù)擬合
|
||||
|
CO變換反應(yīng)，CO+H2O=CO2+H2，用固定床管式積分反應(yīng)器，要擬合動力學方程中的參數(shù)！動力學方程如下圖。轉(zhuǎn)化后的東西學方程為dy/dw=-22.4./v*k(1)*exp(-k(2)/R.*T).*yco^k(3).*(yh2oin-(ycoin-yco))^k(4).*(yco2in+(ycoin-yco))^k(5).*(yh2in+(ycoin-yco))^k(6).*(1-β)，其中β是逆反應(yīng)程度，要擬合指前因子、反應(yīng)活化能和四個反應(yīng)級數(shù)（k1~k6）。整理后的數(shù)據(jù)是已知反應(yīng)器進口一氧化碳、二氧化碳、氫氣、水蒸氣的濕基組成和反應(yīng)器出口一氧化碳濕基組成，主要是一個微分方程dy/dw=f(y),其中f(y)中含有要擬合的6個參數(shù)，涉及到微分方程的初值問題和最小二乘法問題。自己編了一個m文件，可是運行不了，糾結(jié)好長時間也沒有發(fā)現(xiàn)問題，想大家?guī)臀铱匆幌拢〈a如下 function kinetics clear all clc global t T R kp ycoin yh2oin yco2in yh2in v t=[448.69 449.03 450.18 450.71 450.47 430.01 430.31 430.03 430.17 427.29 410.13 410.08 410.37 410.00 409.96 390.00 390.00 389.99 390.00 390.00];%反應(yīng)溫度 T=t+273.15;%絕對溫度K R=8.314;%氣體常數(shù) kp=exp(4575./T-4.33);%平衡常數(shù) ycoin=[0.09865 0.13296 0.18299 0.23883 0.13809 0.13918 0.15477 0.14327 0.16374 0.13077 0.10624 0.11253 0.19904 0.20900 0.19895 0.12235 0.17431 0.14712 0.18064 0.16864];%進口CO含量 yh2oin=[0.52349 0.45073 0.39759 0.33124 0.45811 0.30844 0.39573 0.51067 0.54215 0.32203 0.48650 0.54179 0.32176 0.42416 0.32176 0.41059 0.32107 0.50444 0.48270 0.31474];%進口H2O含量 yco2in=[0.09440 0.08103 0.07302 0.05747 0.09028 0.11332 0.11931 0.04328 0.05115 0.10654 0.06181 0.03400 0.13821 0.09236 0.13688 0.04807 0.07740 0.07265 0.09696 0.08260];%進口CO2含量 yh2in=[0.09363 0.14281 0.17864 0.22005 0.13013 0.20040 0.20556 0.10088 0.10997 0.20593 0.17479 0.13200 0.16644 0.09602 0.17151 0.14214 0.13702 0.18113 0.16127 0.13905];%進口H2含量 v=[21.3952 18.2224 15.0234 13.5257 19.5599 14.5256 18.8683 19.0081 20.6600 16.2300 18.4185 20.6734 13.9253 15.8444 13.9233 16.3679 13.9538 20.7244 17.0155 14.2337];%總原料流量 ycoout=[0.07304 0.08526 0.12500 0.16207 0.08445 0.10695 0.13233 0.10264 0.11417 0.09878 0.09192 0.08657 0.17122 0.17377 0.17204 0.10501 0.14913 0.13370 0.15282 0.14068];%出口CO含量 k0=[3.5e5 8.4e4 0.9 0.6 -0.5 -0.04]; % k(指前因子、活化能、反應(yīng)級數(shù))初始向量 lb=[1e5 3e4 0 0 -1 -0.1]; ub=[1e6 1.5e5 1 1 0.5 0.5];% k上下限 % 使用函數(shù)lsqnonlin()進行參數(shù)估計 [k,resnorm,residual,exitflag]=lsqnonlin(@ObjFunc,k0,lb,ub,[],ycoin,ycoout); % ------------------------------------------------------------------ function f=ObjFunc(k,ycoin,ycoout) % 目標函數(shù) wspan=[0 0.6091]; % 即催化劑裝填量/g ycow=zeros(1,length(ycoout)); [w yco] =ode45(@Euqations,wspan,ycoin,[],k);%4,5階龍格庫塔法求解 ycow=yco(length(yco));%出口CO含量 f=ycoout-ycow; % --------此方程定義微分方程等號右端項 function dydw=Euqations(w,yco,k) % w催化劑質(zhì)量 global T R kp ycoin yh2oin yco2in yh2in v dydw=-22.4./v*k(1)*exp(-k(2)/R.*T).*yco^k(3).*(yh2oin-(ycoin-yco))^k(4).*(yco2in+(ycoin-yco))^k(5).*(yh2in+(ycoin-yco))^k(6).*(1-(yco2in-(ycoin-yco).*yh2in-(ycoin-yco))./(kp.*yco.*yh2oin-(ycoin-yco)));% 反應(yīng)速率 [ 來自科研家族 化工工程工藝 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 matlab |
» 本帖已獲得的紅花（最新10朵）
dingd» 猜你喜歡
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復
-
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有6人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有287人回復
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復
-
求調(diào)劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復
-
光學工程學碩調(diào)劑信息
已經(jīng)有30人回復
-
歡迎加入課題組
已經(jīng)有0人回復
-
散金幣，求好運，祝面上順利！
已經(jīng)有65人回復
-
華東師范大學芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教matlab反應(yīng)動力學參數(shù)估計遇到的問題，謝謝
已經(jīng)有15人回復
- matlab 求指點 動力學方程擬合過程中導數(shù)的獲取
已經(jīng)有10人回復
- matlab擬合擬均相動力學參數(shù)
已經(jīng)有6人回復
- 最小二乘擬合 初值 運算時間
已經(jīng)有12人回復
- 吸附的動力學方程擬合
已經(jīng)有18人回復
- 固定床中氣-固相脫水反應(yīng)本征動力學
已經(jīng)有5人回復
- 【討論】關(guān)于腐蝕速度與動力學參數(shù)
已經(jīng)有15人回復
- origin，sigmaplot等軟件自定義曲線擬合（主要是一些動力學方程的模擬）
已經(jīng)有2人回復
- 【求助】如何從極化曲線求動力學參數(shù)
已經(jīng)有7人回復
» 搶金幣啦！回帖就可以得到:
-
青島高校2026招生碩士研究生
+5/790
-
中國科學院杭州醫(yī)學研究-常皓研究員課題組博士后招聘
+1/283
-
化學化工學院 招收化工、化學、材料等相關(guān)方向研究生（學碩、專碩都有調(diào)劑名額）
+1/85
-
天津理工大學國家杰青王鐵課題組招收2026年博士研究生
+1/84
-
重慶醫(yī)科大學神經(jīng)科學中心/病生與生理系誠邀海內(nèi)外優(yōu)秀博士加盟博士后隊伍（長期）
+1/81
-
加拿大阿爾伯塔大學招收電磁、無線通信、機器學習方向全獎碩士/博士/博士后/訪問學者
+1/53
-
0854電子信息調(diào)劑，閩南師大光電芯片研發(fā)實驗室，柯少穎教授課題組，集成電路
+1/36
-
0854電子信息調(diào)劑，閩南師大光電芯片研發(fā)實驗室
+1/35
-
上海交大化院功能大分子團隊招2026年碩士、聯(lián)培、博士后
+1/32
-
上海交通大學電催化方向博后招聘-有鋰電、高分子、燃料電池背景者優(yōu)先
+1/31
-
湖南大學微生物結(jié)構(gòu)與功能實驗室2026年計劃招收博士研究生
+1/26
-
湖南大學微生物結(jié)構(gòu)與功能實驗室2026年計劃招收博士研究生
+1/19
-
湖南大學微生物結(jié)構(gòu)與功能實驗室2026年計劃招收博士研究生
+1/16
-
西南某211招27屆博士
+1/16
-
南京林業(yè)大學-國家級青年人才團隊 招2026級申請考核制博士（合成化學、植物化學方向）
+1/14
-
5han一本強校招收材料專業(yè)碩士
+1/11
-
香港大學David Srolovitz教授課題組招聘「香江學者」博士后
+1/6
-
軍科某研究院接收調(diào)劑研究生
+1/3
-
國家杰青低維材料與器件力學團隊2026年招收博士研究生
+1/2
-
福建師范大學化學與材料學院杜克釗團隊博士招生
+2/2
csgt0
榮譽版主 (著名寫手)
-
專家經(jīng)驗: +2 - 計算強帖: 2
- 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 帖子: 2809
- 在線: 515.8小時
- 蟲號: 449046
8樓2012-06-01 10:18:18
9樓2012-06-01 12:01:47
11樓2012-06-01 14:47:56
csgt0
榮譽版主 (著名寫手)
-
專家經(jīng)驗: +2 - 計算強帖: 2
- 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 帖子: 2809
- 在線: 515.8小時
- 蟲號: 449046
★ ★
cenwanglai: 金幣+2, 公式輸入辛苦了~ 2012-06-02 09:45:35
cenwanglai: 金幣+2, 公式輸入辛苦了~ 2012-06-02 09:45:35
|
dydw=-22.4./v*k(1)*exp(-k(2)/R.*T).*yco^k(3).*(yh2oin-(ycoin-yco))^k(4).*(yco2in+(ycoin-yco))^k(5).*(yh2in+(ycoin-yco))^k(6).*(1-(yco2in-(ycoin-yco).*yh2in-(ycoin-yco))./(kp.*yco.*yh2oin-(ycoin-yco)));% 反應(yīng)速率 比如：yco^k(3)得用點的吧 ycoin和yco就一個橫一個豎 |
12樓2012-06-01 15:25:45
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +5 | 曼殊2266 2026-03-18 | 5/250 |
|
|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學材料求調(diào)劑，總分315（英一） +3 | sbdksD 2026-03-19 | 3/150 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 085600材料與化工 +5 | 安全上岸！ 2026-03-16 | 5/250 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
|
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 藥學383 求調(diào)劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 277材料科學與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|