| 24小時(shí)熱門(mén)版塊排行榜 |
| 27 | 1/1 | 返回列表 |
| 查看: 4100 | 回復(fù): 26 | ||||
feng198800金蟲(chóng) (正式寫(xiě)手)
|
[交流]
動(dòng)力學(xué)參數(shù)擬合
|
|||
|
CO變換反應(yīng),CO+H2O=CO2+H2,用固定床管式積分反應(yīng)器,要擬合動(dòng)力學(xué)方程中的參數(shù)!動(dòng)力學(xué)方程如下圖。轉(zhuǎn)化后的東西學(xué)方程為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)活化能和四個(gè)反應(yīng)級(jí)數(shù)(k1~k6)。整理后的數(shù)據(jù)是已知反應(yīng)器進(jìn)口一氧化碳、二氧化碳、氫氣、水蒸氣的濕基組成和反應(yīng)器出口一氧化碳濕基組成,主要是一個(gè)微分方程dy/dw=f(y),其中f(y)中含有要擬合的6個(gè)參數(shù),涉及到微分方程的初值問(wèn)題和最小二乘法問(wèn)題。自己編了一個(gè)m文件,可是運(yùn)行不了,糾結(jié)好長(zhǎng)時(shí)間也沒(méi)有發(fā)現(xiàn)問(wèn)題,想大家?guī)臀铱匆幌!代碼如下 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;%絕對(duì)溫度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];%進(jìn)口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];%進(jìn)口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];%進(jìn)口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];%進(jìn)口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)級(jí)數(shù))初始向量 lb=[1e5 3e4 0 0 -1 -0.1]; ub=[1e6 1.5e5 1 1 0.5 0.5];% k上下限 % 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì) [k,resnorm,residual,exitflag]=lsqnonlin(@ObjFunc,k0,lb,ub,[],ycoin,ycoout); % ------------------------------------------------------------------ function f=ObjFunc(k,ycoin,ycoout) % 目標(biāo)函數(shù) wspan=[0 0.6091]; % 即催化劑裝填量/g ycow=zeros(1,length(ycoout)); [w yco] =ode45(@Euqations,wspan,ycoin,[],k);%4,5階龍格庫(kù)塔法求解 ycow=yco(length(yco));%出口CO含量 f=ycoout-ycow; % --------此方程定義微分方程等號(hào)右端項(xiàng) 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)速率 [ 來(lái)自科研家族 化工工程工藝 ] |
matlab |
» 搶金幣啦!回帖就可以得到:
+2/114
+1/87
+1/85
+1/81
+1/69
+1/37
+2/30
+1/20
+1/18
+1/15
+1/11
+1/7
+1/6
+1/5
+2/2
+1/2
+1/2
+1/1
+1/1
+1/1
木蟲(chóng) (著名寫(xiě)手)
榮譽(yù)版主 (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +2 |
鐵桿木蟲(chóng) (職業(yè)作家)
金蟲(chóng) (正式寫(xiě)手)
榮譽(yù)版主 (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +2 |
|
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)得用點(diǎn)的吧 ycoin和yco就一個(gè)橫一個(gè)豎 |
至尊木蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
榮譽(yù)版主 (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +4 |
金蟲(chóng) (正式寫(xiě)手)
榮譽(yù)版主 (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +4 |
|
沒(méi)關(guān)系,我指的是方程形式一樣,在調(diào)用時(shí)分開(kāi),參考 http://www.gaoyang168.com/bbs/viewthread.php?tid=4277851&page=1#pid5 [ Last edited by dbb627 on 2012-6-6 at 16:25 ] |
禁蟲(chóng) (文壇精英)
|
本帖內(nèi)容被屏蔽 |
銅蟲(chóng) (小有名氣)
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
木蟲(chóng) (小有名氣)
|
Optimization terminated: relative function value changing by less than OPTIONS.TolFun. k = 1.0e+005 * Columns 1 through 4 1.000003573678879 0.762988537687594 0.000005843208004 0.000005981921222 Columns 5 through 6 -0.000000744621941 -0.000000996174141 resnorm = 5.657234679833130e-004 residual = -0.006805767672396 0.002960729069534 -0.006244673627901 -0.001713302178392 0.010027828622409 0.004669076015188 -0.006971457219631 -0.002080652235825 0.005813395104906 0.007827253485404 -0.008382737359081 0.000669267176915 -0.001261508186825 -0.002006057821641 -0.001995402177405 -0.001123333955809 0.003628239382343 -0.004182989575642 0.004948094412108 0.007850845335425 exitflag = 3 |
金蟲(chóng) (著名寫(xiě)手)
| 27 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)319求調(diào)劑 +3 | 錦鯉0909 2026-03-17 | 3/150 |
|
|
[考研] 085601專(zhuān)碩,總分342求調(diào)劑,地區(qū)不限 +3 | share_joy 2026-03-16 | 3/150 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[考研] 材料專(zhuān)碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 小熊joy 2026-03-14 | 7/350 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
[考研] 328化工專(zhuān)碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 270求調(diào)劑 085600材料與化工專(zhuān)碩 +3 | YXCT 2026-03-11 | 3/150 |
|
|
[考研] 279求調(diào)劑 +3 | 莫xiao 2026-03-10 | 4/200 |
|