| 查看: 1200 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
哪位大神帮忙看一下MATLAB求解微分方程的问题？ 已有1人参与
|
||
|
clear all clc global L Ct Ac rhoB Cp H1 H2 U dt TJ U TJ us E1 E2 R us L = 4; % 反应管长, m rhoB = 770; % 催化剂堆积密度, kg/m3 rhog =691.82 % 液体混合物密度, kg/m3 CA0=4.136 % 2M1B摩尔流率, moles/hr CB0=37.22 % 2M2B摩尔流率, moles/hr CC0 = 41.36; % MeOH摩尔流率, moles/hr Cp = 131.33; % 比热, kJ/kmol K U = 0.105; % 传热系数, kJ/m2 hr K dt = 0.025; % 管径, m TJ = 333.15; % 冷却温度, K T0 = 343.15; % 物料进口温度(初始温度), K H1 = 33600; % 反应2M1B→T的反应热, kJ/kmol H2 = 26800; % 反应2M2B→T的反应热, kJ/kmol us = 3.11 % 空速，h-1 % 活化能, kJ/kmol E1 =89830 ; E2 =120010 ; R = 8.314; % 理想气体常数, kJ/kmol K Ac = pi*(dt/2)^2; % 反应管的横截面积, m2 [z, y] = ode45(@Equations, [0 L], [CA0 CB0 CC0 T0]) CA = y(:, 1); CB = y(:, 2); CC = y(:, 3); xA = (CA0-CA)./CA0; % A的转化率 xB = (CB0-CB)./CB0; % B的转化率 xB = [0; xB] xC = (CC0-CC)./CC0; % C的转化率 xC = [0; xC] % 图形输出 plot(z, y(:, 4)) % 温度分布 xlabel('z') ylabel('T (K)') figure plot(z, xA, 'r-') % 转化率分布 xlabel('z') ylabel('x_A') figure plot(z, CA, 'r-', z, CB, 'k--', z, CC, 'b:') % 浓度分布 xlabel('z') ylabel('C_A, C_B, C_C') legend('C_A', 'C_B', 'C_C') function [rA, rB, rC, k1, k2,K1,K2] = Rates(y(1), y(2), y(3), T) % 反应动力学 global E1 E2 R CC0 % 速度常数, kmol/kg catalyst hr k1 = exp(-E1/(R*T) + 35.35); k2 = exp(-E2/(R*T) + 42.47); % 反应平衡常数 K1= exp(-8.3881+4041.2/T); K2= exp(-8.2473+3225.3/T); % 反应速度, kmol/kg catalyst hr rA = -k1*CA*CC + k1/K1*(CC0-CC) % 2M1B的消耗速率 rB = -k2*CB*CC + k2/K2*(CC0-CC); % 2M2B的消耗速率 rC = -k1*CA*CC + k1/K1*(CC0-CC)-k2*CB*CC + k2/K2*(CC0-CC); % 甲醇的总反应速度 function dydz = Equations(z, y) % 模型方程组 global Ct Ac rhoB Cp H1 H2 U dt TJ U CA = y(1); % A---2M1B CB = y(2); % B---2M2B CC = y(3); % C---MeOH T = y(4); % 反应速度 [rA, rB, rC, k1, k2,K1,K2] = Rates(y(1), y(2), y(3), T) % 物料平衡 dCAdz = rhoB*rA/us; dCBdz = rhoB*rB/us; dCCdz = rhoB*rC/us; % 热量衡算 dTdz = ( rhoB*((-H1*k1*y(1)*y(3))-(-H2*k2*y(2)*y(3)))-4*U*(T-TJ)/dt )/(us*rhog*Cp); dydz = [dCAdz; dCBdz; dCCdz; dTdz]; |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有201人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于病态方程组求解的疑问
已经有5人回复
- 一个偏微分方程组的求解问题
已经有5人回复
- Matlab求解微分方程(组)及偏微分方程(组)
已经有224人回复
- matlab 非线性微分方程求解
已经有3人回复
- 以下matlab代码出现问题
已经有10人回复
- 偏微分方程求解！
已经有4人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 如何用Matlab求时变微分方程的解？
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- Matlab求解二阶偏微分方程组，希望能给出相应的求解代码
已经有10人回复
- 大神进，数学问题。。求解。
已经有7人回复
- 用matlab求解方程出问题，请帮忙看看
已经有3人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- matlab解微分方程组
已经有15人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
|
算的时候还是出现错误，烦劳您再看看。 ??? Error using ==> feval Undefined function or method 'Equations' for input arguments of type 'double'. Error in ==> odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... |
3楼2014-02-18 10:30:35