| 查看: 465 | 回复: 1 | |||
[交流]
【求助】请帮助解微分方程组 已有1人参与
|
|
% 物料平衡 dxdz = rhoB*r*Mm /(yA0 * MG); % 热量衡算 dTdz = 4*U*(T-TJ)/(dt*MG*Cp)-(rhoB*r*H)/(MG*Cp); 其中r为反应速率， % 反应速度, kmol/kg catalyst hr r=yita*k*kH*PA*PH/(1+kH*PH+KVC*PVC); PA=P0*yA0*(1-x); PH=P0*yH0*(1-x); PVC=P0*x; % 速度常数, kmol/kg catalyst hr if T<403 k=9.21*(10^13)*exp(-27.3*1000/(8.314*T))+0.0476; else k=1.714*(10^14)*exp(-27.55*1000/(8.314*T)); end if T<403 kH=19.32*exp(2.72*1000/(8.314*T)); else kH=1.09*10^20*exp(39.6*1000/(8.314*T)); end KVC=1.783/10^6*exp(12900/(8.314*T))-2.53; yita=0.295*exp(590/T)-120*exp(4.5*x); 对于微分方程组dzdx和dTdx怎么用ODE45来解呀，请高人帮助做这个程序 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有18人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有13人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有22人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有22人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
robert2020(金币+1):多谢解答，辛苦了！ 2010-06-12 20:03:51
robert2020(金币+1):多谢应助，辛苦了！ 2010-06-12 20:06:54
小木虫(金币+0.5):给个红包，谢谢回帖交流
robert2020(金币+1):多谢解答，辛苦了！ 2010-06-12 20:03:51
robert2020(金币+1):多谢应助，辛苦了！ 2010-06-12 20:06:54
|
到下面的地方去下载《实用化工计算机模拟-MATLAB在化学工程中的应用》 你将受益匪浅 http://muchong.com/bbs/viewthread.php?tid=2083386&fpage=1 |
2楼2010-06-05 18:21:11
